In Silico Description of the Direct Inhibition Mechanism of Endothelial Lipase by ANGPTL3.

Angiopoietin-like protein 3 (ANGPTL3) is a plasmatic protein that plays a crucial role in lipoprotein metabolism by inhibiting the lipoprotein lipase (LPL) and the endothelial lipase (EL) responsible for the hydrolysis of phospholipids on high-density lipoprotein (HDL). Interest in developing new pharmacological therapies aimed at inhibiting ANGPTL3 has been growing due to the hypolipidemic and antiatherogenic profile observed in its absence. The goal of this study was the in silico characterization of the interaction between ANGPTL3 and EL. Because of the lack of any structural information on both the trimeric coiled-coil N-terminal domain of ANGPTL3 and the EL homodimer as well as data regarding their interactions, the first step was to obtain the three-dimensional model of these two proteins. The models were then refined via molecular dynamics (MD) simulations and used to investigate the interaction mechanism. The analysis of interactions in different docking poses and their refinement via MD allowed the identification of three specific glutamates of ANGPTL3 that recognize a positively charged patch on the surface of EL. These ANGPTL3 key residues, i.e., Glu154, Glu157, and Glu160, could form a putative molecular recognition site for EL. This study paves the way for future investigations aimed at confirming the recognition site and at designing novel inhibitors of ANGPTL3.

SEMA6A drives GnRH neuron-dependent puberty onset by tuning median eminence vascular permeability.

Innervation of the hypothalamic median eminence by Gonadotropin-Releasing Hormone (GnRH) neurons is vital to ensure puberty onset and successful reproduction. However, the molecular and cellular mechanisms underlying median eminence development and pubertal timing are incompletely understood. Here we show that Semaphorin-6A is strongly expressed by median eminence-resident oligodendrocytes positioned adjacent to GnRH neuron projections and fenestrated capillaries, and that Semaphorin-6A is required for GnRH neuron innervation and puberty onset. In vitro and in vivo experiments reveal an unexpected function for Semaphorin-6A, via its receptor Plexin-A2, in the control of median eminence vascular permeability to maintain neuroendocrine homeostasis. To support the significance of these findings in humans, we identify patients with delayed puberty carrying a novel pathogenic variant of SEMA6A. In all, our data reveal a role for Semaphorin-6A in regulating GnRH neuron patterning by tuning the median eminence vascular barrier and thereby controlling puberty onset.

GM1 structural requirements to mediate neuronal functions.

Since the 1980s, it has been known that the administration of ganglioside GM1 to cultured cells induced or enhanced neuronal differentiation. GM1 mechanism of action relies on its direct interaction and subsequent activation of the membrane tyrosine kinase receptor, TrkA, which naturally serves as NGF receptor. This process is mediated by the sole oligosaccharide portion of GM1, the pentasaccharide ß-Gal-(1-3)-ß-GalNAc-(1-4)-[α-Neu5Ac-(2-3)]-ß-Gal-(1-4)-ß-Glc. Here we detailed the minimum structural requirements of the oligosaccharide portion of GM1 for mediating the TrkA dependent neuritogenic processing. By in vitro and in silico biochemical approaches, we demonstrated that the minimal portion of GM1 required for the TrkA activation is the inner core of the ganglioside's oligosaccharide ß-Gal-(1-3)-ß-GalNAc-(1-4)-[α-Neu5Ac-(2-3)]-ß-Gal. The addition of a sialic acid residue at position 3 of the outer galactose of the GM1 oligosaccharide, which forms the oligosaccharide of GD1a, prevented the interaction with TrkA and the resulting neuritogenesis. On the contrary, the addition of a fucose residue at position 2 of the outer galactose, forming the Fucosyl-GM1 oligosaccharide, did not prevent the TrkA-mediated neuritogenesis.

Missense variants in RPH3A cause defects in excitatory synaptic function and are associated with a clinically variable neurodevelopmental disorder.

PURPOSE: RPH3A encodes a protein involved in the stabilization of GluN2A subunit of N-methyl-D-aspartate (NMDA)-type glutamate receptors at the cell surface, forming a complex essential for synaptic plasticity and cognition. We investigated the effect of variants in RPH3A in patients with neurodevelopmental disorders. METHODS: By using trio-based exome sequencing, GeneMatcher, and screening of 100,000 Genomes Project data, we identified 6 heterozygous variants in RPH3A. In silico and in vitro models, including rat hippocampal neuronal cultures, have been used to characterize the effect of the variants. RESULTS: Four cases had a neurodevelopmental disorder with untreatable epileptic seizures [p.(Gln73His)dn; p.(Arg209Lys); p.(Thr450Ser)dn; p.(Gln508His)], and 2 cases [p.(Arg235Ser); p.(Asn618Ser)dn] showed high-functioning autism spectrum disorder. Using neuronal cultures, we demonstrated that p.(Thr450Ser) and p.(Asn618Ser) reduce the synaptic localization of GluN2A; p.(Thr450Ser) also increased the surface levels of GluN2A. Electrophysiological recordings showed increased GluN2A-dependent NMDA ionotropic glutamate receptor currents for both variants and alteration of postsynaptic calcium levels. Finally, expression of the Rph3AThr450Ser variant in neurons affected dendritic spine morphology. CONCLUSION: Overall, we provide evidence that missense gain-of-function variants in RPH3A increase GluN2A-containing NMDA ionotropic glutamate receptors at extrasynaptic sites, altering synaptic function and leading to a clinically variable neurodevelopmental presentation ranging from untreatable epilepsy to autism spectrum disorder.

Amyposomes, a nanotechnological chaperone with anti-amyloidogenic activity.

AIM: The effect of liposomes bi-functionalized with phosphatidic acid and with a synthetic peptide derived from human apolipoprotein E has been evaluated on the aggregation features of different amyloidogenic proteins: human Amyloid ß1-40 (Aß1-40), transthyretin (TTR) variant S52P, human ß2microglobulin (ß2m) variants ΔN6 and D76N, Serum Amyloid A (SAA). METHODS: The formation of fibrillar aggregates of the proteins was investigated by ThioflavinT fluorescence assay and validated by Atomic Force Microscopy. RESULTS: The results show that liposomes are preventing the transition of non-aggregated forms to the fibrillar state, with stronger effects on Aß1-40, ß2m ΔN6 and SAA. Liposomes also induce disaggregation of the amyloid aggregates of all the proteins investigated, with stronger effects on Aß1-40, ß2 D76N and TTR.SPR assays show that liposomes bind Aß1-40 and SAA aggregates with high affinity (KD in the nanomolar range) whereas binding to TTR aggregates showed a lower affinity (KD in the micromolar range). Aggregates of ß2m variants showed both high and low affinity binding sites. Computed Structural analysis of protein fibrillar aggregates and considerations on the multidentate features of liposomes allow to speculate a common mechanism of action, based on binding the ß-stranded peptide regions responsible for the amyloid formation. CONCLUSION: Thus, multifunctional liposomes perform as pharmacological chaperones with anti-amyloidogenic activity, with a promising potential for the treatment of a number of protein-misfolding diseases.Key messageAmyloidosis is a group of diseases, each due to a specific protein misfolding.Anti-amyloidogenic nanoparticles have been gaining the utmost importance as a potential treatment for protein misfolding disorders.Liposomes bi-functionalized with phosphatidic acid and with a synthetic peptide derived from human apolipoprotein E showed anti-amyloidogenic activity.

Effect of Fc core fucosylation and light chain isotype on IgG1 flexibility.

N-glycosylation plays a key role in modulating the bioactivity of monoclonal antibodies (mAbs), as well as the light chain (LC) isotype can influence their physicochemical properties. However, investigating the impact of such features on mAbs conformational behavior is a big challenge, due to the very high flexibility of these biomolecules. In this work we investigate, by accelerated molecular dynamics (aMD), the conformational behavior of two commercial immunoglobulins G1 (IgG1), representative of κ and λ LCs antibodies, in both their fucosylated and afucosylated forms. Our results show, through the identification of a stable conformation, how the combination of fucosylation and LC isotype modulates the hinge behavior, the Fc conformation and the position of the glycan chains, all factors potentially affecting the binding to the FcγRs. This work also represents a technological enhancement in the conformational exploration of mAbs, making aMD a suitable approach to clarify experimental results.

Exploration of Novel Scaffolds Targeting Cytochrome b of Pyricularia oryzae.

The fulfilment of the European "Farm to Fork" strategy requires a drastic reduction in the use of "at risk" synthetic pesticides; this exposes vulnerable agricultural sectors-among which is the European risiculture-to the lack of efficient means for the management of devastating diseases, thus endangering food security. Therefore, novel scaffolds need to be identified for the synthesis of new and more environmentally friendly fungicides. In the present work, we employed our previously developed 3D model of P. oryzae cytochrome bc1 (cyt bc1) complex to perform a high-throughput virtual screening of two commercially available compound libraries. Three chemotypes were selected, from which a small collection of differently substituted analogues was designed and synthesized. The compounds were tested as inhibitors of the cyt bc1 enzyme function and the mycelium growth of both strobilurin-sensitive (WT) and -resistant (RES) P. oryzae strains. This pipeline has permitted the identification of thirteen compounds active against the RES cyt bc1 and five compounds that inhibited the WT cyt bc1 function while inhibiting the fungal mycelia only minimally. Serendipitously, among the studied compounds we identified a new chemotype that is able to efficiently inhibit the mycelium growth of WT and RES strains by ca. 60%, without inhibiting the cyt bc1 enzymatic function, suggesting a different mechanism of action.

The diachronic trend of female and male stature in Milan over 2000 years.

Stature is a biological trait directly determined by the interaction of genetic and environmental components. As such, it is often evaluated as an indicator for the reconstruction of skeletal biological profiles, past health, and social dynamics of human populations. Based on the analysis of 549 skeletons from the CAL (Collezione Antropologica LABANOF), a study of the diachronic trend of male and female adult stature in Milan (Italy) is being proposed here, covering a time span of about 2000 years, ranging from the Roman era to present-days. The skeletons, from necropolises dedicated to the less wealthy classes of Milanese society, were assigned to one of following five historical periods: Roman Era (first-fifth centuries AD), Early Middle Ages (sixth-tenth centuries AD), Late Middle Ages (eleventh-fifteenth centuries AD), Modern Era (sixteenth-eighteenth centuries AD) and Contemporary Era (nineteenth-twentieth centuries AD), and their stature was estimated according to the regression formulae of Trotter (1970). The collected data were then subjected to statistical analyses with ANOVA using R software. Although stature values showed an ample standard deviation in all periods, statistical analyses showed that stature did not significantly vary across historical periods in Milan for both sexes. This is one of the rare studies showing no diachronic changes in the trend of stature in Europe.

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.

Rabphilin-3A as a novel target to reverse α-synuclein-induced synaptic loss in Parkinson's disease.

Toxic aggregates of α-synuclein (αsyn) are considered key drivers of Parkinson's disease (PD) pathology. In early PD, αsyn induces synaptic dysfunction also modulating the glutamatergic neurotransmission. However, a more detailed understanding of the molecular mechanisms underlying αsyn-triggered synaptic failure is required to design novel therapeutic interventions. Here, we described the role of Rabphilin-3A (Rph3A) as novel target to counteract αsyn-induced synaptic loss in PD. Rph3A is a synaptic protein interacting with αsyn and involved in stabilizing dendritic spines and in promoting the synaptic retention of NMDA-type glutamate receptors. We found that in vivo intrastriatal injection of αsyn-preformed fibrils in mice induces the early loss of striatal synapses associated with decreased synaptic levels of Rph3A and impaired Rph3A/NMDA receptors interaction. Modulating Rph3A striatal expression or interfering with the Rph3A/αsyn complex with a small molecule prevented dendritic spine loss and rescued associated early motor defects in αsyn-injected mice. Notably, the same experimental approaches prevented αsyn-induced synaptic loss in vitro in primary hippocampal neurons. Overall, these findings indicate that approaches aimed at restoring Rph3A synaptic functions can slow down the early synaptic detrimental effects of αsyn aggregates in PD.

Molecular Modelling of NONO and SFPQ Dimerization Process and RNA Recognition Mechanism.

NONO and SFPQ are involved in multiple nuclear processes (e.g., pre-mRNA splicing, DNA repair, and transcriptional regulation). These proteins, along with NEAT1, enable paraspeckle formation, thus promoting multiple myeloma cell survival. In this paper, we investigate NONO and SFPQ dimer stability, highlighting the hetero- and homodimer structural differences, and model their interactions with RNA, simulating their binding to a polyG probe mimicking NEAT1guanine-rich regions. We demonstrated in silico that NONO::SFPQ heterodimerization is a more favorable process than homodimer formation. We also show that NONO and SFPQ RRM2 subunits are primarily required for protein-protein interactions with the other DBHS protomer. Simulation of RNA binding to NONO and SFPQ, beside validating RRM1 RNP signature importance, highlighted the role of ß2 and ß4 strand residues for RNA specific recognition. Moreover, we demonstrated the role of the NOPS region and other protomer's RRM2 ß2/ß3 loop in strengthening the interaction with RNA. Our results, having deepened RNA and DBHS dimer interactions, could contribute to the design of small molecules to modulate the activity of these proteins. RNA-mimetics, able to selectively bind to NONO and/or SFPQ RNA-recognition site, could impair paraspeckle formation, thus representing a first step towards the discovery of drugs for multiple myeloma treatment.

"From Protein Toxins to Applied Toxicological Testing" virtual workshop identifies the need for a bioinformatic framework to assess novel food protein safety.

On October 21-22, 2020 the HESI (Health and Environmental Sciences Institute) Protein Allergens, Toxins, and Bioinformatics Committee, and the Society of Toxicology Food Safety Specialty Section co-hosted a virtual workshop titled "From Protein Toxins to Applied Toxicological Testing". The workshop focused on the safety assessment of novel proteins contained in foods and feeds, was globally represented by over 200 stakeholder attendees, and featured contributions from experts in academia, government and non-government organizations, and agricultural biotechnology developers from the private sector. A range of topics relevant to novel protein safety were discussed, including: the state of protein toxin biology, modes and mechanisms of action, structures and activity, use of bioinformatic analyses to assess the safety of a protein, and ways to leverage computational biology with in silico approaches for protein toxin identification/characterization. Key outcomes of the workshop included the appreciation of the complexity of developing a definition for a protein toxin when viewed from the perspective of food and feed safety, confirming the need for a case-by-case hypothesis-driven interpretation of bioinformatic results that leverages additional metadata rather than an alignment threshold-driven interpretation, and agreement that a "toxin protein database" is not necessary, as the bioinformatic needs for toxin detection may be accomplished by existing databases such as Pfam and UniProtKB/Swiss-Prot. In this paper, a path forward is proposed.

Functional Heterodimerization between the G Protein-Coupled Receptor GPR17 and the Chemokine Receptors 2 and 4: New Evidence.

GPR17, a G protein-coupled receptor, is a pivotal regulator of myelination. Its endogenous ligands trigger receptor desensitization and downregulation allowing oligodendrocyte terminal maturation. In addition to its endogenous agonists, GPR17 could be promiscuously activated by pro-inflammatory oxysterols and chemokines released at demyelinating lesions. Herein, the chemokine receptors CXCR2 and CXCR4 were selected to perform both in silico modelling and in vitro experiments to establish their structural and functional interactions with GPR17. The relative propensity of GPR17 and CXCR2 or CXCR4 to form homo- and hetero-dimers was assessed by homology modelling and molecular dynamics (MD) simulations, and co-immunoprecipitation and immunoenzymatic assay. The interaction between chemokine receptors and GPR17 was investigated by determining receptor-mediated modulation of intracellular cyclic adenosine monophosphate (cAMP). Our data show the GPR17 association with CXCR2 or CXCR4 and the negative regulation of these interactions by CXCR agonists or antagonists. Moreover, GPR17 and CXCR2 heterodimers can functionally influence each other. In contrast, CXCR4 can influence GPR17 functionality, but not vice versa. According to MD simulations, all the dimers reached conformational stability and negative formation energy, confirming the experimental observations. The cross-talk between these receptors could play a role in the development of the neuroinflammatory milieu associated with demyelinating events.

Hemolymph proteins: An overview across marine arthropods and molluscs.

In this compilation we collect information about the main protein components in hemolymph and stress the continued interest in their study. The reasons for such an attention span several areas of biological, veterinarian and medical applications: from the notions for better dealing with the species - belonging to phylum Arthropoda, subphylum Crustacea, and to phylum Mollusca - of economic interest, to the development of 'marine drugs' from the peptides that, in invertebrates, act as antimicrobial, antifungal, antiprotozoal, and/or antiviral agents. Overall, the topic most often on focus is that of innate immunity operated by classes of pattern-recognition proteins. SIGNIFICANCE: The immune response in invertebrates relies on innate rather than on adaptive/acquired effectors. At a difference from the soluble and membrane-bound immunoglobulins and receptors in vertebrates, the antimicrobial, antifungal, antiprotozoal and/or antiviral agents in invertebrates interact with non-self material by targeting some common (rather than some highly specific) structural motifs. Developing this paradigm into (semi) synthetic pharmaceuticals, possibly optimized through the modeling opportunities offered by computational biochemistry, is one of the lessons today's science may learn from the study of marine invertebrates, and specifically of the proteins and peptides in their hemolymph.

Structural Investigation and Molecular Modeling Studies of Strobilurin-Based Fungicides Active against the Rice Blast Pathogen Pyricularia oryzae.

The increasing emergence of fungicide-resistant pathogens requires urgent solutions for crop disease management. Here, we describe a structural investigation of new fungicides obtained by combining strobilurin and succinate dehydrogenase inhibitor pharmacophores. We identified compounds endowed with very good activity against wild-type Pyricularia oryzae, combined in some cases with promising activity against strobilurin-resistant strains. The first three-dimensional model of P. oryzae cytochrome bc1 complex containing azoxystrobin as a ligand was developed. The model was validated with a set of commercially available strobilurins, and it well explains both the resistance mechanism to strobilurins mediated by the mutation G143A and the activity of metyltetraprole against strobilurin-resistant strains. The obtained results shed light on the key recognition determinants of strobilurin-like derivatives in the cytochrome bc1 active site and will guide the further rational design of new fungicides able to overcome resistance caused by G143A mutation in the rice blast pathogen.

rHDL modeling and the anchoring mechanism of LCAT activation.

Lecithin:cholesterol-acyl transferase (LCAT) plays a major role in cholesterol metabolism as it is the only extracellular enzyme able to esterify cholesterol. LCAT activity is required for lipoprotein remodeling and, most specifically, for the growth and maturation of HDLs. In fact, genetic alterations affecting LCAT functionality may cause a severe reduction in plasma levels of HDL-cholesterol with important clinical consequences. Although several hypotheses were formulated, the exact molecular recognition mechanism between LCAT and HDLs is still unknown. We employed a combination of structural bioinformatics procedures to deepen the insights into the HDL-LCAT interplay that promotes LCAT activation and cholesterol esterification. We have generated a data-driven model of reconstituted HDL (rHDL) and studied the dynamics of an assembled rHDL::LCAT supramolecular complex, pinpointing the conformational changes originating from the interaction between LCAT and apolipoprotein A-I (apoA-I) that are necessary for LCAT activation. Specifically, we propose a mechanism in which the anchoring of LCAT lid to apoA-I helices allows the formation of a hydrophobic hood that expands the LCAT active site and shields it from the solvent, allowing the enzyme to process large hydrophobic substrates.

Distinguishing mode of action of compounds inducing craniofacial malformations in zebrafish embryos to support dose-response modeling in combined exposures.

Knowledge on mode-of-action (MOA) is required to understand toxicological effects of compounds, notably in the context of risk assessment of mixtures. Such information is generally scarce, and often complicated by the existence of multiple MOAs per compound. Here, MOAs related to developmental craniofacial malformations were derived from literature, and assembled in a MOA network. A selection of gene expression markers was based on these MOAs. Next, these markers were verified by qPCR in zebrafish embryos, after exposure to reference compounds. These were: triazoles for inhibition of retinoic acid (RA) metabolism, AM580 and CD3254 for selective activation of respectively RA-receptor (RAR) and retinoid-X-receptor (RXR), dithiocarbamates for inhibition of lysyl oxidase, TCDD for activation of the aryl-hydrocarbon-receptor (AhR), VPA for inhibition of histone deacetylase (HDAC), and PFOS for activation of peroxisome proliferator-activated receptor-alpha (PPARα). Next, marker gene profiles for these reference compounds were used to map the profiles of test compounds to known MOAs. In this way, 2,4-dinitrophenol matched with the TCDD and RAR profiles, boric acid with RAR, endosulfan with PFOS, fenpropimorph with dithiocarbamates, PCB126 with AhR, and RA with triazoles and RAR profiles. Prochloraz showed no match. Activities of these compounds in ToxCast assays, and in silico analysis of binding affinity to the respective targets showed limited concordance with the marker gene expression profiles, but still confirmed the complex MOA profiles of reference and test compounds. Ultimately, this approach could be used to support modeling of mixture effects based on upfront knowledge of (dis)similarity of MOAs.

Development of an adverse outcome pathway for cranio-facial malformations: A contribution from in silico simulations and in vitro data.

Mixtures of substances sharing the same molecular initiating event (MIE) are supposed to induce additive effects. The proposed MIE for azole fungicides is CYP26 inhibition with retinoic acid (RA) local increase, triggering key events leading to craniofacial defects. Valproic acid (VPA) is supposed to imbalance RA-regulated gene expression trough histone deacetylases (HDACs) inhibition. The aim was to evaluate effects of molecules sharing the same MIE (azoles) and of such having (hypothetically) different MIEs but which are eventually involved in the same adverse outcome pathway (AOP). An in silico approach (molecular docking) investigated the suggested MIEs. Teratogenicity was evaluated in vitro (WEC). Abnormalities were modelled by PROAST software. The common target was the branchial apparatus. In silico results confirmed azole-related CYP26 inhibition and a weak general VPA inhibition on the tested HDACs. Unexpectedly, VPA showed also a weak, but not marginal, capability to enter the CYP 26A1 and CYP 26C1 catalytic sites, suggesting a possible role of VPA in decreasing RA catabolism, acting as an additional MIE. Our findings suggest a new more complex picture. Consequently two different AOPs, leading to the same AO, can be described. VPA MIEs (HDAC and CYP26 inhibition) impinge on the two converging AOPs.

Some more about dogs: Proteomics of neglected biological fluids.

We report in this manuscript what is known about the protein makeup of a selection of biological fluids in the domestic dog. The samples we review - amniotic and allantoic fluid, seminal fluid, saliva, bile, synovial fluid, tears - are still very poorly characterized in this species. For some of them we can present results from our own, mainly unpublished experiments. SIGNIFICANCE: The dog is one of the most widespread companion animals, and also of medical relevance as model species for some human diseases. Still, investigation of body fluids other than serum and urine is not so commonly undertaken, although - like in humans - also these sample types may have potential for diagnostic purposes. We compile published data about proteomes of fetal fluids, seminal plasma, saliva, bile, synovial fluid and tears, enriched by some yet unpublished data of our own (proteins of amniotic and allantoic fluid, tears). Closing gaps in our knowledge on dog proteins will further our understanding of (patho)physiological processes.

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.