Effect of α-tubulin acetylation on the doublet microtubule structure.

Acetylation of α-tubulin at the lysine 40 residue (αK40) by αTAT1/MEC-17 acetyltransferase modulates microtubule properties and occurs in most eukaryotic cells. Previous literatures suggest that acetylated microtubules are more stable and damage resistant. αK40 acetylation is the only known microtubule luminal post-translational modification site. The luminal location suggests that the modification tunes the lateral interaction of protofilaments inside the microtubule. In this study, we examined the effect of tubulin acetylation on the doublet microtubule (DMT) in the cilia of Tetrahymena thermophila using a combination of cryo-electron microscopy, molecular dynamics, and mass spectrometry. We found that αK40 acetylation exerts a small-scale effect on the DMT structure and stability by influencing the lateral rotational angle. In addition, comparative mass spectrometry revealed a link between αK40 acetylation and phosphorylation in cilia.

The ATPase asymmetry: Novel computational insight into coupling diverse FO motors with tripartite F1.

ATP synthase, a crucial enzyme for cellular bioenergetics, operates via the coordinated coupling of an FO motor, which presents variable symmetry, and a tripartite F1 motor. Despite extensive research, the understanding of their coupling dynamics, especially with non-10-fold symmetrical FO motors, remains incomplete. This study investigates the coupling patterns between eightfold and ninefold FO motors and the constant threefold F1 motor using coarse-grained molecular dynamics simulations. We unveil that in the case of a ninefold FO motor, a 3-3-3 motion is most likely to occur, whereas a 3-3-2 motion predominates with an eightfold FO motor. Furthermore, our findings propose a revised model for the coupling method, elucidating that the pathways' energy usage is primarily influenced by F1 rotation and conformational changes hindered by the b-subunits. Our results present a crucial step toward comprehending the energy landscape and mechanisms governing ATP synthase operation.

Removing the parachuting artifact using two-way scanning data in high-speed atomic force microscopy.

The high-speed atomic force microscopy (HS-AFM) is a unique and prominent method to observe structural dynamics of biomolecules at single molecule level at near-physiological condition. To achieve high temporal resolution, the probe tip scans the stage at high speed which can cause the so-called parachuting artifact in the HS-AFM images. Here, we develop a computational method to detect and remove the parachuting artifact in HS-AFM images using the two-way scanning data. To merge the two-way scanning images, we employed a method to infer the piezo hysteresis effect and to align the forward- and backward-scanning images. We then tested our method for HS-AFM videos of actin filaments, molecular chaperone, and duplex DNA. Together, our method can remove the parachuting artifact from the raw HS-AFM video containing two-way scanning data and make the processed video free from the parachuting artifact. The method is general and fast so that it can easily be applied to any HS-AFM videos with two-way scanning data.

Native doublet microtubules from Tetrahymena thermophila reveal the importance of outer junction proteins.

Cilia are ubiquitous eukaryotic organelles responsible for cellular motility and sensory functions. The ciliary axoneme is a microtubule-based cytoskeleton consisting of two central singlets and nine outer doublet microtubules. Cryo-electron microscopy-based studies have revealed a complex network inside the lumen of both tubules composed of microtubule-inner proteins (MIPs). However, the functions of most MIPs remain unknown. Here, we present single-particle cryo-EM-based analyses of the Tetrahymena thermophila native doublet microtubule and identify 42 MIPs. These data shed light on the evolutionarily conserved and diversified roles of MIPs. In addition, we identified MIPs potentially responsible for the assembly and stability of the doublet outer junction. Knockout of the evolutionarily conserved outer junction component CFAP77 moderately diminishes Tetrahymena swimming speed and beat frequency, indicating the important role of CFAP77 and outer junction stability in cilia beating generation and/or regulation.

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.

FO-F1 coupling and symmetry mismatch in ATP synthase resolved in every FO rotation step.

FOF1 ATP synthase, a ubiquitous enzyme that synthesizes most ATP in living cells, is composed of two rotary motors: a membrane-embedded proton-driven FO motor and a catalytic F1 motor. These motors share both central and peripheral stalks. Although both FO and F1 have pseudo-symmetric structures, their symmetries do not match. How symmetry mismatch is solved remains elusive because of the missing intermediate structures of the rotational steps. Here, for the case of Bacillus PS3 ATP synthases with three- and 10-fold symmetries in F1 and FO, respectively, we uncovered the mechanical couplings between FO and F1 at every 36° rotation step via molecular dynamics simulations and comparative studies of cryoelectron microscopy (cryo-EM) structures from three species. We found that the mismatch could be solved using several elements: 1) the F1 head partially rotates relative to the FO a subunit via elastic distortion of the b subunits, 2) the rotor is twisted, and 3) comparisons of cryo-EM structures further suggest that the c ring rotary angles can deviate from the symmetric ones. In addition, the F1 motor may have non-canonical structures, relieving stronger frustration. Thus, we provide new insights for solving the symmetry mismatch problem.

Rotational Mechanism of FO Motor in the F-Type ATP Synthase Driven by the Proton Motive Force.

In FOF1 ATP synthase, driven by the proton motive force across the membrane, the FO motor rotates the central rotor and induces conformational changes in the F1 motor, resulting in ATP synthesis. Recently, many near-atomic resolution structural models have been obtained using cryo-electron microscopy. Despite high resolution, however, static information alone cannot elucidate how and where the protons pass through the FO and how proton passage is coupled to FO rotation. Here, we review theoretical and computational studies based on FO structure models. All-atom molecular dynamics (MD) simulations elucidated changes in the protonation/deprotonation of glutamate-the protein-carrier residue-during rotation and revealed the protonation states that form the "water wire" required for long-range proton hopping. Coarse-grained MD simulations unveiled a free energy surface based on the protonation state and rotational angle of the rotor. Hybrid Monte Carlo and MD simulations showed how proton transfer is coupled to rotation.

Cooperation among c-subunits of FoF1-ATP synthase in rotation-coupled proton translocation.

In FoF1-ATP synthase, proton translocation through Fo drives rotation of the c-subunit oligomeric ring relative to the a-subunit. Recent studies suggest that in each step of the rotation, key glutamic acid residues in different c-subunits contribute to proton release to and proton uptake from the a-subunit. However, no studies have demonstrated cooperativity among c-subunits toward FoF1-ATP synthase activity. Here, we addressed this using Bacillus PS3 ATP synthase harboring a c-ring with various combinations of wild-type and cE56D, enabled by genetically fused single-chain c-ring. ATP synthesis and proton pump activities were decreased by a single cE56D mutation and further decreased by double cE56D mutations. Moreover, activity further decreased as the two mutation sites were separated, indicating cooperation among c-subunits. Similar results were obtained for proton transfer-coupled molecular simulations. The simulations revealed that prolonged proton uptake in mutated c-subunits is shared between two c-subunits, explaining the cooperation observed in biochemical assays.

Cells need to be able to store and transfer energy to fuel their various activities. To do this, they produce a small molecule called ATP to carry the energy, which is then released when the ATP is broken down. An enzyme found in plants, animals and bacteria, called F_oF₁ ATP synthase, can both create and use ATP. When it does this, protons, or positive hydrogen ions, are transported across cellular boundaries called membranes. The region of the enzyme that is responsible for pumping the protons contains different parts known as the c-ring and the a-subunit. The movement of protons drives the c-ring to rotate relative to the a-subunit, which leads to producing ATP. Previous research using simulations and the protein structures found there are two or three neighbouring amino acids in the c-ring that face the a-subunit, suggesting that these amino acids act together to drive the rotation. To test this hypothesis, Mitome et al. mutated these amino acids to examine the effect on the enzyme's ability to produce ATP. A single mutation reduced the production of ATP, which decreased even further with mutations in two of the amino acids. The extent of this decrease depended on the distance between the two mutations in the c-ring. Simulations of these changes also found similar results. This indicates there is coordination between different parts of the c-ring to increase the rate of ATP production. This study offers new insights into the molecular processes controlling ATP synthesis and confirms previous theoretical research. This will interest specialists in bioenergetics because it addresses a fundamental biological question with broad impact.

Heterogeneous dissociation process of truncated RNAs by oligomerized Vasa helicase.

RNA helicases are enzymes that generally unwind double-stranded RNA using ATP hydrolysis energy, mainly involved in RNA metabolism, transcription, translation, and mRNA splicing. While the helicase core is crucial for RNA unwinding activity, N- and C-terminal extensions of specific helicases may contain an intrinsically disordered region for electrostatic interaction, resulting in the formation of droplets in the cytoplasm. However, how the disordered region of the RNA helicase contributes to RNA unwinding and dissociation remains unclear. Here, we focused on Bombyx mori Vasa, which unwinds truncated target transposon RNAs from the piRNA-induced silencing complex piRISC. In this study, we used single-molecule techniques to visualise how Vasa dynamically interacts with piRISC and investigate how Vasa oligomerization is involved in the process of piRNA amplification, named the ping-pong pathway. We found that Vasa's oligomerization is required during these processes in vitro and in vivo, and that Vasa triggers the dissociation of truncated RNA in heterogeneous pathways. Our single-molecule results suggest that oligomerized Vasa guides the timing of the process regulating overall dissociation efficiency.

Remodeling and activation mechanisms of outer arm dyneins revealed by cryo-EM.

Cilia are thin microtubule-based protrusions of eukaryotic cells. The swimming of ciliated protists and sperm cells is propelled by the beating of cilia. Cilia propagate the flow of mucus in the trachea and protect the human body from viral infections. The main force generators of ciliary beating are the outer dynein arms (ODAs) which attach to the doublet microtubules. The bending of cilia is driven by the ODAs' conformational changes caused by ATP hydrolysis. Here, we report the native ODA complex structure attaching to the doublet microtubule by cryo-electron microscopy. The structure reveals how the ODA complex is attached to the doublet microtubule via the docking complex in its native state. Combined with coarse-grained molecular dynamic simulations, we present a model of how the attachment of the ODA to the doublet microtubule induces remodeling and activation of the ODA complex.

Role of bacterial RNA polymerase gate opening dynamics in DNA loading and antibiotics inhibition elucidated by quasi-Markov State Model.

To initiate transcription, the holoenzyme (RNA polymerase [RNAP] in complex with σ factor) loads the promoter DNA via the flexible loading gate created by the clamp and ß-lobe, yet their roles in DNA loading have not been characterized. We used a quasi-Markov State Model (qMSM) built from extensive molecular dynamics simulations to elucidate the dynamics of Thermus aquaticus holoenzyme's gate opening. We showed that during gate opening, ß-lobe oscillates four orders of magnitude faster than the clamp, whose opening depends on the Switch 2's structure. Myxopyronin, an antibiotic that binds to Switch 2, was shown to undergo a conformational selection mechanism to inhibit clamp opening. Importantly, we reveal a critical but undiscovered role of ß-lobe, whose opening is sufficient for DNA loading even when the clamp is partially closed. These findings open the opportunity for the development of antibiotics targeting ß-lobe of RNAP. Finally, we have shown that our qMSMs, which encode non-Markovian dynamics based on the generalized master equation formalism, hold great potential to be widely applied to study biomolecular dynamics.

Resolving the data asynchronicity in high-speed atomic force microscopy measurement via the Kalman Smoother.

High-speed atomic force microscopy (HS-AFM) is a scanning probe microscopy that can capture structural dynamics of biomolecules in real time at single molecule level near physiological condition. Albeit much improvement, while scanning one frame of HS-AFM movies, biomolecules often change their conformations largely. Thus, the obtained frame images can be hampered by the time-difference, the asynchronicity, in the data acquisition. Here, to resolve this data asynchronicity in the HS-AFM movie, we developed Kalman filter and smoother methods, some of the sequential Bayesian filtering approaches. The Kalman filter/smoother methods use alternative steps of a short time-propagation by a linear dynamical system and a correction by the likelihood of AFM data acquired pixel by pixel. We first tested the method using a toy model of a diffusing cone, showing that the Kalman smoother method outperforms to reproduce the ground-truth movie. We then applied the Kalman smoother to a synthetic movie for conformational change dynamics of a motor protein, i.e., dynein, confirming the superiority of the Kalman smoother. Finally, we applied the Kalman smoother to two real HS-AFM movies, FlhAC and centralspindlin, reducing distortion and noise in the AFM movies. The method is general and can be applied to any HS-AFM movies.

Molecular dynamics simulation of proton-transfer coupled rotations in ATP synthase FO motor.

The FO motor in FOF1 ATP synthase rotates its rotor driven by the proton motive force. While earlier studies elucidated basic mechanisms therein, recent advances in high-resolution cryo-electron microscopy enabled to investigate proton-transfer coupled FO rotary dynamics at structural details. Here, taking a hybrid Monte Carlo/molecular dynamics simulation method, we studied reversible dynamics of a yeast mitochondrial FO. We obtained the 36°-stepwise rotations of FO per one proton transfer in the ATP synthesis mode and the proton pumping in the ATP hydrolysis mode. In both modes, the most prominent path alternatively sampled states with two and three deprotonated glutamates in c-ring, by which the c-ring rotates one step. The free energy transduction efficiency in the model FO motor reached ~ 90% in optimal conditions. Moreover, mutations in key glutamate and a highly conserved arginine increased proton leakage and markedly decreased the coupling, in harmony with previous experiments. This study provides a simple framework of simulations for chemical-reaction coupled molecular dynamics calling for further studies in ATP synthase and others.

How Cytoplasmic Dynein Couples ATP Hydrolysis Cycle to Diverse Stepping Motions: Kinetic Modeling.

Cytoplasmic dynein is a two-headed molecular motor that moves to the minus end of a microtubule by ATP hydrolysis free energy. By employing its two heads (motor domains), cytoplasmic dynein exhibits various bipedal stepping motions: inchworm and hand-over-hand motions, as well as nonalternating steps of one head. However, the molecular basis to achieve such diverse stepping manners remains unclear because of the lack of an experimental method to observe stepping and the ATPase reaction of dynein simultaneously. Here, we propose a kinetic model for bipedal motions of cytoplasmic dynein and perform Gillespie Monte Carlo simulations that qualitatively reproduce most experimental data obtained to date. The model represents the status of each motor domain as five states according to conformation and nucleotide- and microtubule-binding conditions of the domain. In addition, the relative positions of the two domains were approximated by three discrete states. Accompanied by ATP hydrolysis cycles, the model dynein stochastically and processively moved forward in multiple steps via diverse pathways, including inchworm and hand-over-hand motions, similarly to experimental data. The model reproduced key experimental motility-related properties, including velocity and run length, as functions of the ATP concentration and external force, therefore providing a plausible explanation of how dynein achieves various stepping manners with explicit characterization of nucleotide states. Our model highlights the uniqueness of dynein in the coupling of ATPase with its movement during both inchworm and hand-over-hand stepping.

The inner junction complex of the cilia is an interaction hub that involves tubulin post-translational modifications.

Microtubules are cytoskeletal structures involved in stability, transport and organization in the cell. The building blocks, the α- and ß-tubulin heterodimers, form protofilaments that associate laterally into the hollow microtubule. Microtubule also exists as highly stable doublet microtubules in the cilia where stability is needed for ciliary beating and function. The doublet microtubule maintains its stability through interactions at its inner and outer junctions where its A- and B-tubules meet. Here, using cryo-electron microscopy, bioinformatics and mass spectrometry of the doublets of Chlamydomonas reinhardtii and Tetrahymena thermophila, we identified two new inner junction proteins, FAP276 and FAP106, and an inner junction-associated protein, FAP126, thus presenting the complete answer to the inner junction identity and localization. Our structural study of the doublets shows that the inner junction serves as an interaction hub that involves tubulin post-translational modifications. These interactions contribute to the stability of the doublet and hence, normal ciliary motility.

Tubulin lattice in cilia is in a stressed form regulated by microtubule inner proteins.

Cilia, the hair-like protrusions that beat at high frequencies to propel a cell or move fluid around are composed of radially bundled doublet microtubules. In this study, we present a near-atomic resolution map of the Tetrahymena doublet microtubule by cryoelectron microscopy. The map demonstrates that the network of microtubule inner proteins weaves into the tubulin lattice and forms an inner sheath. From mass spectrometry data and de novo modeling, we identified Rib43a proteins as the filamentous microtubule inner proteins in the protofilament ribbon region. The Rib43a-tubulin interaction leads to an elongated tubulin dimer distance every 2 dimers. In addition, the tubulin lattice structure with missing microtubule inner proteins (MIPs) by sarkosyl treatment shows significant longitudinal compaction and lateral angle change between protofilaments. These results are evidence that the MIPs directly affect and stabilize the tubulin lattice. It suggests that the doublet microtubule is an intrinsically stressed filament and that this stress could be manipulated in the regulation of ciliary waveforms.

Allosteric conformational change cascade in cytoplasmic dynein revealed by structure-based molecular simulations.

Cytoplasmic dynein is a giant ATP-driven molecular motor that proceeds to the minus end of the microtubule (MT). Dynein hydrolyzes ATP in a ring-like structure, containing 6 AAA+ (ATPases associated with diverse cellular activities) modules, which is ~15 nm away from the MT binding domain (MTBD). This architecture implies that long-distance allosteric couplings exist between the AAA+ ring and the MTBD in order for dynein to move on the MT, although little is known about the mechanisms involved. Here, we have performed comprehensive molecular simulations of the dynein motor domain based on pre- and post- power-stroke structural information and in doing so we address the allosteric conformational changes that occur during the power-stroke and recovery-stroke processes. In the power-stroke process, the N-terminal linker movement was the prerequisite to the nucleotide-dependent AAA1 transition, from which a transition cascade propagated, on average, in a circular manner on the AAA+ ring until it reached the AAA6/C-terminal module. The recovery-stroke process was initiated by the transition of the AAA6/C-terminal, from which the transition cascade split into the two directions of the AAA+ ring, occurring both clockwise and anti-clockwise. In both processes, the MTBD conformational change was regulated by the AAA4 module and the AAA5/Strut module.