In silico pharmacology: Drug membrane partitioning and crossing

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.

Effect of Force Field Resolution on Membrane Mechanical Response and Mechanoporation Damage under Deformation Simulations

Damage induced by transient disruption and mechanoporation in an intact cell membrane is a vital nanoscale biomechanical mechanism that critically affects cell viability. To complement experimental studies of mechanical membrane damage and disruption, molecular dynamics (MD) simulations have been performed at different force field resolutions, each of which follows different parameterization strategies and thus may influence the properties and dynamics of membrane systems. Therefore, the current study performed tensile deformation MD simulations of bilayer membranes using all-atom (AA), united-atom (UA), and coarse-grained Martini (CG-M) models to investigate how the damage biomechanics differs across atomistic and coarse-grained (CG) simulations. The mechanical response and mechanoporation damage were qualitatively similar but quantitatively different in the three models, including some progressive changes based on the coarse-graining level. The membranes yielded and reached ultimate strength at similar strains; however, the coarser systems exhibited lower average yield stresses and failure strains. The average failure strain in the UA model was approximately 7% lower than the AA, and the CG-M was 20% lower than UA and 27% lower than AA. The CG systems also nucleated a higher number of pores and larger pores, which resulted in higher damage during the deformation process. Overall, the study provides insight on the impact of force field-a critical factor in modeling biomolecular systems and their interactions-in inspecting membrane mechanosensitive responses and serves as a reference for justifying the appropriate force field for future studies of more complex membranes and more diverse biomolecular assemblies.

Liposomes characterization for market approval as pharmaceutical products: Analytical methods, guidelines and standardized protocols

Liposomes are nano-sized lipid-based vesicles widely studied for their drug delivery capabilities. Compared to standard carries they exhibit better properties such as improved site-targeting and drug release, protection of drugs from degradation and clearance, and lower toxic side effects. At present, scientific literature is rich of studies regarding liposomes-based systems, while 14 types of liposomal products have been authorized to the market by EMA and FDA and many others have been approved by national agencies. Although the interest in nanodevices and nanomedicine has steadily increased in the last two decades the development of documentation regulating and standardizing all the phases of their development and quality control still suffers from major inadequacy due to the intrinsic complexity of nano-systems characterization. Many generic documents (Type 1) discussing guidelines for the study of nano-systems (lipidic and not) have been proposed while there is a lack of robust and standardized methods (Type 2 documents). As a result, a widespread of different techniques, approaches and methodologies are being used, generating results of variable quality and hard to compare with each other. Additionally, such documents are often subject to updates and rewriting further complicating the topic. Within this context the aim of this work is focused on bridging the gap in liposome characterization: the most recent standardized methodologies suitable for liposomes characterization are here reported (with the corresponding Type 2 documents) and revised in a short and pragmatical way focused on providing the reader with a practical background of the state of the art. In particular, this paper will put the accent on the methodologies developed to evaluate the main critical quality attributes (CQAs) necessary for liposomes market approval.

Permeability of antioxidants through a lipid bilayer model with coarse-grained simulations

Oxidative stress caused by pollution and lifestyle changes causes an excess of free radicals that react chemically with cell constituents leading to irreversible damage. There are molecules known as antioxidants that reduce the levels of free radicals. Some pigments of fruits and vegetables known as anthocyanins have antioxidant properties. Their interaction with the cell membrane becomes a crucial step in studying these substances. In this research, molecular dynamics simulations, particularly, coarse-grained molecular dynamics (CGMD) were used. This technique aims to replace functional groups with corresponding beads that represent their level of polarity and affinities to other chemical groups. Also, umbrella sampling was carried out to obtain free energy profiles that describe well the orientation and location of antioxidants in a membrane considering Trolox, Cyanidin, Gallic Acid, and Resveratrol molecules to study the structural effects they cause on it. It was concluded in this study that an antioxidant when crossing the membrane does not cause either damage to the structural properties or the loss of packing and stratification of phospholipids. it was also observed that the most reactive part of the molecules could easily approach area A prone to lipid oxidation, which can describe the antioxidant capacity of these molecules.Communicated by Ramaswamy H. Sarma.

Visualization of Molecular Permeation into a Multi-compartment Phospholipid Vesicle

Passive permeation of small molecules into vesicles with multiple compartments is a critical event in many chemical and biological processes. We consider the translocation of the peptide NAF-144-67 labeled with a fluorescent fluorescein dye across membranes of rhodamine-labeled 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) into liposomes with internal vesicles. Time-resolved microscopy revealed a sequential absorbance of the peptide in both the outer and inner micrometer vesicles that developed over a time period of minutes to hours, illustrating the spatial and temporal progress of the permeation. There is minimal perturbation of the membrane structure and no evidence for pore formation. On the basis of molecular dynamics simulations of NAF-144-67, we extended a local defect model to migration processes that include multiple compartments. The model captures the long residence time of the peptide within the membrane and the rate of permeation through the liposome and its internal compartments. Imaging experiments confirm the semi-quantitative description of the permeation of the model by activated diffusion and open the way for studies of more complex systems.

Partitioning into phosphatidylcholine-cholesterol membranes: liposome measurements, coarse-grained simulations, and implications for bioaccumulation

Membrane-water partitioning is an important physical property for the assessment of bioaccumulation and environmental impact. Here, we advance simulation methodology for predicting the partitioning of small molecules into lipid membranes and compare the computational predictions to experimental measurements in liposomes. As a step towards high-throughput screening, we present an automated mapping and parametrization procedure to produce coarse-grained models compatible with the Martini 3 force field. The methodology is general and can also be used for other applications where coarse-grained simulations are appropriate. This article addresses the effect on membrane-water partitioning of adding cholesterol to POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Nine contrasting neutral, zwitterionic and charged solutes are tested. Agreement between experiment and simulation is generally good, with the most challenging cases being permanently charged solutes. For all solutes, partitioning is found to be insensitive to membrane cholesterol concentration up to 25% mole fraction. Hence, for assessment of bioaccumulation into a range of membranes (such as those found in fish), partitioning data measured in pure lipid membranes are still informative.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Substrate binding and lipid-mediated allostery in the human organic anion transporter 1 at the atomic-scale

The Organic Anion Transporter 1 is a membrane transporter known for its central role in drug elimination by the kidney. hOAT1 is an antiporter translocating substrate in exchange for a-ketoglutarate. The understanding of hOAT1 structure and function remains limited due to the absence of resolved structure of hOAT1. Benefiting from conserved structural and functional patterns shared with other Major Facilitator Superfamily transporters, the present study intended to investigate fragments of hOAT1 transport function and modulation of its activity in order to make a step forward the understanding of its transport cycle. μs-long molecular dynamics simulation of hOAT1 were carried out suggesting two plausible binding sites for a typical substrate, adefovir, in line with experimental observations. The well-known B-like motif binding site was observed in line with previous studies. However, we here propose a new inner binding cavity which is expected to be involved in substrate translocation event. Binding modes of hOAT1 co-substrate α-ketoglutarate were also investigated suggesting that it may bind to highly conserved intracellular motifs. We here hypothesise that α-ketoglutarate may disrupt the pseudo-symmetrical intracellular charge-relay system which in turn may participate to the destabilisation of OF conformation. Investigations regarding allosteric communications along hOAT1 also suggest that substrate binding event might modulate the dynamics of intracellular charge relay system, assisted by surrounding lipids as active partners. We here proposed a structural rationalisation of transport impairments observed for two single nucleotide polymorphisms, p.Arg50His and p.Arg454Gln suggesting that the present model may be used to transport dysfunctions arising from hOAT1 mutations.

Current Trends and Changes in Use of Membrane Molecular Dynamics Simulations within Academia and the Pharmaceutical Industry

There has been an almost exponential increase in the use of molecular dynamics simulations in basic research and industry over the last 5 years, with almost a doubling in the number of publications each year. Many of these are focused on neurological membranes, and biological membranes in general, applied to the medical industry. A smaller portion have utilized membrane simulations to answer more basic questions related to the function of specific proteins, chemicals or biological processes. This review covers some newer studies, alongside studies from the last two decades, to determine changes in the field. Some of these are basic, while others are more profound, such as multi-component embedded membrane machinery. It is clear that many facets of the discipline remain the same, while the focus on and uses of the technology are broadening in scope and utilization as a general research tool. Analysis of recent literature provides an overview of the current methodologies, covers some of the recent trends or advances and tries to make predictions of the overall path membrane molecular dynamics will follow in the coming years. In general, the overview presented is geared towards the general scientific community, who may wish to introduce the use of these methodologies in light of these changes, making molecular dynamic simulations more feasible for general scientific or medical research.

A study on attempt for determination of permeation kinetics of coumarin at lipid bilayer of liposomes by using capillary electrophoresis with moment analysis theory

It was tried to develop a moment analysis method for the determination of lipid membrane permeability. The first absolute and second central moments of elution peaks measured by liposome electrokinetic chromatography (LEKC) are analyzed by using moment equations. As a concrete example, elution peak profiles of coumarin in a LEKC system, in which liposomes consisting of 1-palmitoyl-2-oleoyl-sn‑glycero-3-phosphocholine (POPC) and phosphatidylserine (PS) are used as a pseudo-stationary phase, were analyzed. It seems that lipid membrane permeability of coumarin across the lipid bilayer of POPC/PS liposomes was measured by the moment analysis method because previous permeability measurements using parallel artificial membrane permeability assay (PAMPA) and Caco-2 cells indicated that coumarin is permeable across lipid bilayer. However, it was also pointed out that the moment analysis method with LEKC is not effective for the determination of lipid membrane permeability and that it provides information about adsorption/desorption kinetics at lipid bilayer of liposomes. Therefore, different moment equations were also developed for the determination of adsorption/desorption rate constants of coumarin from the LEKC data. It was demonstrated that permeation rate constants at lipid bilayer or adsorption/desorption rate constants can be determined from the LEKC data on the basis of moment analysis theory for the mass transfer phenomena of coumarin at the lipid bilayer of POPC/PS liposomes. Mass transfer kinetics of solutes at lipid bilayer should be determined under the conditions that liposomes originally be because they are self-assembling and dynamic systems formed through weak interactions between phospholipid monomers. The moment analysis method using LEKC is effective for the experimental determination of the mass transfer rate constants at the lipid bilayer of liposomes because neither immobilization nor chemical modification of liposomes is necessary when LEKC data are measured. It is expected that the results of this study contribute to the dissemination of an opportunity for the determination of permeation rate constants or adsorption/desorption rate constants at the lipid bilayer of liposomes to many researchers because capillary electrophoresis is widespread.

Nature-Derived Compounds as Potential Bioactive Leads against CDK9-Induced Cancer: Computational and Network Pharmacology Approaches

Given the importance of cyclin-dependent kinases (CDKs) in the maintenance of cell development, gene transcription, and other essential biological operations, CDK blockers have been generated to manage a variety of disorders resulting from CDK irregularities. Furthermore, CDK9 has a crucial role in transcription by regulating short-lived anti-apoptotic genes necessary for cancer cell persistence. Addressing CDK9 with blockers has consequently emerged as a promising treatment for cancer. This study scrutinizes the effectiveness of nature-derived compounds (geniposidic acid, quercetin, geniposide, curcumin, and withanolide C) against CDK9 through computational approaches. A molecular docking study was performed after preparing the protein and the ligands. The selected blockers of the CDK9 exerted reliable binding affinities (−8.114 kcal/mol to −13.908 kcal/mol) against the selected protein, resulting in promising candidates compared to the co-crystallized ligand (LCI). The binding affinity of geniposidic acid (−13.908 kcal/mol) to CDK9 is higher than quercetin (−10.775 kcal/mol), geniposide (−9.969 kcal/mol), curcumin (−9.898 kcal/mol), withanolide C (−8.114 kcal/mol), and the co-crystallized ligand LCI (−11.425 kcal/mol). Therefore, geniposidic acid is a promising inhibitor of CDK9. Moreover, the molecular dynamics studies assessed the structure–function relationships and protein–ligand interactions. The network pharmacology study for the selected ligands demonstrated the auspicious compound–target–pathway signaling pathways vital in developing tumor, tumor cell growth, differentiation, and promoting tumor cell progression. Moreover, this study concluded by analyzing the computational approaches the natural-derived compounds that have potential interacting activities against CDK9 and, therefore, can be considered promising candidates for CKD9-induced cancer. To substantiate this study’s outcomes, in vivo research is recommended.

Molecular modeling in cardiovascular pharmacology: Current state of the art and perspectives

Molecular modeling in pharmacology is a promising emerging tool for exploring drug interactions with cellular components. Recent advances in molecular simulations, big data analysis, and artificial intelligence (AI) have opened new opportunities for rationalizing drug interactions with their pharmacological targets. Despite the obvious utility and increasing impact of computational approaches, their development is not progressing at the same speed in different fields of pharmacology. Here, we review current in silico techniques used in cardiovascular diseases (CVDs), cardiological drug discovery, and assessment of cardiotoxicity. In silico techniques are paving the way to a new era in cardiovascular medicine, but their use somewhat lags behind that in other fields.

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.

Interactions of chlorogenic acid and isochlorogenic acid A with model lipid bilayer membranes: Insights from molecular dynamics simulations

Because of the negative side-effects of synthetic preservatives, the naturally-occurring polyphenols aroused intense interest of researchers. It has been suggested that chlorogenic acid (CA) and isochlorogenic acid A (iso-CAA) were good candidates to replace the synthetic preservatives. Moreover, the bactericidal activity of iso-CAA was stronger than CA, and the anti-bacterial activities of iso-CAA and CA were highly membrane-dependent. However, the mechanisms were still unclear. Therefore, in the present study, we investigated the mechanisms of the interactions between the two polyphenols and lipid bilayers through molecular dynamics simulations. The results revealed that iso-CAA could be inserted much deeper into POPG lipid bilayer than CA. We also found that hydrophobic interactions and hydrogen bonds both contributed to the insertion of iso-CAA into the POPG lipid bilayer, and the quinic acid moiety was the key structure in iso-CAA to form hydrogen bonds with POPG lipid bilayer. We believed that these findings would provide more useful information to explain the stronger bactericidal activity of iso-CAA than CA at the atomic level.

Automated Coarse-Grained Mapping Algorithm for the Martini Force Field and Benchmarks for Membrane-Water Partitioning

With a view to high-throughput simulations, we present an automated system for mapping and parameterizing organic molecules for use with the coarse-grained Martini force field. The method scales to larger molecules and a broader chemical space than existing schemes. The core of the mapping process is a graph-based analysis of the molecule's bonding network, which has the advantages of being fast, general, and preserving symmetry. The parameterization process pays special attention to coarse-grained beads in aromatic rings. It also includes a method for building efficient and stable frameworks of constraints for molecules with structural rigidity. The performance of the method is tested on a diverse set of 87 neutral organic molecules and the ability of the resulting models to capture octanol-water and membrane-water partition coefficients. In the latter case, we introduce an adaptive method for extracting partition coefficients from free-energy profiles to take into account the interfacial region of the membrane. We also use the models to probe the response of membrane-water partitioning to the cholesterol content of the membrane.

Ligand binding at the protein-lipid interface: strategic considerations for drug design

Many drug targets are embedded within the phospholipid bilayer of cellular membranes, including G protein-coupled receptors, ion channels, transporters and membrane-bound enzymes. Increasing evidence from biophysical and structural studies suggests that many small-molecule drugs commonly associate with these targets at binding sites at the protein-phospholipid interface. Without a direct path from bulk solvent to a binding site, a drug must first partition in the phospholipid membrane before interacting with the protein target. This membrane access mechanism necessarily affects the interpretation of potency data, structure-activity relationships, pharmacokinetics and physicochemical properties for drugs that target these sites. With an increasing number of small-molecule intramembrane binding sites revealed through X-ray crystallography and cryogenic electron microscopy, we suggest that ligand-lipid interactions likely play a larger role in small-molecule drug action than commonly appreciated. This Perspective introduces key concepts and drug design considerations to aid discovery teams operating within this target space, and discusses challenges and future opportunities in the field.

The Extended Eighth-Shell method for periodic boundary conditions with rotational symmetry

The Eighth-Shell method for parallelization of molecular dynamics simulations has previously been shown to be the most optimal for efficiency at large process counts. However, in its current formulation only the P1 space group is supported for periodic boundary conditions (PBC) and thus reflection and/or rotational crystal symmetries are not supported. In this work, we outline the development and implementation of the Extended Eighth-Shell (EES) method that allows rotational symmetry by using an extended import region compared to the ES method. It simulates only the asymmetric unit and communicates coordinates and forces with images that correspond to P2₁ PBC. The P2₁ PBC has application in lipid bilayer simulations as it can be used to allow lipids to switch leaftlets, thus rapidly balancing the chemical potential difference between the two layers. Our results show that the EES method scales efficiently over large number of processes and can be used for simulations with P2₁ symmetry in an orthorhombic crystal.

Mastering the Gram-negative bacterial barrier - Chemical approaches to increase bacterial bioavailability of antibiotics

To win the battle against resistant, pathogenic bacteria, novel classes of anti-infectives and targets are urgently needed. Bacterial uptake, distribution, metabolic and efflux pathways of antibiotics in Gram-negative bacteria determine what we here refer to as bacterial bioavailability. Understanding these mechanisms from a chemical perspective is essential for anti-infective activity and hence, drug discovery as well as drug delivery. A systematic and critical discussion of in bacterio, in vitro and in silico assays reveals that a sufficiently accurate holistic approach is still missing. We expect new findings based on Gram-negative bacterial bioavailability to guide future anti-infective research.

Functional Group Distributions, Partition Coefficients, and Resistance Factors in Lipid Bilayers Using Site Identification by Ligand Competitive Saturation

Small molecules such as metabolites and drugs must pass through the membrane of the cell, a barrier primarily comprising phospholipid bilayers and embedded proteins. To better understand the process of passive diffusion, knowledge of the ability of various functional groups to partition across bilayers and the associated energetics would be of utility. In the present study, the site identification by ligand competitive saturation (SILCS) methodology has been applied to sample the distributions of a diverse set of chemical solutes representing the functional groups of small molecules across phospholipid bilayers composed of 0.9:0.1 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/cholesterol and a mixture of 0.52:0.18:0.3 1,2-dioleoyl-sn-glycero-3-phospho-l-serine/1,2-dioleoyl-sn-glycero-3-phosphocholine/cholesterol used in parallel artificial membrane permeability assay experiments. A combination of oscillating chemical potential grand canonical Monte Carlo and molecular dynamics in the SILCS simulations was applied to achieve solute sampling through the bilayers and surrounding aqueous environment from which the distribution of solutes and the functional groups they represent were obtained. Results show differential distribution of aliphatic versus aromatic groups with the former having increased sampling in the center of the bilayers versus in the region of the glycerol linker for the latter. Variations in the distribution of different polar groups are evident, with large differences between negative acetate and positive methylammonium with accumulation of the polar-neutral and acetate solutes above the bilayer head groups. Conversion of the distributions to absolute free energies allows for a detailed understanding of energetics of functional groups in different regions of the bilayers and for calculation of absolute free-energy profiles of multifunctional drug-like molecules across the bilayers from which partition coefficients and resistance factors suitable for insertion into the homogenous solubility-diffusion equation for calculation of permeability were obtained. Comparisons of the calculated bilayer/solution partition coefficients with 1-octanol/water experimental data for both drug-like molecules and the solutes show overall good agreement, validating the calculated distributions and associated absolute free-energy profiles.

Oxidative Stress Evaluation in Ischemia Reperfusion Models: Characteristics, Limits and Perspectives

Ischemia reperfusion injury is a complex process consisting of a seemingly chaotic but actually organized and compartmentalized shutdown of cell function, of which oxidative stress is a key component. Studying oxidative stress, which results in an imbalance between reactive oxygen species (ROS) production and antioxidant defense activity, is a multi-faceted issue, particularly considering the double function of ROS, assuming roles as physiological intracellular signals and as mediators of cellular component damage. Herein, we propose a comprehensive overview of the tools available to explore oxidative stress, particularly in the study of ischemia reperfusion. Applying chemistry as well as biology, we present the different models currently developed to study oxidative stress, spanning the vitro and the silico, discussing the advantages and the drawbacks of each set-up, including the issues relating to the use of in vitro hypoxia as a surrogate for ischemia. Having identified the limitations of historical models, we shall study new paradigms, including the use of stem cell-derived organoids, as a bridge between the in vitro and the in vivo comprising 3D intercellular interactions in vivo and versatile pathway investigations in vitro. We shall conclude this review by distancing ourselves from "wet" biology and reviewing the in silico, computer-based, mathematical modeling, and numerical simulation options: (a) molecular modeling with quantum chemistry and molecular dynamic algorithms, which facilitates the study of molecule-to-molecule interactions, and the integration of a compound in a dynamic environment (the plasma membrane...); (b) integrative systemic models, which can include many facets of complex mechanisms such as oxidative stress or ischemia reperfusion and help to formulate integrated predictions and to enhance understanding of dynamic interaction between pathways.

A quantitative model for metabolic intervention using gut microbes

As medicine shifts toward precision-based and personalized therapeutics, utilizing more complex biomolecules to treat increasingly difficult and rare conditions, microorganisms provide an avenue for realizing the production and processing necessary for novel drug pipelines. More so, probiotic microbes can be co-opted to deliver therapeutics by oral administration as living drugs, able to survive and safely transit the digestive tract. As living therapeutics are in their nascency, traditional pharmacokinetic-pharmacodynamic (PK-PD) models for evaluating drug candidates are not appropriate for this novel platform. Using a living therapeutic in late-stage clinical development for phenylketonuria (PKU) as a case study, we adapt traditional oral drug delivery models to properly evaluate and inform the engineering of living therapeutics. We develop the adapted for living therapeutics compartmental absorption and transit (ALT-CAT) model to provide metrics for drug efficacy across nine age groups of PKU patients and evaluate model parameters that are influenced by patient physiology, microbe selection and therapeutic production, and dosing formulations. In particular, the ALT-CAT model describes the mathematical framework to model the behavior of orally delivered engineered bacteria that act as living therapeutics by adapting similar methods that have been developed and widely-used for small molecular drug delivery and absorption.

Structure and Dynamics of Dioleoyl-Phosphatidylcholine Bilayers under the Influence of Quercetin and Rutin

Quercetin and rutin, two widely studied flavonoids with applications foreseen in the sectors of pharmaceutical and cosmetic industries, have been chosen as model compounds for a detailed structural and dynamical investigation onto their influence on fluid lipid bilayers. Combining global small angle X-ray scattering analysis with molecular dynamics, various changes in the properties of dioleoyl-phosphatidylcholine (DOPC) bilayers have been determined. The solubility of quercetin in DOPC membranes is assured up to 12 mol %, whereas rutin, with additional glucose and rhamnose groups, are fully soluble only up to 6 mol %. Both flavonoids induce an increase in membrane undulations and thin the bilayers slightly (<1 Å) in a concentration dependent manner, wherein quercetin shows a stronger effect. Concomitantly, in the order of 2-4%, the adjacent bilayer distance increases with the flavonoid's concentration. Partial molecular areas of quercetin and rutin are determined to be 26 and 51 Ų, respectively. Simulated averaged areas per molecule confirm these estimates. A 60° tilted orientation of quercetin is observed with respect to the bilayer normal, whereas the flavonoid moiety of rutin is oriented more perpendicular (α-angle 30°) to the membrane surface. Both flavonoid moieties are located at a depth of 12 and 16 Å for quercetin and rutin, respectively, while their anionic forms display a location closer to the polar interface. Finally, at both simulated concentrations (1.5 and 12 mol %), DOPC-rutin systems induce a stronger packing of the pure DOPC lipid bilayer, mainly due to stronger attractive electrostatic interactions in the polar lipid head region.

Permeation of Biopolymers Across the Cell Membrane: A Computational Comparative Study on Polylactic Acid and Polyhydroxyalkanoate

Polymeric nanoparticles, which by virtue of their size (1-1000 nm) are able to penetrate even into cells, are attracting increasing interest in the emerging field of nanomedicine, as devices for, e.g., drugs or vaccines delivery. Because of the involved dimensional scale in the nanoparticle/cell membrane interactions, modeling approaches at molecular level are the natural choice in order to understand the impact of nanoparticle formulation on cellular uptake mechanisms. In this work, the passive permeation across cell membrane of oligomers made of two employed polymers in the biomedical field [poly-D,L-lactic acid (PDLA) and poly(3-hydroxydecanoate) (P3HD)] is investigated at fundamental atomic scale through molecular dynamics simulations. The free energy profile related to membrane crossing is computed adopting umbrella sampling. Passive permeation is also investigated using a coarse-grained model with MARTINI force field, adopting well-tempered metadynamics. Simulation results showed that P3HD permeation is favored with respect to PDLA by virtue of its higher hydrophobicity. The free energy profiles obtained at full atomistic and coarse-grained scale are in good agreement each for P3HD, while only a qualitative agreement was obtained for PDLA. Results suggest that a reparameterization of non-bonded interactions of the adopted MARTINI beads for the oligomer is needed in order to obtain a better agreement with more accurate simulations at atomic scale.

Cyanine dyes with tail length asymmetry enhance photoselection: A multiscale study on DiD probes in a liquid disordered membrane

Visualization of membrane domains like lipid rafts in natural or artificial membranes is a crucial task for cell biology. For this purpose, fluorescence microscopy is often used. Since fluorescing probes in lipid membranes partition specifically in e.g. local liquid disordered or liquid ordered environments, the consequent changes in their orientation and location are both theoretically and experimentally of interest. Here we focused on a liquid disordered membrane phase and performed molecular dynamics (MD) simulations of the indocarbocyanine DiD probes by varying the length of the attached alkyl tails and also the length of the cyanine backbone. From the probed compounds in a DOPC lipid bilayer at ambient temperature, a varying orientation of the transition dipole moment was observed, which is crucial for fluorescence microscopy and which, through photoselection, was found to be surprisingly more effective for asymmetric probes than for the symmetric ones. Furthermore, we observed that the orientation of the probes was dependent on the tail length; with the methyls or propyls attached, DiD oriented with its tails facing the water, contrary to the ones with longer tails. With advanced hybrid QM/MM calculations we show that the different local environment for differently oriented probes affected the one-photon absorption spectra, that was blue-shifted for the short-tailed DiD with respect to the DiDs with longer tails. We show here that the presented probes can be successfully used for fluorescence microscopy and we believe that the described properties bring further insight for the experimental use of these probes.

Isolation of Bacteriocin-producing Staphylococcus spp. Strains from Human Skin Wounds, Soft Tissue Infections and Bovine Mastitis

A collection of 206 Staphylococcus spp. isolates was investigated for their ability to produce compounds exhibiting antistaphylococcal activity. This group included Staphylococcus aureus, Staphylococcus epidermidis, Staphylococcus xylosus strains recovered from bovine mastitis (n = 158) and human skin wounds and soft tissues infections (n = 48). Production of substances with antimicrobial activity was observed in six strains. Five of them were recovered from bovine mastitis, and one was isolated from the infected human skin wound. Three of the six antimicrobials produced by the different strains showed substantial loss of antimicrobial activity upon treatment with proteolytic enzymes, which suggests their peptidic structure. Additional studies have shown that one of the putative bacteriocins was efficiently secreted to the liquid medium, facilitating its large-scale production and isolation. The peptide produced by the M2B strain exhibited promising activity; however, against narrow spectrum of Staphylococcus spp. clinical and animal isolates. Growth inhibition was observed only in the case of 13 (including nine S. aureus, three S. xylosus and one S. epidermidis strains) out of 206 strains tested. Important advantage of the produced agent was its high thermal stability. Fifteen minutes of incubation at 90°C did not affect its antimicrobial potential. The highest efficiency of production of the agent was demonstrated in TSB medium after 24 hours at 37°C. The researches revealed that ability to production of bacteriocin among staphylococci is not very common. Only one (S. xylosus strain assigned as M2B) out of 206 strains tested produced satisfactory amounts of antistaphylococcal bacteriocin. In spite of that, we would encourage other researchers for investigation of their collections of Staphylococcus spp. isolates towards selection strains producing antimicrobial agents.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Molecular Dynamics Simulations of Membrane Permeability

This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.

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.

Orientational distribution of DPH in lipid membranes: a comparison of molecular dynamics calculations and experimental time-resolved anisotropy experiments

The behavior of the fluorescent probe diphenylhexatriene (DPH) in different lipid phases is investigated. The rotational autocorrelation functions are calculated in order to model the time-resolved fluorescence anisotropy decay. The role of the order parameters is discussed.

Atomistic Picture of Fluorescent Probes with Hydrocarbon Tails in Lipid Bilayer Membranes: An Investigation of Selective Affinities and Fluorescent Anisotropies in Different Environmental Phases

By reverting to spectroscopy, changes in the biological environment of a fluorescent probe can be monitored and the presence of various phases of the surrounding lipid bilayer membranes can be detected. However, it is currently not always clear in which phase the probe resides. The well-known orange 1,1'-dioctadecyl-3,3,3',3'-tetramethylindodicarbo-cyanine perchlorate (DiI-C18(5)) fluorophore, for instance, and the new, blue BODIPY (4,4-difluoro-4-bora-3 a,4 a-diaza- s-indacene) derivative were experimentally seen to target and highlight identical parts of giant unilamellar vesicles of various compositions, comprising mixtures of dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), sphingomyelin (SM), and cholesterol (Chol). However, it was not clear which of the coexisting membrane phases were visualized (Bacalum et al., Langmuir. 2016, 32, 3495). The present study addresses this issue by utilizing large-scale molecular dynamics simulations and the z-constraint method, which allows evaluating Gibbs free-energy profiles. The current calculations give an indication why, at room temperature, both BODIPY and DiI-C18(5) probes prefer the gel (S_o) phase in DOPC/DPPC (2:3 molar ratio) and the liquid-ordered (L_o) phase in DOPC/SM/Chol (1:2:1 molar ratio) mixtures. This study highlights the important differences in orientation and location and therefore in efficiency between the probes when they are used in fluorescence microscopy to screen various lipid bilayer membrane phases. Dependent on the lipid composition, the angle between the transition-state dipole moments of both probes and the normal to the membrane is found to deviate clearly from 90°. It is seen that the DiI-C18(5) probe is located in the headgroup region of the SM/Chol mixture, in close contact with water molecules. A fluorescence anisotropy study also indicates that DiI-C18(5) gives rise to a distinctive behavior in the SM/Chol membrane compared to the other considered membranes. The latter behavior has not been seen for the studied BODIPY probe, which is located deeper in the membrane.

Local low dose curcumin treatment improves functional recovery and remyelination in a rat model of sciatic nerve crush through inhibition of oxidative stress

Traumatic injuries to peripheral nerves are frequent, however, specific pharmacological treatments are currently lacking. Curcumin has antioxidant, anti-inflammatory and neuroprotective properties but high oral doses are required for therapeutic use, particularly due to its low bioavailability. The aim of the present study was to investigate the effects of local and continuous treatment using low curcumin doses on functional recovery and nerve regeneration after rat sciatic nerve crush (SNC). Curcumin was administered by osmotic pumps with a catheter delivering the drug at the injury site (0.2 mg/day for 4 weeks). Functionally, early improvements in mechanical sensitivity, finger spacing of the injured paw, skilful walking and grip strength were observed in curcumin-treated animals. The curcumin treatment increased expression of compact myelin proteins (MPZ and PMP22), myelin sheath thickness and, correspondingly, increased motor and sensitive nerve conduction velocity. Microscopic analysis of gastrocnemius muscle indicated a curcumin-induced decrease in neurogenic lesions. Curcumin treatment reduced the production of reactive oxygen species (ROS) (which were notably produced by macrophages), lipid peroxidation and increased expression of transcription factor Nrf2. In silico analyses indicated that curcumin combines all the characteristics required to be an efficient lipid peroxidation inhibitor at the heart of biological membranes, hence protecting their degradation due to ROS. This antioxidant capacity is likely to contribute to the beneficial effects of curcumin after SNC injury. These results demonstrate that, when administrated locally, low doses of curcumin represent a promising therapy for peripheral nerve regeneration.

Antimicrobial Peptides: Interaction With Model and Biological Membranes and Synergism With Chemical Antibiotics

Antimicrobial peptides (AMPs) are promising novel antibiotics since they have shown antimicrobial activity against a wide range of bacterial species, including multiresistant bacteria; however, toxicity is the major barrier to convert antimicrobial peptides into active drugs. A profound and proper understanding of the complex interactions between these peptides and biological membranes using biophysical tools and model membranes seems to be a key factor in the race to develop a suitable antimicrobial peptide therapy for clinical use. In the search for such therapy, different combined approaches with conventional antibiotics have been evaluated in recent years and demonstrated to improve the therapeutic potential of AMPs. Some of these approaches have revealed promising additive or synergistic activity between AMPs and chemical antibiotics. This review will give an insight into the possibilities that physicochemical tools can give in the AMPs research and also address the state of the art on the current promising combined therapies between AMPs and conventional antibiotics, which appear to be a plausible future opportunity for AMPs treatment.

Digging into Lipid Membrane Permeation for Cardiac Ion Channel Blocker d-Sotalol with All-Atom Simulations

Interactions of drug molecules with lipid membranes play crucial role in their accessibility of cellular targets and can be an important predictor of their therapeutic and safety profiles. Very little is known about spatial localization of various drugs in the lipid bilayers, their active form (ionization state) or translocation rates and therefore potency to bind to different sites in membrane proteins. All-atom molecular simulations may help to map drug partitioning kinetics and thermodynamics, thus providing in-depth assessment of drug lipophilicity. As a proof of principle, we evaluated extensively lipid membrane partitioning of d-sotalol, well-known blocker of a cardiac potassium channel K_v11.1 encoded by the hERG gene, with reported substantial proclivity for arrhythmogenesis. We developed the positively charged (cationic) and neutral d-sotalol models, compatible with the biomolecular CHARMM force field, and subjected them to all-atom molecular dynamics (MD) simulations of drug partitioning through hydrated lipid membranes, aiming to elucidate thermodynamics and kinetics of their translocation and thus putative propensities for hydrophobic and aqueous hERG access. We found that only a neutral form of d-sotalol accumulates in the membrane interior and can move across the bilayer within millisecond time scale, and can be relevant to a lipophilic channel access. The computed water-membrane partitioning coefficient for this form is in good agreement with experiment. There is a large energetic barrier for a cationic form of the drug, dominant in water, to cross the membrane, resulting in slow membrane translocation kinetics. However, this form of the drug can be important for an aqueous access pathway through the intracellular gate of hERG. This route will likely occur after a neutral form of a drug crosses the membrane and subsequently re-protonates. Our study serves to demonstrate a first step toward a framework for multi-scale in silico safety pharmacology, and identifies some of the challenges that lie therein.

Characterization of the interaction forces in a drug carrier complex of doxorubicin with a drug-binding peptide

Polypeptide-based materials are used as building blocks for drug delivery systems aimed at toxicity decrease in chemotherapeutics. A molecular-level approach is adopted for investigating the non-covalent interactions between doxorubicin and a recently synthesized drug-binging peptide as a key part of a system for delivery to neoplastic cells. Molecular dynamics simulations in aqueous solution at room and body temperature are applied to investigate the structure and the binding modes within the drug-peptide complex. The tryptophans are outlined as the main chemotherapeutic adsorption sites, and the importance of their placement in the peptide sequence is highlighted. The drug-peptide binging energy is evaluated by density functional theory calculations. Principal component analysis reveals comparable importance of several types of interaction for the binding strength. π-Stacking is dominant, but other factors are also significant: intercalation, peptide backbone stacking, electrostatics, dispersion, and solvation. Intra- and intermolecular H-bonding also stabilizes the complexes. The influence of solvent molecules on the binding energy is mild. The obtained data characterize the drug-to-peptide attachment as a mainly attractive collective process with interactions spanning a broad range of values. These results explain with atomistic detail the experimentally registered doxorubicin-binging ability of the peptide and outline the complex as a prospective carrying unit that can be employed in design of drug delivery systems.

Anesthetic Agents of Plant Origin: A Review of Phytochemicals with Anesthetic Activity

The majority of currently used anesthetic agents are derived from or associated with natural products, especially plants, as evidenced by cocaine that was isolated from coca (Erythroxylum coca, Erythroxylaceae) and became a prototype of modern local anesthetics and by thymol and eugenol contained in thyme (Thymus vulgaris, Lamiaceae) and clove (Syzygium aromaticum, Myrtaceae), respectively, both of which are structurally and mechanistically similar to intravenous phenolic anesthetics. This paper reviews different classes of phytochemicals with the anesthetic activity and their characteristic molecular structures that could be lead compounds for anesthetics and anesthesia-related drugs. Phytochemicals in research papers published between 1996 and 2016 were retrieved from the point of view of well-known modes of anesthetic action, that is, the mechanistic interactions with Na⁺ channels, γ-aminobutyric acid type A receptors, N-methyl-d-aspartate receptors and lipid membranes. The searched phytochemicals include terpenoids, alkaloids and flavonoids because they have been frequently reported to possess local anesthetic, general anesthetic, antinociceptive, analgesic or sedative property. Clinical applicability of phytochemicals to local and general anesthesia is discussed by referring to animal in vivo experiments and human pre-clinical trials. This review will give structural suggestions for novel anesthetic agents of plant origin.

Differentiating Physicochemical Properties between Addictive and Nonaddictive ADHD Drugs Revealed by Molecular Dynamics Simulation Studies

Attention-deficit/hyperactivity disorder (ADHD) is the most commonly diagnosed mental disorder of children and adolescents. Although psychostimulants are currently the first-line drugs for ADHD, their highly addictive profile raises great abuse concerns. It is known that psychostimulants' addictiveness is largely attributed to their interaction with dopamine transporter (DAT) and their binding modes in DAT can thus facilitate the understanding of the mechanism underlining drugs' addictiveness. However, no DAT residue able to discriminate ADHD drugs' addictiveness is identified, and the way how different drug structures affect their abuse liability is still elusive. In this study, multiple computational methods were integrated to differentiate binding modes between approved psychostimulants and ADHD drugs of little addictiveness. As a result, variation in energy contribution of 8 residues between addictive and nonaddictive drugs was observed, and a reduction in hydrophobicity of drugs' 2 functional groups was identified as the indicator of drugs' addictiveness. This finding agreed well with the physicochemical properties of 8 officially reported controlled substances. The identified variations in binding mode can shed light on the mechanism underlining drugs' addictiveness, which may thus facilitate the discovery of improved ADHD therapeutics with reduced addictive profile.

Models for drug absorption from the small intestine: where are we and where are we going?

The small intestine is a complex organ with movements, flora, mucus and flows. Despite this, the most widely used absorption models consider the organ a cylindrical monoepithelial tube. This review presents the recent evolution of models to take into consideration the complex nature of gut physiology. The most commonly encountered issues are ethical (in vivo models) and differences in drug transport as a result of a modified expression of drug transporters or metabolic enzymes compared with human (in vitro and in vivo models). Finally, this review discusses the way forward to reach an ideal equilibrium between reproducibility, predictability and efficiency for predicting permeability. The features of an ideal model are listed as a guideline for future development.