ProBiS-ligands: a web server for prediction of ligands by examination of protein binding sites

Development of a novel representation of drug 3D structures and enhancement of the TSR-based method for probing drug and target interactions

Understanding the mechanisms underlying interactions between drugs and target proteins is critical for drug discovery. In our earlier studies, we introduced the Triangular Spatial Relationship (TSR)-based algorithm, which enables the representation of a protein's 3D structure as a vector of integers (TSR keys). These TSR keys correspond to substructures of the 3D structure of a protein and are computed based on the triangles constructed by all possible triples of Cα atoms within the protein. In this study, we report on a new TSR-based algorithm for probing drug and target interactions. Specifically, we have extended the previous algorithm in three novel directions: TSR keys for representing the 3D structure of a drug or a ligand, cross TSR keys between drugs and their targets and intra-residual TSR keys for phosphorylated amino acids. The outcomes illustrate the key contributions as follows: (i) The TSR-based method, which uses the TSR keys as features, is unique in its capability to interpret hierarchical relationships of drugs as well as drug - target complexes using common and specific TSR keys. (ii) The method can distinguish not only the binding sites from the rest of the protein structures, but also the binding sites of primary targets from those of off-targets. (iii) The method has the potential to correlate the 3D structures of drugs with their functions. (iv) Representation of 3D structures by TSR keys has its unique advantage in terms of ease of making searching for similar substructures across structure datasets easier. In summary, this study presents a novel computational methodology, with significant advantages, for providing insights into the mechanism underlying drug and target interactions.

Structure Prediction and Genome Mining-Aided Discovery of the Bacterial C-Terminal Tryptophan Prenyltransferase PalQ

Post-translational prenylations, found in eukaryotic primary metabolites and bacterial secondary metabolites, play crucial roles in biomolecular interactions. Employing genome mining methods combined with AlphaFold2-based predictions of protein interactions, PalQ , a prenyltransferase responsible for the tryptophan prenylation of RiPPs produced by Paenibacillus alvei, is identified. PalQ differs from cyanobactin prenyltransferases because of its evolutionary relationship to isoprene synthases, which enables PalQ to transfer extended prenyl chains to the indole C3 position. This prenylation introduces structural diversity to the tryptophan side chain and also leads to conformational dynamics in the peptide backbone, attributed to the cis/trans isomerization that arises from the formation of a pyrrolidine ring. Additionally, PalQ exhibited pronounced positional selectivity for the C-terminal tryptophan. Such enzymatic characteristics offer a toolkit for peptide therapeutic lipidation.

PAS domain of flagellar histidine kinase FlrB has a unique architecture and binds heme as a sensory ligand in an unconventional fashion

Phosphorylation of the σ54-dependent transcription activator FlrC by the sensor histidine kinase FlrB is essential for flagellar synthesis of Vibrio cholerae. Despite that, the structure, sensory signal, and mechanistic basis of function of FlrB were elusive. Here, we report the crystal structure of the sensory PAS domain of FlrB in its functional dimeric state that exhibits a unique architecture. Series of biochemical/biophysical experiments on different constructs and mutants established that heme binds hydrophobically as sensory ligand in the shallow ligand-binding cleft of FlrB-PAS without axial coordination. Intriguingly, ATP binding to the C-terminal ATP-binding (CA) domain assists PAS domain to bind heme, vis-à-vis, heme binding to the PAS facilitates ATP binding to the CA domain. We hypothesize that synergistic binding of heme and ATP triggers conformational signaling in FlrB, leading to downstream flagellar gene transcription. Enhanced swimming motility of V. cholerae with increased heme uptake supports this proposition.

Identifying Metal Binding Sites in Proteins Using Homologous Structures, the MADE Approach

In order to identify the locations of metal ions in the binding sites of proteins, we have developed a method named the MADE (MAcromolecular DEnsity and Structure Analysis) approach. The MADE approach represents an evolution of our previous toolset, the ProBiS H2O (MD) methodology, for the identification of conserved water molecules. Our method uses experimental structures of proteins homologous to a query, which are subsequently superimposed upon it. Areas with a particular species present in a similar location among many homologous protein structures are identified using a clustering algorithm. Dense clusters likely represent positions containing species important to the query protein structure or function. We analyze well-characterized apo protein structures and show that the MADE approach can identify clusters corresponding to the expected positions of metal ions in their binding sites. The greatest advantage of our method lies in its generality. It can in principle be applied to any species found in protein records; it is not only limited to metal ions. We additionally demonstrate that the MADE approach can be successfully applied to predict the location of cofactors in computer-modeled structures, e.g., via AlphaFold. We also conduct a careful protein superposition method comparison and find our methodology robust and the results largely independent of the selected protein superposition algorithm. We postulate that with increasing structural data availability, additional applications of the MADE approach will be possible such as non-protein systems, water network identification, protein binding site elaboration, and analysis of binding events, all in a dynamic manner. We have implemented the MADE approach as a plugin for the PyMOL molecular visualization tool. The MADE plugin is available free of charge at https://gitlab.com/Jukic/made_software.

Identified Isosteric Replacements of Ligands' Glycosyl Domain by Data Mining

Biologically equivalent replacements of key moieties in molecules rationalize scaffold hopping, patent busting, or R-group enumeration. Yet, this information may depend upon the expert-defined space, and might be subjective and biased toward the chemistries they get used to. Most importantly, these practices are often informatively incomplete since they are often compromised by a try-and-error cycle, and although they depict what kind of substructures are suitable for the replacement occurrence, they fail to explain the driving forces to support such interchanges. The protein data bank (PDB) encodes a receptor-ligand interaction pattern and could be an optional source to mine structural surrogates. However, manual decoding of PDB has become almost impossible and redundant to excavate the bioisosteric know-how. Therefore, a text parsing workflow has been developed to automatically extract the local structural replacement of a specific structure from PDB by finding spatial and steric interaction overlaps between the fragments in endogenous ligands and particular ligand fragments. Taking the glycosyl domain for instance, a total of 49 520 replacements that overlap on nucleotide ribose were identified and categorized based on their SMILE codes. A predominately ring system, such as aliphatic and aromatic rings, was observed; yet, amide and sulfonamide replacements also occur. We believe these findings may enlighten medicinal chemists on the structure design and optimization of ligands using the bioisosteric replacement strategy.

Structure-based discovery and in vitro validation of inhibitors of chloride intracellular channel 4 protein

The use of computer-aided methods have continued to propel accelerated drug discovery across various disease models, interestingly allowing the specific inhibition of pathogenic targets. Chloride Intracellular Channel Protein 4 (CLIC4) is a novel class of intracellular ion channel highly implicated in tumor and vascular biology. It regulates cell proliferation, apoptosis and angiogenesis; and is involved in multiple pathologic signaling pathways. Absence of specific inhibitors however impedes its advancement to translational research. Here, we integrate structural bioinformatics and experimental research approaches for the discovery and validation of small-molecule inhibitors of CLIC4. High-affinity allosteric binders were identified from a library of 1615 Food and Drug Administration (FDA)-approved drugs via a high-performance computing-powered blind-docking approach, resulting in the selection of amphotericin B and rapamycin. NMR assays confirmed the binding and conformational disruptive effects of both drugs while they also reversed stress-induced membrane translocation of CLIC4 and inhibited endothelial cell migration. Structural and dynamics simulation studies further revealed that the inhibitory mechanisms of these compounds were hinged on the allosteric modulation of the catalytic glutathione (GSH)-like site loop and the extended catalytic β loop which may elicit interference with the catalytic activities of CLIC4. Structure-based insights from this study provide the basis for the selective targeting of CLIC4 to treat the associated pathologies.

Protein binding sites for drug design

Drug development is a lengthy and challenging process that can be accelerated at early stages by new mathematical approaches and modern computers. To address this important issue, we are developing new mathematical solutions for the detection and characterization of protein binding sites that are important for new drug development. In this review, we present algorithms based on graph theory combined with molecular dynamics simulations that we have developed for studying biological target proteins to provide important data for optimizing the early stages of new drug development. A particular focus is the development of new protein binding site prediction algorithms (ProBiS) and new web tools for modeling pharmaceutically interesting molecules-ProBiS Tools (algorithm, database, web server), which have evolved into a full-fledged graphical tool for studying proteins in the proteome. ProBiS differs from other structural algorithms in that it can align proteins with different folds without prior knowledge of the binding sites. It allows detection of similar binding sites and can predict molecular ligands of various types of pharmaceutical interest that could be advanced to drugs to treat a disease, based on the entire Protein Data Bank (PDB) and AlphaFold database, including proteins not yet in the PDB. All ProBiS Tools are freely available to the academic community at http://insilab.org and https://probis.nih.gov.

ProBiS-Fold Approach for Annotation of Human Structures from the AlphaFold Database with No Corresponding Structure in the PDB to Discover New Druggable Binding Sites

ProBiS (Protein Binding Sites), a local structure-based comparison algorithm, is used in the new ProBiS-Fold web server to annotate human structures from the AlphaFold database without a corresponding structure in the Protein Data Bank (PDB) to discover new druggable binding sites. The ProBiS algorithm is used to compare each query protein structure predicted by the AlphaFold approach with the protein structures from the PDB to identify similarities between known binding sites found in the PDB and yet unknown binding sites in the AlphaFold database. Ligands bound in these identified similar PDB sites are then transferred to each query protein from the AlphaFold database, and binding sites are identified as ligand clusters on an AlphaFold protein. Small molecule binding sites are assigned druggability scores based on the similarity of their ligands to known drugs, allowing them to be ranked according to their perceived and actual potential for drug development. ProBiS-Fold provides interactive and downloadable binding sites for the entire human structural proteome, including more than 3000 new druggable binding sites that have no corresponding structure in the PDB, taking into account AlphaFold's model quality, to enable protein structure-function relationship studies and pharmaceutical drug discovery research. The web server is freely accessible to academic users at http://probis-fold.insilab.org.

GraphSite: Ligand Binding Site Classification with Deep Graph Learning

The binding of small organic molecules to protein targets is fundamental to a wide array of cellular functions. It is also routinely exploited to develop new therapeutic strategies against a variety of diseases. On that account, the ability to effectively detect and classify ligand binding sites in proteins is of paramount importance to modern structure-based drug discovery. These complex and non-trivial tasks require sophisticated algorithms from the field of artificial intelligence to achieve a high prediction accuracy. In this communication, we describe GraphSite, a deep learning-based method utilizing a graph representation of local protein structures and a state-of-the-art graph neural network to classify ligand binding sites. Using neural weighted message passing layers to effectively capture the structural, physicochemical, and evolutionary characteristics of binding pockets mitigates model overfitting and improves the classification accuracy. Indeed, comprehensive cross-validation benchmarks against a large dataset of binding pockets belonging to 14 diverse functional classes demonstrate that GraphSite yields the class-weighted F1-score of 81.7%, outperforming other approaches such as molecular docking and binding site matching. Further, it also generalizes well to unseen data with the F1-score of 70.7%, which is the expected performance in real-world applications. We also discuss new directions to improve and extend GraphSite in the future.

Virtual screening and docking analysis of novel ligands for selective enhancement of tea (Camellia sinensis) flavonoids

•

Tea-specific flavonoid biosynthetic pathway (FBP) was retrieved from KEGG.

•

Putative ligands were predicted to enhance enzymes-substrate binding affinity.

•

FBP genes showed moderately higher expression & relatively strong codon adaptation.

•

Most of the genes were AT-rich and biased to A/U-ending synonymous codons.

•

Mutational selection was determining the selective constraints on codon bias.

Flavour of tea is mainly contributed by a group of polyphenols – flavonoids. However, the content of flavonoid fluctuates seasonally and is found to be higher in the first flush of tea, when compared to the second flush. This disparity in the flavonoid content, and hence taste, incurs heavy economic losses to the tea plantation industry each harvest season. For our present study, four key product-specific enzymes (PAL, FNS, FLS and ANS) of the tea-specific flavonoid pathway were selected to perform molecular docking studies with specific virtually screened allosteric modulators. Results of docking analyses showed Naringenin, 2-Morpholin-4-ium-4-ylethanesulfonate, 6-C-Glucosylquercetin, 2-Oxoglutaric acid, 3,5,7,3′,4′-pentahydroxyflavone to be capable of improving the spontaneity of the enzyme-substrate reactions in terms of docking score, RMSD values, and non-covalent interactions (H-bond,hydrophobic interaction, Π-stacking, salt bridge, etc.). Further, the evolutionary relationship of tea flavonoid pathway enzymes was constructed and compared with related taxa. The codon usage-based of tea flavonoid biosynthetic genes indicated the non-biasness of their nucleotide composition. Overall this study will provide a direction towards putative ligand-dependent enhancement of flavonoid content, irrespective of seasonal variation.

ProBiS-Dock: A Hybrid Multitemplate Homology Flexible Docking Algorithm Enabled by Protein Binding Site Comparison

The protein data bank (PDB) is a rich source of protein ligand structures, but ligands are not explicitly used in current docking algorithms. We have developed ProBiS-Dock, a docking algorithm complementary to the ProBiS-Dock Database (J. Chem. Inf. Model. 2021, 61, 4097-4107) that treats small molecules and proteins as fully flexible entities and allows conformational changes in both after ligand binding. A new scoring function is described that consists of a binding site-specific scoring function (ProBiS-Score) and a general statistical scoring function. ProBiS-Dock enables rapid docking of small molecules to proteins and has been successfully validated in silico against standard benchmarks. It enables rapid search for new active ligands by leveraging existing knowledge in the PDB. The potential of the software for drug development has been confirmed in vitro by the discovery of new inhibitors of human indoleamine 2,3-dioxygenase 1, an enzyme that is an attractive target for cancer therapy and catalyzes the first rate-determining step of l-tryptophan metabolism via the kynurenine pathway. The software is freely available to academic users at http://insilab.org/probisdock.

Drug repurposing for cancer therapy in the era of precision medicine

Drug repurposing refers to the identification of clinically approved drugs, with the known safety profiles and defined pharmacokinetic properties, to new indications. Despite the advances in oncology research, cancers are still associated with the most unmet medical needs. Drug repurposing has emerged as a useful approach for the search for effective and durable cancer treatment. It may also represent a promising strategy to facilitate precision cancer treatment and to overcome drug resistance. The repurposing of non-cancer drugs for precision oncology effectively extends the inventory of actionable molecular targets and thus increases the number of patients who may benefit from precision cancer treatment. In cancer types where genetic heterogeneity is so high that it is not feasible to identify strong repurposed drug candidates for standard treatment, the precision oncology approach offers individual patients access to novel treatment options. For repurposed candidates with low potency, a combination of multiple repurposed drugs may produce a synergistic therapeutic effect. Precautions should be taken when combining repurposed drugs with anticancer agents to avoid detrimental drug-drug interactions and unwanted side effects. New multifactorial data analysis and artificial intelligence methods are needed to untangle the complex association of molecular signatures influencing specific cancer subtypes to facilitate drug repurposing in precision oncology.

CIFDock: A novel CHARMM-based flexible receptor-flexible ligand docking protocol

Docking studies play a critical role in the current workflow of drug discovery. However, limitations may often arise through factors including inadequate ligand sampling, a lack of protein flexibility, scoring function inadequacies (e.g., due to metals, co-factors, etc.), and difficulty in retaining explicit water molecules. Herein, we present a novel CHARMM-based induced fit docking (CIFDock) workflow that can circumvent these limitations by employing all-atom force fields coupled to enhanced sampling molecular dynamics procedures. Self-guided Langevin dynamics simulations are used to effectively sample relevant ligand conformations, side chain orientations, crystal water positions, and active site residue motion. Protein flexibility is further enhanced by dynamic sampling of side chain orientations using an expandable rotamer library. Steps in the procedure consisting of fixing individual components (e.g., the ligand) while sampling the other components (e.g., the residues in the active site of the protein) allow for the complex to adapt to conformational changes. Ultimately, all components of the complex-the protein, ligand, and waters-are sampled simultaneously and unrestrained with SGLD to capture any induced fit effects. This modular flexible docking procedure is automated using CHARMM scripting, interfaced with SLURM array processing, and parallelized to use the desired number of processors. We validated the CIFDock procedure by performing cross-docking studies using a data set comprised of 21 pharmaceutically relevant proteins. Five variants of the CHARMM-based SWISSDOCK scoring functions were created to quantify the results of the final generated poses. Results obtained were comparable to, or in some cases improved upon, commercial docking program data.

MolADI: A Web Server for Automatic Analysis of Protein-Small Molecule Dynamic Interactions

Protein-ligand interaction analysis is important for drug discovery and rational protein design. The existing online tools adopt only a single conformation of the complex structure for calculating and displaying the interactions, whereas both protein residues and ligand molecules are flexible to some extent. The interactions evolved with time in the trajectories are of greater interest. MolADI is a user-friendly online tool which analyzes the protein-ligand interactions in detail for either a single structure or a trajectory. Interactions can be viewed easily with both 2D graphs and 3D representations. MolADI is available as a web application.

ProBiS-Dock Database: A Web Server and Interactive Web Repository of Small Ligand-Protein Binding Sites for Drug Design

We have developed a new system, ProBiS-Dock, which can be used to determine the different types of protein binding sites for small ligands. The binding sites identified this way are then used to construct a new binding site database, the ProBiS-Dock Database, that allows for the ranking of binding sites according to their utility for drug development. The newly constructed database currently has more than 1.4 million binding sites and offers the possibility to investigate potential drug targets originating from different biological species. The interactive ProBiS-Dock Database, a web server and repository that consists of all small-molecule ligand binding sites in all of the protein structures in the Protein Data Bank, is freely available at http://probis-dock-database.insilab.org. The ProBiS-Dock Database will be regularly updated to keep pace with the growth of the Protein Data Bank, and our anticipation is that it will be useful in drug discovery.

Ligand Binding Site Comparison - LiBiSCo - a web-based tool for analyzing interactions between proteins and ligands to explore amino acid specificity within active sites

Interaction between protein and ligands are ubiquitous in a biological cell, and understanding these interactions at the atom level in protein-ligand complexes is crucial for structural bioinformatics and drug discovery. Here, we present a web-based protein-ligand interaction application named Ligand Binding Site Comparison (LiBiSCo) for comparing the amino acid residues interacting with atoms of a ligand molecule between different protein-ligand complexes available in the Protein Data Bank (PDB) database. The comparison is performed at the ligand atom level irrespectively of having binding site similarity or not between the protein structures of interest. The input used in LiBiSCo is one or several PDB IDs of protein-ligand complex(es) and the tool returns a list of identified interactions at ligand atom level including both bonded and non-bonded interactions. A sequence profile for the interaction for each ligand atoms is provided as a WebLogo. The LiBiSco is useful in understanding ligand binding specificity and structural promiscuity among families that are structurally unrelated. The LiBiSCo tool can be accessed through https://albiorix.bioenv.gu.se/LiBiSCo/HomePage.py.

NOD: a web server to predict New use of Old Drugs to facilitate drug repurposing

Computational methods accelerate the drug repurposing pipelines that are a quicker and cost-effective alternative to discovering new molecules. However, there is a paucity of web servers to conduct fast, focussed, and customized investigations for identifying new uses of old drugs. We present the NOD web server, which has the mentioned characteristics. NOD uses a sensitive sequence-guided approach to identify close and distant homologs of a protein of interest. NOD then exploits this evolutionary information to suggest potential compounds from the DrugBank database that can be repurposed against the input protein. NOD also allows expansion of the chemical space of the potential candidates through similarity searches. We have validated the performance of NOD against available experimental and/or clinical reports. In 65.6% of the investigated cases in a control study, NOD is able to identify drugs more effectively than the searches made in DrugBank. NOD is freely-available at http://pauling.mbu.iisc.ac.in/NOD/NOD/ .

Exploring the effectiveness of the TSR-based protein 3-D structural comparison method for protein clustering, and structural motif identification and discovery of protein kinases, hydrolases, and SARS-CoV-2's protein via the application of amino acid grouping

Development of protein 3-D structural comparison methods is essential for understanding protein functions. Some amino acids share structural similarities while others vary considerably. These structures determine the chemical and physical properties of amino acids. Grouping amino acids with similar structures potentially improves the ability to identify structurally conserved regions and increases the global structural similarity between proteins. We systematically studied the effects of amino acid grouping on the numbers of Specific/specific, Common/common, and statistically different keys to achieve a better understanding of protein structure relations. Common keys represent substructures found in all types of proteins and Specific keys represent substructures exclusively belonging to a certain type of proteins in a data set. Our results show that applying amino acid grouping to the Triangular Spatial Relationship (TSR)-based method, while computing structural similarity among proteins, improves the accuracy of protein clustering in certain cases. In addition, applying amino acid grouping facilitates the process of identification or discovery of conserved structural motifs. The results from the principal component analysis (PCA) demonstrate that applying amino acid grouping captures slightly more structural variation than when amino acid grouping is not used, indicating that amino acid grouping reduces structure diversity as predicted. The TSR-based method uniquely identifies and discovers binding sites for drugs or interacting proteins. The binding sites of nsp16 of SARS-CoV-2, SARS-CoV and MERS-CoV that we have defined will aid future antiviral drug design for improving therapeutic outcome. This approach for incorporating the amino acid grouping feature into our structural comparison method is promising and provides a deeper insight into understanding of structural relations of proteins.

Mechanistic Insights into Side Effects of Troglitazone and Rosiglitazone Using a Novel Inverse Molecular Docking Protocol

Thiazolidinediones form drugs that treat insulin resistance in type 2 diabetes mellitus. Troglitazone represents the first drug from this family, which was removed from use by the FDA due to its hepatotoxicity. As an alternative, rosiglitazone was developed, but it was under the careful watch of FDA for a long time due to suspicion, that it causes cardiovascular diseases, such as heart failure and stroke. We applied a novel inverse molecular docking protocol to discern the potential protein targets of both drugs. Troglitazone and rosiglitazone were docked into predicted binding sites of >67,000 protein structures from the Protein Data Bank and examined. Several new potential protein targets with successfully docked troglitazone and rosiglitazone were identified. The focus was devoted to human proteins so that existing or new potential side effects could be explained or proposed. Certain targets of troglitazone such as 3-oxo-5-beta-steroid 4-dehydrogenase, neutrophil collagenase, stromelysin-1, and VLCAD were pinpointed, which could explain its hepatoxicity, with additional ones indicating that its application could lead to the treatment/development of cancer. Results for rosiglitazone discerned its interaction with members of the matrix metalloproteinase family, which could lead to cancer and neurodegenerative disorders. The concerning cardiovascular side effects of rosiglitazone could also be explained. We firmly believe that our results deepen the mechanistic understanding of the side effects of both drugs, and potentially with further development and research maybe even help to minimize them. On the other hand, the novel inverse molecular docking protocol on the other hand carries the potential to develop into a standard tool to predict possible cross-interactions of drug candidates potentially leading to adverse side effects.

Proteins and Their Interacting Partners: An Introduction to Protein-Ligand Binding Site Prediction Methods with a Focus on FunFOLD3

Proteins are essential molecules with a diverse range of functions; elucidating their biological and biochemical characteristics can be difficult and time consuming using in vitro and/or in vivo methods. Additionally, in vivo protein-ligand binding site elucidation is unable to keep place with current growth in sequencing, leaving the majority of new protein sequences without known functions. Therefore, the development of new methods, which aim to predict the protein-ligand interactions and ligand-binding site residues directly from amino acid sequences, is becoming increasingly important. In silico prediction can utilise either sequence information, structural information or a combination of both. In this chapter, we will discuss the broad range of methods for ligand-binding site prediction from protein structure and we will describe our method, FunFOLD3, for the prediction of protein-ligand interactions and ligand-binding sites based on template-based modelling. Additionally, we will describe the step-by-step instructions using the FunFOLD3 downloadable application along with examples from the Critical Assessment of Techniques for Protein Structure Prediction (CASP) where FunFOLD3 has been used to aid ligand and ligand-binding site prediction. Finally, we will introduce our newer method, FunFOLD3-D, a version of FunFOLD3 which aims to improve template-based protein-ligand binding site prediction through the integration of docking, using AutoDock Vina.

The latest automated docking technologies for novel drug discovery

INTRODUCTION

Molecular docking has been consolidated as one of the most important methods in the molecular modeling field. It has been recognized as a prominent tool in the study of protein-ligand complexes, to describe intermolecular interactions, to accurately predict poses of multiple ligands, to discover novel promising bioactive compounds. Molecular docking methods have evolved in terms of their accuracy and reliability; but there are pending issues to solve for improving the connection between the docking results and the experimental evidence.

AREAS COVERED

In this article, the author reviews very recent innovative molecular docking applications with special emphasis on reverse docking, treatment of protein flexibility, the use of experimental data to guide the selection of docking poses, the application of Quantum mechanics(QM) in docking, and covalent docking.

EXPERT OPINION

There are several issues being worked on in recent years that will lead to important breakthroughs in molecular docking methods in the near future These developments are related to more efficient exploration of large datasets and receptor conformations, advances in electronic description, and the use of structural information for guiding the selection of results.

High resolution structure of human apolipoprotein (a) kringle IV type 2: beyond the lysine binding site

Lipoprotein (a) [Lp(a)] is characterized by an LDL-like composition in terms of lipids and apoB100, and by one copy of a unique glycoprotein, apo(a). The apo(a) structure is mainly based on the repetition of tandem kringle domains with high homology to plasminogen kringles 4 and 5. Among them, kringle IV type 2 (KIV-2) is present in a highly variable number of genetically encoded repeats, whose length is inversely related to Lp(a) plasma concentration and cardiovascular risk. Despite it being the major component of apo(a), the actual function of KIV-2 is still unclear. Here, we describe the first high-resolution crystallographic structure of this domain. It shows a general fold very similar to other KIV domains with high and intermediate affinity for the lysine analog, ε-aminocaproic acid. Interestingly, KIV-2 presents a lysine binding site (LBS) with a unique shape and charge distribution. KIV-2 affinity for predicted small molecule binders was found to be negligible in surface plasmon resonance experiments; and with the LBS being nonfunctional, we propose to rename it "pseudo-LBS". Further investigation of the protein by computational small-molecule docking allowed us to identify a possible heparin-binding site away from the LBS, which was confirmed by specific reverse charge mutations abolishing heparin binding. This study opens new possibilities to define the pathogenesis of Lp(a)-related diseases and to facilitate the design of specific therapeutic drugs.

Computational Drug Repositioning: Current Progress and Challenges

Novel drug discovery is time-consuming, costly, and a high-investment process due to the high attrition rate. Therefore, many trials are conducted to reuse existing drugs to treat pressing conditions and diseases, since their safety profiles and pharmacokinetics are already available. Drug repositioning is a strategy to identify a new indication of existing or already approved drugs, beyond the scope of their original use. Various computational and experimental approaches to incorporate available resources have been suggested for gaining a better understanding of disease mechanisms and the identification of repurposed drug candidates for personalized pharmacotherapy. In this review, we introduce publicly available databases for drug repositioning and summarize the approaches taken for drug repositioning. We also highlight and compare their characteristics and challenges, which should be addressed for the future realization of drug repositioning.

Loop Grafting between Similar Local Environments for Fc-Silent Antibodies

Reduction of the affinity of the fragment crystallizable (Fc) region with immune receptors by substitution of one or a few amino acids, known as Fc-silencing, is an established approach to reduce the immune effector functions of monoclonal antibody therapeutics. This approach to Fc-silencing, however, is problematic as it can lead to instability and immunogenicity of the developed antibodies. We evaluated loop grafting as a novel approach to Fc-silencing in which the Fc loops responsible for immune receptor binding were replaced by loops of up to 20 amino acids from similar local environments in other human and mouse antibodies. Molecular dynamics simulations of the designed variants of an Fc region in a complex with the immune receptor FcγIIIa confirmed that loop grafting potentially leads to a significant reduction in the binding of the antibody variants to the receptor, while retaining their stability. In comparison, standard variants with less than eight substituted amino acids showed possible instability and a lower degree of Fc-silencing due to the occurrence of compensatory interactions. The presented approach to Fc-silencing is general and could be used to modulate undesirable side effects of other antibody therapeutics without affecting their stability or increasing their immunogenicity.

Drug ReposER: a web server for predicting similar amino acid arrangements to known drug binding interfaces for potential drug repositioning

A common drug repositioning strategy is the re-application of an existing drug to address alternative targets. A crucial aspect to enable such repurposing is that the drug's binding site on the original target is similar to that on the alternative target. Based on the assumption that proteins with similar binding sites may bind to similar drugs, the 3D substructure similarity data can be used to identify similar sites in other proteins that are not known targets. The Drug ReposER (DRug REPOSitioning Exploration Resource) web server is designed to identify potential targets for drug repurposing based on sub-structural similarity to the binding interfaces of known drug binding sites. The application has pre-computed amino acid arrangements from protein structures in the Protein Data Bank that are similar to the 3D arrangements of known drug binding sites thus allowing users to explore them as alternative targets. Users can annotate new structures for sites that are similarly arranged to the residues found in known drug binding interfaces. The search results are presented as mappings of matched sidechain superpositions. The results of the searches can be visualized using an integrated NGL viewer. The Drug ReposER server has no access restrictions and is available at http://mfrlab.org/drugreposer/.

In Silico Laboratory: Tools for Similarity-Based Drug Discovery

Computational methods that predict and evaluate binding of ligands to receptors implicated in different pathologies have become crucial in modern drug design and discovery. Here, we describe protocols for using the recently developed package of computational tools for similarity-based drug discovery. The ProBiS stand-alone program and web server allow superimposition of protein structures against large protein databases and predict ligands based on detected binding site similarities. GenProBiS allows mapping of human somatic missense mutations related to cancer and non-synonymous single nucleotide polymorphisms and subsequent visual exploration of specific interactions in connection to these mutations. We describe protocols for using LiSiCA, a fast ligand-based virtual screening software that enables easy screening of large databases containing billions of small molecules. Finally, we show the use of BoBER, a web interface that enables user-friendly access to a large database of bioisosteric and scaffold hopping replacements.

Identification of novel analgesics through a drug repurposing strategy

The identification of new indications for approved or failed drugs is a process called drug repositioning or drug repurposing. The motivation includes overcoming the productivity gap that exists in drug development, which is a high-cost–high-risk process. Repositioning also includes rescuing drugs that have safely entered the market but have failed to demonstrate sufficient efficiency for the initial clinical indication. Considering the high prevalence of chronic pain, the lack of sufficient efficacy and the safety issues of current analgesics, repositioning seems to be an attractive approach. This review presents example of drugs that already have been repositioned and highlights new technologies that are available for the identification of additional compounds to stimulate the curiosity of readers for further exploration.

GenProBiS: web server for mapping of sequence variants to protein binding sites

Discovery of potentially deleterious sequence variants is important and has wide implications for research and generation of new hypotheses in human and veterinary medicine, and drug discovery. The GenProBiS web server maps sequence variants to protein structures from the Protein Data Bank (PDB), and further to protein–protein, protein–nucleic acid, protein–compound, and protein–metal ion binding sites. The concept of a protein–compound binding site is understood in the broadest sense, which includes glycosylation and other post-translational modification sites. Binding sites were defined by local structural comparisons of whole protein structures using the Protein Binding Sites (ProBiS) algorithm and transposition of ligands from the similar binding sites found to the query protein using the ProBiS-ligands approach with new improvements introduced in GenProBiS. Binding site surfaces were generated as three-dimensional grids encompassing the space occupied by predicted ligands. The server allows intuitive visual exploration of comprehensively mapped variants, such as human somatic mis-sense mutations related to cancer and non-synonymous single nucleotide polymorphisms from 21 species, within the predicted binding sites regions for about 80 000 PDB protein structures using fast WebGL graphics. The GenProBiS web server is open and free to all users at http://genprobis.insilab.org.

Discovery of Novel Potential Human Targets of Resveratrol by Inverse Molecular Docking

Resveratrol is a polyphenol known for its antioxidant and anti-inflammatory properties, which support its use as a treatment for variety of diseases. There are already known connections of resveratrol to chemoprevention of cancer because of its ability to prevent tumor initiation and inhibit tumor promotion and progression. Resveratrol is also believed to be important in cardiovascular diseases and neurological disorders, such as Alzheimer's disease. Using an inverse molecular docking approach, we sought to find new potential targets of resveratrol. Docking of resveratrol into each ProBiS predicted binding site of >38 000 protein structures from the Protein Data Bank was examined, and a number of novel potential targets into which resveratrol was docked successfully were found. These explain known actions or predict new effects of resveratrol. The results included three human proteins that are already known to bind resveratrol. A majority of proteins discovered however have no already described connections with resveratrol. We report new potential target human proteins and proteins connected with different organisms into which resveratrol can dock. Our results reveal previously unknown potential target human proteins, whose connection with cardiovascular and neurological disorders could lead to new potential treatments for variety of diseases. We believe that our research could help in future experimental studies on revestratol bioactivity in humans.

Binding site comparisons for target-centered drug discovery

INTRODUCTION

The success of binding site comparisons in drug discovery is based on the recognized fact that many different proteins have similar binding sites. Indeed, binding site comparisons have found many uses in drug development and have the potential to dramatically cut the cost and shorten the time necessary for the development of new drugs. Areas covered: The authors review recent methods for comparing protein binding sites and their use in drug repurposing and polypharmacology. They examine emerging fields including the use of binding site comparisons in precision medicine, the prediction of structured water molecules, the search for targets of natural compounds, and their application in the development of protein-based drugs by loop modeling and for comparison of RNA binding sites. Expert opinion: Binding site comparisons have produced many interesting results in drug development, but relatively little work has been done on protein-protein interaction sites, which are particularly relevant in view of the success of biological drugs. Growth of protein loop modeling for modulating biological drugs is anticipated. The fusion of currently distinct methods for the comparison of RNA and protein binding sites into a single comprehensive approach could allow the search for new selective ribosomal antibiotics and initiate pharmaceutical research into other nucleoproteins.

Functional characterization of β2-adrenergic and insulin receptor heteromers

This study aimed to functionally characterize β₂-adrenergic (β₂AR) and insulin receptor (IR) heteromers in regard to β-arrestin 2 (βarr2) recruitment and cAMP signaling and to examine the involvement of the cytoplasmic portion of the IR β chain in heteromerization with β₂AR. Evidence for β₂AR:IR:βarr2 complex formation and the specificity of the IR:βarr2 interaction was first provided by bioinfomatics analysis. Receptor-heteromer investigation technology (HIT) then provided functional evidence of β₂AR:IR heterodimerization by showing isoproterenol-induced but not insulin-induced GFP²-βarr2 recruitment to the β₂AR:IR complex; the IR:βarr2 interaction was found to only be constitutive. The constitutive IR:βarr2 BRET signal (BRET_const) was significantly smaller in cells coexpressing IR-RLuc8 and a GFP²-βarr2 1-185 mutant lacking the proposed IR binding domain. β₂AR:IR heteromerization also influenced the pharmacological phenotype of β₂AR, i.e., its efficacy in recruiting βarr2 and activating cAMP signaling. Evidence suggesting involvement of the cytoplasmic portion of the IR β chain in the interaction with β₂AR was provided by BRET² saturation and HIT assays using an IR 1-1271 stop mutant lacking the IR C-terminal tail region. For the complex consisting of IR 1-1271-RLuc8:β₂AR-GFP², saturation was not reached, most likely reflecting random collisions between IR 1-1271 and β₂AR. Furthermore, in the HIT assay, no substantial agonist-induced increase in the BRET² signal was detected that would be indicative of βarr2 recruitment to the IR 1-1271:β₂AR heteromer. Complementary 3D visualization of β₂AR:IR provided supporting evidence for stability of the heterotetramer complex and identified amino acid residues involved in β₂AR:IR heteromerization. This article is part of the Special Issue entitled 'Receptor heteromers and their allosteric receptor-receptor interactions'.

Glioblastoma-specific anti-TUFM nanobody for in-vitro immunoimaging and cancer stem cell targeting

Glioblastoma multiforme (GBM) is the most common and lethal form of brain tumor. The prognosis for patients remains poor, despite the combination of new preoperative and intraoperative neuroimaging, radical surgery, and recent advances in radiotherapy and chemotherapy. To improve GBM therapy and patient outcome, sustained drug delivery to glioma cells is needed, while minimizing toxicity to adjacent neurons and glia cells. This might be achieved through an anti-proteomic approach based on nanobodies, the single-domain antigen-binding fragments of heavy-chain antibodies of the camelid adaptive immune system. We report here on the validation and quantification of a nanobody raised against mitochondrial translation elongation factor (TUFM). Differential expression of TUFM was examined in different GBM cell lines and GBM tissue at the protein and mRNA levels, as compared to their expression in neural stem cells and normal brain tissue. We further used in-silico modelling and immunocytochemistry to define the specificity of anti-TUFM nanobody (Nb206) towards GBM stem cells, as compared to GBM cell lines (U251MG and U87MG cells). Due to its specificity and pronounced inhibitory effect on GBM stem cell growth, we propose the use of this anti-TUFM nanobody for GBM in vitro immunoimaging and potentially also cancer stem cell targeting.

Ranking Enzyme Structures in the PDB by Bound Ligand Similarity to Biological Substrates

There are numerous applications that use the structures of protein-ligand complexes from the PDB, such as 3D pharmacophore identification, virtual screening, and fragment-based drug design. The structures underlying these applications are potentially much more informative if they contain biologically relevant bound ligands, with high similarity to the cognate ligands. We present a study of ligand-enzyme complexes that compares the similarity of bound and cognate ligands, enabling the best matches to be identified. We calculate the molecular similarity scores using a method called PARITY (proportion of atoms residing in identical topology), which can conveniently be combined to give a similarity score for all cognate reactants or products in the reaction. Thus, we generate a rank-ordered list of related PDB structures, according to the biological similarity of the ligands bound in the structures.

•

We present PARITY, matching atoms in identical topology to gauge ligand similarity

•

Bound-cognate ligand similarity is a useful metric for ranking PDB structures

•

Only 26% of enzyme structures in the PDB have bound-cognate ligand similarity ≥0.7

•

We provide rank-ordered lists of PDBs with the most biologically relevant ligands

Insights from Ion Binding Site Network Analysis into Evolution and Functions of Proteins

Many biological phenomena can be represented as complex networks. Using a protein binding site comparison approach, we generated a network of ion binding sites on the scale of all known protein structures from the Protein Data Bank. We found that this ion binding site similarity network is scale-free, indicating a network in which a few ion binding site scaffolds are the network hubs, and these are connected to hundreds of nodes, whereas the vast majority of nodes have only a few neighbors. Enrichment and statistical analysis of the network components and communities yielded insights into underlying processes from the functional and the structural perspective. Largest components and communities were observed to be closely related to basic metabolic processes and some of the most common structural folds, which, from the evolutionary point of view, indicates that they may be the oldest ones. Further, we derived the first comprehensive map of ion interchangeability, based on binding site similarity. Several highly interchangeable protein-ion binding site pairs emerged (e.g., Ca²⁺ and Mg²⁺ ), as well as structurally distinct ones. The constructed network of ion binding site similarities will aid in understanding the general principles of protein-ion binding sites structure, function and evolution. We demonstrate potential uses of the network on proteins involved in cancer development and immune response, where individual ions play prominent roles in disease development.

BoBER: web interface to the base of bioisosterically exchangeable replacements

We describe a novel freely available web server Base of Bioisosterically Exchangeable Replacements (BoBER), which implements an interface to a database of bioisosteric and scaffold hopping replacements. Bioisosterism and scaffold hopping are key concepts in drug design and optimization, and can be defined as replacements of biologically active compound's fragments with other fragments to improve activity, reduce toxicity, change bioavailability or to diversify the scaffold space. Our web server enables fast and user-friendly searches for bioisosteric and scaffold replacements which were obtained by mining the whole Protein Data Bank. The working of the web server is presented on an existing MurF inhibitor as example. BoBER web server enables medicinal chemists to quickly search for and get new and unique ideas about possible bioisosteric or scaffold hopping replacements that could be used to improve hit or lead drug-like compounds.

Global organization of a binding site network gives insight into evolution and structure-function relationships of proteins

The global organization of protein binding sites is analyzed by constructing a weighted network of binding sites based on their structural similarities and detecting communities of structurally similar binding sites based on the minimum description length principle. The analysis reveals that there are two central binding site communities that play the roles of the network hubs of smaller peripheral communities. The sizes of communities follow a power-law distribution, which indicates that the binding sites included in larger communities may be older and have been evolutionary structural scaffolds of more recent ones. Structurally similar binding sites in the same community bind to diverse ligands promiscuously and they are also embedded in diverse domain structures. Understanding the general principles of binding site interplay will pave the way for improved drug design and protein design.

ProBiS tools (algorithm, database, and web servers) for predicting and modeling of biologically interesting proteins

ProBiS (Protein Binding Sites) Tools consist of algorithm, database, and web servers for prediction of binding sites and protein ligands based on the detection of structurally similar binding sites in the Protein Data Bank. In this article, we review the operations that ProBiS Tools perform, provide comments on the evolution of the tools, and give some implementation details. We review some of its applications to biologically interesting proteins. ProBiS Tools are freely available at http://probis.cmm.ki.si and http://probis.nih.gov.

Bioinformatics Tools and Resources for Analyzing Protein Structures

The dramatic increase in the number of protein sequences and structures deposited in biological databases has led to the development of many bioinformatics tools and programs to manage, validate, compare, and interpret this large volume of data. In addition, powerful tools are being developed to use this sequence and structural data to facilitate protein classification and infer biological function of newly identified proteins. This chapter covers freely available bioinformatics resources on the World Wide Web that are commonly used for protein structure analysis.

Discovery of Mycobacterium tuberculosis InhA Inhibitors by Binding Sites Comparison and Ligands Prediction

Drug discovery is usually focused on a single protein target; in this process, existing compounds that bind to related proteins are often ignored. We describe ProBiS plugin, extension of our earlier ProBiS-ligands approach, which for a given protein structure allows prediction of its binding sites and, for each binding site, the ligands from similar binding sites in the Protein Data Bank. We developed a new database of precalculated binding site comparisons of about 290000 proteins to allow fast prediction of binding sites in existing proteins. The plugin enables advanced viewing of predicted binding sites, ligands' poses, and their interactions in three-dimensional graphics. Using the InhA query protein, an enoyl reductase enzyme in the Mycobacterium tuberculosis fatty acid biosynthesis pathway, we predicted its possible ligands and assessed their inhibitory activity experimentally. This resulted in three previously unrecognized inhibitors with novel scaffolds, demonstrating the plugin's utility in the early drug discovery process.

Computational identification of non-synonymous polymorphisms within regions corresponding to protein interaction sites

BACKGROUND

Protein-protein interactions (PPI) play an important role in function of all organisms and enable understanding of underlying metabolic processes. Computational predictions of PPIs are an important aspect in proteomics, as experimental methods may result in high degree of false positive results and are more expensive. Although there are many databases collecting predicted PPIs, exploration of genetics information underlying PPI interactions has not been investigated thoroughly. The aim of the present study was to identify genomic locations corresponding to regions involved in predicted PPIs and to collect non-synonymous polymorphisms (nsSNPs) located within those regions; which we termed PPI-SNPs.

METHODS

Predicted PPIs were obtained from PiSITE database (http://pisite.hgc.jp). Non-synonymous SNPs mapped on protein structural data (PDBs) were obtained from the UCSC server. Polymorphism locations on protein structures were mapped to predicted PPI regions. DAVID tool was used for pathway enrichment and gene cluster analysis (https://david.ncifcrf.gov/).

RESULTS

We collected 544 polymorphisms located within predicted PPI sites that map to 197 genes. We identified 9 SNPs, previously associated with diseases, but not yet associated with PPI sites. We also found examples in which polymorphisms located within predicted PPI regions are also occurring within previously experimentally validated PPIs and within experimentally determined functional domains.

CONCLUSIONS

Our study provides the first catalog of nsSNPs located within predicted PPIs. These prioritized SNPs present the basis for planning experimental validation of SNPs that cause gain or loss of PPIs. Our implementation is expandable, as datasets used are constantly updated.

Protein ligand-specific binding residue predictions by an ensemble classifier

BACKGROUND

Prediction of ligand binding sites is important to elucidate protein functions and is helpful for drug design. Although much progress has been made, many challenges still need to be addressed. Prediction methods need to be carefully developed to account for chemical and structural differences between ligands.

RESULTS

In this study, we present ligand-specific methods to predict the binding sites of protein-ligand interactions. First, a sequence-based method is proposed that only extracts features from protein sequence information, including evolutionary conservation scores and predicted structure properties. An improved AdaBoost algorithm is applied to address the serious imbalance problem between the binding and non-binding residues. Then, a combined method is proposed that combines the current template-free method and four other well-established template-based methods. The above two methods predict the ligand binding sites along the sequences using a ligand-specific strategy that contains metal ions, acid radical ions, nucleotides and ferroheme. Testing on a well-established dataset showed that the proposed sequence-based method outperformed the profile-based method by 4-19% in terms of the Matthews correlation coefficient on different ligands. The combined method outperformed each of the individual methods, with an improvement in the average Matthews correlation coefficients of 5.55% over all ligands. The results also show that the ligand-specific methods significantly outperform the general-purpose methods, which confirms the necessity of developing elaborate ligand-specific methods for ligand binding site prediction.

CONCLUSIONS

Two efficient ligand-specific binding site predictors are presented. The standalone package is freely available for academic usage at http://dase.ecnu.edu.cn/qwdong/TargetCom/TargetCom_standalone.tar.gz  or request upon the corresponding author.

Comparison of in silico tools for binding site prediction applied for structure-based design of autolysin inhibitors

Autolysin E (AtlE) is a bacteriolytic enzyme which plays an important role in division and growth of bacterial cells and therefore represents a promising potential drug target. Its 3D structure has been recently elucidated. We used in silico prediction tools to study substrate or ligand (inhibitor) binding regions of AtlE. We applied several freely available tools and a commercial tool for binding site identification and compared results of the prediction. Calculation time, number of predictions and output data provided by specific software vary according to the different approaches utilized by specific method categories. Despite different approaches, binding sites in similar locations on the protein were predicted. Specific amino acid residues that form these binding sites were predicted as binding residues. The predicted residues, especially those with predicted highest conservation score, could theoretically have catalytic and binding properties. According to our results, we assume that E138, which has the highest conservation score, is the catalytic residue and F161, G162 and Y224, which are also highly conserved, are responsible for substrate binding. Ligands developed with binding specificity towards these residues could inhibit the catalysis and binding of the substrate of AtlE. The molecules with inhibitory potency could therefore represent potential new antibacterial agents.

Identification of Ecdysone Hormone Receptor Agonists as a Therapeutic Approach for Treating Filarial Infections

Background

A homologue of the ecdysone receptor has previously been identified in human filarial parasites. As the ecdysone receptor is not found in vertebrates, it and the regulatory pathways it controls represent attractive potential chemotherapeutic targets.

Methodology/ Principal Findings

Administration of 20-hydroxyecdysone to gerbils infected with B. malayi infective larvae disrupted their development to adult stage parasites. A stable mammalian cell line was created incorporating the B. malayi ecdysone receptor ligand-binding domain, its heterodimer partner and a secreted luciferase reporter in HEK293 cells. This was employed to screen a series of ecdysone agonist, identifying seven agonists active at sub-micromolar concentrations. A B. malayi ecdysone receptor ligand-binding domain was developed and used to study the ligand-receptor interactions of these agonists. An excellent correlation between the virtual screening results and the screening assay was observed. Based on both of these approaches, steroidal ecdysone agonists and the diacylhydrazine family of compounds were identified as a fruitful source of potential receptor agonists. In further confirmation of the modeling and screening results, Ponasterone A and Muristerone A, two compounds predicted to be strong ecdysone agonists stimulated expulsion of microfilaria and immature stages from adult parasites.

Conclusions

The studies validate the potential of the B. malayi ecdysone receptor as a drug target and provide a means to rapidly evaluate compounds for development of a new class of drugs against the human filarial parasites.

The human filarial parasites are the causative agents of two neglected tropical diseases targeted for elimination by the international community. The current elimination programs rely upon the mass distribution of a limited number of drugs, leaving the programs open to failure in the event that resistance develops. Thus, there is a critical need for novel chemotherapeutic agents to supplement the current arsenal. The filarial parasites are ecdysozoans, whose developmental processes are controlled by a master regulator, the ecdysone receptor. Here we validate the potential of the filarial ecdysone receptor as a chemotherapeutic target and report the development of high throughput and virtual screening assays that may be used to compounds that target it.

Minimal Functional Sites in Metalloproteins and Their Usage in Structural Bioinformatics

Metal ions play a functional role in numerous biochemical processes and cellular pathways. Indeed, about 40% of all enzymes of known 3D structure require a metal ion to be able to perform catalysis. The interactions of the metals with the macromolecular framework determine their chemical properties and reactivity. The relevant interactions involve both the coordination sphere of the metal ion and the more distant interactions of the so-called second sphere, i.e., the non-bonded interactions between the macromolecule and the residues coordinating the metal (metal ligands). The metal ligands and the residues in their close spatial proximity define what we call a minimal functional site (MFS). MFSs can be automatically extracted from the 3D structures of metal-binding biological macromolecules deposited in the Protein Data Bank (PDB). They are 3D templates that describe the local environment around a metal ion or metal cofactor and do not depend on the overall macromolecular structure. MFSs provide a different view on metal-binding proteins and nucleic acids, completely focused on the metal. Here we present different protocols and tools based upon the concept of MFS to obtain deeper insight into the structural and functional properties of metal-binding macromolecules. We also show that structure conservation of MFSs in metalloproteins relates to local sequence similarity more strongly than to overall protein similarity.

Mechanistic Study of Human Glucose Transport Mediated by GLUT1

The glucose transporter 1 (GLUT1) belongs to the major facilitator superfamily (MFS) and is responsible for the constant uptake of glucose. However, the molecular mechanism of sugar transport remains obscure. In this study, homology modeling and molecular dynamics (MD) simulations in lipid bilayers were performed to investigate the combination of the alternate and multisite transport mechanism of glucose with GLUT1 in atomic detail. To explore the substrate recognition mechanism, the outward-open state human GLUT1 homology model was generated based on the template of xylose transporter XylE (PDB ID: 4GBZ), which shares up to 29% sequence identity and 49% similarity with GLUT1. Through the MD simulation study of glucose across lipid bilayer with both the outward-open GLUT1 and the GLUT1 inward-open crystal structure, we investigated six different conformational states and identified four key binding sites in both exofacial and endofacial loops that are essential for glucose recognition and transport. The study further revealed that four flexible gates consisting of W65/Y292/Y293-M420/TM10b-W388 might play important roles in the transport cycle. The study showed that some side chains close to the central ligand binding site underwent larger position changes. These conformational interchanges formed gated networks within an S-shaped central channel that permitted staged ligand diffusion across the transporter. This study provides new inroads for the understanding of GLUT1 ligand recognition paradigm and configurational features which are important for molecular, structural, and physiological research of the MFS members, especially for GLUT1-targeted drug design and discovery.