Parsing the roles of neck-linker docking and tethered head diffusion in the stepping dynamics of kinesin

Competition between physical search and a weak-to-strong transition rate-limits kinesin binding times

The self-organization of cells relies on the profound complexity of protein-protein interactions. Challenges in directly observing these events have hindered progress toward understanding their diverse behaviors. One notable example is the interaction between molecular motors and cytoskeletal systems that combine to perform a variety of cellular functions. In this work, we leverage theory and experiments to identify and quantify the rate-limiting mechanism of the initial association between a cargo-bound kinesin motor and a microtubule track. Recent advances in optical tweezers provide binding times for several lengths of kinesin motors trapped at varying distances from a microtubule, empowering the investigation of competing models. We first explore a diffusion-limited model of binding. Through Brownian dynamics simulations and simulation-based inference, we find this simple diffusion model fails to explain the experimental binding times, but an extended model that accounts for the ADP state of the molecular motor agrees closely with the data, even under the scrutiny of penalizing for additional model complexity. We provide quantification of both kinetic rates and biophysical parameters underlying the proposed binding process. Our model suggests that a typical binding event is limited by ADP state rather than physical search. Lastly, we predict how these association rates can be modulated in distinct ways through variation of environmental concentrations and physical properties.

ADP release can explain spatially-dependent kinesin binding times

The self-organization of cells relies on the profound complexity of protein-protein interactions. Challenges in directly observing these events have hindered progress toward understanding their diverse behaviors. One notable example is the interaction between molecular motors and cytoskeletal systems that combine to perform a variety of cellular functions. In this work, we leverage theory and experiments to identify and quantify the rate-limiting mechanism of the initial association between a cargo-bound kinesin motor and a microtubule track. Recent advances in optical tweezers provide binding times for several lengths of kinesin motors trapped at varying distances from a microtubule, empowering the investigation of competing models. We first explore a diffusion-limited model of binding. Through Brownian dynamics simulations and simulation-based inference, we find this simple diffusion model fails to explain the experimental binding times, but an extended model that accounts for the ADP state of the molecular motor agrees closely with the data, even under the scrutiny of penalizing for additional model complexity. We provide quantification of both kinetic rates and biophysical parameters underlying the proposed binding process. Our model suggests that most but not every motor binding event is limited by their ADP state. Lastly, we predict how these association rates can be modulated in distinct ways through variation of environmental concentrations and spatial distances.

Regulatory mechanisms of the dynein-2 motility by post-translational modification revealed by MD simulation

Intraflagellar transport for ciliary assembly and maintenance is driven by dynein and kinesins specific to the cilia. It has been shown that anterograde and retrograde transports run on different regions of the doublet microtubule, i.e., separate train tracks. However, little is known about the regulatory mechanism of this selective process. Since the doublet microtubule is known to display specific post-translational modifications of tubulins, i.e., "tubulin code", for molecular motor regulations, we investigated the motility of ciliary specific dynein-2 under different post-translational modification by coarse-grained molecular dynamics. Our setup allows us to simulate the landing behaviors of dynein-2 on un-modified, detyrosinated, poly-glutamylated and poly-glycylated microtubules in silico. Our study revealed that poly-glutamylation can play an inhibitory effect on dynein-2 motility. Our result indicates that poly-glutamylation of the B-tubule of the doublet microtubule can be used as an efficient means to target retrograde intraflagellar transport onto the A-tubule.

Road-blocker HSP disease mutation disrupts pre-organization for ATP hydrolysis in kinesin through a second sphere control

Kinesin motor proteins perform several essential cellular functions powered by the adenosine triphosphate (ATP) hydrolysis reaction. Several single-point mutations in the kinesin motor protein KIF5A have been implicated to hereditary spastic paraplegia disease (HSP), a lethal neurodegenerative disease in humans. In earlier studies, we have shown that a series of HSP-related mutations can impair the kinesin's long-distance displacement or processivity by modulating the order-disorder transition of the linker connecting the heads to the coiled coil. On the other hand, the reduction of kinesin's ATP hydrolysis reaction rate by a distal asparagine-to-serine mutation is also known to cause HSP disease. However, the molecular mechanism of the ATP hydrolysis reaction in kinesin by this distal mutation is still not fully understood. Using classical molecular dynamics simulations combined with quantum mechanics/molecular mechanics calculations, the pre-organization geometry required for optimal hydrolysis in kinesin motor bound to α/β-tubulin is determined. This optimal geometry has only a single salt-bridge (of the possible two) between Arg203-Glu236, putting a reactive water molecule at a perfect position for hydrolysis. Such geometry is also needed to create the appropriate configuration for proton translocation during ATP hydrolysis. The distal asparagine-to-serine mutation is found to disrupt this optimal geometry. Therefore, the current study along with our previous one demonstrates how two different effects on kinesin dynamics (processivity and ATP hydrolysis), caused by a different set of genotypes, can give rise to the same phenotype leading to HSP disease.

Plus and minus ends of microtubules respond asymmetrically to kinesin binding by a long-range directionally driven allosteric mechanism

Although it is known that majority of kinesin motors walk predominantly toward the plus end of microtubules (MTs) in a hand-over-hand manner, the structural origin of the stepping directionality is not understood. To resolve this issue, we modeled the structures of kinesin-1 (Kin1), MT, and the Kin1-MT complex using the elastic network model and calculated the residue-dependent responses to a local perturbation in the constructs. Kin1 binding elicits an asymmetric response that is pronounced in α/β-tubulin dimers in the plus end of the MT. Kin1 opens the clefts of multiple plus end α/β-tubulin dimers, creating binding-competent conformations, which is required for processivity. Reciprocally, MT induces correlations between switches I and II in the motor and enhances fluctuations in adenosine 5'-diphosphate and the residues in the binding pocket. Our findings explain both the directionality of stepping and MT effects on a key step in the catalytic cycle of kinesin.

Theory and simulations of condensin mediated loop extrusion in DNA

Condensation of hundreds of mega-base-pair-long human chromosomes in a small nuclear volume is a spectacular biological phenomenon. This process is driven by the formation of chromosome loops. The ATP consuming motor, condensin, interacts with chromatin segments to actively extrude loops. Motivated by real-time imaging of loop extrusion (LE), we created an analytically solvable model, predicting the LE velocity and step size distribution as a function of external load. The theory fits the available experimental data quantitatively, and suggests that condensin must undergo a large conformational change, induced by ATP binding, bringing distant parts of the motor to proximity. Simulations using a simple model confirm that the motor transitions between an open and a closed state in order to extrude loops by a scrunching mechanism, similar to that proposed in DNA bubble formation during bacterial transcription. Changes in the orientation of the motor domains are transmitted over ~50 nm, connecting the motor head and the hinge, thus providing an allosteric basis for LE.

Effects of Gold Nanoparticles on the Stepping Trajectories of Kinesin

A substantial increase in the temporal resolution of the stepping of dimeric molecular motors is possible by tracking the position of a large gold nanoparticle (GNP) attached to a labeled site on one of the heads. This technique was employed to measure the stepping trajectories of conventional kinesin (Kin1) using the time-dependent position of the GNP as a proxy. The trajectories revealed that the detached head always passes to the right of the head that is tightly bound to the microtubule (MT) during a step. In interpreting the results of such experiments, it is assumed that the GNP does not significantly alter the diffusive motion of the detached head. We used coarse-grained simulations of a system consisting of the MT-Kin1 complex with and without attached GNP to investigate how the stepping trajectories are affected. The two significant findings are: (1) The GNP does not faithfully track the position of the stepping head, and (2) the rightward bias is typically exaggerated by the GNP. Both these findings depend on the precise residue position to which the GNP is attached. Surprisingly, the stepping trajectories of kinesin are not significantly affected if, in addition to the GNP, a 1 μm diameter cargo is attached to the coiled coil. Our simulations suggest the effects of the large probe have to be considered when inferring the stepping mechanisms using GNP tracking experiments.

Multiscale Coarse-Grained Model for the Stepping of Molecular Motors with Application to Kinesin

Conventional kinesin, a motor protein that transports cargo within cells, walks by taking multiple steps toward the plus end of the microtubule (MT). While significant progress has been made in understanding the details of the walking mechanism of kinesin, there are many unresolved issues. From a computational perspective, a central challenge is the large size of the system, which limits the scope of time scales accessible in standard computer simulations. Here, we create a general multiscale coarse-grained model for motors that enables us to simulate the stepping process of motors on polar tracks (actin and MT) with a focus on kinesin. Our approach greatly shortens the computation times without a significant loss in detail, thus allowing us to better describe the molecular basis of the stepping kinetics. The small number of parameters, which are determined by fitting to experimental data, allows us to develop an accurate method that may be adopted to simulate stepping in other molecular motors. The model enables us to simulate a large number of steps, which was not possible previously. We show in agreement with experiments that due to the docking of the neck linker (NL) of kinesin, sometimes deemed as the power stroke, the space explored diffusively by the tethered head is severely restricted, allowing the step to be completed in tens of microseconds. We predict that increasing the interaction strength between the NL and the motor head, achievable by mutations in the NL, decreases the stepping time but reaches a saturation value. Furthermore, the full three-dimensional dynamics of the cargo are fully resolved in our model, contributing to the predictive power and allowing us to study the important aspects of cargo-motor interactions.

Computational modeling of dynein motor proteins at work

Along with various experimental methods, a combination of theoretical and computational methods is essential to explore different length-scale and time-scale processes in the biological system. The functional mechanism of a dynein, an ATP-fueled motor protein, working in a multiprotein complex, involves a wide range of length/time-scale events. It generates mechanical force from chemical energy and moves on microtubules towards the minus end direction while performing a large number of biological processes including ciliary beating, intracellular material transport, and cell division. Like in the cases of other conventional motor proteins, a combination of experimental techniques including X-crystallography, cryo-electron microscopy, and single molecular assay have provided a wealth of information about the mechanochemical cycle of a dynein. Dyneins have a large and complex structural architecture and therefore, computational modeling of different aspects of a dynein is extremely challenging. As the process of dynein movement involves varying length and timescales, it demands, like in experiments, a combination of computational methods covering such a wide range of processes for the comprehensive investigation of the mechanochemical cycle. In this review article, we will summarize how the use of state-of-the-art computational methods can provide a detailed molecular understanding of the mechanochemical cycle of the dynein. We implemented all-atom molecular dynamics simulations and hybrid quantum-mechanics/molecular-mechanics simulations to explore the ATP hydrolysis mechanisms at the primary ATPase site (AAA1) of dynein. To investigate the large-scale conformational changes we employed coarse-grained structure-based molecular dynamics simulations to capture the domain motions. Here we explored the conformational changes upon binding of ATP at AAA1, nucleotide state-dependent regulation of the mechanochemical cycle, and inter-head coordination by inter-head tension. Additionally, implementing a phenomenological theoretical model we explore the force-dependent detachment rate of a motorhead from the microtubule and the principle of multi-dynein cooperation during cargo transport.

Mechanistic basis of propofol-induced disruption of kinesin processivity

Propofol is a widely used general anesthetic to induce and maintain anesthesia, and its effects are thought to occur through impact on the ligand-gated channels including the GABA_A receptor. Propofol also interacts with a large number of proteins including molecular motors and inhibits kinesin processivity, resulting in significant decrease in the run length for conventional kinesin-1 and kinesin-2. However, the molecular mechanism by which propofol achieves this outcome is not known. The structural transition in the kinesin neck-linker region is crucial for its processivity. In this study, we analyzed the effect of propofol and its fluorine derivative (fropofol) on the transition in the neck-linker region of kinesin. Propofol binds at two crucial surfaces in the leading head: one at the microtubule-binding interface and the other in the neck-linker region. We observed in both the cases the order-disorder transition of the neck-linker was disrupted and kinesin lost its signal for forward movement. In contrast, there was not an effect on the neck-linker transition with propofol binding at the trailing head. Free-energy calculations show that propofol at the microtubule-binding surface significantly reduces the microtubule-binding affinity of the kinesin head. While propofol makes pi-pi stacking and H-bond interactions with the propofol binding cavity, fropofol is unable to make a suitable interaction at this binding surface. Therefore, the binding affinity of fropofol is much lower compared to propofol. Hence, this study provides a mechanism by which propofol disrupts kinesin processivity and identifies transitions in the ATPase stepping cycle likely affected.

Protein Diffusion on Charged Biopolymers: DNA versus Microtubule

Protein diffusion in lower-dimensional spaces is used for various cellular functions. For example, sliding on DNA is essential for proteins searching for their target sites, and protein diffusion on microtubules is important for proper cell division and neuronal development. On the one hand, these linear diffusion processes are mediated by long-range electrostatic interactions between positively charged proteins and negatively charged biopolymers and have similar characteristic diffusion coefficients. On the other hand, DNA and microtubules have different structural properties. Here, using computational approaches, we studied the mechanism of protein diffusion along DNA and microtubules by exploring the diffusion of both protein types on both biopolymers. We found that DNA-binding and microtubule-binding proteins can diffuse on each other's substrates; however, the adopted diffusion mechanism depends on the molecular properties of the diffusing proteins and the biopolymers. On the protein side, only DNA-binding proteins can perform rotation-coupled diffusion along DNA, with this being due to their higher net charge and its spatial organization at the DNA recognition helix. By contrast, the lower net charge on microtubule-binding proteins enables them to diffuse more quickly than DNA-binding proteins on both biopolymers. On the biopolymer side, microtubules possess intrinsically disordered, negatively charged C-terminal tails that interact with microtubule-binding proteins, thus supporting their diffusion. Thus, although both DNA-binding and microtubule-binding proteins can diffuse on the negatively charged biopolymers, the unique molecular features of the biopolymers and of their natural substrates are essential for function.

Tubulin tails and their modifications regulate protein diffusion on microtubules

Microtubules (MTs) are essential components of the eukaryotic cytoskeleton that serve as "highways" for intracellular trafficking. In addition to the well-known active transport of cargo by motor proteins, many MT-binding proteins seem to adopt diffusional motility as a transportation mechanism. However, because of the limited spatial resolution of current experimental techniques, the detailed mechanism of protein diffusion has not been elucidated. In particular, the precise role of tubulin tails and tail modifications in the diffusion process is unclear. Here, using coarse-grained molecular dynamics simulations validated against atomistic simulations, we explore the molecular mechanism of protein diffusion along MTs. We found that electrostatic interactions play a central role in protein diffusion; the disordered tubulin tails enhance affinity but slow down diffusion, and diffusion occurs in discrete steps. While diffusion along wild-type MT is performed in steps of dimeric tubulin, the removal of the tails results in a step of monomeric tubulin. We found that the energy barrier for diffusion is larger when diffusion on MTs is mediated primarily by the MT tails rather than the MT body. In addition, globular proteins (EB1 and PRC1) diffuse more slowly than an intrinsically disordered protein (Tau) on MTs. Finally, we found that polyglutamylation and polyglycylation of tubulin tails lead to slower protein diffusion along MTs, although polyglycylation leads to faster diffusion across MT protofilaments. Taken together, our results explain experimentally observed data and shed light on the roles played by disordered tubulin tails and tail modifications in the molecular mechanism of protein diffusion along MTs.

Influenza as a molecular walker

The surface of the influenza virus is decorated with the receptor-binding protein hemagglutinin (HA) and the receptor-cleaving enzyme neuraminidase (NA). HA is responsible for host cell recognition, while NA prevents aggregation and entrapment, but the intricate mechanism of how the functions of these glycoproteins cooperate and how they are regulated by mutational responses to environmental pressures remains unclear. Recently, several groups have described the motion of influenza over surfaces and reported that this motion is inhibited by NA inhibitors. We argue that the motion of influenza resembles the motility of artificial receptor-cleaving particles called "molecular spiders". The cleaving of receptors by this type of molecular walkers leads to self-avoiding motion across a surface. When the binding and cleaving rates of molecular spiders are balanced, they move both rapidly and efficiently. The studies of molecular spiders offer new insights into the functional balance of HA and NA, but they do not address the asymmetric distribution of HA and NA on the surface of influenza. We propose that receptor-cleaving molecular walkers could play an important role in the further investigation of the motility of influenza viruses.

Molecular dynamics simulation of the nanosecond pulsed electric field effect on kinesin nanomotor

Kinesin is a biological molecular nanomotor which converts chemical energy into mechanical work. To fulfill various nanotechnological tasks in engineered environments, the function of biological molecular motors can be altered by artificial chemical modifications. The drawback of this approach is the necessity of designing and creating a new motor construct for every new task. We propose that intense nanosecond-scale pulsed electric field could modify the function of nanomotors. To explore this hypothesis, we performed molecular dynamics simulation of a kinesin motor domain docked on a subunit of its microtubule track - a single tubulin heterodimer. In the simulation, we exposed the kinesin motor domain to intense (100 MV/m) electric field up to 30 ns. We found that both the magnitude and angle of the kinesin dipole moment are affected. Furthermore, we found that the electric field affects contact surface area between kinesin and tubulin, the structure and dynamics of the functionally important kinesin segments, including microtubule binding motifs as well as nucleotide hydrolysis site which power the nanomotor. These findings indicate that external intense nanosecond-scale electric field could alter kinesin behavior. Our results contribute to developing novel electromagnetic methods for modulating the function of biomolecular matter at the nanoscale.

How kinesin waits for ATP affects the nucleotide and load dependence of the stepping kinetics

Conventional kinesin, responsible for directional transport of cellular vesicles, takes multiple nearly uniform 8.2-nm steps by consuming one ATP molecule per step as it walks toward the plus end of the microtubule (MT). Despite decades of intensive experimental and theoretical studies, there are gaps in the elucidation of key steps in the catalytic cycle of kinesin. How the motor waits for ATP to bind to the leading head is controversial. Two experiments using a similar protocol have arrived at different conclusions. One asserts that kinesin waits for ATP in a state with both the heads bound to the MT, whereas the other shows that ATP binds to the leading head after the trailing head detaches. To discriminate between the 2 scenarios, we developed a minimal model, which analytically predicts the outcomes of a number of experimental observable quantities (the distribution of run length, the distribution of velocity [[Formula: see text]], and the randomness parameter) as a function of an external resistive force (F) and ATP concentration ([T]). The differences in the predicted bimodality in [Formula: see text] as a function of F between the 2 models may be amenable to experimental testing. Most importantly, we predict that the F and [T] dependence of the randomness parameters differ qualitatively depending on the waiting states. The randomness parameters as a function of F and [T] can be quantitatively measured from stepping trajectories with very little prejudice in data analysis. Therefore, an accurate measurement of the randomness parameter and the velocity distribution as a function of load and nucleotide concentration could resolve the apparent controversy.

All-atom molecular dynamics simulations reveal how kinesin transits from one-head-bound to two-heads-bound state

Kinesin dimer walks processively along a microtubule (MT) protofilament in a hand-over-hand manner, transiting alternately between one-head-bound (1HB) and two-heads-bound (2HB) states. In 1HB state, one head bound by adenosine diphosphate (ADP) is detached from MT and the other head is bound to MT. Here, using all-atom molecular dynamics simulations we determined the position and orientation of the detached ADP-head relative to the MT-bound head in 1HB state. We showed that in 1HB state when the MT-bound head is in ADP or nucleotide-free state, with its neck linker being undocked, the detached ADP-head and the MT-bound head have the high binding energy, and after adenosine triphosphate (ATP) binds to the MT-bound head, with its neck linker being docked, the binding energy between the two heads is reduced greatly. These results reveal how the kinesin dimer retains 1HB state before ATP binding and how the dimer transits from 1HB to 2HB state after ATP binding. Key residues involved in the head-head interaction in 1HB state were identified.

Mouse embryo geometry drives formation of robust signaling gradients through receptor localization

Morphogen signals are essential for cell fate specification during embryogenesis. Some receptors that sense these morphogens are known to localize to only the apical or basolateral membrane of polarized cell lines in vitro. How such localization affects morphogen sensing and patterning in the developing embryo remains unknown. Here, we show that the formation of a robust BMP signaling gradient in the early mouse embryo depends on the restricted, basolateral localization of BMP receptors. The mis-localization of receptors to the apical membrane results in ectopic BMP signaling in the mouse epiblast in vivo. With evidence from mathematical modeling, human embryonic stem cells in vitro, and mouse embryos in vivo, we find that the geometric compartmentalization of BMP receptors and ligands creates a signaling gradient that is buffered against fluctuations. Our results demonstrate the importance of receptor localization and embryo geometry in shaping morphogen signaling during embryogenesis.

Structural basis for power stroke vs. Brownian ratchet mechanisms of motor proteins

Two mechanisms have been proposed for the function of motor proteins: The power stroke and the Brownian ratchet. The former refers to generation of a large downhill free energy gradient over which the motor protein moves nearly irreversibly in making a step, whereas the latter refers to biasing or rectifying the diffusive motion of the motor. Both mechanisms require input of free energy, which generally involves the processing of an ATP (adenosine 5'-triphosphate) molecule. Recent advances in experiments that reveal the details of the stepping motion of motor proteins, together with computer simulations of atomistic structures, have provided greater insights into the mechanisms. Here, we compare the various models of the power stroke and the Brownian ratchet that have been proposed. The 2 mechanisms are not mutually exclusive, and various motor proteins employ them to different extents to perform their biological function. As examples, we discuss linear motor proteins Kinesin-1 and myosin-V, and the rotary motor F₁-ATPase, all of which involve a power stroke as the essential element of their stepping mechanism.

Theory of Nonequilibrium Free Energy Transduction by Molecular Machines

Biomolecular machines are protein complexes that convert between different forms of free energy. They are utilized in nature to accomplish many cellular tasks. As isothermal nonequilibrium stochastic objects at low Reynolds number, they face a distinct set of challenges compared with more familiar human-engineered macroscopic machines. Here we review central questions in their performance as free energy transducers, outline theoretical and modeling approaches to understand these questions, identify both physical limits on their operational characteristics and design principles for improving performance, and discuss emerging areas of research.

Insights into Kinesin-1 Stepping from Simulations and Tracking of Gold Nanoparticle-Labeled Motors

High-resolution tracking of gold nanoparticle-labeled proteins has emerged as a powerful technique for measuring the structural kinetics of processive enzymes and other biomacromolecules. These techniques use point spread function (PSF) fitting methods borrowed from single-molecule fluorescence imaging to determine molecular positions below the diffraction limit. However, compared to fluorescence, gold nanoparticle tracking experiments are performed at significantly higher frame rates and utilize much larger probes. In the current work, we use Brownian dynamics simulations of nanoparticle-labeled proteins to investigate the regimes in which the fundamental assumptions of PSF fitting hold and where they begin to break down. We find that because gold nanoparticles undergo tethered diffusion around their anchor point, PSF fitting cannot be extended to arbitrarily fast frame rates. Instead, camera exposure times that allow the nanoparticle to fully populate its stationary positional distribution achieve a spatial averaging that increases fitting precision. We furthermore find that changes in the rotational freedom of the tagged protein can lead to artifactual translations in the fitted particle position. Finally, we apply these lessons to dissect a standing controversy in the kinesin field over the structure of a dimer in the ATP waiting state. Combining new experiments with simulations, we determine that the rear kinesin head in the ATP waiting state is unbound but not displaced from its previous microtubule binding site and that apparent differences in separately published reports were simply due to differences in the gold nanoparticle attachment position. Our results highlight the importance of gold conjugation decisions and imaging parameters to high-resolution tracking results and will serve as a useful guide for the design of future gold nanoparticle tracking experiments.

Structural consequences of hereditary spastic paraplegia disease-related mutations in kinesin

A wide range of mutations in the kinesin motor Kif5A have been linked to a neuronal disorder called hereditary spastic paraplegia (HSP). The position of these mutations can vary, and a range of different motile behaviors have been observed, indicating that the HSP mutants can alter distinct aspects of kinesin mechanochemistry. While focusing on four key HSP-associated mutants, this study examined the structural and dynamic perturbations that arise from these mutations using a series of different computational methods, ranging from bioinformatics analyses to all-atom simulations, that account for solvent effects explicitly. We show that two catalytic domain mutations (R280S and K253N) reduce the microtubule (MT) binding affinity of the kinesin head domains appreciably, while N256S has a much smaller impact. Bioinformatics analysis suggests that the stalk mutation A361V perturbs motor dimerization. Subsequent integration of these effects into a coarse-grained structure-based model of dimeric kinesin revealed that the order-disorder transition of the neck linker is substantially affected, indicating a hampered directionality and processivity of kinesin. The present analyses therefore suggest that, in addition to kinesin-MT binding and coiled-coil dimerization, HSP mutations affecting motor stepping transitions and processivity can lead to disease.

Antitubulin sulfonamides: The successful combination of an established drug class and a multifaceted target

Tubulin, the microtubules and their dynamic behavior are amongst the most successful antitumor, antifungal, antiparasitic, and herbicidal drug targets. Sulfonamides are exemplary drugs with applications in the clinic, in veterinary and in the agrochemical industry. This review summarizes the actual state and recent progress of both fields looking from the double point of view of the target and its drugs, with special focus onto the structural aspects. The article starts with a brief description of tubulin structure and its dynamic assembly and disassembly into microtubules and other polymers. Posttranslational modifications and the many cellular means of regulating and modulating tubulin's biology are briefly presented in the tubulin code. Next, the structurally characterized drug binding sites, their occupying drugs and the effects they induce are described, emphasizing on the structural requirements for high potency, selectivity, and low toxicity. The second part starts with a summary of the favorable and highly tunable combination of physical-chemical and biological properties that render sulfonamides a prototypical example of privileged scaffolds with representatives in many therapeutic areas. A complete description of tubulin-binding sulfonamides is provided, covering the different species and drug sites. Some of the antimitotic sulfonamides have met with very successful applications and others less so, thus illustrating the advances, limitations, and future perspectives of the field. All of them combine in a mechanism of action and a clinical outcome that conform efficient drugs.

Mechanochemical Model of the Power Stroke of the Single-Headed Motor Protein KIF1A

During the process of ATP binding to an apo-kinesin microtubule (MT), the kinesin core rotates on the MT, and the neck linker (NL) of the kinesin undergoes an undocked to docked transition. This has been suggested to be a power stroke of kinesin, on the basis of the structural analysis. Here, we developed a mesoscopic structure-based model and studied the power stroke of KIF1A. We quantified the underlying free energy landscape and showed the emergence of several states for the power stroke of KIF1A: UB-UR-UD (unbound, unrotating, undock), B-IR-UD (bound, initial rotating, undock), B-PR-UD (bound, partial rotating, undock), and B-R-D (bound, rotating, dock). We found that ATP binding triggered conformational fluctuations of key elements. We also explored the conformational change of key structural elements during the rotation of KIF1A and docking of the NL. In addition, we semiquantitatively and qualitatively estimated the free energy released by the ATP binding, and how much of this remains for the docking of the NL during the power stroke process at different temperatures. Finally, based on results from the thermodynamics landscape and conformational change of structural key elements, we proposed a mechanochemical model of the power stroke of KIF1A.

Interface Residues That Drive Allosteric Transitions Also Control the Assembly of l-Lactate Dehydrogenase

The allosteric enzyme, l-lactate dehydrogenase (LDH), is activated by fructose 1,6-metaphosphate (FBP) to reduce pyruvate to lactate. The molecular details of the FBP-driven transition from the low affinity T state to the high affinity R state in LDH, a tetramer composed of identical subunits, are not known. The dynamics of the T → R allosteric transition, investigated using Brownian dynamics (BD) simulations of the self-organized polymer (SOP) model, revealed that coordinated rotations of the subunits drive the T → R transition. We used the structural perturbation method (SPM), which requires only the static structure, to identify the allostery wiring diagram (AWD), a network of residues that transmits signals across the tetramer, as LDH undergoes the T → R transition. Interestingly, the residues that play a major role in the dynamics, which are predominantly localized at the interfaces, coincide with the AWD identified using the SPM. Although the allosteric pathways are highly heterogeneous, on the basis of our simulations, we surmise that predominantly the conformational changes in the T → R transition start from the region near the active site, comprised of helix αC, helix α1/2G, helix α3G, and helix α2F, and proceed to other structural units, thus completing the global motion. Brownian dynamics simulations of the tetramer assembly, triggered by a temperature quench from the fully disrupted conformations, show that the bottleneck for assembly is the formation of the correct orientational registry between the subunits, requiring contacts between the interface residues. Surprisingly, these residues are part of the AWD, which was identified using the SPM. Taken together, our results show that LDH, and perhaps other multidomain proteins, may have evolved to stabilize distinct states of allosteric enzymes using precisely the same AWD that also controls the functionally relevant allosteric transitions.

Directionally biased sidestepping of Kip3/kinesin-8 is regulated by ATP waiting time and motor-microtubule interaction strength

Kinesin-8 motors, which move in a highly processive manner toward microtubule plus ends where they act as depolymerases, are essential regulators of microtubule dynamics in cells. To understand their navigation strategy on the microtubule lattice, we studied the 3D motion of single yeast kinesin-8 motors, Kip3, on freely suspended microtubules in vitro. We observed short-pitch, left-handed helical trajectories indicating that kinesin-8 motors frequently switch protofilaments in a directionally biased manner. Intriguingly, sidestepping was not directly coupled to forward stepping but rather depended on the average dwell time per forward step under limiting ATP concentrations. Based on our experimental findings and numerical simulations we propose that effective sidestepping toward the left is regulated by a bifurcation in the Kip3 step cycle, involving a transition from a two-head-bound to a one-head-bound conformation in the ATP-waiting state. Results from a kinesin-1 mutant with extended neck linker hint toward a generic sidestepping mechanism for processive kinesins, facilitating the circumvention of intracellular obstacles on the microtubule surface.

Theoretical Investigations of the Role of Mutations in Dynamics of Kinesin Motor Proteins

Motor proteins are active enzymatic molecules that are critically important for a variety of biological phenomena. It is known that some neurodegenerative diseases are caused by specific mutations in motor proteins that lead to their malfunctioning. Hereditary spastic paraplegia is one of such diseases, and it is associated with the mutations in the neuronal conventional kinesin gene, producing the decreased speed and processivity of this motor protein. Despite the importance of this problem, there is no clear understanding on the role of mutations in modifying dynamic properties of motor proteins. In this work, we investigate theoretically the molecular basis for negative effects of two specific mutations, N256S and R280S, on the dynamics of kinesin motor proteins. We hypothesize that these mutations might accelerate the adenosine triphosphate (ATP) release by increasing the probability of open conformations for the ATP-binding pocket. Our approach is based on the use of coarse-grained structure-based molecular dynamics simulations to analyze the conformational changes and chemical transitions in the kinesin molecule, which is also supplemented by investigation of a mesoscopic discrete-state stochastic model. Computer simulations suggest that mutations N256S and R280S can decrease the free energy difference between open and closed biochemical states, making the open conformation more stable and the ATP release faster, which is in agreement with our hypothesis. Furthermore, we show that in the case of N256S mutation, this effect is caused by disruption of interactions between α helix and switch I and loop L11 structural elements. Our computational results are qualitatively supported by the explicit analysis of the discrete-state stochastic model.