Hoechst 33342 (H33342) is a fluorescent probe that is commonly used to stain the DNA of living cells. To do so, it needs to interact with and permeate through cell membranes, despite its high overall charge at physiological pH values. In this work, we address the effect of pH in the association of H33342 with lipid bilayers using a combined experimental and computational approach. The partition of H33342 to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid membranes was experimentally quantified using fluorescence spectroscopy and isothermal titration calorimetry (ITC) measurements. Quantum chemical calculations were performed to select the most stable isomer of H33342 for the overall charges 0, +1, and +2, expected to predominate across the 5 < pH < 10 range. The interaction of these isomers with POPC bilayers was then studied by both unrestrained and umbrella sampling molecular dynamics (MD) simulations. Both experimental results and computational free energy profiles indicate that the partition coefficient of H33342 displays a small variation over a wide pH range, not exceeding one order of magnitude. The enthalpy variation upon partition to the membrane suggests efficient hydrogen bonding between the probe and the lipid, namely, for the protonated +2 form, which was confirmed in the MD simulation studies. The relatively high lipophilicity obtained for the charged species contrasts with the decrease in their general hydrophobicity as estimated from octanol/water partition. This highlights the distinction between lipophilicity and hydrophobicity, as well as the importance of considering the association with lipid bilayers when predicting the affinity for biomembranes.

Rifampicin is a broad-spectrum antibiotic, active against several bacterial infections such as tuberculosis. It is a relatively large and structurally complex molecule, including numerous polar groups. Although violating several of Lipinski’s rules for efficient intestinal absorption, rifampicin shows good oral bioavailability, permeating through cell membranes in the absorption pathway and those of the target organisms. Some hypotheses have been proposed for its efficient membrane permeation, but the details are mostly unknown. In this work, the interaction of rifampicin with POPC lipid bilayers is studied using experimental biophysics methodologies and atomistic molecular dynamics simulations considering the two most prevalent ionic species at physiological pH, the anionic and the zwitterionic forms. The results show that both ionization forms of rifampicin establish favorable interactions with the membrane lipids, in agreement with the relatively high partition coefficient obtained experimentally. The results from MD simulations and isothermal titration calorimetry using different pH buffers show that the piperazine group inserts deeply in the membrane and is accompanied by a stabilization of its neutral form. The bulky nature of rifampicin and its deep insertion in the membrane lead to a strong perturbation in the lipids local order, decreasing the membrane barrier properties as evaluated from the rate of carboxyfluorescein leaching. Altogether, the comparison between the experimental and MD simulations results provides important insight regarding the rifampicin molecular features responsible for its relatively fast membrane permeation. The lipid POPC used in this study was selected as a simple membrane with relevance for different organisms across all kingdoms. Further studies using more complex lipid compositions will provide details on eventual specificities for rifampicin interaction with the membranes of distinct organisms.

Hoechst 33342 (H33342) is a fluorescent probe that is commonly used to stain DNA of living cells. To do so, it needs to interact with and permeate through cell membranes. The pH of biological media may vary quite significantly and this influences H33342 ionization state, and may change its rate of permeation through biomembranes. In this work, we address the effect of pH in the association of H33342 with lipid bilayers using a combined experimental and computational approach. The partition of H33342 to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid membranes was experimentally quantified using fluorescence spectroscopy and isothermal titration calorimetry (ITC) measurements. Quantum chemical calculations were performed to select the most stable isomer of H33342 for the overall charges 0, +1, and +2, expected to predominate across the 5 < pH < 10 range. The interaction of these isomers with POPC bilayers was then studied by both unrestrained and umbrella sampling molecular dynamics simulations. While probe insertion occurs for all ionization states, with similar overall equilibrium location, orientation was found to depend on the probe charge, with the neutral form lying mostly parallel and the doubly charged displaying a significant tilt relative to the membrane plane. Both experimental results and computational free energy profiles indicate that the partition coefficient of H33342 displays a small variation over a wide pH range, not exceeding one order of magnitude. The enthalpy of interaction is negative and is larger at lower pH values. This suggests efficient hydrogen bonding between the probe and the lipid, namely for the protonated +2 form. The relatively high lipophilicity obtained for the charged species contrasts with the decrease in their general hydrophobicity as estimated from water and octanol partition. This highlights the distinction between lipophilicity and hydrophobicity, and the importance of considering association with lipid bilayers when predicting the affinity for biomembranes.

Predicting the rate at which substances permeate membrane barriers in vivo is crucial for drug development. Permeability coefficients obtained from in vitro studies are valuable for this goal. These are normally determined by following the dynamics of solute equilibration between two membrane-separated compartments. However, the correct calculation of permeability coefficients from such data is not always straightforward. To address these problems, here we develop a kinetic model for solute permeation through lipid membrane barriers that includes the two membrane leaflets as compartments in a four-compartment model. Accounting for solute association with the membrane allows assessing various methods in a wide variety of conditions. The results showed that the often-used expression P_app = β × r/3 is inapplicable to very large or very small vesicles, to moderately or highly lipophilic solutes, or when the development of a significant pH gradient opposes the solute’s flux. We establish useful relationships that overcome these limitations and allow predicting permeability in compartmentalised in vitro or in vivo systems with specific properties. Finally, from the parameters for the interaction of the solute with the membrane barrier, we defined an intrinsic permeability coefficient that facilitates quantitative comparisons between solutes.

Rifampicin is a broad-spectrum antibiotic, active against several bacterial infections such as tuberculosis. It is a relatively large and structurally complex molecule, including numerous polar groups. Although violating several of Lipinski’s rules for efficient intestinal absorption, rifampicin permeates through cell membranes of its target organisms and shows good oral bioavailability. Some hypotheses have been proposed for its efficient membrane permeation, but the details are mostly unknown. In this work, the interaction of rifampicin with POPC lipid bilayers is studied using experimental biophysics methodologies and atomistic molecular dynamics simulations considering the two most prevalent ionic species at physiological pH, the anionic and the zwitterionic forms. The results show that both ionization forms of rifampicin establish favorable interactions with the membrane lipids, in agreement with the relatively high partition coefficient obtained experimentally. Results from MD simulations and Isothermal Titration Calorimetry using different pH buffers show that the piperazine group inserts deeply in the membrane and is accompanied by a stabilization of its neutral form. The bulky nature of rifampicin and its deep insertion in the membrane lead to a strong perturbation in the lipids local order, decreasing the membrane barrier properties as evaluated from the rate of carboxyfluorescein leaching. Altogether, the comparison between the experimental and MD simulations results provides important insight regarding the rifampicin molecular features responsible for its relatively fast membrane permeation.