Proton Control of Transitions in an Amino Acid Transporter

Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives

Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.

Investigation of molecular details of a bacterial cationic amino acid transporter (GkApcT) during arginine transportation using molecular dynamics simulation and umbrella sampling techniques

CONTEXT

Cationic amino acid transporters (CATs) facilitate arginine transport across membranes and maintain its levels in various tissues and organs, but their overexpression has been associated with severe cancers. A recent study identified the alternating access mechanism and critical residues involved in arginine transportation in a cationic amino acid transporter from Geobacillus kaustophilus (GkApcT). Here, we used molecular dynamics (MD) simulation methods to investigate the transportation mechanism of arginine (Arg) through GkApcT. The results revealed that arginine strongly interacts with specific binding site residues (Thr43, Asp111, Glu115, Lys191, Phe231, Ile234, and Asp237). Based on the umbrella sampling, the main driving force for arginine transport is the polar interactions of the arginine with channel-lining residues. An in-depth description of the dissociation mechanism and binding energy analysis brings valuable insight into the interactions between arginine and transporter residues, facilitating the design of effective CAT inhibitors in cancer cells.

METHODS

The membrane-protein system was constructed by uploading the prokaryotic CAT (PDB ID: 6F34) to the CHARMM-GUI web server. Molecular dynamics simulations were done using the GROMACS package, version 5.1.4, with the CHARMM36 force field and TIP3P water model. The MM-PBSA approach was performed for determining the arginine binding free energy. Furthermore, the hotspot residues were identified through per-residue decomposition analysis. The characteristics of the channel such as bottleneck radius and channel length were analyzed using the CaverWeb 1.1 web server. The proton wire inside the transporter was investigated based on the classic Grotthuss mechanism. We also investigated the atomistic details of arginine transportation using the path-based free energy umbrella sampling technique (US).

Arginine transportation mechanism through cationic amino acid transporter 1: insights from molecular dynamics studies

Metabolic and signaling mechanisms in mammalian cells are facilitated by the transportation of L-arginine (Arg) across the plasma membrane through cationic amino acid transporter (CAT) proteins. Due to a lack of argininosuccinate synthase (ASS) activity in various tumor cells such as acute myeloid leukemia, acute lymphocytic leukemia, and chronic lymphocytic leukemia, these tumor entities are arginine-auxotrophic and therefore depend on the uptake of the amino acid arginine. Cationic amino acid transporter-1 (CAT-1) is the leading arginine importer expressed in the aforementioned tumor entities. Hence, in the present study, to investigate the transportation mechanism of arginine in CAT-1, we performed molecular dynamics (MD) simulation methods on the modeled human CAT-1. The MM-PBSA approach was conducted to determine the critical residues interacting with arginine within the corresponding binding site of CAT-1. In addition, we found out that the water molecules have the leading role in forming the transportation channel within CAT-1. The conductive structure of CAT-1 was formed only when the water molecules were continuously distributed across the channel. Steered molecular dynamics (SMD) simulation approach showed various energy barriers against arginine transportation through CAT-1, especially while crossing the bottlenecks of the related channel. These findings at the molecular level might shed light on identifying the crucial amino acids in the binding of arginine to eukaryotic CATs and also provide fundamental insights into the arginine transportation mechanisms through CAT-1. Understanding the transportation mechanism of arginine is essential to developing CAT-1 blockers, which can be potential medications for some types of cancers.Communicated by Ramaswamy H. Sarma.

Correction Schemes for Absolute Binding Free Energies Involving Lipid Bilayers

Absolute binding free-energy (ABFE) calculations are playing an increasing role in drug design, especially as they can be performed on a range of disparate compounds and direct comparisons between them can be made. It is, however, especially important to ensure that they are as accurate as possible, as unlike relative binding free-energy (RBFE) calculations, one does not benefit as much from a cancellation of errors during the calculations. In most modern implementations of ABFE calculations, a particle mesh Ewald scheme is typically used to treat the electrostatic contribution to the free energy. A central requirement of such schemes is that the box preserves neutrality throughout the calculation. There are many ways to deal with this problem that have been discussed over the years ranging from a neutralizing plasma with a post hoc correction term through to a simple co-alchemical ion within the same box. The post hoc correction approach is the most widespread. However, the vast majority of these studies have been applied to a soluble protein in a homogeneous solvent (water or salt solution). In this work, we explore which of the more common approaches would be the most suitable for a simulation box with a lipid bilayer within it. We further develop the idea of the so-called Rocklin correction for lipid-bilayer systems and show how such a correction could work. However, we also show that it will be difficult to make this generalizable in a practical way and thus we conclude that the use of a “co-alchemical ion” is the most useful approach for simulations involving lipid membrane systems.

Structural basis of antifolate recognition and transport by PCFT

Folates (also known as vitamin B9) have a critical role in cellular metabolism as the starting point in the synthesis of nucleic acids, amino acids and the universal methylating agent S-adenylsmethionine1,2. Folate deficiency is associated with a number of developmental, immune and neurological disorders3-5. Mammals cannot synthesize folates de novo; several systems have therefore evolved to take up folates from the diet and distribute them within the body3,6. The proton-coupled folate transporter (PCFT) (also known as SLC46A1) mediates folate uptake across the intestinal brush border membrane and the choroid plexus4,7, and is an important route for the delivery of antifolate drugs in cancer chemotherapy8-10. How PCFT recognizes folates or antifolate agents is currently unclear. Here we present cryo-electron microscopy structures of PCFT in a substrate-free state and in complex with a new-generation antifolate drug (pemetrexed). Our results provide a structural basis for understanding antifolate recognition and provide insights into the pH-regulated mechanism of folate transport mediated by PCFT.

A signal capture and proofreading mechanism for the KDEL-receptor explains selectivity and dynamic range in ER retrieval

ER proteins of widely differing abundance are retrieved from the Golgi by the KDEL-receptor. Abundant ER proteins tend to have KDEL rather than HDEL signals, whereas ADEL and DDEL are not used in most organisms. Here, we explore the mechanism of selective retrieval signal capture by the KDEL-receptor and how HDEL binds with 10-fold higher affinity than KDEL. Our results show the carboxyl-terminus of the retrieval signal moves along a ladder of arginine residues as it enters the binding pocket of the receptor. Gatekeeper residues D50 and E117 at the entrance of this pocket exclude ADEL and DDEL sequences. D50N/E117Q mutation of human KDEL-receptors changes the selectivity to ADEL and DDEL. However, further analysis of HDEL, KDEL, and RDEL-bound receptor structures shows that affinity differences are explained by interactions between the variable −4 H/K/R position of the signal and W120, rather than D50 or E117. Together, these findings explain KDEL-receptor selectivity, and how signal variants increase dynamic range to support efficient ER retrieval of low and high abundance proteins.

Sequence analysis and spatiotemporal developmental distribution of the Cat-1-type transporter slc7a1a in zebrafish (Danio rerio)

Cationic amino acid transporter 1 (Cat-1 alias Slc7a1) is a Na⁺-independent carrier system involved in transport and absorption of the cationic amino acids lysine, arginine, histidine, and ornithine and has also been shown to be indispensable in a large variety of biological processes. Starting from isolated full-length zebrafish (Danio rerio) cDNA for slc7a1a, we performed comparative and phylogenetic sequence analysis, investigated the conservation of the gene during vertebrate evolution, and defined tissue expression during zebrafish development. Whole mount in situ hybridization first detected slc7a1a transcripts in somites, eyes, and brain at 14 h post-fertilization (hpf) with additional expression in the distal nephron at 24 hpf and in branchial arches at 3 days post-fertilization (dpf), with significant increase by 5 dpf. Taken together, the expression analysis of the zebrafish Cat-1 system gene slc7a1a suggests a functional role(s) during the early development of the central nervous system, muscle, gills, and kidney. Graphical abstract.

The KDEL trafficking receptor exploits pH to tune the strength of an unusual short hydrogen bond

The endoplasmic reticulum (ER) is the main site of protein synthesis in eukaryotic cells and requires a high concentration of luminal chaperones to function. During protein synthesis, ER luminal chaperones are swept along the secretory pathway and must be retrieved to maintain cell viability. ER protein retrieval is achieved by the KDEL receptor, which recognises a C-terminal Lys-Asp-Glu-Leu (KDEL) sequence. Recognition of ER proteins by the KDEL receptor is pH dependent, with binding occurring under acidic conditions in the Golgi and release under conditions of higher pH in the ER. Recent crystal structures of the KDEL receptor in the apo and peptide bound state suggested that peptide binding drives the formation of a short-hydrogen bond that locks the KDEL sequence in the receptor and activates the receptor for COPI binding in the cytoplasm. Using quantum mechanical calculations we demonstrate that the strength of this short hydrogen bond is reinforced following protonation of a nearby histidine, providing a conceptual link between receptor protonation and KDEL peptide binding. Protonation also controls the water networks adjacent to the peptide binding site, leading to a conformational change that ultimately allows the receptor-complex to be recognized by the COPI system.