A force-sensitive mutation reveals a non-canonical role for dynein in anaphase progression.

The diverse roles of the dynein motor in shaping microtubule networks and cargo transport complicate in vivo analysis of its functions significantly. To address this issue, we have generated a series of missense mutations in Drosophila Dynein heavy chain. We show that mutations associated with human neurological disease cause a range of defects, including impaired cargo trafficking in neurons. We also describe a novel microtubule-binding domain mutation that specifically blocks the metaphase-anaphase transition during mitosis in the embryo. This effect is independent from dynein's canonical role in silencing the spindle assembly checkpoint. Optical trapping of purified dynein complexes reveals that this mutation only compromises motor performance under load, a finding rationalized by the results of all-atom molecular dynamics simulations. We propose that dynein has a novel function in anaphase progression that depends on it operating in a specific load regime. More broadly, our work illustrates how in vivo functions of motors can be dissected by manipulating their mechanical properties.

The mechanism and energetics of the dynein priming stroke.

Dyneins are an AAA+ motor responsible for motility and force generation toward the minus end of microtubules. Dynein motility is powered by nucleotide-dependent transitions of its linker domain, which transitions between straight (post-powerstroke) and bent (pre-powerstroke) conformations. To understand the dynamics and energetics of the linker, we performed all-atom molecular dynamics simulations of human dynein-2 primed for its power stroke. Simulations revealed that the linker can adopt either a bent conformation or a semi-bent conformation, separated by a 5.7 kT energy barrier. The linker cannot switch back to its straight conformation in the pre-powerstroke state due to a steric clash with the AAA+ ring. Simulations also showed that an isolated linker has a free energy minimum near the semi-bent conformation in the absence of the AAA+ ring, indicating that the linker stores energy as it bends and releases this energy during the powerstroke.

A force-sensitive mutation reveals a spindle assembly checkpoint-independent role for dynein in anaphase progression.

The cytoplasmic dynein-1 (dynein) motor organizes cells by shaping microtubule networks and moving a large variety of cargoes along them. However, dynein's diverse roles complicate in vivo studies of its functions significantly. To address this issue, we have used gene editing to generate a series of missense mutations in Drosophila Dynein heavy chain (Dhc). We find that mutations associated with human neurological disease cause a range of defects in larval and adult flies, including impaired cargo trafficking in neurons. We also describe a novel mutation in the microtubule-binding domain (MTBD) of Dhc that, remarkably, causes metaphase arrest of mitotic spindles in the embryo but does not impair other dynein-dependent processes. We demonstrate that the mitotic arrest is independent of dynein's well-established roles in silencing the spindle assembly checkpoint. In vitro reconstitution and optical trapping assays reveal that the mutation only impairs the performance of dynein under load. In silico all-atom molecular dynamics simulations show that this effect correlates with increased flexibility of the MTBD, as well as an altered orientation of the stalk domain, with respect to the microtubule. Collectively, our data point to a novel role of dynein in anaphase progression that depends on the motor operating in a specific load regime. More broadly, our work illustrates how cytoskeletal transport processes can be dissected in vivo by manipulating mechanical properties of motors.

Omicron BA.1 and BA.2 variants increase the interactions of SARS-CoV-2 spike glycoprotein with ACE2.

SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. First the BA.1 and later the BA.2 variant became the most dominant strains of the Omicron variant. To investigate how the mutations of these strains affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water, and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ∼34% and ∼51% stronger binding free energies to PD for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells.

SARS-CoV-2 Delta Variant Decreases Nanobody Binding and ACE2 Blocking Effectivity.

The Delta variant spreads more rapidly than previous variants of SARS-CoV-2. This variant comprises several mutations on the receptor-binding domain (RBDDelta) of its spike glycoprotein, which binds to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. The RBD-PD interaction has been targeted by antibodies and nanobodies to prevent viral infection, but their effectiveness against the Delta variant remains unclear. Here, we investigated RBDDelta-PD interactions in the presence and absence of nanobodies H11-H4, H11-D4, and Ty1 by performing 21.8 µs of all-atom molecular dynamics simulations. Unbiased simulations revealed that Delta variant mutations strengthen RBD binding to ACE2 by increasing the hydrophobic interactions and salt bridge formation, but weaken interactions with H11-H4, H11-D4, and Ty1. Among these nanobodies H11-H4 and H11-D4 bind RBD without overlapping ACE2. They were unable to dislocate ACE2 from RBDDelta when bound side by side with ACE2 on RBD. Steered molecular dynamics simulations at comparable loading rates to high-speed atomic force microscopy (AFM) experiments estimated lower rupture forces of the nanobodies from RBDDelta compared to ACE2. Our results suggest that existing nanobodies are less effective to inhibit RBDDelta-PD interactions and a new generation of nanobodies is needed to neutralize the Delta variant.

Structural and functional insight into regulation of kinesin-1 by microtubule-associated protein MAP7.

Microtubule (MT)-associated protein 7 (MAP7) is a required cofactor for kinesin-1-driven transport of intracellular cargoes. Using cryo-electron microscopy and single-molecule imaging, we investigated how MAP7 binds MTs and facilitates kinesin-1 motility. The MT-binding domain (MTBD) of MAP7 bound MTs as an extended α helix between the protofilament ridge and the site of lateral contact. Unexpectedly, the MTBD partially overlapped with the binding site of kinesin-1 and inhibited its motility. However, by tethering kinesin-1 to the MT, the projection domain of MAP7 prevented dissociation of the motor and facilitated its binding to available neighboring sites. The inhibitory effect of the MTBD dominated as MTs became saturated with MAP7. Our results reveal biphasic regulation of kinesin-1 by MAP7 in the context of their competitive binding to MTs.

Binding Mechanism of Neutralizing Nanobodies Targeting SARS-CoV-2 Spike Glycoprotein.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters human cells upon binding of its spike (S) glycoproteins to ACE2 receptors. Several nanobodies neutralize SARS-CoV-2 infection by binding to the receptor-binding domain (RBD) of the S protein, but how their binding antagonizes S-ACE2 interactions is not well understood. Here, we identified interactions between the RBD and nanobodies H11-H4, H11-D4, and Ty1 by performing all-atom molecular dynamics simulations. H11-H4 and H11-D4 can bind to RBD without overlapping with ACE2. H11-H4, and to a lesser extent H11-D4, binding dislocates ACE2 from its binding site due to electrostatic repulsion. In comparison, Ty1 overlaps with ACE2 on RBD and has a similar binding strength to ACE2. Mutations in the Alpha variant of SARS-CoV-2 had a minor effect in RBD binding strengths of ACE2 and nanobodies, but reduced the ability of H11-H4 and H11-D4 to dislocate ACE2 from RBD. In comparison, the Beta variant weakened the RBD binding strengths of H11-H4 and H11-D4, which were less effective to dislocate ACE2 binding. Unexpectedly, mutations in Beta strengthened Ty1 binding to RBD, suggesting that this nanobody may be more effective to neutralize the Beta variant of SARS-CoV-2.

Critical Interactions Between the SARS-CoV-2 Spike Glycoprotein and the Human ACE2 Receptor.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects human cells by binding its spike (S) glycoproteins to angiotensin-converting enzyme 2 (ACE2) receptors and causes the coronavirus disease 2019 (COVID-19). Therapeutic approaches to prevent SARS-CoV-2 infection are mostly focused on blocking S-ACE2 binding, but critical residues that stabilize this interaction are not well understood. By performing all-atom molecular dynamics (MD) simulations, we identified an extended network of salt bridges, hydrophobic and electrostatic interactions, and hydrogen bonds between the receptor-binding domain (RBD) of the S protein and ACE2. Mutagenesis of these residues on the RBD was not sufficient to destabilize binding but reduced the average work to unbind the S protein from ACE2. In particular, the hydrophobic end of RBD serves as the main anchor site and is the last to unbind from ACE2 under force. We propose that blocking the hydrophobic surface of RBD via neutralizing antibodies could prove to be an effective strategy to inhibit S-ACE2 interactions.

Conformational transition of SARS-CoV-2 spike glycoprotein between its closed and open states.

In 2020, the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected millions of people worldwide and caused the coronavirus disease 2019 (COVID-19). Spike (S) glycoproteins on the viral membrane bind to ACE2 receptors on the host cell membrane and initiate fusion, and S protein is currently among the primary drug target to inhibit viral entry. The S protein can be in a receptor inaccessible (closed) or accessible (open) state based on down and up positions of its receptor-binding domain (RBD), respectively. However, conformational dynamics and the transition pathway between closed to open states remain unexplored. Here, we performed all-atom molecular dynamics (MD) simulations starting from closed and open states of the S protein trimer in the presence of explicit water and ions. MD simulations showed that RBD forms a higher number of interdomain interactions and exhibits lower mobility in its down position than its up position. MD simulations starting from intermediate conformations between the open and closed states indicated that RBD switches to the up position through a semi-open intermediate that potentially reduces the free energy barrier between the closed and open states. Free energy landscapes were constructed, and a minimum energy pathway connecting the closed and open states was proposed. Because RBD-ACE2 binding is compatible with the semi-open state, but not with the closed state of the S protein, we propose that the formation of the intermediate state is a prerequisite for the host cell recognition.

Molecular dynamics simulations provide molecular insights into the role of HLA-B51 in Behçet's disease pathogenesis.

Behçet's disease is an inflammatory disorder of unknown etiology. Genetic tendency has an important role in its pathogenesis, and HLA-B51, a class I MHC antigen, has been recognized as the strongest susceptibility factor for Behçet's disease. Despite the confirmation of the association of HLA-B51 with Behçet's disease in different populations, its pathogenic mechanisms remain elusive. HLA-B51 differs in only two amino acids from HLA-B52, other split antigen of HLA-B5, which is not associated with Behçet's disease. These two amino acids are located in the B pocket of the antigen-binding groove, which occupies the second amino acids of the bound peptides. To understand the nature of the HLA-peptide interactions, differences in structure and dynamics of two HLA alleles were investigated by molecular dynamics simulations using YAYDGKDYI, LPRSTVINI, and IPYQDLPHL peptides. For HLA-B51, all bound peptides fluctuated to larger extent than HLA-B52. Free energy profiles of unbinding process for YAYDGKDYI by steered molecular dynamics simulations showed that unbinding from HLA-B52 results in greater free energy differences than HLA-B51. These results suggest the possibility of an instability of HLA-B51 associated with the repertoire of peptides, and this finding may provide significant insight to its pathogenic role in Behçet's disease.

Thermodynamic first law efficiency of membrane proteins.

Proteins are nature's biomolecular machines. Proteins, such as transporters, pumps and motors, have complex function/operating-machinery/mechanisms, comparable to the macro-scaled machines that we encounter in our daily life. These proteins, as it is for their macro-scaled counterparts, convert (part of) other/various forms of energy into work. In this study, we are performing the first law analysis on a set of proteins, including the dopamine transporter, glycine transporters I and II, glutamate transporter, sodium-potassium pump and Ca2+ ATPase. Each of these proteins operates on a thermodynamic/mechanic cycle to perform their function. In each of these cycles, they receive energy from a source, convert part of this energy into work and reject the remaining part of the energy to the environment. Conservation of energy principle was applied to the thermodynamic/mechanic cycle of each protein, and thermodynamic first law efficiency was evaluated for each cycle, which shows how much of the energy input per cycle was converted into useful work. Interestingly, calculations based on experimental data indicate that proteins can operate under a range of efficiencies, which vary based on the extracellular and intracellular ion and substrate concentrations. The lowest observed first law efficiency was 50%, which is a very high value if compared to the efficiency of the macro-scaled heat engines we encounter in our daily lives.Communicated by Ramaswamy H. Sarma.