Predictions from First-Principles of Membrane Permeability to Small Molecules: How Useful Are They in Practice?

Elucidating Daptomycin's Antibacterial Efficacy: Insights into the Tripartite Complex with Lipid II and Phospholipids in Bacterial Septum Membrane

This study elucidated the mechanism of formation of a tripartite complex containing daptomycin (Dap), lipid II, and phospholipid phosphatidylglycerol in the bacterial septum membrane, which was previously reported as the cause of the antibacterial action of Dap against gram-positive bacteria via molecular dynamics and enhanced sampling methods. Others have suggested that this transient complex ushers in the inhibition of cell wall synthesis by obstructing the downstream polymerization and cross-linking processes involving lipid II, which is absent in the presence of cardiolipin lipid in the membrane. In this work, we observed that the complex was stabilized by Ca2+-mediated electrostatic interactions between Dap and lipid head groups, hydrophobic interaction, hydrogen bonds, and salt bridges between the lipopeptide and lipids and was associated with Dap concentration-dependent membrane depolarization, thinning of the bilayer, and increased lipid tail disorder. Residues Orn6 and Kyn13, along with the DXDG motif, made simultaneous contact with constituent lipids, hence playing a crucial role in the formation of the complex. Incorporating cardiolipin into the membrane model led to its competitively displacing lipid II away from the Dap, reducing the lifetime of the complex and the nonexistence of lipid tail disorder and membrane depolarization. No evidence of water permeation inside the membrane hydrophobic interior was noted in all of the systems studied. Additionally, it was shown that using hydrophobic contacts between Dap and lipids as collective variables for enhanced sampling gave rise to a free energy barrier for the translocation of the lipopeptide. A better understanding of Dap's antibacterial mechanism, as studied through this work, will help develop lipopeptide-based antibiotics for rising Dap-resistant bacteria.

Impact of the Unstirred Water Layer on the Permeation of Small-Molecule Drugs

Over the last two decades, numerous molecular dynamics (MD) simulation-based investigations have attempted to predict the membrane permeability to small-molecule drugs as indicators of their bioavailability, a majority of which utilize the inhomogeneous solubility diffusion (ISD) model. However, MD-based membrane permeability is routinely 3-4 orders of magnitude larger than the values measured with the intestinal perfusion technique. There have been contentious discussions on the sources of the large discrepancies, and the two indisputable, potentially dominant ones are the fixed protonation state of the permeant and the neglect of the unstirred water layer (UWL). Employing six small-molecule drugs of different biopharmaceutical classification system classes, the current MD study relies on the ISD model but introduces the (de)protonation of the permeant by characterizing the permeation free energy of both neutral and charged states. In addition, the role of the UWL as a potential resistance against permeation is explored. The new MD protocol closely mimics the nature of small-molecule permeation and yields estimates that agree well with in vivo intestinal permeability.

How is Membrane Permeation of Small Ionizable Molecules Affected by Protonation Kinetics?

According to the pH-partition hypothesis, the aqueous solution adjacent to a membrane is a mixture of the ionization states of the permeating molecule at fixed Henderson-Hasselbalch concentrations, such that each state passes through the membrane in parallel with its own specific permeability. An alternative view, based on the assumption that the rate of switching ionization states is instantaneous, represents the permeation of ionizable molecules via an effective Boltzmann-weighted average potential (BWAP). Such an assumption is used in constant-pH molecular dynamics simulations. The inhomogeneous solubility-diffusion framework can be used to compute the pH-dependent membrane permeability for each of these two limiting treatments. With biased WTM-eABF molecular dynamics simulations, we computed the potential of mean force and diffusivity of each ionization state of two weakly basic small molecules: nicotine, an addictive drug, and varenicline, a therapeutic for treating nicotine addiction. At pH = 7, the BWAP effective permeability is greater than that determined by pH-partitioning by a factor of 2.5 for nicotine and 5 for varenicline. To assess the importance of ionization kinetics, we present a Smoluchowski master equation that includes explicitly the protonation and deprotonation processes coupled with the diffusive motion across the membrane. At pH = 7, the increase in permeability due to the explicit ionization kinetics is negligible for both nicotine and varenicline. This finding is reaffirmed by combined Brownian dynamics and Markov state model simulations for estimating the permeability of nicotine while allowing changes in its ionization state. We conclude that for these molecules the pH-partition hypothesis correctly captures the physics of the permeation process. The small free energy barriers for the permeation of nicotine and varenicline in their deprotonated neutral forms play a crucial role in establishing the validity of the pH-partitioning mechanism. Essentially, BWAP fails because ionization kinetics are too slow on the time scale of membrane crossing to affect the permeation of small ionizable molecules such as nicotine and varenicline. For the singly protonated state of nicotine, the computational results agree well with experimental measurements (P1 = 1.29 × 10-7 cm/s), but the agreement for neutral (P0 = 6.12 cm/s) and doubly protonated nicotine (P2 = 3.70 × 10-13 cm/s) is slightly worse, likely due to factors associated with the aqueous boundary layer (neutral form) or leaks through paracellular pathways (doubly protonated form).