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

Isomer-sourced structure iteration methods for in silico development of inhibitors: Inducing GTP-bound NRAS-Q61 oncogenic mutations to an "off-like" state

The NRAS-mutant subset of melanoma represent some of the most aggressive and deadliest types associated with poor overall survival. Unfortunately, for more than 40 years, no therapeutic agent directly targeting NRAS mutations has been clinically approved. In this work, based on microsecond scale molecular dynamics simulations, the effect of Q61 mutations on NRAS conformational characteristics is revealed at the atomic level. The GTP-bound NRAS-Q61R and Q61K mutations show a specific targetable pocket between Switch-II and α-helix 3 whereas the NRAS-Q61L non-polar mutation category shows a different targetable pocket. Moreover, a new isomer-sourced structure iteration method has been developed for the in silico design of potential inhibitor prototypes for oncogenes. We show the possibility of a designed prototype HM-387 to target activated NRAS-Q61R and that it can gradually induce the transition from the activated NRAS-Q61R to an "off-like" state.

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.

Structural and Dynamical Studies of a Lipid-Nanoclay Composite Layer at the Air-Water Interface

We modulate the adsorption affinities of nanoclay particles for the air-water interface by changing the cationic surface charge composition of the lipid monolayer and thereby tune the attractive electrostatic interaction between the positively charged lipid layer and the zwitterionic nanoclay particles in the water subphase. Our findings emphasize the significance of electrostatic interaction between lipids and the nanoclay, as well as its impact on the structural and viscoelastic features of the composite layer. We use surface pressure (Π)-mean molecular area (A) isotherms, atomic force microscope (AFM), Brewster angle microscopy (BAM), and energy dispersive X-ray spectrsocopy (EDXS) measurements to analyze the structure phases of lipid and lipid-nanoclay composite interfacial layer. The Π-A isotherm curve shows that the lipid-nanoclay composite layer has a larger lift-off area than the neat lipid layer, indicating that nanoparticles adsorb at the lipid layer via electrostatic interaction between lipid and nanoclay molecules. The surface density of the adsorbed nanoclay particles increases with an increase in the composition of the cationic lipid molecules. The stress relaxation response of the composite layer, measured using step compression measurements, exhibits exponential decay and ubiquitous dependence on the cationic dimyristoy-trimethylammonium propane (DMTAP) composition in the lipid layer with crossover to faster relaxation dynamics at DMTAP > 0.75. The power-law study of the frequency-dependent dynamic viscoelastic responses of the interfacial layer, measured using the barrier oscillation method, reveals a transition from glass-like response from neat lipid layer to gel-like dynamic response for the lipid-nanoclay composite layer. A solid-like behavior is evident for all the interface layers with dilation elastic modulus (E') > dilational viscous modulus (E″); however, the dynamic response of the neat layer is largely frequency-independent, whereas lipid-nanoclay composite layers with DMTAP > 0.75 reveal a frequency-dependent dynamic responses. The frequency-dependent power-law exponent of E', E″ increases on increasing the fractional composition of cationic DMTAP from 0.1 to 1.0, which forms a saturated interface of laponite particles and behaves as a viscoelastic gel in 2D.

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 Å.

Free-Energy Simulations Support a Lipophilic Binding Route for Melatonin Receptors

The effects of the neurohormone melatonin are mediated by the activation of the GPCRs MT₁ and MT₂ in a variety of tissues. Crystal structures suggest ligand access to the orthosteric binding site of MT₁ and MT₂ receptors through a lateral channel between transmembrane (TM) helices IV and V. We investigated the feasibility of this lipophilic entry route for 2-iodomelatonin, a nonselective agonist with a slower dissociation rate from the MT₂ receptor, applying enhanced sampling simulations and free-energy calculations. 2-Iodomelatonin unbinding was investigated with steered molecular dynamics simulations which revealed different trajectories passing through the gap between TM helices IV and V for both receptors. For one of these unbinding trajectories from the MT₁ receptor, an umbrella-sampling protocol with path-collective variables provided a calculated energy barrier consistent with the experimental dissociation rate. The side-chain flexibility of Tyr5.38 was significantly different in the two receptor subtypes, as assessed by metadynamics simulations, and during ligand unbinding it frequently assumes an open conformation in the MT₁ but not in the MT₂ receptor, favoring 2-iodomelatonin egress. Taken together, our simulations are consistent with the possibility that the gap between TM IV and V is a way of connecting the orthosteric binding site and the membrane core for lipophilic melatonin receptor ligands. Our simulations also suggest that the open state of Tyr5.38 generates a small pocket on the surface of MT₁ receptor, which could participate in the recognition of MT₁-selective ligands and may be exploited in the design of new selective compounds.

Energy landscapes and dynamics of ion translocation through membrane transporters: a meeting ground for physics, chemistry, and biology

The dynamics of ion translocation through membrane transporters is visualized from a comprehensive point of view by a Gibbs energy landscape approach. The ΔG calculations have been performed with the Kirkwood-Tanford-Warshel (KTW) electrostatic theory that properly takes into account the self-energies of the ions. The Gibbs energy landscapes for translocation of a single charge and an ion pair are calculated, compared, and contrasted as a function of the order parameter, and the characteristics of the frustrated system with bistability for the ion pair are described and quantified in considerable detail. These calculations have been compared with experimental data on the ΔG of ion pairs in proteins. It is shown that, under suitable conditions, the adverse Gibbs energy barrier can be almost completely compensated by the sum of the electrostatic energy of the charge-charge interactions and the solvation energy of the ion pair. The maxima in ΔG_KTW with interionic distance in the bound H⁺ - A^- charge pair on the enzyme is interpreted in thermodynamic and molecular mechanistic terms, and biological implications for molecular mechanisms of ATP synthesis are discussed. The timescale at which the order parameter moves between two stable states has been estimated by solving the dynamical equations of motion, and a wealth of novel insights into energy transduction during ATP synthesis by the membrane-bound F_OF₁-ATP synthase transporter is offered. In summary, a unifying analytical framework that integrates physics, chemistry, and biology has been developed for ion translocation by membrane transporters for the first time by means of a Gibbs energy landscape approach.

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.

Long-lasting Salt Bridges Provide the Anchoring Mechanism of Oncogenic Kirsten Rat Sarcoma Proteins at Cell Membranes

RAS proteins work as GDP-GTP binary switches and regulate cytoplasmic signaling networks that are able to control several cellular processes, playing an essential role in signal transduction pathways involved in cell growth, differentiation, and survival, so that overacting RAS signaling can lead to cancer. One of the hardest challenges to face is the design of mutation-selective therapeutic strategies. In this work, a G12D-mutated farnesylated GTP-bound Kirsten RAt sarcoma (KRAS) protein has been simulated at the interface of a DOPC/DOPS/cholesterol model anionic cell membrane. A specific long-lasting salt bridge connection between farnesyl and the hypervariable region of the protein has been identified as the main mechanism responsible for the binding of oncogenic farnesylated KRAS-4B to the cell membrane. Free-energy landscapes allowed us to characterize local and global minima of KRAS-4B binding to the cell membrane, revealing the main pathways between anchored and released states.

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.

A Role of Newly Found Auxiliary Site in Phospholipase A1 from Thai Banded Tiger Wasp (Vespa affinis) in Its Enzymatic Enhancement: In Silico Homology Modeling and Molecular Dynamics Insights

Phospholipase A₁ from Thai banded tiger wasp (Vespa affinis) venom also known as Ves a 1 plays an essential role in fatal vespid allergy. Ves a 1 becomes an important therapeutic target for toxin remedy. However, established Ves a 1 structure or a mechanism of Ves a 1 function were not well documented. This circumstance has prevented efficient design of a potential phospholipase A₁ inhibitor. In our study, we successfully recruited homology modeling and molecular dynamic (MD) simulation to model Ves a 1 three-dimensional structure. The Ves a 1 structure along with dynamic behaviors were visualized and explained. In addition, we performed molecular docking of Ves a 1 with 1,2-Dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC) lipid to assess a possible lipid binding site. Interestingly, molecular docking predicted another lipid binding region apart from its corresponding catalytic site, suggesting an auxiliary role of the alternative site at the Ves a 1 surface. The new molecular mechanism related to the surface lipid binding site (auxiliary site) provided better understanding of how phospholipase A₁ structure facilitates its enzymatic function. This auxiliary site, conserved among Hymenoptera species as well as some mammalian lipases, could be a guide for interaction-based design of a novel phospholipase A₁ inhibitor.