Allosteric conversation in the androgen receptor ligand-binding domain surfaces

Androgen receptor (AR) is a major therapeutic target that plays pivotal roles in prostate cancer (PCa) and androgen insensitivity syndromes. We previously proposed that compounds recruited to ligand-binding domain (LBD) surfaces could regulate AR activity in hormone-refractory PCa and discovered several surface modulators of AR function. Surprisingly, the most effective compounds bound preferentially to a surface of unknown function [binding function 3 (BF-3)] instead of the coactivator-binding site [activation function 2 (AF-2)]. Different BF-3 mutations have been identified in PCa or androgen insensitivity syndrome patients, and they can strongly affect AR activity. Further, comparison of AR x-ray structures with and without bound ligands at BF-3 and AF-2 showed structural coupling between both pockets. Here, we combine experimental evidence and molecular dynamic simulations to investigate whether BF-3 mutations affect AR LBD function and dynamics possibly via allosteric conversation between surface sites. Our data indicate that AF-2 conformation is indeed closely coupled to BF-3 and provide mechanistic proof of their structural interconnection. BF-3 mutations may function as allosteric elicitors, probably shifting the AR LBD conformational ensemble toward conformations that alter AF-2 propensity to reorganize into subpockets that accommodate N-terminal domain and coactivator peptides. The induced conformation may result in either increased or decreased AR activity. Activating BF-3 mutations also favor the formation of another pocket (BF-4) in the vicinity of AF-2 and BF-3, which we also previously identified as a hot spot for a small compound. We discuss the possibility that BF-3 may be a protein-docking site that binds to the N-terminal domain and corepressors. AR surface sites are attractive pharmacological targets to develop allosteric modulators that might be alternative lead compounds for drug design.

Prediction of 3D structure and ligand-binding properties : ABC transporters and flavodiiron proteins from Synechocystis sp. strain PCC 6803

Cyanobacteria are unicellular, non-nitrogen-fixing prokaryotes, which perform photosynthesis similarly as higher plants. The cyanobacterium Synechocystis sp. strain PCC 6803 is used as a model organism in photosynthesis research. My research described herein aims at understanding the function of the photosynthetic machinery and how it responds to changes in the environment. Detailed knowledge of the regulation of photosynthesis in cyanobacteria can be utilized for biotechnological purposes, for example in the harnessing of solar energy for biofuel production. In photosynthesis, iron participates in electron transfer. Here, we focused on iron transport in Synechocystis sp. strain PCC 6803 and particularly on the environmental regulation of the genes encoding the FutA2BC ferric iron transporter, which belongs to the ABC transporter family. A homology model built for the ATP-binding subunit FutC indicates that it has a functional ATPbinding site as well as conserved interactions with the channel-forming subunit FutB in the transporter complex. Polyamines are important for the cell proliferation, differentiation and apoptosis in prokaryotic and eukaryotic cells. In plants, polyamines have special roles in stress response and in plant survival. The polyamine metabolism in cyanobacteria in response to environmental stress is of interest in research on stress tolerance of higher plants. In this thesis, the potd gene encoding an polyamine transporter subunit from Synechocystis sp. strain PCC 6803 was characterized for the first time. A homology model built for PotD protein indicated that it has capability of binding polyamines, with the preference for spermidine. Furthermore, in order to investigate the structural features of the substrate specificity, polyamines were docked into the binding site. Spermidine was positioned very similarly in Synechocystis PotD as in the template structure and had most favorable interactions of the docked polyamines. Based on the homology model, experimental work was conducted, which confirmed the binding preference. Flavodiiron proteins (Flv) are enzymes, which protect the cell against toxicity of oxygen and/or nitric oxide by reduction. In this thesis, we present a novel type of photoprotection mechanism in cyanobacteria by the heterodimer of Flv2/Flv4. The constructed homology model of Flv2/Flv4 suggests a functional heterodimer capable of rapid electron transfer. The unknown protein sll0218, encoded by the flv2-flv4 operon, is assumed to facilitate the interaction of the Flv2/Flv4 heterodimer and energy transfer between the phycobilisome and PSII. Flv2/Flv4 provides an alternative electron transfer pathway and functions as an electron sink in PSII electron transfer.

Binding sites and actions of Tx1, a neurotoxin from the venom of the spider Phoneutria nigriventer, in guinea pig ileum

Tx1, a neurotoxin isolated from the venom of the South American spider Phoneutria nigriventer, produces tail elevation, behavioral excitation and spastic paralysis of the hind limbs after intracerebroventricular injection in mice. Since Tx1 contracts isolated guinea pig ileum, we have investigated the effect of this toxin on acetylcholine release, as well as its binding to myenteric plexus-longitudinal muscle membranes from the guinea pig ileum. [125I]-Tx1 binds specifically and with high affinity (Kd = 0.36 ± 0.02 nM) to a single, non-interacting (nH = 1.1), low capacity (Bmax 1.1 pmol/mg protein) binding site. In competition experiments using several compounds (including ion channel ligands), only PhTx2 and PhTx3 competed with [125I]-Tx1 for specific binding sites (K0.5 apparent = 7.50 x 10-4 g/l and 1.85 x 10-5 g/l, respectively). PhTx2 and PhTx3, fractions from P. nigriventer venom, contain toxins acting on sodium and calcium channels, respectively. However, the neurotoxin PhTx2-6, one of the isoforms found in the PhTx2 pool, did not affect [125I]-Tx1 binding. Tx1 reduced the [3H]-ACh release evoked by the PhTx2 pool by 33%, but did not affect basal or KCl-induced [3H]-ACh release. Based on these results, as well as on the homology of Tx1 with toxins acting on calcium channels (w-Aga IA and IB) and its competition with [125I]-w-Cono GVIA in the central nervous system, we suggest that the target site for Tx1 may be calcium channels.

Studies on ATP-diphosphohydrolase nucleotide-binding sites by intrinsic fluorescence

Potato apyrase, a soluble ATP-diphosphohydrolase, was purified to homogeneity from several clonal varieties of Solanum tuberosum. Depending on the source of the enzyme, differences in kinetic and physicochemical properties have been described, which cannot be explained by the amino acid residues present in the active site. In order to understand the different kinetic behavior of the Pimpernel (ATPase/ADPase = 10) and Desirée (ATPase/ADPase = 1) isoenzymes, the nucleotide-binding site of these apyrases was explored using the intrinsic fluorescence of tryptophan. The intrinsic fluorescence of the two apyrases was slightly different. The maximum emission wavelengths of the Desirée and Pimpernel enzymes were 336 and 340 nm, respectively, suggesting small differences in the microenvironment of Trp residues. The Pimpernel enzyme emitted more fluorescence than the Desirée apyrase at the same concentration although both enzymes have the same number of Trp residues. The binding of the nonhydrolyzable substrate analogs decreased the fluorescence emission of both apyrases, indicating the presence of conformational changes in the neighborhood of Trp residues. Experiments with quenchers of different polarities, such as acrylamide, Cs+ and I- indicated the existence of differences in the nucleotide-binding site, as further shown by quenching experiments in the presence of nonhydrolyzable substrate analogs. Differences in the nucleotide-binding site may explain, at least in part, the kinetic differences of the Pimpernel and Desirée isoapyrases.

Modulation allostérique de la fonction des récepteurs FP et PAF

Les récepteurs couplés aux protéines-G (RCPGs) constituent la première étape d’une série de cascades signalétiques menant à la régulation d’une multitude de processus physiologiques. Dans le modèle classique connu, la liaison du ligand induit un changement de conformation du récepteur qui mène à sa forme active. Une fois activés, les RCPGs vont réguler l’activité d’une protéine membranaire cible qui peut être tant une enzyme qu’un canal ionique. L’interaction entre le récepteur et la cible nécessite l’intermédiaire d’une protéine hétérotrimérique appelée « protéine G », qui est activée pour favoriser l’échange du GDP (guanosine diphosphate) pour un GTP (guanosine triphosphate) et assurer la transduction du signal du récepteur à l’effecteur. Les mécanismes moléculaires menant à l’activation des effecteurs spécifiques via l’activation des RCPGs par les protéines G hétérotrimériques sont encore plutôt méconnus. Dans notre étude nous nous sommes intéressés aux récepteurs FP et PAF, à leurs ligands naturels, la PGF2α et le Carbamyl-PAF respectivement, et à des ligands à action antagoniste sur ces récepteurs. Des ligands considérés comme agonistes, sont des molécules qui interagissent avec le récepteur et induisent les mêmes effets que le ligand naturel. Les antagonistes, par contre, sont des molécules qui interagissent avec le récepteur et bloquent l’action du ligand naturel en prévenant le changement conformationnel du complexe, et ils peuvent avoir une action compétitive ou non-compétitive. Nous avons étudié aussi des ligands orthostériques et allostériques du récepteur FP des prostaglandines et du récepteur PAF. Un ligand orthostérique peut se comporter comme agoniste ou antagoniste en se fixant au site de liaison du ligand (agoniste) naturel. Un ligand allostérique est un agoniste ou antagoniste se fixant à un site autre que celui du ligand naturel entraînant un changement de conformation ayant pour conséquence soit une augmentation (effecteur positif), soit une diminution (effecteur négatif) de l'activité du ligand naturel.

Caractérisation moléculaire du rôle de facteurs accessoires ArgR et PepA au niveau de la recombinaison spécifique sur le site cer

Mon projet de recherche avait pour but de caractériser le rôle de deux protéines, ArgR et PepA, qui agissent en tant que facteurs accessoires de la recombinaison au niveau de deux sites cer du plasmide ColE1 présent dans la bactérie Escherichia coli. Ces deux protéines, couplées aux deux recombinases à tyrosine XerC et XerD, permettent la catalyse de la recombinaison site spécifique au niveau de la séquence cer, convertissant les multimères instables de ColE1 en monomères stables. Cette étude a principalement porté sur la région C-terminale de la protéine ArgR. Cette région de la protéine ArgR possède une séquence en acides-aminés et une structure similaire à celle de la protéine AhrC de Bacillus subtilis. De plus, AhrC, le répresseur de l’arginine de cette bactérie, est capable de complémenter des Escherichia coli mutantes déficientes en ArgR. Les régions C-terminales de ces protéines, montrent une forte similarité. De précédents travaux dans notre laboratoire ont démontré que des mutants d’ArgR comprenant des mutations dans cette région, en particulier les mutants ArgR149, une version tronquée d’ArgR de 149 acides-aminés, et ArgR5aa, une version comprenant une insertion de cinq acides-aminés dans la partie C-terminale, perdaient la capacité de permettre la recombinaison au niveau de deux sites cer présents dans le plasmide pCS210. Malgré cette incapacité à promouvoir la réaction de recombinaison en cer, ces deux mutants étaient toujours capables de se lier spécifiquement à l’ADN et de réprimer une fusion argA :: lacZ. Dans ce travail, les versions mutantes et sauvages d’ArgR furent surexprimées en tant que protéines de fusion 6-histidine. Des analyses crosslinking ont montré que la version sauvage et ArgR5aa pouvaient former des hexamères in-vitro de manière efficace, alors qu’ArgR149 formait des multimères de plus faible poids moléculaire. Des formes tronquées d’ArgR qui comportaient 150 acides-aminés ou plus, étaient encore capables de permettre la recombinaison en cer. Les mutants par substitution ArgRL149A et ArgRL151A ont tous montré que les substitutions d’un seul acide-aminé au sein de cette région avaient peu d’effets sur la recombinaison en cer. Les expériences de crosslinking protéine-à-protéine ont montré que le type sauvage et les formes mutantes d’ArgR étaient capables d’interagir avec la protéine accessoire PepA, également impliquée dans la recombinaison en cer. Les expériences de recombinaison in-vitro utilisant la forme sauvage et les formes mutantes d’ArgR combinées avec les protéines PepA, XerC et XerD purifiées, ont montré que le mutant ArgR149 ne soutenait pas la recombinaison, mais que le mutant ArgR5aa permettait la formation d’une jonction d’Holliday. Des expériences de topologie ont montré que PepA était capable de protéger l’ADN de la topoisomérase 1, et d’empêcher ArgRWt de se lier à l’ADN. Les deux mutants ArgR149 et ArgR5aa protègent aussi l’ADN avec plus de surenroulements. Quand on ajoute PepA, les profils de migration montrent un problème de liaison des deux mutants avec PepA. D’autres expériences impliquant le triplet LEL (leucine-acide glutamique-leucine) et les acides-aminés alentour devraient être réalisés dans le but de connaitre l’existence d’un site de liaison potentiel pour PepA.

Destins des S-RNases et interactions moléculaires dans le tube pollinique dans le cadre de l’auto-incompatibilité gamétophytique chez Solanum chacoense

L’auto-incompatibilité (AI) est une barrière reproductive prézygotique qui permet aux pistils d’une fleur de rejeter leur propre pollen. Les systèmes d’AI peuvent prévenir l’autofertilisation et ainsi limiter l’inbreeding. Dans l’AI gamétophytique, le génotype du pollen détermine son propre phénotype d’incompatibilité, et dans ce système, les déterminants mâles et femelles de l’AI sont codés par un locus multigénique et multi-allélique désigné le locus S. Chez les Solanaceae, le déterminant femelle de l’AI est une glycoprotéine stylaire extracellulaire fortement polymorphique possédant une activité ribonucléase et désignée S-RNase. Les S-RNases montrent un patron caractéristique de deux régions hypervariables (HVa et HVb), responsables de leur détermination allélique, et cinq régions hautement conservées (C1 à C5) impliquées dans l’activité catalytique ou la stabilisation structurelle de ces protéines. Dans ce travail, nous avons investigué plusieurs caractéristiques des S-RNases et identifié un nouveau ligand potentiel aux S-RNases chez Solanum chacoense. L’objectif de notre première étude était l’élucidation du rôle de la région C4 des S-RNases. Afin de tester l’hypothèse selon laquelle la région C4 serait impliquée dans le repliement ou la stabilité des S-RNases, nous avons généré un mutant dans lequel les quatre résidus chargés présents en région C4 furent remplacés par des résidus glycine. Cette protéine mutante ne s’accumulant pas à des niveaux détectables, la région C4 semble bien avoir un rôle structurel. Afin de vérifier si C4 est impliquée dans une liaison avec une autre protéine, nous avons généré le mutant R115G, dans lequel un acide aminé chargé fût éliminé afin de réduire les affinités de liaison dans cette région. Ce mutant n’affectant pas le phénotype de rejet pollinique, il est peu probable que la région C4 soit impliquée dans la liaison des S-RNases avec un ligand ou leur pénétration à l’intérieur des tubes polliniques. Enfin, le mutant K113R, dans lequel le seul résidu lysine conservé parmi toutes les S-RNases fût remplacé par un résidu arginine, fût généré afin de vérifier si cette lysine était un site potentiel d’ubiquitination des S-RNases. Toutefois, la dégradation des S-RNases ne fût pas inhibée. Ces résultats indiquent que C4 joue probablement un rôle structurel de stabilisation des S-RNases. Dans une seconde étude, nous avons analysé le rôle de la glycosylation des S-RNases, dont un site, en région C2, est conservé parmi toutes les S-RNases. Afin d’évaluer la possibilité que les sucres conjugués constituent une cible potentielle d’ubiquitination, nous avons généré une S11-RNase dont l‘unique site de glycosylation en C2 fût éliminé. Ce mutant se comporte de manière semblable à une S11-RNase de type sauvage, démontrant que l’absence de glycosylation ne confère pas un phénotype de rejet constitutif du pollen. Afin de déterminer si l’introduction d’un sucre dans la région HVa de la S11-RNase pourrait affecter le rejet pollinique, nous avons généré un second mutant comportant un site additionnel de glycosylation dans la région HVa et une troisième construction qui comporte elle aussi ce nouveau site mais dont le site en région C2 fût éliminé. Le mutant comportant deux sites de glycosylation se comporte de manière semblable à une S11-RNase de type sauvage mais, de manière surprenante, le mutant uniquement glycosylé en région HVa peut aussi rejeter le pollen d’haplotype S13. Nous proposons que la forme non glycosylée de ce mutant constitue un allèle à double spécificité, semblable à un autre allèle à double spécificité préalablement décrit. Il est intéressant de noter que puisque ce phénotype n’est pas observé dans le mutant comportant deux sites de glycosylation, cela suggère que les S-RNases ne sont pas déglycosylées à l’intérieur du pollen. Dans la dernière étude, nous avons réalisé plusieurs expériences d’interactions protéine-protéine afin d’identifier de potentiels interactants polliniques avec les S-RNases. Nous avons démontré que eEF1A, un composant de la machinerie de traduction chez les eucaryotes, peut lier une S11-RNase immobilisée sur résine concanavaline A. Des analyses de type pull-down utilisant la protéine eEF1A de S. chacoense étiquetée avec GST confirment cette interaction. Nous avons aussi montré que la liaison, préalablement constatée, entre eEF1A et l’actine est stimulée en présence de la S11-RNase, bien que cette dernière ne puisse directement lier l’actine. Enfin, nous avons constaté que dans les tubes polliniques incompatibles, l’actine adopte une structure agrégée qui co-localise avec les S-RNases. Ces résultats suggèrent que la liaison entre eEF1A et les S-RNases pourrait constituer un potentiel lien fonctionnel entre les S-RNases et l’altération du cytosquelette d’actine observée lors des réactions d’AI. Par ailleurs, si cette liaison est en mesure de titrer les S-RNases disponibles à l’intérieur du tube pollinique, ce mécanisme pourrait expliquer pourquoi des quantités minimales ou « seuils » de S-RNases sont nécessaires au déclenchement des réactions d’AI.

Nuclear magnetic resonance structure of the nucleic acid-binding domain of severe acute respiratory syndrome coronavirus nonstructural protein 3

The nuclear magnetic resonance (NMR) structure of a globular domain of residues 1071 to 1178 within the previously annotated nucleic acid-binding region (NAB) of severe acute respiratory syndrome coronavirus nonstructural protein 3 (nsp3) has been determined, and N- and C-terminally adjoining polypeptide segments of 37 and 25 residues, respectively, have been shown to form flexibly extended linkers to the preceding globular domain and to the following, as yet uncharacterized domain. This extension of the structural coverage of nsp3 was obtained from NMR studies with an nsp3 construct comprising residues 1066 to 1181 [ nsp3(1066-1181)] and the constructs nsp3(1066-1203) and nsp3(1035-1181). A search of the protein structure database indicates that the globular domain of the NAB represents a new fold, with a parallel four-strand beta-sheet holding two alpha-helices of three and four turns that are oriented antiparallel to the beta-strands. Two antiparallel two-strand beta-sheets and two 3(10)-helices are anchored against the surface of this barrel-like molecular core. Chemical shift changes upon the addition of single-stranded RNAs (ssRNAs) identified a group of residues that form a positively charged patch on the protein surface as the binding site responsible for the previously reported affinity for nucleic acids. This binding site is similar to the ssRNA-binding site of the sterile alpha motif domain of the Saccharomyces cerevisiae Vts1p protein, although the two proteins do not share a common globular fold.

Thermodynamic analysis of ferrous ion binding to Escherichia coli ferritin EcFtnA

Iron oxidation in the bacterial ferritin EcFtnA from Escherichia coli shows marked differences from its homologue human H-chain ferritin (HuHF). While the amino acid residues that constitute the dinuclear center in these proteins are highly conserved, EcFtnA has a third iron-binding site (C site) in close proximity to the dinuclear center that is seemingly responsible for these differences. Here, we describe the first thermodynamic study of Fe2+ binding to EcFtnA and its variants to determine the location of the primary ferrous ion-binding sites on the protein and to better understand the role of the third C site in iron binding. Isothermal titration calorimetric analyses of the wild-type protein reveal the presence of two main classes of binding sites in the pH range of 6.5-7.5, ascribed to Fe2+ binding, first at the A and then the B sites. Site-directed mutagenesis of ligands in the A, B, or C sites affects the apparent Fe2+-binding stoichiometries at the unaltered sites. The data imply some degree of inter- and intrasubunit negative cooperative interaction between sites. Unlike HuHF where only the A site initially binds Fe2+, both A and B sites in EcFtnA bind Fe2+, implying a role for the C site in influencing the binding of Fe2+ at the B site of the di-iron center of EcFtnA. The ITC equations describing a binding model for three classes of independent binding sites are reported here for the first time.

Mass spectroscopic characterization of the coronavirus infectious bronchitis virus nucleoprotein and elucidation of the role of phosphorylation in RNA binding by using surface plasmon resonance

Phosphorylation of the coronavirus nucleoprotein (N protein) has been predicted to play a role in RNA binding. To investigate this hypothesis, we examined the kinetics of RNA binding between nonphosphorylated and phosphorylated infectious bronchitis virus N protein with nonviral and viral RNA by surface plasmon resonance (Biacore). Mass spectroscopic analysis of N protein identified phosphorylation sites that were proximal to RNA binding domains. Kinetic analysis, by surface plasmon resonance, indicated that nonphospborylated N protein bound with the same affinity to viral RNA as phosphorylated N protein. However, phosphorylated N protein bound to viral RNA with a higher binding affinity than nonviral RNA, suggesting that phosphorylation of N protein determined the recognition of virus RNA. The data also indicated that a known N protein binding site (involved in transcriptional regulation) consisting of a conserved core sequence present near the 5' end of the genome (in the leader sequence) functioned by promoting high association rates of N protein binding. Further analysis of the leader sequence indicated that the core element was not the only binding site for N protein and that other regions functioned to promote high-affinity binding.

Oligomers of D-2 dopamine receptors - Evidence from ligand binding

There is increasing evidence that G protein-coupled receptors form oligomers and that this might be important for their function. We have studied this phenomenon for the D-2 dopamine receptor and have shown-using a variety of biochemical and biophysical techniques-that this receptor forms dimers or higher-order oligomers. Using ligand-binding studies, we have also found evidence that this oligomer formation has functional relevance. Thus, for the receptor expressed in either CHO cells or Sf 9 insect cells, the binding properties of several radioligands (in saturation, competition, and dissociation assays) do not conform to those expected for a monomeric receptor with a single binding site. We propose that the receptors exist in oligomers with homotropic and heterotropic negatively cooperative interactions between ligands.

Insights into the effects on metal binding of the systematic substitution of five key glutamate ligands in the ferritin of Escherichia coli

Ferritins are nearly ubiquitous iron storage proteins playing a fundamental role in iron metabolism. They are composed of 24 subunits forming a spherical protein shell encompassing a central iron storage cavity. The iron storage mechanism involves the initial binding and subsequent O-2-dependent oxidation of two Fe2+ ions located at sites A and B within the highly conserved dinuclear "ferroxidase center" in individual subunits. Unlike animal ferritins and the heme-containing bacterioferritins, the Escherichia coli ferritin possesses an additional iron-binding site (site C) located on the inner surface of the protein shell close to the ferroxidase center. We report the structures of five E. coli ferritin variants and their Fe3+ and Zn2+ (a redox-stable alternative for Fe2+) derivatives. Single carboxyl ligand replacements in sites A, B, and C gave unique effects on metal binding, which explain the observed changes in Fe2+ oxidation rates. Binding of Fe2+ at both A and B sites is clearly essential for rapid Fe2+ oxidation, and the linking of Fe-B(2+) to Fe-C(2+) enables the oxidation of three Fe2+ ions. The transient binding of Fe2+ at one of three newly observed Zn2+ sites may allow the oxidation of four Fe2+ by one dioxygen molecule.

Oligomers of D-2 dopamine receptors: evidence from ligand binding

There is increasing evidence that G protein-coupled receptors form oligomers and that this might be important for their function. We have studied this phenomenon for the D-2 dopamine receptor and have shown-using a variety of biochemical and biophysical techniques-that this receptor forms dimers or higher-order oligomers. Using ligand-binding studies, we have also found evidence that this oligomer formation has functional relevance. Thus, for the receptor expressed in either CHO cells or Sf 9 insect cells, the binding properties of several radioligands (in saturation, competition, and dissociation assays) do not conform to those expected for a monomeric receptor with a single binding site. We propose that the receptors exist in oligomers with homotropic and heterotropic negatively cooperative interactions between ligands

Titanocene anticancer complexes and their binding mode of action to human serum albumin: a computational study

Due to the pivotal role played by human serum albumin (HSA) in the transport and cytotoxicity of titanocene complexes, a docking study has been performed on a selected set of titanocene complexes to aid in the current understanding of the potential mode of action of these titanocenes upon binding HSA. Analysis of the docking results has revealed potential binding at the known drug binding sites in HSA and has provided some explanation for the specificity and subsequent cytotoxicity of these titanocenes. Additionally, a new alternative binding site for these titanocenes has been postulated.

In silico identification and characterization of protein-ligand binding sites

Protein–ligand binding site prediction methods aim to predict, from amino acid sequence, protein–ligand interactions, putative ligands, and ligand binding site residues using either sequence information, structural information, or a combination of both. In silico characterization of protein–ligand interactions has become extremely important to help determine a protein’s functionality, as in vivo-based functional elucidation is unable to keep pace with the current growth of sequence databases. Additionally, in vitro biochemical functional elucidation is time-consuming, costly, and may not be feasible for large-scale analysis, such as drug discovery. Thus, in silico prediction of protein–ligand interactions must be utilized to aid in functional elucidation. Here, we briefly discuss protein function prediction, prediction of protein–ligand interactions, the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated EvaluatiOn (CAMEO) competitions, along with their role in shaping the field. We also discuss, in detail, our cutting-edge web-server method, FunFOLD for the structurally informed prediction of protein–ligand interactions. Furthermore, we provide a step-by-step guide on using the FunFOLD web server and FunFOLD3 downloadable application, along with some real world examples, where the FunFOLD methods have been used to aid functional elucidation.