Interaction of MRI Contrast Agent [Gd(DOTA)]- with Lipid Membranes: A Molecular Dynamics Study.

Contrast agents are important imaging probes in clinical MRI, allowing the identification of anatomic changes that otherwise would not be possible. Intensive research on the development of new contrast agents is being made to image specific pathological markers or sense local biochemical changes. The most widely used MRI contrast agents are based on gadolinium(III) complexes. Due to their very high charge density, they have low permeability through tight biological barriers such as the blood-brain barrier, hampering their application in the diagnosis of neurological disorders. In this study, we explore the interaction between the widely used contrast agent [Gd(DOTA)]- (Dotarem) and POPC lipid bilayers by means of molecular dynamics simulations. This metal complex is a standard reference where several chemical modifications have been introduced to improve key properties such as bioavailability and targeting. The simulations unveil detailed insights into the agent's interaction with the lipid bilayer, offering perspectives beyond experimental methods. Various properties, including the impact on global and local bilayer properties, were analyzed. As expected, the results indicate a low partition coefficient (KP) and high permeation barrier for this reference compound. Nevertheless, favorable interactions are established with the membrane leading to moderately long residence times. While coordination of one inner-sphere water molecule is maintained for the membrane-associated chelate, the physical-chemical attributes of [Gd(DOTA)]- as a MRI contrast agent are affected. Namely, increases in the rotational correlation times and in the residence time of the inner-sphere water are observed, with the former expected to significantly increase the water proton relaxivity. This work establishes a reference framework for the use of simulations to guide the rational design of new contrast agents with improved relaxivity and bioavailability and for the development of liposome-based formulations for use as imaging probes or theranostic agents.

Fluorescent Probes cis- and trans-Parinaric Acids in Fluid and Gel Lipid Bilayers: A Molecular Dynamics Study.

Fluorescence probes are indispensable tools in biochemical and biophysical membrane studies. Most of them possess extrinsic fluorophores, which often constitute a source of uncertainty and potential perturbation to the host system. In this regard, the few available intrinsically fluorescent membrane probes acquire increased importance. Among them, cis- and trans-parinaric acids (c-PnA and t-PnA, respectively) stand out as probes of membrane order and dynamics. These two compounds are long-chained fatty acids, differing solely in the configurations of two double bonds of their conjugated tetraene fluorophore. In this work, we employed all-atom and coarse-grained molecular dynamics simulations to study the behavior of c-PnA and t-PnA in lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), representative of the liquid disordered and solid ordered lipid phases, respectively. All-atom simulations indicate that the two probes show similar location and orientation in the simulated systems, with the carboxylate facing the water/lipid interface and the tail spanning the membrane leaflet. The two probes establish interactions with the solvent and lipids to a similar degree in POPC. However, the almost linear t-PnA molecules have tighter lipid packing around them, especially in DPPC, where they also interact more with positively charged lipid choline groups. Probably for these reasons, while both probes show similar partition (assessed from computed free energy profiles across bilayers) to POPC, t-PnA clearly partitions more extensively than c-PnA to the gel phase. t-PnA also displays more hindered fluorophore rotation, especially in DPPC. Our results agree very well with experimental fluorescence data from the literature and allow deeper understanding of the behavior of these two reporters of membrane organization.

Interaction of Hoechst 33342 with POPC Membranes at Different pH Values.

Hoechst 33342 (H33342) is a fluorescent probe that is commonly used to stain the DNA of living cells. To do so, it needs to interact with and permeate through cell membranes, despite its high overall charge at physiological pH values. In this work, we address the effect of pH in the association of H33342 with lipid bilayers using a combined experimental and computational approach. The partition of H33342 to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid membranes was experimentally quantified using fluorescence spectroscopy and isothermal titration calorimetry (ITC) measurements. Quantum chemical calculations were performed to select the most stable isomer of H33342 for the overall charges 0, +1, and +2, expected to predominate across the 5 < pH < 10 range. The interaction of these isomers with POPC bilayers was then studied by both unrestrained and umbrella sampling molecular dynamics (MD) simulations. Both experimental results and computational free energy profiles indicate that the partition coefficient of H33342 displays a small variation over a wide pH range, not exceeding one order of magnitude. The enthalpy variation upon partition to the membrane suggests efficient hydrogen bonding between the probe and the lipid, namely, for the protonated +2 form, which was confirmed in the MD simulation studies. The relatively high lipophilicity obtained for the charged species contrasts with the decrease in their general hydrophobicity as estimated from octanol/water partition. This highlights the distinction between lipophilicity and hydrophobicity, as well as the importance of considering the association with lipid bilayers when predicting the affinity for biomembranes.

Modeling Gd3+ Complexes for Molecular Dynamics Simulations: Toward a Rational Optimization of MRI Contrast Agents.

The correct parametrization of lanthanide complexes is of the utmost importance for their characterization using computational tools such as molecular dynamics simulations. This allows the optimization of their properties for a wide range of applications, including medical imaging. Here we present a systematic study to establish the best strategies for the correct parametrization of lanthanide complexes using [Gd(DOTA)]- as a reference, which is used as a contrast agent in MRI. We chose the bonded model to parametrize the lanthanide complexes, which is especially important when considering the study of the complex as a whole (e.g., for the study of the dynamics of its interaction with proteins or membranes). We followed two strategies: a so-called heuristic approach employing strategies already published by other authors and another based on the more recent MCPB.py tool. Adjustment of the Lennard-Jones parameters of the metal was required. The final topologies obtained with both strategies were able to reproduce the experimental ion to oxygen distance, vibrational frequencies, and other structural properties. We report a new strategy to adjust the Lennard-Jones parameters of the metal ion in order to capture dynamic properties such as the residence time of the capping water (τm). For the first time, the correct assessment of the τm value for Gd-based complexes was possible by recording the dissociative events over up to 10 µs all-atom simulations. The MCPB.py tool allowed the accurate parametrization of [Gd(DOTA)]- in a simpler procedure, and in this case, the dynamics of the water molecules in the outer hydration sphere was also characterized. This sphere was divided into the first hydration layer, an intermediate region, and an outer hydration layer, with a residence time of 18, 10 and 19 ps, respectively, independent of the nonbonded parameters chosen for Gd3+. The Lennard-Jones parameters of Gd3+ obtained here for [Gd(DOTA)]- may be used with similarly structured gadolinium MRI contrast agents. This allows the use of molecular dynamics simulations to characterize and optimize the contrast agent properties. The characterization of their interaction with membranes and proteins will permit the design of new targeted contrast agents with improved pharmacokinetics.

Molecular Dynamics Simulations: Advances and Applications.

Molecular dynamics (MD) simulations have led to great advances in many scientific disciplines, such as chemical physics, materials science, and biophysics [...].

Interactions between Rhodamine Dyes and Model Membrane Systems-Insights from Molecular Dynamics Simulations.

BACKGROUND: rhodamines are dyes widely used as fluorescent tags in cell imaging, probing of mitochondrial membrane potential, and as P-glycoprotein model substrates. In all these applications, detailed understanding of the interaction between rhodamines and biomembranes is fundamental. METHODS: we combined atomistic molecular dynamics (MD) simulations and fluorescence spectroscopy to characterize the interaction between rhodamines 123 and B (Rh123 and RhB, respectively) and POPC bilayers. RESULTS: while the xanthene moiety orients roughly parallel to the membrane plane in unrestrained MD simulations, variations on the relative position of the benzoic ring (below the xanthene for Rh123, above it for RhB) were observed, and related to the structure of the two dyes and their interactions with water and lipids. Subtle distinctions were found among different ionization forms of the probes. Experimentally, RhB displayed a lipid/water partition coefficient more than two orders of magnitude higher than Rh123, in agreement with free energy profiles obtained from umbrella sampling MD. CONCLUSIONS: this work provided detailed insights on the similarities and differences in the behavior of bilayer-inserted Rh123 and RhB, related to the structure of the probes. The much higher affinity of RhB for the membranes increases the local concentration and explains its higher apparent affinity for P-glycoprotein reconstituted in model membranes.

Dancing with oils - the interaction of lipases with different oil/water interfaces.

The use of enzymes as biocatalysts in industrial applications has received much attention during the last few years. Lipases are widely employed in the food and cosmetic industry, for the synthesis of novel biomaterials and as a greener solution for the treatment of waste cooking oils (WCO). The latter topic has been widely explored with the use of enzymes from several origins and types, for the treatment of different used and non-used cooking oils. The experimental conditions of such works are also quite broad, hampering the detailed understanding of the process. In this work we present a detailed characterization of the interaction of several commonly used lipases with different types of vegetal oils and food fats through coarse-grained molecular dynamics simulations. First, the molecular details of the oil/water (O/W) mixtures, namely at the O/W interface, are described. The O/W interface was found to be enriched in triglyceride molecules with higher polarity. Then, the interaction of lipases with oil mixtures is characterized from different perspectives, including the identification of the most important protein residues for this process. The lipases from Thermomyces lanuginosus (TLL), Rhizomucor miehei (RML) and Candida antarctica (CALB) were found to bind to the O/W interface in a manner that makes the protein binding site more available for the oil molecules. These enzymes were also found to efficiently bind to the O/W interface of all oil mixtures, which in addition to reactivity factors, may explain the efficient applicability of these enzymes to a large variety of edible oils and WCO.

The Secret Lives of Fluorescent Membrane Probes as Revealed by Molecular Dynamics Simulations.

Fluorescent probes have been employed for more than half a century to study the structure and dynamics of model and biological membranes, using spectroscopic and/or microscopic experimental approaches. While their utilization has led to tremendous progress in our knowledge of membrane biophysics and physiology, in some respects the behavior of bilayer-inserted membrane probes has long remained inscrutable. The location, orientation and interaction of fluorophores with lipid and/or water molecules are often not well known, and they are crucial for understanding what the probe is actually reporting. Moreover, because the probe is an extraneous inclusion, it may perturb the properties of the host membrane system, altering the very properties it is supposed to measure. For these reasons, the need for independent methodologies to assess the behavior of bilayer-inserted fluorescence probes has been recognized for a long time. Because of recent improvements in computational tools, molecular dynamics (MD) simulations have become a popular means of obtaining this important information. The present review addresses MD studies of all major classes of fluorescent membrane probes, focusing in the period between 2011 and 2020, during which such work has undergone a dramatic surge in both the number of studies and the variety of probes and properties accessed.

Orientation of nitro-group governs the fluorescence lifetime of nitrobenzoxadiazole (NBD)-labeled lipids in lipid bilayers.

Nitrobenzoxadiazole (NBD) labeled lipids are popular fluorescent probes of membrane structure and dynamics, and have been widely used in both model systems and living cells. Irrespective of attachment to the lipid head group or hydrocarbon chains, the NBD fluorophore generally adopts a transverse bilayer location near the host lipid carbonyl/glycerol moieties. Still, considerable variability is observed in the measured fluorescence lifetimes, indicating that overall fluorophore location is not the determinant of NBD fluorescence properties. Combining fluorescence experiments and molecular dynamics simulations, we show that for two almost identical NBD probes, significant differences in fluorophore orientation and fluorescence lifetime are observed. Integrating these findings with literature data, we demonstrate a correlation between NBD orientation and fluorescence lifetime. The latter is longer when the NBD nitro group is predominantly oriented towards the bilayer interior, compared to probes for which it points to the water medium.

Fluorescence of nitrobenzoxadiazole (NBD)-labeled lipids in model membranes is connected not to lipid mobility but to probe location.

Nitrobenzoxadiazole (NBD)-labeled lipids are popular fluorescent membrane probes. However, the understanding of important aspects of the photophysics of NBD remains incomplete, including the observed shift in the emission spectrum of NBD-lipids to longer wavelengths following excitation at the red edge of the absorption spectrum (red-edge excitation shift or REES). REES of NBD-lipids in membrane environments has been previously interpreted as reflecting restricted mobility of solvent surrounding the fluorophore. However, this requires a large change in the dipole moment (Δµ) of NBD upon excitation. Previous calculations of the value of Δµ of NBD in the literature have been carried out using outdated semi-empirical methods, leading to conflicting values. Using up-to-date density functional theory methods, we recalculated the value of Δµ and verified that it is rather small (∼2 D). Fluorescence measurements confirmed that the value of REES is ∼16 nm for 1,2-dioleoyl-sn-glycero-3-phospho-l-serine-N-(NBD) (NBD-PS) in dioleoylphosphatidylcholine vesicles. However, the observed shift is independent of both the temperature and the presence of cholesterol and is therefore insensitive to the mobility and hydration of the membrane. Moreover, red-edge excitation leads to an increased contribution of the decay component with a shorter lifetime, whereas time-resolved emission spectra of NBD-PS displayed an atypical blue shift following excitation. This excludes restrictions to solvent relaxation as the cause of the measured REES and TRES of NBD, pointing instead to the heterogeneous transverse location of probes as the origin of these effects. The latter hypothesis was confirmed by molecular dynamics simulations, from which the calculated heterogeneity of the hydration and location of NBD correlated with the measured fluorescence lifetimes/REES. Globally, our combination of theoretical and experiment-based techniques has led to a considerably improved understanding of the photophysics of NBD and a reinterpretation of its REES in particular.

Behaviour of NBD-head group labelled phosphatidylethanolamines in POPC bilayers: a molecular dynamics study.

A complete homologous series of fluorescent phosphatidylethanolamines (diCnPE), labelled at the head group with a 7-nitrobenz-2-oxa-1,3-diazo-4-yl(NBD) fluorophore and inserted in 1-palmitoyl, 2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers, was studied using atomistic molecular dynamics simulations. The longer-chained derivatives of NBD-diCnPE, with n = 14, 16, and 18, are commercially available, and widely used as fluorescent membrane probes. Properties such as location of atomic groups and acyl chain order parameters of both POPC and NBD-diCnPE, fluorophore orientation and hydrogen bonding, membrane electrostatic potential and lateral diffusion were calculated for all derivatives in the series. Most of these probes induce local disordering of POPC acyl chains, which is on the whole counterbalanced by ordering resulting from binding of sodium ions to lipid carbonyl/glycerol oxygen atoms. An exception is found for NBD-diC16PE, which displays optimal matching with POPC acyl chain length and induces a slight local ordering of phospholipid acyl chains. Compared to previously studied fatty amines, acyl chain-labelled phosphatidylcholines, and sterols bearing the same fluorescent tag, the chromophore in NBD-diCnPE locates in a similar region of the membrane (near the glycerol backbone/carbonyl region) but adopts a different orientation (with the NO2 group facing the interior of the bilayer). This modification leads to an inverted orientation of the P-N axis in the labelled lipid, which affects the interface properties, such as the membrane electrostatic potential and hydrogen bonding to lipid head group atoms. The implications of this study for the interpretation of the photophysical properties of NBD-diCnPE (complex fluorescence emission kinetics, differences with other NBD lipid probes) are discussed.

Interaction of NBD-labelled fatty amines with liquid-ordered membranes: a combined molecular dynamics simulation and fluorescence spectroscopy study.

A complete homologous series of fluorescent 7-nitrobenz-2-oxa-1,3-diazol-4-yl-(NBD) labelled fatty amines of varying alkyl chain lengths, NBD-Cn, inserted in 1-palmitoyl, 2-oleoyl-sn-glycero-3-phosphocholine (POPC) or N-palmitoyl sphingomyelin (SpM) bilayers, with 50 mol% and 40 mol% cholesterol (Chol), respectively, was studied using atomistic molecular dynamics simulations. For all amphiphiles in both bilayers, the NBD fluorophore locates at the interface, in a more external position than that previously observed for pure POPC bilayers. This shallower location of the NBD group agrees with the lower fluorescent quantum yield, shorter fluorescence lifetime, and higher ionisation constants (smaller pKa) determined experimentally. The more external location is also consistent with the changes measured in steady-state fluorescence anisotropy from POPC to POPC/Chol (1 : 1) vesicles. Accordingly, the equilibrium location of the NBD group within the various bilayers is mainly dictated by bilayer compositions, and is mostly unaffected by the length of the attached alkyl chain. Similarly to the behaviour observed in POPC bilayers, the longer-chained NBD-Cn amphiphiles show significant mass density near the mixed bilayers' midplanes, and the alkyl chains of the longer derivatives, mainly NBD-C16, penetrate the opposite bilayer leaflet to some extent. However, this effect is quantitatively less pronounced in these ordered bilayers than in POPC. Similarly to POPC bilayers, the effects of these amphiphiles on the structure and dynamics of the host lipid were found to be relatively mild, in comparison with acyl-chain phospholipid analogues.

Homeostasis of free cholesterol in the blood: a preliminary evaluation and modeling of its passive transport.

The rate of noncatalyzed transfer of cholesterol (Chol) among lipoproteins and cells in the blood is of fundamental importance as a baseline to assess the role of active transport mechanisms, but remains unknown. Here we address this gap by characterizing the associa-tion of the Chol analog, ergosta-5,7,9(11),22-tetraen-3ß-ol (DHE), with the lipoproteins VLDL, LDL, HDL2, and HDL3 Combining these results with data for the association of DHE with liposomes, we elaborated a kinetic model for the noncatalyzed exchange of free Chol among blood compartments. The computational results are in good agreement with experimental values. The small deviations are explained by the nonequilibrium distribution of unesterified Chol in vivo, due to esterification and entry of new unesterified Chol, and eventual effects introduced by incubations at low temperatures. The kinetic profile of the homeostasis of unesterified Chol in the blood predicted by the model developed in this work is in good agreement with the observations in vivo, highlighting the importance of passive processes.

Permeation of a Homologous Series of NBD-Labeled Fatty Amines through Lipid Bilayers: A Molecular Dynamics Study.

Permeation through biomembranes is ubiquitous for drugs to reach their active sites. Asymmetry of the cell plasma membrane (PM) has been described as having an important role in this process. Here we describe the interaction of a homologous series of 7-nitrobenz-2-oxa-1,3-diazol-4-yl (NBD)-labeled amphiphiles (NBD-Cn, n = 4 to 16) with lipid bilayers of different compositions (1-palmitoyl, 2-oleoyl-sn-glycero-3-phosphocholine (POPC):cholesterol (1:1) and palmitoylated sphingomyelin (SpM):cholesterol (6:4)), including an asymmetric bilayer. Both unrestrained and umbrella sampling (US) simulations (at varying distances to the bilayer center) were carried out. The free energy profile of NBD-Cn at different depths in the membrane was obtained from the US simulations. The behavior of the amphiphiles during the permeation process was described regarding their orientation, chain elongation, and H-bonding to lipid and water molecules. Permeability coefficients were also calculated for the different amphiphiles of the series, using the inhomogeneous solubility-diffusion model (ISDM). Quantitative agreement with values obtained from kinetic modeling of the permeation process could not be obtained. However, for the longer, and more hydrophobic amphiphiles, the variation trend along the homologous series was qualitatively better matched by the ISDM when the equilibrium location of each amphiphile was taken as reference (ΔG = 0), compared to the usual choice of bulk water.

Interaction of Near-Infrared (NIR)-Light Responsive Probes with Lipid Membranes: A Combined Simulation and Experimental Study.

Cancer is considered a major societal challenge for the next decade worldwide. Developing strategies for simultaneous diagnosis and treatment has been considered a promising tool for fighting cancer. For this, the development of nanomaterials incorporating prototypic near-infrared (NIR)-light responsive probes, such as heptamethine cyanines, has been showing very promising results. The heptamethine cyanine-incorporating nanomaterials can be used for a tumor's visualization and, upon interaction with NIR light, can also produce a photothermal/photodynamic effect with a high spatio-temporal resolution and minimal side effects, leading to an improved therapeutic outcome. In this work, we studied the interaction of 12 NIR-light responsive probes with lipid membrane models by molecular dynamics simulations. We performed a detailed characterization of the location, orientation, and local perturbation effects of these molecules on the lipid bilayer. Based on this information, the probes were divided into two groups, predicting a lower and higher perturbation of the lipid bilayer. From each group, one molecule was selected for testing in a membrane leakage assay. The experimental data validate the hypothesis that molecules with charged substituents, which function as two polar anchors for the aqueous phase while spanning the membrane thickness, are more likely to disturb the membrane by the formation of defects and pores, increasing the membrane leakage. The obtained results are expected to contribute to the selection of the most suitable molecules for the desired application or eventually guiding the design of probe modifications for achieving an optimal interaction with tumor cell membranes.

Animal-derived products in science and current alternatives.

More than fifty years after the 3Rs definition and despite the continuous implementation of regulatory measures, animals continue to be widely used in basic research. Their use comprises not only in vivo experiments with animal models, but also the production of a variety of supplements and products of animal origin for cell and tissue culture, cell-based assays, and therapeutics. The animal-derived products most used in basic research are fetal bovine serum (FBS), extracellular matrix proteins such as Matrigel™, and antibodies. However, their production raises several ethical issues regarding animal welfare. Additionally, their biological origin is associated with a high risk of contamination, resulting, frequently, in poor scientific data for clinical translation. These issues support the search for new animal-free products able to replace FBS, Matrigel™, and antibodies in basic research. In addition, in silico methodologies play an important role in the reduction of animal use in research by refining the data previously to in vitro and in vivo experiments. In this review, we depicted the current available animal-free alternatives in in vitro research.

Interaction of a Homologous Series of Amphiphiles with P-glycoprotein in a Membrane Environment-Contributions of Polar and Non-Polar Interactions.

The transport of drugs by efflux transporters in biomembranes limits their bioavailability and is a major determinant of drug resistance development by cancer cells and pathogens. A large number of chemically dissimilar drugs are transported, and despite extensive studies, the molecular determinants of substrate specificity are still not well understood. In this work, we explore the role of polar and non-polar interactions on the interaction of a homologous series of fluorescent amphiphiles with the efflux transporter P-glycoprotein. The interaction of the amphiphiles with P-glycoprotein is evaluated through effects on ATPase activity, efficiency in inhibition of [125I]-IAAP binding, and partition to the whole native membranes containing the transporter. The results were complemented with partition to model membranes with a representative lipid composition, and details on the interactions established were obtained from MD simulations. We show that when the total concentration of amphiphile is considered, the binding parameters obtained are apparent and do not reflect the affinity for P-gp. A new formalism is proposed that includes sequestration of the amphiphiles in the lipid bilayer and the possible binding of several molecules in P-gp's substrate-binding pocket. The intrinsic binding affinity thus obtained is essentially independent of amphiphile hydrophobicity, highlighting the importance of polar interactions. An increase in the lipophilicity and amphiphilicity led to a more efficient association with the lipid bilayer, which maintains the non-polar groups of the amphiphiles in the bilayer, while the polar groups interact with P-gp's binding pocket. The presence of several amphiphiles in this orientation is proposed as a mechanism for inhibition of P-pg function.

Antiviral and antibacterial activity of hand sanitizer and surface disinfectant formulations.

Emergent diseases caused by viral and bacterial infections have proven to be a current and future challenge. The occurrence of these diseases is usually accompanied by the lack of vaccines and dedicated therapies leaving prevention as the best strategy to adopt. In that context, and apart from confinement and physical distancing measures, an increase in hygiene actions, namely hand and surface cleaning and disinfection can reduce the infection spread originated from our day-to-day routines. However, during crisis situations the high disinfectants demand can very likely lead to having them running out of stock. This impels many individuals and companies to produce their own disinfectants. Here, we explore the main components of a disinfection formulation, both for hand-rub and surface cleaning. Alcohol and non-alcohol based formulations are described, including the possibility to fine tune the properties of the final product in order to increase public acceptance while maintaining product efficacy. The action mechanisms of the main active principles are also described conjugating information from experimental and theoretical data. Overall, the main aspects to develop a disinfectant formulation are addressed, as well as their function, helping formulation developers to better understand the impact of their choices.

Effect of Protein Flexibility from Coarse-Grained Elastic Network Parameterizations on the Calculation of Free Energy Profiles of Ligand Binding.

The characterization of the affinity and binding mechanism of specific molecules to a protein active site is scientifically and industrially relevant for many applications. In principle, this information can be obtained using molecular dynamics (MD) simulations by calculating the free energy profile of the process. However, this is a computationally demanding calculation. Currently, coarse-grained (CG) force fields are very well implemented for MD simulations of biomolecular systems. These computationally efficient force fields are a major advantage to the study of large model systems and/or those requiring long simulation times. The Martini model is currently one of the most popular CG force fields for these systems. For the specific case of protein simulations, to correctly maintain the macromolecular three-dimensional structure, the Martini model needs to include an elastic network (EN). In this work, the effect of protein flexibility, as induced by three EN models compatible with the Martini force field, was tested on the calculation of free energy profiles for protein-ligand binding. The EN models used were ElNeDyn, GoMartini, and GEN. The binding of triolein (TOG) and triacetin (TAG) to a lipase protein (thermomyces lanuginosa lipase-TLL) was used as a case study. The results show that inclusion of greater flexibility in the CG parameterization of proteins is of high importance in the calculation of the free energy profiles of protein-ligand systems. However, care must be taken in order to avoid unjustified large protein deformations. In addition, due to molecular flexibility there may be no absolute need for the center of the ligand to reach the center of the protein-binding site. The calculation of the energy profile to a distance of about 0.5 nm from the active site center can be sufficient to differentiate the affinity of different ligands to a protein.

Quercetin dual interaction at the membrane level.

Quercetin impact on lipid bilayers suggests a dual action mechanism at cell membranes. This widespread polyphenol displays high partition with low interference in the more fluid membrane domains, more vulnerable to oxidative attack, but strong perturbation of cholesterol/sphingolipid enriched domains, where signalling platforms are expected to assemble.