Search

Förster resonance energy transfer (FRET) and electron paramagnetic resonance (EPR) spectroscopy are complementary techniques for quantifying distances in the nanometer range. Both approaches are commonly employed for probing the conformations and conformational changes of biological macromolecules based on site-directed fluorescent or paramagnetic labeling. FRET can be applied in solution at ambient temperature and thus provides direct access to dynamics, especially if used at the single-molecule level, whereas EPR requires immobilization or work at cryogenic temperatures but provides data that can be more reliably used to extract distance distributions. However, a combined analysis of the complementary data from the two techniques has been complicated by the lack of a common modeling framework. Here, we demonstrate a systematic analysis approach based on rotamer libraries for both FRET and EPR labels to predict distance distributions between two labels from a structural model. Dynamics of the fluorophores within these distance distributions are taken into account by diffusional averaging, which improves the agreement with experiment. Benchmarking this methodology with a series of surface-exposed pairs of sites in a structured protein domain reveals that the lowest resolved distance differences can be as small as ∼0.25 nm for both techniques, with quantitative agreement between experimental and simulated transfer efficiencies within a range of ±0.045. Rotamer library analysis thus establishes a coherent way of treating experimental data from EPR and FRET and provides a basis for integrative structural modeling, including studies of conformational distributions and dynamics of biological macromolecules using both techniques., Biophysical Journal, 120 (21), ISSN:0006-3495, ISSN:1542-0086

Damage to arterial vessel walls leads to the formation of platelet aggregate, which acts as a physical obstacle for bleeding. An arterial thrombus is heterogeneous; it has a dense inner part (core) and an unstable outer part (shell). The thrombus shell is very dynamic, being composed of loosely connected discoid platelets. The mechanisms underlying the observed mobility of the shell and its (patho)physiological implications are unclear. To investigate arterial thrombus mechanics, we developed a novel, to our knowledge, two-dimensional particle-based computational model of microvessel thrombosis. The model considers two types of interplatelet interactions: primary reversible (glycoprotein Ib (GPIb)-mediated) and stronger integrin-mediated interaction, which intensifies with platelet activation. At high shear rates, the former interaction leads to adhesion, and the latter is primarily responsible for stable platelet aggregation. Using a stochastic model of GPIb-mediated interaction, we initially reproduced experimental curves that characterize individual platelet interactions with a von Willebrand factor-coated surface. The addition of the second stabilizing interaction results in thrombus formation. The comparison of thrombus dynamics with experimental data allowed us to estimate the magnitude of critical interplatelet forces in the thrombus shell and the characteristic time of platelet activation. The model predicts moderate dependence of maximal thrombus height on the injury size in the absence of thrombin activity. We demonstrate that the developed stochastic model reproduces the observed highly dynamic behavior of the thrombus shell. The presence of primary stochastic interaction between platelets leads to the properties of thrombus consistent with in vivo findings; it does not grow upstream of the injury site and covers the whole injury from the first seconds of the formation. А simplified model, in which GPIb-mediated interaction is deterministic, does not reproduce these features. Thus, the stochasticity of platelet interactions is critical for thrombus plasticity, suggesting that interaction via a small number of bonds drives the dynamics of arterial thrombus shell.

N-methyl-D-aspartate (NMDA) receptors mediate synaptic excitatory signaling in the mammalian central nervous system by forming calcium-permeable transmembrane channels upon binding glutamate and coagonist glycine. Ca2+ influx through NMDA receptors leads to channel inactivation through a process mediated by resident calmodulin bound to the intracellular C-terminal segment of the GluN1 subunit of the receptor. Using single-molecule FRET investigations, we show that in the presence of calcium-calmodulin, the distance across the two GluN1 subunits at the entrance of the first transmembrane segment is shorter and the bilobed cleft of the glycine-binding domain in GluN1 is more closed when bound to glycine and glutamate relative to what is observed in the presence of barium-calmodulin. Consistent with these observations, the glycine deactivation rate is slower in the presence of calcium-calmodulin. Taken together, these results show that the binding of calcium-calmodulin to the C-terminus has long-range allosteric effects on the extracellular segments of the receptor that may contribute to the calcium-dependent inactivation.

Bacterial cell division is tightly coupled to the dynamic behavior of FtsZ, a tubulin homolog. Recent experimental work in vitro and in vivo has attributed FtsZ's assembly dynamics to treadmilling, in which subunits add to the bottom and dissociate from the top of protofilaments. However, the molecular mechanisms producing treadmilling have yet to be characterized and quantified. We have developed a Monte Carlo model for FtsZ assembly that explains treadmilling, and also explains assembly nucleation by the same mechanisms. A key element of the model is a conformational change from R (relaxed), which is highly favored for monomers, to T (tense), which is favored for subunits in a protofilament. This model was created in MATLAB. Kinetic parameters were converted to probabilities of execution during a single, small time step. These were used to stochastically determine FtsZ dynamics. Our model is able to accurately describe the results of several in vitro and in vivo studies for a variety of FtsZ flavors. With standard conditions, the model FtsZ polymerized and produced protofilaments that treadmilled at 23 nm/s, hydrolyzed GTP at 3.6-3.7 GTP min-1 FtsZ-1, and had an average length of 30-40 subunits, all similar to experimental results. Adding a bottom capper resulted in shorter protofilaments and higher GTPase, similar to the effect of the known bottom capper protein MciZ. The model could match nucleation kinetics of several flavors of FtsZ using the same parameters as treadmilling and varying only the R to T transition of monomers.

Precise quantification of the energetics and interactions that stabilize membrane proteins in a lipid bilayer is a long-sought goal. Toward this end, atomic force microscopy has been used to unfold individual membrane proteins embedded in their native lipid bilayer, typically by retracting the cantilever at a constant velocity. Recently, unfolding intermediates separated by as few as two amino acids were detected using focused-ion-beam-modified ultrashort cantilevers. However, unambiguously discriminating between such closely spaced states remains challenging, in part because any individual unfolding trajectory only occupies a subset of the total number of intermediates. Moreover, structural assignment of these intermediates via worm-like-chain analysis is hindered by brief dwell times compounded with thermal and instrumental noise. To overcome these issues, we moved the cantilever in a sawtooth pattern of 6–12 nm, offset by 0.25–1 nm per cycle, generating a “zigzag” force ramp of alternating positive and negative loading rates. We applied this protocol to the model membrane protein bacteriorhodopsin (bR). In contrast to conventional studies that extract bR’s photoactive retinal along with the first transmembrane helix, we unfolded bR in the presence of its retinal. To do so, we introduced a previously developed enzymatic-cleavage site between helices E and F and pulled from the top of the E helix using a site-specific, covalent attachment. The resulting zigzag unfolding trajectories occupied 40% more states per trajectory and occupied those states for longer times than traditional constant-velocity records. In total, we identified 31 intermediates during the unfolding of five helices of EF-cleaved bR. These included a previously reported, mechanically robust intermediate located between helices C and B that, with our enhanced resolution, is now shown to be two distinct states separated by three amino acids. Interestingly, another intermediate directly interacted with the retinal, an interaction confirmed by removing the retinal.

The pathology of Plasmodium falciparum malaria is largely defined by the cytoadhesion of infected erythrocytes to the microvascular endothelial lining. The complexity of the endothelial surface and the large range of interactions available for the infected erythrocyte via parasite-encoded adhesins make analysis of critical contributions during cytoadherence challenging to define. Here, we have explored supported membranes functionalized with two important adhesion receptors, ICAM1 or CD36, as a quantitative biomimetic surface to help understand the processes involved in cytoadherence. Parasitized erythrocytes bound to the receptor-functionalized membranes with high efficiency and selectivity under both static and flow conditions, with infected wild-type erythrocytes displaying a higher binding capacity than do parasitized heterozygous sickle cells. We further show that the binding efficiency decreased with increasing intermolecular receptor distance and that the cell-surface contacts were highly dynamic and increased with rising wall shear stress as the cell underwent a shape transition. Computer simulations using a deformable cell model explained the wall-shear-stress-induced dynamic changes in cell shape and contact area via the specific physical properties of erythrocytes, the density of adhesins presenting knobs, and the lateral movement of receptors in the supported membrane.

The coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is responsible for the coronavirus disease 2019 pandemic, and the closely related SARS-CoV coronavirus enter cells by binding at the human angiotensin converting enzyme 2 (hACE2). The stronger hACE2 affinity of SARS-CoV-2 has been connected with its higher infectivity. In this work, we study hACE2 complexes with the receptor-binding domains (RBDs) of the human SARS-CoV-2 and human SARS-CoV viruses, using all-atom molecular dynamics simulations and computational protein design with a physics-based energy function. The molecular dynamics simulations identify charge-modifying substitutions between the CoV-2 and CoV RBDs, which either increase or decrease the hACE2 affinity of the SARS-CoV-2 RBD. The combined effect of these mutations is small, and the relative affinity is mainly determined by substitutions at residues in contact with hACE2. Many of these findings are in line and interpret recent experiments. Our computational protein design calculations redesign positions 455, 493, 494, and 501 of the SARS-CoV-2 receptor binding motif, which contact hACE2 in the complex and are important for ACE2 recognition. Sampling is enhanced by an adaptive importance sampling Monte Carlo method. Sequences with increased affinity replace CoV-2 glutamine by a negative residue at position 493; serine by a nonpolar or aromatic residue or an asparagine at position 494; and asparagine by valine or threonine at position 501. Substitutions at positions 455 and 501 have a smaller effect on affinity. Substitutions suggested by our design are seen in viral sequences encountered in other species, including bat and pangolin. Our results might be used to identify potential virus strains with higher human infectivity and assist in the design of peptide-based or peptidomimetic compounds with the potential to inhibit SARS-CoV-2 binding at hACE2.

The ability of proteins to assemble at sites of high membrane curvature is essential to diverse membrane remodeling processes, including clathrin-mediated endocytosis. Multiple adaptor proteins within the clathrin pathway have been shown to sense regions of high membrane curvature, leading to local recruitment of the clathrin coat. Because clathrin triskelia do not bind to the membrane directly, it has remained unclear whether the clathrin coat plays an active role in sensing membrane curvature or is passively recruited by adaptor proteins. Using a synthetic tag to assemble clathrin directly on membrane surfaces, here we show that clathrin is a strong sensor of membrane curvature, comparable with previously studied adaptor proteins. Interestingly, this sensitivity arises from clathrin assembly rather than from the properties of unassembled triskelia, suggesting that triskelia have preferred angles of interaction, as predicted by earlier structural data. Furthermore, when clathrin is recruited by adaptors, its curvature sensitivity is amplified by 2- to 10-fold, such that the resulting protein complex is up to 100 times more likely to assemble on a highly curved surface compared with a flatter one. This exquisite sensitivity points to a synergistic relationship between the coat and its adaptor proteins, which enables clathrin to pinpoint sites of high membrane curvature, an essential step in ensuring robust membrane traffic. More broadly, these findings suggest that protein networks, rather than individual protein domains, are likely the most potent drivers of membrane curvature sensing.

Proteins carry out a wide range of functions that are tightly regulated in space and time. Protein phosphorylation is the most common post-translation modification of proteins and plays a key role in the regulation of many biological processes. The finding that many phosphorylated residues are not solvent exposed in the unphosphorylated state opens several questions for understanding the mechanism that underlies phosphorylation and how phosphorylation may affect protein structures. First, because kinases need access to the phosphorylated residue, how do such buried residues become modified? Second, once phosphorylated, what are the structural effects of phosphorylation of buried residues, and do they lead to changed conformational dynamics? We have used the ternary complex between p27Kip1 (p27), Cdk2, and cyclin A to study these questions using enhanced sampling molecular dynamics simulations. In line with previous NMR and single-molecule fluorescence experiments, we observe transient exposure of Tyr88 in p27, even in its unphosphorylated state. Once Tyr88 is phosphorylated, we observe a coupling to a second site, thus making Tyr74 more easily exposed and thereby the target for a second phosphorylation step. Our observations provide atomic details on how protein dynamics plays a role in modulating multisite phosphorylation in p27, thus supplementing previous experimental observations. More generally, we discuss how the observed phenomenon of transient exposure of buried residues may play a more general role in regulating protein function.

Zebrafish (Danio rerio) swim within days of fertilization, powered by muscles of the axial myotomes. Forces generated by these muscles can be measured rapidly in whole, intact larval tails by adapting protocols developed forex vivomuscle mechanics. But it is not known how well these measurements reflect the function of the underlying muscle fibers and sarcomeres. Here we consider the anatomy of the 5-day-old, wild-type larval tail, and implement technical modifications to measuring muscle physiology in intact tails. Specifically, we quantify fundamental relationships between force, length, and shortening velocity, and capture the extreme contractile speeds required to swim with tail-beat frequencies of 80-100 Hz. Therefore, we analyze 1000 frames/second movies to track the movement of structures, visible in the transparent tail, which correlate with sarcomere length. We also characterize the passive viscoelastic properties of the preparation to isolate forces contributed by non-muscle structures within the tail. Myotomal muscles generate more than 95% of their maximum isometric stress (76±3 mN/mm2) over the range of muscle lengths usedin vivo. They have rapid twitch kinetics (full width at half-maximum stress: 11±1 msec) and a high twitch to tetanus ratio (0.91±0.05), indicating adaptations for fast excitation-contraction coupling. Although contractile stress is relatively low, myotomal muscles develop high net power (134±20 W/kg at 80 Hz) in cyclical work loop experiments designed to simulate thein vivodynamics of muscle fibers during swimming. When shortening at a constant speed of 7±1 muscle lengths/second, muscles develop 86±2% of isometric stress, while peak instantaneous power during 100Hz work loops occurs at 18±2 muscle lengths/second. These approaches can improve the usefulness of zebrafish as a model system for muscle research by providing a rapid and sensitive functional readout for experimental interventions.Statement of significanceThe zebrafish (Danio rerio) may prove a uniquely efficient model system for characterizing vertebrate muscle physiology. Transparent, drug-permeable larva – each, in essence, a fully functional muscle – can be generated rapidly, inexpensively, and in large numbers. Critically, the zebrafish genome contains homologs of major muscle genes and is highly amenable to manipulation. To reach its potential, reliable (and preferably rapid) means are needed to observe the effects of experimental interventions on larval muscle function. In the present study we show how mechanical measurements made on whole, intact larval tails can provide a readout of fundamental muscle-mechanical properties. Additionally, we show that these muscles are among the fastest ever measured, and therefore worthy of study in their own right.

In-cell NMR spectroscopy is a powerful tool to determine the properties of proteins and nucleic acids within living cells. In-cell NMR can give site-specific measurements of interactions, modifications, and dynamics as well as their modulation by the cellular environment. In-cell NMR requires selective incorporation of heavy isotopes into a protein of interest, either through the introduction of exogenously produced protein to a cell’s interior or the selective overexpression of a protein. We developed conditions to allow the use of Saccharomyces cerevisiae, which was chosen because of its genetic tractability, as a eukaryotic expression system for in-cell NMR. We demonstrate this technique using a fragment of S. cerevisiae Nsp1, an FG Nup. FG Nups are intrinsically disordered proteins containing phenylalanine (F)-glycine (G) repeats and form the selective barrier within the nuclear pore complex. Yeast FG Nups have previously been shown to be maintained in a highly dynamic state within living bacteria as measured by in-cell NMR. Interactions thought to stabilize this dynamic state are also present in the protein’s native organism, although site specificity of interaction is different between the two cytosols.

Interactions between EphB4 receptor tyrosine kinases and their membrane-bound ephrin-B2 ligands on apposed cells play a regulatory role in neural stem cell differentiation. With both receptor and ligand constrained to move within the membranes of their respective cells, this signaling system inevitably experiences spatial confinement and mechanical forces in conjunction with receptor-ligand binding. In this study, we reconstitute the EphB4-ephrin-B2 juxtacrine signaling geometry using a supported-lipid-bilayer system presenting laterally mobile and monomeric ephrin-B2 ligands to live neural stem cells. This experimental platform successfully reconstitutes EphB4-ephrin-B2 binding, lateral clustering, downstream signaling activation, and neuronal differentiation, all in a configuration that preserves the spatiomechanical aspects of the natural juxtacrine signaling geometry. Additionally, the supported bilayer system allows control of lateral movement and clustering of the receptor-ligand complexes through patterns of physical barriers to lateral diffusion fabricated onto the underlying substrate. The results from this study reveal a distinct spatiomechanical effect on the ability of EphB4-ephrin-B2 signaling to induce neuronal differentiation. These observations parallel similar studies of the EphA2-ephrin-A1 system in a very different biological context, suggesting that such spatiomechanical regulation may be a common feature of Eph-ephrin signaling.

The mechanical properties of biological tissues are increasingly recognized as important factors in developmental and pathological processes. Most existing mechanical measurement techniques either necessitate destruction of the tissue for access or provide insufficient spatial resolution. Here, we show for the first time to our knowledge a systematic application of confocal Brillouin microscopy to quantitatively map the mechanical properties of spinal cord tissues during biologically relevant processes in a contact-free and nondestructive manner. Living zebrafish larvae were mechanically imaged in all anatomical planes during development and after spinal cord injury. These experiments revealed that Brillouin microscopy is capable of detecting the mechanical properties of distinct anatomical structures without interfering with the animal's natural development. The Brillouin shift within the spinal cord remained comparable during development and transiently decreased during the repair processes after spinal cord transection. By taking into account the refractive index distribution, we explicitly determined the apparent longitudinal modulus and viscosity of different larval zebrafish tissues. Importantly, mechanical properties differed between tissues in situ and in excised slices. The presented work constitutes the first step toward an in vivo assessment of spinal cord tissue mechanics during regeneration, provides a methodical basis to identify key determinants of mechanical tissue properties, and allows us to test their relative importance in combination with biochemical and genetic factors during developmental and regenerative processes.

Leiomodins (Lmods) are a family of actin filament nucleators related to tropomodulins (Tmods), which are pointed end-capping proteins. Whereas Tmods have alternating tropomyosin- and actin-binding sites (TMBS1, ABS1, TMBS2, ABS2), Lmods lack TMBS2 and half of ABS1, and present a C-terminal extension containing a proline-rich domain and an actin-binding Wiskott-Aldrich syndrome protein homology 2 (WH2) domain that is absent in Tmods. Most of the nucleation activity of Lmods resides within a fragment encompassing ABS2 and the C-terminal extension. This fragment recruits actin monomers into a polymerization nucleus. Here, we revise a recently reported structure of this region of Lmod2 in complex with actin and provide biochemical validation for the newly revised structure. We find that instead of two actin subunits connected by a single Lmod2 polypeptide, as reported in the original structure, the P1 unit cell contains two nearly identical copies of actin monomers, each bound to Lmod2's ABS2 and WH2 domain, with no electron density connecting these two domains. Moreover, we show that the two actin molecules in the unit cell are related to each other by a local twofold noncrystallographic symmetry axis, a conformation clearly distinct from that of actin subunits in the helical filament. We further find that a proposed actin-binding site within the missing connecting region of Lmod2, termed helix h1, does not bind actin in vitro and that the electron density assigned to it in the original structure corresponds instead to a WH2 domain with opposite backbone directionality. Polymerization assays using Lmod2 mutants of helix h1 and the WH2 domain support this conclusion. Finally, we find that deleting the C-terminal extension of Lmod1 and Lmod2 results in an approximately threefold decrease in the nucleation activity, which is only partially accounted for by the lack of the WH2 domain.

The highly charged RNA molecules, with each phosphate carrying a single negative charge, cannot fold into well-defined architectures with tertiary interactions, in the absence of ions. For ribozymes, divalent cations are known to be more efficient than monovalent ions in driving them to a compact state although Mg$^{2+}$ ions are needed for catalytic activities. Therefore, how ions interact with RNA is relevant in understanding RNA folding. It is often thought that most of the ions are territorially and non-specifically bound to the RNA, as predicted by the counterion condensation (CIC) theory. Here, we show using simulations of ${\it Azoarcus}$ ribozyme, based on an accurate coarse-grained Three Site Interaction (TIS) model, with explicit divalent and monovalent cations, that ion condensation is highly specific and depends on the nucleotide position. The regions with high coordination between the phosphate groups and the divalent cations are discernible even at very low Mg$^{2+}$ concentrations when the ribozyme does not form tertiary interactions. Surprisingly, these regions also contain the secondary structural elements that nucleate subsequently in the self-assembly of RNA, implying that ion condensation is determined by the architecture of the folded state. These results are in sharp contrast to interactions of ions (monovalent and divalent) with rigid charged rods in which ion condensation is uniform and position independent. The differences are explained in terms of the dramatic non-monotonic shape fluctuations in the ribozyme as it folds with increasing Mg$^{2+}$ or Ca$^{2+}$ concentration., Comment: 22 pages including 9 figures, 5 SI figures, and 1 SI table

Escherichia coli single-strand (ss) DNA-binding protein (SSB) is an essential protein that binds ssDNA intermediates formed during genome maintenance. SSB homotetramers bind ssDNA in several modes differing in occluded site size and cooperativity. The 35-site-size ((SSB)35) mode favored at low [NaCl] and high SSB/DNA ratios displays high “unlimited” nearest-neighbor cooperativity (ω35), forming long protein clusters, whereas the 65-site-size ((SSB)65) mode in which ssDNA wraps completely around the tetramer is favored at higher [NaCl] (>200 mM) and displays “limited” cooperativity (ω65), forming only dimers of tetramers. In addition, a non-nearest-neighbor high cooperativity can also occur in the (SSB)65 mode on long ssDNA even at physiological salt concentrations in the presence of glutamate and requires its intrinsically disordered C-terminal linker (IDL) region. However, whether cooperativity exists between the different modes and the role of the IDL in nearest-neighbor cooperativity has not been probed. Here, we combine sedimentation velocity and fluorescence titration studies to examine nearest-neighbor cooperativity in each binding mode and between binding modes using (dT)70 and (dT)140. We find that the (SSB)35 mode always shows extremely high “unlimited” cooperativity that requires the IDL. At high salt, wild-type SSB and a variant without the IDL, SSB-ΔL, bind in the (SSB)65 mode but show little cooperativity, although cooperativity increases at lower [NaCl] for wild-type SSB. We also find significant intermode nearest-neighbor cooperativity (ω65/35), with ω65 ≪ ω65/35

The arrangement of receptors in the plasma membrane strongly affects the ability of a cell to sense its environment both in terms of sensitivity and in terms of spatial resolution. The spatial and temporal arrangement of the receptors is affected in turn by the mechanical properties and the structure of the cell membrane. Here we focus on characterizing the flow of the membrane in response to the motion of a protein embedded in it. We do so by measuring the correlated diffusion of extracellularly tagged transmembrane neurotrophin receptors TrkB and p75 on transfected neuronal cells. In accord with previous reports, we find that the motion of single receptors exhibits transient confinement to sub-micron domains. We confirm predictions based on hydrodynamics of fluid membranes, finding long-range correlations in the motion of the receptors in the plasma membrane. However, we discover that these correlations do not persist for long ranges, as predicted, but decay exponentially, with a typical decay length on the scale of the average confining domain size.

The TMEM16 family comprises Ca 2+ -activated Cl − channels and phospholipid scramblases. The crystal structure of a fungal homolog, nhTMEM16, revealed an important architectural feature of this protein family in the form of a bilayer-spanning hydrophilic groove that is directly exposed to the membrane. This groove likely provides a pathway for lipid translocation. As mutations that alter ion channel activity of the TMEM16 proteins localize around the groove, it was suggested that the ion and lipid pathways coincide such that the ion pore is partly lined by phospholipids. However, this proposal was not supported by the observation that nhTMEM16 does not mediate ion transport. Here we show that nhTMEM16 mediates both ion and lipid transport and that its properties closely resemble those of a previously characterized fungal homolog, afTMEM16. We show that the reported lack of ion transport activity of nhTMEM16 is due to the lipid composition of the reconstitution membranes and to the presence of a GFP tag. Thus, nhTMEM16, like afTMEM16 and the mammalian TMEM16F, mediates simultaneous lipid scrambling and nonspecific ion transport. This supports the hypothesis that these two processes are tightly correlated and likely to be a general functional feature of the TMEM16 scramblases and therefore of general importance in understanding their biological roles.

Hydra, a simple freshwater animal famous for its regenerative capabilities, must tear a hole through its epithelial tissue each time it opens its mouth. The feeding response of Hydra has been well-characterized physiologically and is regarded as a classical model system for environmental chemical biology. However, due to a lack of in vivo labeling and imaging tools, the biomechanics of mouth opening have remained completely unexplored. We take advantage of the availability of transgenic Hydra lines to perform the first dynamical analysis, to our knowledge, of Hydra mouth opening and test existing hypotheses regarding the underlying cellular mechanisms. Through cell position and shape tracking, we show that mouth opening is accompanied by changes in cell shape, but not cellular rearrangements as previously suggested. Treatment with a muscle relaxant impairs mouth opening, supporting the hypothesis that mouth opening is an active process driven by radial contractile processes (myonemes) in the ectoderm. Furthermore, we find that all events exhibit the same relative rate of opening. Because one individual can open consecutively to different amounts, this suggests that the degree of mouth opening is controlled through neuronal signaling. Finally, from the opening dynamics and independent measurements of the elastic properties of the tissues, we estimate the forces exerted by the myonemes to be on the order of a few nanoNewtons. Our study provides the first dynamical framework, to our knowledge, for understanding the remarkable plasticity of the Hydra mouth and illustrates that Hydra is a powerful system for quantitative biomechanical studies of cell and tissue behaviors in vivo.

The human protein TDP-43 is a major component of the cellular aggregates found in amyotrophic lateral sclerosis and other neurodegenerative diseases. Insoluble cytoplasmic aggregates isolated from the brain of amyotrophic lateral sclerosis and frontotemporal lobar degeneration patients contain ubiquitinated, hyperphosphorylated, and N-terminally truncated TDP-43. Truncated fragments of TDP-43 identified from patient tissues contain part of the second RNA recognition motif (RRM2) and the disordered C-terminus, indicating that both domains can be involved in aggregation and toxicity. Here, we focus on RRM2. Using all-atom replica-averaged metadynamics simulations with NMR chemical shift restraints, we characterized the atomic structure of non-native states of RRM2, sparsely populated under native conditions. These structures reveal the exposure to the solvent of aggregation-prone peptide regions, normally buried in the native state, supporting a role in aggregation for the partially folded states of RRM2.

Cell sorting, whereby a heterogeneous cell mixture organizes into distinct tissues, is a fundamental patterning process in development. Hydra is a powerful model system for carrying out studies of cell sorting in three dimensions, because of its unique ability to regenerate after complete dissociation into individual cells. The physicists Alfred Gierer and Hans Meinhardt recognized Hydra's self-organizing properties more than 40 years ago. However, what drives cell sorting during regeneration of Hydra from cell aggregates is still debated. Differential motility and differential adhesion have been proposed as driving mechanisms, but the available experimental data are insufficient to distinguish between these two. Here, we answer this longstanding question by using transgenic Hydra expressing fluorescent proteins and a multiscale experimental and numerical approach. By quantifying the kinematics of single cell and whole aggregate behaviors, we show that no differences in cell motility exist among cell types and that sorting dynamics follow a power law with an exponent of ∼0.5. Additionally, we measure the physical properties of separated tissues and quantify their viscosities and surface tensions. Based on our experimental results and numerical simulations, we conclude that tissue interfacial tensions are sufficient to explain cell sorting in aggregates of Hydra cells. Furthermore, we demonstrate that the aggregate's geometry during sorting is key to understanding the sorting dynamics and explains the exponent of the power law behavior. Our results answer the long standing question of the physical mechanisms driving cell sorting in Hydra cell aggregates. In addition, they demonstrate how powerful this organism is for biophysical studies of self-organization and pattern formation.

Intrinsically disordered proteins and intrinsically disordered regions are implicated in many biological functions and associated with many diseases, but their conformational characterizations are challenging. The disordered β6/β7 loop of Staphylococcus aureus sortase A is involved in the binding of both sorting signals and calcium. Calcium binding allosterically activates the enzyme, but the detailed mechanism has been unclear. Here we adapted the replica exchange with solute tempering method to sample the conformations of the β6/β7 loop, in apo form and in three liganded forms (bound with a sorting signal or calcium or both). The extensive molecular dynamics simulations yield atomic details of the disordered-to-order transition of the loop and suggest a mechanism for allosteric activation: calcium binding results in partial closure and ordering of the loop, thereby leading to preorganization of the binding pocket for the sorting signal. The approach has general applicability to the study of intrinsically disordered regions.

Magic-angle spinning nuclear magnetic resonance is well suited for the study of membrane proteins in the nativelike lipid environment. However, the natural cellular membrane is invariably more complex than the proteoliposomes most often used for solid-state NMR (SSNMR) studies, and differences may affect the structure and dynamics of the proteins under examination. In this work we use SSNMR and other biochemical and biophysical methods to probe the structure of a seven-transmembrane helical photoreceptor, Anabaena sensory rhodopsin (ASR), prepared in the Escherichia coli inner membrane, and compare it to that in a bilayer formed by DMPC/DMPA lipids. We find that ASR is organized into trimers in both environments but forms two-dimensional crystal lattices of different symmetries. It favors hexagonal packing in liposomes, but may form a square lattice in the E. coli membrane. To examine possible changes in structure site-specifically, we perform two- and three-dimensional SSNMR experiments and analyze the differences in chemical shifts and peak intensities. Overall, this analysis reveals that the structure of ASR is largely conserved in the inner membrane of E. coli, with many of the important structural features of rhodopsins previously observed in ASR in proteoliposomes being preserved. Small, site-specific perturbations in protein structure that occur as a result of the membrane changes indicate that the protein can subtly adapt to its environment without large structural rearrangement.

Fast amoeboid migration requires cells to apply mechanical forces on their surroundings via transient adhesions. However, the role these forces play in controlling cell migration speed remains largely unknown. We used three-dimensional force microscopy to measure the three-dimensional forces exerted by chemotaxing Dictyostelium cells, and examined wild-type cells as well as mutants with defects in contractility, internal F-actin crosslinking, and cortical integrity. We showed that cells pull on their substrate adhesions using two distinct, yet interconnected mechanisms: axial actomyosin contractility and cortical tension. We found that the migration speed increases when axial contractility overcomes cortical tension to produce the cell shape changes needed for locomotion. We demonstrated that the three-dimensional pulling forces generated by both mechanisms are internally balanced by an increase in cytoplasmic pressure that allows cells to push on their substrate without adhering to it, and which may be relevant for amoeboid migration in complex three-dimensional environments.

The CaV2.2 (N-type) and CaV2.1 (P/Q-type) voltage-dependent calcium channels are prevalent throughout the nervous system where they mediate synaptic transmission, but the basis for the selective presence at individual synapses still remains an open question. The CaV2.1 channels have been proposed to respond more effectively to brief action potentials (APs), an idea supported by computational modeling. However, the side-by-side comparison of CaV2.1 and CaV2.2 kinetics in intact neurons failed to reveal differences. As an alternative means for direct functional comparison we expressed zebrafish CaV2.1 and CaV2.2 α-subunits, along with their accessory subunits, in HEK293 cells. HEK cells lack calcium currents, thereby circumventing the need for pharmacological inhibition of mixed calcium channel isoforms present in neurons. HEK cells also have a simplified morphology compared to neurons, which improves voltage control. Our measurements revealed faster kinetics and shallower voltage-dependence of activation and deactivation for CaV2.1. Additionally, recordings of calcium current in response to a command waveform based on the motorneuron AP show, directly, more effective activation of CaV2.1. Analysis of calcium currents associated with the AP waveform indicate an approximately fourfold greater open probability (PO) for CaV2.1. The efficient activation of CaV2.1 channels during APs may contribute to the highly reliable transmission at zebrafish neuromuscular junctions.

Mechanotransduction pathways convert forces that stress and strain structures within cells into gene expression levels that impact development, homeostasis, and disease. The levels of some key structural proteins in the nucleus, cytoskeleton, or extracellular matrix have been recently reported to scale with tissue- and cell-level forces or mechanical properties such as stiffness, and so the mathematics of mechanotransduction becomes important to understand. Here, we show that if a given structural protein positively regulates its own gene expression, then stresses need only inhibit degradation of that protein to achieve stable, mechanosensitive gene expression. This basic use-it-or-lose-it module is illustrated by application to meshworks of nuclear lamin A, minifilaments of myosin II, and extracellular matrix collagen fibers—all of which possess filamentous coiled-coil/supercoiled structures. Past experiments not only suggest that tension suppresses protein degradation mediated and/or initiated by various enzymes but also that transcript levels vary with protein levels because key transcription factors are regulated by these structural proteins. Coupling between modules occurs within single cells and between cells in tissue, as illustrated during embryonic heart development where cardiac fibroblasts make collagen that cardiomyocytes contract. With few additional assumptions, the basic module has sufficient physics to control key structural genes in both development and disease.

Voltage-dependence of Na(+)-coupled phosphate cotransporters of the SLC34 family arises from displacement of charges intrinsic to the protein and the binding/release of one Na(+) ion in response to changes in the transmembrane electric field. Candidate coordination residues for the cation at the Na1 site were previously predicted by structural modeling using the x-ray structure of dicarboxylate transporter VcINDY as template and confirmed by functional studies. Mutations at Na1 resulted in altered steady-state and presteady-state characteristics that should be mirrored in the conformational changes induced by membrane potential changes. To test this hypothesis by functional analysis, double mutants of the flounder SLC34A2 protein were constructed that contain one of the Na1-site perturbing mutations together with a substituted cysteine for fluorophore labeling, as expressed in Xenopus oocytes. The locations of the mutations were mapped onto a homology model of the flounder protein. The effects of the mutagenesis were characterized by steady-state, presteady-state, and fluorometric assays. Changes in fluorescence intensity (ΔF) in response to membrane potential steps were resolved at three previously identified positions. These fluorescence data corroborated the altered presteady-state kinetics upon perturbation of Na1, and furthermore indicated concomitant changes in the microenvironment of the respective fluorophores, as evidenced by changes in the voltage dependence and time course of ΔF. Moreover, iodide quenching experiments indicated that the aqueous nature of the fluorophore microenvironment depended on the membrane potential. These findings provide compelling evidence that membrane potential and cation interactions induce significant large-scale structural rearrangements of the protein.

The molecular cascade that controls switching of the direction of rotation of Escherichia coli flagellar motors is well known, but the conformational changes that allow the rotor to switch are still unclear. The signaling molecule CheY, when phosphorylated, binds to the C-ring at the base of the rotor, raising the probability that the motor spins clockwise. When the concentration of CheY-P is so low that the motor rotates exclusively counterclockwise (CCW), the C-ring recruits more monomers of FliM and tetramers of FliN, the proteins to which CheY-P binds, thus increasing the motor's sensitivity to CheY-P and allowing it to switch once again. Motors that rotate exclusively CCW have more FliM and FliN subunits in their C-rings than motors that rotate exclusively clockwise. How are the new subunits accommodated? Does the diameter of the C-ring increase, or do FliM and FliN get packed in a different pattern, keeping the overall diameter of the C-ring constant? Here, by measuring fluorescence anisotropy of yellow fluorescent protein-labeled motors, we show that the CCW C-rings accommodate more FliM monomers without changing the spacing between them, and more FliN monomers at the same time as increasing their effective spacing and/or changing their orientation within the tetrameric structure.

Extracellular vesicles (EVs) are proposed to play important roles in intercellular communication. Two classes of EVs can be distinguished based on their intracellular origin. Exosomes are generated within endosomes and released when these fuse with the plasma membrane, whereas ectosomes bud directly from the plasma membrane. Studies of EV function have been hindered by limited understanding of their biogenesis. Components of the endosomal sorting complex required for transport (ESCRT) machinery play essential roles in topologically equivalent processes at both the endosome and the plasma membrane and are consistently recovered in EVs, but whether they are generally required to produce EVs is still debated. Here, we study the effects of inhibiting the ESCRT-associated AAA+ ATPase VPS4 on EV release from cultured cells using two methods for EV recovery, differential centrifugation and polyethylene glycol precipitation followed by lectin affinity chromatography. We find that inhibiting VPS4 in HEK293 cells decreases release of EV-associated proteins and miRNA as well as the overall number of EV particles. The tetraspanins CD63 and CD9 are among the most frequently monitored EV proteins, but they differ in their subcellular localization, with CD63 primarily in endosomes and CD9 on the plasma membrane. We find that CD63 and CD9 are enriched in separable populations of EVs that are both sensitive to VPS4 inhibition. Serum stimulation increases release of both types of EVs and is also reduced by inhibiting VPS4. Taken together, our data indicate that VPS4 activity is important for generating exosomes and ectosomes, thereby generally implicating the ESCRT machinery in EV biogenesis.

We previously used an orthogonal high-throughput screen to select peptides that spontaneously cross synthetic lipid bilayers without bilayer disruption. Many of the 12-residue spontaneous membrane translocating peptides (SMTPs) selected from the library contained a 5-residue consensus motif, LRLLR in positions 5-9. We hypothesized that the conserved motif could be a necessary and sufficient minimal motif for translocation. To test this and to explore the mechanism of spontaneous membrane translocation, we synthesized seven arginine placement variants of LRLLRWC and compared their membrane partitioning, translocation, and perturbation to one of the parent SMTPs, called "TP2". Several motif variant peptides translocate into synthetic vesicles with rates that are similar to TP2. However, the peptide containing the selected motif, LRLLRWC, was not the fastest; sequence context is also important for translocation efficiency. Although none of these peptides permeabilize bilayers, the motif peptides translocate faster at higher peptide to lipid ratios, suggesting that bilayer perturbation and/or cooperative interactions are important for their translocation. On the other hand, TP2 translocates slower as its concentration is increased, suggesting that TP2 translocates as a monomer and is inhibited by lateral interactions in the membrane. TP2 and the LRLLR motif peptide induce lipid translocation, suggesting that lipids chaperone them across the bilayer. The other motif peptides do not induce lipid flip-flop, suggesting an alternate mechanism. Concatenated motifs translocate slower than the motifs alone. Variants of TP2 with shorter and longer arginine side-chain analogs translocate slower than TP2. In summary, these results suggest that multiple patterns of leucine and arginine can support spontaneous membrane translocation, and that sequence context is important for the contribution of the motifs. Because motifs do not make simple, additive contributions to spontaneous translocation, rational engineering of novel SMTPs will remain difficult, providing even more reason to pursue SMTP discovery with synthetic molecular evolution.

Spatiotemporal control of the cAMP signaling pathway is governed by both hormonal stimulation of cAMP generation by adenylyl cyclases (activation phase) and cAMP hydrolysis by phosphodiesterases (PDEs) (termination phase). The termination phase is initiated by PDEs actively targeting the protein kinase A (PKA) R-subunit through formation of a PDE-PKAR-cyclic adenosine monophosphate (cAMP) complex (the termination complex). Our results using PDE8 as a model PDE, reveal that PDEs mediate active hydrolysis of cAMP bound to its receptor RIα by enhancing the enzymatic activity. This accelerated cAMP turnover occurs via formation of a stable PDE8-RIα complex, where the protein-protein interface forms peripheral contacts and the central ligand cements this ternary interaction. The basis for enhanced catalysis is active translocation of cAMP from its binding site on RIα to the hydrolysis site on PDE8 through direct “channeling.” Our results reveal cAMP channeling in the PDE8-RIα complex and a molecular description of how this channel facilitates processive hydrolysis of unbound cAMP. Thus, unbound cAMP maintains the PDE8-RIα complex while being hydrolyzed, revealing an undiscovered mode for amplification of PKA activity by cAMP-mediated sequestration of the R-subunit by PDEs. This novel regulatory mode explains the paradox of cAMP signal amplification by accelerated PDE-mediated cAMP turnover. This highlights how target effector proteins of small-molecule ligands can promote enzyme-mediated ligand hydrolysis by scaffolding effects. Enhanced activity of the PDE8-RIα complex facilitates robust desensitization, allowing the cell to respond to dynamic levels of cAMP rather than steady-state levels. The PDE8-RIα complex represents a new class of PDE-based complexes for specific drug discovery targeting the cAMP signaling pathway.

Mechanical forces are critical but poorly understood inputs for organogenesis and wound healing. Calcium ions (Ca2+) are critical second messengers in cells for integrating environmental and mechanical cues, but the regulation of Ca2+ signaling is poorly understood in developing epithelial tissues. Here we report a chip-based regulated environment for microorgans that enables systematic investigations of the crosstalk between an organ’s mechanical stress environment and biochemical signaling under genetic and chemical perturbations. This method enabled us to define the essential conditions for generating organ-scale intercellular Ca2+ waves in Drosophila wing discs that are also observed in vivo during organ development. We discovered that mechanically induced intercellular Ca2+ waves require fly extract growth serum as a chemical stimulus. Using the chip-based regulated environment for microorgans, we demonstrate that not the initial application but instead the release of mechanical loading is sufficient, but not necessary, to initiate intercellular Ca2+ waves. The Ca2+ response depends on the prestress intercellular Ca2+ activity and not on the magnitude or duration of the mechanical stimulation applied. Mechanically induced intercellular Ca2+ waves rely on IP3R-mediated Ca2+-induced Ca2+ release and propagation through gap junctions. Thus, intercellular Ca2+ waves in developing epithelia may be a consequence of stress dissipation during organ growth.

The lack of physiological recordings from Caenorhabditis elegans embryos stands in stark contrast to the comprehensive anatomical and gene expression datasets already available. Using light-sheet fluorescence microscopy to address the challenges associated with functional imaging at this developmental stage, we recorded calcium dynamics in muscles and neurons and developed analysis strategies to relate activity and movement. In muscles, we found that the initiation of twitching was associated with a spreading calcium wave in a dorsal muscle bundle. Correlated activity in muscle bundles was linked with early twitching and eventual coordinated movement. To identify neuronal correlates of behavior, we monitored brainwide activity with subcellular resolution and identified a particularly active cell associated with muscle contractions. Finally, imaging neurons of a well-defined adult motor circuit, we found that reversals in the eggshell correlated with calcium transients in AVA interneurons.

The ability of cells to sense and respond to mechanical cues from the surrounding environment has been implicated as a key regulator of cell differentiation, migration, and proliferation. The extracellular matrix (ECM) is an oft-overlooked component of the interface between cells and their surroundings. Cells assemble soluble ECM proteins into insoluble fibrils with unique mechanical properties that can alter the mechanical cues a cell receives. In this study, we construct a model that predicts the dynamics of cellular traction force generation and subsequent assembly of fibrils of the ECM protein fibronectin (FN). FN fibrils are the primary component in primordial ECM and, as such, FN assembly is a critical component in the cellular mechanical response. The model consists of a network of Hookean springs, each representing an extensible domain within an assembling FN fibril. As actomyosin forces stretch the spring network, simulations predict the resulting traction force and FN fibril formation. The model accurately predicts FN fibril morphometry and demonstrates a mechanism by which FN fibril assembly regulates traction force dynamics in response to mechanical stimuli and varying surrounding substrate stiffness.

Cells can alter the lipid content of their plasma membranes upon changes in their environment to maintain and adjust membrane function. Recent work suggests that some membrane functions arise because cellular plasma membranes are poised close to a miscibility transition under growth conditions. Here we report experiments utilizing giant plasma membrane vesicles (GPMVs) to explore how membrane transition temperature varies with growth temperature in a zebrafish cell line (ZF4) that can be adapted for growth between 20 and 32°C. We find that GPMV transition temperatures adjust to be 16.7 ± 1.2°C below growth temperature for four growth temperatures investigated and that adjustment occurs over roughly 2 days when temperature is abruptly lowered from 28 to 20°C. We also find that GPMVs have slightly different lipidomes when isolated from cells adapted for growth at 28 and 20°C. Similar to past work in vesicles derived from mammalian cells, fluctuating domains are observed in ZF4-derived GPMVs, consistent with their having critical membrane compositions. Taken together, these experimental results suggest that cells in culture biologically tune their membrane composition in a way that maintains specific proximity to a critical miscibility transition.

To study the impact of neural activity on cellular physiology, one would like to combine precise control of firing patterns with highly sensitive probes of cellular physiology. Light-gated ion channels, e.g., Channelrhodopsin-2, enable precise control of firing patterns; green fluorescent protein-based reporters, e.g., the GCaMP6f Ca2+ reporter, enable highly sensitive probing of cellular physiology. However, for most actuator-reporter combinations, spectral overlap prevents straightforward combination within a single cell. Here we explore multiwavelength control of channelrhodopsins to circumvent this limitation. The “stoplight” technique described in this article uses channelrhodopsin variants that are opened by blue light and closed by orange light. Cells are illuminated with constant blue light to excite fluorescence of a green fluorescent protein-based reporter. Modulated illumination with orange light negatively regulates activation of the channelrhodopsin. We performed detailed photophysical characterization and kinetic modeling of four candidate stoplight channelrhodopsins. The variant with the highest contrast, sdChR(C138S,E154A), enabled all-optical measurements of activity-induced calcium transients in cultured rat hippocampal neurons, although cell-to-cell variation in expression levels presents a challenge for quantification.

Soft x-ray tomography (SXT) is increasingly being recognized as a valuable method for visualizing and quantifying the ultrastructure of cryopreserved cells. Here, we describe the combination of SXT with cryogenic confocal fluorescence tomography (CFT). This correlative approach allows the incorporation of molecular localization data, with isotropic precision, into high-resolution three-dimensional (3-D) SXT reconstructions of the cell. CFT data are acquired first using a cryogenically adapted confocal light microscope in which the specimen is coupled to a high numerical aperture objective lens by an immersion fluid. The specimen is then cryo-transferred to a soft x-ray microscope (SXM) for SXT data acquisition. Fiducial markers visible in both types of data act as common landmarks, enabling accurate coalignment of the two complementary tomographic reconstructions. We used this method to identify the inactive X chromosome (Xi) in female v-abl transformed thymic lymphoma cells by localizing enhanced green fluorescent protein-labeled macroH2A with CFT. The molecular localization data were used to guide segmentation of Xi in the SXT reconstructions, allowing characterization of the Xi topological arrangement in near-native state cells. Xi was seen to adopt a number of different topologies with no particular arrangement being dominant.

Fluorescently labeled proteins that are found both in the cytoplasm and at the plasma membrane, such as peripheral membrane proteins, create stratified fluorescent layers that present a challenging environment for brightness studies with fluorescence fluctuation spectroscopy. The geometry of each layer along with fluorescence and brightness contributions from adjacent layers generates a convoluted raw brightness that conceals the underlying brightness of each individual layer. Because the brightness at a layer establishes the oligomeric state of the fluorescently labeled protein at said layer, we developed a method that connects the experimental raw brightness with the physical brightness at each layered compartment. The technique determines the oligomerization in each compartment from an axial intensity scan through the sample, followed by a fluorescence fluctuation spectroscopy measurement at each layer. We experimentally verify the technique with H-Ras-EGFP as a model system and determine its oligomeric state at both the plasma membrane and in the cytoplasm. Furthermore, we study the oligomerization of the Gag matrix domain of Human T-lymphotropic virus Type 1. The matrix domain targets the Gag polyprotein to the plasma membrane where, subsequently, viral assembly occurs. We determine the oligomerization of matrix in the cytoplasm and observe the onset of protein-protein interactions at the membrane. These observations shed light on the early assembly steps of the retrovirus.

Sodium-glucose transporters (SGLTs) facilitate the movement of water across the cell membrane, playing a central role in cellular homeostasis. Here, we present a detailed analysis of the mechanism of water permeation through the inward-facing state of vSGLT based on nearly 10 μs of molecular dynamics simulations. These simulations reveal the transient formation of a continuous water channel through the transporter that permits water to permeate the protein. Trajectories in which spontaneous release of galactose is observed, as well as those in which galactose remains in the binding site, show that the permeation rate, although modulated by substrate occupancy, is not tightly coupled to substrate release. Using a, to our knowledge, novel channel-detection algorithm, we identify the key residues that control water flow through the transporter and show that solvent gating is regulated by side-chain motions in a small number of residues on the extracellular face. A sequence alignment reveals the presence of two insertion sites in mammalian SGLTs that flank these outer-gate residues. We hypothesize that the absence of these sites in vSGLT may account for the high water permeability values for vSGLT determined via simulation compared to the lower experimental estimates for mammalian SGLT1.

Plasma membrane voltage is a fundamentally important property of a living cell; its value is tightly coupled to membrane transport, the dynamics of transmembrane proteins, and to intercellular communication. Accurate measurement of the membrane voltage could elucidate subtle changes in cellular physiology, but existing genetically encoded fluorescent voltage reporters are better at reporting relative changes than absolute numbers. We developed an Archaerhodopsin-based fluorescent voltage sensor whose time-domain response to a stepwise change in illumination encodes the absolute membrane voltage. We validated this sensor in human embryonic kidney cells. Measurements were robust to variation in imaging parameters and in gene expression levels, and reported voltage with an absolute accuracy of 10 mV. With further improvements in membrane trafficking and signal amplitude, time-domain encoding of absolute voltage could be applied to investigate many important and previously intractable bioelectric phenomena.

© 2017 Biophysical Society Cells employing amoeboid motility exhibit repetitive cycles of rapid expansion and contraction and apply coordinated traction forces to their environment. Although aspects of this process are well studied, it is unclear how the cell controls the coordination of cell length changes with adhesion to the surface. Here, we develop a simple model to mechanistically explain the emergence of periodic changes in length and spatiotemporal dynamics of traction forces measured in chemotaxing unicellular amoeba, Dictyostelium discoideum. In contrast to the biochemical mechanisms that have been implicated in the coordination of some cellular processes, we show that many features of amoeboid locomotion emerge from a simple mechanochemical model. The mechanism for interaction with the environment in Dictyostelium is unknown and thus, we explore different cell-environment interaction models to reveal that mechanosensitive adhesions are necessary to reproduce the spatiotemporal adhesion patterns. In this modeling framework, we find that the other motility modes, such as smooth gliding, arise naturally with variations in the physical properties of the surface. Thus, our work highlights the prominent role of biomechanics in determining the emergent features of amoeboid locomotion.

Myofilament extensibility is a key structural parameter for interpreting myosin cross-bridge kinetics in striated muscle. Previous studies reported much higher thick-filament extensibility at low tension than the better-known and commonly used values at high tension, but in interpreting mechanical studies of muscle, a single value for thick-filament extensibility has usually been assumed. Here, we established the complete thick-filament force-extension curve from actively contracting, intact vertebrate skeletal muscle. To access a wide range of tetanic forces, the myosin inhibitor blebbistatin was used to induce low tetanic forces in addition to the higher tensions obtained from tetanic contractions of the untreated muscle. We show that the force/extensibility curve of the thick filament is nonlinear, so assuming a single value for thick-filament extensibility at all force levels is not justified. We also show that independent of whether tension is generated passively by sarcomere stretch or actively by cross-bridges, the thick-filament extensibility is nonlinear. Myosin head periodicity, however, only changes when active tension is generated under calcium-activated conditions. The nonlinear thick-filament force-extension curve in skeletal muscle, therefore, reflects a purely passive response to either titin-based force or actomyosin-based force, and it does not include a thick-filament activation mechanism. In contrast, the transition of myosin head periodicity to an active configuration appears to only occur in response to increased active force when calcium is present.

Motors that move DNA, or that move along DNA, play essential roles in DNA replication, transcription, recombination, and chromosome segregation. The mechanisms by which these DNA translocases operate remain largely unknown. Some double-stranded DNA (dsDNA) viruses use an ATP-dependent motor to drive DNA into preformed capsids. These include several human pathogens as well as dsDNA bacteriophages—viruses that infect bacteria. We previously proposed that DNA is not a passive substrate of bacteriophage packaging motors but is instead an active component of the machinery. We carried out computational studies on dsDNA in the channels of viral portal proteins, and they reveal DNA conformational changes consistent with that hypothesis. dsDNA becomes longer (“stretched”) in regions of high negative electrostatic potential and shorter (“scrunched”) in regions of high positive potential. These results suggest a mechanism that electrostatically couples the energy released by ATP hydrolysis to DNA translocation: The chemical cycle of ATP binding, hydrolysis, and product release drives a cycle of protein conformational changes. This produces changes in the electrostatic potential in the channel through the portal, and these drive cyclic changes in the length of dsDNA as the phosphate groups respond to the protein’s electrostatic potential. The DNA motions are captured by a coordinated protein-DNA grip-and-release cycle to produce DNA translocation. In short, the ATPase, portal, and dsDNA work synergistically to promote genome packaging.

Superresolution fluorescence microscopy is used to locate single copies of RNA polymerase (RNAP) in live Escherichia coli and track their diffusive motion. On a timescale of 0.1–1 s, most copies separate remarkably cleanly into two diffusive states. The “slow” RNAPs, which move indistinguishably from DNA loci, are assigned to specifically bound copies (with fractional population ftrxn) that are initiating transcription, elongating, pausing, or awaiting termination. The “mixed-state” RNAP copies, with effective diffusion constant Dmixed = 0.21 μm2 s−1, are assigned as a rapidly exchanging mixture of nonspecifically bound copies (fns) and copies undergoing free, three-dimensional diffusion within the nucleoids (ffree). Longer trajectories of 7-s duration reveal transitions between the slow and mixed states, corroborating the assignments. Short trajectories of 20-ms duration enable direct observation of the freely diffusing RNAP copies, yielding Dfree = 0.7 μm2 s−1. Analysis of single-particle trajectories provides quantitative estimates of the partitioning of RNAP into different states of activity: ftrxn = 0.54 ± 0.07, fns = 0.28 ± 0.05, ffree = 0.12 ± 0.03, and fnb = 0.06 ± 0.05 (fraction unable to bind to DNA on a 1-s timescale). These fractions disagree with earlier estimates.

In 1965, Monod et al. (1) published a plausible model on the nature and role of allosteric transitions in proteins. That article triggered the intense interest of a vast number of talented scientists starting with Perutz (2), the father of structural biology. The competition with the induced-fit model of Koshland et al. (3), started an exciting scientific debate. In science, the general validity of an idea can best be appreciated in the long run. After half a century, allostery is still alive and flourishing; and interesting, novel applications of the theory have led to original pharmacological approaches for the therapy of central-nervous system diseases and cancer.

Autonomous circadian oscillations arise from transcriptional-translational feedback loops of core clock components. The period of a circadian oscillator is relatively insensitive to changes in nutrients (e.g., glucose), which is referred to as "nutrient compensation". Recently, a transcription repressor, CSP-1, was identified as a component of the circadian system in Neurospora crassa. The transcription of csp-1 is under the circadian regulation. Intriguingly, CSP-1 represses the circadian transcription factor, WC-1, forming a negative feedback loop that can influence the core oscillator. This feedback mechanism is suggested to maintain the circadian period in a wide range of glucose concentrations. In this report, we constructed a mathematical model of the Neurospora circadian clock incorporating the above WC-1/CSP-1 feedback loop, and investigated molecular mechanisms of glucose compensation. Our model shows that glucose compensation exists within a narrow range of parameter space where the activation rates of csp-1 and wc-1 are balanced with each other, and simulates loss of glucose compensation in csp-1 mutants. More importantly, we experimentally validated rhythmic oscillations of the wc-1 gene expression and loss of glucose compensation in the wc-1(ov) mutant as predicted in the model. Furthermore, our stochastic simulations demonstrate that the CSP-1-dependent negative feedback loop functions in glucose compensation, but does not enhance the overall robustness of oscillations against molecular noise. Our work highlights predictive modeling of circadian clock machinery and experimental validations employing Neurospora and brings a deeper understanding of molecular mechanisms of glucose compensation.

Efficient and accurate mapping of transition pathways is a challenging problem in allosteric proteins. We propose here a to our knowledge new methodology called collective molecular dynamics (coMD). coMD takes advantage of the collective modes of motions encoded by the fold, simultaneously evaluating the interactions and energetics via a full-atomic MD simulation protocol. The basic approach is to deform the structure collectively along the modes predicted by the anisotropic network model, upon selecting them via a Monte Carlo/Metropolis algorithm from among the complete pool of all accessible modes. Application to adenylate kinase, an allosteric enzyme composed of three domains, CORE, LID, and NMP, shows that both open-to-closed and closed-to-open transitions are readily sampled by coMD, with large-scale motions of the LID dominating. An energy-barrier crossing occurs during the NMP movements. The energy barrier originates from a switch between the salt bridges K136-D118 at the LID-CORE interface and K57-E170 and D33-R156 at the CORE-NMP and LID-NMP interfaces, respectively. Despite its simplicity and computing efficiency, coMD yields ensembles of transition pathways in close accord with detailed full atomic simulations, lending support to its utility as a multiscale hybrid method for efficiently exploring the allosteric transitions of multidomain or multimeric proteins.

The cytoskeletal protein MreB is an essential component of the bacterial cell-shape generation system. Using a superresolution variant of total internal reflection microscopy with structured illumination, as well as three-dimensional stacks of deconvolved epifluorescence microscopy, we found that inside living Bacillus subtilis cells, MreB forms filamentous structures of variable lengths, typically not longer than 1 μm. These filaments move along their orientation and mainly perpendicular to the long bacterial axis, revealing a maximal velocity at an intermediate length and a decreasing velocity with increasing filament length. Filaments move along straight trajectories but can reverse or alter their direction of propagation. Based on our measurements, we provide a mechanistic model that is consistent with all observations. In this model, MreB filaments mechanically couple several motors that putatively synthesize the cell wall, whereas the filaments’ traces mirror the trajectories of the motors. On the basis of our mechanistic model, we developed a mathematical model that can explain the nonlinear velocity length dependence. We deduce that the coupling of cell wall synthesis motors determines the MreB filament transport velocity, and the filament mechanically controls a concerted synthesis of parallel peptidoglycan strands to improve cell wall stability.