A novel fatty acid analogue triggers CD36-GPR120 interaction and exerts anti-inflammatory action in endotoxemia.

Inflammation is a mediator of a number of chronic pathologies. We synthesized the diethyl (9Z,12Z)-octadeca-9,12-dien-1-ylphosphonate, called NKS3, which decreased lipopolysaccharide (LPS)-induced mRNA upregulation of proinflammatory cytokines (IL-1ß, IL-6 and TNF-α) not only in primary intraperitoneal and lung alveolar macrophages, but also in freshly isolated mice lung slices. The in-silico studies suggested that NKS3, being CD36 agonist, will bind to GPR120. Co-immunoprecipitation and proximity ligation assays demonstrated that NKS3 induced protein-protein interaction of CD36 with GPR120in RAW 264.7 macrophage cell line. Furthermore, NKS3, via GPR120, decreased LPS-induced activation of TAB1/TAK1/JNK pathway and the LPS-induced mRNA expression of inflammatory markers in RAW 264.7 cells. In the acute lung injury model, NKS3 decreased lung fibrosis and inflammatory cytokines (IL-1ß, IL-6 and TNF-α) and nitric oxide (NO) production in broncho-alveolar lavage fluid. NKS3 exerted a protective effect on LPS-induced remodeling of kidney and liver, and reduced circulating IL-1ß, IL-6 and TNF-α concentrations. In a septic shock model, NKS3 gavage decreased significantly the LPS-induced mortality in mice. In the last, NKS3 decreased neuroinflammation in diet-induced obese mice. Altogether, these results suggest that NKS3 is a novel anti-inflammatory agent that could be used, in the future, for the treatment of inflammation-associated pathologies.

Adsorption of water in Na-LTA zeolites: an ab initio molecular dynamics investigation.

The very wide range of applications of LTA zeolites, including the storage of tritiated water, implies that a detailed and accurate atomic-scale description of the adsorption processes taking place in their structure is crucial. To unravel with an unprecedented accuracy the mechanisms behind the water filling in NaA, we have conducted a systematic ab initio molecular dynamics investigation. Two LTA structural models, the conventional Z4A and the reduced one ZK4, have been used for static and dynamic ab initio calculations, respectively. After assessing this reduced model with comparative static DFT calculations, we start the filling of the α and ß cages by water, molecule by molecule. This allowed us to thoroughly study the interaction of water molecules with the zeolite structure and between water molecules, progressively forming H-bond chains and ring patterns as the cage is being filled. The adsorption energies could then be calculated with an unprecedented accuracy, which showed that the interaction of the molecules with the zeolite weakens as their number increases. By these methods, we have been able to highlight the primary role of Na+ cations in the interaction of water with zeolite, and inversely, the role of water in the displacement of cations when it is sufficiently solvated, allowing the passage between the α and ß cages. This phenomenon is possible thanks to the inhomogeneous distribution of water molecules on the cationic sites, as shown by our AIMD simulations, which allows the formation of water clusters. These results are important because they help in understanding how the coverage of cationic sites by water will affect the adsorption of other molecules inside the Na-LTA zeolite.

Influence of Antiplatelet Agents on the Lipid Composition of Platelet Plasma Membrane: A Lipidomics Approach with Ticagrelor and Its Active Metabolite.

Lipids contained in the plasma membrane of platelets play an important role in platelet function. Modifications in the lipid composition can fluidify or rigidify the environment around embedded receptors, in order to facilitate the access of the receptor by the drug. However, data concerning the lipid composition of platelet plasma membrane need to be updated. In addition, data on the impact of drugs on plasma membrane composition, in particular antiplatelet agents, remain sparse. After isolation of platelet plasma membrane, we assessed, using lipidomics, the effect of ticagrelor, a P2Y12 antagonist, and its active metabolite on the lipid composition of these plasma membranes. We describe the exact lipid composition of plasma membrane, including all sub-species. Ticagrelor and its active metabolite significantly increased cholesterol and phosphatidylcholine ether with short saturated acyl chains 16:0/16:0, and decreased phosphatidylcholine, suggesting overall rigidification of the membrane. Furthermore, ticagrelor and its active metabolite decreased some arachidonylated plasmalogens, suggesting a decrease in availability of arachidonic acid from the membrane phospholipids for synthesis of biologically active mediators. To conclude, ticagrelor and its active metabolite seem to influence the lipid environment of receptors embedded in the lipid bilayer and modify the behavior of the plasma membrane.

Low-Density Lipoproteins and Human Serum Albumin as Carriers of Squalenoylated Drugs: Insights from Molecular Simulations.

We have studied the interaction of three clinically promising squalenoylated drugs (gemcitabine-squalene, adenine-squalene, and doxorubicin-squalene) with low-density lipoproteins (LDL) by means of atomistic molecular dynamics simulations. It is shown that all studied squalenoylated drugs accumulate inside the LDL particles. This effect is promoted by the squalene moiety, which acts as an anchor and drives the hydrophilic drugs into the hydrophobic core of the LDL lipid droplet. Our data suggest that LDL particles could be a universal carriers of squalenoylated drugs in the bloodstream. Interaction of gemcitabine-squalene with human serum albumin (HSA) was also studied by ensemble of docking simulations. It is shown that HSA could also act as a passive carrier of this bioconjugate. It should be noted that the binding of squalene moiety to HSA was unspecific and did not occur in the binding pockets devoted to fatty acids.

Selective Inhibition of STAT3 with Respect to STAT1: Insights from Molecular Dynamics and Ensemble Docking Simulations.

STAT3 protein, which is known to be involved in cancer development, is a promising target for anticancer therapy. Successful inhibitors of STAT3 should not affect an activity of closely related protein STAT1, which makes their development challenging. The mechanisms of selectivity of several existing STAT3 inhibitors are not clear. In this work, we studied molecular mechanisms of selectivity of 13 experimentally tested STAT3 inhibitors by means of extensive molecular dynamics and ensemble docking simulations. It is shown that all studied inhibitors bind to the large part of the protein surface in an unspecific statistical manner. The binding to the dimerization interface of the SH2 domain, which is usually considered as the main target region, is not energetically preferable. Binding in this region is remarkably similar for STAT1 and STAT3 proteins and cannot explain experimentally observed selectivity toward STAT3. We propose a new mechanism of selectivity called "selectivity by distraction" for existing STAT3 inhibitors. This mechanism is based on equilibrium statistical partitioning of inhibitor molecules between protein domains. The unspecific binding of inhibitors to the DNA-binding and the coil-coil domains is stronger in STAT1 in comparison to STAT3 while the energies of their binding to SH2 domains are comparable. This "distracts" inhibitor molecules from the SH2 domain of STAT1 and leads to higher effective concentration of inhibitors in the vicinity of the SH2 domain of STAT3.

Interaction of C60 fullerenes with asymmetric and curved lipid membranes: a molecular dynamics study.

Interaction of fullerenes with asymmetric and curved DOPC/DOPS bicelles is studied by means of coarse-grained molecular dynamics simulations. The effects caused by asymmetric lipid composition of the membrane leaflets and the curvature of the membrane are analyzed. It is shown that the aggregates of fullerenes prefer to penetrate into the membrane in the regions of the moderately positive mean curvature. Upon penetration into the hydrophobic core of the membrane fullerenes avoid the regions of the extreme positive or the negative curvature. Fullerenes increase the ordering of lipid tails, which are in direct contact with them, but do not influence other lipids significantly. Our data suggest that the effects of the membrane curvature should be taken into account in the studies concerning permeability of the membranes to fullerenes and fullerene-based drug delivery systems.

Theoretical insights into the mechanism of redox switch in heat shock protein Hsp33.

Heat shock protein 33 (Hsp33) is activated in the presence of H2O2 by a very interesting redox switch based on a tetra-coordinated zinc-cysteine complex present in the fully reduced and inactive protein form. The oxidation of this zinc center by H2O2 induces formation of two S-S bridges and the zinc release followed by the protein unfolding. We report here a theoretical study of the step-by-step sequence of the overall process starting with the oxidation of the first cysteine residue and ending with the zinc release. Each reaction step is characterized by its Gibbs free energy barrier (∆G (‡)). It is predicted that the first reaction step consists in the oxidation of Cys263 by H2O2 which is by far the most reactive cysteine (∆G (‡) = 15.4 kcal mol(-1)). The next two reaction steps are the formation of the first S-S bridge between Cys263 and Cys266 (∆G (‡) = 13.6 kcal mol(-1)) and the oxidation of Cys231 by H2O2 (∆G (‡) = 20.4 kcal mol(-1)). It is then shown that the formation of the second S-S bridge (Cys231-Cys233) before the zinc release is most unlikely (∆G (‡) = 34.8 kcal mol(-1)). Instead, the release of zinc just after the oxidation of the third cysteine (Cys231) is shown to be thermodynamically (dissociation Gibbs free energy ∆G d = 6.0 kcal mol(-1)) and kinetically (reaction rate constant k d ≈ 10(6) s(-1)) favored. This result is in good agreement with the experimental data on the oxidation mechanism of Hsp33 zinc center available to date.

A molecular dynamics study of catestatin docked on nicotinic acetylcholine receptors to identify amino acids potentially involved in the binding of chromogranin A fragments.

Catestatin, a cationic and hydrophobic 21-amino acid fragment of chromogranin A, is known to be a non-competitive nicotinic antagonist acting through nicotinic acetylcholine receptors (nAChRs) to inhibit catecholamine release. Since this receptor is the target of several neuronal and non-neuronal disorder prophylaxes and treatments, this study aims at the elucidation of the binding of human catestatin to the entire nAChR reconstructed in lipid bilayers by means of docking followed by full atomistic molecular dynamics simulations. The obtained results show that the minimum free energy for the binding of the peptide and the receptor attains minimal values for locations at the pore site and in the outer beta subunit. This result is consistent with previous studies showing that catestatin occludes the pore opening. A new finding is an additional even stronger binding seat at the beta subunit and that membrane presence could be an important factor. Specific amino acids involved in catestatin binding have been identified, indicating targets for point mutation studies. In addition to improving the understanding of the interaction between the peptide and muscle-type and even other nAChR subtypes, the results of this study provide directions for future peptidomimetic research.

A principal component analysis of the dynamics of subdomains and binding sites in human serum albumin.

The conformational dynamics of human serum albumin (HSA) was investigated by principal component analysis (PCA) applied to three molecular dynamics trajectories of 200 ns each. The overlap of the essential subspaces spanned by the first 10 principal components (PC) of different trajectories was about 0.3 showing that the PCA based on a trajectory length of 200 ns is not completely convergent for this protein. The contributions of the relative motion of subdomains and of the subdomains (internal) distortion to the first 10 PCs were found to be comparable. Based on the distribution of the first 3 PC, 10 protein conformers are identified showing relative root mean square deviations (RMSD) between 2.3 and 4.6 Å. The main PCs are found to be delocalized over the whole protein structure indicating that the motions of different protein subdomains are coupled. This coupling is considered as being related to the allosteric effects observed upon ligand binding to HSA. On the other hand, the first PC of one of the three trajectories describes a conformational transition of the protein domain I that is close to that experimentally observed upon myristate binding. This is a theoretical support for the older hypothesis stating that changes of the protein onformation favorable to binding can precede the ligand complexation. A detailed all atoms PCA performed on the primary Sites 1 and 2 confirms the multiconformational character of the HSA binding sites as well as the significant coupling of their motions.

Efficiency of the monofunctionalized C60 fullerenes as membrane targeting agents studied by all-atom molecular dynamics simulations.

Transmembrane translocation of C60 fullerenes functionalized by the single amino-derivative in neutral and charged forms was studies by extensive all-atom molecular dynamics simulations. It is shown that these complexes exhibit very strong affinity to the membrane core, but their spontaneous translocation through the membrane is not possible at practical time scale. In contrast, free amino derivatives translocate through the membrane much easier than their complexes with fullerenes, but do not have pronounced affinity to the membrane interior. Our results suggest that monofunctionalized C60 could be extremely efficient membrane targeting agents, which facilitate accumulation of the water-soluble compounds in the hydrophobic core of lipid bilayer.

Ticagrelor increases its own potency at the P2Y12 receptor by directly changing the plasma membrane lipid order in platelets.

BACKGROUND AND PURPOSE: Although the amphiphilic nature of the widely used antithrombotic drug Ticagrelor is well known, it was never considered as a membranotropic agent capable of interacting with the lipid bilayer in a receptor-independent way. In this study, we investigated the influence of Ticagrelor on plasma membrane lipid order in platelets and if this modulates the potency of Ticagrelor at the P2Y12 receptor. EXPERIMENTAL APPROACH: We combined fluorescent in situ, in vitro and in silico approaches to probe the interactions between the plasma membrane of platelets and Ticagrelor. The influence of Ticagrelor on the lipid order of the platelet plasma membrane and large unilamellar vesicles was studied using the advanced fluorescent probe NR12S. Furthermore, the properties of model lipid bilayers in the presence of Ticagrelor were characterized by molecular dynamics simulations. Finally, the influence of an increased lipid order on the dose-response of platelets to Ticagrelor was studied. KEY RESULTS: Ticagrelor incorporates spontaneously into lipid bilayers and affects the lipid order of the membranes of model vesicles and isolated platelets, in a nontrivial composition and concentration-dependent manner. We showed that higher plasma membrane lipid order in platelets leads to a lower IC50 value for Ticagrelor. It is shown that membrane incorporation of Ticagrelor increases its potency at the P2Y12 receptor, by increasing the order of the platelet plasma membrane. CONCLUSION AND IMPLICATIONS: A novel dual mechanism of Ticagrelor action is suggested that combines direct binding to P2Y12 receptor with simultaneous modulation of receptor-lipid microenvironment.

A colorimetric assay adapted to fragment screening revealing aurones and chalcones as new arginase inhibitors.

Arginase, a difficult-to-target metalloenzyme, is implicated in a wide range of diseases, including cancer, infectious, and cardiovascular diseases. Despite the medical need, existing inhibitors have limited structural diversity, consisting predominantly of amino acids and their derivatives. The search for innovative arginase inhibitors has now extended to screening approaches. Due to the small and narrow active site of arginase, screening must meet the criteria of fragment-based screening. However, the limited binding capacity of fragments requires working at high concentrations, which increases the risk of interference and false positives. In this study, we investigated three colorimetric assays and selected one based on interference for screening under these challenging conditions. The subsequent adaptation and application to the screening a library of metal chelator fragments resulted in the identification of four compounds with moderate activity. The synthesis and evaluation of a series of compounds from one of the hits led to compound 21a with an IC50 value of 91.1 µM close to the reference compound piceatannol. Finally, molecular modelling supports the potential binding of aurones and chalcones to the active site of arginase, suggesting them as new candidates for the development of novel arginase inhibitors.

Oxidation reactivity of zinc-cysteine clusters in metallothionein.

Evaluating the reactivity of the metal-thiolate clusters in metallothionein (MT) is a key step in understanding the biological functions of this protein. The effects of the metal clustering and protein environment on the thiolate reactivity with hydrogen peroxide (H(2)O(2)) were investigated by performing quantum theory calculations with chemical accuracy at two levels of complexity. At the first level, the reactivity with H(2)O(2) of a model system ([(Zn)(3)(MeS)(9)](3-), MeS is methanethiolate) of the ß domain cluster of MT was evaluated using density functional theory (DFT) with the mPW1PW91 functional. At the second level of complexity, the protein environment was included in the reactant system and the calculations were performed with the hybrid ONIOM method combining the DFT-mPW1PW91 and the semiempirical PM6 levels of theory. In these conditions, the energy barrier for the oxidation of the most reactive terminal thiolate was 21.5 kcal mol(-1). This is 3 kcal mol(-1) higher than that calculated for the terminal thiolate in the model system [(Zn)(3)(MeS)(9)](3-) and about 7 kcal mol(-1) higher than that obtained for the free thiolate. In spite of this rise of the energy barrier induced by the protein environment, the thiolate oxidation by H(2)O(2) is confirmed as a possible way for metal release from MT. On the other hand, the results suggest that the antioxidant role of MT in the living cell cannot be as important as that of glutathione (which bears a free thiol).

Will C-Laurdan dethrone Laurdan in fluorescent solvent relaxation techniques for lipid membrane studies?

Development of fluorescence methods involves the necessity of understanding the fluorescent probes behavior in their ground and excited states. In the case of biological membranes, the position of the dyes in the lipid bilayer and their response after excitation are necessary parameters to correctly analyze the experiments. In the present work, we focus on two fluorescent markers, Laurdan (6-lauroyl-2-(N,N-dimethylamino)naphthalene) and its derivative C-Laurdan (6-dodecanoyl-2-[N-methyl-N-(carboxymethyl)amino]naphthalene), recently proposed for lipid raft visualization [Kim, H. M.; et al. ChemBioChem 2007, 8, 553]. C-Laurdan, by the presence of an additional carboxyl group, has an advantage over Laurdan since it has a higher sensitivity to the membrane polarity at the lipid headgroup region and a higher water solubility. This theoretical study, based on quantum mechanical (QM) and molecular dynamics (MD) simulations in a fully hydrated lipid membrane model, compare the equilibrium and dynamic properties of both probes taking into account their fluorescence lifetimes. C-Laurdan is found to be more stable than Laurdan in the headgroup region of the membrane and also much more aligned with the lipids. This study suggests that, besides the lipid raft imaging, the C-Laurdan marker can considerably extend the capabilities of fluorescent solvent relaxation method.

Enhanced potassium selectivity in a bioinspired solid nanopore.

Biological ion channels present unique ionic properties. They can be highly permeable to ions, while selecting only one type of ions, without external energy supply. An important research field has been developed to transfer these properties to solid state nanoporous membranes in order to develop artificial biomimetic nanofilters. One of the promising ways to develop biomimetic structures is based on the direct insertion of the gramicidin A, i.e. an ionic channel, inside a nanopore. Experiments have recently proved the feasibility of such a hybrid membrane with very interesting results regarding the ionic selectivity. Here, we propose to interpret these experiments using theoretical molecular dynamic simulations which allow us to analyze more profoundly the structures of the proteins confined inside the nanopore and the relation between their conformation and the observed ionic properties.

Plasma membrane lipid bilayer is druggable: Selective delivery of gemcitabine-squalene nano-medicine to cancer cells.

Up to now the lipid bilayers were rarely considered as targets in cancer therapy despite pronounced differences in lipid composition between plasma membranes of benign and malignant cells. In this study we demonstrate that the lipid bilayer of the plasma membrane is druggable and suitable for facilitating selective delivery of amphiphilic gemcitabine-squalene nanomedicines to cancer cells. Data from radioactive assays, fluorescent membrane probes and molecular dynamics simulations provide evidence of selective accumulation of gemcitabine-squalene in the plasma membranes with disrupted lipid asymmetry and its subsequent preferential uptake by malignant cells. This causes pronounced cytotoxicity on cancer cells in comparison to their benign counterparts originating from the same tissue.

Novel Fat Taste Receptor Agonists Curtail Progressive Weight Gain in Obese Male Mice.

BACKGROUND & AIMS: The spontaneous preference for dietary lipids is principally regulated by 2 lingual fat taste receptors, CD36 and GPR120. Obese animals and most of human subjects exhibit low orosensory perception of dietary fat because of malfunctioning of these taste receptors. Our aim was to target the 2 fat taste receptors by newly synthesized high affinity fatty acid agonists to decrease fat-rich food intake and obesity. METHODS: We synthesized 2 fat taste receptor agonists (FTA), NKS-3 (CD36 agonist) and NKS-5 (CD36 and GPR120 agonist). We determined their molecular dynamic interactions with fat taste receptors and the effect on Ca2+ signaling in mouse and human taste bud cells (TBC). In C57Bl/6 male mice, we assessed their gustatory perception and effects of their lingual application on activation of tongue-gut loop. We elucidated their effects on obesity and its related parameters in male mice fed a high-fat diet. RESULTS: The two FTA, NKS-3 and NKS-5, triggered higher Ca2+ signaling than a dietary long-chain fatty acid in human and mouse TBC. Mice exhibited a gustatory attraction for these compounds. In conscious mice, the application of FTA onto the tongue papillae induced activation of tongue-gut loop, marked by the release of pancreato-bile juice into collecting duct and cholecystokinin and peptide YY into blood stream. Daily intake of NKS-3 or NKS-5 via feeding bottles decreased food intake and progressive weight gain in obese mice but not in control mice. CONCLUSIONS: Our results show that targeting fat sensors in the tongue by novel chemical fat taste agonists might represent a new strategy to reduce obesity.

New bioinspired membrane made of a biological ion channel confined into the cylindrical nanopore of a solid-state polymer.

A hybrid nanoporous membrane made of a solid-state polymeric thin film in which an ion channel is confined is realized. The primary and extremely encouraging results obtained by confocal fluorescence spectroscopy and ion diffusion measurement demonstrate respectively that (i) the considered ion channel, that is, Gramicidin-A, can be confined selectively inside the nanopores and (ii) the ionic permeability of the membrane is enhanced. Atomistic molecular simulations are also reported and fruitfully compared to the experimental findings.

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.

Numerical studies of the membrane fluorescent dyes dynamics in ground and excited states.

Fluorescence methods are widely used in studies of biological and model membranes. The dynamics of membrane fluorescent markers in their ground and excited electronic states and correlations with their molecular surrounding within the fully hydrated phospholipid bilayer are still not well understood. In the present work, Quantum Mechanical (QM) calculations and Molecular Dynamics (MD) simulations are used to characterize location and interactions of two membrane polarity probes (Prodan; 6-propionyl-2-dimethylaminonaphthalene and its derivative Laurdan; 2-dimethylamino-6-lauroylnaphthalene) with the dioleoylphosphatidylcholine (DOPC) lipid bilayer model. MD simulations with fluorophores in ground and excited states are found to be a useful tool to analyze the fluorescent dye dynamics and their immediate vicinity. The results of QM calculations and MD simulations are in excellent agreement with available experimental data. The calculation shows that the two amphiphilic dyes initially placed in bulk water diffuse within 10 ns towards their final location in the lipid bilayer. Analysis of solvent relaxation process in the aqueous phase occurs on the picoseconds timescale whereas it takes nanoseconds at the lipid/water interface. Four different relaxation time constants, corresponding to different relaxation processes, where observed when the dyes were embedded into the membrane.