Hydration of biologically relevant tetramethylammonium cation by neutron scattering and molecular dynamics.

Neutron scattering and molecular dynamics studies were performed on a concentrated aqueous tetramethylammonium (TMA) chloride solution to gain insight into the hydration shell structure of TMA, which is relevant for understanding its behavior in biological contexts of, e.g., properties of phospholipid membrane headgroups or interactions between DNA and histones. Specifically, neutron diffraction with isotopic substitution experiments were performed on TMA and water hydrogens to extract the specific correlation between hydrogens in TMA (HTMA) and hydrogens in water (HW). Classical molecular dynamics simulations were performed to help interpret the experimental neutron scattering data. Comparison of the hydration structure and simulated neutron signals obtained with various force field flavors (e.g. overall charge, charge distribution, polarity of the CH bonds and geometry) allowed us to gain insight into how sensitive the TMA hydration structure is to such changes and how much the neutron signal can capture them. We show that certain aspects of the hydration, such as the correlation of the hydrogen on TMA to hydrogen on water, showed little dependence on the force field. In contrast, other correlations, such as the ion-ion interactions, showed more marked changes. Strikingly, the neutron scattering signal cannot discriminate between different hydration patterns. Finally, ab initio molecular dynamics was used to examine the three-dimensional hydration structure and thus to benchmark force field simulations. Overall, while neutron scattering has been previously successfully used to improve force fields, in the particular case of TMA we show that it has only limited value to fully determine the hydration structure, with other techniques such as ab initio MD being of a significant help.

Structural Effects of Cation Binding to DPPC Monolayers.

Ions at the two sides of the plasma membrane maintain the transmembrane potential, participate in signaling, and affect the properties of the membrane itself. The extracellular leaflet is particularly enriched in phosphatidylcholine lipids and under the influence of Na+, Ca2+, and Cl- ions. In this work, we combined molecular dynamics simulations performed using state-of-the-art models with vibrational sum frequency generation (VSFG) spectroscopy to study the effects of these key ions on the structure of dipalmitoylphosphatidylcholine. We used lipid monolayers as a proxy for membranes, as this approach enabled a direct comparison between simulation and experiment. We find that the effects of Na+ are minor. Ca2+, on the other hand, strongly affects the lipid headgroup conformations and induces a tighter packing of lipids, thus promoting the liquid condensed phase. It does so by binding to both the phosphate and carbonyl oxygens via direct and water-mediated binding modes, the ratios of which depend on the monolayer packing. Clustering analysis performed on simulation data revealed that changes in area per lipid or CaCl2 concentration both affect the headgroup conformations, yet their effects are anticorrelated. Cations at the monolayer surface also attract Cl-, which at large CaCl2 concentrations penetrates deep to the monolayer. This phenomenon coincides with a radical change in the VSFG spectra of the phosphate group, thus indicating the emergence of a new binding mode.

From the outer ear to the nerve: A complete computer model of the peripheral auditory system.

Computer models of the individual components of the peripheral auditory system - the outer, middle, and inner ears and the auditory nerve - have been developed in the past, with varying level of detail, breadth, and faithfulness of the underlying parameters. Building on previous work, we advance the modeling of the ear by presenting a complete, physiologically justified, bottom-up computer model based on up-to-date experimental data that integrates all of these parts together seamlessly. The detailed bottom-up design of the present model allows for the investigation of partial hearing mechanisms and their defects, including genetic, molecular, and microscopic factors. Also, thanks to the completeness of the model, one can study microscopic effects in the context of their implications on hearing as a whole, enabling the correlation with neural recordings and non-invasive psychoacoustic methods. Such a model is instrumental for advancing quantitative understanding of the mechanism of hearing, for investigating various forms of hearing impairment, as well as for devising next generation hearing aids and cochlear implants.

Collective catalysis under spatial constraints: Phase separation and size-scaling effects on mass action kinetics.

Chemical reactions are usually studied under the assumption that both substrates and catalysts are well-mixed (WM) throughout the system. Although this is often applicable to test-tube experimental conditions, it is not realistic in cellular environments, where biomolecules can undergo liquid-liquid phase separation (LLPS) and form condensates, leading to important functional outcomes, including the modulation of catalytic action. Similar processes may also play a role in protocellular systems, like primitive coacervates, or in membrane-assisted prebiotic pathways. Here we explore whether the demixing of catalysts could lead to the formation of microenvironments that influence the kinetics of a linear (multistep) reaction pathway, as compared to a WM system. We implemented a general lattice model to simulate LLPS of a collection of different catalysts and extended it to include diffusion and a sequence of reactions of small substrates. We carried out a quantitative analysis of how the phase separation of the catalysts affects reaction times depending on the affinity between substrates and catalysts, the length of the reaction pathway, the system size, and the degree of homogeneity of the condensate. A key aspect underlying the differences reported between the two scenarios is that the scale invariance observed in the WM system is broken by condensation processes. The main theoretical implications of our results for mean-field chemistry are drawn, extending the mass action kinetics scheme to include substrate initial "hitting times" to reach the catalysts condensate. We finally test this approach by considering open nonlinear conditions, where we successfully predict, through microscopic simulations, that phase separation inhibits chemical oscillatory behavior, providing a possible explanation for the marginal role that this complex dynamic behavior plays in real metabolisms.

Statistical mechanics of biomolecular condensates via cavity methods.

Physical mechanisms of phase separation in living systems play key physiological roles and have recently been the focus of intensive studies. The strongly heterogeneous nature of such phenomena poses difficult modeling challenges that require going beyond mean-field approaches based on postulating a free energy landscape. The pathway we take here is to calculate the partition function starting from microscopic interactions by means of cavity methods, based on a tree approximation for the interaction graph. We illustrate them on the binary case and then apply them successfully to ternary systems, in which simpler one-factor approximations are proved inadequate. We demonstrate the agreement with lattice simulations and contrast our theory with coacervation experiments of associative de-mixing of nucleotides and poly-lysine. Different types of evidence are provided to support cavity methods as ideal tools for modeling biomolecular condensation, giving an optimal balance between the consideration of spatial aspects and fast computational results.

Resolving the Equal Number Density Puzzle: Molecular Picture from Simulations of LiCl(aq) and NaCl(aq).

The change in number densities of aqueous solutions of alkali chlorides should be qualitatively predictable. Typically, as cations get larger, the number density of the solution decreases. However, aqueous solutions of lithium and sodium chloride exhibit at ambient conditions practically identical number densities at equal molalities despite different ionic sizes. Here, we provide an atomistic interpretation of this experimentally observed anomalous behavior using molecular dynamics simulations. The obtained results show that the rigidity of the Li+ first and second solvation shells and the associated compromised hydrogen bonding result in practically equal average water densities in the local hydration regions for Li+ and Na+ despite different sizes of the cations. In addition, in more distant regions from the cations, the water densities of these two solutions also coincide. These findings thus provide an atomistic interpretation for matching number densities of LiCl and NaCl solutions. In contrast, the number density differences between NaCl and KCl solutions as well as between LiCl and KCl solutions behave in a regular fashion with lower number densities of solutions observed for larger cations.

Quantitative linear dichroism imaging of molecular processes in living cells made simple by open software tools.

Fluorescence-detected linear dichroism microscopy allows observing various molecular processes in living cells, as well as obtaining quantitative information on orientation of fluorescent molecules associated with cellular features. Such information can provide insights into protein structure, aid in development of genetically encoded probes, and allow determinations of lipid membrane properties. However, quantitating and interpreting linear dichroism in biological systems has been laborious and unreliable. Here we present a set of open source ImageJ-based software tools that allow fast and easy linear dichroism visualization and quantitation, as well as extraction of quantitative information on molecular orientations, even in living systems. The tools were tested on model synthetic lipid vesicles and applied to a variety of biological systems, including observations of conformational changes during G-protein signaling in living cells, using fluorescent proteins. Our results show that our tools and model systems are applicable to a wide range of molecules and polarization-resolved microscopy techniques, and represent a significant step towards making polarization microscopy a mainstream tool of biological imaging.