Strain Mediated Adaptation Is Key for Myosin Mechanochemistry: Discovering General Rules for Motor Activity

How fingers affect folding of a thumb: Inter-subdomain cooperation in the folding of SARS-CoV-2 RdRp protein

The RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 is a critical enzyme essential for the virus's replication and transcription, making it a key therapeutic target. The RdRp protein exhibits a characteristic cupped right-hand shaped structure with two vital subdomains: the fingers and the thumb. Despite being distinct, biophysical experiments suggest that these subdomains cooperate to facilitate RNA accommodation, ensuring RdRp functionality. To investigate the structure-based mechanisms underlying the fingers-thumb interaction in both apo and RNA-bound RdRp, we constructed a coarse-grained structure-based model based on recent cryo-electron microscopy data. The simulations reveal frequent open-to-closed conformational transitions in apo RdRp, akin to a breathing-like motion. These conformational changes are regulated by the fingers-thumb association and the folding dynamics of the thumb subdomain. The thumb adopts a stable fold only when tethered by the fingers-thumb interface; when these subdomains are disconnected, the thumb transitions into an open state. A significant number of open-to-closed transition events were analyzed to generate a transition contact probability map, which highlights a few specific residues at the thumb-fingers interface, distant from the RNA accommodation sites, as essential for inducing the thumb's folding process. Given that thumb subdomain folding is critical for RNA binding and viral replication, the study proposes that these interfacial residues may function as remote regulatory switches and could be targeted for the development of allosteric drugs against SARS-CoV-2 and similar RNA viruses.

On the Control of Directionality of Myosin

The origin of the unique directionality of myosin has been a problem of fundamental and practical importance. This work establishes in a conclusive way that the directionality is controlled by tuning the barrier for the rate-determining step, namely, the ADP release step. This conclusion is based on exploring the molecular origin behind the reverse directionality of myosins V and VI and the determination of the origin of the change in the barriers of the ADP release for the forward and backward motions. Our investigation is performed by combining different simulation methods such as steer molecular dynamics (SMD), umbrella sampling, renormalization method, and automated path searching method. It is found that in the case of myosin V, the ADP release from the postrigor (trailing head) state overcomes a lower barrier than the prepowerstroke (leading head) state, which is also evident from experimental observation. In the case of myosin VI, we noticed a different trend when compared to myosin V. Since the directionality of myosins V and VI follows a reverse trend, we conclude that such differences in the directionality are controlled by the free energy barrier for the ADP release. Overall, the proof that the directionality of myosin is determined by the activation barrier of the rate-determining step in the cycle, rather than by some unspecified dynamical effects, has general importance.

Intrinsically Disordered Regions Function as a Cervical Collar to Remotely Regulate the Nodding Dynamics of SARS-CoV-2 Prefusion Spike Heads

The SARS-CoV-2 prefusion spike heads (receptor binding domains, RBDs) frequently nod down and up to interact with host cell receptors. As the spike protein is a trimeric unit of significant size, to understand its large-scale structural dynamics associated with the nodding mechanism and the mutational impact on the same, we develop a topological symmetry-information-loaded coarse-grained structure-based model of a spike trimer using recent cryo-EM structural data. Our study reveals the control of two distant intrinsically disordered regions (IDRs), namely, 630 and FPPR loops, over the nodding dynamics of spike heads. We find that the order-disorder transition of IDRs becomes more evident in the variants of concern (VOCs) that are associated with the characteristic mutation, D614G, in the proximity of these IDRs. In some VOCs, the two other mutations A570D and S982A also show an integral effect. The driver mutation D614G instigates a salt-bridge disruption, altering the order-disorder dynamics of both 630 and FPPR loops and their interaction with the C-terminal domains (CTD1/CTD2). This altered connectivity in these mutants allows the two IDRs to act collectively as a "cervical collar" for the RBD, supporting various spike head postures, consistent with cryo-EM results available for specific cases. The IDRs' control over the spike structure and dynamics presents an exciting opportunity where they can be targeted as remote operational switches to artificially maneuver the nod for effective therapeutic interventions.

A VDAC1-mediated NEET protein chain transfers [2Fe-2S] clusters between the mitochondria and the cytosol and impacts mitochondrial dynamics

Here we address the important question of cross-talk between the mitochondria and cytosol. We show that the inner mitochondrial protein, MiNT, interacts with a protein on the outer mitochondrial membrane (mNT). This interaction occurs within the major outer membrane protein VDAC1. Inside the inner space of VDAC1, MiNT transfers its [2Fe-2S] clusters to mNT, which was shown to be a [2Fe-2S] cluster donor protein that donates its cluster(s) to apo-acceptor proteins residing in the cytosol. Hence, we suggest a pathway for transferring [2Fe-2S] clusters from inside the mitochondria to the cytosol.

Mitochondrial inner NEET (MiNT) and the outer mitochondrial membrane (OMM) mitoNEET (mNT) proteins belong to the NEET protein family. This family plays a key role in mitochondrial labile iron and reactive oxygen species (ROS) homeostasis. NEET proteins contain labile [2Fe-2S] clusters which can be transferred to apo-acceptor proteins. In eukaryotes, the biogenesis of [2Fe-2S] clusters occurs within the mitochondria by the iron–sulfur cluster (ISC) system; the clusters are then transferred to [2Fe-2S] proteins within the mitochondria or exported to cytosolic proteins and the cytosolic iron–sulfur cluster assembly (CIA) system. The last step of export of the [2Fe-2S] is not yet fully characterized. Here we show that MiNT interacts with voltage-dependent anion channel 1 (VDAC1), a major OMM protein that connects the intermembrane space with the cytosol and participates in regulating the levels of different ions including mitochondrial labile iron (mLI). We further show that VDAC1 is mediating the interaction between MiNT and mNT, in which MiNT transfers its [2Fe-2S] clusters from inside the mitochondria to mNT that is facing the cytosol. This MiNT–VDAC1–mNT interaction is shown both experimentally and by computational calculations. Additionally, we show that modifying MiNT expression in breast cancer cells affects the dynamics of mitochondrial structure and morphology, mitochondrial function, and breast cancer tumor growth. Our findings reveal a pathway for the transfer of [2Fe-2S] clusters, which are assembled inside the mitochondria, to the cytosol.

Quantification and demonstration of the collective constriction-by-ratchet mechanism in the dynamin molecular motor

Dynamin oligomerizes into helical filaments on tubular membrane templates and, through constriction, cleaves them in a GTPase-driven way. Structural observations of GTP-dependent cross-bridges between neighboring filament turns have led to the suggestion that dynamin operates as a molecular ratchet motor. However, the proof of such mechanism remains absent. Particularly, it is not known whether a powerful enough stroke is produced and how the motor modules would cooperate in the constriction process. Here, we characterized the dynamin motor modules by single-molecule Förster resonance energy transfer (smFRET) and found strong nucleotide-dependent conformational preferences. Integrating smFRET with molecular dynamics simulations allowed us to estimate the forces generated in a power stroke. Subsequently, the quantitative force data and the measured kinetics of the GTPase cycle were incorporated into a model including both a dynamin filament, with explicit motor cross-bridges, and a realistic deformable membrane template. In our simulations, collective constriction of the membrane by dynamin motor modules, based on the ratchet mechanism, is directly reproduced and analyzed. Functional parallels between the dynamin system and actomyosin in the muscle are seen. Through concerted action of the motors, tight membrane constriction to the hemifission radius can be reached. Our experimental and computational study provides an example of how collective motor action in megadalton molecular assemblies can be approached and explicitly resolved.

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.

Dynamic Asymmetry Exposes 2019-nCoV Prefusion Spike

The novel coronavirus (2019-nCoV) spike protein is a smart molecular machine that instigates the entry of coronavirus to the host cell causing the COVID-19 pandemic. In this study, a symmetry-information-loaded structure-based Hamiltonian is developed using recent Cryo-EM structural data to explore the complete conformational energy landscape of the full-length prefusion spike protein. The study finds the 2019-nCoV prefusion spike to adopt a unique strategy by undertaking a dynamic conformational asymmetry that results in two prevalent asymmetric structures of spike where one or two spike heads rotate up to provide better exposure to the host-cell receptor. A few unique interchain interactions are identified at the interface of closely associated N-terminal domain (NTD) and receptor binding domain (RBD) playing a crucial role in the thermodynamic stabilization of the up conformation of the RBD in the case of the 2019-nCoV spike. The interaction-level information decoded in this study may provide deep insight into developing effective therapeutic targets.

Role of AAA3 Domain in Allosteric Communication of Dynein Motor Proteins

Cytoplasmic dynein, an AAA+ motif containing motor, generates force and movement along the microtubule to execute important biological functions including intracellular material transport and cell division by hydrolyzing ATP. Among the six AAA+ domains, AAA1 is the primary ATPase site where a single ATP hydrolysis generates a single step. Nucleotide states in AAA3 gate dynein's activity, suggesting that AAA3 acts as a regulatory switch. However, the comprehensive structural perspective of AAA3 in dynein's mechanochemical cycle remains unclear. Here, we explored the allosteric transition path of dynein involving AAA3 using a coarse-grained structure-based model. ATP binding to the AAA1 domain creates a cascade of conformational changes through the other domains of the ring, which leads to the pre-power stroke formation. The linker domain, which is the mechanical element of dynein, shifts from a straight to a bent conformation during this process. In our present study, we observe that AAA3 gates the allosteric communication from AAA1 to the microtubule binding domain (MTBD) through AAA4 and AAA5. The MTBD is linked to the AAA+ ring via a coiled-coil stalk and a buttress domain, which are extended from AAA4 and AAA5, respectively. Further analysis also uncovers the role of AAA3 in the linker movement. The free energy calculation shows that the linker prefers the straight conformation when AAA3 remains in the ATP-bound condition. As AAA3 restricts the motion of AAA4 and AAA5, the linker/AAA5 interactions get stabilized, and the linker cannot move to the pre-power stroke state that halts the complete structural transition required for the mechanochemical cycle. Therefore, we suggest that AAA3 governs dynein's mechanochemical cycle and motility by controlling the AAA4 and AAA5 domains that further regulate the linker movement and the power stroke formation.

Atomistic Models from Orientation and Distance Constraints Using EPR of a Bifunctional Spin Label

We have used high-resolution orientation and distance measurements derived from electron paramagnetic resonance of a bifunctional spin label (BSL) to build and refine atomistic models of protein structure. We demonstrate this approach by investigating the effects of nucleotide binding on the structure of myosin's catalytic domain while myosin is in complex with actin. Constraints for orientation of individual helices were obtained in a previous study from continuous-wave electron paramagnetic resonance of myosin labeled at specific sites with BSLs in oriented muscle fibers. In this study, new distance constraints were derived from double electron-electron resonance on myosin constructs labeled with a BSL specifically at two sites. Using these complementary constraints together, we thoroughly characterize the BSL's rigid, highly stereoselective attachment to protein α-helices, which permits accurate measurements of orientation and distance. We also leverage these measurements to derive a novel, to our knowledge, structural model for myosin-II in complex with actin and MgADP and compare our model to other recent actomyosin structures. The described approach is applicable to any orientable complex (e.g., membranes or filaments) in which site-specific di-Cys mutation is feasible.

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.

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.

Kinematics of the lever arm swing in myosin VI

Myosin VI (MVI) is the only known member of the myosin superfamily that, upon dimerization, walks processively toward the pointed end of the actin filament. The leading head of the dimer directs the trailing head forward with a power stroke, a conformational change of the motor domain exaggerated by the lever arm. Using a unique coarse-grained model for the power stroke of a single MVI, we provide the molecular basis for its motility. We show that the power stroke occurs in two major steps. First, the motor domain attains the poststroke conformation without directing the lever arm forward; and second, the lever arm reaches the poststroke orientation by undergoing a rotational diffusion. From the analysis of the trajectories, we discover that the potential that directs the rotating lever arm toward the poststroke conformation is almost flat, implying that the lever arm rotation is mostly uncoupled from the motor domain. Because a backward load comparable to the largest interhead tension in a MVI dimer prevents the rotation of the lever arm, our model suggests that the leading-head lever arm of a MVI dimer is uncoupled, in accord with the inference drawn from polarized total internal reflection fluorescence (polTIRF) experiments. Without any adjustable parameter, our simulations lead to quantitative agreement with polTIRF experiments, which validates the structural insights. Finally, in addition to making testable predictions, we also discuss the implications of our model in explaining the broad step-size distribution of the MVI stepping pattern.

Exploring the mechanochemical cycle of dynein motor proteins: structural evidence of crucial intermediates

Exploration of the biologically relevant pathways of dynein's mechanochemical cycle using structure based models.