Quantifying how step-wise fluorination tunes local solute hydrophobicity, hydration shell thermodynamics and the quantum mechanical contributions of solute-water interactions.

The ability to locally tune solute-water interactions and thus control the hydrophilic/hydrophobic character of a solute is key to control molecular self-assembly and to develop new drugs and biocatalysts; it has been a holy grail in synthetic chemistry and biology. To date, the connection between (i) the hydrophobicity of a functional group; (ii) the local structure and thermodynamics of its hydration shell; and (iii) the relative influence of van der Waals (dispersion) and electrostatic interactions on hydration remains unclear. We investigate this connection using spectroscopic, classical simulation and ab initio methods by following the transition from hydrophile to hydrophobe induced by the step-wise fluorination of methyl groups. Along the transition, we find that water-solute hydrogen bonds are progressively transformed into dangling hydroxy groups. Each structure has a distinct thermodynamic, spectroscopic and quantum-mechanical signature connected to the associated local solute hydrophobicity and correlating with the relative contribution of electrostatics and dispersion to the solute-water interactions.

Hydrophobic but Water-Friendly: Favorable Water-Perfluoromethyl Interactions Promote Hydration Shell Defects.

Although perfluorination is known to enhance hydrophobicity and change protein activity, its influence on hydration-shell structure and thermodynamics remains an open question. Here we address that question by combining experimental Raman multivariate curve resolution spectroscopy with theoretical classical simulations and quantum mechanical calculations. Perfluorination of the terminal methyl group of ethanol is found to enhance the disruption of its hydration-shell hydrogen bond network. Our results reveal that this disruption is not due to the associated volume change but rather to the electrostatic stabilization of the water dangling OH···F interaction. Thus, the hydration shell structure of fluorinated methyl groups results from a delicate balance of solute-water interactions that is intrinsically different from that associated with a methyl group.

Unexpected trends in the hydrophobicity of fluorinated amino acids reflect competing changes in polarity and conformation.

Fluorination can dramatically improve the thermal and proteolytic stability of proteins and their enzymatic activity. Key to the impact of fluorination on protein properties is the hydrophobicity of fluorinated amino acids. We use molecular dynamics simulations, together with a new fixed-charge, atomistic force field, to quantify the changes in hydration free energy, ΔGHyd, for amino acids with alkyl side chains and with 1 to 6 -CH → -CF side chain substitutions. Fluorination changes ΔGHyd by -1.5 to +2 kcal mol-1, but the number of fluorines is a poor predictor of hydrophobicity. Changes in ΔGHyd reflect two main contributions: (i) fluorination alters side chain-water interactions; we identify a crossover point from hydrophilic to hydrophobic fluoromethyl groups which may be used to estimate the hydrophobicity of fluorinated alkyl side-chains; (ii) fluorination alters the number of backbone-water hydrogen bonds via changes in the relative side chain-backbone conformation. Our results offer a road map to mechanistically understand how fluorination alters hydrophobicity of (bio)polymers.

Diphenylhexatriene membrane probes DPH and TMA-DPH: A comparative molecular dynamics simulation study.

Fluorescence spectroscopy and microscopy have been utilized as tools in membrane biophysics for decades now. Because phospholipids are non-fluorescent, the use of extrinsic membrane probes in this context is commonplace. Among the latter, 1,6-diphenylhexatriene (DPH) and its trimethylammonium derivative (TMA-DPH) have been extensively used. It is widely believed that, owing to its additional charged group, TMA-DPH is anchored at the lipid/water interface and reports on a bilayer region that is distinct from that of the hydrophobic DPH. In this study, we employ atomistic MD simulations to characterize the behavior of DPH and TMA-DPH in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and POPC/cholesterol (4:1) bilayers. We show that although the dynamics of TMA-DPH in these membranes is noticeably more hindered than that of DPH, the location of the average fluorophore of TMA-DPH is only ~3-4Å more shallow than that of DPH. The hindrance observed in the translational and rotational motions of TMA-DPH compared to DPH is mainly not due to significant differences in depth, but to the favorable electrostatic interactions of the former with electronegative lipid atoms instead. By revealing detailed insights on the behavior of these two probes, our results are useful both in the interpretation of past work and in the planning of future experiments using them as membrane reporters.

Influence of the sterol aliphatic side chain on membrane properties: a molecular dynamics study.

Following a recent experimental investigation of the effect of the length of the alkyl side chain in a series of cholesterol analogues (Angew. Chem., Int. Ed., 2013, 52, 12848-12851), we report here an atomistic molecular dynamics characterization of the behaviour of methyl-branched side chain sterols (iso series) in POPC bilayers. The studied sterols included androstenol (i-C0-sterol) and cholesterol (i-C8-sterol), as well as four other derivatives (i-C5, i-C10, i-C12 and i-C14-sterol). For each sterol, both subtle local effects and more substantial differential alterations of membrane properties along the iso series were investigated. The location and orientation of the tetracyclic ring system is almost identical in all compounds. Among all the studied sterols, cholesterol is the sterol that presents the best matching with the hydrophobic length of POPC acyl chains, whereas longer-chained sterols interdigitate into the opposing membrane leaflet. In accordance with the experimental observations, a maximal ordering effect is observed for intermediate sterol chain length (i-C5, cholesterol, i-C10). Only for these sterols a preferential interaction with the saturated sn-1 chain of POPC (compared to the unsaturated sn-2 chain) was observed, but not for either shorter or longer-chained derivatives. This work highlights the importance of the sterol alkyl chain in the modulation of membrane properties and lateral organization in biological membranes.

Influence of Methylene Fluorination and Chain Length on the Hydration Shell Structure and Thermodynamics of Linear Diols.

The interplay between the local hydration shell structure, the length of hydrophobic solutes, and their identity (perfluorinated or not) remains poorly understood. We address this issue by combining Raman-multivariate curve resolution (Raman-MCR) spectroscopy, simulation, and quantum-mechanical calculations to quantify the thermodynamics and the first principle interactions behind the formation of defects in the hydration shell of alkyl-diol and perfluoroalkyl-diol chains. The hydration shell of the fluorinated diols contains substantially more defects than that of the nonfluorinated diols; these defects are water hydroxy groups that do not donate hydrogen bonds and which either point to the solute (radial-dangling OH) or not (nonradial-dangling OH). The number of radial-dangling OH defects per carbon decreases for longer chains and toward the interior of the fluorinated diols, mainly due to less favorable electrostatics and exchange interactions; nonradial-dangling OH defects per carbon increase with chain length. In contrast, the hydration shell of the nonfluorinated diols only contains radial-dangling defects, which become more abundant toward the center of the chain and for larger chains, predominantly because of more favorable dispersion interactions. These results have implications for how the folding of macromolecules, ligand binding to biomacromolecules, and chemical reactions at water-oil interfaces could be modified through the introduction of fluorinated groups or solvents.

NBD-labeled cholesterol analogues in phospholipid bilayers: insights from molecular dynamics.

Nitrobenzoxadiazole (NBD)-labeled sterols are commonly used as fluorescent cholesterol analogues in membrane biophysics. However, some experimental reports have questioned their ability to emulate the behavior of cholesterol in phospholipid bilayers. For the purpose of a detailed clarification of this matter, atomistic molecular dynamics simulations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers, containing either cholesterol or one of two fluorescent cholesterol analogues, 22-NBD-cholesterol or 25-NBD-cholesterol, were carried out. It is found that these sterol probes tend to adopt conformations in which their tail-labeled fluorophore is oriented toward the lipid/water interface, with a location similar to that observed in molecular dynamics simulations of other NBD probes. This implies that in these molecules the long sterol axis is no longer aligned with the membrane normal, and preferentially adopts orientations approximately parallel to the bilayer plane. In turn, these stretched conformations, together with NBD-POPC atomic interactions, lead to slowed-down lateral diffusion of both fluorescent sterols, compared to cholesterol. From computation of the deuterium order parameter and acyl chain tilts of POPC chains for varying POPC-sterol distance, it is observed that the local ordering effect of sterol is altered in both fluorescent derivatives. In agreement with reported experimental data, both fluorescent sterols are able to increase the order of POPC at 20 mol % concentration (as some molecules adopt an upright conformation, possibly related to formation of transbilayer aggregates), albeit to a smaller extent to that of cholesterol. Altogether, this study indicates that both 22- and 25-NBD-cholesterol are unable to mimic the most important features of cholesterol's behavior in lipid bilayers.

Behavior of fluorescent cholesterol analogues dehydroergosterol and cholestatrienol in lipid bilayers: a molecular dynamics study.

Molecular dynamics simulations of bilayer systems consisting of varying proportions of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), cholesterol (Chol), and intrinsically fluorescent Chol analogues dehydroergosterol (DHE) or cholestatrienol (CTL) were carried out to study in detail the extent to which these fluorescent probes mimic Chol's behavior (location, orientation, dynamics) in membranes as well as their effect on host bilayer structure and dynamics (namely their ability to induce membrane ordering in comparison with Chol). Control properties of POPC and POPC/Chol bilayers agree well with published experimental and simulation work. Both probes and Chol share similar structural and dynamical properties within the bilayers. Additionally, the fluorescent sterols induce membrane ordering to a similar (slightly lower) extent to that of Chol. These findings combined demonstrate that the two studied fluorescent sterols are adequate analogues of Chol, and may be used with advantage over side-chain labeled sterols. The small structural differences between the three studied sterols are responsible for the slight variations in the calculated properties, with CTL presenting a more similar behavior to Chol (correlating with its larger structural similarity to Chol) compared to DHE.