In silico Description of LAT1 Transport Mechanism at an Atomistic Level

Atomistic description of the OCTN1 recognition mechanism via in silico methods

The Organic Cation Transporter Novel 1 (OCTN1), also known as SLC22A4, is widely expressed in various human tissues, and involved in numerous physiological and pathological processes remains. It facilitates the transport of organic cations, zwitterions, with selectivity for positively charged solutes. Ergothioneine, an antioxidant compound, and acetylcholine (Ach) are among its substrates. Given the lack of experimentally solved structures of this protein, this study aimed at generating a reliable 3D model of OCTN1 to shed light on its substrate-binding preferences and the role of sodium in substrate recognition and transport. A chimeric model was built by grafting the large extracellular loop 1 (EL1) from an AlphaFold-generated model onto a homology model. Molecular dynamics simulations revealed domain-specific mobility, with EL1 exhibiting the highest impact on overall stability. Molecular docking simulations identified cytarabine and verapamil as highest affinity ligands, consistent with their known inhibitory effects on OCTN1. Furthermore, MM/GBSA analysis allowed the categorization of substrates into weak, good, and strong binders, with molecular weight strongly correlating with binding affinity to the recognition site. Key recognition residues, including Tyr211, Glu381, and Arg469, were identified through interaction analysis. Ach demonstrated a low interaction energy, supporting the hypothesis of its one-directional transport towards to outside of the membrane. Regarding the role of sodium, our model suggested the involvement of Glu381 in sodium binding. Molecular dynamics simulations of systems at increasing levels of Na+ concentrations revealed increased sodium occupancy around Glu381, supporting experimental data associating Na+ concentration to molecule transport. In conclusion, this study provides valuable insights into the 3D structure of OCTN1, its substrate-binding preferences, and the role of sodium in the recognition. These findings contribute to the understanding of OCTN1 involvement in various physiological and pathological processes and may have implications for drug development and disease management.

In silico investigation on structure-function relationship of members belonging to the human SLC52 transporter family

Riboflavin is an essential water-soluble vitamin that needs to be provided through the diet because of the conversion into flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN), important cofactors in hundreds of flavoenzymes. The adsorption and distribution of riboflavin is mediated by transmembrane transporters of the SLC52 family, namely RFVT1-3, whose mutations are mainly associated with two diseases, MADD and the Brown-Vialetto-Van Laere syndrome. Interest in RFVTs as pharmacological targets has increased in the last few years due to their overexpression in several cancer cells, which can be exploited both by blocking the uptake of riboflavin into the cancerous cells, and by performing cancer targeted delivery of drugs with a high affinity for RFVTs. In this work, we propose three-dimensional structural models for all three human riboflavin transporters obtained by state-of-the-art artificial intelligence-based methods, which were then further refined with molecular dynamics simulations. Furthermore, two of the most notable mutations concerning RFVT2 and RFVT3 (W31S and N21S, respectively) were investigated studying the interactions between the wild-type and mutated transporters with riboflavin.

Molecular simulations of fluoxetine in hydrated lipid bilayers, as well as in aqueous solutions containing β-cyclodextrin

Fluoxetine, which is a well-known antidepressant drug, is studied in hydrated cholesterol-free and cholesterol-containing lipid bilayers through unbiased and biased atomistic molecular dynamics simulations. The latter are conducted for the calculation of the potential of mean force (PMF) of fluoxetine along an axis perpendicular to the two leaflets of the bilayer. The PMF indicates that the drug prefers to reside inside the lipid phase and allows us to calculate important thermodynamic properties, such as the Gibbs energy difference of partitioning from the water to the lipid phase and the Gibbs energy barrier for hopping events between the two leaflets of the bilayer. The results from the biased simulations are in accord with the mass density profiles calculated from the unbiased simulations. Moreover, we estimate the effect of fluoxetine mole fraction on the order parameters of the lipid alkyl chains and on the area per lipid. It is also found that fluoxetine forms a hydrogen bond network with lipids and water molecules penetrating into the lipid phase. In addition, fluoxoetine is studied in detail in aqueous solutions containing β-cyclodextrin. It is observed from unbiased molecular dynamics simulations that the two aforementioned molecules form a noncovalent complex spontaneously and the calculated binding free energy is in agreement with the literature.

LAT1, a novel pharmacological target for the treatment of glioblastoma

The L-Type Amino Acid transporter, LAT1 (SLC7A5), has a crucial role in mediating amino acid uptake into the cells, thus modulating cell growth and proliferation as well as other intracellular functions. Different studies have reported a central role of LAT1 in glioblastoma development and progression, suggesting that the modulation of its activity could be a novel therapeutic strategy. LAT1 also has an important role in the peripheral immune system, by regulating the activation status of several immune cells through modulation of the mechanistic target of rapamycin kinase. In glioblastoma (GBM), the blood-brain barrier is disrupted, which allows the recruitment of peripheral immune cells to the tumour site. These cells, together with resident microglia, contribute to cancer growth and progression. Currently, little is known about the function of LAT1 in the reprogramming of the immune component of the tumour microenvironment in the context of GBM. In this article, we review the available data on the role of LAT1 in the regulation of GBM biology, including its potential role in the tumour microenvironment, particularly in infiltrating-peripheral immune cells and resident microglial cells. In addition, we review the available data on the main pharmacological inhibitors of LAT1, aiming to evaluate their possible role as novel therapeutics for GBM.

Comparison of Experimental Strategies to Study l-Type Amino Acid Transporter 1 (LAT1) Utilization by Ligands

l-Type amino acid transporter 1 (LAT1), expressed abundantly in the brain and placenta and overexpressed in several cancer cell types, has gained a lot of interest in drug research and development, as it can be utilized for brain-targeted drug delivery, as well as inhibiting the essential amino acid supply to cancer cells. The structure of LAT1 is today very well-known and the interactions of ligands at the binding site of LAT1 can be modeled and explained. However, less is known of LAT1′s life cycle within the cells. Moreover, the functionality of LAT1 can be measured by several different methods, which may vary between the laboratories and make the comparison of the results challenging. In the present study, the usefulness of indirect cis-inhibition methods and direct cellular uptake methods and their variations to interpret the interactions of LAT1-ligands were evaluated. Moreover, this study also highlights the importance of understanding the intracellular kinetics of LAT1-ligands, and how they can affect the regular function of LAT1 in critical tissues, such as the brain. Hence, it is discussed herein how the selected methodology influences the outcome and created knowledge of LAT1-utilizing compounds.

Photo-Switchable Sulfonylureas Binding to ATP-Sensitive Potassium Channel Reveal the Mechanism of Light-Controlled Insulin Release

ATP-sensitive potassium (KATP) channels are present in numerous organs, including the heart, brain, and pancreas. Physiological opening and closing of KATPs present in pancreatic β-cells, in response to changes in the ATP/ADP concentration ratio, are correlated with insulin release into the bloodstream. Sulfonylurea drugs, commonly used in type 2 diabetes mellitus treatment, bind to the octamer KATP channels composed of four pore-forming Kir6.2 and four SUR1 subunits and increase the probability of insulin release. Azobenzene-based derivatives of sulfonylureas, such as JB253 inspired by well-established antidiabetic drug glimepiride, allow for control of this process by light. The mechanism of that phenomenon was not known until now. In this paper, we use molecular docking, molecular dynamics, and metadynamics to reveal structural determinants explaining light-controlled insulin release. We show that both trans- and cis-JB253 bind to the same SUR1 cavity as antidiabetic sulfonylurea glibenclamide (GBM). Simulations indicate that, in contrast to trans-JB253, the cis-JB253 structure generated by blue light absorption promotes open structures of SUR1, in close similarity to the GBM effect. We postulate that in the open SUR1 structures, the N-terminal tail from Kir6.2 protruding into the SUR1 pocket is stabilized by flexible enough sulfonylureas. Therefore, the adjacent Kir6.2 pore is more often closed, which in turn facilitates insulin release. Thus, KATP conductance is regulated by peptide linkers between its Kir6.2 and SUR1 subunits, a phenomenon present in other biological signaling pathways. Our data explain the observed light-modulated activity of photoactive sulfonylureas and widen a way to develop new antidiabetic drugs having reduced adverse effects.

In silico study of levodopa in hydrated lipid bilayers at the atomistic level

This article presents atomistic molecular dynamics and umbrella sampling simulations of levodopa at various concentrations in hydrated cholesterol-free 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and cholesterol-containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. Levodopa is the standard medication for Parkinson's disease and is marketed under various trade names; in the context of this article, the levodopa molecule is mostly studied in its zwitterionic form but some results concerning the neutral levodopa are presented as well for comparison purposes. The motivation is to study in detail how levodopa behaves in different hydrated lipid membranes, primarily from the thermodynamic point of view, and reveal aspects of mechanism of its permeation through them. Dependencies of properties on the levodopa concentration are also investigated. Special attention is paid to the calculation of mass density profiles, order parameters and self-diffusion coefficients. Levodopa zwitterions, which form a hydrogen bond network with water and phospholipid molecules, are found to be preferentially located at the water/lipid interface, as well as in the aqueous phase surrounding the cholesterol-free and cholesterol-containing bilayers. This is concluded from the potentials of mean force calculated by umbrella sampling simulations as levodopa is transferred from the lipid to the aqueous phase along an axis perpendicular to the two leaflets of the membranes.

Molecular characteristics supporting l-Type amino acid transporter 1 (LAT1)-mediated translocation

l-Type amino acid transporter 1 (LAT1) is an interesting protein due to its peculiar expression profile. It can be utilized not only as a carrier for improved or targeted drug delivery, e.g., into the brain but also as a target protein by which amino acid supply can be restricted, e.g., from the cancer cells. The recognition and binding processes of LAT1-ligands, such as amino acids and clinically used small molecules, including l-dopa, gabapentin, and melphalan, are today well-known. Binding to LAT1 is crucial, particularly when designing the LAT1-inhibitors. However, it will not guarantee effective translocation across the cell membrane via LAT1, which is a definite requirement for LAT1-substrates, such as drugs that elicit their pharmacological effects inside the cells. Therefore, in the present study, the accumulation of known LAT1-utilizing compounds into the selected LAT1-expressing cancer cells (MCF-7) was explored experimentally over a time period. The differences found among the transport efficiency and affinity of the studied compounds for LAT1 were subsequently explained by docking the ligands into the human LAT1 model (based on the recent cryo-electron microscopy structure). Thus, the findings of this study clarify the favorable structural requirements of the size, shape, and polarity of the ligands that support the translocation and effective transport across the cell membrane via LAT1. This knowledge can be applied in future drug design to attain improved or targeted drug delivery and hence, successful LAT1-utilizing drugs with increased therapeutic effects.

ATP modulates SLC7A5 (LAT1) synergistically with cholesterol

The plasma membrane transporter hLAT1 is responsible for providing cells with essential amino acids. hLAT1 is over-expressed in virtually all human cancers making the protein a hot-spot in the fields of cancer and pharmacology research. However, regulatory aspects of hLAT1 biology are still poorly understood. A remarkable stimulation of transport activity was observed in the presence of physiological levels of cholesterol together with a selective increase of the affinity for the substrate on the internal site, suggesting a stabilization of the inward open conformation of hLAT1. A synergistic effect by ATP was also observed only in the presence of cholesterol. The same phenomenon was detected with the native protein. Altogether, the biochemical assays suggested that cholesterol and ATP binding sites are close to each other. The computational analysis identified two neighboring regions, one hydrophobic and one hydrophilic, to which cholesterol and ATP were docked, respectively. The computational data predicted interaction of the ϒ-phosphate of ATP with Lys 204, which was confirmed by site-directed mutagenesis. The hLAT1-K204Q mutant showed an impaired function and response to ATP. Interestingly, this residue is conserved in several members of the SLC7 family.

Demystifying the Molecular Basis of Pyrazoloquinolinones Recognition at the Extracellular α1+/β3- Interface of the GABAA Receptor by Molecular Modeling

GABA_A receptors are pentameric ligand-gated ion channels that serve as major inhibitory neurotransmitter receptors in the mammalian brain and the target of numerous clinically relevant drugs interacting with different ligand binding sites. Here, we report an in silico approach to investigate the binding of pyrazoloquinolinones (PQs) that mediate allosteric effects through the extracellular α+/β- interface of GABA_A receptors. First, we docked a potent prototype of PQs into the α1+/β3- site of a homology model of the human α1β3γ2 subtype of the GABA_A receptor. Next, for each docking pose, we computationally derived protein-ligand complexes for 18 PQ analogs with known experimental potency. Subsequently, binding energy was calculated for all complexes using the molecular mechanics-generalized Born surface area method. Finally, docking poses were quantitatively assessed in the light of experimental data to derive a binding hypothesis. Collectively, the results indicate that PQs at the α1+/β3- site likely exhibit a common binding mode that can be characterized by a hydrogen bond interaction with β3Q64 and hydrophobic interactions involving residues α1F99, β3Y62, β3M115, α1Y159, and α1Y209. Importantly, our results are in good agreement with the recently resolved cryo-Electron Microscopy structures of the human α1β3γ2 and α1β2γ2 subtypes of GABA_A receptors.

SLC6A14, a Pivotal Actor on Cancer Stage: When Function Meets Structure

SLC6A14 (ATB0,+) is a sodium- and chloride-dependent neutral and dibasic amino acid transporter that regulates the distribution of amino acids across cell membranes. The transporter is overexpressed in many human cancers characterized by an increased demand for amino acids; as such, it was recently acknowledged as a novel target for cancer therapy. The knowledge on the molecular mechanism of SLC6A14 transport is still limited, but some elegant studies on related transporters report the involvement of the 12 transmembrane α-helices in the transport mechanism, and describe structural rearrangements mediated by electrostatic interactions with some pivotal gating residues.

In the present work, we constructed a SLC6A14 model in outward-facing conformation via homology modeling and used molecular dynamics simulations to predict amino acid residues critical for substrate recognition and translocation. We docked the proteinogenic amino acids and other known substrates in the SLC6A14 binding site to study both gating regions and the exposed residues involved in transport. Interestingly, some of these residues correspond to those previously identified in other LeuT-fold transporters; however, we could also identify a novel relevant residue with such function.

For the first time, by combined approaches of molecular docking and molecular dynamics simulations, we highlight the potential role of these residues in neutral amino acid transport. This novel information unravels new aspects of the human SLC6A14 structure–function relationship and may have important outcomes for cancer treatment through the design of novel inhibitors of SLC6A14-mediated transport.