Binding and dynamics of melatonin at the interface of phosphatidylcholine-cholesterol membranes

Insight into molecular structures and dynamical properties of niosome bilayers containing melatonin molecules: a molecular dynamics simulation approach

Niosomes represent vesicular carriers capable of encapsulating both hydrophobic and hydrophilic drugs within their inner core or bilayer shell. They are typically composed of non-ionic synthetic surfactants such as sorbitan monostearate (Span60) with the addition of cholesterol (Chol). The physical properties and stability of niosomal vesicles strongly depend on the composition of their bilayers, which plays a significant role in determining the efficiency of drug encapsulation and release in drug delivery systems. In this study, we have explored the interactions between melatonin (Mel) molecules and the niosome bilayer, as well as their resulting physical properties. Molecular dynamics simulations were employed to investigate melatonin-inserted niosome bilayers, both with and without the inclusion of cholesterol. The simulation results revealed that cholesterol notably influences the location of melatonin molecules within the niosome bilayers. In the absence of cholesterol, melatonin tends to occupy the region around the Span60 tail groups. However, in the presence of cholesterol, melatonin is found in the vicinity of the Span60 head groups. Melatonin molecules in niosome bilayers without cholesterol exhibit a more ordered orientation when compared to those in bilayers containing 50 mol% cholesterol. The bilayer structure of the Span60/Mel and Span60/Chol/Mel systems exhibited a liquid-disordered phase (Ld). In contrast, the Span60/Chol bilayer system displays a liquid-ordered phase (Lo) with less fluidity. This study reveals that melatonin induces a disorderly bilayer structure and greater lateral expansion, whereas cholesterol induces an orderly bilayer structure and a more condensed effect. Cholesterol plays a crucial role in condensing the bilayer structure with stronger interactions between Span60 and cholesterol. The addition of 50 mol% cholesterol in the Span60 bilayers not only enhances the stability and rigidity of niosomes but also facilitates the easier release of melatonin from the bilayer membranes. This finding is particularly valuable in the context of preparing niosomes for drug delivery systems.

Computational Advances of Protein/Neurotransmitter-membrane Interactions Involved in Vesicle Fusion and Neurotransmitter Release

Vesicle fusion is of crucial importance to neuronal communication at neuron terminals. The exquisite but complex fusion machinery for neurotransmitter release is tightly controlled and regulated by protein/neurotransmitter-membrane interactions. Computational 'microscopies', in particular molecular dynamics simulations and related techniques, have provided notable insight into the physiological process over the past decades, and have made enormous contributions to fields such as neurology, pharmacology and pathophysiology. Here we review the computational advances of protein/neurotransmitter-membrane interactions related to presynaptic vesicle-membrane fusion and neurotransmitter release, and outline the in silico challenges ahead for understanding this important physiological process.

In Silico Drug Design of Benzothiadiazine Derivatives Interacting with Phospholipid Cell Membranes

The use of drugs derived from benzothiadiazine, a bicyclic heterocyclic benzene derivative, has become a widespread treatment for diseases such as hypertension, low blood sugar or the human immunodeficiency virus, among others. In this work we have investigated the interactions of benzothiadiazine and four of its derivatives designed in silico with model zwitterionic cell membranes formed by dioleoylphosphatidylcholine, 1,2-dioleoyl-sn-glycero-3-phosphoserine and cholesterol at the liquid–crystal phase inside aqueous potassium chloride solution. We have elucidated the local structure of benzothiadiazine by means of microsecond molecular dynamics simulations of systems including a benzothiadiazine molecule or one of its derivatives. Such derivatives were obtained by the substitution of a single hydrogen site of benzothiadiazine by two different classes of chemical groups, one of them electron-donating groups (methyl and ethyl) and another one by electron-accepting groups (fluorine and trifluoromethyl). Our data have revealed that benzothiadiazine derivatives have a strong affinity to stay at the cell membrane interface although their solvation characteristics can vary significantly—they can be fully solvated by water in short periods of time or continuously attached to specific lipid sites during intervals of 10–70 ns. Furthermore, benzothiadiazines are able to bind lipids and cholesterol chains by means of single and double hydrogen-bonds of characteristic lengths between 1.6 and 2.1 Å.

Predicting the conformational variability of oncogenic GTP-bound G12D mutated KRas-4B proteins at zwitterionic model cell membranes

KRas proteins are the largest family of mutated Ras isoforms, participating in a wide variety of cancers. Due to their importance, large effort is being carried out on drug development by small-molecule inhibitors. However, understanding protein conformational variability remains a challenge in drug discovery. In the case of the Ras family, their multiple conformational states can affect the binding of potential drug inhibitors. To overcome this challenge, we propose a computational framework based on combined all-atom Molecular Dynamics and Metadynamics simulations in order to accurately access conformational variants of the target protein. We tested the methodology using a G12D mutated GTP bound oncogenic KRas-4B protein located at the interface of a DOPC/DOPS/cholesterol model anionic cell membrane. Two main orientations of KRas-4B at the anionic membrane have been determined. The corresponding torsional angles are taken as reliable reaction coordinates so that free-energy landscapes are obtained by well-tempered metadynamics simulations, revealing local and global minima of the free-energy hypersurface and unveiling reactive paths of the system between the two preferential orientations. We have observed that GTP-binding to KRas-4B has huge influence on the stabilisation of the protein and it can potentially help to open Switch I/II druggable pockets, lowering energy barriers between stable states and resulting in cumulative conformers of KRas-4B. This may highlight new opportunities for targeting the unique meta-stable states through the design of new efficient drugs.

Structure of benzothiadiazine at zwitterionic phospholipid cell membranes

The use of drugs derived from benzothiadiazine, which is a bicyclic heterocyclic benzene derivative, has become a widespread treatment for diseases such as hypertension (treated with diuretics such as bendroflumethiazide or chlorothiazide), low blood sugar (treated with non-diuretic diazoxide), or the human immunodeficiency virus, among others. In this work, we have investigated the interactions of benzothiadiazine with the basic components of cell membranes and solvents, such as phospholipids, cholesterol, ions, and water. The analysis of the mutual microscopic interactions is of central importance to elucidate the local structure of benzothiadiazine as well as the mechanisms responsible for the access of benzothiadiazine to the interior of the cell. We have performed molecular dynamics simulations of benzothiadiazine embedded in three different model zwitterionic bilayer membranes made by dimyristoylphosphatidylcholine, dioleoylphosphatidylcholine, 1,2-dioleoyl-sn-glycero-3-phosphoserine, and cholesterol inside aqueous sodium-chloride solution in order to systematically examine microscopic interactions of benzothiadiazine with the cell membrane at liquid-crystalline phase conditions. From data obtained through radial distribution functions, hydrogen-bonding lengths, and potentials of mean force based on reversible work calculations, we have observed that benzothiadiazine has a strong affinity to stay at the cell membrane interface although it can be fully solvated by water in short periods of time. Furthermore, benzothiadiazine is able to bind lipids and cholesterol chains by means of single and double hydrogen-bonds of different characteristic lengths.

Melatonin: Regulation of Biomolecular Condensates in Neurodegenerative Disorders

Biomolecular condensates are membraneless organelles (MLOs) that form dynamic, chemically distinct subcellular compartments organizing macromolecules such as proteins, RNA, and DNA in unicellular prokaryotic bacteria and complex eukaryotic cells. Separated from surrounding environments, MLOs in the nucleoplasm, cytoplasm, and mitochondria assemble by liquid-liquid phase separation (LLPS) into transient, non-static, liquid-like droplets that regulate essential molecular functions. LLPS is primarily controlled by post-translational modifications (PTMs) that fine-tune the balance between attractive and repulsive charge states and/or binding motifs of proteins. Aberrant phase separation due to dysregulated membrane lipid rafts and/or PTMs, as well as the absence of adequate hydrotropic small molecules such as ATP, or the presence of specific RNA proteins can cause pathological protein aggregation in neurodegenerative disorders. Melatonin may exert a dominant influence over phase separation in biomolecular condensates by optimizing membrane and MLO interdependent reactions through stabilizing lipid raft domains, reducing line tension, and maintaining negative membrane curvature and fluidity. As a potent antioxidant, melatonin protects cardiolipin and other membrane lipids from peroxidation cascades, supporting protein trafficking, signaling, ion channel activities, and ATPase functionality during condensate coacervation or dissolution. Melatonin may even control condensate LLPS through PTM and balance mRNA- and RNA-binding protein composition by regulating N6-methyladenosine (m6A) modifications. There is currently a lack of pharmaceuticals targeting neurodegenerative disorders via the regulation of phase separation. The potential of melatonin in the modulation of biomolecular condensate in the attenuation of aberrant condensate aggregation in neurodegenerative disorders is discussed in this review.

Microscopic Interactions of Melatonin, Serotonin and Tryptophan with Zwitterionic Phospholipid Membranes

The interactions at the atomic level between small molecules and the main components of cellular plasma membranes are crucial for elucidating the mechanisms allowing for the entrance of such small species inside the cell. We have performed molecular dynamics and metadynamics simulations of tryptophan, serotonin, and melatonin at the interface of zwitterionic phospholipid bilayers. In this work, we will review recent computer simulation developments and report microscopic properties, such as the area per lipid and thickness of the membranes, atomic radial distribution functions, angular orientations, and free energy landscapes of small molecule binding to the membrane. Cholesterol affects the behaviour of the small molecules, which are mainly buried in the interfacial regions. We have observed a competition between the binding of small molecules to phospholipids and cholesterol through lipidic hydrogen-bonds. Free energy barriers that are associated to translational and orientational changes of melatonin have been found to be between 10-20 kJ/mol for distances of 1 nm between melatonin and the center of the membrane. Corresponding barriers for tryptophan and serotonin that are obtained from reversible work methods are of the order of 10 kJ/mol and reveal strong hydrogen bonding between such species and specific phospholipid sites. The diffusion of tryptophan and melatonin is of the order of 10-7 cm2/s for the cholesterol-free and cholesterol-rich setups.

Influence of Cholesterol on the Orientation of the Farnesylated GTP-Bound KRas-4B Binding with Anionic Model Membranes

The Ras family of proteins is tethered to the inner leaflet of the cell membranes which plays an essential role in signal transduction pathways that promote cellular proliferation, survival, growth, and differentiation. KRas-4B, the most mutated Ras isoform in different cancers, has been under extensive study for more than two decades. Here we have focused our interest on the influence of cholesterol on the orientations that KRas-4B adopts with respect to the plane of the anionic model membranes. How cholesterol in the bilayer might modulate preferences for specific orientation states is far from clear. Herein, after analyzing data from in total 4000 ns-long molecular dynamics (MD) simulations for four KRas-4B systems, properties such as the area per lipid and thickness of the membrane as well as selected radial distribution functions, penetration of different moieties of KRas-4B, and internal conformational fluctuations of flexible moieties in KRas-4B have been calculated. It has been shown that high cholesterol content in the plasma membrane (PM) favors one orientation state (OS1), exposing the effector-binding loop for signal transduction in the cell from the atomic level. We confirm that high cholesterol in the PM helps KRas-4B mutant stay in its constitutively active state, which suggests that high cholesterol intake can increase mortality and may promote cancer progression for cancer patients. We propose that during the treatment of KRas-4B-related cancers, reducing the cholesterol level in the PM and sustaining cancer progression by controlling the plasma cholesterol intake might be taken into account in anti-cancer therapies.

Chiral Recognition of Flexible Melatonin Receptor Ligands Induced by Conformational Equilibria

N-anilinoethylamides are a class of melatoninergic agents with the aniline portion mimicking the indole ring of the natural ligand and the ethylamide chain reproducing that of melatonin. The simplest compound in this class, N-{2-[(3-methoxyphenyl)methylamino]ethyl}acetamide (UCM793), has nanomolar binding affinity for MT₁ and MT₂ membrane receptors. To explore the effect of chain conformation on receptor binding, a methyl group was inserted on the methylene alpha or beta to the amide nitrogen and conformational equilibria were investigated by NMR spectroscopy and molecular dynamics simulations. Receptor affinity was conserved only for the beta-methyl derivative, which also showed significant stereoselectivity, with the (S) enantiomer being the eutomer. Molecular dynamics simulations, validated by NMR spectroscopy, showed that the beta-methyl group affects the conformational preferences of the ethylamide chain. Docking into the receptor crystal structure provides a rationale for the observed chiral recognition, suggesting that the (S)-beta-methyl group favors the conformation that better fits the receptor binding site.

Cellular absorption of small molecules: free energy landscapes of melatonin binding at phospholipid membranes

Free energy calculations are essential to unveil mechanisms at the atomic scale such as binding of small solutes and their translocation across cell membranes, eventually producing cellular absorption. Melatonin regulates biological rhythms and is directly related to carcinogenesis and neurodegenerative disorders. Free energy landscapes obtained from well-tempered metadynamics simulations precisely describe the characteristics of melatonin binding to specific sites in the membrane and reveal the role of cholesterol in free energy barrier crossing. A specific molecular torsional angle and the distance between melatonin and the center of the membrane along the normal to the membrane Z-axis have been considered as suitable reaction coordinates. Free energy barriers between two particular orientations of the molecular structure (folded and extended) have been found to be of about 18 kJ/mol for z-distances of about 1–2 nm. The ability of cholesterol to expel melatonin out of the internal regions of the membrane towards the interface and the external solvent is explained from a free energy perspective. The calculations reported here offer detailed free energy landscapes of melatonin embedded in model cell membranes and reveal microscopic information on its transition between free energy minima, including the location of relevant transition states, and provide clues on the role of cholesterol in the cellular absorption of small molecules.