Different Binding Modes of Small and Large Binders of GAT1

Well‐known inhibitors of the γ‐aminobutyric acid (GABA) transporter GAT1 share a common scaffold of a small cyclic amino acid linked by an alkyl chain to a moiety with two aromatic rings. Tiagabine, the only FDA‐approved GAT1 inhibitor, is a typical example. Some small amino acids such as (R)‐nipecotic acid are medium‐to‐strong binders of GAT1, but similar compounds, such as proline, are very weak binders. When substituted with 4,4‐diphenylbut‐3‐en‐1‐yl (DPB) or 4,4‐bis(3‐methylthiophen‐2‐yl)but‐3‐en‐1‐yl (BTB) groups, the resulting compounds have similar pK_i and pIC₅₀ values, even though the pure amino acids have very different values. To investigate if small amino acids and their substituted counterparts share a similar binding mode, we synthesized butyl‐, DPB‐, and BTB‐substituted derivatives of small amino acids. Supported by the results of docking studies, we propose different binding modes not only for unsubstituted und substituted, but also for strong‐ and weak‐binding amino acids. These data lead to the conclusion that following a fragment‐based approach, not pure but N‐butyl‐substituted amino acids should be used as starting points, giving a better estimate of the activity when a BTB or DPB substituent is added.     

Homology modeling of the serotonin transporter: insights into the primary escitalopram-binding site

The serotonin transporter (SERT) is one of the neurotransmitter transporters that plays a critical role in the regulation of endogenous amine concentrations and therefore is an important target for therapeutic agents affecting the central nervous system. The recently published, high resolution X-ray structure of the closely related amino acid transporter, Aquifex aeolicus leucine transporter (LeuT), provides an opportunity to develop a three-dimensional model of the structure of SERT. We present herein a homology model of SERT using LeuT as the template and containing escitalopram as a bound ligand. Our model explains selectivities known from mutational studies and varying ligand data, which are discussed and illustrated in the paper.     

A Toxoplasma gondii Oxopurine Transporter Binds Nucleobases and Nucleosides Using Different Binding Modes

Toxoplasma gondii is unable to synthesize purines de novo, instead salvages them from its environment, inside the host cell, for which they need high affinity carriers. Here, we report the expression of a T. gondii Equilibrative Nucleoside Transporter, Tg244440, in a Trypanosoma brucei strain from which nucleobase transporters have been deleted. Tg244440 transported hypoxanthine and guanine with similar affinity (K_m ~1 µM), while inosine and guanosine displayed K_i values of 4.05 and 3.30 µM, respectively. Low affinity was observed for adenosine, adenine, and pyrimidines, classifying Tg244440 as a high affinity oxopurine transporter. Purine analogues were used to probe the substrate-transporter binding interactions, culminating in quantitative models showing different binding modes for oxopurine bases, oxopurine nucleosides, and adenosine. Hypoxanthine and guanine interacted through protonated N1 and N9, and through unprotonated N3 and N7 of the purine ring, whereas inosine and guanosine mostly employed the ribose hydroxy groups for binding, in addition to N1H of the nucleobase. Conversely, the ribose moiety of adenosine barely made any contribution to binding. Tg244440 is the first gene identified to encode a high affinity oxopurine transporter in T. gondii and, to the best of our knowledge, the first purine transporter to employ different binding modes for nucleosides and nucleobases.     

Mutagenic Analysis of the Putative ABCC6 Substrate-Binding Cavity Using a New Homology Model

Inactivating mutations in ABCC6 underlie the rare hereditary mineralization disorder pseudoxanthoma elasticum. ABCC6 is an ATP-binding cassette (ABC) integral membrane protein that mediates the release of ATP from hepatocytes into the bloodstream. The released ATP is extracellularly converted into pyrophosphate, a key mineralization inhibitor. Although ABCC6 is firmly linked to cellular ATP release, the molecular details of ABCC6-mediated ATP release remain elusive. Most of the currently available data support the hypothesis that ABCC6 is an ATP-dependent ATP efflux pump, an un-precedented function for an ABC transporter. This hypothesis implies the presence of an ATP-binding site in the substrate-binding cavity of ABCC6. We performed an extensive mutagenesis study using a new homology model based on recently published structures of its close homolog, bovine Abcc1, to characterize the substrate-binding cavity of ABCC6. Leukotriene C4 (LTC₄), is a high-affinity substrate of ABCC1. We mutagenized fourteen amino acid residues in the rat ortholog of ABCC6, rAbcc6, that corresponded to the residues in ABCC1 found in the LTC₄ binding cavity. Our functional characterization revealed that most of the amino acids in rAbcc6 corresponding to those found in the LTC₄ binding pocket in bovine Abcc1 are not critical for ATP efflux. We conclude that the putative ATP binding site in the substrate-binding cavity of ABCC6/rAbcc6 is distinct from the bovine Abcc1 LTC₄-binding site.     

Recent Advances on P-Glycoprotein (ABCB1) Transporter Modelling with In Silico Methods

ABC transporters play a critical role in both drug bioavailability and toxicity, and with the discovery of the P-glycoprotein (P-gp), this became even more evident, as it plays an important role in preventing intracellular accumulation of toxic compounds. Over the past 30 years, intensive studies have been conducted to find new therapeutic molecules to reverse the phenomenon of multidrug resistance (MDR) ), that research has found is often associated with overexpression of P-gp, the most extensively studied drug efflux transporter; in MDR, therapeutic drugs are prevented from reaching their targets due to active efflux from the cell. The development of P-gp inhibitors is recognized as a good way to reverse this type of MDR, which has been the subject of extensive studies over the past few decades. Despite the progress made, no effective P-gp inhibitors to reverse multidrug resistance are yet on the market, mainly because of their toxic effects. Computational studies can accelerate this process, and in silico models such as QSAR models that predict the activity of compounds associated with P-gp (or analogous transporters) are of great value in the early stages of drug development, along with molecular modelling methods, which provide a way to explain how these molecules interact with the ABC transporter. This review highlights recent advances in computational P-gp research, spanning the last five years to 2022. Particular attention is given to the use of machine-learning approaches, drug–transporter interactions, and recent discoveries of potential P-gp inhibitors that could act as modulators of multidrug resistance.     

Identification of putative binding sites of P-glycoprotein based on its homology model

A homology model of P-glycoprotein based on the crystal structure of the multidrug transporter Sav1866 is developed, incorporated into a membrane environment, and optimized. The resulting model is analyzed in relation to the functional state and potential binding sites. The comparison of modeled distances to distances reported in experimental studies between particular residues suggests that the model corresponds most closely to the first ATP hydrolysis step of the protein transport cycle. Comparison to the protein 3D structure confirms this suggestion. Using SiteID and Site Finder programs three membrane related binding regions are identified: a region at the interface between the membrane and cytosol and two regions located in the transmembrane domains. The regions contain binding pockets of different size, orientation, and amino acids. A binding pocket located inside the membrane cavity is also identified. The pockets are analyzed in relation to amino acids shown experimentally to influence the protein function. The results suggest that the protein has multiple binding sites and may bind and/or release substrates in multiple pathways.     

Inhibitory Action of Antidepressants on Mouse Betaine/GABA Transporter (BGT1) Heterologously Expressed in Cell Cultures

Betaine/γ-aminobutyric acid (GABA) transporter (BGT1, SLC6A12) is a member of the Na(+)- and Cl(-)-dependent neurotransmitter transporter gene family with a homology to the GABA transporters (GATs), GAT1 (SLC6A1), GAT2 (SLC6A13) and GAT3 (SLC6A11) (HUGO nomenclature). Since antidepressants have been reported to inhibit GABA uptake, we examined those effects on mouse BGT1 (mBGT1) in comparison with other mouse GAT (mGAT) subtypes in the heterologously expressed cell cultures. All antidepressants tested here inhibited the [(3)H]GABA uptake through mBGT1 and mGATs in a rank order of potency with mBGT1 > mGAT1-3. Kinetic analyses for maprotilline, mianserine and trimipramine revealed that they inhibited mBGT1 and mGAT1 noncompetitively, except that mianserine competitively inhibited mBGT1. These results provided a clue to investigate the structure-function relationship of mBGT1 using antidepressants as a tool, leading to the identification of potential candidates for selective and specific inhibitors of mBGT1.     

Binding Mode Prediction of Evodiamine within Vanilloid Receptor TRPV1

Accurate assessment of the potential binding mode of drugs is crucial to computer-aided drug design paradigms. It has been reported that evodiamine acts as an agonist of the vanilloid receptor Transient receptor potential vanilloid-1 (TRPV1). However, the precise interaction between evodiamine and TRPV1 was still not fully understood. In this perspective, the homology models of TRPV1 were generated using the crystal structure of the voltage-dependent shaker family K(+) channel as a template. We then performed docking and molecular dynamics simulation to gain a better understanding of the probable binding modes of evodiamine within the TRPV1 binding pocket. There are no significant interspecies differences in evodiamine binding in rat, human and rabbit TRPV1 models. Pharmacophore modeling further provided confidence for the validity of the docking studies. This study is the first to shed light on the structural determinants required for the interaction between TRPV1 and evodiamine, and gives new suggestions for the rational design of novel TRPV1 ligands.     

Structural Model for the Binding Sites of Allosterically Potentiating Ligands on Nicotinic Acetylcholine Receptors

Current treatments of Alzheimer's disease include the allosteric potentiation of nicotinic acetylcholine receptor (nAChR) response. The location of the binding site for allosteric potentiating ligands (APLs) within the receptor is not yet fully understood. Based on homology models for the ligand binding domain of human α7, human α4β2, and chicken α7 receptors, as well as blind docking experiments with galanthamine, physostigmine, codeine, and 5HT, we identified T197 as an essential element of the APL binding site at the outer surface of the ligand binding domain (LBD) of nAChR. We also found the previously known galanthamine binding site in the region of K123 at the inside of the receptor funnel, which, however, was shown to not be part of the APL site. Our results are verified by site‐directed mutagenesis and electrophysiological experiments, and suggest that APL and ACh bind to different sites on nicotinic receptors and that allosteric potentiation may arise from a direct interplay between both these sites.     

Homology Model Adjustment and Ligand Screening with a Pseudoreceptor of the Human Histamine H4 Receptor

A new pseudoreceptor modeling method (PRPS) was applied to the refinement of a homology model of the human histamine H₄ receptor (H₄R), the prediction of a ligand binding site, and virtual screening. Retrieval of two new H₄R ligands demonstrates the biological relevance of the pseudoreceptor model and provides a means for finding new hits and leads in the early phases of drug discovery.

     

Homology modeling of NR2B modulatory domain of NMDA receptor and analysis of ifenprodil binding

NMDA receptors are glutamate-gated ion channels (iGluRs) that are involved in several important physiological functions such as neuronal development, synaptic plasticity, learning, and memory. Among iGluRs, NMDA receptors have been perhaps the most actively investigated for their role in chronic neurodegeneration such as Alzheimer's, Parkinson's, and Huntington's diseases. Recent studies have shown that the NTD of subunit NR2B modulates ion channel gating through the binding of allosteric modulators such as the prototypical compound ifenprodil. In the present paper, the construction of a three-dimensional model for the NR2B modulatory domain is described and docking calculations allow, for the first time, definition of the ifenprodil binding pose at an atomic level and fully explain all the available structure-activity relationships. Moreover, in an attempt to add further insight into the ifenprodil mechanism of action, as it is not completely clear if it binds and stabilizes an open or a closed conformation of the NR2B modulatory domain, a matter, which is fundamental for the rational design of NMDA antagonists, MD simulations followed by an MM-PBSA analysis were performed. These calculations reveal that the closed conformation of the R1-R2 domain, rather than the open, constitutes the high affinity binding site for ifenprodil and that a profound stabilization of the closed conformation upon ifenprodil binding occurs. Thus, for a rational design and/or for virtual screening experiments, the closed conformation of the R1-R2 domain should be taken into account and our 3D model can provide valuable hints for the design of NR2B-selective antagonists.     

Insight into the inhibition of human choline kinase: homology modeling and molecular dynamics simulations

A homology model of human choline kinase (CK-alpha) based on the X-ray crystallographic structure of C. elegans choline kinase (CKA-2) is presented. Molecular dynamics simulations performed on CK-alpha confirm the quality of the model, and also support the putative ATP and choline binding sites. A good correlation between the MD results and reported CKA-2 mutagenesis assays has been found for the main residues involved in catalytic activity. Preliminary docking studies performed on the CK-alpha model indicate that inhibitors can bind to the binding sites of both substrates (ATP and choline). A possible reason for inhibition of choline kinase by Ca(2+) ion is also proposed.     

Screening for New Inhibitors of Glycine Transporter 1 and 2 by Means of MS Binding Assays

A straightforward screening of a compound library comprising 2439 substances for the identification of new inhibitors for the neurotransmitter transporters GlyT1 and GlyT2 is described. Screening and full-scale competition experiments were performed using recently developed GlyT1 and GlyT2 MS Binding Assays. That way for both targets, GlyT1 and GlyT2, ligands were identified, which exhibited affinities (pK_i values) in the low micromolar to sub-micromolar range. The majority of these binders exhibit new chemical scaffolds in the class of GlyT1 and GlyT2 inhibitors, which could be of interest for the development of new ligands with improved affinities for the target proteins. Additionally, compounds with excellent fluorescent properties were found for GlyT2, which renders them promising compounds for future fluorescence-based techniques. All in all, this study demonstrates that MS Binding Assays represent a powerful technology platform also well suited for the screening of compound libraries in a highly reliable and effective manner.     

基于MD方法的HMX黏结剂优选研究

从优选高聚物黏结剂角度考虑,用分子动力学(MD)方法,对含少量液态黏结剂的奥克托今(HMX)炸药进行了MD模拟。选择端羟基聚丁二烯(HTPB)、聚叠氮缩水甘油醚(GAP)、硅橡胶和乙烯基酯(VE)4种黏结剂。通过计算HMX与黏结剂之间的结合能,选择结合能最大的一组的黏结剂作为最优黏结剂。结果表明,HTPB可作为最优黏结剂。     

Homology modeling and molecular dynamics simulations of the mu opioid receptor in a membrane-aqueous system

Three types of opioid receptors-mu, delta, and kappa-belong to the rhodopsin subfamily in the G protein-coupled receptor superfamily. With the recent characterization of the high-resolution X-ray crystal structure of bovine rhodopsin, considerable attention has been focused on molecular modeling of these transmembrane proteins. In this study, a homology model of the mu opioid receptor was constructed based on the X-ray crystal structure of bovine rhodopsin. A phospholipid bilayer was built around the receptor, and two water layers were placed on both surfaces of the lipid bilayer. Molecular-dynamics simulations were carried out by using CHARMM for the entire system, which consisted of 316 amino acid residues, 92 phospholipid molecules, 8327 water molecules, and 11 chloride counter ions-40 931 atoms altogether. The whole system was equilibrated for 250 ps followed by another 2 ns dynamic simulation. The opioid ligand naltrexone was docked into the optimized model, and the critical amino acid residues for binding were identified. The mu opioid receptor homology model optimized in a complete membrane-aqueous system should provide a good starting point for further characterization of the binding modes for opioid ligands. Furthermore, the method developed herein will be applicable to molecular model building to other opioid receptors as well as other GPCRs.     

Target Molecular Simulations of RecA Family Protein Filaments

Modeling of the RadA family mechanism is crucial to understanding the DNA SOS repair process. In a 2007 report, the archaeal RadA proteins function as rotary motors (linker region: I71-K88) such as shown in Figure 1. Molecular simulations approaches help to shed further light onto this phenomenon. We find 11 rotary residues (R72, T75-K81, M84, V86 and K87) and five zero rotary residues (I71, K74, E82, R83 and K88) in the simulations. Inclusion of our simulations may help to understand the RadA family mechanism.     

Exploring the Molecular Basis for Selective Binding of Homoserine Dehydrogenase from Mycobacterium leprae TN toward Inhibitors: A Virtual Screening Study

Homoserine dehydrogenase (HSD) from Mycobacterium leprae TN is an antifungal target for antifungal properties including efficacy against the human pathogen. The 3D structure of HSD has been firmly established by homology modeling methods. Using the template, homoserine dehydrogenase from Thiobacillus denitrificans (PDB Id 3MTJ), a sequence identity of 40% was found and molecular dynamics simulation was used to optimize a reliable structure. The substrate and co-factor-binding regions in HSD were identified. In order to determine the important residues of the substrate (l-aspartate semialdehyde (l-ASA)) binding, the ASA was docked to the protein; Thr163, Asp198, and Glu192 may be important because they form a hydrogen bond with HSD through AutoDock 4.2 software. After use of a virtual screening technique of HSD, the four top-scoring docking hits all seemed to cation–π ion pair with the key recognition residue Lys107, and Lys207. These ligands therefore seemed to be new chemotypes for HSD. Our results may be helpful for further experimental investigations.     

Analysis of Binding Determinants for Different Classes of Competitive and Noncompetitive Inhibitors of Glycine Transporters

Glycine transporters are interesting therapeutic targets as they play significant roles in glycinergic and glutamatergic systems. The search for new selective inhibitors of particular types of glycine transporters (GlyT-1 and GlyT-2) with beneficial kinetics is hampered by limited knowledge about the spatial structure of these proteins. In this study, a pool of homology models of GlyT-1 and GlyT-2 in different conformational states was constructed using the crystal structures of related transporters from the SLC6 family and the recently revealed structure of GlyT-1 in the inward-open state, in order to investigate their binding sites. The binding mode of the known GlyT-1 and GlyT-2 inhibitors was determined using molecular docking studies, molecular dynamics simulations, and MM-GBSA free energy calculations. The results of this study indicate that two amino acids, Gly373 and Leu476 in GlyT-1 and the corresponding Ser479 and Thr582 in GlyT-2, are mainly responsible for the selective binding of ligands within the S1 site. Apart from these, one pocket of the S2 site, which lies between TM3 and TM10, may also be important. Moreover, selective binding of noncompetitive GlyT-1 inhibitors in the intracellular release pathway is affected by hydrophobic interactions with Ile399, Met382, and Leu158. These results can be useful in the rational design of new glycine transporter inhibitors with desired selectivity and properties in the future.     

Comparative analysis of putative agonist-binding modes in the human A1 adenosine receptor

A recent study reported a model of the human A(1) adenosine receptor and its agonist binding site, proposing two putative binding modes in the same binding site for the natural agonist, adenosine. The present work investigates the flexibility of this binding site by exhaustive exploration with the natural agonist and with three other adenosine derivatives: N6-cyclopentyladenosine (CPA), 2-chloro-N6-cyclopentyladenosine (CCPA), and 5'-N-ethylcarboxamidoadenosine (NECA). Our aim was to find a common binding mode for agonists that would explain the role in the binding process of the different substitutions allowed at the 2, N6, and 5' positions of adenosine. This problem was addressed through docking simulations, molecular dynamics studies, and estimations of the ligand-binding free energy with both the AUTODOCK scoring function and the linear interaction energy (LIE) approach. The results point to a single receptor-binding position that explains the effects of the different chemical modifications on the adenosine derivatives considered here.