Interfacial water on collagen nanoribbons by 3D AFM.

Collagen is the most abundant structural protein in mammals. Type I collagen in its fibril form has a characteristic pattern structure that alternates two regions called gap and overlap. The structure and properties of collagens are highly dependent on the water and mineral content of the environment. Here, we apply 3D AFM to characterize at angstrom-scale resolution the interfacial water structure of collagen nanoribbons. For a neutral tip, the interfacial water structure is characterized by the oscillation of the water particle density distribution with a value of 0.3 nm (hydration layers). The interfacial structure does not depend on the collagen region. For a negatively charged tip, the interfacial structure might depend on the collagen region. Hydration layers are observed in overlap regions, while in gap regions, the interfacial solvent structure is dominated by electrostatic interactions. These interactions generate interlayer distances of 0.2 nm.

Layering and capillary waves in the structure factor of liquid surfaces.

Within the extended Capillary Wave Theory (ECWT), to extract the bending modulus of a liquid surface from the total structure factor of the interfacial region requires to separate the capillary waves (CW) signal from a non-CW background. Some years ago, Höfling and Dietrich (HD), working in the strict grazing incidence limit qz = 0, proposed a background that combines the liquid and vapor bulk structure factors in the amounts set by Gibbs's plane. We contrast that proposal with Molecular Dynamics (MD) simulations of the Lennard-Jones model analyzed with the Intrinsic Sampling Method (ISM). The study is extended to qz ≠ 0, to test the stronger consistency requirements of the ECWT and the experimental conditions; it shows a good MD-ECWT matching although we need some fine tuning over HD proposal. Then, the agreement with the ISM result for the surface bending modulus is good and that provides an interpretation, in terms of the molecular layering at the liquid edge, for the fluctuating surface represented by the CW signal in the surface structure factor, both for MD simulations and surface diffraction experiments.

A coarse-grained approach to model the dynamics of the actomyosin cortex.

BACKGROUND: The dynamics of the actomyosin machinery is at the core of many important biological processes. Several relevant cellular responses such as the rhythmic compression of the cell cortex are governed, at a mesoscopic level, by the nonlinear interaction between actin monomers, actin crosslinkers, and myosin motors. Coarse-grained models are an optimal tool to study actomyosin systems, since they can include processes that occur at long time and space scales, while maintaining the most relevant features of the molecular interactions. RESULTS: Here, we present a coarse-grained model of a two-dimensional actomyosin cortex, adjacent to a three-dimensional cytoplasm. Our simplified model incorporates only well-characterized interactions between actin monomers, actin crosslinkers and myosin, and it is able to reproduce many of the most important aspects of actin filament and actomyosin network formation, such as dynamics of polymerization and depolymerization, treadmilling, network formation, and the autonomous oscillatory dynamics of actomyosin. CONCLUSIONS: We believe that the present model can be used to study the in vivo response of actomyosin networks to changes in key parameters of the system, such as alterations in the attachment of actin filaments to the cell cortex.

Bending Modulus of Lipid Membranes from Density Correlation Functions.

The bending modulus κ quantifies the elasticity of biological membranes in terms of the free energy cost of increasing the membrane corrugation. Molecular dynamics (MD) simulations provide a powerful approach to quantify κ by analyzing the thermal fluctuations of the lipid bilayer. However, existing methods require the identification and filtering of non-mesoscopic fluctuation modes. State of the art methods rely on identifying a smooth surface to describe the membrane shape. These methods introduce uncertainties in calculating κ since they rely on different criteria to select the relevant fluctuation modes. Here, we present a method to compute κ using molecular simulations. Our approach circumvents the need to define a mesoscopic surface or an orientation field for the lipid tails explicitly. The bending and tilt moduli can be extracted from the analysis of the density correlation function (DCF). The method introduced here builds on the Bedeaux and Weeks (BW) theory for the DCF of fluctuating interfaces and on the coupled undulatory (CU) mode introduced by us in previous work. We test the BW-DCF method by computing the elastic properties of lipid membranes with different system sizes (from 500 to 6000 lipid molecules) and using coarse-grained (for POPC and DPPC lipids) and fully atomistic models (for DPPC). Further, we quantify the impact of cholesterol on the bending modulus of DPPC bilayers. We compare our results with bending moduli obtained with X-ray diffraction data and different computer simulation methods.

Intrinsic structure perspective for MIPS interfaces in two-dimensional systems of active Brownian particles.

Suspensions of active Brownian particles (ABPs) undergo motility-induced phase separation (MIPS) over a wide range of mean density and activity strength, which implies the spontaneous aggregation of particles due to the persistence of their direction of motion, even in the absence of an explicit attraction. Both similarities and qualitative differences have been obtained when the MIPS is analysed in the same terms as a liquid-gas phase coexistence in an equilibrium attractive system. Negative values of the mechanical surface tension have been reported, from the total forces across the interface, while stable fluctuations of the interfacial line could be interpreted as a positive capillary surface tension; in equilibrium liquid surfaces, these two magnitudes are equal. We present here the analysis of 2D-ABP interfaces in terms of the intrinsic density and force profiles, calculated with the particle distance to the instantaneous interfacial line. Our results provide new insight into the origin of MIPS from the local rectification of the random active force on the particles near the interface. As has been reported, this effect acts as an external potential that produces a pressure gradient across the interface, such that the mechanical surface tension of the MIPS cannot be described as that of equilibrium coexisting phases; however, our analysis shows that most of that effect comes from the tightly caged particles at the dense (inner) side of the MIPS interface, rather than from the free moving particles at the outer side that collide with the dense cluster. Moreover, a clear correlation appears between the decay of the hexatic order parameter at the dense slab and the end of the MIPS as the strength of the active force is lowered. We show that, using the strong active forces required for MIPS, the interfacial structure and properties are very similar for ABPs with purely repulsive interactions (the Weeks-Chandler-Andersen-Lennard-Jones (WCA-LJ) model truncated at its minimum) and when the interaction includes a range of the LJ attractive force.

Tip Charge Dependence of Three-Dimensional AFM Mapping of Concentrated Ionic Solutions.

A molecular scale understanding of the organization and structure of a liquid near a solid surface is currently a major challenge in surface science. It has implications across different fields from electrochemistry and energy storage to molecular biology. Three-dimensional AFM generates atomically resolved maps of solid-liquid interfaces. The imaging mechanism behind those maps is under debate, in particular, for concentrated ionic solutions. Theory predicts that the observed contrast should depend on the tip's charged state. Here, by using neutrally, negatively, and positively charged tips, we demonstrate that the 3D maps depend on the tip's polarization. A neutral tip will explore the total particle density distribution (water and ions) while a charged tip will reveal the charge density distribution. The experimental data reproduce the key findings of the theory.

Density functional analysis of atomic force microscopy in a dense fluid.

We present a density functional (DF) analysis for the entropic force in Atomic Force Microscopy (AFM) across the layers of a dense fluid. Previous theoretical analysis, based on the ideal gas entropy, was apparently supported by the similarity in the oscillatory decay for the force and density profile. We point out that such similarity is a generic DF result, which carries no information on the interface, since the decaying mode is characteristic of the bulk fluid correlation. The truly interfacial information, from the layering measured by AFM, comes in its amplitude and not in the decay mode. With our rigorous study of a simple hard sphere model, we find semiempirical clues to disentangle the role of the tip radius and to relate the amplitude of the molecular layering to the oscillatory force on the AFM tip.

A MOF@COF Composite with Enhanced Uptake through Interfacial Pore Generation.

Herein, we describe a new class of porous composites comprising metal-organic framework (MOF) crystals confined in single spherical matrices made of packed covalent-organic framework (COF) nanocrystals. These MOF@COF composites are synthesized through a two-step method of spray-drying and subsequent amorphous (imine-based polymer)-to-crystalline (imine-based COF) transformation. This transformation around the MOF crystals generates micro- and mesopores at the MOF/COF interface that provide far superior porosity compared to that of the constituent MOF and COF components added together. We report that water sorption in these new pores occurs within the same pressure window as in the COF pores. Our new MOF@COF composites, with their additional pores at the MOF/COF interface, should have implications for the development of new composites.

Bacterial cell division: modeling FtsZ assembly and force generation from single filament experimental data.

The bacterial cytoskeletal protein FtsZ binds and hydrolyzes GTP, self-aggregates into dynamic filaments and guides the assembly of the septal ring on the inner side of the membrane at midcell. This ring constricts the cell during division and is present in most bacteria. Despite exhaustive studies undertaken in the last 25 years after its discovery, we do not yet know the mechanism by which this GTP-dependent self-aggregating protein exerts force on the underlying membrane. This paper reviews recent experiments and theoretical models proposed to explain FtsZ filament dynamic assembly and force generation. It highlights how recent observations of single filaments on reconstituted model systems and computational modeling are contributing to develop new multiscale models that stress the importance of previously overlooked elements as monomer internal flexibility, filament twist and flexible anchoring to the cell membrane. These elements contribute to understand the rich behavior of these GTP consuming dynamic filaments on surfaces. The aim of this review is 2-fold: (1) to summarize recent multiscale models and their implications to understand the molecular mechanism of FtsZ assembly and force generation and (2) to update theoreticians with recent experimental results.

Density correlation in liquid surfaces: Bedeaux-Weeks high order terms and non capillary wave background.

We present Molecular Dynamics (MD) simulations of liquid-vapor surfaces, and their Intrinsic Sampling Method analysis, to get a quantitative test for the theoretical prediction of the capillary wave (CW) effects on density correlation done by Bedeaux and Weeks (BW) in 1985. The results are contrasted with Wertheim's proposal which is the first term in BW series and are complemented with a (formally defined and computational accessible) proposal for the background of non-CW fluctuations. Our conclusion is that BW theory is both accurate and needed since it may differ significantly from Wertheim's proposal. We discuss the implications for the analysis of experimental X-ray surface diffraction data and MD simulations.

Capillary waves as eigenmodes of the density correlation at liquid surfaces.

We analyze the density correlations in a liquid-vapor surface to establish a quantitative connection between the Density Functional (DF) formalism, Molecular Dynamic (MD) simulations, and the Capillary Wave (CW) theory. Instead of the integrated structure factor, we identify the CW fluctuations as eigenmodes of the correlation function. The square-gradient DF approximation appears as fully consistent with the use of the thermodynamic surface tension to describe the surface fluctuations for any wavevector because it misses the upper cutoff in the surface Hamiltonian from the merging of the CW mode with the non-CW band. This mesoscopic cutoff may be accurately predicted from the main peak in the structure factor of the bulk liquid. We explore the difference between the full density-density correlation mode and the bare CW that represents the correlation between the corrugation of the intrinsic surface and the density at the interfacial region. The non-local decay of the CW effects, predicted from DF analysis and observed in MD simulations with the intrinsic sampling method, is found to characterize the bare CW fluctuations, which also require a wavevector-dependent surface tension.

Deconstructing Temperature Gradients across Fluid Interfaces: The Structural Origin of the Thermal Resistance of Liquid-Vapor Interfaces.

The interfacial thermal resistance determines condensation-evaporation processes and thermal transport across material-fluid interfaces. Despite its importance in transport processes, the interfacial structure responsible for the thermal resistance is still unknown. By combining nonequilibrium molecular dynamics simulations and interfacial analyses that remove the interfacial thermal fluctuations we show that the thermal resistance of liquid-vapor interfaces is connected to a low density fluid layer that is adsorbed at the liquid surface. This thermal resistance layer (TRL) defines the boundary where the thermal transport mechanism changes from that of gases (ballistic) to that characteristic of dense liquids, dominated by frequent particle collisions involving very short mean free paths. We show that the thermal conductance is proportional to the number of atoms adsorbed in the TRL, and hence we explain the structural origin of the thermal resistance in liquid-vapor interfaces.

Atomically resolved three-dimensional structures of electrolyte aqueous solutions near a solid surface.

Interfacial liquid layers play a central role in a variety of phenomena ranging from friction to molecular recognition. Liquids near a solid surface form an interfacial layer where the molecular structure is different from that of the bulk. Here we report atomic resolution three-dimensional images of electrolyte solutions near a mica surface that demonstrate the existence of three types of interfacial structures. At low concentrations (0.01-1 M), cations are adsorbed onto the mica. The cation layer is topped by a few hydration layers. At higher concentrations, the interfacial layer extends several nanometres into the liquid. It involves the alternation of cation and anion planes. Fluid Density Functional calculations show that water molecules are a critical factor for stabilizing the structure of the interfacial layer. The interfacial layer stabilizes a crystal-like structure compatible with liquid-like ion and solvent mobilities. At saturation, some ions precipitate and small crystals are formed on the mica.

Collective effects of torsion in FtsZ filaments.

Recent evidence points to the presence of torsion in FtsZ bonds. In addition, experiments with FtsZ mutants on surfaces resulted in new aggregates that cannot be explained by older models for FtsZ dynamics. We use an interaction model for FtsZ derived from molecular dynamics simulations and expand a fine-grained lattice model used to describe FtsZ aggregates on a surface. This new model includes different anchoring angles for the monomers and allows bond twist, two ingredients that oppose each other resulting in a more dynamic and interesting system. We study the role and importance of these conflicting elements and how the aggregates are characterized by the different interaction parameters.

Capillary wave Hamiltonian for the Landau-Ginzburg-Wilson density functional.

We study the link between the density functional (DF) formalism and the capillary wave theory (CWT) for liquid surfaces, focused on the Landau-Ginzburg-Wilson (LGW) model, or square gradient DF expansion, with a symmetric double parabola free energy, which has been extensively used in theoretical studies of this problem. We show the equivalence between the non-local DF results of Parry and coworkers and the direct evaluation of the mean square fluctuations of the intrinsic surface, as is done in the intrinsic sampling method for computer simulations. The definition of effective wave-vector dependent surface tensions is reviewed and we obtain new proposals for the LGW model. The surface weight proposed by Blokhuis and the surface mode analysis proposed by Stecki provide consistent and optimal effective definitions for the extended CWT Hamiltonian associated to the DF model. A non-local, or coarse-grained, definition of the intrinsic surface provides the missing element to get the mesoscopic surface Hamiltonian from the molecular DF description, as had been proposed a long time ago by Dietrich and coworkers.

Capillary waves and the decay of density correlations at liquid surfaces.

Wertheim predicted strong density-density correlations at free liquid surfaces, produced by capillary wave fluctuations of the interface [M. S. Wertheim, J. Chem. Phys. 65, 2377 (1976)JCPSA60021-960610.1063/1.433352]. That prediction has been used to search for a link between capillary wave (CW) theory and density functional (DF) formalism for classical fluids. In particular, Parry et al. have recently analyzed the decaying tails of these CW effects moving away from the interface as a clue for the extended CW theory [A. O. Parry et al., J. Phys.: Condens. Matter 28, 244013 (2016)JCOMEL0953-898410.1088/0953-8984/28/24/244013], beyond the strict long-wavelength limit studied by Wertheim. Some apparently fundamental inconsistencies between the CW and the DF theoretical views of the fluid interfaces arose from the asymptotic analysis of the CW signal. In this paper we revisit the problem of the CW asymptotic decay with a separation of local non-CW surface correlation effects from those that are a truly nonlocal propagation of the CW fluctuations from the surface towards the liquid bulk.

A computer simulation approach to quantify the true area and true area compressibility modulus of biological membranes.

We present a new computational approach to quantify the area per lipid and the area compressibility modulus of biological membranes. Our method relies on the analysis of the membrane fluctuations using our recently introduced coupled undulatory (CU) mode [Tarazona et al., J. Chem. Phys. 139, 094902 (2013)], which provides excellent estimates of the bending modulus of model membranes. Unlike the projected area, widely used in computer simulations of membranes, the CU area is thermodynamically consistent. This new area definition makes it possible to accurately estimate the area of the undulating bilayer, and the area per lipid, by excluding any contributions related to the phospholipid protrusions. We find that the area per phospholipid and the area compressibility modulus features a negligible dependence with system size, making possible their computation using truly small bilayers, involving a few hundred lipids. The area compressibility modulus obtained from the analysis of the CU area fluctuations is fully consistent with the Hooke's law route. Unlike existing methods, our approach relies on a single simulation, and no a priori knowledge of the bending modulus is required. We illustrate our method by analyzing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers using the coarse grained MARTINI force-field. The area per lipid and area compressibility modulus obtained with our method and the MARTINI forcefield are consistent with previous studies of these bilayers.

Mesoscopic Hamiltonian for the fluctuations of adsorbed Lennard-Jones liquid films.

We use Monte Carlo simulations of a Lennard-Jones fluid adsorbed on a short-range planar wall substrate to study the fluctuations in the thickness of the wetting layer, and we get a quantitative and consistent characterization of their mesoscopic Hamiltonian, H[ξ]. We have observed important finite-size effects, which were hampering the analysis of previous results obtained with smaller systems. The results presented here support an appealing simple functional form for H[ξ], close but not exactly equal to the theoretical nonlocal proposal made on the basis a generic density-functional analysis by Parry and coworkers. We have analyzed systems under different wetting conditions, as a proof of principle for a method that provides a quantitative bridge between the molecular interactions and the phenomenology of wetting films at mesoscopic scales.

Modeling the interplay between protein and lipid aggregation in supported membranes.

We present a theoretical model that deals with the complex interplay between lipid segregation and the self-aggregation of lipid-attached proteins. The model, in contrast to previous ones that consider proteins only as passive elements affecting the lipid distribution, describes the system including three terms: the dynamic interactions between protein monomers, the interactions between lipid components, and a mixed term considering protein-lipid interactions. It is used to explain experimental results performed on a well-defined system in which a self-aggregating soluble bacterial cytoskeletal protein polymerizes on a lipid bilayer containing two lipid components. All the elements considered in a previously described protein model, including torsion of the monomers within the filament, are needed to account for the observed filament shapes. The model also points out that lipid segregation can affect the length and curvature of the filaments and that the dynamic behavior of the lipids and proteins can have different time scales, giving rise to memory effects. This simple model that considers a dynamic protein assembly on a fluid and active lipid surface can be easily extended to other biologically relevant situations in which the interplay between protein and lipid aggregation is needed to fully describe the system.

Effect of dispersion forces on the capillary-wave fluctuations of liquid surfaces.

We present molecular dynamics evidence for the nonanalytic effects of the long-range dispersion forces on the capillary waves fluctuations of a Lennard-Jones liquid surface. The results of the intrinsic sampling method, for the analysis of the instantaneous interfacial shape, are obtained in large systems for several cut-off distances of the potential tail, and they show good agreement with the theoretical prediction by Napiórkowski and Dietrich, based on a density functional analysis. The enhancement of the capillary waves is quantified to be within 1% for a simple liquid near its triple point.