Outcomes of the EMDataResource cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

Bending-torsional elasticity and energetics of the plus-end microtubule tip.

SignificanceThe mechanochemical basis of microtubule growth, which is essential for the normal function and division of eukaryotic cells, has remained elusive and controversial, despite extensive work. In particular, recent findings have created the paradox that the microtubule plus-end tips look very similar during both growing and shrinking phases, thereby challenging the traditional textbook picture. Our large-scale atomistic simulations resolve this paradox and explain microtubule growth and shrinkage dynamics as a process governed by energy barriers between protofilament conformations, the heights of which are in turn fine-tuned by different nucleotide states, thus implementing an information-driven Brownian ratchet.

Cryo-EM model validation recommendations based on outcomes of the 2019 EMDataResource challenge.

This paper describes outcomes of the 2019 Cryo-EM Model Challenge. The goals were to (1) assess the quality of models that can be produced from cryogenic electron microscopy (cryo-EM) maps using current modeling software, (2) evaluate reproducibility of modeling results from different software developers and users and (3) compare performance of current metrics used for model evaluation, particularly Fit-to-Map metrics, with focus on near-atomic resolution. Our findings demonstrate the relatively high accuracy and reproducibility of cryo-EM models derived by 13 participating teams from four benchmark maps, including three forming a resolution series (1.8 to 3.1 Å). The results permit specific recommendations to be made about validating near-atomic cryo-EM structures both in the context of individual experiments and structure data archives such as the Protein Data Bank. We recommend the adoption of multiple scoring parameters to provide full and objective annotation and assessment of the model, reflective of the observed cryo-EM map density.

Choice of fluorophore affects dynamic DNA nanostructures.

The ability to dynamically remodel DNA origami structures or functional nanodevices is highly desired in the field of DNA nanotechnology. Concomitantly, the use of fluorophores to track and validate the dynamics of such DNA-based architectures is commonplace and often unavoidable. It is therefore crucial to be aware of the side effects of popular fluorophores, which are often exchanged without considering the potential impact on the system. Here, we show that the choice of fluorophore can strongly affect the reconfiguration of DNA nanostructures. To this end, we encapsulate a triple-stranded DNA (tsDNA) into water-in-oil compartments and functionalize their periphery with a single-stranded DNA handle (ssDNA). Thus, the tsDNA can bind and unbind from the periphery by reversible opening of the triplex and subsequent strand displacement. Using a combination of experiments, molecular dynamics (MD) simulations, and reaction-diffusion modelling, we demonstrate for 12 different fluorophore combinations that it is possible to alter or even inhibit the DNA nanostructure formation-without changing the DNA sequence. Besides its immediate importance for the design of pH-responsive switches and fluorophore labelling, our work presents a strategy to precisely tune the energy landscape of dynamic DNA nanodevices.

Microtubule instability driven by longitudinal and lateral strain propagation.

Tubulin dimers associate longitudinally and laterally to form metastable microtubules (MTs). MT disassembly is preceded by subtle structural changes in tubulin fueled by GTP hydrolysis. These changes render the MT lattice unstable, but it is unclear exactly how they affect lattice energetics and strain. We performed long-time atomistic simulations to interrogate the impacts of GTP hydrolysis on tubulin lattice conformation, lateral inter-dimer interactions, and (non-)local lateral coordination of dimer motions. The simulations suggest that most of the hydrolysis energy is stored in the lattice in the form of longitudinal strain. While not significantly affecting lateral bond stability, the stored elastic energy results in more strongly confined and correlated dynamics of GDP-tubulins, thereby entropically destabilizing the MT lattice.

A refined reaction-diffusion model of tau-microtubule dynamics and its application in FDAP analysis.

Fluorescence decay after photoactivation (FDAP) and fluorescence recovery after photobleaching (FRAP) are well established approaches for studying the interaction of the microtubule (MT)-associated protein tau with MTs in neuronal cells. Previous interpretations of FDAP/FRAP data have revealed dwell times of tau on MTs in the range of several seconds. However, this is difficult to reconcile with a dwell time recently measured by single-molecule analysis in neuronal processes that was shorter by two orders of magnitude. Questioning the validity of previously used phenomenological interpretations of FDAP/FRAP data, we have generalized the standard two-state reaction-diffusion equations by 1), accounting for the parallel and discrete arrangement of MTs in cell processes (i.e., homogeneous versus heterogeneous distribution of tau-binding sites); and 2), explicitly considering both active (diffusion upon MTs) and passive (piggybacking upon MTs at rates of slow axonal transport) motion of bound tau. For some idealized cases, analytical solutions were derived. By comparing them with the full numerical solution and Monte Carlo simulations, the respective validity domains were mapped. Interpretation of our FDAP data (from processes of neuronally differentiated PC12 cells) in light of the heterogeneous formalism yielded independent estimates for the association (∼2 ms) and dwell (∼100 ms) times of tau to/on a single MT rather than in an MT array. The dwell time was shorter by orders of magnitude than that in a previous report where a homogeneous topology of MTs was assumed. We found that the diffusion of bound tau was negligible in vivo, in contrast to an earlier report that tau diffuses along the MT lattice in vitro. Methodologically, our results demonstrate that the heterogeneity of binding sites cannot be ignored when dealing with reaction-diffusion of cytoskeleton-associated proteins. Physiologically, the results reveal the behavior of tau in cellular processes, which is noticeably different from that in vitro.

Single-particle Cryo-EM and molecular dynamics simulations: A perfect match.

Knowledge of the structure and dynamics of biomolecules is key to understanding the mechanisms underlying their biological functions. Single-particle cryo-electron microscopy (cryo-EM) is a powerful structural biology technique to characterize complex biomolecular systems. Here, we review recent advances of how Molecular Dynamics (MD) simulations are being used to increase and enhance the information extracted from cryo-EM experiments. We will particularly focus on the physics underlying these experiments, how MD facilitates structure refinement, in particular for heterogeneous and non-isotropic resolution, and how thermodynamic and kinetic information can be extracted from cryo-EM data.

Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge.

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

SESCA: Predicting Circular Dichroism Spectra from Protein Molecular Structures.

Circular dichroism (CD) spectroscopy is a highly sensitive but low-resolution technique to study the structure of proteins. Combined with molecular modeling or other complementary techniques, CD spectroscopy can provide essential information at higher resolution. To this end, we introduce a new computational method to calculate the electronic circular dichroism spectra of proteins from a structural model or ensemble using the average secondary structure composition and a precalculated set of basis spectra. The method is designed for model validation to estimate the error of a given protein structural model based on the measured CD spectrum. We compared the predictive power of our method to that of existing algorithms, namely, DichroCalc and PDB2CD, and found that it predicts CD spectra more accurately. Our results indicate that the derived basis sets are robust to both experimental errors in the reference spectra and the choice of the secondary structure classification algorithm. For over 80% of the globular reference proteins, our basis sets accurately predict the experimental spectrum solely from their secondary structure composition. For the remaining 20%, correcting for intensity normalization considerably improves the prediction power. Additionally, we show that the predictions for short peptides and an example complex of intrinsically disordered proteins strongly benefit from accounting for side-chain contributions and structural flexibility.

Automated cryo-EM structure refinement using correlation-driven molecular dynamics.

We present a correlation-driven molecular dynamics (CDMD) method for automated refinement of atomistic models into cryo-electron microscopy (cryo-EM) maps at resolutions ranging from near-atomic to subnanometer. It utilizes a chemically accurate force field and thermodynamic sampling to improve the real-space correlation between the modeled structure and the cryo-EM map. Our framework employs a gradual increase in resolution and map-model agreement as well as simulated annealing, and allows fully automated refinement without manual intervention or any additional rotamer- and backbone-specific restraints. Using multiple challenging systems covering a wide range of map resolutions, system sizes, starting model geometries and distances from the target state, we assess the quality of generated models in terms of both model accuracy and potential of overfitting. To provide an objective comparison, we apply several well-established methods across all examples and demonstrate that CDMD performs best in most cases.

Microtubule assembly governed by tubulin allosteric gain in flexibility and lattice induced fit.

Microtubules (MTs) are key components of the cytoskeleton and play a central role in cell division and development. MT assembly is known to be associated with a structural change in [Formula: see text]-tubulin dimers from kinked to straight conformations. How GTP binding renders individual dimers polymerization-competent, however, is still unclear. Here, we have characterized the conformational dynamics and energetics of unassembled tubulin using atomistic molecular dynamics and free energy calculations. Contrary to existing allosteric and lattice models, we find that GTP-tubulin favors a broad range of almost isoenergetic curvatures, whereas GDP-tubulin has a much lower bending flexibility. Moreover, irrespective of the bound nucleotide and curvature, two conformational states exist differing in location of the anchor point connecting the monomers that affects tubulin bending, with one state being strongly favored in solution. Our findings suggest a new combined model in which MTs incorporate and stabilize flexible GTP-dimers with a specific anchor point state.

Single-molecule imaging reveals dynamic biphasic partition of RNA-binding proteins in stress granules.

Stress granules (SGs) are cytosolic, nonmembranous RNA-protein complexes. In vitro experiments suggested that they are formed by liquid-liquid phase separation; however, their properties in mammalian cells remain unclear. We analyzed the distribution and dynamics of two paradigmatic RNA-binding proteins (RBPs), Ras GTPase-activating protein SH3-domain-binding protein (G3BP1) and insulin-like growth factor II mRNA-binding protein 1 (IMP1), with single-molecule resolution in living neuronal cells. Both RBPs exhibited different exchange kinetics between SGs. Within SGs, single-molecule localization microscopy revealed distributed hotspots of immobilized G3BP1 and IMP1 that reflect the presence of relatively immobile nanometer-sized nanocores. We demonstrate alternating binding in nanocores and anomalous diffusion in the liquid phase with similar characteristics for both RBPs. Reduction of low-complexity regions in G3BP1 resulted in less detectable mobile molecules in the liquid phase without change in binding in nanocores. The data provide direct support for liquid droplet behavior of SGs in living cells and reveal transient binding of RBPs in nanocores. Our study uncovers a surprising disconnect between SG partitioning and internal diffusion and interactions of RBPs.

Presence of a carboxy-terminal pseudorepeat and disease-like pseudohyperphosphorylation critically influence tau's interaction with microtubules in axon-like processes.

A current challenge of cell biology is to investigate molecular interactions in subcellular compartments of living cells to overcome the artificial character of in vitro studies. To dissect the interaction of the neuronal microtubule (MT)-associated protein tau with MTs in axon-like processes, we used a refined fluorescence decay after photoactivation approach and single-molecule tracking. We found that isoform variation had only a minor influence on the tau-MT interaction, whereas the presence of a C-terminal pseudorepeat region (PRR) greatly increased MT binding by a greater-than-sixfold reduction of the dissociation rate. Bioinformatic analysis revealed that the PRR contained a highly conserved motif of 18 amino acids. Disease-associated tau mutations in the PRR (K369I, G389R) did not influence apparent MT binding but increased its dynamicity. Simulation of disease-like tau hyperphosphorylation dramatically diminished the tau-MT interaction by a greater-than-fivefold decrease of the association rate with no major change in the dissociation rate. Apparent binding of tau to MTs was similar in axons and dendrites but more sensitive to increased phosphorylation in axons. Our data indicate that under the conditions of high MT density that prevail in the axon, tau's MT binding and localization are crucially affected by the presence of the PRR and tau hyperphosphorylation.

Single-molecule tracking of tau reveals fast kiss-and-hop interaction with microtubules in living neurons.

The microtubule-associated phosphoprotein tau regulates microtubule dynamics and is involved in neurodegenerative diseases collectively called tauopathies. It is generally believed that the vast majority of tau molecules decorate axonal microtubules, thereby stabilizing them. However, it is an open question how tau can regulate microtubule dynamics without impeding microtubule-dependent transport and how tau is also available for interactions other than those with microtubules. Here we address this apparent paradox by fast single-molecule tracking of tau in living neurons and Monte Carlo simulations of tau dynamics. We find that tau dwells on a single microtubule for an unexpectedly short time of ∼40 ms before it hops to the next. This dwell time is 100-fold shorter than previously reported by ensemble measurements. Furthermore, we observed by quantitative imaging using fluorescence decay after photoactivation recordings of photoactivatable GFP-tagged tubulin that, despite this rapid dynamics, tau is capable of regulating the tubulin-microtubule balance. This indicates that tau's dwell time on microtubules is sufficiently long to influence the lifetime of a tubulin subunit in a GTP cap. Our data imply a novel kiss-and-hop mechanism by which tau promotes neuronal microtubule assembly. The rapid kiss-and-hop interaction explains why tau, although binding to microtubules, does not interfere with axonal transport.

Interplay between phosphorylation and palmitoylation mediates plasma membrane targeting and sorting of GAP43.

Phosphorylation and lipidation provide posttranslational mechanisms that contribute to the distribution of cytosolic proteins in growing nerve cells. The growth-associated protein GAP43 is susceptible to both phosphorylation and S-palmitoylation and is enriched in the tips of extending neurites. However, how phosphorylation and lipidation interplay to mediate sorting of GAP43 is unclear. Using a combination of biochemical, genetic, and imaging approaches, we show that palmitoylation is required for membrane association and that phosphorylation at Ser-41 directs palmitoylated GAP43 to the plasma membrane. Plasma membrane association decreased the diffusion constant fourfold in neuritic shafts. Sorting to the neuritic tip required palmitoylation and active transport and was increased by phosphorylation-mediated plasma membrane interaction. Vesicle tracking revealed transient association of a fraction of GAP43 with exocytic vesicles and motion at a fast axonal transport rate. Simulations confirmed that a combination of diffusion, dynamic plasma membrane interaction and active transport of a small fraction of GAP43 suffices for efficient sorting to growth cones. Our data demonstrate a complex interplay between phosphorylation and lipidation in mediating the localization of GAP43 in neuronal cells. Palmitoylation tags GAP43 for global sorting by piggybacking on exocytic vesicles, whereas phosphorylation locally regulates protein mobility and plasma membrane targeting of palmitoylated GAP43.