Rigorous sampling of docking poses unveils binding hypothesis for the halogenated ligands of L-type Amino acid Transporter 1 (LAT1)

Insight into telomere regulation: road to discovery and intervention in plasma drug-protein targets

BACKGROUND

Telomere length is a critical metric linked to aging, health, and disease. Currently, the exploration of target proteins related to telomere length is usually limited to the context of aging and specific diseases, which limits the discovery of more relevant drug targets. This study integrated large-scale plasma cis-pQTLs data and telomere length GWAS datasets. We used Mendelian randomization(MR) to identify drug target proteins for telomere length, providing essential clues for future precision therapy and targeted drug development.

METHODS

Using plasma cis-pQTLs data from a previous GWAS study (3,606 Pqtls associated with 2,656 proteins) and a GWAS dataset of telomere length (sample size: 472,174; GWAS ID: ieu-b-4879) from UK Biobank, using MR, external validation, and reverse causality testing, we identified essential drug target proteins for telomere length. We also performed co-localization, Phenome-wide association studies and enrichment analysis, protein-protein interaction network construction, search for existing intervening drugs, and potential drug/compound prediction for these critical targets to strengthen and expand our findings.

RESULTS

After Bonferron correction (p < 0.05/734), RPN1 (OR: 0.96; 95%CI: (0.95, 0.97)), GDI2 (OR: 0.94; 95%CI: (0.92, 0.96)), NT5C (OR: 0.97; 95%CI: (0.95, 0.98)) had a significant negative causal association with telomere length; TYRO3 (OR: 1.11; 95%CI: (1.09, 1.15)) had a significant positive causal association with telomere length. GDI2 shared the same genetic variants with telomere length (coloc.abf-PPH 4 > 0.8).

CONCLUSION

Genetically determined plasma RPN1, GDI2, NT5C, and TYRO3 have significant causal effects on telomere length and can potentially be drug targets. Further exploration of the role and mechanism of these proteins/genes in regulating telomere length is needed.

Design strategies and recent development of bioactive modulators for glutamine transporters

Glutamine transporters are integral to the metabolism of glutamine in both healthy tissues and cancerous cells, playing a pivotal role in maintaining amino acid balance, synthesizing biomolecules, and regulating redox equilibrium. Their critical functions in cellular metabolism make them promising targets for oncological therapies. Recent years have witnessed substantial progress in the field of glutamine transporters, marked by breakthroughs in understanding of their protein structures and the discovery of novel inhibitors, prodrugs, and radiotracers. This review provides a comprehensive update on the latest advancements in modulators targeting the glutamine transporter, with special attention given to LAT1 and ASCT2. It also discusses innovative approaches in drug design aimed at these transporters.

Descriptor-free QSAR: effectiveness in screening for putative inhibitors of FGFR1

The long short-term memory (LSTM) algorithm has provided solutions to the limitations of the descriptors-utilizing QSAR models in drug design. However, the direct application of LSTM remains scarce. The effectiveness of a descriptor-free QSAR (LSTM-SM) in modeling the FGFR1 inhibitors dataset while comparing with two conventional QSAR using descriptors (126 bits Morgan fingerprint and 2 D descriptors respectively) as a baseline model was investigated in this study. The validated descriptor-free QSAR model was thereafter used to screen for active FGFR1 inhibitors in the ChemDiv database and subjected to molecular docking, induced-fit docking, QM-MM optimization, and molecular dynamics simulations to filter for compounds with high binding affinity and suggest the putative mechanism of inhibition and specificity. The LSTM-SM model performed better than conventional QSAR; having accuracy, specificity, and sensitivity of 0.92, model loss of 0.025, and AUC of 0.95. Fifteen thousand compounds were predicted as actives from the ChemDiv database and four compounds were finally selected. Of the four, two showed putatively effective binding interactions with key active site residues. Molecular dynamics simulations on these compounds in complex with the receptor further give insight into the conformational dynamics of each compound bounded to the receptor. The complexes formed are stable and exhibit a similar degree of compactness. Our findings predicted the advent of self-feature extracting machine learning algorithms of compounds, and have provided the possibility of better predictive model quality that is not necessarily limited by compound descriptors. The putative FGFR1 inhibitors, with their mechanism of inhibition and specificity, were elucidated using this approachCommunicated by Ramaswamy H. Sarma.

In Vitro, In Vivo, and In Silico Analyses of Molecular Anti-Pigmentation Mechanisms of Selected Thai Rejuvenating Remedy and Bioactive Metabolites

Thai rejuvenating remedies are mixed herbal formulas promoting longevity. Due to the complexity, the biological activities of these remedies are minimal. Therefore, in this study, the authors evaluated the anti-pigmentation effect at the molecular level of the selected Thai rejuvenating remedy to fulfill the knowledge gap. First, the authors found that the selected remedy showed promising activity against the tyrosinase enzyme with an IC₅₀ value of 9.41 µg/mL. In the comparison, kojic acid (positive control) exhibited an IC₅₀ value of 3.92 µg/mL against the same enzyme. Later, the authors identified glabridin as a bioactive molecule against tyrosinase with an IC₅₀ value of 0.08 µg/mL. However, ethyl p-methoxycinnamate was the most abundant metabolite found in the remedy. The authors also found that the selected remedy and glabridin reduced the melanin content in the cell-based assay (B16F1) but not in the zebrafish larvae experiment. Finally, the authors conducted a computational investigation through molecular docking proposing a theoretical molecular interplay between glabridin, ethyl p-methoxycinnamate, and target proteins (tyrosinase and melanocortin-1 receptor, MC1R). Hence, in this study, the authors reported the molecular anti-pigmentation mechanism of the selected Thai rejuvenating remedy for the first time by combining the results from in silico, in vitro, and in vivo experiments.

Docking-Based Prediction of Peptide Binding to MHC Proteins

Major histocompatibility complex (MHC) proteins are the most polymorphic and polygenic proteins in humans. They bind peptides, derived from cleavage of different pathogenic antigens, and are responsible for presenting them to T cells. The peptides recognized by the T cell receptors are denoted as epitopes and they trigger an immune response.In this chapter, we describe a docking protocol for predicting the peptide binding to a given MHC protein using the software tool GOLD. The protocol starts with the construction of a combinatorial peptide library used in the docking and ends with the derivation of a quantitative matrix (QM) accounting for the contribution of each amino acid at each peptide position.

Accuracy and precision of binding free energy prediction for a tacrine related lead inhibitor of acetylcholinesterase with an arsenal of supercomputerized molecular modelling methods: a comparative study

Nowadays, advanced computational chemistry methods offer various strategies for revealing prospective hit structures in drug development essentially through accurate binding free energy predictions. After the era of molecular docking and quantitative structure-activity relationships, much interest has been lately oriented to perturbed molecular dynamic approaches like replica exchange with solute tempering and free energy perturbation (REST/FEP) and the potential of the mean force with adaptive biasing and accelerated weight histograms (PMF/AWH). Both of these receptor-based techniques can exploit exascale CPU&GPU supercomputers to achieve high throughput performance. In this fundamental study, we have compared the predictive power of a panel of supercomputerized molecular modelling methods to distinguish the major binding modes and the corresponding binding free energies of a promising tacrine related potential antialzheimerics in human acetylcholinesterase. The binding free energies were estimated using flexible molecular docking, molecular mechanics/generalized Born surface area/Poisson-Boltzmann surface area (MM/GBSA/PBSA), transmutation REST/FEP with 12 x 5 ns/λ windows, annihilation FEP with 20 x 5 ns/λ steps, PMF with weight histogram analysis method (WHAM) and 40 x 5 ns samples, and PMF/AWH with 10 x 100 ns replicas. Confrontation of the classical approaches such as canonical molecular dynamics and molecular docking with alchemical calculations and steered molecular dynamics enabled us to show how large errors in ΔG predictions can be expected if these in silico methods are employed in the elucidation of a common case of enzyme inhibition.Communicated by Ramaswamy H. Sarma.

Discovery of new potent lysine specific histone demythelase-1 inhibitors (LSD-1) using structure based and ligand based molecular modelling and machine learning

Lysine-specific histone demethylase 1 (LSD-1) is an epigenetic enzyme that oxidatively cleaves methyl groups from monomethyl and dimethyl Lys4 of histone H3 and is highly overexpressed in different types of cancer. Therefore, it has been widely recognized as a promising therapeutic target for cancer therapy. Towards this end, we employed various Computer Aided Drug Design (CADD) approaches including pharmacophore modelling and machine learning. Pharmacophores generated by structure-based (SB) (either crystallographic-based or docking-based) and ligand-based (LB) (either supervised or unsupervised) modelling methods were allowed to compete within the context of genetic algorithm/machine learning and were assessed by Shapley additive explanation values (SHAP) to end up with three successful pharmacophores that were used to screen the National Cancer Institute (NCI) database. Seventy-five NCI hits were tested for their LSD-1 inhibitory properties against neuroblastoma SH-SY5Y cells, pancreatic carcinoma Panc-1 cells, glioblastoma U-87 MG cells and in vitro enzymatic assay, culminating in 3 nanomolar LSD-1 inhibitors of novel chemotypes.

Describing inhibitor specificity for the amino acid transporter LAT1 from metainference simulations

The human L-type amino acid transporter 1 (LAT1; SLC7A5) is a membrane transporter of amino acids, thyroid hormones, and drugs such as the Parkinson's disease drug levodopa (L-Dopa). LAT1 is found in the blood-brain barrier, testis, bone marrow, and placenta, and its dysregulation has been associated with various neurological diseases, such as autism and epilepsy, as well as cancer. In this study, we combine metainference molecular dynamics simulations, molecular docking, and experimental testing, to characterize LAT1-inhibitor interactions. We first conducted a series of molecular docking experiments to identify the most relevant interactions between LAT1's substrate-binding site and ligands, including both inhibitors and substrates. We then performed metainference molecular dynamics simulations using cryoelectron microscopy structures in different conformations of LAT1 with the electron density map as a spatial restraint, to explore the inherent heterogeneity in the structures. We analyzed the LAT1 substrate-binding site to map important LAT1-ligand interactions as well as newly described druggable pockets. Finally, this analysis guided the discovery of previously unknown LAT1 ligands using virtual screening and cellular uptake experiments. Our results improve our understanding of LAT1-inhibitor recognition, providing a framework for rational design of future lead compounds targeting this key drug target.

Homology Modeling and Molecular Docking Approaches for the Proposal of Novel Insecticides against the African Malaria Mosquito (Anopheles gambiae)

Vector-borne infectious diseases are responsible for the deaths of over 700,000 people annually, than 400,000 of them resulting from malaria. The mosquito Anopheles gambiae is one of the dominant vector species of human malaria transmission. A significant issue of the conventional insecticides which target the arthropod borne infectious diseases is their induced resistance. To overcome this inconvenience, insecticides with new modes of action are required. One of the most promising targets for the development of new potential insecticides as evidenced by current studies is the D1-like dopamine receptor (DAR). To get a deeper understanding of the structural information of this receptor, the 3D homology model was built. The possible sites within the protein were identified and the most probable binding site was highlighted. The homology model along with a series of DAR antagonists with known activity against Anopheles gambiae larvae were used in docking experiments to gain insight into their intermolecular interactions. Furthermore, virtual screening of the natural compounds from the SPECS database led to the prediction of toxicity and environmental hazards for one potential new insecticide against the Anopheles gambiae mosquito.

Optimization of covalent docking for organophosphates interaction with Anopheles acetylcholinesterase

Organophosphates (OPs) used as potent insecticides for malaria vector control, covalently phosphorylate the catalytic serine residue of Anopheles gambiae AChE (AgAChE) in a reaction that liberates their leaving groups. In the recent 10-year insecticide use assessment, OPs were the most frequently used World Health Organization prequalified insecticides. Molecular modelling programs are best suited to display molecular interactions between ligands and the target proteins. The docking modes that generate ligand poses closer to the binding site show high accuracy in predicting the ligand binding mode. The implicit solvation approach such as molecular mechanics-generalized born surface area (MM-GBSA) is a more reliable method to predict ligand onformations and binding affinities. Apart from covalent docking studies being scarce, current molecular docking programs do not adequately possess the covalent docking reaction algorithm to display the molecular mechanism of OPs at the AgAChE catalytic site. This results into OP docking studies commonly being conducted through noncovalent pannels. The aim of this study was to establish the optimim covalent docking system for OPs through manual customization of Schrödinger's Glide covalent docking reaction algorithm. To achieve this, a newly customized covalent reaction algorithm was assessed on a set of ligands covering aromatic, non-aromatic and hydrophobic OPs and compared to the noncovalent docking results in terms of reliability based on the reported X-ray diffraction molecular interactions and crystal poses. The study established that by virtue of omitting the well-known OP hydrolysis, noncovalent mode suggested molecular interactions that were further from the catalytic triad and could not otherwise occur when the molecule is hydrolyzed as in the customized covalent docking mode. Moreover, the MM-GBSA concurred with the optimized covalent docking in eliminating such inaccurate molecular interactions. Additionally, the covalent docking mode confined the interactions and ligand poses to the catalytic site indicating relatively high accuracy and reliability. This study reports the optimized covalent docking panel that effectively confirmed the molecular mechanisms of OPs, as well as indentifying the corresponding amino acid residues required to stabilize the aromatic, non-aromatic and hydrophobic OPs at the AgAChE catalytic site in line with the reported X-ray diffraction studies. As such, the proposed manual customization of the Schrödinger's Glide covalent docking platform can be used to reliably predict molecular interactions between OPs and AgAChE target.

Highly Specific L-Type Amino Acid Transporter 1 Inhibition by JPH203 as a Potential Pan-Cancer Treatment

Accelerated cancer cell growth requires a massive intake of amino acids. Overexpression of L-type (large) amino acid transporter 1 (LAT1) on the cancer cell membrane facilitates such a demand, which is limited in normal organs. Therefore, LAT1 overexpression is ideal as a molecular cancer therapeutic target. JPH203, a LAT1-selective non-transportable blocker, had demonstrated LAT1 inhibition in <10 µM IC50 values and effectively suppressed cancer cell growth in studies involving several types of cancer cell lines and tumor xenograft models. A limited phase I clinical trial was performed on five different solid tumors and showed that JPH203 is well-tolerated and has a promising activity for the treatment of bile duct cancer. This review details the development and prospect of JPH203 as a LAT1-targeting cancer therapy.

Molecular mechanism of mitoquinol mesylate in mitigating the progression of hepatocellular carcinoma-in silico and in vivo studies

The safety and efficacy of mitoquinol mesylate (MitoQ) in attenuating the progression of hepatocellular carcinoma (HCC) in Wistar rats has been reported. However, the binding modes for MitoQ as well as its molecular mechanisms in cirrhosis and liver cancer have not been fully investigated. This study sought to understand the structural and molecular mechanisms of MitoQ in modulating the expression of nuclear factor erythroid 2-related factor 2 (Nrf2) and mitochondrial succinate dehydrogenase (SDH) in cirrhotic-HCC rats. The research indicates that the upregulated Nrf2 expression in cirrhotic-HCC rats was significantly (p < 0.05) reduced by MitoQ while the activity of SDH was significantly (p < 0.05) increased. Analysis of binding modes revealed MitoQ interacts with amino acid residues in the active pocket of tramtrack and bric-a-brac (BTB) and KELCH domains of KEAP1 with average binding affinities of -66.46 and -74.74 kcal/mol, respectively. Also, MitoQ interacted with the key amino acid residues at the active site of mitochondrial complex II with a higher average binding affinity of -75.76 kcal/mol compared to co-crystallized ligand of complex II (-62.31 kcal/mol). Molecular dynamics simulations data showed the binding of MitoQ to be stable with low eigenvalues while the quantum mechanics calculations suggest MitoQ to be very reactive with its mechanism of chemical reactivity to be via electrophilic reactions. Thus, MitoQ modulates expression of Nrf2 and enhances activity of mitochondrial SDH in cirrhotic-HCC rats via its interaction with key amino acid residues in the active pocket of BTB and KELCH domains of KEAP1 as well as amino residues at the active site of SDH. These findings are significant in demonstrating the potential of Nrf2 and SDH as possible biomarkers for the diagnosis and/or prognosis of hepatocellular carcinoma in patients. This study also supports repurposing of mitoQ for the treatment/management of liver cirrhosis and HCC.

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.

Mechanism of substrate transport and inhibition of the human LAT1-4F2hc amino acid transporter

LAT1 (SLC7A5) is one of the representative light chain proteins of heteromeric amino acid transporters, forming a heterodimer with its heavy chain partner 4F2hc (SLC3A2). LAT1 is overexpressed in many types of tumors and mediates the transfer of drugs and hormones across the blood-brain barrier. Thus, LAT1 is considered as a drug target for cancer treatment and may be exploited for drug delivery into the brain. Here, we synthesized three potent inhibitors of human LAT1, which inhibit transport of leucine with IC₅₀ values between 100 and 250 nM, and solved the cryo-EM structures of the corresponding LAT1-4F2hc complexes with these inhibitors bound at resolution of up to 2.7 or 2.8 Å. The protein assumes an outward-facing occluded conformation, with the inhibitors bound in the classical substrate binding pocket, but with their tails wedged between the substrate binding site and TM10 of LAT1. We also solved the complex structure of LAT1-4F2hc with 3,5-diiodo-l-tyrosine (Diiodo-Tyr) at 3.4 Å overall resolution, which revealed a different inhibition mechanism and might represent an intermediate conformation between the outward-facing occluded state mentioned above and the outward-open state. To our knowledge, this is the first time that the outward-facing conformation is revealed for the HAT family. Our results unveil more important insights into the working mechanisms of HATs and provide a structural basis for future drug design.

Novel, natural allosteric inhibitors and enhancers of Candida rugosa lipase activity

Candida rugosa lipase (CRL) is an enzyme commonly used in medicinal and biotechnological applications. Allosteric modulators of CRL could aid in modifying lipase-related diseases as well as improving biotechnological processes. Thus, a combinatorial approach of computational in-silico and high-throughput in-vitro screening was used to identify allosteric modulators of CRL. The screening of natural product libraries resulted in 132 compounds of which 53 were tested in-vitro. Subsequently, four inhibitors and three enhancers were identified of which rutin and cynaroside represented the strongest inhibitors of CRL activity (IC50: 227 ± 26 µM and 446 ± 15 µM, respectively) and NP-008496 the strongest enhancer (EC50: 425 ± 18 µM). All three compounds were predicted to bind the same allosteric site suggesting a common mechanism. Therefore, the present study demonstrated a reliable work-flow, identified an allosteric site of CRL and determined inhibitors and enhancers with numerous potential medical and biotechnological applications.

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.

L-Type amino acid transporter 1 as a target for drug delivery

Our growing understanding of membrane transporters and their substrate specificity has opened a new avenue in the field of targeted drug delivery. The L-type amino acid transporter 1 (LAT1) has been one of the most extensively investigated transporters for delivering drugs across biological barriers. The transporter is predominantly expressed in cerebral cortex, blood-brain barrier, blood-retina barrier, testis, placenta, bone marrow and several types of cancer. Its physiological function is to mediate Na⁺ and pH independent exchange of essential amino acids: leucine, phenylalanine, etc. Several drugs and prodrugs designed as LAT1 substrates have been developed to improve targeted delivery into the brain and cancer cells. Thus, the anti-parkinsonian drug, L-Dopa, the anti-cancer drug, melphalan and the anti-epileptic drug gabapentin, all used in clinical practice, utilize LAT1 to reach their target site. These examples provide supporting evidence for the utility of the LAT1-mediated targeted delivery of the (pro)drug. This review comprehensively summarizes recent advances in LAT1-mediated targeted drug delivery. In addition, the use of LAT1 is critically evaluated and limitations of the approach are discussed.

Glutamine transporters as pharmacological targets: From function to drug design

Among the different targets of administered drugs, there are membrane transporters that play also a role in drug delivery and disposition. Moreover, drug-transporter interactions are responsible for off-target effects of drugs underlying their toxicity. The improvement of the drug design process is subjected to the identification of those membrane transporters mostly relevant for drug absorption, delivery and side effect production. A peculiar group of proteins with great relevance to pharmacology is constituted by the membrane transporters responsible for managing glutamine traffic in different body districts. The interest around glutamine metabolism lies in its physio-pathological role; glutamine is considered a conditionally essential amino acid because highly proliferative cells have an increased request of glutamine that cannot be satisfied only by endogenous synthesis. Then, glutamine transporters provide cells with this special nutrient. Among the glutamine transporters, SLC1A5, SLC6A14, SLC6A19, SLC7A5, SLC7A8 and some members of SLC38 family are the best characterized, so far, in both physiological and pathological conditions. Few 3D structures have been solved by CryoEM; other structural data on these transporters have been obtained by computational analysis. Interactions with drugs have been described for several transporters of this group. For some of them, the studies are at an advanced stage, for others, the studies are still in nuce and novel biochemical findings open intriguing perspectives.