TolC is the trimeric outer membrane component of the efflux pump system in Escherichia coli that is responsible for antibiotic efflux from bacterial cells. Overexpression of efflux pumps has been reported to decrease susceptibility to antibiotics in a variety of bacterial pathogens. Reliable production of membrane proteins allows for the biophysical and structural characterization needed to better understand efflux and for the development of therapeutics. Preparation of recombinant protein for biochemical/structural studies often involves the production of proteins as inclusion body aggregates from which active proteins are recovered. Here, we find that the in vitro folding of TolC into its functional trimeric state from inclusion bodies is dependent on the headgroup composition of detergent micelles used. Nonionic detergent favors the formation of functional trimeric TolC, whereas zwitterionic detergents induce the formation of a non-native, oligomeric TolC fold. We also find that nonionic detergents with shorter alkyl lengths facilitate TolC folding. It remains to be seen whether the charges in lipid headgroups have similar effects on membrane insertion and folding in biological systems.

The polymerase chain reaction (PCR) is a central technique in biotechnology. Its ability to amplify a specific target region of a DNA sequence has led to prominent applications, including virus tests, DNA sequencing, genotyping, and genome cloning. These applications rely on the specificity of the primer hybridization and therefore require effective suppression of hybridization errors. A simple and effective method to achieve that is to add blocker strands, also called clamps, to the PCR mixture. These strands bind to the unwanted target sequence, thereby blocking the primer mishybridization. Because of its simplicity, this method is applicable to a broad nucleic-acid-based biotechnology. However, the precise mechanism by which blocker strands suppress PCR errors remains to be understood, limiting the applicability of this technique. Here, we combine experiments and theoretical modeling to reveal this mechanism. We find that the blocker strands both energetically destabilize the mishybridized complex and sculpt a kinetic barrier to suppress mishybridization. This combination of energetic and kinetic biasing extends the viable range of annealing temperatures, which reduces design constraint of the primer sequence and extends the applicability of PCR.

When multivesicular endosomes (MVEs) fuse with the plasma membrane, exosomes are released into the extracellular space where they can affect other cells. The ability of exosomes to regulate cells nearby or further away depends on whether they remain attached to the secreting cell membrane. The regulation and kinetics of exosome secretion are not well characterized, but probes for directly imaging single MVE fusion events have allowed for visualization of the fusion and release process. In particular, the design of an exosome marker with a pH-sensitive dye in the middle of the tetraspanin protein CD63 has facilitated studies of individual MVE fusion events. Using TIRF microscopy, single fusion events were measured in A549 cells held at 23-37°C and events were identified using an automated detection algorithm. Stable docking precedes fusion almost always and a decrease in temperature was accompanied by decrease in the rate of content loss and in the frequency of fusion events. The loss of CD63-pHluorin fluorescence was measured at fusion sites and fit with a single or double exponential decay, with most events requiring two components and a plateau because the loss of fluorescence was typically incomplete. To interpret the kinetics, fusion events were simulated as a localized release of tethered/untethered exosomes coupled with the membrane diffusion of CD63. The experimentally observed decay required three components in the simulation: 1) free exosomes, 2) CD63 membrane diffusion from the endosomal membrane into the plasma membrane, and 3) tethered exosomes. Modeling with slow diffusion of the tethered exosomes (0.0015-0.004 μm2/s) accurately fits the experimental data for all temperatures. However, simulating with immobile tethers or the absence of tethers fails to replicate the data. Our model suggests that exosome release from the fusion site is incomplete due to postfusion, membrane attachment.

Osmolytes are ubiquitous in the cell and play an important role in controlling protein stability under stress. The natural osmolyte trimethylamine N-oxide (TMAO) is used by marine animals to counteract the effect of pressure denaturation at large depths. The molecular mechanism of TMAO stabilization against pressure and urea denaturation has been extensively studied, but unlike the case of other osmolytes, the ability of TMAO to protect proteins from high temperature has not been quantified. To reveal the effect of TMAO on folded and unfolded protein ensembles and the hydration shell at different temperatures, we study a mutant of the well-characterized, fast-folding model protein B (PRB). We carried out, in total, >190 μs all-atom simulations of thermal folding/unfolding of PRB at multiple temperatures and concentrations of TMAO. The simulations show increased thermal stability of PRB in the presence of TMAO. Partly structured, compact ensembles are favored over the unfolded state. TMAO forms two shells near the protein: an outer shell away from the protein surface has altered H-bond lifetimes of water molecules and increases hydration of the protein to help stabilize it; a less-populated inner shell with an opposite TMAO orientation closer to the protein surface binds exclusively to basic side chains. The cooperative cosolute effect of the inner and outer shell TMAO has a small number of TMAO molecules "herding" water molecules into two hydration shells at or near the protein surface. The stabilizing effect of TMAO on our protein saturates at 1 M despite higher TMAO solubility, so there may be little evolutionary pressure for extremophiles to produce higher intracellular TMAO concentrations, if true in general.

Smoothened (SMO) is a membrane protein of the class F subfamily of G protein-coupled receptors (GPCRs) and maintains homeostasis of cellular differentiation. SMO undergoes conformational change during activation, transmitting the signal across the membrane, making it amenable to bind to its intracellular signaling partner. Receptor activation has been studied at length for class A receptors, but the mechanism of class F receptor activation remains unknown. Agonists and antagonists bound to SMO at sites in the transmembrane domain (TMD) and the cysteine-rich domain have been characterized, giving a static view of the various conformations SMO adopts. Although the structures of the inactive and active SMO outline the residue-level transitions, a kinetic view of the overall activation process remains unexplored for class F receptors. We describe SMO's activation process in atomistic detail by performing 300 μs of molecular dynamics simulations and combining it with Markov state model theory. A molecular switch, conserved across class F and analogous to the activation-mediating D-R-Y motif in class A receptors, is observed to break during activation. We also show that this transition occurs in a stage-wise movement of the transmembrane helices: TM6 first, followed by TM5. To see how modulators affect SMO activity, we simulated agonist and antagonist-bound SMO. We observed that agonist-bound SMO has an expanded hydrophobic tunnel in SMO's core TMD, whereas antagonist-bound SMO shrinks this tunnel, further supporting the hypothesis that cholesterol travels through a tunnel inside Smoothened to activate it. In summary, this study elucidates the distinct activation mechanism of class F GPCRs and shows that SMO's activation process rearranges the core TMD to open a hydrophobic conduit for cholesterol transport.

Fluorescently labeled oligonucleotides are powerful tools for characterizing DNA processes; however, their use is limited by the cost and sequence requirements of current labeling technologies. Here, we develop an easy, inexpensive, and sequence-independent method for site-specifically labeling DNA oligonucleotides. We utilize commercially synthesized oligonucleotides containing phosphorothioate diester(s) in which a nonbridging oxygen is replaced with a sulfur (PS-DNA). The increased nucleophilicity of the thiophosphoryl sulfur relative to the phosphoryl oxygen permits selective reactivity with iodoacetamide compounds. As such, we leverage a long-existing bifunctional linker, N,N'-bis(α-iodoacetyl)-2-2'-dithiobis(ethylamine) (BIDBE), that reacts with PS-DNAs to leave a free thiol, allowing conjugation of the wide variety of commercial maleimide-functionalized compounds. We optimized BIDBE synthesis and its attachment to PS-DNA and then fluorescently labeled the BIDBE-PS-DNA using standard protocols for labeling cysteines. We purified the individual epimers, and using single-molecule Förster resonance energy transfer (FRET), we show that the FRET efficiency is independent of the epimeric attachment. Subsequently, we demonstrate that an epimeric mixture of double-labeled Holliday junctions (HJs) can be used to characterize their conformational properties in the absence and presence of the structure-specific endonuclease Drosophila melanogaster Gen. Finally, we use a biochemical activity assay to show that this double-labeled HJ is functional for cleavage by Gen and that the double-labeled HJ allows multiple DNA species to be identified in a single experiment. In conclusion, our results indicate that dye-labeled BIDBE-PS-DNAs are comparable to commercially labeled DNAs at a significantly reduced cost. Notably, this technology could be applied to other maleimide-functionalized compounds, such as spin labels, biotin, and proteins. The sequence independence of labeling, coupled with its ease and low cost, enables unrestricted exploration of dye placement and choice, providing the potential for creation of differentially labeled DNA libraries and opening previously inaccessible experimental avenues.

Essential cellular processes such as metabolism, protein synthesis, and autophagy require the intracellular transport of membrane-bound vesicles. The importance of the cytoskeleton and associated molecular motors for transport is well documented. Recent research has suggested that the endoplasmic reticulum (ER) may also play a role in vesicle transport through a tethering of vesicles to the ER. We use single-particle tracking fluorescence microscopy and a Bayesian change-point algorithm to characterize vesicle motility in response to the disruption of the ER, actin, and microtubules. This high-throughput change-point algorithm allows us to efficiently analyze thousands of trajectory segments. We find that palmitate-mediated disruption of the ER leads to a significant decrease in vesicle motility. A comparison with the disruption of actin and microtubules shows that disruption of the ER has a significant impact on vesicle motility, greater than the disruption of actin. Vesicle motility was dependent on cellular region, with greater motility in the cell periphery than the perinuclear region, possibly due to regional differences in actin and the ER. Overall, these results suggest that the ER is an important factor in vesicle transport.

NADH and NADPH play key roles in the regulation of metabolism. Their endogenous fluorescence is sensitive to enzyme binding, allowing changes in cellular metabolic state to be determined using fluorescence lifetime imaging microscopy (FLIM). However, to fully uncover the underlying biochemistry, the relationships between their fluorescence and binding dynamics require greater understanding. Here we accomplish this through time- and polarization-resolved fluorescence and polarized two-photon absorption measurements. Two lifetimes result from binding of both NADH to lactate dehydrogenase and NADPH to isocitrate dehydrogenase. The composite fluorescence anisotropy indicates the shorter (1.3-1.6 ns) decay component to be accompanied by local motion of the nicotinamide ring, pointing to attachment solely via the adenine moiety. For the longer lifetime (3.2-4.4 ns), the nicotinamide conformational freedom is found to be fully restricted. As full and partial nicotinamide binding are recognized steps in dehydrogenase catalysis, our results unify photophysical, structural, and functional aspects of NADH and NADPH binding and clarify the biochemical processes that underlie their contrasting intracellular lifetimes.

D-Dopachrome tautomerase (D-DT; or MIF-2) is a multifunctional protein with immunomodulatory properties and a documented pathogenic role in inflammation and cancer that is associated with activation of the cell surface receptor CD74. Alongside D-DT, macrophage migration inhibitory factor (MIF) is also known to activate CD74, promoting pathogenesis. While the role of the MIF/CD74 axis has been extensively studied in various disease models, the late discovery of the D-DT/CD74 axis has led to a poor investigation into the D-DT-induced activation mechanism of CD74. A previous study has identified 4-(3-carboxyphenyl)-2,5-pyridinedicarboxylic acid (4-CPPC) as the first selective and reversible inhibitor of D-DT and reported its potency to block the D-DT-induced activation of CD74 in a cell-based model. In this study, we employ molecular dynamics simulations and nuclear magnetic resonance experiments to study 4-CPPC-induced changes to the dynamic profile of D-DT. We found that binding of the inhibitor remarkably promotes the conformational flexibility of C-terminal without impacting the structural stability of the biological assembly. Consequently, long-range intrasubunit (>11 Å) and intersubunit (>30 Å) communications are enabled between distal regions. Communication across the three subunits is accomplished via 4-CPPC, which serves as a communication bridge after Val113 is displaced from its hydrophobic pocket. This previously unrecognized structural property of D-DT is not shared with its human homolog, MIF, which exhibits an impressive C-terminal rigidity even in the presence of an inhibitor. Considering the previously reported role of MIF's C-terminal in the activation of CD74, our results break new ground for understanding the functionality of D-DT in health and disease.

Proliferating bacterial cells exhibit stochastic growth and size dynamics, but the regulation of noise in bacterial growth and morphogenesis remains poorly understood. A quantitative understanding of morphogenetic noise control, and how it changes under different growth conditions, would provide better insights into cell-to-cell variability and intergenerational fluctuations in cell physiology. Using multigenerational growth and width data of single Escherichia coli and Caulobacter crescentus cells, we deduce the equations governing growth and size dynamics of rod-like bacterial cells. Interestingly, we find that both E. coli and C. crescentus cells deviate from exponential growth within the cell cycle. In particular, the exponential growth rate increases during the cell cycle irrespective of nutrient or temperature conditions. We propose a mechanistic model that explains the emergence of super-exponential growth from autocatalytic production of ribosomes coupled to the rate of cell elongation and surface area synthesis. Using this new model and statistical inference on large datasets, we construct the Langevin equations governing cell growth and size dynamics of E. coli cells in different nutrient conditions. The single-cell level model predicts how noise in intragenerational and intergenerational processes regulate variability in cell morphology and generation times, revealing quantitative strategies for cellular resource allocation and morphogenetic noise control in different growth conditions.

We present a method for extracting temperature-dependent thermodynamic and photophysical properties of SYTO-13 dye bound to DNA from fluorescence measurements. Together, mathematical modeling, control experiments, and numerical optimization enable dye binding strength, dye brightness, and experimental noise (or error) to be discriminated from one another. By focusing on the low-dye-coverage regime, the model avoids bias and can simplify quantification. Utilizing the temperature-cycling capabilities and multi-reaction chambers of a real-time PCR machine increases throughput. Significant well-to-well and plate-to-plate variation is quantified by using total least squares to account for error in both fluorescence and nominal dye concentration. Properties computed independently for single-stranded DNA and double-stranded DNA by numerical optimization are consistent with intuition and explain the advantageous performance of SYTO-13 in high-resolution melting and real-time PCR assays. Distinguishing between binding, brightness, and noise also clarifies the mechanism for the increased fluorescence of dye in a solution of double-stranded DNA compared to single-stranded DNA; in fact, the explanation changes with temperature.

Although cells with distinct apical areas have been widely observed in epithelial tissues, how the size of cells affects their behavior during tissue deformation and morphogenesis as well as key physical factors modulating such influence remains elusive. Here, we showed that the elongation of cells within the monolayer under anisotropic biaxial stretching increases with their size because the strain released by local cell rearrangement (i.e., T1 transition) is more significant for small cells that possess higher contractility. On the other hand, by incorporating the nucleation, peeling, merging, and breakage dynamics of subcellular stress fibers into classical vertex formulation, we found that stress fibers with orientations predominantly aligned with the main stretching direction will be formed at tricellular junctions, in good agreement with recent experiments. The contractile forces generated by stress fibers help cells to resist imposed stretching, reduce the occurrence of T1 transitions, and, consequently, modulate their size-dependent elongation. Our findings demonstrate that epithelial cells could utilize their size and internal structure to regulate their physical and related biological behaviors. The theoretical framework proposed here can also be extended to investigate the roles of cell geometry and intracellular contraction in processes such as collective cell migration and embryo development.

Single-channel patch-clamp recordings allow observing the action of a single protein complex in real time and hence the deduction of the underlying conformational changes in the ion-channel protein. Commonly, recordings are modeled using hidden Markov chains, connecting open and closed states in the experimental data with protein conformations. The rates between states denote transition probabilities that could be modified by membrane voltage or ligand binding. Modeling algorithms have to deal with limited recording bandwidth and a very noisy background. It was previously shown that the fit of two-dimensional (2D)-dwell-time histograms with simulations is very robust in that regard. Errors introduced by the low-pass filter or noise cancel out to a certain degree when comparing experimental and simulated data. In addition, the topology of models (that is, the chain of open and closed states) could be inferred from 2D-histograms. However, the 2D-fit was never applied to its full potential. A major reason may be the extremely time-consuming and often unreliable fitting process, due to the stochastic variability in the simulations. We have now solved these issues by introducing a message-passing interface (MPI) allowing massive parallel computing on a high-performance computing (HPC) cluster and obtaining ensemble solutions. With ensembles, we have demonstrated how important ranked solutions are for difficult tasks related to a noisy background, fast gating events beyond the corner frequency of the low-pass filter, and topology estimation of the underlying Markov model. Finally, we have shown that, by combining the objective function of the 2D-fit with the deviation of the current amplitude distributions, automatic determination of the current level of the conducting state is possible, even with an apparent current reduction due to low-pass filtering. Making use of an HPC cluster, the power of 2D-dwell-time analysis can be used to its fullest with minor input of the experimenter.

The four-point-one ezrin-radixin-moesin homology (FERM) protein domain is a multifunctional protein-lipid binding site, constituting an integral part of numerous membrane-associated proteins. Its interaction with the lipid phosphatidylinositol-4,5-bisphosphate (PIP2), located at the inner leaflet of eukaryotic plasma membranes, is important for localization, anchorage, and activation of FERM-containing proteins. FERM-PIP2 complexes structurally determined so far exclusively feature a 1:1 binding stoichiometry of protein and lipid, with a few basic FERM residues neutralizing the -4 charge of the bound PIP2. Whether this picture from static crystal structures also applies to the dynamic interaction of FERM domains on PIP2 membranes is unknown. We here quantified the stoichiometry of FERM-PIP2 binding in a lipid bilayer using atomistic molecular dynamics simulations and experiments on solid supported membranes for the FERM domains of focal adhesion kinase and ezrin. In contrast to the structural data, we find much higher average stoichiometries of FERM-PIP2 binding, amounting to 1:3 or 1:4 ratios, respectively. In simulations, the full set of basic residues at the membrane interface, 7 and 15 residues for focal adhesion kinase and ezrin, respectively, engages in PIP2 interactions. In addition, Na ions enter the FERM-membrane binding interface, compensating negative PIP2 charges in case of high charge surpluses from bound PIP2. We propose the multivalent binding of FERM domains to PIP2 in lipid bilayers to significantly enhance the stability of FERM-membrane binding and to render the FERM-membrane linkage highly adjustable.

The 3D organization of eukaryotic genomes plays an important role in genome function. While significant progress has been made in deciphering the folding mechanisms of individual chromosomes, the principles of the dynamic large-scale spatial arrangement of all chromosomes inside the nucleus are poorly understood. We use polymer simulations to model the diploid human genome compartmentalization relative to nuclear bodies such as nuclear lamina, nucleoli, and speckles. We show that a self-organization process based on a cophase separation between chromosomes and nuclear bodies can capture various features of genome organization, including the formation of chromosome territories, phase separation of A/B compartments, and the liquid property of nuclear bodies. The simulated 3D structures quantitatively reproduce both sequencing-based genomic mapping and imaging assays that probe chromatin interaction with nuclear bodies. Importantly, our model captures the heterogeneous distribution of chromosome positioning across cells while simultaneously producing well-defined distances between active chromatin and nuclear speckles. Such heterogeneity and preciseness of genome organization can coexist due to the nonspecificity of phase separation and the slow chromosome dynamics. Together, our work reveals that the cophase separation provides a robust mechanism for us to produce functionally important 3D contacts without requiring thermodynamic equilibration that can be difficult to achieve.

Small-conductance (SK) calcium-activated potassium channels are a promising treatment target in atrial fibrillation. However, the functional properties that differentiate SK inhibitors remain poorly understood. The objective of this study was to determine how two unrelated SK channel inhibitors, apamin and AP14145, impact SK channel function in excised inside-out single-channel recordings. Surprisingly, both apamin and AP14145 exert much of their inhibition by inducing a class of very-long-lived channel closures (apamin: τc,vl = 11.8 ± 7.1 s, and AP14145: τc,vl = 10.3 ± 7.2 s), which were never observed under control conditions. Both inhibitors also induced changes to the three closed and two open durations typical of normal SK channel gating. AP14145 shifted the open duration distribution to favor longer open durations, whereas apamin did not alter open-state kinetics. AP14145 also prolonged the two shortest channel closed durations (AP14145: τc,s = 3.50 ± 0.81 ms, and τc,i = 32.0 ± 6.76 ms versus control: τc,s = 1.59 ± 0.19 ms, and τc,i = 13.5 ± 1.17 ms), thus slowing overall gating kinetics within bursts of channel activity. In contrast, apamin accelerated intraburst gating kinetics by shortening the two shortest closed durations (τc,s = 0.75 ± 0.10 ms and τc,i = 5.08 ± 0.49 ms) and inducing periods of flickery activity. Finally, AP14145 introduced a unique form of inhibition by decreasing unitary current amplitude. SK channels exhibited two clearly distinguishable amplitudes (control: Ahigh = 0.76 ± 0.03 pA, and Alow = 0.54 ± 0.03 pA). AP14145 both reduced the fraction of patches exhibiting the higher amplitude (AP14145: 4/9 patches versus control: 16/16 patches) and reduced the mean low amplitude (0.38 ± 0.03 pA). Here, we have demonstrated that both inhibitors introduce very long channel closures but that each also exhibits unique effects on other components of SK gating kinetics and unitary current. The combination of these effects is likely to be critical for understanding the functional differences of each inhibitor in the context of cyclical Ca2+-dependent channel activation in vivo.

Huntingtin (HTT) is a scaffolding protein that recruits motor proteins to vesicular cargoes, enabling it to regulate kinesin-1, dynein, and myosin-VI-dependent transport. To maintain the native stoichiometry of HTT with its interacting partners, we used CRISPR/Cas9 to induce a phosphomimetic mutation of the endogenous HTT at S421 (HTT-S421D). Using single-particle tracking, optical tweezers, and immunofluorescence, we examined the effects of this mutation on the motility of early endosomes and lysosomes. In HTT-S421D cells, lysosomes exhibit longer displacements and higher processive fractions compared with wild-type (HTT-WT) cells. Kinesins and dyneins exert greater forces on early endosomes and lysosomes in cells expressing HTT-S421D. In addition, endosomes bind to microtubules faster and are more resistant to detachment under load. The recruitment of kinesins and dyneins to microtubules is enhanced in HTT-S421D cells. In contrast, overexpression of HTT had variable effects on the processivity, displacement, and directional bias of both early endosomes and lysosomes. These data indicate that phosphorylation of the endogenous HTT causes early endosomes and lysosomes to move longer distances and more processively by recruiting and activating both kinesin-1 and dynein.

Noncoding RNAs (ncRNAs) are an emerging epigenetic factor and have been recognized as playing a key role in many gene expression pathways. Structurally, binding of ncRNAs to isolated DNA is strongly dependent on sequence complementary and results in the formation of an RNA.DNA-DNA (RDD) triple helix. However, in vivo DNA is not isolated but is rather packed in chromatin fibers, the fundamental unit of which is the nucleosome. Biochemical experiments have shown that ncRNA binding to nucleosomal DNA is elevated at DNA entry and exit sites and is dependent on the presence of the H3 N-terminal tails. However, the structural and dynamical bases for these mechanisms remain unknown. Here, we have examined the mechanisms and effects of RDD formation in the context of the nucleosome using a series of all-atom molecular dynamics simulations. Results highlight the importance of DNA sequence on complex stability, elucidate the effects of the H3 tails on RDD structures, show how RDD formation impacts the structure and dynamics of the H3 tails, and show how RNA alters the local and global DNA double-helical structure. Together, our results suggest ncRNAs can modify nucleosome, and potentially higher-order chromatin, structures and dynamics as a means of exerting epigenetic control.

The parasitic weed Striga has led to billions of dollars' worth of agricultural productivity loss worldwide. Striga detects host plants using compounds of the strigolactone class of phytohormones. Early steps in the strigolactone signaling pathway involve substrate binding and hydrolysis followed by a conformational change to an "active" or "closed" state, after which it associates with a MAX2-family downstream signaling partner. The structures of the inactive and active states of strigolactone receptors are known through X-ray crystallography, and the transition pathway from the inactive to active state in apo receptors has previously been characterized using molecular dynamics simulations. However, it also has been suggested that a covalent butenolide modification of the receptor on the catalytic histidine through substrate hydrolysis promotes formation of the active state. Using molecular dynamics simulations, we show that the presence of the covalent butenolide enhances activation in both AtD14, a receptor found in Arabidopsis, and ShHTL7, a receptor found in Striga, but the enhancement is ∼50 times greater in ShHTL7. We also show that several conserved interactions with the covalent butenolide modification promote transition to the active state in both AtD14 (non-parasite) and ShHTL7 (parasite). Finally, we demonstrate that the enhanced activation of ShHTL7 likely results from disruption of ShHTL7-specific histidine interactions that inhibited activation in the apo case. These results provide a possible explanation for difference in strigolactone sensitivity seen between different strigolactone-sensitive proteins and can be used to aid the design of selective modulators to control Striga parasites.

Focal adhesions (FAs) mediate the interaction of the cytoskeleton with the extracellular matrix in a highly dynamic fashion. Talin is a central regulator, adaptor protein, and mechano-sensor of FA complexes. For recruitment and firm attachment at FAs, talin's N-terminal FERM domain binds to phosphatidylinositol 4,5-bisphosphate (PIP2)-enriched membranes. A newly published autoinhibitory structure of talin-1, where the known PIP2 interaction sites are covered up, lead us to hypothesize that a hitherto less examined loop insertion of the FERM domain acts as an additional and initial site of contact. We evaluated direct interactions of talin-1 with a PIP2 membrane by means of atomistic molecular dynamics simulations. We show that this unstructured, 33-residue-long loop strongly interacts with PIP2 and can facilitate further membrane contacts, including the canonical PIP2 interactions, by serving as a flexible membrane anchor. Under force as present at FAs, the extensible FERM loop ensures talin maintains membrane contacts when pulled away from the membrane by up to 7 nm. We identify key basic residues of the anchor mediating the highly dynamic talin-membrane interaction. Our results put forward an intrinsically disordered loop as a key and highly adaptable PIP2 recognition site of talin and potentially other PIP2-binding mechano-proteins.

Hydrophobic gating is an emerging mechanism in regulation of protein ion channels where the pore remains physically open but becomes dewetted to block ion permeation. Atomistic molecular dynamics simulations have played a crucial role in understanding hydrophobic gating by providing the molecular details to complement mutagenesis and structural studies. However, existing studies rely on direct simulations and do not quantitatively describe how the sequence and structural changes may control the delicate liquid-vapor equilibrium of confined water in the pore of the channel protein. To address this limitation, we explore two enhanced sampling methods, namely metadynamics and umbrella sampling, to derive free-energy profiles of pore hydration in both the closed and open states of big potassium (BK) channels, which are important in cardiovascular and neural systems. It was found that metadynamics required substantially longer sampling times and struggled to generate stably converged free-energy profiles due to the slow dynamics of cooperative pore water diffusion even in the barrierless limit. Using umbrella sampling, well-converged free-energy profiles can be readily generated for the wild-type BK channels as well as three mutants with pore-lining mutations experimentally known to dramatically perturb the channel gating voltage. The results show that the free energy of pore hydration faithfully reports the gating voltage of the channel, providing further support for hydrophobic gating in BK channels. Free-energy analysis of pore hydration should provide a powerful approach for quantitative studies of how protein sequence, structure, solution conditions, and/or drug binding may modulate hydrophobic gating in ion channels.

Optical methods for examining cellular structure based on endogenous contrast rely on analysis of refractive index changes to discriminate cell phenotype. These changes can be visualized using techniques such as phase contrast microscopy, detected by light scattering, or analyzed numerically using quantitative phase imaging. The statistical variations of refractive index at the nanoscale can be quantified using disorder strength, a metric seen to increase with neoplastic change. In contrast, the spatial organization of these variations is typically characterized using a fractal dimension, which is also seen to increase with cancer progression. Here, we seek to link these two measurements using multiscale measurements of optical phase to calculate disorder strength and in turn to determine the fractal dimension of the structures. First, quantitative phase images are analyzed to show that the disorder strength metric changes with resolution. The trend of disorder strength with length scales is analyzed to determine the fractal dimension of the cellular structures. Comparison of these metrics is presented for different cell lines with varying phenotypes including MCF10A, MCF7, BT474, HT-29, A431, and A549 cell lines, in addition to three cell populations with modified phenotypes. Our results show that disorder strength and fractal dimension can both be obtained with quantitative phase imaging and that these metrics can independently distinguish between different cell lines. Furthermore, their combined use presents a new approach for better understanding cellular restructuring during different pathways.

Transforming growth factor-β1, -β2, and -β3 (TGF-β1, -β2, and -β3) are secreted signaling ligands that play essential roles in tissue development, tissue maintenance, immune response, and wound healing. TGF-β ligands form homodimers and signal by assembling a heterotetrameric receptor complex comprised of two type I receptor (TβRI):type II receptor (TβRII) pairs. TGF-β1 and TGF-β3 ligands signal with high potency due to their high affinity for TβRII, which engenders high-affinity binding of TβRI through a composite TGF-β:TβRII binding interface. However, TGF-β2 binds TβRII 200-500 more weakly than TGF-β1 and TGF-β3 and signals with lower potency compared with these ligands. Remarkably, the presence of an additional membrane-bound coreceptor, known as betaglycan, increases TGF-β2 signaling potency to levels similar to TGF-β1 and -β3. The mediating effect of betaglycan occurs even though it is displaced from and not present in the heterotetrameric receptor complex through which TGF-β2 signals. Published biophysics studies have experimentally established the kinetic rates of the individual ligand-receptor and receptor-receptor interactions that initiate heterotetrameric receptor complex assembly and signaling in the TGF-β system; however, current experimental approaches are not able to directly measure kinetic rates for the intermediate and latter steps of assembly. To characterize these steps in the TGF-β system and determine the mechanism of betaglycan in the potentiation of TGF-β2 signaling, we developed deterministic computational models with different modes of betaglycan binding and varying cooperativity between receptor subtypes. The models identified conditions for selective enhancement of TGF-β2 signaling. The models provide support for additional receptor binding cooperativity that has been hypothesized but not evaluated in the literature. The models further showed that betaglycan binding to the TGF-β2 ligand through two domains provides an effective mechanism for transfer to the signaling receptors that has been tuned to efficiently promote assembly of the TGF-β2(TβRII)2(TβRI)2 signaling complex.

Mitochondrial inner membrane potentials in cardiomyocytes may oscillate in cycles of depolarization/repolarization when the mitochondrial network is exposed to metabolic or oxidative stress. The frequencies of such oscillations are dynamically changing while clusters of weakly coupled mitochondrial oscillators adjust to a common phase and frequency. Across the cardiac myocyte, the averaged signal of the mitochondrial population follows self-similar or fractal dynamics; however, fractal properties of individual mitochondrial oscillators have not yet been examined. We show that the largest synchronously oscillating cluster exhibits a fractal dimension, D, that is indicative of self-similar behavior with D=1.27±0.11, in contrast to the remaining network mitochondria whose fractal dimension is close to that of Brownian noise, D=1.58±0.10. We further demonstrate that fractal behavior is correlated with local coupling mechanisms, whereas it is only weakly linked to measures of functional connections between mitochondria. Our findings suggest that individual mitochondrial fractal dimensions may serve as a simple measure of local mitochondrial coupling.

[FeFe] hydrogenases are enzymes that have acquired a unique capacity to synthesize or consume molecular hydrogen (H2). This function relies on a complex catalytic mechanism involving the active site and two distinct electron and proton transfer networks working in concert. By an analysis based on terahertz vibrations of [FeFe] hydrogenase structure, we are able to predict and identify the existence of rate-promoting vibrations at the catalytic site and the coupling with functional residues involved in reported electron and proton transfer networks. Our findings suggest that the positioning of the cluster is influenced by the response of the scaffold to thermal fluctuations, which in turn drives the formation of networks for electron transfer through phonon-assisted mechanisms. Thus, we address the problem of linking the molecular structure to the catalytic function through picosecond dynamics, while raising the functional gain brought by the cofactors or clusters, using the concept of fold-encoded localized vibrations.

The functional properties of some biological ion channels and membrane transport proteins are proposed to exploit anion-hydrophobic interactions. Here, we investigate a chloride-pumping rhodopsin as an example of a membrane protein known to contain a defined anion binding site composed predominantly of hydrophobic residues. Using molecular dynamics simulations, we explore Cl- binding to this hydrophobic site and compare the dynamics arising when electronic polarization is neglected (CHARMM36 [c36] fixed-charge force field), included implicitly (via the prosECCo force field), or included explicitly (through the polarizable force field, AMOEBA). Free energy landscapes of Cl- moving out of the binding site and into bulk solution demonstrate that the inclusion of polarization results in stronger ion binding and a second metastable binding site in chloride-pumping rhodopsin. Simulations focused on this hydrophobic binding site also indicate longer binding durations and closer ion proximity when polarization is included. Furthermore, simulations reveal that Cl- within this binding site interacts with an adjacent loop to facilitate rebinding events that are not observed when polarization is neglected. These results demonstrate how the inclusion of polarization can influence the behavior of anions within protein binding sites and can yield results comparable with more accurate and computationally demanding methods.

Residual force enhancement (RFE), an increase in isometric force after active stretching of a muscle compared with the purely isometric force at the corresponding length, has been consistently observed throughout the structural hierarchy of skeletal muscle. Similar to RFE, passive force enhancement (PFE) is also observable in skeletal muscle and is defined as an increase in passive force when a muscle is deactivated after it has been actively stretched compared with the passive force following deactivation of a purely isometric contraction. These history-dependent properties have been investigated abundantly in skeletal muscle, but their presence in cardiac muscle remains unresolved and controversial. The purpose of this study was to investigate whether RFE and PFE exist in cardiac myofibrils and whether the magnitudes of RFE and PFE increase with increasing stretch magnitudes. Cardiac myofibrils were prepared from the left ventricles of New Zealand White rabbits, and the history-dependent properties were tested at three different final average sarcomere lengths (n = 8 for each), 1.8, 2, and 2.2 μm, while the stretch magnitude was kept at 0.2 μm/sarcomere. The same experiment was repeated with a final average sarcomere length of 2.2 μm and a stretching magnitude of 0.4 μm/sarcomere (n = 8). All 32 cardiac myofibrils exhibited increased forces after active stretching compared with the corresponding purely isometric reference conditions (p < 0.05). Furthermore, the magnitude of RFE was greater when myofibrils were stretched by 0.4 compared with 0.2 μm/sarcomere (p < 0.05). We conclude that, like in skeletal muscle, RFE and PFE are properties of cardiac myofibrils and are dependent on stretch magnitude.

RNA pseudoknots are a kind of important tertiary motif, and the structures and stabilities of pseudoknots are generally critical to the biological functions of RNAs with the motifs. In this work, we have carefully refined our previously developed coarse-grained model with salt effect through involving a new coarse-grained force field and a replica-exchange Monte Carlo algorithm, and employed the model to predict structures and stabilities of complex RNA pseudoknots in ion solutions beyond minimal H-type pseudoknots. Compared with available experimental data, the newly refined model can successfully predict 3D structures from sequences for the complex RNA pseudoknots including SARS-CoV-2 programming-1 ribosomal frameshifting element and Zika virus xrRNA, and can reliably predict the thermal stabilities of RNA pseudoknots with various sequences and lengths over broad ranges of monovalent/divalent salts. In addition, for complex pseudoknots including SARS-CoV-2 frameshifting element, our analyses show that their thermally unfolding pathways are mainly dependent on the relative stabilities of unfolded intermediate states, in analogy to those of minimal H-type pseudoknots.

Understanding how cells remember previous mechanical environments to influence their fate, or mechanical memory, informs the design of biomaterials and therapies in medicine. Current regeneration therapies, such as cartilage regeneration procedures, require 2D cell expansion processes to achieve large cell populations critical for the repair of damaged tissues. However, the limit of mechanical priming for cartilage regeneration procedures before inducing long-term mechanical memory following expansion processes is unknown, and mechanisms defining how physical environments influence the therapeutic potential of cells remain poorly understood. Here, we identify a threshold to mechanical priming separating reversible and irreversible effects of mechanical memory. After 16 population doublings in 2D culture, expression levels of tissue-identifying genes in primary cartilage cells (chondrocytes) are not recovered when transferred to 3D hydrogels, while expression levels of these genes were recovered for cells only expanded for eight population doublings. Additionally, we show that the loss and recovery of the chondrocyte phenotype correlates with a change in chromatin architecture, as shown by structural remodeling of the trimethylation of H3K9. Efforts to disrupt the chromatin architecture by suppressing or increasing levels of H3K9me3 reveal that only with increased levels of H3K9me3 did the chromatin architecture of the native chondrocyte phenotype partially return, along with increased levels of chondrogenic gene expression. These results further support the connection between the chondrocyte phenotype and chromatin architecture, and also reveal the therapeutic potential of inhibitors of epigenetic modifiers as disruptors of mechanical memory when large numbers of phenotypically suitable cells are required for regeneration procedures.

The distribution of red blood cells (RBCs) in the microcirculation determines the oxygen delivery and solute transport to tissues. This process relies on the partitioning of RBCs at successive bifurcations throughout the microvascular network, and it has been known since the last century that RBCs partition disproportionately to the fractional blood flow rate, therefore leading to heterogeneity of the hematocrit (i.e., volume fraction of RBCs in blood) in microvessels. Usually, downstream of a microvascular bifurcation, the vessel branch with a higher fraction of blood flow receives an even higher fraction of RBC flux. However, both temporal and time-average deviations from this phase-separation law have been observed in recent studies. Here, we quantify how the microscopic behavior of RBC lingering (i.e., RBCs temporarily residing near the bifurcation apex with diminished velocity) influences their partitioning, through combined in vivo experiments and in silico simulations. We developed an approach to quantify the cell lingering at highly confined capillary-level bifurcations and demonstrate that it correlates with deviations of the phase-separation process from established empirical predictions by Pries et al. Furthermore, we shed light on how the bifurcation geometry and cell membrane rigidity can affect the lingering behavior of RBCs; e.g., rigid cells tend to linger less than softer ones. Taken together, RBC lingering is an important mechanism that should be considered when studying how abnormal RBC rigidity in diseases such as malaria and sickle-cell disease could hinder the microcirculatory blood flow or how the vascular networks are altered under pathological conditions (e.g., thrombosis, tumors, aneurysm).

The stress-free state (SFS) of red blood cells (RBCs) is a fundamental reference configuration for the calibration of computational models, yet it remains unknown. Current experimental methods cannot measure the SFS of cells without affecting their mechanical properties, whereas computational postulates are the subject of controversial discussions. Here, we introduce data-driven estimates of the SFS shape and the visco-elastic properties of RBCs. We employ data from single-cell experiments that include measurements of the equilibrium shape of stretched cells and relaxation times of initially stretched RBCs. A hierarchical Bayesian model accounts for these experimental and data heterogeneities. We quantify, for the first time, the SFS of RBCs and use it to introduce a transferable RBC (t-RBC) model. The effectiveness of the proposed model is shown on predictions of unseen experimental conditions during the inference, including the critical stress of transitions between tumbling and tank-treading cells in shear flow. Our findings demonstrate that the proposed t-RBC model provides predictions of blood flows with unprecedented accuracy and quantified uncertainties.

Increased blood viscosity in type 2 diabetes mellitus (T2DM) is a risk factor for the development of insulin resistance and diabetes-related vascular complications; however, individuals with T2DM exhibit heterogeneous hemorheological properties, including cell deformation and aggregation. Using a multiscale red blood cell (RBC) model with key parameters derived from patient-specific data, we present a computational study of the rheological properties of blood from individual patients with T2DM. Specifically, one key model parameter, which determines the shear stiffness of the RBC membrane (μ) is informed by the high-shear-rate blood viscosity of patients with T2DM. At the same time, the other, which contributes to the strength of the RBC aggregation interaction (D0), is derived from the low-shear-rate blood viscosity of patients with T2DM. The T2DM RBC suspensions are simulated at different shear rates, and the predicted blood viscosity is compared with clinical laboratory-measured data. The results show that the blood viscosity obtained from clinical laboratories and computational simulations are in agreement at both low and high shear rates. These quantitative simulation results demonstrate that the patient-specific model has truly learned the rheological behavior of T2DM blood by unifying the mechanical and aggregation factors of the RBCs, which provides an effective way to extract quantitative predictions of the rheological properties of the blood of individual patients with T2DM.

α-Helical secondary structures impart specific mechanical and physiochemical properties to peptides and proteins, enabling them to perform a vast array of molecular tasks ranging from membrane insertion to molecular allostery. Loss of α-helical content in specific regions can inhibit native protein function or induce new, potentially toxic, biological activities. Thus, identifying specific residues that exhibit loss or gain of helicity is critical for understanding the molecular basis of function. Two-dimensional infrared (2D IR) spectroscopy coupled with isotope labeling is capable of capturing detailed structural changes in polypeptides. Yet, questions remain regarding the inherent sensitivity of isotope-labeled modes to local changes in α-helicity, such as terminal fraying; the origin of spectral shifts (hydrogen-bonding versus vibrational coupling); and the ability to definitively detect coupled isotopic signals in the presence of overlapping side chains. Here, we address each of these points individually by characterizing a short, model α-helix (DPAEAAKAAAGR-NH2) with 2D IR and isotope labeling. These results demonstrate that pairs of 13C18O probes placed three residues apart can detect subtle structural changes and variations along the length of the model peptide as the α-helicity is systematically tuned. Comparison of singly and doubly labeled peptides affirm that frequency shifts arise primarily from hydrogen-bonding, while vibrational coupling between paired isotopes leads to increased peak areas that can be clearly differentiated from underlying side-chain modes or uncoupled isotope labels not participating in helical structures. These results demonstrate that 2D IR in tandem with i,i+3 isotope-labeling schemes can capture residue-specific molecular interactions within a single turn of an α-helix.

Despite the molecular evidence that a nearly linear steady-state current-voltage relationship in mammalian astrocytes reflects a total current resulting from more than one differentially regulated K+ conductance, detailed ordinary differential equation (ODE) models of membrane voltage Vm are still lacking. Various experimental results reporting altered rectification of the major Kir currents in glia, dominated by Kir4.1, have motivated us to develop a detailed model of Vm dynamics incorporating the weaker potassium K2P-TREK1 current in addition to Kir4.1, and study the stability of the resting state Vr. The main question is whether, with the loss of monotonicity in glial I-V curve resulting from altered Kir rectification, the nominal resting state Vr remains stable, and the cell retains the trivial, potassium electrode behavior with Vm after EK. The minimal two-dimensional model of Vm near Vr showed that an N-shape deformed Kir I-V curve induces multistability of Vm in a model that incorporates K2P activation kinetics, and nonspecific K+ leak currents. More specifically, an asymmetrical, nonlinear decrease of outward Kir4.1 conductance, turning the channels into inward rectifiers, introduces instability of Vr. That happens through a robust bifurcation giving birth to a second, more depolarized stable resting state Vdr > -10 mV. Realistic recordings from electrographic seizures were used to perturb the model. Simulations of the model perturbed by constant current through gap junctions and seizure-like discharges as local field potentials led to depolarization and switching of Vm between the two stable states, in a downstate-upstate manner. In the event of prolonged depolarizations near Vdr, such catastrophic instability would affect all aspects of the glial function, from metabolic support to membrane transport, and practically all neuromodulatory roles assigned to glia.