Electron transfer through 3D monolayers on Au25 clusters.

The monolayer protecting small gold nanoparticles (monolayer-protected clusters, MPCs) is generally represented as the 3D equivalent of 2D self-assembled monolayers (SAMs) on extended gold surfaces. However, despite the growing relevance of MPCs in important applied areas, such as catalysis and nanomedicine, our knowledge of the structure of 3D SAMs in solution is still extremely limited. We prepared a large series of monodisperse Au25(SCnH2n+1)18 clusters (n=2, 4, 6, 8, 10, 12, 14, 16, 18) and studied how electrons tunnel through these monolayers. Electron transfer results, nicely supported by 1H NMR spectroscopy, IR absorption spectroscopy, and molecular dynamics results, show that there is a critical ligand length marking the transition between short ligands, which form a quite fluid monolayer structure, and longer alkyl chains, which self-organize into bundles. At variance with the truly protecting 2D SAMs, efficient electronic communication of the Au25 core with the outer environment is thus possible even for long alkyl chains. These conclusions provide a different picture of how an ultrasmall gold core talks with the environment through/with its protecting but not-so-shielding monolayer.

Passive Membrane Permeability: Beyond the Standard Solubility-Diffusion Model.

The spontaneous diffusion of solutes through lipid bilayers is still a challenge for theoretical predictions. Since permeation processes remain beyond the capabilities of unbiased molecular dynamics simulations, an alternative strategy is currently adopted to gain insight into their mechanism and time scale. This is based on a monodimensional description of the translocation process only in terms of the position of the solute along the normal to the lipid bilayer, which is formally expressed in the solubility-diffusion model. Actually, a role of orientational and conformational motions has been pointed out, and the use of advanced simulation techniques has been proposed to take into account their effect. Here, we discuss the limitations of the standard solubility-diffusion approach and propose a more general description of membrane translocation as a diffusion process on a free energy surface, which is a function of the translational and rotational degrees of freedom of the molecule. Simple expressions for the permeability coefficient are obtained under suitable conditions. For fast solute reorientation, the classical solubility-diffusion equation is recovered. Under the assumption that well-defined minima can be identified on the free energy landscape, a mechanistic interpretation of the permeability coefficient in terms of all possible permeation paths is given.

Flip-flop of steroids in phospholipid bilayers: effects of the chemical structure on transbilayer diffusion.

The transverse motion of molecules from one leaflet to the other of a lipid bilayer, or flip-flop, represents a putative mechanism for their transmembrane transport and may contribute to the asymmetric distribution of components in biomembranes. However, a clear understanding of this process is still missing. The scarce knowledge derives from the difficulty of experimental determination. Because of its slow rate on the molecular time scale, flip-flop is challenging also for computational techniques. Here, we report a study of the passive transbilayer diffusion of steroids, based on a kinetic model derived from the analysis of their free energy surface, as a function of their position and orientation in the bilayer. An implicit membrane description is used, where the anisotropy and the nonuniformity of the bilayer environment are taken into account in terms of the gradients of density, dielectric permittivity, acyl chain order parameters, and lateral pressure. The flip-flop rates are determined by solving the Master Equation that governs the time evolution of the system, with transition rates between free energy minima evaluated according to the Kramers theory. Considering various steroids (cholesterol, lanosterol, ketosterone, 5-cholestene, 25-hydroxycholesterol, and testosterone), we can discuss how differences in molecular shape and polarity affect the pathway and the rate of flip-flop in a liquid crystalline 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) bilayer, at low steroid concentration. We predict time scales ranging from microseconds to milliseconds, strongly affected by the presence of polar substituents and by their position in the molecular skeleton.

Monomeric fullerenes in lipid membranes: effects of molecular shape and polarity.

We report a combined theoretical and experimental study on the single-molecule interaction of fullerenes with phospholipid membranes. We studied pristine C(60) (1) and two N-substituted fulleropyrrolidines (2 and 3), one of which (3) bore a paramagnetic nitroxide group. Theoretical predictions of fullerene distribution and permeability across lipid bilayers were combined with electron paramagnetic resonance (EPR) experiments in aligned DMPC/DHPC bicelles containing the paramagnetic fulleropyrrolidine 3 or either one of the diamagnetic fullerenes together with spin-labeled lipids. We found that, at low concentrations, fullerenes are present in the bilayer as single molecules. Their preferred location in the membrane is only slightly influenced by the derivatization: all derivatives were confined just below the hydrophilic/hydrophobic interface, because of the key role played by dispersion interactions between the highly polarizable fullerene cage and the hydrocarbon chains, which are especially tight within this region. However, the deviation from spherical shape is sufficient to induce a preferential orientation of 2 and 3 in the membrane. We predict that monomeric fullerenes spontaneously penetrate the bilayer, in agreement with the results of molecular dynamics simulations, but we point out the limits of the currently used permeability model when applied to hydrophobic solutes.

Polarity-sensitive fluorescent probes in lipid bilayers: bridging spectroscopic behavior and microenvironment properties.

We have studied the emission features of the fluorescent polarity-sensitive probes known as Prodan and Laurdan in a liquid-crystalline DPPC bilayer. To this purpose, we have combined high-level quantum mechanical electronic structure calculations with a molecular field theory for the positional-orientational-conformational distribution of the probes, in their ground and excited states, inside of the lipid bilayer, taking into account at both levels the nonuniformity and anisotropy of the environment. Thus, we can interpret the features of the fluorescence spectra of Prodan and Laurdan in relation to the position and orientation of their chromophore in the bilayer. We have found that the environment polarity is not sufficient to explain the large red shifts experimentally observed and that specific effects due to hydrogen bonding must be considered. We show that the orientation of the probe is important in determining the accessibility to water of the H-bond-acceptor group; in the case of Laurdan interesting conformational effects are highlighted.

Solute Partitioning into Lipid Bilayers: An Implicit Model for Nonuniform and Ordered Environment.

We have developed a theoretical and computational methodology to evaluate the coupled orientational-positional distribution of solutes in lipid bilayers. Four different contributions to the solute free energy are considered, which can be traced back to (i) electrostatic and (ii) dispersion interactions between the solute and environment, (iii) work for the formation of a solute-shaped cavity, and (iv) anisotropic interactions with the ordered acyl chains. An atomistic representation of the solute is adopted, which includes the conformational degrees of freedom, whereas an implicit model is used for the water/bilayer environment. The highly nonuniform and anisotropic nature of this is introduced through the profiles of density, dielectric permittivity, lateral pressure, and acyl chain order parameters, which can be derived from experiments or simulations. Effects of chemical composition and physical state of the bilayer can be accounted for by a proper form of these profiles. The methodology which we propose is suitable for the integrated calculation of spectroscopic observables for probes in membranes, for the estimate of partition and permeability coefficients of solutes, and for the implicit modeling of the membrane environment in molecular dynamics and Monte Carlo simulations. Here, the method is presented, and the underlying assumptions are discussed. Cholesterol in the liquid crystalline DPPC bilayer is taken as a case study, to illustrate the capabilities of the proposed approach. Free energy maps, distribution profiles, and orientational properties are shown; they compare well with those obtained from all-atom molecular dynamics simulations, as well as with available experimental data, suggesting that the model used is able to capture the subtle effects of the interplay between intermolecular interaction and nanoscale architecture of the lipid bilayer. The detailed picture provided by our calculations appears suitable to investigate the determinants of the behavior of solutes in lipid membranes, highlighting even nonstraightforward issues, which may have biophysical implications.