Polymer-like Model to Study the Dynamics of Dynamin Filaments on Deformable Membrane Tubes

Mesoscopic elasticity controls dynamin-driven fission of lipid tubules

Mesoscale physics bridges the gap between the microscopic degrees of freedom of a system and its large-scale continuous behavior and highlights the role of a few key quantities in complex and multiscale phenomena, like dynamin-driven fission of lipid membranes. The dynamin protein wraps the neck formed during clathrin-mediated endocytosis, for instance, and constricts it until severing occurs. Although ubiquitous and fundamental for life, the cooperation between the GTP-consuming conformational changes within the protein and the full-scale response of the underlying lipid substrate is yet to be unraveled. In this work, we build an effective mesoscopic model from constriction to fission of lipid tubules based on continuum membrane elasticity and implicitly accounting for ratchet-like power strokes of dynamins. Localization of the fission event, the overall geometry, and the energy expenditure we predict comply with the major experimental findings. This bolsters the idea that a continuous picture emerges soon enough to relate dynamin polymerization length and membrane rigidity and tension with the optimal pathway to fission. We therefore suggest that dynamins found in in vivo processes may optimize their structure accordingly. Ultimately, we shed light on real-time conductance measurements available in literature and predict the fission time dependency on elastic parameters.

Integrative modeling of guanylate binding protein dimers

Guanylate-binding proteins (GBPs) are essential interferon-γ-activated large GTPases that play a crucial role in host defense against intracellular bacteria and parasites. While their protective functions rely on protein polymerization, our understanding of the structural intricacies of these multimerized states remains limited. To bridge this knowledge gap, we present dimer models for human GBP1 (hGBP1) and murine GBP2 and 7 (mGBP2 and mGBP7) using an integrative approach, incorporating the crystal structure of hGBP1's GTPase domain dimer, crosslinking mass spectrometry, small-angle X-ray scattering, protein-protein docking, and molecular dynamics simulations. Our investigation begins by comparing the protein dynamics of hGBP1, mGBP2, and mGBP7. We observe that the M/E domain in all three proteins exhibits significant mobility and hinge motion, with mGBP7 displaying a slightly less pronounced motion but greater flexibility in its GTPase domain. These dynamic distinctions can be attributed to variations in the sequences of mGBP7 and hGBP1/mGBP2, resulting in different dimerization modes. Unlike hGBP1 and its close ortholog mGBP2, which exclusively dimerize through their GTPase domains, we find that mGBP7 exhibits three equally probable alternative dimer structures. The GTPase domain of mGBP7 is only partially involved in its dimerization, primarily due to an accumulation of negative charge, allowing mGBP7 to dimerize independently of GTP. Instead, mGBP7 exhibits a strong tendency to dimerize in an antiparallel arrangement across its stalks. The results of this work go beyond the sequence-structure-function relationship, as the sequence differences in mGBP7 and mGBP2/hGBP1 do not lead to different structures, but to different protein dynamics and dimerization. The distinct GBP dimer structures are expected to encode specific functions crucial for disrupting pathogen membranes.

Structural insights into the mechanism of GTP initiation of microtubule assembly

In eukaryotes, the dynamic assembly of microtubules (MT) plays an important role in numerous cellular processes. The underlying mechanism of GTP triggering MT assembly is still unknown. Here, we present cryo-EM structures of tubulin heterodimer at their GTP- and GDP-bound states, intermediate assembly states of GTP-tubulin, and final assembly stages of MT. Both GTP- and GDP-tubulin heterodimers adopt similar curved conformations with subtle flexibility differences. In head-to-tail oligomers of tubulin heterodimers, the inter-dimer interface of GDP-tubulin exhibits greater flexibility, particularly in tangential bending. Cryo-EM of the intermediate assembly states reveals two types of tubulin lateral contacts, "Tube-bond" and "MT-bond". Further, molecular dynamics (MD) simulations show that GTP triggers lateral contact formation in MT assembly in multiple sequential steps, gradually straightening the curved tubulin heterodimers. Therefore, we propose a flexible model of GTP-initiated MT assembly, including the formation of longitudinal and lateral contacts, to explain the nucleation and assembly of MT.

Cryo-electron tomography reveals structural insights into the membrane remodeling mode of dynamin-like EHD filaments

Eps15-homology domain containing proteins (EHDs) are eukaryotic, dynamin-related ATPases involved in cellular membrane trafficking. They oligomerize on membranes into filaments that induce membrane tubulation. While EHD crystal structures in open and closed conformations were previously reported, little structural information is available for the membrane-bound oligomeric form. Consequently, mechanistic insights into the membrane remodeling mechanism have remained sparse. Here, by using cryo-electron tomography and subtomogram averaging, we determined structures of nucleotide-bound EHD4 filaments on membrane tubes of various diameters at an average resolution of 7.6 Å. Assembly of EHD4 is mediated via interfaces in the G-domain and the helical domain. The oligomerized EHD4 structure resembles the closed conformation, where the tips of the helical domains protrude into the membrane. The variation in filament geometry and tube radius suggests a spontaneous filament curvature of approximately 1/70 nm^-1. Combining the available structural and functional data, we suggest a model for EHD-mediated membrane remodeling.

SMOG 2 and OpenSMOG: Extending the limits of structure-based models

Applying simulations with structure-based G o ¯ - like models has proven to be an effective strategy for investigating the factors that control biomolecular dynamics. The common element of these models is that some (or all) of the intra/inter-molecular interactions are explicitly defined to stabilize an experimentally determined structure. To facilitate the development and application of this broad class of models, we previously released the SMOG 2 software package. This suite allows one to easily customize and distribute structure-based (i.e., SMOG) models for any type of polymer-ligand system. The force fields generated by SMOG 2 may then be used to perform simulations in highly optimized MD packages, such as Gromacs, NAMD, LAMMPS, and OpenMM. Here, we describe extensions to the software and demonstrate the capabilities of the most recent version (SMOG v2.4.2). Changes include new tools that aid user-defined customization of force fields, as well as an interface with the OpenMM simulation libraries (OpenSMOG v1.1.0). The OpenSMOG module allows for arbitrary user-defined contact potentials and non-bonded potentials to be employed in SMOG models, without source-code modifications. To illustrate the utility of these advances, we present applications to systems with millions of atoms, long polymers and explicit ions, as well as models that include non-structure-based (e.g., AMBER-based) energetic terms. Examples include large-scale rearrangements of the SARS-CoV-2 Spike protein, the HIV-1 capsid with explicit ions, and crystallographic lattices of ribosomes and proteins. In summary, SMOG 2 and OpenSMOG provide robust support for researchers who seek to develop and apply structure-based models to large and/or intricate biomolecular systems.

Hydrodynamic coupling for particle-based solvent-free membrane models

The great challenge with biological membrane systems is the wide range of scales involved, from nanometers and picoseconds for individual lipids to the micrometers and beyond millisecond for cellular signaling processes. While solvent-free coarse-grained membrane models are convenient for large-scale simulations and promising to provide insight into slow processes involving membranes, these models usually have unrealistic kinetics. One major obstacle is the lack of an equally convenient way of introducing hydrodynamic coupling without significantly increasing the computational cost of the model. To address this, we introduce a framework based on anisotropic Langevin dynamics, for which major in-plane and out-of-plane hydrodynamic effects are modeled via friction and diffusion tensors from analytical or semi-analytical solutions to Stokes hydrodynamic equations. Using this framework, in conjunction with our recently developed membrane model, we obtain accurate dispersion relations for planar membrane patches, both free-standing and in the vicinity of a wall. We briefly discuss how non-equilibrium dynamics is affected by hydrodynamic interactions. We also measure the surface viscosity of the model membrane and discuss the affecting dissipative mechanisms.

CryoEM structure of the super-constricted two-start dynamin 1 filament

Dynamin belongs to the large GTPase superfamily, and mediates the fission of vesicles during endocytosis. Dynamin molecules are recruited to the neck of budding vesicles to assemble into a helical collar and to constrict the underlying membrane. Two helical forms were observed: the one-start helix in the constricted state and the two-start helix in the super-constricted state. Here we report the cryoEM structure of a super-constricted two-start dynamin 1 filament at 3.74 Å resolution. The two strands are joined by the conserved GTPase dimeric interface. In comparison with the one-start structure, a rotation around Hinge 1 is observed, essential for communicating the chemical power of the GTPase domain and the mechanical force of the Stalk and PH domain onto the underlying membrane. The Stalk interfaces are well conserved and serve as fulcrums for adapting to changing curvatures. Relative to one-start, small rotations per interface accumulate to bring a drastic change in the helical pitch. Elasticity theory rationalizes the diversity of dynamin helical symmetries and suggests corresponding functional significance.

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.