Light chain 2 is a Tctex-type related axonemal dynein light chain that regulates directional ciliary motility in Trypanosoma brucei.

Flagellar motility is essential for the cell morphology, viability, and virulence of pathogenic kinetoplastids. Trypanosoma brucei flagella beat with a bending wave that propagates from the flagellum's tip to its base, rather than base-to-tip as in other eukaryotes. Thousands of dynein motor proteins coordinate their activity to drive ciliary bending wave propagation. Dynein-associated light and intermediate chains regulate the biophysical mechanisms of axonemal dynein. Tctex-type outer arm dynein light chain 2 (LC2) regulates flagellar bending wave propagation direction, amplitude, and frequency in Chlamydomonas reinhardtii. However, the role of Tctex-type light chains in regulating T. brucei motility is unknown. Here, we used a combination of bioinformatics, in-situ molecular tagging, and immunofluorescence microscopy to identify a Tctex-type light chain in the procyclic form of T. brucei (TbLC2). We knocked down TbLC2 expression using RNAi in both wild-type and FLAM3, a flagellar attachment zone protein, knockdown cells and quantified TbLC2's effects on trypanosome cell biology and biophysics. We found that TbLC2 knockdown reduced the directional persistence of trypanosome cell swimming, induced an asymmetric ciliary bending waveform, modulated the bias between the base-to-tip and tip-to-base beating modes, and increased the beating frequency. Together, our findings are consistent with a model of TbLC2 as a down-regulator of axonemal dynein activity that stabilizes the forward tip-to-base beating ciliary waveform characteristic of trypanosome cells. Our work sheds light on axonemal dynein regulation mechanisms that contribute to pathogenic kinetoplastids' unique tip-to-base ciliary beating nature and how those mechanisms underlie dynein-driven ciliary motility more generally.

Long-range electrostatic interactions significantly modulate the affinity of dynein for microtubules.

The dynein family of microtubule minus-end-directed motor proteins drives diverse functions in eukaryotic cells, including cell division, intracellular transport, and flagellar beating. Motor protein processivity, which characterizes how far a motor walks before detaching from its filament, depends on the interaction between its microtubule-binding domain (MTBD) and the microtubule. Dynein's MTBD switches between high- and low-binding affinity states as it steps. Significant structural and functional data show that specific salt bridges within the MTBD and between the MTBD and the microtubule govern these affinity state shifts. However, recent computational work suggests that nonspecific, long-range electrostatic interactions between the MTBD and the microtubule may also play an important role in the processivity of dynein. To investigate this hypothesis, we mutated negatively charged amino acids remote from the dynein MTBD-microtubule-binding interface to neutral residues and measured the binding affinity using microscale thermophoresis and optical tweezers. We found a significant increase in the binding affinity of the mutated MTBDs for microtubules. Furthermore, we found that charge screening by free ions in solution differentially affected the binding and unbinding rates of MTBDs to microtubules. Together, these results demonstrate a significant role for long-range electrostatic interactions in regulating dynein-microtubule affinity. Moreover, these results provide insight into the principles that potentially underlie the biophysical differences between molecular motors with various processivities and protein-protein interactions more generally.

Calculating the force-dependent unbinding rate of biological macromolecular bonds from force-ramp optical trapping assays.

The non-covalent biological bonds that constitute protein-protein or protein-ligand interactions play crucial roles in many cellular functions, including mitosis, motility, and cell-cell adhesion. The effect of external force ([Formula: see text]) on the unbinding rate ([Formula: see text]) of macromolecular interactions is a crucial parameter to understanding the mechanisms behind these functions. Optical tweezer-based single-molecule force spectroscopy is frequently used to obtain quantitative force-dependent dissociation data on slip, catch, and ideal bonds. However, analyses of this data using dissociation time or dissociation force histograms often quantitatively compare bonds without fully characterizing their underlying biophysical properties. Additionally, the results of histogram-based analyses can depend on the rate at which force was applied during the experiment and the experiment's sensitivity. Here, we present an analytically derived cumulative distribution function-like approach to analyzing force-dependent dissociation force spectroscopy data. We demonstrate the benefits and limitations of the technique using stochastic simulations of various bond types. We show that it can be used to obtain the detachment rate and force sensitivity of biological macromolecular bonds from force spectroscopy experiments by explicitly accounting for loading rate and noisy data. We also discuss the implications of our results on using optical tweezers to collect force-dependent dissociation data.

Out-of-Equilibrium Biophysical Chemistry: The Case for Multidimensional, Integrated Single-Molecule Approaches.

Out-of-equilibrium processes are ubiquitous across living organisms and all structural hierarchies of life. At the molecular scale, out-of-equilibrium processes (for example, enzyme catalysis, gene regulation, and motor protein functions) cause biological macromolecules to sample an ensemble of conformations over a wide range of time scales. Quantifying and conceptualizing the structure-dynamics to function relationship is challenging because continuously evolving multidimensional energy landscapes are necessary to describe nonequilibrium biological processes in biological macromolecules. In this perspective, we explore the challenges associated with state-of-the-art experimental techniques to understanding biological macromolecular function. We argue that it is time to revisit how we probe and model functional out-of-equilibrium biomolecular dynamics. We suggest that developing integrated single-molecule multiparametric force-fluorescence instruments and using advanced molecular dynamics simulations to study out-of-equilibrium biomolecules will provide a path towards understanding the principles of and mechanisms behind the structure-dynamics to function paradigm in biological macromolecules.

Ensemble Switching Unveils a Kinetic Rheostat Mechanism of the Eukaryotic Thiamine Pyrophosphate Riboswitch.

Thiamine pyrophosphate (TPP) riboswitches regulate thiamine metabolism by inhibiting the translation of enzymes essential to thiamine synthesis pathways upon binding to thiamine pyrophosphate in cells across all domains of life. Recent work on the Arabidopsis thaliana TPP riboswitch suggests a multi-step TPP binding process involving multiple riboswitch configurational ensembles and that Mg2+ dependence underlies the mechanism of TPP recognition and subsequent transition to the expression-inhibiting state of the aptamer domain followed by changes in the expression platform. However, details of the relationship between TPP riboswitch conformational changes and interactions with TPP and Mg2+ ¬¬in the aptamer domain constituting this mechanism are unknown. Therefore, we integrated single-molecule multiparameter fluorescence and force spectroscopy with atomistic molecular dynamics simulations and found that conformational transitions within the aptamer domain's sensor helices associated with TPP and Mg2+ ligand binding occurred between at least five different ensembles on timescales ranging from µs to ms. These dynamics are orders of magnitude faster than the 10 second-timescale folding kinetics associated with expression-state switching in the switch sequence. Together, our results show that a TPP and Mg2+ dependent mechanism determines dynamic configurational state ensemble switching of the aptamer domain's sensor helices that regulates the stability of the switch helix, which ultimately may lead to the expression-inhibiting state of the riboswitch. Additionally, we propose that two pathways exist for ligand recognition and that this mechanism underlies a kinetic rheostat-like behavior of the Arabidopsis thaliana TPP riboswitch.

Reporting on the future of integrative structural biology ORAU workshop.

Integrative and hybrid methods have the potential to bridge long-standing knowledge gaps in structural biology. These methods will have a prominent role in the future of the field as we make advances toward a complete, unified representation of biology that spans the molecular and cellular scales. The Department of Physics and Astronomy at Clemson University hosted The Future of Integrative Structural Biology workshop on April 29, 2017 and partially sponsored by partially sponsored by a program of the Oak Ridge Associated Universities (ORAU). The workshop brought experts from multiple structural biology disciplines together to discuss near-term steps toward the goal of a molecular atlas of the cell. The discussion focused on the types of structural data that should be represented, how this data should be represented, and how the time domain might be incorporated into such an atlas. The consensus was that an explorable, map-like Virtual Cell, containing both spatial and temporal data bridging the atomic and cellular length scales obtained by multiple experimental methods, represents the best path toward a complete atlas of the cell.

From isolated structures to continuous networks: A categorization of cytoskeleton-based motile engineered biological microstructures.

As technology at the small scale is advancing, motile engineered microstructures are becoming useful in drug delivery, biomedicine, and lab-on-a-chip devices. However, traditional engineering methods and materials can be inefficient or functionally inadequate for small-scale applications. Increasingly, researchers are turning to the biology of the cytoskeleton, including microtubules, actin filaments, kinesins, dyneins, myosins, and associated proteins, for both inspiration and solutions. They are engineering structures with components that range from being entirely biological to being entirely synthetic mimics of biology and on scales that range from isotropic continuous networks to single isolated structures. Motile biological microstructures trace their origins from the development of assays used to study the cytoskeleton to the array of structures currently available today. We define 12 types of motile biological microstructures, based on four categories: entirely biological, modular, hybrid, and synthetic, and three scales: networks, clusters, and isolated structures. We highlight some key examples, the unique functionalities, and the potential applications of each microstructure type, and we summarize the quantitative models that enable engineering them. By categorizing the diversity of motile biological microstructures in this way, we aim to establish a framework to classify these structures, define the gaps in current research, and spur ideas to fill those gaps. This article is categorized under: Nanotechnology Approaches to Biology > Nanoscale Systems in Biology Nanotechnology Approaches to Biology > Cells at the Nanoscale Biology-Inspired Nanomaterials > Protein and Virus-Based Structures Therapeutic Approaches and Drug Discovery > Emerging Technologies.

Resistance is futile: Centering forces yield for asymmetric cell division.

Asymmetric cell division relies on microtubule-based forces to asymmetrically position the mitotic apparatus. In this issue, Sallé et al. (2019. J. Cell Biol. https://doi.org/10.1083/jcb.201807102) use magnetic tweezers to induce asymmetric division in sea urchin zygotes, demonstrating that asymmetry could arise from a time-dependent weakening of centering forces.

E-hooks provide guidance and a soft landing for the microtubule binding domain of dynein.

Macromolecular binding is a complex process that involves sensing and approaching the binding partner, adopting the proper orientation, and performing the physical binding. We computationally investigated the role of E-hooks, which are intrinsically disordered regions (IDRs) at the C-terminus of tubulin, on dynein microtubule binding domain (MTBD) binding to the microtubule as a function of the distance between the MTBD and its binding site on the microtubule. Our results demonstrated that the contacts between E-hooks and the MTBD are dynamical; multiple negatively charted patches of amino acids on the E-hooks grab and release the same positively charged patches on the MTBD as it approaches the microtubule. Even when the distance between the MTBD and the microtubule was greater than the E-hook length, the E-hooks sensed and guided MTBD via long-range electrostatic interactions in our simulations. Moreover, we found that E-hooks exerted electrostatic forces on the MTBD that were distance dependent; the force pulls the MTBD toward the microtubule at long distances but opposes binding at short distances. This mechanism provides a "soft-landing" for the MTBD as it binds to the microtubule. Finally, our analysis of the conformational states of E-hooks in presence and absence of the MTBD indicates that the binding process is a mixture of the induced-fit and lock-and-key macromolecular binding hypotheses. Overall, this novel binding mechanism is termed "guided-soft-binding" and could have broad-reaching impacts on the understanding of how IDRs dock to structured proteins.

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.

Mechanism of microtubule lumen entry for the α-tubulin acetyltransferase enzyme αTAT1.

Microtubules are structural polymers inside of cells that are subject to posttranslational modifications. These posttranslational modifications create functionally distinct subsets of microtubule networks in the cell, and acetylation is the only modification that takes place in the hollow lumen of the microtubule. Although it is known that the α-tubulin acetyltransferase (αTAT1) is the primary enzyme responsible for microtubule acetylation, the mechanism for how αTAT1 enters the microtubule lumen to access its acetylation sites is not well understood. By performing biochemical assays, fluorescence and electron microscopy experiments, and computational simulations, we found that αTAT1 enters the microtubule lumen through the microtubule ends, and through bends or breaks in the lattice. Thus, microtubule structure is an important determinant in the acetylation process. In addition, once αTAT1 enters the microtubule lumen, the mobility of αTAT1 within the lumen is controlled by the affinity of αTAT1 for its acetylation sites, due to the rapid rebinding of αTAT1 onto highly concentrated α-tubulin acetylation sites. These results have important implications for how acetylation could gradually accumulate on stable subsets of microtubules inside of the cell.

Cytoplasmic dynein binding, run length, and velocity are guided by long-range electrostatic interactions.

Dyneins are important molecular motors involved in many essential biological processes, including cargo transport along microtubules, mitosis, and in cilia. Dynein motility involves the coupling of microtubule binding and unbinding to a change in the configuration of the linker domain induced by ATP hydrolysis, which occur some 25 nm apart. This leaves the accuracy of dynein stepping relatively inaccurate and susceptible to thermal noise. Using multi-scale modeling with a computational focusing technique, we demonstrate that the microtubule forms an electrostatic funnel that guides the dynein's microtubule binding domain (MTBD) as it finally docks to the precise, keyed binding location on the microtubule. Furthermore, we demonstrate that electrostatic component of the MTBD's binding free energy is linearly correlated with the velocity and run length of dynein, and we use this linearity to predict the effect of mutating each glutamic and aspartic acid located in MTBD domain to alanine. Lastly, we show that the binding of dynein to the microtubule is associated with conformational changes involving several helices, and we localize flexible hinge points within the stalk helices. Taken all together, we demonstrate that long range electrostatic interactions bring a level of precision to an otherwise noisy dynein stepping process.

Multiscale method for modeling binding phenomena involving large objects: application to kinesin motor domains motion along microtubules.

Many biological phenomena involve the binding of proteins to a large object. Because the electrostatic forces that guide binding act over large distances, truncating the size of the system to facilitate computational modeling frequently yields inaccurate results. Our multiscale approach implements a computational focusing method that permits computation of large systems without truncating the electrostatic potential and achieves the high resolution required for modeling macromolecular interactions, all while keeping the computational time reasonable. We tested our approach on the motility of various kinesin motor domains. We found that electrostatics help guide kinesins as they walk: N-kinesins towards the plus-end, and C-kinesins towards the minus-end of microtubules. Our methodology enables computation in similar, large systems including protein binding to DNA, viruses, and membranes.

The motility of axonemal dynein is regulated by the tubulin code.

Microtubule diversity, arising from the utilization of different tubulin genes and from posttranslational modifications, regulates many cellular processes including cell division, neuronal differentiation and growth, and centriole assembly. In the case of cilia and flagella, multiple cell biological studies show that microtubule diversity is important for axonemal assembly and motility. However, it is not known whether microtubule diversity directly influences the activity of the axonemal dyneins, the motors that drive the beating of the axoneme, nor whether the effects on motility are indirect, perhaps through regulatory pathways upstream of the motors, such as the central pair, radial spokes, or dynein regulatory complex. To test whether microtubule diversity can directly regulate the activity of axonemal dyneins, we asked whether in vitro acetylation or deacetylation of lysine 40 (K40), a major posttranslational modification of α-tubulin, or whether proteolytic cleavage of the C-terminal tail (CTT) of α- and ß-tubulin, the location of detyrosination, polyglutamylation, and polyglycylation modifications as well as most of the genetic diversity, can influence the activity of outer arm axonemal dynein in motility assays using purified proteins. By quantifying the motility with displacement-weighted velocity analysis and mathematically modeling the results, we found that K40 acetylation increases and CTTs decrease axonemal dynein motility. These results show that axonemal dynein directly deciphers the tubulin code, which has important implications for eukaryotic ciliary beat regulation.

Displacement-weighted velocity analysis of gliding assays reveals that Chlamydomonas axonemal dynein preferentially moves conspecific microtubules.

In vitro gliding assays, in which microtubules are observed to glide over surfaces coated with motor proteins, are important tools for studying the biophysics of motility. Gliding assays with axonemal dyneins have the unusual feature that the microtubules exhibit large variations in gliding speed despite measures taken to eliminate unsteadiness. Because axonemal dynein gliding assays are usually done using heterologous proteins, i.e., dynein and tubulin from different organisms, we asked whether the source of tubulin could underlie the unsteadiness. By comparing gliding assays with microtubules polymerized from Chlamydomonas axonemal tubulin with those from porcine brain tubulin, we found that the unsteadiness is present despite matching the source of tubulin to the source of dynein. We developed a novel, to our knowledge, displacement-weighted velocity analysis to quantify both the velocity and the unsteadiness of gliding assays systematically and without introducing bias toward low motility. We found that the quantified unsteadiness is independent of tubulin source. In addition, we found that the short Chlamydomonas microtubules translocate significantly faster than their porcine counterparts. By modeling the effect of length on velocity, we propose that the observed effect may be due to a higher rate of binding of Chlamydomonas axonemal dynein to Chlamydomonas microtubules than to porcine microtubules.

Reconstitution of flagellar sliding.

The motile structure within eukaryotic cilia and flagella is the axoneme. This structure typically consists of nine doublet microtubules arranged around a pair of singlet microtubules. The axoneme contains more than 650 different proteins that have structural, force-generating, and regulatory functions. Early studies on sea urchin sperm identified the force-generating components, the dynein motors. It was shown that dynein can slide adjacent doublet microtubules in the presence of ATP. How this sliding gives rise to the beating of the axoneme is still unknown. Reconstitution assays provide a clean system, free from cellular effects, to elucidate the underlying beating mechanisms. These assays can be used to identify the components that are both necessary and sufficient for the generation of flagellar beating.

One-step purification of assembly-competent tubulin from diverse eukaryotic sources.

We have developed a protocol that allows rapid and efficient purification of native, active tubulin from a variety of species and tissue sources by affinity chromatography. The affinity matrix comprises a bacterially expressed, recombinant protein, the TOG1/2 domains from Saccharomyces cerevisiae Stu2, covalently coupled to a Sepharose support. The resin has a high capacity to specifically bind tubulin from clarified crude cell extracts, and, after washing, highly purified tubulin can be eluted under mild conditions. The eluted tubulin is fully functional and can be efficiently assembled into microtubules. The method eliminates the need to use heterologous systems for the study of microtubule-associated proteins and motor proteins, which has been a major issue in microtubule-related research.

Effect of ligands on thermal dissipation from gold nanorods.

Thermal interface conductance was measured for soluble gold nanorods (NRs) coated with mercaptocarboxylic acids (HS-(CH(2))(n)COOH, n = 5, 10, 15), thiolated polyethylene glycols (MW = 356, 1000, 5000), and HS-(CH(2))(15)-COOH-coated NRs further coated with alternating layers of poly(diallyldimethylammonium chloride) and poly(sodium styrenesulfonate). Ferguson analysis determined ligand thickness. The thermal-diffusion-dominated regime of transient absorption spectra was fit to a continuum heat diffusion finite element model to obtain the thermal interface conductance, G, which varied with ligand chemistry but not molecule length. The results suggest that the ability to exclude water from the NR surface governs ligand G values.