Search

Synthetic ion channels based on benzo(crown-ether) compounds have been previously reported to function as ion-selective channels in planar lipid bilayers, with hydrogen bonding networks implicated in the formation of self-aggregated complexes. Herein, we report the synthesis and characterization of two new families of benzo(crown-ether) compounds, termed monoacylated and monoalkylated benzo(crown-ethers) (MABCE), both of which lack hydrogen bond donors. Depending on the length of alkyl chain substituent and the size of macrocycle, MABCE compounds inhibit bacterial growth and transport ions across biological membranes. Single-channel recordings show that the activity is higher in the presence of K

The conserved fold of thioredoxin (Trx)-like thiol/disulfide oxidoreductases contains an invariant cis-proline residue (P76 in Escherichia coli Trx) that is essential for Trx function and that is responsible for the folding rate-limiting step. E. coli Trx contains four additional prolines, which are all in the trans conformation in the native state. Notably, a recent study revealed that replacement of all four trans prolines in Trx by alanines (Trx variant Trx1P) further slowed the rate-limiting step 25-fold, indicating that one or several of the four trans prolines accelerate the trans-to-cis transition of P76 in Trx wild-type (wt). Here, we characterized the folding kinetics of Trx variants containing cisP76 and one or several of the natural trans prolines of Trx wt with NMR spectroscopy. First, we demonstrate that the isomerization reaction in Trx1P is a pure two-state transition between two distinct tertiary structures, in which all observed NMR resonances changes follow the same first-order kinetics. Moreover, we show that trans-P68 is the critical residue responsible for the faster folding of wt Trx relative to the single-proline (P76) variant Trx1P, as the two-proline variant Trx2P(P76P68) already folds seven times faster than Trx1P. trans-P34 also accelerates trans-to-cis isomerization of P76, albeit to a smaller extent. Overall, the results demonstrate that trans prolines can significantly modulate the kinetics of rate-limiting trans-to-cis proline isomerization in protein folding. Finally, we discuss possible mechanisms of acceleration and the potential significance of a protein-internal folding acceleration mechanism for Trx in a living cell.

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

International audience; Evanescent-wave scattering is a topic in classical electrodynamics and in the study of colloidal particles near a boundary. However, how such near-surface scattering at subcellular refractive-index heterogeneities degrades the excitation confinement in biological total internal reflection fluorescence microscopy has not been well studied. An elegant theoretical work by Axelrod and Axelrod now addresses this very relevant question and reveals that-even when scattered-evanescent light preserves some of its surprising optical properties.

Upon contact with antigen presenting cells (APCs) cytotoxic T lymphocytes (T cells) establish a highly organized contact zone denoted as immunological synapse (IS). The formation of the IS implies relocation of the microtubule organizing center (MTOC) towards the contact zone, which necessitates a proper connection between MTOC and IS via dynamic microtubules (MTs). The efficiency of the MTs finding the IS within relevant time scale is, however, still illusive. We investigate how MTs search the three-dimensional constrained cellular volume for the IS and bind upon encounter to dynein anchored at the IS cortex. The search efficiency is estimated by calculating the time required for the MTs to reach the dynein-enriched region of the IS. In this study, we develop simple mathematical and numerical models incorporating relevant components of a cell and propose an optimal search strategy. Using the mathematical model, we have quantified the average search time for a wide range of model parameters and proposed an optimized set of values leading to the minimum capture time. Our results show that search times are minimal when the IS formed at the nearest or at the farthest sites on the cell surface, with respect to the perinuclear MTOC. The search time increases monotonically away from these two specific sites and are maximal at an intermediate position near the equator of the cell. We observed that search time strongly depends on the number of searching MTs and distance of the MTOC from the nuclear surface., 35 pages, 9 Figures, and 2 Tables

DNA strand displacement, where a single-stranded nucleic acid invades a DNA duplex, is pervasive in genomic processes and DNA engineering applications. The kinetics of strand displacement have been studied in bulk; however, the kinetics of the underlying strand exchange were obfuscated by a slow bimolecular association step. Here, we use a novel single-molecule Fluorescence Resonance Energy Transfer (smFRET) approach termed the "fission" assay to obtain the full distribution of first passage times of unimolecular strand displacement. At a frame time of 4.4 ms, the first passage time distribution for a 14-nt displacement domain exhibited a nearly monotonic decay with little delay. Among the eight different sequences we tested, the mean displacement time was on average 35 ms and varied by up to a factor of 13. The measured displacement kinetics also varied between complementary invaders and between RNA and DNA invaders of the same base sequence except for T$\rightarrow$U substitution. However, displacement times were largely insensitive to the monovalent salt concentration in the range of 0.25 M to 1 M. Using a one-dimensional random walk model, we infer that the single-step displacement time is in the range of $\sim 30 \mu s$ to $\sim 300 \mu s$ depending on the base identity. The framework presented here is broadly applicable to the kinetic analysis of multistep processes investigated at the single-molecule level., Comment: To be published in Biophysical Journal

Mechanical properties of the cell are important biomarkers for probing its architectural changes caused by cellular processes and/or pathologies. The development of microfluidic technologies has enabled measuring the cell’s mechanical properties at high throughput so that mechanical phenotyping can be applied to large samples in reasonable timescales. These studies typically measure the stiffness of the cell as the only mechanical biomarker and do not disentangle the rheological contributions of different structural components of the cell, including the cell cortex, the interior cytoplasm and its immersed cytoskeletal structures, and the nucleus. Recent advancements in high-speed fluorescent imaging have enabled probing the deformations of the cell cortex while also tracking different intracellular components in rates applicable to microfluidic platforms. We present a, to our knowledge, novel method to decouple the mechanics of the cell cortex and the cytoplasm by analyzing the correlation between the cortical deformations that are induced by external microfluidic flows and the nucleus displacements, induced by those cortical deformations, i.e., we use the nucleus as a high-throughput microrheological probe to study the rheology of the cytoplasm, independent of the cell cortex mechanics. To demonstrate the applicability of this method, we consider a proof-of-concept model consisting of a rigid spherical nucleus centered in a spherical cell. We obtain analytical expressions for the time-dependent nucleus velocity as a function of the cell deformations when the interior cytoplasm is modeled as a viscous, viscoelastic, porous, and poroelastic material and demonstrate how the nucleus velocity can be used to characterize the linear rheology of the cytoplasm over a wide range of forces and timescales/frequencies.

The COVID-19 pandemic has swept over the world in the past months, causing significant loss of life and consequences to human health. Although numerous drug and vaccine developments efforts are underway, many questions remain outstanding on the mechanism of SARS-CoV-2 viral association to angiotensin-converting enzyme 2 (ACE2), its main host receptor, and entry in the cell. Structural and biophysical studies indicate some degree of flexibility in the viral extracellular Spike glycoprotein and at the receptor binding domain-receptor interface, suggesting a role in infection. Here, we perform all-atom molecular dynamics simulations of the glycosylated, full-length membrane-bound ACE2 receptor, in both an apo and spike receptor binding domain (RBD) bound state, in order to probe the intrinsic dynamics of the ACE2 receptor in the context of the cell surface. A large degree of fluctuation in the full length structure is observed, indicating hinge bending motions at the linker region connecting the head to the transmembrane helix, while still not disrupting the ACE2 homodimer or ACE2-RBD interfaces. This flexibility translates into an ensemble of ACE2 homodimer conformations that could sterically accommodate binding of the spike trimer to more than one ACE2 homodimer, and suggests a mechanical contribution of the host receptor towards the large spike conformational changes required for cell fusion. This work presents further structural and functional insights into the role of ACE2 in viral infection that can be exploited for the rational design of effective SARS-CoV-2 therapeutics.Statement of SignificanceAs the host receptor of SARS-CoV-2, ACE2 has been the subject of extensive structural and antibody design efforts in aims to curtail COVID-19 spread. Here, we perform molecular dynamics simulations of the homodimer ACE2 full-length structure to study the dynamics of this protein in the context of the cellular membrane. The simulations evidence exceptional plasticity in the protein structure due to flexible hinge motions in the head-transmembrane domain linker region and helix mobility in the membrane, resulting in a varied ensemble of conformations distinct from the experimental structures. Our findings suggest a dynamical contribution of ACE2 to the spike glycoprotein shedding required for infection, and contribute to the question of stoichiometry of the Spike-ACE2 complex.

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.

International audience; Antimicrobial peptides are promising therapeutic agents to mitigate the global rise of antibiotic resistance. They generally act by perturbing the bacterial cell membrane and are thus less likely to induce resistance. Because they are membrane-active molecules, it is critical to verify and understand their potential action toward eukaryotic cells to help design effective and safe drugs. In this work, we studied the interaction of two antimicrobial peptides, aurein 1.2 and caerin 1.1, with red blood cell (RBC) membranes using in situ 31P and 2H solid-state NMR (SS-NMR). We established a protocol to integrate up to 25% of deuterated fatty acids in the membranes of ghosts, which are obtained when hemoglobin is removed from RBCs. Fatty acid incorporation and the integrity of the lipid bilayer were confirmed by SS-NMR and fluorescence confocal microscopy. Leakage assays were performed to assess the lytic power of the antimicrobial peptides. The in situ perturbation of the ghost membranes by aurein 1.2 and caerin 1.1 revealed by 31P and 2H SS-NMR is consistent with membrane perturbation through a carpet mechanism for aurein 1.2, whereas caerin 1.1 acts on RBCs via pore formation. These results are compatible with fluorescence microscopy images of the ghosts. The peptides interact with eukaryotic membranes following similar mechanisms that take place in bacteria, highlighting the importance of hydrophobicity when determining such interactions. Our work bridges model membranes and in vitro studies and provides an analytical toolbox to assess drug toxicity toward eukaryotic cells.

The sarcoplasmic reticulum Ca2+-ATPase (SERCA) transports two Ca2+ions from the cytoplasm to the reticulum lumen at the expense of ATP hydrolysis. In addition to transporting Ca2+, SERCA facilitates bidirectional proton transport across the sarcoplasmic reticulum to maintain the charge balance of the transport sites and to balance the charge deficit generated by the exchange of Ca2+. Previous studies have shown the existence of a transient water-filled pore in SERCA that connects the Ca2+-binding sites with the lumen, but the capacity of this pathway to sustain passive proton transport has remained unknown. In this study, we used the multiscale reactive molecular dynamics (MS-RMD) method and free energy sampling to quantify the free energy profile and timescale of the proton transport across this pathway while also explicitly accounting for the dynamically coupled hydration changes of the pore. We find that proton transport from the central binding site to the lumen has a microsecond timescale, revealing a novel passive cytoplasm-to-lumen proton flow beside the well-known inverse proton countertransport occurring in active Ca2+transport. We propose that this proton transport mechanism is operational and serves as a functional conduit for passive proton transport across the sarcoplasmic reticulum.SIGNIFICANCEMultiscale reactive molecular dynamics combined with free energy sampling was applied to study proton transport through a transient water pore connecting the Ca2+-binding site to the lumen in SERCA. This is the first computational study of this large biomolecular system that treats the hydrated excess proton and its transport through water structures and amino acids explicitly. When also correctly accounting for the hydration fluctuations of the pore, it is found that a transiently hydrated channel can transport protons on a microsecond timescale. These results quantitatively support the hypothesis of the proton intake into the sarcoplasm via SERCA, in addition to the well-known proton pumping by SERCA to the cytoplasm along with Ca2+transport.

Cell membranes mainly consist of lipid bilayers with an actively regulated composition. The underlying processes are still poorly understood, in particular, how the hundreds of components are controlled. Cholesterol has been found to correlate with phospholipid saturation for reasons that remain unclear. To better understand the link between cell membrane regulation and chemical composition, we establish a computational framework based on chemical reaction networks, resulting in multiple semigrand canonical ensembles. By running computer simulations, we show that regulating the chemical potential of lipid species is sufficient to reproduce the experimentally observed increase in acyl tail saturation with added cholesterol. Our model proposes a different picture of lipid regulation in which components can be regulated passively instead of actively. In this picture, phospholipid acyl tail composition naturally adapts to added molecules such as cholesterol or proteins. A comparison between regulated membranes with commonly studied ternary model membranes shows stark differences: for instance, correlation lengths and viscosities observed are independent of lipid chemical affinity.

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.

Voltage-dependent potassium (K(v)) channels play a fundamental role in neuronal and cardiac excitability and are potential therapeutic targets. They assemble as tetramers with a centrally located pore domain surrounded by a voltage-sensing domain (VSD), which is critical for sensing transmembrane potential and subsequent gating. Although the sensor is supposed to be in “Up” conformation in both n-octylglucoside (OG) micelles and phospholipid membranes in the absence of membrane potential, toxins that bind VSD and modulate the gating behavior of K(v) channels exhibit dramatic affinity differences in these membrane-mimetic systems. In this study, we have monitored the structural dynamics of the S3b-S4 loop of the paddle motif in activated conformation of KvAP-VSD by site-directed fluorescence approaches, using the environment-sensitive fluorescent probe 7-nitrobenz-2-oxa-1,3-diazol-4-yl-ethylenediamine (NBD). Emission maximum of NBD-labeled loop region of KvAP-VSD (residues 110–117) suggests a significant change in the polarity of local environment in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1′-rac-glycerol) membranes compared to OG micelles. This indicates that S3b-S4 loop residues might be partitioning to membrane interface, which is supported by an overall increased mean fluorescence lifetimes and significantly reduced water accessibility in membranes. Further, the magnitude of red edge excitation shift (REES) supports the presence of restricted/bound water molecules in the loop region of the VSD in micelles and membranes. Quantitative analysis of REES data using Gaussian probability distribution function clearly indicates that the sensor loop has fewer discrete equilibrium conformational states when reconstituted in membranes. Interestingly, this reduced molecular heterogeneity is consistent with the site-specific NBD polarization results, which suggest that the membrane environment offers a relaxed/dynamic organization for most of the S3b-S4 loop residues of the sensor. Overall, our results are relevant for understanding toxin-VSD interaction and gating mechanisms of K(v) channels in membranes.

Lipid bilayers can exhibit asymmetric states, in which the physical characteristics of one leaflet differ from those of the other. This most visibly manifests in a different lipid composition, but it can also involve opposing lateral stresses in each leaflet that combine to an overall vanishing membrane tension. Here, we use theoretical modeling and coarse-grained simulation to explore the interplay between a compositional asymmetry and a nonvanishing differential stress. Minimizing the total elastic energy leads to a preferred spontaneous curvature that balances torques due to both bending moments and differential stress, with sometimes unexpected consequences. For instance, asymmetric flat bilayers, whose specific areas in each leaflet are matched to those of corresponding tensionless symmetric flat membranes, still exhibit a residual differential stress because the conditions of vanishing area strain and vanishing bending moment differ. We also measure the curvature rigidity of asymmetric bilayers and find that a sufficiently strong differential stress, but not compositional asymmetry alone, can increase the bending modulus. The likely cause is a stiffening of the compressed leaflet, which appears to be related to its gel transition but not identical with it. We finally show that the impact of cholesterol on differential stress depends on the relative strength of elastic and thermodynamic driving forces: if cholesterol solvates equally well in both leaflets, it will redistribute to cancel both leaflet tensions almost completely, but if its partitioning free energy prefers one leaflet over the other, the resulting distribution bias may even create differential stress. Because cells keep most of their lipid bilayers in an asymmetric nonequilibrium steady state, our findings suggest that biomembranes are elastically more complex than previously thought: besides a spontaneous curvature, they might also exhibit significant differential stress, which could strongly affect their curvature energetics.

The nucleus is delineated by the nuclear envelope (NE), which is a double membrane barrier composed of the inner and outer nuclear membranes as well as a ∼40-nm wide lumen. In addition to its barrier function, the NE acts as a critical signaling node for a variety of cellular processes, which are mediated by protein complexes within this subcellular compartment. Although fluorescence fluctuation spectroscopy is a powerful tool for characterizing protein complexes in living cells, it was recently demonstrated that conventional fluorescence fluctuation spectroscopy methods are not suitable for applications in the NE because of the presence of slow nuclear membrane undulations. We previously addressed this challenge by developing time-shifted mean-segmented Q (tsMSQ) analysis and applied it to successfully characterize protein homo-oligomerization in the NE. However, many NE complexes, such as the linker of the nucleoskeleton and cytoskeleton complex, are formed by heterotypic interactions, which single-color tsMSQ is unable to characterize. Here, we describe the development of dual-color (DC) tsMSQ to analyze NE heteroprotein complexes built from proteins that carry two spectrally distinct fluorescent labels. Experiments performed on model systems demonstrate that DC tsMSQ properly identifies heteroprotein complexes and their stoichiometry in the NE by accounting for spectral cross talk and local volume fluctuations. Finally, we applied DC tsMSQ to study the assembly of the linker of the nucleoskeleton and cytoskeleton complex, a heteroprotein complex composed of Klarsicht/ANC-1/SYNE homology and Sad1/UNC-84 (SUN) proteins, in the NE of living cells. Using DC tsMSQ, we demonstrate the ability of the SUN protein SUN2 and the Klarsicht/ANC-1/SYNE homology protein nesprin-2 to form a heterocomplex in vivo. Our results are consistent with previously published in vitro studies and demonstrate the utility of the DC tsMSQ technique for characterizing NE heteroprotein complexes.

Neural precursor cells expressed developmentally downregulated protein 4-2 (Nedd4-2), a homologous to the E6-AP carboxyl terminus (HECT) ubiquitin ligase, triggers the endocytosis and degradation of its downstream target molecules by regulating signal transduction through interactions with other targets, including 14-3-3 proteins. In our previous study, we found that 14-3-3 binding induces a structural rearrangement of Nedd4-2 by inhibiting interactions between its structured domains. Here, we used time-resolved fluorescence intensity and anisotropy decay measurements, together with fluorescence quenching and mass spectrometry, to further characterize interactions between Nedd4-2 and 14-3-3 proteins. The results showed that 14-3-3 binding affects the emission properties of AEDANS-labeled WW3, WW4, and, to a lesser extent, WW2 domains, and reduces their mobility, but not those of the WW1 domain, which remains mobile. In contrast, 14-3-3 binding has the opposite effect on the active site of the HECT domain, which is more solvent exposed and mobile in the complexed form than in the apo form of Nedd4-2. Overall, our results suggest that steric hindrance of the WW3 and WW4 domains combined with conformational changes in the catalytic domain may account for the 14-3-3 binding-mediated regulation of Nedd4-2.

Tissue stiffness (Young's modulus) is a key control parameter in cell behavior and bioengineered gels where defined mechanical properties have become an essential part of the toolkit for interrogating mechanotransduction. Here, we show using a mechanical cell model that the effective substrate stiffness experienced by a cell depends, not just on the engineered mechanical properties of the substrate but critically also on the particular arrangement of adhesions between cell and substrate. In particular, we find that cells with different adhesion patterns can experience two different gel stiffnesses as equivalent and will generate the same mean cell deformations. In considering small patches of adhesion, which mimic focal adhesion complexes, we show how the experimentally observed focal adhesion growth and elongation on stiff substrates can be explained by energy considerations. Relatedly, energy arguments also provide a reason why nascent adhesions do not establish into focal adhesions on soft substrates, as has been commonly observed. Fewer and larger adhesions are predicted to be preferred over more and smaller, an effect enhanced by random spot placing with the simulations predicting qualitatively realistic cell shapes in this case.

Wounds can be produced when cells and tissues are subjected to excessive forces, for instance, under pathological conditions or nonphysiological loading. However, the cellular behaviors in the wound formation process are not clear. Here we tested the behaviors of wound formation in the epithelial layer with an in-suit uniaxial stretching device. We found that the wound often nucleates at the position where the cells are dividing. The polarization direction of cells near the wound is preferentially along the wound edge, whereas the cells far from the wound are preferentially perpendicular to the stretching direction. The larger the wound area is, the higher is the aspect ratio of the cells around the wound. Increasing the cell density will strengthen the cell layer. The higher the cell density is, the smaller is the area of the wounds, and the weaker is the effect of stretching on the polarization of the cells. Furthermore, we built a coarse-grained cell model that can explicitly consider the elasticity and viscoelasticity of cells, cell-cell interaction, and cell active stress, by which we simulated the wound formation process and quantitatively analyzed the force and stress fields in the cell layer, particularly around the wound. These analyses reveal the cellular mechanisms of wound formation behaviors in the cell layer under stretching and shed useful light on tissue engineering and regenerative medicine for biomedical applications.

Proteins often interconvert between different conformations in ways critical to their function. Although manipulating such equilibria for biophysical study is often challenging, the application of pressure is a potential route to achieve such control by favoring the population of lower volume states. Here, we use this feature to study the interconversion of ARNT PAS-B Y456T, which undergoes a dramatic +3 slip in the β-strand register as it switches between two stably folded conformations. Using high-pressure biomolecular NMR approaches, we obtained the first, to our knowledge, quantitative data testing two key hypotheses of this process: the slipped conformation is both smaller and less compressible than the wild-type equivalent, and the interconversion proceeds through a chiefly unfolded intermediate state. Data collected in steady-state pressure and time-resolved pressure-jump modes, including observed pressure-dependent changes in the populations of the two conformers and increased rate of interconversion between conformers, support both hypotheses. Our work exemplifies how these approaches, which can be generally applied to protein conformational switches, can provide unique information that is not easily accessible through other techniques.

Plant cell size and shape are tuned to their function and specified primarily by cellulose microfibril (CMF) patterning of the cell wall. Arabidopsis thaliana leaf trichomes are unicellular structures that act as a physical defense to deter insect feeding. This highly polarized cell type employs a strongly anisotropic cellulose wall to extend and taper, generating sharply pointed branches. During elongation, the mechanisms by which shifts in fiber orientation generate cells with predictable sizes and shapes are unknown. Specifically, the axisymmetric growth of trichome branches is often thought to result from axisymmetric CMF patterning. Here, we analyzed the direction and degree of twist of branches after desiccation to reveal the presence of an asymmetric cell wall organization with a left-hand bias. CMF organization, quantified using computational modeling, suggests a limited reorientation of microfibrils during growth and a maximum branch length limited by the wall axial stiffness. The model provides a mechanism for CMF asymmetry, which occurs after the branch bending stiffness becomes low enough that ambient bending affects the principal stresses. After this stage, the CMF synthesis results in a constant bending stiffness for longer branches. The bending vibration natural frequencies of branches with respect to their length are also discussed.

Recent evidence has demonstrated that, when cultured on micro-patterned surfaces, living cells can move in a coordinated manner and form distinct migration patterns, including flowing chain, suspended propagating bridge, rotating vortex, etc. However, the fundamental question of exactly how and why cells migrate in these fashions remains elusive. Here, we present a theoretical investigation to show that the tight interplay between internal cellular activities, such as chemo-mechanical feedbacks and polarization, and external geometrical constraints are behind these intriguing experimental observations. In particular, on narrow strip patterns, strongly force-dependent cellular contractility and intercellular adhesion were found to be critical for reinforcing the leading edge of the migrating cell monolayer and eventually result in the formation of suspended cell bridges flying over nonadhesive regions. On the other hand, a weak force-contractility feedback led to the movement of cells like a flowing chain along the adhesive strip. Finally, we also showed that the random polarity forces generated in migrating cells are responsible for driving them into rotating vortices on strips with width above a threshold value (~10 times the size of the cell).

The active isometric force-length relation (FLR) of striated muscle sarcomeres is central to understanding and modeling muscle function. The mechanistic basis of the descending arm of the FLR is well explained by the decreasing thin:thick filament overlap that occurs at long sarcomere lengths. The mechanistic basis of the ascending arm of the FLR (the decrease in force that occurs at short sarcomere lengths), alternatively, has never been well explained. Because muscle is a constant-volume system, interfilament lattice distances must increase as sarcomere length shortens. This increase would decrease thin and thick-filament electrostatic interactions independently of thin:thick filament overlap. To examine this effect, we present here a fundamental, physics-based model of the sarcomere that includes filament molecular properties, calcium binding, sarcomere geometry including both thin:thick filament overlap and interfilament radial distance, and electrostatics. The model gives extremely good fits to existing FLR data from a large number of different muscles across their entire range of measured activity levels, with the optimized parameter values in all cases lying within anatomically and physically reasonable ranges. A local first-order sensitivity analysis (varying individual parameters while holding the values of all others constant) shows that model output is most sensitive to a subset of model parameters, most of which are related to sarcomere geometry, with model output being most sensitive to interfilament radial distance. This conclusion is supported by re-running the fits with only this parameter subset being allowed to vary, which increases fit errors only moderately. These results show that the model well reproduces existing experimental data, and indicate that changes in interfilament spacing play as central a role as changes in filament overlap in determining the FLR, particularly on its ascending arm.

Fluorescence microscopy is an excellent tool to gain knowledge on cellular structures and biochemical processes. Stimulated emission depletion (STED) microscopy provides a resolution in the range of a few 10 nm at relatively fast data acquisition. As cellular structures can be oriented in any direction, it is of great benefit if the microscope exhibits an isotropic resolution. Here, we present an isoSTED microscope that utilizes water-immersion objective lenses and enables imaging of cellular structures with an isotropic resolution of better than 60 nm in living samples at room temperature and without CO(2) supply or another pH control. This corresponds to a reduction of the focal volume by far more than two orders of magnitude as compared to confocal microscopy. The imaging speed is in the range of 0.8 s/μm(3). Because fluorescence signal can only be detected from a diffraction-limited volume, a background signal is inevitably observed at resolutions well beyond the diffraction limit. Therefore, we additionally present a method that allows us to identify this unspecific background signal and to remove it from the image.

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 ability to accurately predict RNA hairpin structure and stability for different loop sequences and salt conditions is important for understanding, modeling, and designing larger RNA folds. However, traditional RNA secondary structure models cannot treat loop-sequence and ionic effects on RNA hairpin folding. Here, we describe a general, three-dimensional (3D) conformation-based computational method for modeling salt concentration-dependent conformational distributions and the detailed 3D structures for a set of three RNA hairpins that contain a variable, 15-nucleotide loop sequence. For a given RNA sequence, the new, to our knowledge, method integrates a Vfold2D two-dimensional structure folding model with IsRNA coarse-grained molecular dynamics 3D folding simulations and Monte Carlo tightly bound ion estimations of ion-mediated electrostatic interactions. The model predicts free-energy landscapes for the different RNA hairpin-forming sequences with variable salt conditions. The theoretically predicted results agree with the experimental fluorescence measurements, validating the strategy. Furthermore, the theoretical model goes beyond the experimental results by enabling in-depth 3D structural analysis, revealing energetic mechanisms for the sequence- and salt-dependent folding stability. Although the computational framework presented here is developed for RNA hairpin systems, the general method may be applied to investigate other RNA systems, such as multiway junctions or pseudoknots in mixed metal ion solutions.

Spatially discordant alternans (SDA) of action potential duration (APD) has been widely observed in cardiac tissue and is linked to cardiac arrhythmogenesis. Theoretical studies have shown that conduction velocity restitution (CVR) is required for the formation of SDA. However, this theory is not completely supported by experiments, indicating that other mechanisms may exist. In this study, we carried out computer simulations using mathematical models of action potentials to investigate the mechanisms of SDA in cardiac tissue. We show that when CVR is present and engaged, such as fast pacing from one side of the tissue, the spatial pattern of APD in the tissue undergoes either spatially concordant alternans or SDA, independent of initial conditions or tissue heterogeneities. When CVR is not engaged, such as simultaneous pacing of the whole tissue or under normal/slow heart rates, the spatial pattern of APD in the tissue can have multiple solutions, including spatially concordant alternans and different SDA patterns, depending on heterogeneous initial conditions or pre-existing repolarization heterogeneities. In homogeneous tissue, curved nodal lines are not stable, which either evolve into straight lines or disappear. However, in heterogeneous itssue, curved nodal lines can be stable, depending on their initial locations and shapes relative to the structure of the heterogeneity. Therefore, CVR-induced SDA and non-CVR-induced SDA exhibit different dynamical properties, which may be responsible for the different SDA properties observed in experimental studies and arrhythmogenesis in different clinical settings.

Fluorine incorporation is ideally suited to many NMR techniques, and incorporation of fluorine into proteins and fragment libraries for drug discovery has become increasingly common. Here, we use one-dimensional

In cell membranes, the functional constituents such as peptides, proteins, and polysaccharides diffuse in a sea of lipids as single molecules and molecular aggregates. Thus, the fluidity of the heterogeneous multicomponent membrane is important for understanding the roles of the membrane in cell functionality. Recently, Henle and Levine described the hydrodynamics of molecular diffusion in a spherical membrane. A tangential point force at the north pole induces a pair of vortices whose centers lie on a line perpendicular to the point force and are symmetrical with respect to the point force. The position of the vortex center depends on η(m)/Rη(w), where R is the radius of the spherical membrane, and η(m) and η(w) are the viscosities of the membrane and the surrounding medium, respectively. Based on this theoretical prediction, we applied a point force to a phase-separated spherical vesicle composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dioleoyl-sn-glycero-3-phosphocholine/cholesterol by means of a microinjection technique. The pathlines were visualized by trajectories of microdomains. We determined the position of the vortex center and estimated the membrane viscosity using the dependence of the position of the vortex center on η(m)/Rη(w). The obtained apparent membrane viscosities for various compositions are mapped on the phase diagram. The membrane viscosity is almost constant in the range of 0 < ϕ(Lo) ≤ 0.5 (ϕ(Lo): area fraction of the liquid ordered phase), whereas that in the range of 0.5 ≤ ϕ(Lo) < 1.0 exponentially increases with increase of ϕ(Lo). The obtained viscosity landscape provides a basic understanding of the fluidity of heterogeneous multicomponent membranes.

DNA under torsional strain undergoes a buckling transition that is the fundamental step in plectoneme nucleation and supercoil dynamics, which are critical for the processing of genomic information. Despite its importance, quantitative models of the buckling transition, in particular to also explain the surprising two-orders-of-magnitude difference between the buckling times for RNA and DNA revealed by single-molecule tweezers experiments, are currently lacking. Additionally, little is known about the configurations of the DNA during the buckling transition because they are not directly observable experimentally. Here, we use a discrete worm-like chain model and Brownian dynamics to simulate the DNA/RNA buckling transition. Our simulations are in good agreement with experimentally determined parameters of the buckling transition. The simulations show that the buckling time strongly and exponentially depends on the bending stiffness, which accounts for more than half the measured difference between DNA and RNA. Analyzing the microscopic conformations of the chain revealed by our simulations, we find clear evidence for a solenoid-shaped transition state and a curl intermediate. The curl intermediate features a single loop and becomes increasingly populated at low forces. Taken together, the simulations suggest that the worm-like chain model can account semiquantitatively for the buckling dynamics of both DNA and RNA.

Proteins interact with one another across a broad spectrum of affinities. Our understanding of the low end of this spectrum, as characterized by millimolar dissociation constants, relies on a handful of cases in which weak encounters have experimentally been identified. These weak interactions away from the specific target binding site can lead toward a higher-affinity complex. Recently, we detected weak encounters between two paralogous phosphotransferase pathways of Escherichia coli, which regulate various metabolic processes and stress responses. In addition to encounters that are known to occur between cognate proteins, i.e., those that can exchange phosphate groups with each other, surprisingly, encounters involving noncognates were also observed. It is not clear whether these “futile” encounters have a cooperative or competitive role. Using agent-based simulations, we find that the encounter complexes can be cooperative or competitive so as to increase or lower the effective binding affinity of the specific complex under different circumstances. This finding invites further questions into how organisms might exploit such low affinities to connect their signaling components.

Lateral segregation and the formation of lateral domains are well-known phenomena in ternary lipid bilayers composed of an unsaturated (low gel-to-liquid phase transition temperature (Tm)) phospholipid, a saturated (high-Tm) phospholipid, and cholesterol. The formation of lateral domains has been shown to be influenced by differences in phospholipid acyl chain unsaturation and length. Recently, we also showed that differential interactions of cholesterol with low- and high-Tm phospholipids in the bilayer can facilitate phospholipid segregation. Now, we have investigated phospholipid-cholesterol interactions and their role in lateral segregation in ternary bilayers composed of different unsaturated phosphatidylcholines (PCs) with varying acyl chain lengths, N-palmitoyl-D-erythro-sphingomyelin (PSM), and cholesterol. Using deuterium NMR spectroscopy, we determined how PSM was influenced by the acyl chain composition in surrounding PC environments and correlated this with the affinity of cholestatrienol (a fluorescent cholesterol analog) for PSM in the different PC environments. Results from a combination of time-resolved fluorescence measurements of trans-parinaric acid and Forster resonance energy transfer experiments showed that the relative affinity of cholesterol for phospholipids determined the degree to which the sterol promoted domain formation. From Forster resonance energy transfer, deuterium NMR, and differential scanning calorimetry results, it was clear that cholesterol also influenced both the thermostability of the domains and the degree of order in and outside the PSM-rich domains. The results of this study have shown that the affinity of cholesterol for both low-Tm and high-Tm phospholipids and the effects of low- and high-Tm phospholipids on each other influence both lateral structure and domain properties in complex bilayers. We envision that similar effects also contribute to lateral heterogeneity in even more complex biological membranes.

Cells contain numerous membraneless organelles that assemble by intracellular liquid-liquid phase separation. The viscous properties and associated biomolecular mobility within these condensed phase droplets, or condensates, are increasingly recognized as important for cellular function and also dysfunction, for example, in protein aggregation pathologies. Fluorescence recovery after photobleaching (FRAP) is widely used to assess condensate fluidity and to estimate protein diffusion coefficients. However, the models and assumptions utilized in FRAP analysis of protein condensates are often not carefully considered. Here, we combine FRAP experiments on both in vitro reconstituted droplets and intracellular condensates with systematic examination of different models that can be used to fit the data and evaluate the impact of model choice on measured values. A key finding is that model boundary conditions can give rise to widely divergent measured values. This has important implications, for example, in experiments that bleach subregions versus the entire condensate, two commonly employed experimental approaches. We suggest guidelines for determining the appropriate modeling framework and highlight emerging questions about the molecular dynamics at the droplet interface. The ability to accurately determine biomolecular mobility both in the condensate interior and at the interface is important for obtaining quantitative insights into condensate function, a key area for future research.

We consider the plasma membrane that contains a cholesterol molar fraction of 0.4 and ask how that cholesterol is distributed between the two leaves. Because of the rapid flip-flop of cholesterol between leaves, we assume that its distribution is determined by the equality of its chemical potentials in the two leaves. When we consider only the contributions of entropy and interactions to the cholesterol chemical potential in our model system, we find, not surprisingly, that the cholesterol is mostly in the outer leaf because of the strong attraction between cholesterol and sphingomyelin (SM), which is predominantly in that leaf. We find 72% there. We then include the contribution from the bending energy in each leaf that must be overcome to join the leaves in a flat bilayer. The product of bending modulus and spontaneous curvature is obtained from simulation. We find that the addition of cholesterol to the outer leaf reduces the spontaneous curvature, which is initially positive, until it passes through zero when the molar fraction of cholesterol in the outer leaf is 0.28. Additional cholesterol is driven toward the inner leaf by the sphingomyelin phosphatidylcholine mixture. This is resisted by the bending energy contribution to the inner leaf. We find, again by simulation, that the addition of cholesterol monotonically increases the magnitude of the spontaneous curvature of the inner leaf, which is negative. This increases its bending energy. We conclude that, as a result of these competing effects, the percentage of cholesterol in the outer leaf is reduced to ∼63 ± 6%.

Proteins are dynamic molecules that undergo conformational changes to a broad spectrum of different excited states. Unfortunately, the small populations of these states make it difficult to determine their structures or functional implications. Computer simulations are an increasingly powerful means to identify and characterize functionally relevant excited states. However, this advance has uncovered a further challenge: it can be extremely difficult to identify the most salient features of large simulation data sets. We reasoned that many functionally relevant conformational changes are likely to involve large, cooperative changes to the surfaces that are available to interact with potential binding partners. To examine this hypothesis, we introduce a method that returns a prioritized list of potentially functional conformational changes by segmenting protein structures into clusters of residues that undergo cooperative changes in their solvent exposure, along with the hierarchy of interactions between these groups. We term these groups exposons to distinguish them from other types of clusters that arise in this analysis and others. We demonstrate, using three different model systems, that this method identifies experimentally validated and functionally relevant conformational changes, including conformational switches, allosteric coupling, and cryptic pockets. Our results suggest that key functional sites are hubs in the network of exposons. As a further test of the predictive power of this approach, we apply it to discover cryptic allosteric sites in two different β-lactamase enzymes that are widespread sources of antibiotic resistance. Experimental tests confirm our predictions for both systems. Importantly, we provide the first evidence, to our knowledge, for a cryptic allosteric site in CTX-M-9 β-lactamase. Experimentally testing this prediction did not require any mutations and revealed that this site exerts the most potent allosteric control over activity of any pockets found in β-lactamases to date. Discovery of a similar pocket that was previously overlooked in the well-studied TEM-1 β-lactamase demonstrates the utility of exposons.

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.