Structure-based prediction of drug distribution across the headgroup and core strata of a phospholipid bilayer using surrogate phases

Application of Generative Artificial Intelligence in Predicting Membrane Partitioning of Drugs: Combining Denoising Diffusion Probabilistic Models and MD Simulations Reduces the Computational Cost to One-Third

The optimal interaction of drugs with plasma membranes and membranes of subcellular organelles is a prerequisite for desirable pharmacology. Importantly, for drugs targeting the transmembrane lipid-facing sites of integral membrane proteins, the relative affinity of a drug to the bilayer lipids compared to the surrounding aqueous phase affects the partitioning, access, and binding of the drug to the target site. Molecular dynamics (MD) simulations, including enhanced sampling techniques such as steered MD, umbrella sampling (US), and metadynamics, offer valuable insights into the interactions of drugs with the membrane lipids and water in atomistic detail. However, these methods are computationally prohibitive for the high-throughput screening of drug candidates. This study shows that applying denoising diffusion probabilistic models (DDPMs), a generative AI method, to US simulation data reduces the computational cost significantly. Specifically, the models used only partial (one-third) data from the US simulations and reproduced the complete potential of mean force (PMF) profiles for three FDA-approved drugs (β2-adrenergic agonists) and ∼20 biologically relevant chemicals with known experimentally characterized bilayer locations. Intriguingly, the model can predict the solvation-free energies for partitioning and crossing the bilayer, preferred bilayer locations (low-energy well), and orientations of the ligands with high accuracy. The results indicate that DDPMs can be used to characterize the complete membrane partitioning profile of drug molecules using fewer umbrella sampling simulations at select positions along the bilayer normal (z-axis), irrespective of their amphiphilic-lipophilic-cephalophilic characteristics.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Antileishmanial activity of synthetic analogs of the naturally occurring quinolone alkaloid N-methyl-8-methoxyflindersin

Leishmaniasis is a neglected, parasitic tropical disease caused by an intracellular protozoan from the genus Leishmania. Quinoline alkaloids, secondary metabolites found in plants from the Rutaceae family, have antiparasitic activity against Leishmania sp. N-methyl-8-methoxyflindersin (1), isolated from the leaves of Raputia heptaphylla and also known as 7-methoxy-2,2-dimethyl-2H,5H,6H-pyran[3,2-c]quinolin-5-one, shows antiparasitic activity against Leishmania promastigotes and amastigotes. This study used in silico tools to identify synthetic quinoline alkaloids having structure similar to that of compound 1 and then tested these quinoline alkaloids for their in vitro antiparasitic activity against Leishmania (Viannia) panamensis, in vivo therapeutic response in hamsters suffering from experimental cutaneous leishmaniasis (CL), and ex vivo immunomodulatory potential in healthy donors' human peripheral blood (monocyte)-derived macrophages (hMDMs). Compounds 1 (natural), 2 (synthetic), and 8 (synthetic) were effective against intracellular promastigotes (9.9, 3.4, and 1.6 μg/mL medial effective concentration [EC50], respectively) and amastigotes (5.07, 7.94, and 1.91 μg/mL EC50, respectively). Compound 1 increased nitric oxide production in infected hMDMs and triggered necrosis-related ultrastructural alterations in intracellular amastigotes, while compound 2 stimulated oxidative breakdown in hMDMs and caused ultrastructural alterations in the parasite 4 h posttreatment, and compound 8 failed to induce macrophage modulation but selectively induced apoptosis of infected hMDMs and alterations in the intracellular parasite ultrastructure. In addition, synthetic compounds 2 and 8 improved the health of hamsters suffering from experimental CL, without evidence of treatment-associated adverse toxic effects. Therefore, synthetic compounds 2 and 8 are potential therapeutic candidates for topical treatment of CL.

A relationship between membrane permeation and partitioning of nitroaromatic explosives and their functional groups. A computational study

Nitroaromatic explosives, such as 2,4,6-trinitrotoluene, are representative aromatic compounds, which are generally highly toxic.

Polarity of Hydrated Phosphatidylcholine Headgroups

The headgroup (H) stratum (sometimes called the polar region) of membrane bilayers is a relevant yet poorly understood solvation phase for small molecules and macromolecules interacting with the membranes. Solvation of compounds in bilayer strata is characterized experimentally by wide- and small-angle X-ray scattering, neutron diffraction, and various NMR techniques. The quantification is tedious and only available for a limited set of small molecules. Our recently published model of liposome partitioning of small molecules shows that solvation of compounds in the H-stratum of fluid phosphatidylcholine (PC) bilayers correlates well with their solvation in hydrated diacetyl phosphatidylcholine (DAcPC), and solvation in the core (C) depends in a similar way on that in n-hexadecane. These two correlations became a basis for a model describing the location of compounds in the H- and C-strata and at the connecting interface as a nonlinear function of the fragment solvation characteristics of the compounds. In this study, refractivity of hydrated DAcPC phases with varying water contents was measured and polarity was determined using the steady-state fluorescence of indole and Nile Red. The results were compared with the published data obtained by other techniques for PC bilayers in liposomes or on solid supports. The demonstrated qualitative agreement, as well as the polarity and refractivity dependencies on the DAcPC concentration, supports the suitability of hydrated DAcPC as the H-stratum surrogate. Interestingly, depending on hydrations typical for the H-strata of fluid PC bilayers, the dielectric constant could decrease significantly from 31.0 to 7.3 for 16 and 8 water molecules per headgroup, respectively. Although additional experiments are needed for confirmation, this observation could help set proper dielectric constant magnitudes in continuum-based computational models of accumulation and crossing of the PC bilayers with varying hydration levels thanks to the temperature or the structure of fatty acid chains.

Does the Lipid Bilayer Orchestrate Access and Binding of Ligands to Transmembrane Orthosteric/Allosteric Sites of G Protein-Coupled Receptors?

The ligand-binding sites of many G protein-coupled receptors (GPCRs) are situated around and deeply embedded within the central pocket formed by their seven transmembrane-spanning α-helical domains. Generally, these binding sites are assumed accessible to endogenous ligands from the aqueous phase. Recent advances in the structural biology of GPCRs, along with biophysical and computational studies, suggest that amphiphilic and lipophilic molecules may gain access to these receptors by first partitioning into the membrane and then reaching the binding site via lateral diffusion through the lipid bilayer. In addition, several crystal structures of class A and class B GPCRs bound to their ligands offer unprecedented details on the existence of lipid-facing allosteric binding sites outside the transmembrane helices that can only be reached via lipid pathways. The highly organized structure of the lipid bilayer may direct lipophilic or amphiphilic drugs to a specific depth within the bilayer, changing local concentration of the drug near the binding site and affecting its binding kinetics. Additionally, the constraints of the lipid bilayer, including its composition and biophysical properties, may play a critical role in "pre-organizing" ligand molecules in an optimal orientation and conformation to facilitate receptor binding. Despite its clear involvement in molecular recognition processes, the critical role of the membrane in binding ligands to lipid-exposed transmembrane binding sites remains poorly understood and warrants comprehensive investigation. Understanding the mechanistic basis of the structure-membrane interaction relationship of drugs will not only provide useful insights about receptor binding kinetics but will also enhance our ability to take advantage of the apparent membrane contributions when designing drugs that target transmembrane proteins with improved efficacy and safety. In this minireview, we summarize recent structural and computational studies on membrane contributions to binding processes, elucidating both lipid pathways of ligand access and binding mechanisms for several orthosteric and allosteric ligands of class A and class B GPCRs.

Scalable Synthesis and Purification of Acetylated Phosphatidyl Choline Headgroup

The acetylated headgroup of the most abundant mammalian phospholipid, 1,2-diacetyl-3-sn-phosphatidyl choline (DAcPC), has several important applications in research. For instance, it can be dissolved in the same amount of water as in the fluid PC bilayer, to create a surrogate of a PC headgroup stratum, for studying solvation of small molecules and the influence of their structure on the process. In contrast to PC derivatives with longer acyl chains, DAcPC does not self-aggregate, rendering the aqueous solution homogeneous and suitable for simplified analyses of interactions of molecules with the headgroups. Several studies have been published where DAcPC was used in a crudely purified form. Here we describe a one-step preparation of DAcPC from commercially available bulk chemicals and purification of the product by crystallization and washing. The process gives a good yield and is easily scalable. The availability of enantiopure, crystalline DAcPC could open the door to more extensive biochemical, pharmacological, and nutritional studies of this interesting chemical.

Assessing gastric toxicity of xanthone derivatives of anti-inflammatory activity using simulation and experimental approaches

Xanthones are tricyclic compounds of natural or synthetic origin exhibiting a broad spectrum of therapeutic activities. Three synthetic xanthone derivatives (KS1, KS2, and KS3) with properties typical for nonsteroidal anti-inflammatory drugs (NSAID) were objects of the presented model study. NSAIDs are in common use however; several of them exhibit gastric toxicity predominantly resulting from their direct interactions with the outermost lipid layer of the gastric mucosa that impair its hydrophobic barrier property. Among the studied xanthones, gastric toxicity of only KS2 has been determined in previous pharmacological studies, and it is low. In this study, carried out using X-ray diffraction and computer simulation, a palmitoyloleoylphosphatidylcholine-cholesterol bilayer (POPC-Chol) was used as a model of a hydrophobic layer of lipids protecting gastric mucosa as POPC and Chol are the main lipids in human mucus. X-ray diffraction data were used to validate the computer model. The aim of the study was to assess potential gastric toxicity of the xanthones by analysing their atomic level interactions with lipids, ions, and water in the lipid bilayer and their effect on the bilayer physicochemical properties. The results show that xanthones have small effect on the bilayer properties except for its rigidity whereas their interactions with water, ions, and lipids depend on their protonation state and for a given state, are similar for all the xanthones. As gastric toxicity of KS2 is low, based on MD simulations one can predict that toxicity of KS1 and KS3 is also low.

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.