Dynamical and structural properties of charged and uncharged lidocaine in a lipid bilayer

Micro-pharmacokinetics of lidocaine and bupivacaine transfer across a myelinated nerve fiber

BACKGROUND

The aim of the present study was to predict the time to onset and duration of action of two local anesthetics (lidocaine and bupivacaine) based on experimental dimensions of a typical nerve and experimental octanol/water partition coefficients.

METHODS

We began our compilation of experimental data with a numerical solution of the Smoluchowski equation for the transfer of lidocaine and bupivacaine across the axon membrane in the region of the node of Ranvier (axolemma) and across the Schwann cell. The difference between the aqueous and lipid environments of the neuron was simulated by including the coordinate-dependent chemical potential. In the second step, the permeation rates calculated using the diffusion equation were used to solve a system of four ordinary differential equations. This approach allowed us to simulate the cellular environment for a longer time and to compare our model with pharmacokinetic properties (time to onset and duration of action) of local anesthetics from the literature. The behavior of local anesthetics under physiological conditions and in case of local acidosis was also simulated.

RESULTS

We demonstrated that local anesthetics cross the axolemma in a time span of less than 1 μs. The time to onset of action, controlled by diffusion from the epineurium to an axon with a typical distance of 500 μm, was 167 s and 186 s for lidocaine and bupivacaine, respectively. The calculated half-life, which is a measure of the duration of action, was 41 min and 328 min for lidocaine and bupivacaine, respectively.

CONCLUSIONS

Duration of action is controlled by the storage capacity of lipophilic compartments around the axon, which is higher for bupivacaine but lower in local acidosis. For the latter case, the literature, including textbooks, provides a misinterpretation, namely that protonated species cannot penetrate the membrane.

Protonation and orientation: a computational approach to cocaine diffusion through a model membrane

We study the diffusion of cocaine through a DMPC lipid bilayer as an example of a protonable, amphiphilic molecule passing a biological membrane. Using classical molecular dynamics simulations, the free...

Phosphodiesterase Type 5 Inhibitors Greatly Affect Physicochemical Properties of Model Lipid Membranes

Although phosphodiesterase type 5 inhibitors are widely used and well-studied drugs, the potential benefits of their application in the treatment of various diseases and new drug delivery systems, including liposome forms, are still being discussed. In this regard, the role of the lipid matrix of cell membranes in the pharmacological action of the inhibitors is of special interest. It was shown that sildenafil, vardenafil, and tadalafil caused a significant decrease in the boundary potential of model membranes composed of palmitoyloleoylphosphatidylcholine or its mixture with cholesterol, by 70–80 mV. The reduction in the membrane dipole potential induced by inhibitors led to a 20–25% increase in the conductance of cation-selective pores formed by the antimicrobial peptide gramicidin A. The addition of sildenafil or vardenafil also led to a significant decrease in the temperature of the main phase transition of dipalmytoylphosphatidylcholine, by about 1.5 °C, while tadalafil did not change the melting temperature. Sildenafil, vardenafil, and tadalafil enhanced the pore-forming activity of the antifungal polyene antibiotic nystatin by 11, 13, and 2 times, respectively. This fact might indicate the induction of membrane curvature stress by the inhibitors. The data obtained might be of special interest for the development of lipid-mediated forms of drugs.

Determination of bupivacaine tissue concentration in human biopsy samples using high-performance liquid chromatography with mass spectrometry

In the present study, we developed a simple and rapid analytical method for the quantification of bupivacaine hydrochloride in human biopsy samples of adipose, muscle, neural, connective and cartilage tissue using liquid chromatography-mass spectrometry. Anesthetics were extracted from the tissue samples using 0.1% formic acid in acetonitrile for protein denaturation and hexane for removal of lipophilic impurities. Analytes were separated adequately on Phenomenex Luna Omega polar C18 column using a gradient mobile phase 0.1% formic acid in water and 0.1% formic acid in acetonitrile. The lower limits of quantification were ≤ 97 ng g^-1 tissue for all studied tissues. Intra-day recoveries were between 48.2% and 82.1% with relative standard deviations (RSDs) between 1.47% and 14.28%, whereas inter-day recoveries were between 52.2% and 77.6% with RSDs between 2.98% and 14.79%. The calibration curve showed a linear fit with R² higher than 0.99 in the concentration range from 0.16 to 100 μg g^-1 . Lidocaine hydrochloride was tested as internal standard because its recoveries and matrix effects were comparable to bupivacaine hydrochloride. Post-analytical corrections of measured bupivacaine tissue concentrations can accordingly be made based on recovery of lidocaine as internal standard, with recoveries between 51.2% and 86.9% and RSDs between 1.99% and 16.88%. The developed method could be used to study time-dependent spread of bupivacaine locally or to more distant locations across tissue barriers.

Alkaloids Modulate the Functioning of Ion Channels Produced by Antimicrobial Agents via an Influence on the Lipid Host

It is widely recognized that an alteration in membrane physical properties induced by the adsorption of various drugs and biologically active compounds might greatly affect the functioning of peptides and proteins embedded in the membrane, in particular various ion channels. This study aimed to obtain deep insight into the diversity of the molecular mechanisms of membrane action of one of the most numerous and extremely important class of phytochemicals, the alkaloids. Protoalkaloids (derivatives of β-phenylethylamine, benzylamines, and colchicines), heterocyclic alkaloids (derivatives of purine, quinolysidine, piperidine, pyridine, quinoline, and isoquinoline), and steroid alkaloids were tested. We evaluated the effects of 22 compounds on lipid packing by investigating the thermotropic behavior of membrane lipids and the leakage of a fluorescent marker from unilamellar lipid vesicles. The alteration in the transmembrane distribution of the electrical potential was estimated by measuring the alkaloid induced changes in the boundary potential of planar lipid bilayers. We found that benzylamines, the chili pepper active components, capsaicin and dihydrocapsaicin, strongly affect not only the elastic properties of the lipid host, but also its electrostatics by dramatic decrease in membrane dipole potential. We concluded that the increase in the conductance and lifetime of gramicidin A channels induced by benzylamines was related to alteration in membrane dipole potential not to decrease in membrane stiffness. A sharp decrease in the lifetime of single ion pores induced by the antifungal lipopeptide syringomycin E, after addition of benzylamines and black pepper alkaloid piperine, was also mainly due to the reduction in dipole potential. At the same time, we showed that the disordering of membrane lipids in the presence of benzylamines and piperine plays a decisive role in the regulation of the conductance induced by the antifungal polyene macrolide antibiotic nystatin, while the inhibition of steady-state transmembrane current produced by the antimicrobial peptide cecropin A was attributed to both the dipole potential drop and membrane lipid disordering in the presence of pepper alkaloids. These data might lead to a better understanding of the biological activity of alkaloids, especially their action on voltage-gated and mechanosensitive ion channels in cell membranes.

Atomistic Molecular Dynamics Simulations of Propofol and Fentanyl in Phosphatidylcholine Lipid Bilayers

Atomistic molecular dynamics (MD) and steered MD simulations in combination with umbrella sampling methodology were utilized to study the general anesthetic propofol and the opioid analgesic fentanyl and their interaction with lipid bilayers, which is not yet fully understood. These molecules were inserted into two different fully hydrated phospholipid bilayers, namely, dioleoylphosphatidylcholine (DOPC) and dipalmitoylphosphatidylcholine (DPPC), to investigate the effects that these drugs have on the bilayer. We determined the role of the lipid chain length and saturation on the behavior of the two drugs. Pure, fully hydrated DOPC and DPPC bilayers were also simulated, and the results were in excellent agreement with the experimental values. Various structural and mechanical properties of each system, such as the area per lipid, area compressibility modulus, order parameter, lateral lipid diffusion, hydrogen bonds, and radial distribution functions, have been calculated to assess how the drug molecules affect the different bilayers. From the calculated results, we show that fentanyl and propofol generally follow similar trends in each bilayer but adopt different favorable positions close to the headgroup/chain interface at the carbonyl groups. Propofol was shown to selectively form hydrogen bonds at the carbonyl carbon in each bilayer, whereas fentanyl interacts with water molecules at the headgroup interface. From the calculated free-energy profiles, we determined that both molecules show a preference for the low-density, low-order acyl chain region of the bilayers and both significantly preferred the DOPC bilayer with propofol and fentanyl having energy minima at -6.66 and -43.07 kcal mol^-1, respectively. This study suggests that different chain lengths and levels of saturation directly affect the properties of these two important molecules, which are seen to work together to control anesthesia in surgical applications.

Appropriate amounts proportions of lidocaine gel, indigo carmine and lipiodol mixture for preoperative marking in video-assisted thoracic surgery

Preoperative marking is necessary in thoracoscopic wedge resections of lung nodules expected to be invisible or nonpalpable during surgery. Recently, lidocaine gel was added to a dye solution containing indigo carmine and lipiodol to promote micelle formation, but the optimal mixing ratio was not determined. Thus, an in vitro experiment was performed to identify the optimal mixing ratio of lidocaine gel, indigo carmine and lipiodol. To fixed volumes of indigo carmine and lipiodol of 0.5 ml each, 0.1, 0.2, 0.3, 0.4 and 0.5 ml of lidocaine gel was added. Changes were examined every 2 h. No changes were seen in the mixtures containing 0.1 and 0.2 ml of lidocaine gel, whereas those with 0.3, 0.4 and 0.5 ml had clearly separated after 2 h. Our findings suggest that the ideal proportion of indigo carmine, lipiodol and lidocaine gel for use in preoperative marking is 5:5:2.

Diffusional interaction behavior of NSAIDs in lipid bilayer membrane using molecular dynamics (MD) simulation: Aspirin and Ibuprofen

In this research, for the first time, molecular dynamics (MD) method was used to simulate aspirin and ibuprofen at various concentrations and in neutral and charged states. Effects of the concentration (dosage), charge state, and existence of an integral protein in the membrane on the diffusion rate of drug molecules into lipid bilayer membrane were investigated on 11 systems, for which the parameters indicating diffusion rate and those affecting the rate were evaluated. Considering the diffusion rate, a suitable score was assigned to each system, based on which, analysis of variance (ANOVA) was performed. By calculating the effect size of the indicative parameters and total scores, an optimum system with the highest diffusion rate was determined. Consequently, diffusion rate controlling parameters were obtained: the drug-water hydrogen bond in protein-free systems and protein-drug hydrogen bond in the systems containing protein.

Protonation state of inhibitors determines interaction sites within voltage-gated sodium channels

Voltage-gated sodium channels are essential for carrying electrical signals throughout the body, and mutations in these proteins are responsible for a variety of disorders, including epilepsy and pain syndromes. As such, they are the target of a number of drugs used for reducing pain or combatting arrhythmias and seizures. However, these drugs affect all sodium channel subtypes found in the body. Designing compounds to target select sodium channel subtypes will provide a new therapeutic pathway and would maximize treatment efficacy while minimizing side effects. Here, we examine the binding preferences of nine compounds known to be sodium channel pore blockers in molecular dynamics simulations. We use the approach of replica exchange solute tempering (REST) to gain a more complete understanding of the inhibitors' behavior inside the pore of NavMs, a bacterial sodium channel, and NavPas, a eukaryotic sodium channel. Using these simulations, we are able to show that both charged and neutral compounds partition into the bilayer, but neutral forms more readily cross it. We show that there are two possible binding sites for the compounds: (i) a site on helix 6, which has been previously determined by many experimental and computational studies, and (ii) an additional site, occupied by protonated compounds in which the positively charged part of the drug is attracted into the selectivity filter. Distinguishing distinct binding poses for neutral and charged compounds is essential for understanding the nature of pore block and will aid the design of subtype-selective sodium channel inhibitors.

Experimental and theoretical studies of emodin interacting with a lipid bilayer of DMPC

Emodin is one of the most abundant anthraquinone derivatives found in nature. It is the active principle of some traditional herbal medicines with known biological activities. In this work, we combined experimental and theoretical studies to reveal information about location, orientation, interaction and perturbing effects of Emodin on lipid bilayers, where we have taken into account the neutral form of the Emodin (EMH) and its anionic/deprotonated form (EM^-). Using both UV/Visible spectrophotometric techniques and molecular dynamics (MD) simulations, we showed that both EMH and EM^- are located in a lipid membrane. Additionally, using MD simulations, we revealed that both forms of Emodin are very close to glycerol groups of the lipid molecules, with the EMH inserted more deeply into the bilayer and more disoriented relative to the normal of the membrane when compared with the EM^-, which is more exposed to interfacial water. Analysis of several structural properties of acyl chains of the lipids in a hydrated pure DMPC bilayer and in the presence of Emodin revealed that both EMH and EM^- affect the lipid bilayer, resulting in a remarkable disorder of the bilayer in the vicinity of the Emodin. However, the disorder caused by EMH is weaker than that caused by EM^-. Our results suggest that these disorders caused by Emodin might lead to distinct effects on lipid bilayers including its disruption which are reported in the literature.

Study of procaine and tetracaine in the lipid bilayer using molecular dynamics simulation

Despite available experimental results, the molecular mechanism of action of local anesthetics upon the nervous system and contribution of the cell membrane to the process are still controversial. In this work, molecular dynamics simulations were performed to investigate the effect of two clinically used local anesthetics, procaine and tetracaine, on the structure and dynamics of a fully hydrated dimyristoylphosphatidylcholine lipid bilayer. We focused on comparing the main effects of uncharged and charged drugs on various properties of the lipid membrane: mass density distribution, diffusion coefficient, order parameter, radial distribution function, hydrogen bonding, electrostatic potential, headgroup angle, and water dipole orientation. To compare the diffusive nature of anesthetic through the lipid membrane quantitatively, we investigated the hexadecane/water partition coefficient using expanded ensemble simulation. We predicted the permeability coefficient of anesthetics in the following order: uncharged tetracaine > uncharged procaine > charged tetracaine > charged procaine. We also shown that the charged forms of drugs are more potent in hydrogen bonding, disturbing the lipid headgroups, changing the orientation of water dipoles, and increasing the headgroup electrostatic potential more than uncharged drugs, while the uncharged drugs make the lipid diffusion faster and increase the tail order parameter. The results of these simulation studies suggest that the different forms of anesthetics induce different structural modifications in the lipid bilayer, which provides new insights into their molecular mechanism.

Naratriptan aggregation in lipid bilayers: perspectives from molecular dynamics simulations

In order to understand the interaction between naratriptan and a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline (POPC), we carried out molecular dynamics simulations. The simulations were performed considering neutral and protonated ionization states, starting from different initial conditions. At physiological pH, the protonated state of naratriptan is predominant. It is expected that neutral compounds could have larger membrane partition than charged compounds. However, for the specific case of triptans, it is difficult to study neutral species in membranes experimentally, making computer simulations an interesting tool. When the naratriptan molecules were originally placed in water, they partitioned between the bilayer/water interface and water phase, as has been described for similar compounds. From this condition, the drugs displayed low access to the hydrophobic environment, with no significant effects on bilayer organization. The molecules anchored in the interface, due mainly to the barrier function of the polar and oriented lipid heads. On the other hand, when placed inside the bilayer, both neutral and protonated naratriptan showed self-aggregation in the lipid tail environment. In particular, the protonated species exhibited a pore-like structure, dragging water through this environment. Graphical Abstract Different behaviour of Naratriptan and Sumatriptan, when the drugs were originally placed in the lipid core.

In silico pharmacology: Drug membrane partitioning and crossing

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.

From molecular modelling to photophysics of neutral oligo- and polyfluorenes incorporated into phospholipid bilayers

The combination of various experimental techniques with theoretical simulations has allowed elucidation of the mode of incorporation of fluorene based derivatives into phospholipid bilayers. Molecular dynamics (MD) simulations on a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) bilayer, with benzene (B), biphenyl (BP), fluorene (F) and tri-(9,9-di-n-octylfluorenyl-2,7-diyl), TF, have provided insights into the topography of these molecules when they are present in the phospholipid bilayer, and suggest marked differences between the behavior of the small molecules and the oligomer. Further information on the interaction of neutral fluorenes within the phospholipid bilayer was obtained by an infrared (IR) spectroscopic study of films of DMPC and of the phospholipid with PFO deuterated specifically on its alkyl chains (DMPC-PFO-d34). This was complemented by measurements of the effect of F, TF and two neutral polymers: polyfluorene poly(9,9-di-n-octylfluorenyl-2,7-diyl), PFO, and poly(9,9-di-n-dodecylfluorenyl-2,7-diyl), PFD, on the phospholipid phase transition temperature using differential scanning calorimetry (DSC). Changes in liposome size upon addition of F and PFO were followed by dynamic light scattering. In addition, the spectroscopic properties of F, TF, PFO and PFD solubilised in DMPC liposomes (absorption, steady-state and time-resolved fluorescence) were compared with those of the same probes in typical organic solvents (chloroform, cyclohexane and ethanol). Combining the insight from MD simulations with the results at the molecular level from the various experimental techniques suggests that while the small molecules have a tendency to be located in the phospholipid head group region, the polymers are incorporated within the lipid bilayers, with the backbone predominantly orthogonal to the phospholipid alkyl chains and with interdigitation of them and the PFO alkyl chains.

Anesthetic diffusion through lipid membranes depends on the protonation rate

Hundreds of substances possess anesthetic action. However, despite decades of research and tests, a golden rule is required to reconcile the diverse hypothesis behind anesthesia. What makes an anesthetic to be local or general in the first place? The specific targets on proteins, the solubility in lipids, the diffusivity, potency, action time? Here we show that there could be a new player equally or even more important to disentangle the riddle: the protonation rate. Indeed, such rate modulates the diffusion speed of anesthetics into lipid membranes; low protonation rates enhance the diffusion for local anesthetics while high ones reduce it. We show also that there is a pH and membrane phase dependence on the local anesthetic diffusion across multiple lipid bilayers. Based on our findings we incorporate a new clue that may advance our understanding of the anesthetic phenomenon.

Bupivacaine-induced cellular entry of QX-314 and its contribution to differential nerve block

BACKGROUND AND PURPOSE

Selective nociceptor fibre block is achieved by introducing the cell membrane impermeant sodium channel blocker lidocaine N-ethyl bromide (QX-314) through transient receptor potential V1 (TRPV1) channels into nociceptors. We screened local anaesthetics for their capacity to activate TRP channels, and characterized the nerve block obtained by combination with QX-314.

EXPERIMENTAL APPROACH

We investigated TRP channel activation in dorsal root ganglion (DRG) neurons by calcium imaging and patch-clamp recordings, and cellular QX-314 uptake by MS. To characterize nerve block, compound action potential (CAP) recordings from isolated nerves and behavioural responses were analysed.

KEY RESULTS

Of the 12 compounds tested, bupivacaine was the most potent activator of ruthenium red-sensitive calcium entry in DRG neurons and activated heterologously expressed TRPA1 channels. QX-314 permeated through TRPA1 channels and accumulated intracellularly after activation of these channels. Upon sciatic injections, QX-314 markedly prolonged bupivacaine's nociceptive block and also extended (to a lesser degree) its motor block. Bupivacaine's blockade of C-, but not A-fibre, CAPs in sciatic nerves was extended by co-application of QX-314. Surprisingly, however, this action was the same in wild-type, TRPA1-knockout and TRPV1/TRPA1-double knockout mice, suggesting a TRP-channel independent entry pathway. Consistent with this, high doses of bupivacaine promoted a non-selective, cellular uptake of QX-314.

CONCLUSIONS AND IMPLICATIONS

Bupivacaine, combined with QX-314, produced a long-lasting sensory nerve block. This did not require QX-314 permeation through TRPA1, although bupivacaine activated these channels. Regardless of entry pathway, the greatly extended duration of block produced by QX-314 and bupivacaine may be clinically useful.

The interaction of local anesthetics with lipid membranes

Molecular Dynamic Simulations are performed to evaluate the interaction of lidocaine, procaine and tetracaine with a lipid membrane. The main interest is to evaluate the structural changes produced by these local anesthetics in the bilayers. Penetration trajectories, interaction energies, entropy changes and an order parameter are calculated to quantify the destabilization of the lipid configurations. We show that such structural parameters give important information to understand how anesthetic agents influence the structure of plasma membranes. Graphic processing units (GPUs) are used in our simulations.

Conformational dynamics of dry lamellar crystals of sugar based lipids: an atomistic simulation study

The rational design of a glycolipid application (e.g. drug delivery) with a tailored property depends on the detailed understanding of its structure and dynamics. Because of the complexity of sugar stereochemistry, we have undertaken a simulation study on the conformational dynamics of a set of synthetic glycosides with different sugar groups and chain design, namely dodecyl β-maltoside, dodecyl β-cellobioside, dodecyl β-isomaltoside and a C12C10 branched β-maltoside under anhydrous conditions. We examined the chain structure in detail, including the chain packing, gauche/trans conformations and chain tilting. In addition, we also investigated the rotational dynamics of the headgroup and alkyl chains. Monoalkylated glycosides possess a small amount of gauche conformers (∼20%) in the hydrophobic region of the lamellar crystal (LC) phase. In contrast, the branched chain glycolipid in the fluid Lα phase has a high gauche population of up to ∼40%. Rotational diffusion analysis reveals that the carbons closest to the headgroup have the highest correlation times. Furthermore, its value depends on sugar type, where the rotational dynamics of an isomaltose was found to be 11-15% and more restrained near the sugar, possibly due to the chain disorder and partial inter-digitation compared to the other monoalkylated lipids. Intriguingly, the present simulation demonstrates the chain from the branched glycolipid bilayer has the ability to enter into the hydrophilic region. This interesting feature of the anhydrous glycolipid bilayer simulation appears to arise from a combination of lipid crowding and the amphoteric nature of the sugar headgroups.

Molecular dynamics simulation of the partitioning of benzocaine and phenytoin into a lipid bilayer

Molecular dynamics simulations were used to examine the partitioning behaviour of the local anaesthetic benzocaine and the anti-epileptic phenytoin into lipid bilayers, a factor that is critical to their mode of action. Free energy methods are used to quantify the thermodynamics of drug movement between water and octanol as well as for permeation across a POPC membrane. Both drugs are shown to favourably partition into the lipid bilayer from water and are likely to accumulate just inside the lipid headgroups where they may alter bilayer properties or interact with target proteins. Phenytoin experiences a large barrier to cross the centre of the bilayer due to less favourable energetic interactions in this less dense region of the bilayer. Remarkably, in our simulations both drugs are able to pull water into the bilayer, creating water chains that extend back to bulk, and which may modify the local bilayer properties. We find that the choice of atomic partial charges can have a significant impact on the quantitative results, meaning that careful validation of parameters for new drugs, such as performed here, should be performed prior to their use in biomolecular simulations.

Membrane permeable local anesthetics modulate Na(V)1.5 mechanosensitivity

Voltage-gated sodium selective ion channel Na_V1.5 is expressed in the heart and the gastrointestinal tract, which are mechanically active organs. Na_V1.5 is mechanosensitive at stimuli that gate other mechanosensitive ion channels. Local anesthetic and antiarrhythmic drugs act upon Na_V1.5 to modulate activity by multiple mechanisms. This study examined whether Na_V1.5 mechanosensitivity is modulated by local anesthetics. Na_V1.5 channels wereexpressed in HEK-293 cells, and mechanosensitivity was tested in cell-attached and excised inside-out configurations. Using a novel protocol with paired voltage ladders and short pressure pulses, negative patch pressure (-30 mmHg) in both configurations produced a hyperpolarizing shift in the half-point of the voltage-dependence of activation (V_1/2a) and inactivation (V_1/2i) by about -10 mV. Lidocaine (50 µM) inhibited the pressure-induced shift of V_1/2a but not V_1/2i. Lidocaine inhibited the tonic increase in pressure-induced peak current in a use-dependence protocol, but it did not otherwise affect use-dependent block. The local anesthetic benzocaine, which does not show use-dependent block, also effectively blocked a pressure-induced shift in V_1/2a. Lidocaine inhibited mechanosensitivity in Na_V1.5 at the local anesthetic binding site mutated (F1760A). However, a membrane impermeable lidocaine analog QX-314 did not affect mechanosensitivity of F1760A Na_V1.5 when applied from either side of the membrane. These data suggest that the mechanism of lidocaine inhibition of the pressure-induced shift in the half-point of voltage-dependence of activation is separate from the mechanisms of use-dependent block. Modulation of Na_V1.5 mechanosensitivity by the membrane permeable local anesthetics may require hydrophobic access and may involve membrane-protein interactions.

Intramolecular hydrogen bonding in articaine can be related to superior bone tissue penetration: a molecular dynamics study

Local anesthetics (LAs) are drugs that cause reversible loss of nociception during surgical procedures. Articaine is a commonly used LA in dentistry that has proven to be exceptionally effective in penetrating bone tissue and induce anesthesia on posterior teeth in maxilla and mandibula. In the present study, our aim was to gain a deeper understanding of the penetration of articaine through biological membranes by studying the interactions of articaine with a phospholipid membrane. Our approach involves Langmuir monolayer experiments combined with molecular dynamics simulations. Membrane permeability of LAs can be modulated by pH due to a titratable amine group with a pKa value close to physiological pH. A change in protonation state is thus known to act as a lipophilicity switch in LAs. Our study shows that articaine has an additional unique lipophilicity switch in its ability to form an intramolecular hydrogen bond. We suggest this intramolecular hydrogen bond as a novel and additional solvent-dependent mechanism for modulation of lipophilicity of articaine which may enhance its diffusion through membranes and connective tissue.

Partitioning and localization of environment-sensitive 2-(2'-pyridyl)- and 2-(2'-pyrimidyl)-indoles in lipid membranes: a joint refinement using fluorescence measurements and molecular dynamics simulations

Fluorescence of environment-sensitive dyes is widely applied to monitor local structure and solvation dynamics of biomolecules. It has been shown that, in comparison with a parent indole fluorophore, fluorescence of 2-(2'-pyridyl)-5-methylindole (5M-PyIn-0) and 2-[2'-(4',6'-dimethylpyrimidyl)]-indole (DMPmIn-0) is remarkably sensitive to hydrogen bonding with protic partners. Strong fluorescence, observed for these compounds in nonpolar and polar aprotic solvents, is efficiently quenched in aqueous solution. This study demonstrates that 5M-PyIn-0 and DMPmIn-0, which are almost nonemitting in aqueous solution, become highly fluorescent upon titrating with phospholipid vesicles. The fluorescence enhancement is accompanied by a significant blue shift of emission maximum. The Gibbs free energy of membrane partitioning, measured by the increase in the steady-state fluorescence intensities during transfer from an aqueous environment to a lipid bilayer, is very favorable for both compounds, being in a range from -7.1 to -8.0 kcal/mol and depending only slightly on lipid composition of the membrane. The fluorescence enhancement upon membrane partitioning is indicative of the loss of the specific hydrogen-bonding interactions between the excited fluorophore and water molecules, causing efficient fluorescence quenching in bulk water. This conclusion is supported by atomistic molecular dynamics (MD) simulations, demonstrating that both 5M-PyIn-0 and DMPmIn-0 bind rapidly and partition deeply into a lipid bilayer. MD simulations also show a rapid, nanosecond-scale decrease in the probability of solute-solvent hydrogen bonding during passive diffusion of the probe molecules from bulk water into a lipid bilayer. At equilibrium conditions, both 5M-PyIn-0 and DMPmIn-0 prefer deep localization within the hydrophobic, water-free region of the bilayer. A free energy profile of penetration across a bilayer estimated using MD umbrella sampling shows that both indole derivatives favor residence in a rather wide potential energy well located 10-15 Å from the bilayer center.

Distribution and dynamics of adamantanes in a lipid bilayer

The adamantanes are a class of compounds that have found use in the treatment of influenza A and Parkinson's disease, among others. The mode of action for influenza A is based on the adamantanes' interaction with the transmembrane M2 channel, whereas the treatment of Parkinson's disease is thought to relate to a channel block of N-methyl-D-aspartate receptors. An understanding of how these compounds interact with the lipid bilayer is thus of great interest. We used molecular-dynamics simulations to calculate the potential of mean force of adamantanes in a lipid bilayer. Our results demonstrate a preference for the interfacial region of the lipid bilayer for both protonated and deprotonated species, with the protonated species proving significantly more favorable. However, the protonated species have a large free-energy barrier in the center of the membrane. In contrast, there is no barrier (compared with aqueous solution) at the center of the bilayer for deprotonated species, suggesting that the permeant species is indeed the neutral form, as commonly assumed. We discuss the results with respect to proposed mechanisms of action and implications for drug-delivery in general.

Membrane interactivity of charged local anesthetic derivative and stereoselectivity in membrane interaction of local anesthetic enantiomers

With respect to the membrane lipid theory as a molecular mechanism for local anesthetics, two critical subjects, the negligible effects of charged drugs when applied extracellularly and the stereoselective effects of enantiomers, were verified by paying particular attention to membrane components, phospholipids with the anionic property, and cholesterol with several chiral carbons. The membrane interactivities of structurally-different anesthetics were determined by their induced fluidity changes of liposomal membranes. Lidocaine (3.0 μmol/mL) fluidized phosphatidylcholine membranes, but not its quaternary derivative QX-314 (3.0 μmol/mL). Similarly to the mother molecule lidocaine, however, QX-314 fluidized phosphatidylserine-containing nerve cell model membranes and acidic phospholipids-constituting membranes depending on the acidity of membrane lipids. Positively charged local anesthetics are able to act on lipid bilayers by ion-pairing with anionic (acidic) phospholipids. Bupivacaine (0.75 mol/mL) and ropivacaine (0.75 and 1.0 μmol/mL) fluidized nerve cell model membranes with the potency being S(−)-enantiomer < racemate < R(+)-enantiomer (P < 0.01, vs antipode and racemate) and cardiac cell model membranes with the potency being S(−)-ropivacaine < S(−)-bupivacaine < R(+)-bupivacaine (P < 0.01). However, their membrane effects were not different when removing cholesterol from the model membranes. Stereoselectivity is producible by cholesterol which increases the chirality of lipid bilayers and enables to discriminate anesthetic enantiomers. The membrane lipid interaction should be reevaluated as the mode of action of local anesthetics.

NMR investigations of interactions between anesthetics and lipid bilayers

Interactions between anesthetics (lidocaine and short chain alcohols) and lipid membranes formed by dimyristoylphosphatidylcholine (DMPC) were studied using NMR spectroscopy. The orientational order of lidocaine was investigated using deuterium NMR on a selectively labelled compound whereas segmental ordering in the lipids was probed by two-dimensional 1H-13C separated local field experiments under magic-angle spinning conditions. In addition, trajectories generated in molecular dynamics (MD) computer simulations were used for interpretation of the experimental results. Separate simulations were carried out with charged and uncharged lidocaine molecules. Reasonable agreement between experimental dipolar interactions and the calculated counterparts was observed. Our results clearly show that charged lidocaine affects significantly the lipid headgroup. In particular the ordering of the lipids is increased accompanied by drastic changes in the orientation of the P-N vector in the choline group.

Distribution and favorable binding sites of pyrroloquinoline and its analogues in a lipid bilayer studied by molecular dynamics simulations

The distribution of 1H-pyrrolo[3,2-h]quinoline (PQ), 11H-dipyrido[2,3-a]carbazole (PC) and 7-azaindole (7AI) at a water/membrane interface has been investigated by molecular dynamics (MD) simulations. The MD study focused on favorable binding sites of the azaaromatic probes across a dipalmitoylphosphatidylcholine (DPPC) bilayer. Our simulations show that PQ and PC are preferably accommodated at the hydrocarbon core of the bilayer below the glycerol moiety. In addition, it is found that the hydrophobic aromatic parts of the probes are located inside a more ordered region of DPPC, consisting of hydrophobic lipid chains. In contrast to PQ and PC, 7AI is characterized by a broad distribution between a DPPC interface and water, so that the three preferable binding sites are found across a water/membrane interface. It is found that in the sequence 7AI-PQ-PC, due to the increase of the number of aromatic rings and, hence, the hydrophobic character of the probes, the depth of the probe localization is gradually shifted deeper inside the hydrocarbon core of the bilayer. We found that the probe-lipid hydrogen-bonding contributes weakly to the favorable localizations of the azaaromatic probes inside the DPPC bilayer, so that the probe localization is mainly driven by electrostatic dipole-dipole and van der Waals interactions.

Effect of local anesthetic lidocaine on electrostatic properties of a lipid bilayer

The influence of the local anesthetic lidocaine on electrostatic properties of a lipid membrane bilayer was studied by molecular dynamics simulations. The electrostatic dipole potential, charge densities, and orientations of the headgroup angle have been examined in the presence of different amounts of charged or uncharged forms of lidocaine. Important changes in the membrane properties caused by the presence of both forms of lidocaine are presented and discussed. Our simulations have shown that both charged and uncharged lidocaine cause almost the same increase in the electrostatic potential in the middle of the membrane, although for different reasons. The increase, approximately 90 mV for 9 mol % of lidocaine and 220 mV for 28 mol % of lidocaine, is of a size that may affect the functioning of voltage-gated ion channels.