Sequence Determines the Switch in the Fibril Forming Regions in the Low-Complexity FUS Protein and Its Variants.

Residues spanning distinct regions of the low-complexity domain of the RNA-binding protein, Fused in Sarcoma (FUS-LC), form fibril structures with different core morphologies. Solid-state NMR experiments show that the 214-residue FUS-LC forms a fibril with an S-bend (core-1, residues 39-95), while the rest of the protein is disordered. In contrast, the fibrils of the C-terminal variant (FUS-LC-C; residues 111-214) have a U-bend topology (core-2, residues 112-150). Absence of the U-bend in FUS-LC implies that the two fibril cores do not coexist. Computer simulations show that these perplexing findings could be understood in terms of the population of sparsely populated fibril-like excited states in the monomer. The propensity to form core-1 is higher compared to core-2. We predict that core-2 forms only in truncated variants that do not contain the core-1 sequence. At the monomer level, sequence-dependent enthalpic effects determine the relative stabilities of the core-1 and core-2 topologies.

Molecular Transfer Model for pH Effects on Intrinsically Disordered Proteins: Theory and Applications.

We present a theoretical method to study how changes in pH shape the heterogeneous conformational ensemble explored by intrinsically disordered proteins (IDPs). The theory is developed in the context of coarse-grained models, which enable a fast, accurate, and extensive exploration of conformational space at a given protonation state. In order to account for pH effects, we generalize the molecular transfer model (MTM), in which conformations are re-weighted using the transfer free energy, which is the free energy necessary for bringing to equilibrium in a new environment a "frozen" conformation of the system. Using the semi-grand ensemble, we derive an exact expression of the transfer free energy, which amounts to the appropriate summation over all the protonation states. Because the exact result is computationally too demanding to be useful for large polyelectrolytes or IDPs, we introduce a mean-field (MF) approximation of the transfer free energy. Using a lattice model, we compare the exact and MF results for the transfer free energy and a variety of observables associated with the model IDP. We find that the precise location of the charged groups (the sequence), and not merely the net charge, determines the structural properties. We demonstrate that some of the limitations previously noted for MF theory in the context of globular proteins are mitigated when disordered polymers are studied. The excellent agreement between the exact and MF results poises us to use the method presented here as a computational tool to study the properties of IDPs and other biological systems as a function of pH.

Step-Wise Hydration of Magnesium by Four Water Molecules Precedes Phosphate Release in a Myosin Motor.

Molecular motors, such as myosin, kinesin, and dynein, convert the energy released by the hydrolysis of ATP into mechanical work, thus allowing them to undergo directional motion on cytoskeletal tracks. A pivotal step in the chemomechanical transduction in myosin motors occurs after they bind to the actin filament, which triggers the release of phosphate (Pi, product of ATP hydrolysis) and the rotation of the lever arm. Here, we investigate the mechanism of phosphate release in myosin VI using extensive molecular dynamics simulations involving multiple trajectories of several µs. Because the escape of phosphate is expected to occur on time-scales on the order of milliseconds or more in myosin VI, we observed Pi release only if the trajectories were initiated with a rotated phosphate inside the nucleotide binding pocket. We discovered that although Pi populates the traditional "back door" route, phosphate exits through various other gateways, thus establishing the heterogeneity in the escape routes. Remarkably, we observed that the release of phosphate is preceded by a stepwise hydration of the ADP-bound magnesium ion. The release of the anion occurred only after four water molecules hydrated the cation (Mg2+). By performing comparative structural analyses, we show that hydration of magnesium is the key step in the phosphate release in a number of ATPases and GTPases. Nature may have evolved hydration of Mg2+ as a general molecular switch for Pi release, which is a universal step in the catalytic cycle of many machines that share little sequence or structural similarity.

Dynamics of Allosteric Transitions in Dynein.

Cytoplasmic dynein, whose motor domain belongs to the AAA+ family, walks on microtubules toward the minus end. Using the available structures in different nucleotide states, we performed simulations of a coarse-grained model to elucidate the dynamics of allosteric transitions. Binding of ATP closes the cleft between the AAA1 and AAA2 domains, triggering conformational changes in the rest of the motor domain, thus forming the pre-power stroke state. Interactions with the microtubule, modeled implicitly, enhance ADP release rate, and the formation of the post-power stroke state. The dynamics of the linker (LN), which reversibly changes from a straight to a bent state, is heterogeneous. Persistent interactions between the LN and the insert loops in the AAA2 domain prevent the formation of pre-power stroke state when ATP is bound to AAA3, thus locking dynein in a repressed non-functional state. Application of mechanical force to the LN restores motility in the repressed state.