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

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.

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.

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

Asn25 Deamidation as an Allosteric Tool to Increase IFNß-1a Biological Activity.

Interferon beta (IFNß) is a well-known cytokine, belonging to the type I family, that exerts antiviral, immunomodulatory, and antiproliferative activity. It has been reported that the artificially deamidated form of recombinant IFNß-1a at Asn25 position shows an increased biological activity. As a deepening of the previous study, the molecular mechanism underlying this biological effect was investigated in this work by combining experimental and computational techniques. Specifically, the binding to IFNAR1 and IFNAR2 receptors and the canonical pathway of artificially deamidated IFNß-1a molecule were analyzed in comparison to the native form. As a result, a change in receptor affinity of deamidated IFNß-1a with respect to the native form was observed, and to better explore this molecular interaction, molecular dynamics simulations were carried out. Results confirmed, as previously hypothesized, that the N25D mutation can locally change the interaction network of the mutated residue but also that this effect can be propagated throughout the molecule. In fact, many residues not involved in the interaction with IFNAR1 in the native form participate to the recognition in the deamidated molecule, enhancing the binding to IFNAR1 receptor and consequently an increase of signaling cascade activation. In particular, a higher STAT1 phosphorylation and interferon-stimulated gene expression was observed under deamidated IFNß-1a cell treatment. In conclusion, this study increases the scientific knowledge of deamidated IFNß-1a, deciphering its molecular mechanism, and opens new perspectives to novel therapeutic strategies.

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.

IgG1 conformational behavior: elucidation of the N-glycosylation role via molecular dynamics.

Currently, monoclonal antibodies (mAbs) are the most used biopharmaceuticals for human therapy. One of the key aspects in their development is the control of effector functions mediated by the interaction between fragment crystallizable (Fc) and Fcγ receptors, which is a secondary mechanism of the action of biotherapeutics. N-glycosylation at the Fc portion can regulate these mechanisms, and much experimental evidence suggests that modifications of glycosidic chains can affect antibody binding to FcγRIIIa, consequently impacting the immune response. In this work, we try to elucidate via in silico procedures the structural role exhibited by glycans, particularly fucose, in mAb conformational freedom that can potentially affect the receptor recognition. By using adalimumab, a marketed IgG1, as a general template, after rebuilding its three-dimensional (3D) structure through homology modeling approaches, we carried out molecular dynamics simulations of three differently glycosylated species: aglycosylated, afucosylated, and fucosylated antibody. Trajectory analysis showed different dynamical behaviors and pointed out that sugars can influence the overall 3D structure of the antibody. As a result, we propose a putative structural mechanism by which the presence of fucose introduces conformational constraints in the whole antibody and not only in the Fc domain, preventing a conformation suitable for the interaction with the receptor. As secondary evidence, we observed a high flexibility of the antibodies that is translated into an asymmetric behavior of Fab portions shown by all the simulated biopolymers, making the dynamical asymmetry a new, to our knowledge, molecular aspect that may be further investigated. In conclusion, these findings can help understand the contribution of sugars on the structural architecture of mAbs, paving the way to novel strategies of pharmaceutical development.

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.

Surface Plasmon Resonance as a Tool for Ligand Binding Investigation of Engineered GPR17 Receptor, a G Protein Coupled Receptor Involved in Myelination.

The aim of this study was to investigate the potential of surface plasmon resonance (SPR) spectroscopy for the measurement of real-time ligand-binding affinities and kinetic parameters for GPR17, a G protein-coupled receptor (GPCR) of major interest in medicinal chemistry as potential target in demyelinating diseases. The receptor was directly captured, in a single-step, from solubilized membrane extracts on the sensor chip through a covalently bound anti-6x-His-antibody and retained its ligand binding activity for over 24 h. Furthermore, our experimental setup made possible, after a mild regeneration step, to remove the bound receptor without damaging the antibody, and thus to reuse many times the same chip. Two engineered variants of GPR17, designed for crystallographic studies, were expressed in insect cells, extracted from crude membranes and analyzed for their binding with two high affinity ligands: the antagonist Cangrelor and the agonist Asinex 1. The calculated kinetic parameters and binding constants of ligands were in good agreement with those reported from activity assays and highlighted a possible functional role of the N-terminal residues of the receptor in ligand recognition and binding. Validation of SPR results was obtained by docking and molecular dynamics of GPR17-ligands interactions and by functional in vitro studies. The latter allowed us to confirm that Asinex 1 behaves as GPR17 receptor agonist, inhibits forskolin-stimulated adenylyl cyclase pathway and promotes oligodendrocyte precursor cell maturation and myelinating ability.