The popular MARTINI coarse-grained model is used as a test case to analyze the adherence of top-down coarse-grained molecular dynamics models (i.e., models primarily parametrized to match experimental results) to the known features of statistical mechanics for the underlying all-atom representations. Specifically, the temperature dependence of various pair distribution functions, and hence their underlying potentials of mean force via the reversible work theorem, are compared between MARTINI 2.0, Dry MARTINI, and all-atom simulations mapped onto equivalent coarse-grained sites for certain lipid bilayers. It is found that the MARTINI models do not completely capture the lipid structure seen in atomistic simulations as projected onto the coarse-grained mappings and that issues of accuracy and temperature transferability arise due to an incorrect enthalpy–entropy decomposition of these potentials of mean force. The potential of mean force for the association of two amphipathic helices in a lipid bilayer is also calculated, and especially at shorter ranges, the MARTINI and all-atom projection results differ substantially. The former is much less repulsive and hence will lead to a higher probability of MARTINI helix association in the MARTINI bilayer than occurs in the actual all-atom case. Additionally, the bilayer height fluctuation spectra are calculated for the MARTINI model, and compared to the all-atom results, it is found that the magnitude of thermally averaged amplitudes at intermediate length scales are quite different, pointing to a number of possible consequences for realistic modeling of membrane processes. Taken as a whole, the results presented here show disagreement in the enthalpic and entropic driving forces driving lateral structure in lipid bilayers as well as quantitative differences in association of embedded amphipathic helices, which can help direct future efforts to parametrize CG models with better agreement to the all-atom systems they aspire to represent.

Protein-mediated membrane remodeling is a ubiquitous and critical process for proper cellular function. Inverse Bin/Amphiphysin/Rvs (I-BAR) domains drive local membrane deformation as a precursor to large-scale membrane remodeling. We employ a multiscale approach to provide the molecular mechanism of unusual I-BAR domain-driven membrane remodeling at a low protein surface concentration with near-atomistic detail. We generate a bottom-up coarse-grained model that demonstrates similar membrane-bound I-BAR domain aggregation behavior as our recent Mesoscopic Membrane with Explicit Proteins model. Together, these models bridge several length scales and reveal an aggregation behavior of I-BAR domains. We find that at low surface coverage (i.e., low bound protein density), I-BAR domains form transient, tip-to-tip strings on periodic flat membrane sheets. Inside of lipid bilayer tubules, we find linear aggregates parallel to the axis of the tubule. Finally, we find that I-BAR domains form tip-to-tip aggregates around the edges of membrane domes. These results are supported by in vitro experiments showing low curvature bulges surrounded by I-BAR domains on giant unilamellar vesicles. Overall, our models reveal new I-BAR domain aggregation behavior in membrane tubules and on the surface of vesicles at low surface concentration that add insight into how I-BAR domain proteins may contribute to certain aspects of membrane remodeling in cells.

Filopodia are actin-rich membrane protrusions essential for cell morphogenesis, motility, and cancer invasion. How cells control filopodium initiation on the plasma membrane remains elusive. We performed experiments in cellulo, in vitro, and in silico to unravel the mechanism of filopodium initiation driven by the membrane curvature sensor IRSp53 (insulin receptor substrate protein of 53 kDa). We showed that full-length IRSp53 self-assembles into clusters on membranes depending on PIP 2 . Using well-controlled in vitro reconstitution systems, we demonstrated that IRSp53 clusters recruit the actin polymerase VASP (vasodilator-stimulated phosphoprotein) to assemble actin filaments locally on membranes, leading to the generation of actin-filled membrane protrusions reminiscent of filopodia. By pulling membrane nanotubes from live cells, we observed that IRSp53 can only be enriched and trigger actin assembly in nanotubes at highly dynamic membrane regions. Our work supports a regulation mechanism of IRSp53 in its attributes of curvature sensation and partner recruitment to ensure a precise spatial-temporal control of filopodium initiation.

We investigated the conformational free energy landscape of glucose solvated in water and in the ionic liquids (ILs) 1-butyl-3-methylimidazolium chloride ([Bmim][Cl]) and 1-butyl-3-methylimidazoulim boron tetrafluoride ([Bmim][BF4]). To quantify equilibrium thermodynamic solvent effects, molecular dynamics simulations in conjunction with enhanced sampling based on the metadynamics framework were used. The results show that the solvent choice induces significant differences in the equilibrium ring structures, which may help further resolve the molecular mechanism governing IL-mediated cellulose dissolution.

Filopodia are actin-rich membrane protrusions essential for cell morphogenesis, motility, and cancer invasion. How cells control filopodia initiation on the plasma membrane remains elusive. We performed experiments in cellulo, in vitro and in silico to unravel the mechanism of filopodia initiation driven by the membrane curvature sensor IRSp53. We showed that full-length IRSp53 self-assembles into clusters on membranes depending on PIP2. Using well-controlled in vitro reconstitution systems, we demonstrated that IRSp53 clusters recruit the actin polymerase VASP to assemble actin filaments locally on membranes, leading to the generation of actin-filled membrane protrusions reminiscent of filopodia. By pulling membrane nanotubes from live cells, we observed that IRSp53 can only be enriched and trigger actin assembly in nanotubes at highly dynamic membrane regions. Our work supports a regulation mechanism of IRSp53 in its attributes of curvature sensation and partner recruitment to ensure a precise spatial-temporal control of filopodia initiation.