Search

Many methods to accelerate sampling of molecular configurations are based on the idea that temperature can be used to accelerate rare transitions. These methods typically compute equilibrium properties at a target temperature using reweighting or through Monte Carlo exchanges between replicas at higher temperatures. A recent paper [G. M. Rotskoff and E. Vanden-Eijnden, Phys. Rev. Lett. 122, 150602 (2019)] demonstrated that accurate equilibrium densities of states can also be computed through a nonequilibrium "quench" process, where sampling is performed at a higher temperature to encourage rapid mixing and then quenched to lower energy states with dissipative dynamics. Here, we provide an implementation of the quench dynamics in LAMMPS and evaluate a new formulation of nonequilibrium estimators for the computation of partition functions or free energy surfaces (FESs) of molecular systems. We show that the method is exact for a minimal model of N-independent harmonic springs and use these analytical results to develop heuristics for the amount of quenching required to obtain accurate sampling. We then test the quench approach on alanine dipeptide, where we show that it gives an FES that is accurate near the most stable configurations using the quench approach but disagrees with a reference umbrella sampling calculation in high FE regions. We then show that combining quenching with umbrella sampling allows the efficient calculation of the free energy in all regions. Moreover, by using this combined scheme, we obtain the FES across a range of temperatures at no additional cost, making it much more efficient than standard umbrella sampling if this information is required. Finally, we discuss how this approach can be extended to solute tempering and demonstrate that it is highly accurate for the case of solvated alanine dipeptide without any additional modifications. [ABSTRACT FROM AUTHOR]

We evaluate neural network (NN) coarse-grained (CG) force fields compared to traditional CG molecular mechanics force fields. We conclude that NN force fields are able to extrapolate and sample from unseen regions of the free energy surface when trained with limited data. Our results come from 88 NN force fields trained on different combinations of clustered free energy surfaces from four protein mapped trajectories. We used a statistical measure named total variation similarity to assess the agreement between reference free energy surfaces from mapped atomistic simulations and CG simulations from trained NN force fields. Our conclusions support the hypothesis that NN CG force fields trained with samples from one region of the proteins' free energy surface can, indeed, extrapolate to unseen regions. Additionally, the force matching error was found to only be weakly correlated with a force field's ability to reconstruct the correct free energy surface. [ABSTRACT FROM AUTHOR]

From rock salt to nanoparticle superlattices, complex structure can emerge from simple building blocks that attract each other through Coulombic forces.sup.1-4. On the micrometre scale, however, colloids in water defy the intuitively simple idea of forming crystals from oppositely charged partners, instead forming non-equilibrium structures such as clusters and gels.sup.5-7. Although various systems have been engineered to grow binary crystals.sup.8-11, native surface charge in aqueous conditions has not been used to assemble crystalline materials. Here we form ionic colloidal crystals in water through an approach that we refer to as polymer-attenuated Coulombic self-assembly. The key to crystallization is the use of a neutral polymer to keep particles separated by well defined distances, allowing us to tune the attractive overlap of electrical double layers, directing particles to disperse, crystallize or become permanently fixed on demand. The nucleation and growth of macroscopic single crystals is demonstrated by using the Debye screening length to fine-tune assembly. Using a variety of colloidal particles and commercial polymers, ionic colloidal crystals isostructural to caesium chloride, sodium chloride, aluminium diboride and K.sub.4C.sub.60 are selected according to particle size ratios. Once fixed by simply diluting out solution salts, crystals are pulled out of the water for further manipulation, demonstrating an accurate translation from solution-phase assembly to dried solid structures. In contrast to other assembly approaches, in which particles must be carefully engineered to encode binding information.sup.12-18, polymer-attenuated Coulombic self-assembly enables conventional colloids to be used as model colloidal ions, primed for crystallization. Oppositely charged colloidal particles are assembled in water through an approach that allows electrostatic interactions to be precisely tuned to generate macroscopic single crystals., Author(s): Theodore Hueckel [sup.1] , Glen M. Hocky [sup.1] , Jeremie Palacci [sup.2] , Stefano Sacanna [sup.1] Author Affiliations: (1) Department of Chemistry, New York University, New York, USA (2) [...]

Protein–ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent "bonds" being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed—slip bonds, where the unbinding rate increases, and catch bonds, where the unbinding rate decreases. Despite much theoretical work, we cannot predict for which protein–ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin–streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite the use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predicts an overall slip bond for larger forces consistent with prior experimental results although there are substantial deviations at small forces that require further investigation. This work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future. [ABSTRACT FROM AUTHOR]

Many proteins in cells are capable of sensing and responding to piconewton-scale forces, a regime in which conformational changes are small but significant for biological processes. In order to efficiently and effectively sample the response of these proteins to small forces, enhanced sampling techniques will be required. In this work, we derive, implement, and evaluate an efficient method to simultaneously sample the result of applying any constant pulling force within a specified range to a molecular system of interest. We start from simulated tempering in force, whereby force is added as a linear bias on a collective variable to the system's Hamiltonian, and the coefficient is taken as a continuous auxiliary degree of freedom. We derive a formula for an average collective-variable-dependent force, which depends on a set of weights learned on-the-fly throughout a simulation, that reflect the limit where force varies infinitely quickly. Simulation data can then be used to retroactively compute averages of any observable at any force within the specified range. This technique is based on recent work deriving similar equations for infinite switch simulated tempering in temperature, which showed that the infinite switch limit is the most efficient for sampling. Here, we demonstrate that our method accurately samples molecular systems at all forces within a user defined force range simultaneously and show how it can serve as an enhanced sampling tool for cases where the pulling direction destabilizes states that have low free-energy at zero-force. This method is implemented in and freely distributed with the PLUMED open-source sampling library, and hence can be readily applied to problems using a wide range of molecular dynamics software packages. [ABSTRACT FROM AUTHOR]

Accounting for electrons and nuclei simultaneously is a powerful capability of ab initio molecular dynamics (AIMD). However, AIMD is often unable to accurately reproduce properties of systems such as water due to inaccuracies in the underlying electronic density functionals. This shortcoming is often addressed by added empirical corrections and/or increasing the simulation temperature. We present here a maximum-entropy approach to directly incorporate limited experimental data via a minimal bias. Biased AIMD simulations of water and an excess proton in water are shown to give significantly improved properties both for observables which were biased to match experimental data and for unbiased observables. This approach also yields new physical insight into inaccuracies in the underlying density functional theory as utilized in the unbiased AIMD. [ABSTRACT FROM AUTHOR]

Arp2/3 complex nucleates branched actin filaments that provide pushing forces to drive cellular processes such as lamellipodial protrusion and endocytosis. Arp2/3 complex is intrinsically inactive, and multiple classes of nucleation promoting factors (NPFs) stimulate its nucleation activity. When activated by WASP family NPFs, the complex must bind to the side of a preexisting (mother) filament of actin to complete the nucleation process, ensuring that WASP-mediated activation creates branched rather than linear actin filaments. How actin filaments contribute to activation is currently not understood, largely due to the lack of high-resolution structures of activated Arp2/3 complex bound to the side of a filament. Here, we present the 3.9-Å cryo-electron microscopy structure of the Arp2/3 complex at a branch junction. The structure reveals contacts between Arp2/3 complex and the side of the mother actin filament that likely stimulate subunit flattening, a conformational change that allows the actin-related protein subunits in the complex (Arp2 and Arp3) to mimic filamentous actin subunits. In contrast, limited contact between the bottom half of the complex and the mother filament suggests that clamp twisting, a second major conformational change observed in the active state, is not stimulated by actin filaments, potentially explaining why actin filaments are required but insufficient to trigger nucleation during WASP-mediated activation. Along with biochemical and live-cell imaging data and molecular dynamics simulations, the structure reveals features critical for the interaction of Arp2/3 complex with actin filaments and regulated assembly of branched actin filament networks in cells. [ABSTRACT FROM AUTHOR]

Understanding the role of nonequilibrium driving in self-organization is crucial for developing a predictive description of biological systems, yet it is impeded by their complexity. The actin cytoskeleton serves as a paradigm for how equilibrium and nonequilibrium forces combine to give rise to self-organization. Motivated by recent experiments that show that actin filament growth rates can tune the morphology of a growing actin bundle cross-linked by two competing types of actin-binding proteins [S. L. Freedman et al., Proc. Natl. Acad. Sci. U.S.A. 116, 16192-16197 (2019)], we construct a minimal model for such a system and show that the dynamics of a growing actin bundle are subject to a set of thermodynamic constraints that relate its nonequilibrium driving, morphology, and molecular fluxes. The thermodynamic constraints reveal the importance of correlations between these molecular fluxes and offer a route to estimating microscopic driving forces from microscopy experiments. [ABSTRACT FROM AUTHOR]

Ultrastable glasses have risen to prominence due to their potentially useful material properties and the tantalizing possibility of a general method of preparation via vapor deposition. Despite the importance of this novel class of amorphous materials, numerical studies have been scarce because achieving ultrastability in atomistic simulations is an enormous challenge. Here, we bypass this difficulty and establish that randomly pinning the position of a small fraction of particles inside an equilibrated supercooled liquid generates ultrastable configurations at essentially no numerical cost, while avoiding undesired structural changes due to the preparation protocol. Building on the analogy with vapor-deposited ultrastable glasses, we study the melting kinetics of these configurations following a sudden temperature jump into the liquid phase. In homogeneous geometries, we find that enhanced kinetic stability is accompanied by large scale dynamic heterogeneity, while a competition between homogeneous and heterogeneous melting is observed when a liquid boundary invades the glass at constant velocity. Our work demonstrates the feasibility of large-scale, atomistically resolved, and experimentally relevant simulations of the kinetics of ultrastable glasses. [ABSTRACT FROM AUTHOR]

Members of the ADF/cofilin family of regulatory proteins bind actin filaments cooperatively, locally change actin subunit conformation and orientation, and sever filaments at "boundaries" between bare and cofilin-occupied segments. A cluster of bound cofilin introduces two distinct classes of boundaries due to the intrinsic polarity of actin filaments, one at the "pointed" end side and the other at the "barbed" end-side of the cluster; severing occurs more readily at the pointed end side of the cluster ("fast-severing" boundary) than the barbed end side ("slow-severing" boundary). A recent electron-cryomicroscopy (cryo-EM) model of the slow-severing boundary revealed structural "defects" at the interface that potentially contribute to severing. However, the structure of the fast-severing boundary remains uncertain. Here, we use extensive molecular dynamics simulations to produce atomic resolution models of both severing boundaries. Our equilibrated simulation model of the slow-severing boundary is consistent with the cryo-EM structural model. Simulations indicate that actin subunits at both boundaries adopt structures intermediate between those of bare and cofilin-bound actin subunits. These "intermediate" states have compromised intersubunit contacts, but those at the slow-severing boundary are stabilized by cofilin bridging interactions, accounting for its lower fragmentation probability. Simulations where cofilin proteins are removed from cofilactin filaments favor a mechanism in which a cluster of two contiguously bound cofilins is needed to fully stabilize the cofilactin conformation, promote cooperative binding interactions, and accelerate filament severing. Together, these studies provide a molecular-scale foundation for developing coarse-grained and theoretical descriptions of cofilin-mediated actin filament severing. [ABSTRACT FROM AUTHOR]

Many cytoskeletal proteins perform fundamental biological processes and are evolutionarily ancient. For example, the superfamily of actin-related proteins (Arps) specialized early in eukaryotic evolution for diverse cellular roles in the cytoplasm and the nucleus. Despite its strict conservation across eukaryotes, we find that the Arp superfamily has undergone dramatic lineage-specific diversification in Drosophila. Our phylogenomic analyses reveal four independent Arp gene duplications that occurred in the common ancestor of the obscura group of Drosophila and have been mostly preserved in this lineage. All four obscura -specific Arp paralogs are predominantly expressed in the male germline and have evolved under positive selection. We focus our analyses on the divergent Arp2D paralog, which arose via a retroduplication event from Arp2 , a component of the Arp2/3 complex that polymerizes branched actin networks. Computational modeling analyses suggest that Arp2D can replace Arp2 in the Arp2/3 complex and bind actin monomers. Together with the signature of positive selection, our findings suggest that Arp2D may augment Arp2's functions in the male germline. Indeed, we find that Arp2D is expressed during and following male meiosis, where it localizes to distinct locations such as actin cones—specialized cytoskeletal structures that separate bundled spermatids into individual mature sperm. We hypothesize that this unprecedented burst of genetic innovation in cytoskeletal proteins may have been driven by the evolution of sperm heteromorphism in the obscura group of Drosophila. [ABSTRACT FROM AUTHOR]

Cytokinesis physically separates dividing cells by forming a contractile actomyosin ring. The fission yeast contractile ring has been proposed to assemble by Search-Capture-Pull-Release from cytokinesis precursor nodes that include the molecular motor type-II myosin Myo2 and the actin assembly factor formin Cdc12. By successfully reconstituting Search-Capture-Pull in vitro, we discovered that formin Cdc12 is a mechanosensor, whereby myosin pulling on formin-bound actin filaments inhibits Cdc12-mediated actin assembly. We mapped Cdc12 mechanoregulation to its formin homology 1 domain, which facilitates delivery of new actin subunits to the elongating actin filament. Quantitative modeling suggests that the pulling force of the myosin propagates through the actin filament, which behaves as an entropic spring, and thereby may stretch the disordered formin homology 1 domain and impede formin-mediated actin filament elongation. Finally, live cell imaging of mechano-insensitive formin mutant cells established that mechanoregulation of formin Cdc12 is required for efficient contractile ring assembly in vivo. [ABSTRACT FROM AUTHOR]

The fission yeast actin cytoskeleton is an ideal, simplified system to investigate fundamental mechanisms behind cellular self-organization. By focusing on the stabilizing protein tropomyosin Cdc8, bundling protein fimbrin Fim1, and severing protein coffin Adf1, we examined how their pairwise and collective interactions with actin filaments regulate their activity and segregation to functionally diverse F-actin networks. Utilizing multi-color TIRF microscopy of in vitro reconstituted F-actin networks, we observed and characterized two distinct Cdc8 cables loading and spreading cooperatively on individual actin filaments. Furthermore, Cdc8, Fim1, and Adf1 all compete for association with F-actin by different mechanisms, and their cooperative association with actin filaments affects their ability to compete. Finally, competition between Fim1 and Adf1 for F-actin synergizes their activities, promoting rapid displacement of Cdc8 from a dense F-actin network. Our findings reveal that competitive and cooperative interactions between actin binding proteins help define their associations with different F-actin networks. [ABSTRACT FROM AUTHOR]

Ions regulate the assembly and mechanical properties of actin filaments. Recent work using structural bioinformatics and site-specific mutagenesis favors the existence of two discrete and specific divalent cation binding sites on actin filaments, positioned in the long axis between actin subunits. Cation binding at one site drives polymerization, while the other modulates filament stiffness and plays a role in filament severing by the regulatory protein, cofilin. Existing structural methods have not been able to resolve filament-associated cations, and so in this work we turn to molecular dynamics simulations to suggest a candidate binding pocket geometry for each site and to elucidate the mechanism by which occupancy of the "stiffness site" affects filament mechanical properties. Incorporating a magnesium ion in the "polymerization site" does not seem to require any large-scale change to an actin subunit's conformation. Binding of a magnesium ion in the "stiffness site" adheres the actin DNase-binding loop (D-loop) to its long-axis neighbor, which increases the filament torsional stiffness and bending persistence length. Our analysis shows that bound D-loops occupy a smaller region of accessible conformational space. Cation occupancy buries key conserved residues of the D-loop, restricting accessibility to regulatory proteins and enzymes that target these amino acids. [ABSTRACT FROM AUTHOR]

We investigate the connection between local structure and dynamical heterogeneity in supercooled liquids. Through the study of four different models, we show that the correlation between a particle's mobility and the degree of local order in nearby regions is highly system dependent. Our results suggest that the correlation between local structure and dynamics is weak or absent in systems that conform well to the mean-field picture of glassy dynamics and strong in those that deviate from this paradigm. Finally, we investigate the role of order-agnostic point-to-set correlations and reveal that they provide similar information content to local structure measures, at least in the system where local order is most pronounced. [ABSTRACT FROM AUTHOR]

We study the relaxation dynamics of a binary Lennard-Jones liquid in the presence of an amorphous wall generated from equilibrium particle configurations. In qualitative agreement with the results presented by Kob et al. [Nat. Phys. 8. 164 (2012).] for a liquid of harmonic spheres, we find that our binary mixture shows a saturation of the dynamical length scale close to the mode-coupling temperature Tc. Furthermore we show that, due to the broken symmetry imposed by the wall, signatures of an additional change in dynamics become apparent at a temperature well above Tc. We provide evidence that this modification in the relaxation dynamics occurs at a recently proposed dynamical crossover temperature Ts > Tc, which is related to the breakdown of the Stokes-Einstein relation. We find that this dynamical crossover at Ts is also observed for the harmonic spheres as well as a WCA liquid, showing that it may be a general feature of glass-forming systems. [ABSTRACT FROM AUTHOR]

It has been demonstrated recently that supercooled liquids sharing simple structural features (e.g. pair distribution functions) may exhibit strikingly distinct dynamical behavior. Here we show that a more subtle structural feature correlates with relaxation times in three simulated systems that have nearly identical radial distribution functions but starkly different dynamical behavior. In particular, for the first time we determine the thermodynamic "point-to-set" length scale in several canonical model systems and demonstrate the quantitative connection between this length scale and the growth of relaxation times. Our results provide clues necessary for distinguishing competing theories of the glass transition. [ABSTRACT FROM AUTHOR]

In this work, we study the nature of transitions between inherent structures of a two-dimensional model supercooled liquid. We demonstrate that these transitions occur predominately along a small number of directions on the energy landscape. Moreover, we show that the number of such directions decreases as the temperature of the liquid is decreased in the mildly supercooled regime, in concert with earlier studies on an athermal jamming system. We show that this decrease happens in parallel with a change in character of the transitions as dynamics in the system become more heterogeneous and localized. We investigate the origin of these trends, which suggests interesting connections between jamming and thermal glassy phenomena. [ABSTRACT FROM AUTHOR]

Summary Cells assemble and maintain functionally distinct actin cytoskeleton networks with various actin filament organizations and dynamics through the coordinated action of different sets of actin-binding proteins. The biochemical and functional properties of diverse actin-binding proteins, both alone and in combination, have been increasingly well studied. Conversely, how different sets of actin-binding proteins properly sort to distinct actin filament networks in the first place is not nearly as well understood. Actin-binding protein sorting is critical for the self-organization of diverse dynamic actin cytoskeleton networks within a common cytoplasm. Using in vitro reconstitution techniques including biomimetic assays and single-molecule multi-color total internal reflection fluorescence microscopy, we discovered that sorting of the prominent actin-bundling proteins fascin and α-actinin to distinct networks is an intrinsic behavior, free of complicated cellular signaling cascades. When mixed, fascin and α-actinin mutually exclude each other by promoting their own recruitment and inhibiting recruitment of the other, resulting in the formation of distinct fascin- or α-actinin-bundled domains. Subdiffraction-resolution light microscopy and negative-staining electron microscopy revealed that fascin domains are densely packed, whereas α-actinin domains consist of widely spaced parallel actin filaments. Importantly, other actin-binding proteins such as fimbrin and espin show high specificity between these two bundle types within the same reaction. Here we directly observe that fascin and α-actinin intrinsically segregate to discrete bundled domains that are specifically recognized by other actin-binding proteins. [ABSTRACT FROM AUTHOR]

The cell interior is packed with macromolecules of mesoscale size, and this crowded milieu significantly influences cellular physiology. Cellular stress responses almost universally lead to inhibition of translation, resulting in polysome collapse and release of mRNA. The released mRNA molecules condense with RNA-binding proteins to form ribonucleoprotein (RNP) condensates known as processing bodies and stress granules. Here, we show that polysome collapse and condensation of RNA transiently fluidize the cytoplasm, and coarse-grained molecular dynamic simulations support this as a minimal mechanism for the observed biophysical changes. Increased mesoscale diffusivity correlates with the efficient formation of quality control bodies (Q-bodies), membraneless organelles that compartmentalize misfolded peptides during stress. Synthetic, light-induced RNA condensation also fluidizes the cytoplasm. Together, our study reveals a functional role for stress-induced translation inhibition and formation of RNP condensates in modulating the physical properties of the cytoplasm to enable efficient response of cells to stress conditions. [Display omitted] • Diverse stresses cause transient fluidization of the cytoplasm • Polysome collapse is required for cytoplasmic fluidization • mRNA condensation into stress granules or P-bodies is also required • Cytoplasmic fluidization facilitates the formation of new mesoscale structures Cells must reorganize when exposed to stress. Xie et al. discover that the cytoplasm is initially fluidized in all stress conditions they test. This biophysical change requires polysomes to disassemble and mRNA to condense into P-bodies or stress granules. Fluidization allows for the efficient formation of new stress-induced structures. [ABSTRACT FROM AUTHOR]

Many biological processes, including cell division, growth, and motility, rely on rapid remodeling of the actin cytoskeleton and on actin filament severing by the regulatory protein cofilin. Phosphorylation of vertebrate cofilin at Ser-3 regulates both actin binding and severing. Substitution of serine with aspartate at position 3 (S3D) is widely used to mimic cofilin phosphorylation in cells and in vitro. The S3D substitution weakens cofilin binding to filaments, and it is presumed that subsequent reduction in cofilin occupancy inhibits filament severing, but this hypothesis has remained untested. Here, using time-resolved phosphorescence anisotropy, electron cryomicroscopy, and allatom molecular dynamics simulations, we show that S3D cofilin indeed binds filaments with lower affinity, but also with a higher cooperativity than wild-type cofilin, and severs actin weakly across a broad range of occupancies.Wefound that three factors contribute to the severing deficiency of S3D cofilin. First, the high cooperativity of S3D cofilin generates fewer boundaries between bare and decorated actin segments where severing occurs preferentially. Second, S3D cofilin only weakly alters filament bending and twisting dynamics and therefore does not introduce the mechanical discontinuities required for efficient filament severing at boundaries. Third, Ser-3 modification (i.e. substitution with Asp or phosphorylation) "undocks" and repositions the cofilinNterminus away from the filament axis, which compromises S3D cofilin's ability to weaken longitudinal filament subunit interactions. Collectively, our results demonstrate that, in addition to inhibiting actin binding, Ser-3 modification favors formation of a cofilin-binding mode that is unable to sufficiently alter filament mechanical properties and promote severing. [ABSTRACT FROM AUTHOR]

Vinculin forms a catch bond with the cytoskeletal polymer actin, displaying an increased bond lifetime upon force application. Notably, this behavior depends on the direction of the applied force, which has significant implications for cellular mechanotransduction. In this study, we present a comprehensive molecular dynamics simulation study, employing enhanced sampling techniques to investigate the thermodynamic, kinetic, and mechanistic aspects of this phenomenon at physiologically relevant forces. We dissect a catch bond mechanism in which force shifts vinculin between either a weakly- or strongly-bound state. Our results demonstrate that models for these states have unbinding times consistent with those from single-molecule studies, and suggest that both have some intrinsic catch bonding behavior. We provide atomistic insight into this behavior, and show how a directional pulling force can promote the strong or weak state. Crucially, our strategy can be extended to capture the difficult-to-capture effects of small mechanical forces on biomolecular systems in general, and those involved in mechanotransduction more specifically.

The outer membrane is the defining structure of Gram-negative bacteria. We previously demonstrated that it is critical for the mechanical integrity of the cell envelope and therefore to the robustness of the bacterial cell as a whole. Here, to determine the key molecules and moieties within the outer membrane that underlie its contribution to cell envelope mechanics, we measured cell-envelope stiffness across several sets of mutants with altered outer-membrane sugar content, protein content, and electric charge. To decouple outer membrane stiffness from total cell envelope stiffness, we developed a novel microfluidics-based "osmotic force extension" assay. In tandem, we developed a simple method to increase throughput of microfluidics experiments by performing them on color-coded pools of mutants. Using Escherichia coli as a model Gram-negative bacterium, we found that truncating the core oligosaccharide, deleting the β-barrel protein OmpA, or deleting lipoprotein outer membrane-cell wall linkers all had the same modest, convergent effect on total cell-envelope stiffness but had large, varying effects on the ability of the cell wall to transfer tension to the outer membrane during large hyperosmotic shocks. Surprisingly, altering lipid A charge had little effect on the mechanical properties of the envelope. Importantly, the presence or absence of OmpA determined whether truncating the core oligosaccharide decreased or increased envelope stiffness (respectively), revealing sign epistasis between these components. Based on these data we propose a specific structural model in which the chemical interactions between lipopolysaccharides, β-barrel proteins, and phospholipids coordinately determine cell envelope stiffness, and the ability of the outer membrane to functionally share mechanical loads with the cell wall., Statement of Significance: The outer membrane is the defining cellular structure of Gram-negative bacteria, a group that contains many important pathogens like Escherichia coli , Vibrio cholerae , and Pseudomonas aeruginosa . One role of the outer membrane is to block the uptake of small molecules like antibiotics. However, it is becoming increasingly clear that it also functions as a structural exoskeleton that is critical for the cell's ability to cope with internal and external mechanical forces. Here, we carefully dissect the molecular basis for the load-bearing capacity of the outer membrane by screening a set of mutants with a new cell biophysics assay., Competing Interests: Declaration of interests. The authors declare no competing interest.

In this work, we investigate the question: do code-generating large language models know chemistry? Our results indicate, mostly yes. To evaluate this, we introduce an expandable framework for evaluating chemistry knowledge in these models, through prompting models to solve chemistry problems posed as coding tasks. To do so, we produce a benchmark set of problems, and evaluate these models based on correctness of code by automated testing and evaluation by experts. We find that recent LLMs are able to write correct code across a variety of topics in chemistry and their accuracy can be increased by 30 percentage points via prompt engineering strategies, like putting copyright notices at the top of files. Our dataset and evaluation tools are open source which can be contributed to or built upon by future researchers, and will serve as a community resource for evaluating the performance of new models as they emerge. We also describe some good practices for employing LLMs in chemistry. The general success of these models demonstrates that their impact on chemistry teaching and research is poised to be enormous., Competing Interests: After submission of this manuscript, A. D. W. worked as a paid consultant for OpenAI, the developers of some of the models presented in this work., (This journal is © The Royal Society of Chemistry.)