Run length distribution of dimerized kinesin-3 molecular motors: comparison with dimeric kinesin-1.

Kinesin-3 and kinesin-1 molecular motors are two families of the kinesin superfamily. It has been experimentally revealed that in monomeric state kinesin-3 is inactive in motility and cargo-mediated dimerization results in superprocessive motion, with an average run length being more than 10-fold longer than that of kinesin-1. In contrast to kinesin-1 showing normally single-exponential distribution of run lengths, dimerized kinesin-3 shows puzzlingly Gaussian distribution of run lengths. Here, based on our proposed model, we studied computationally the dynamics of kinesin-3 and compared with that of kinesin-1, explaining quantitatively the available experimental data and revealing the origin of superprocessivity and Gaussian run length distribution of kinesin-3. Moreover, predicted results are provided on ATP-concentration dependence of run length distribution and force dependence of mean run length and dissociation rate of kinesin-3.

Reexamining the origin of the directionality of myosin V.

The nature of the conversion of chemical energy to directional motion in myosin V is examined by careful simulations that include two complementary methods: direct Langevin Dynamics (LD) simulations with a scaled-down potential that provided a detailed time-resolved mechanism, and kinetic equations solution for the ensemble long-time propagation (based on information collected for segments of the landscape using LD simulations and experimental information). It is found that the directionality is due to the rate-limiting ADP release step rather than the potential energy of the lever arm angle. We show that the energy of the power stroke and the barriers involved in it are of minor consequence to the selectivity of forward over backward steps and instead suggest that the selective release of ADP from a postrigor myosin motor head promotes highly selective and processive myosin V. Our model is supported by different computational methods-LD simulations, Monte Carlo simulations, and kinetic equations solution-as well as by structure-based binding energy calculations.

Simulating the dynamics of the mechanochemical cycle of myosin-V.

The detailed dynamics of the cycle of myosin-V are explored by simulation approaches, examining the nature of the energy-driven motion. Our study started with Langevin dynamics (LD) simulations on a very coarse landscape with a single rate-limiting barrier and reproduced the stall force and the hand-over-hand dynamics. We then considered a more realistic landscape and used time-dependent Monte Carlo (MC) simulations that allowed trajectories long enough to reproduce the force/velocity characteristic sigmoidal correlation, while also reproducing the hand-over-hand motion. Overall, our study indicated that the notion of a downhill lever-up to lever-down process (popularly known as the powerstroke mechanism) is the result of the energetics of the complete myosin-V cycle and is not the source of directional motion or force generation on its own. The present work further emphasizes the need to use well-defined energy landscapes in studying molecular motors in general and myosin in particular.

Assessment of neonatal platelet adhesion, activation, and aggregation.

BACKGROUND: Acquired and inherited bleeding disorders may present in the neonatal period with devastating lifelong effects. Diagnosing bleeding disorders in the neonatal population could aid in preventing and treating the associated complications. However, currently available platelet function testing is limited in neonates, owing to difficulties in obtaining an adequate blood volume, a lack of normal reference ranges, and an incomplete understanding of the neonatal platelet functional phenotype. OBJECTIVE: To develop small-volume, whole blood platelet function assays in order to quantify and compare neonatal and adult platelet function. METHODS AND RESULTS: Peripheral blood was obtained from healthy, full-term neonates at 24 h of life. Platelet activation, secretion and aggregation were measured via flow cytometry. Platelet adhesion and aggregation were assessed under static and flow conditions. As compared with adult platelets, peripheral neonatal platelet P-selectin expression and integrin glycoprotein IIbIIIa activation were significantly reduced in response to the G-protein-coupled receptor (GPCR) agonists thrombin receptor activator peptide-6 (TRAP-6), ADP, and U46619, and the immunoreceptor tyrosine-based activation motif (ITAM) signaling pathway agonists collagen-related peptide (CRP) and rhodocytin. Neonatal platelet aggregation was markedly reduced in response to TRAP-6, ADP, U46619, CRP and rhodocytin as compared with adult platelets. The extents of neonatal and adult platelet adhesion and aggregate formation under static and shear conditions on collagen and von Willebrand factor were similar. CONCLUSIONS: As compared with adult platelets, we found that neonatal platelet activation and secretion were blunted in response to GPCR or ITAM agonists, whereas the extent of neonatal platelet adhesion and aggregate formation was similar to that of adult platelets.

Unique double-ring structure of the peroxisomal Pex1/Pex6 ATPase complex revealed by cryo-electron microscopy.

Members of the AAA family of ATPases assemble into hexameric double rings and perform vital functions, yet their molecular mechanisms remain poorly understood. Here, we report structures of the Pex1/Pex6 complex; mutations in these proteins frequently cause peroxisomal diseases. The structures were determined in the presence of different nucleotides by cryo-electron microscopy. Models were generated using a computational approach that combines Monte Carlo placement of structurally homologous domains into density maps with energy minimization and refinement protocols. Pex1 and Pex6 alternate in an unprecedented hexameric double ring. Each protein has two N-terminal domains, N1 and N2, structurally related to the single N domains in p97 and N-ethylmaleimide sensitive factor (NSF); N1 of Pex1 is mobile, but the others are packed against the double ring. The N-terminal ATPase domains are inactive, forming a symmetric D1 ring, whereas the C-terminal domains are active, likely in different nucleotide states, and form an asymmetric D2 ring. These results suggest how subunit activity is coordinated and indicate striking similarities between Pex1/Pex6 and p97, supporting the hypothesis that the Pex1/Pex6 complex has a role in peroxisomal protein import analogous to p97 in ER-associated protein degradation.

Detachment, futile cycling, and nucleotide pocket collapse in myosin-V stepping.

Myosin-V is a highly processive dimeric protein that walks with 36-nm steps along actin tracks, powered by coordinated adenosine triphosphate (ATP) hydrolysis reactions in the two myosin heads. No previous theoretical models of the myosin-V walk reproduce all the observed trends of velocity and run length with adenosine diphosphate (ADP), ATP and external forcing. In particular, a result that has eluded all theoretical studies based upon rigorous physical chemistry is that run length decreases with both increasing [ADP] and [ATP]. We systematically analyze which mechanisms in existing models reproduce which experimental trends and use this information to guide the development of models that can reproduce them all. We formulate models as reaction networks between distinct mechanochemical states with energetically determined transition rates. For each network architecture, we compare predictions for velocity and run length to a subset of experimentally measured values, and fit unknown parameters using a bespoke Monte Carlo simulated annealing optimization routine. Finally we determine which experimental trends are replicated by the best-fit model for each architecture. Only two models capture them all: one involving [ADP]-dependent mechanical detachment, and another including [ADP]-dependent futile cycling and nucleotide pocket collapse. Comparing model-predicted and experimentally observed kinetic transition rates favors the latter.

A stochastic model for kinesin bidirectional stepping.

In this paper, a hand-over-hand stochastic model for the dynamics of the conventional kinesin is constructed. In the model, both forward and backward motions are taken into consideration. First passage time distributions, average velocities, dwell times, and forward/backward step ratios are investigated based on the model. A good agreement between the results of the model and experimental data is achieved under a variety of external loads.

Fluctuation-preserving coarse graining for biochemical systems.

Finite stochastic Markov models play a major role in modeling biological systems. Such models are a coarse-grained description of the underlying microscopic dynamics and can be considered mesoscopic. The level of coarse-graining is to a certain extent arbitrary since it depends on the resolution of accommodating measurements. Here we present a systematic way to simplify such stochastic descriptions which preserves both the meso-micro and the meso-macro connections. The former is achieved by demanding locality, the latter by considering cycles on the network of states. Our method preserves fluctuations of observables much better than naïve approaches.

Load and Pi control flux through the branched kinetic cycle of myosin V.

Myosin V is a processive actin-based motor protein that takes multiple 36-nm steps to deliver intracellular cargo to its destination. In the laser trap, applied load slows myosin V heavy meromyosin stepping and increases the probability of backsteps. In the presence of 40 mm phosphate (P(i)), both forward and backward steps become less load-dependent. From these data, we infer that P(i) release commits myosin V to undergo a highly load-dependent transition from a state in which ADP is bound to both heads and its lead head trapped in a pre-powerstroke conformation. Increasing the residence time in this state by applying load increases the probability of backstepping or detachment. The kinetics of detachment indicate that myosin V can detach from actin at two distinct points in the cycle, one of which is turned off by the presence of P(i). We propose a branched kinetic model to explain these data. Our model includes P(i) release prior to the most load-dependent step in the cycle, implying that P(i) release and load both act as checkpoints that control the flux through two parallel pathways.

Perturbation-based Markovian Transmission Model for macromolecular machinery in cell.

The study of the dynamics of a complex system is an important problem that includes large macromolecular complexes, molecular interaction networks, and cell functional modules. Large macromolecular complexes in cellular machinery can be modeled as a connected network, as in the elastic or Gaussian network models as demonstrated by Bahar and colleagues. Here we propose the Perturbation-based Markovian Transmission Model for studying the dynamics of signal transmission in macromolecular machinery. The initial perturbation is transmitted by a Markovian processes, and the dynamics of the probability flow is analytically solved using the master equation. Due to the large size of macromolecular complexes, it is very difficult to obtain analytical time-dependent Markovian dynamics of all atoms from the first perturbation until stationary state. To overcome it, we decrease the level of complexity of the transition matrix using a Krylov subspace method. This method is equivalent to integrating all eigen modes, and we show it can provide a globally accurate solution to the dynamics problem of signal transmission for very large macromolecular complexes with reasonable computational time. We give results of the dynamics of the GroEL-GroES chaperone system by applying uniform perturbation to all residues. We are able to identify experimentally found important residues and provide a set of predicted pivot, messenger, and effector residues, each with distinct dynamic behavior. Further results of selective perturbation on the surface of ATP binding pocket identifies the path of maximal probability flow of signal.

Did protein kinase regulatory mechanisms evolve through elaboration of a simple structural component?

Statistical analysis of the functional constraints acting on eukaryotic protein kinases (EPKs) and on distantly related kinases suggests that EPK regulatory mechanisms evolved around an ancient structural component whose most distinctive features include the HxD-motif adjoining the catalytic loop, the F-helix, an F-helix aspartate, and the DFG-motif adjoined to the activation loop. The HxD-histidine constitutes a convergence point for signal integration, as conserved interactions link it to key catalytic residues, to the F-helix aspartate, and to both ends of the DFG-motif. These and other conserved features appear to be associated with DFG conformational changes and with coordinated movements possibly associated with phosphate transfer and ADP release. The EPKs have acquired structural features that link this core component to likely substrate-interacting regions at either end of the F-helix (most notably involving an F-helix tryptophan) and to three regions undergoing conformational changes upon kinase activation: the activation segment, the C-helix, and the nucleotide-binding pocket.

Disparity in allosteric interactions of monastrol with Eg5 in the presence of ADP and ATP: a difference FT-IR investigation.

Eg5 is a kinesin-like motor protein required for mitotic progression in higher eukaryotes. It is thought to cross-link antiparallel microtubules, and provides a force required for the formation of a bipolar spindle. Monastrol causes the catastrophic collapse of the mitotic spindle through the allosteric inhibition of Eg5. Utilizing a truncated Eg5 protein, we employ difference infrared spectroscopy to probe structural changes that occur in the motor protein with monastrol in the presence of either ADP or ATP. Difference FT-IR spectra of Eg5-monastrol-nucleotide complexes demonstrate that there are triggered conformational changes corresponding to an interconversion of secondary structural elements in the motor upon interaction with nucleotides. Notably, conformational changes elicited in the presence of ADP are different from those in the presence of ATP. In Eg5-monastrol complexes, exchange of ADP is associated with a decrease in random structure and an increase in alpha-helical content. In contrast, formation of the Eg5-monastrol-ATP complex is associated with a decrease in alpha-helical content and a concomitant increase in beta-sheet content. One or more carboxylic acid residues in Eg5 undergo unique changes when ATP, but not ADP, interacts with the motor domain in the presence of monastrol. This first direct dissection of inhibitor-protein interactions, using these methods, demonstrates a clear disparity in the structural consequences of monastrol in the presence of ADP versus ATP.