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.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

The effect of lysobisphosphatidic acid (LBPA) on α-tocopherol transfer protein (α-TTP) binding to lipid membranes

The a-tocopherol transfer protein (a-TTP) is responsible for the retention of the atocopherol form of vitamin E in living organisms. The detailed ligand transfer mechanism by a-TTP is still yet to be fully elucidated. To date, studies show that a-TTP transfers a-tocopherol from late endosomes in liver cells to the plasma membrane where it is repackaged into very low density lipoprotein (VLDL) and released into the circulation. Late endosomes have been shown to contain a lipid known as lysobisphosphatidic acid (LBP A) that is unique to this cellular compartment. LBPA plays a role in intracellular trafficking and controlling membrane curvature. Taking these observations into account plus the fact that certain proteins are recruited to membranes based on membrane curvature, the specific aim of this project was to examine the effect of LBP A on a-TTP binding to lipid membranes. To achieve this objective, dual polarization interferometry (DPI) and a vesicle binding assay were employed. Whilst DPI allows protein binding affinity to be measured on a flat lipid surface, the vesicle binding assay determines protein binding affinity to lipid vesicles mimicking curved membranes. DPI analysis revealed that the amount of a-TTP bound to lipid membranes is higher when LBPA is present. Using the vesicle binding assay, a similar result was seen where a greater amount of protein is bound to large unilamellar vesicles (LUV s) containing LBP A. However, the effect of LBP A was attenuated when small unilamellar vesicles (SUVs) were replaced with LUVs. The outcome of this project suggests that aTTP binding to membranes is influenced by membrane curvature, which in turn is induced by the presence of LBP A.

FunFOLDQA: a quality assessment tool for protein-ligand binding site residue predictions

The estimation of prediction quality is important because without quality measures, it is difficult to determine the usefulness of a prediction. Currently, methods for ligand binding site residue predictions are assessed in the function prediction category of the biennial Critical Assessment of Techniques for Protein Structure Prediction (CASP) experiment, utilizing the Matthews Correlation Coefficient (MCC) and Binding-site Distance Test (BDT) metrics. However, the assessment of ligand binding site predictions using such metrics requires the availability of solved structures with bound ligands. Thus, we have developed a ligand binding site quality assessment tool, FunFOLDQA, which utilizes protein feature analysis to predict ligand binding site quality prior to the experimental solution of the protein structures and their ligand interactions. The FunFOLDQA feature scores were combined using: simple linear combinations, multiple linear regression and a neural network. The neural network produced significantly better results for correlations to both the MCC and BDT scores, according to Kendall’s τ, Spearman’s ρ and Pearson’s r correlation coefficients, when tested on both the CASP8 and CASP9 datasets. The neural network also produced the largest Area Under the Curve score (AUC) when Receiver Operator Characteristic (ROC) analysis was undertaken for the CASP8 dataset. Furthermore, the FunFOLDQA algorithm incorporating the neural network, is shown to add value to FunFOLD, when both methods are employed in combination. This results in a statistically significant improvement over all of the best server methods, the FunFOLD method (6.43%), and one of the top manual groups (FN293) tested on the CASP8 dataset. The FunFOLDQA method was also found to be competitive with the top server methods when tested on the CASP9 dataset. To the best of our knowledge, FunFOLDQA is the first attempt to develop a method that can be used to assess ligand binding site prediction quality, in the absence of experimental data.

Protein-Ion Binding Process on Finite Macromolecular Concentration. A Poisson-Boltzmann and Monte Carlo Study

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Characterization of suramin binding sites on the human group IIA secreted phospholipase A(2) by site-directed mutagenesis and molecular dynamics simulation

Suramin is a polysulphonated naphthylurea with inhibitory activity against the human secreted group IIA phospholipase A(2) (hsPLA2GIIA), and we have investigated suramin binding to recombinant hsPLA2GIIA using site-directed mutagenesis and molecular dynamics (MD) simulations. The changes in suramin binding affinity of 13 cationic residue mutants of the hsPLA2GIIA was strongly correlated with alterations in the inhibition of membrane damaging activity of the protein. Suramin binding to hsPLA2GIIA was also studied by MD simulations, which demonstrated that altered intermolecular potential energy of the suramin/mutant complexes was a reliable indicator of affinity change. Although residues in the C-terminal region play a major role in the stabilization of the hsPLA2GIIA/suramin complex, attractive and repulsive hydrophobic and electrostatic interactions with residues throughout the protein together with the adoption of a bent suramin conformation, all contribute to the stability of the complex. Analysis of the h5PLA2GIIA/suramin interactions allows the prediction of the properties of suramin analogues with improved binding and higher affinities which may be candidates for novel phospholipase A(2) inhibitors. (C) 2012 Elsevier Inc. All rights reserved.

Human leukocyte antigens (HLA) associated drug hypersensitivity: consequences of drug binding to HLA

Recent publications have shown that certain human leukocyte antigen (HLA) alleles are strongly associated with hypersensitivity to particular drugs. As HLA molecules are a critical element in T-cell stimulation, it is no surprise that particular HLA alleles have a direct functional role in the pathogenesis of drug hypersensitivity. In this context, a direct interaction of the relevant drug with HLA molecules as described by the p-i concept appears to be more relevant than presentation of hapten-modified peptides. In some HLA-associated drug hypersensitivity reactions, the presence of a risk allele is a necessary but incomplete factor for disease development. In carbamazepine and HLA-B*15:02, certain T-cell receptor (TCR) repertoires are required for immune activation. This additional requirement may be one of the 'missing links' in explaining why most individuals carrying this allele can tolerate the drug. In contrast, abacavir generates polyclonal T-cell response by interacting specifically with HLA-B*57:01 molecules. T cell stimulation may be due to presentation of abacavir or of altered peptides. While the presence of HLA-B*58:01 allele substantially increases the risk of allopurinol hypersensitivity, it is not an absolute requirement, suggesting that other factors also play an important role. In summary, drug hypersensitivity is the end result of a drug interaction with certain HLA molecules and TCRs, the sum of which determines whether the ensuing immune response is going to be harmful or not.

U7 snRNP-specific Lsm11 protein: dual binding contacts with the 100 kDa zinc finger processing factor (ZFP100) and a ZFP100-independent function in histone RNA 3' end processing

The 3' cleavage generating non-polyadenylated animal histone mRNAs depends on the base pairing between U7 snRNA and a conserved histone pre-mRNA downstream element. This interaction is enhanced by a 100 kDa zinc finger protein (ZFP100) that forms a bridge between an RNA hairpin element upstream of the processing site and the U7 small nuclear ribonucleoprotein (snRNP). The N-terminus of Lsm11, a U7-specific Sm-like protein, was shown to be crucial for histone RNA processing and to bind ZFP100. By further analysing these two functions of Lsm11, we find that Lsm11 and ZFP100 can undergo two interactions, i.e. between the Lsm11 N-terminus and the zinc finger repeats of ZFP100, and between the N-terminus of ZFP100 and the Sm domain of Lsm11, respectively. Both interactions are not specific for the two proteins in vitro, but the second interaction is sufficient for a specific recognition of the U7 snRNP by ZFP100 in cell extracts. Furthermore, clustered point mutations in three phylogenetically conserved regions of the Lsm11 N-terminus impair or abolish histone RNA processing. As these mutations have no effect on the two interactions with ZFP100, these protein regions must play other roles in histone RNA processing, e.g. by contacting the pre-mRNA or additional processing factors.

Development of a Novel Green Fluorescent Protein-Based Binding Assay to Study the Association of Plakins with Intermediate Filament Proteins.

Protein-protein interactions are fundamental for most biological processes, such as the formation of cellular structures and enzymatic complexes or in signaling pathways. The identification and characterization of protein-protein interactions are therefore essential for understanding the mechanisms and regulation of biological systems. The organization and dynamics of the cytoskeleton, as well as its anchorage to specific sites in the plasma membrane and organelles, are regulated by the plakins. These structurally related proteins anchor different cytoskeletal networks to each other and/or to other cellular structures. The association of several plakins with intermediate filaments (IFs) is critical for maintenance of the cytoarchitecture. Pathogenic mutations in the genes encoding different plakins can lead to dramatic manifestations, occurring principally in the skin, striated muscle, and/or nervous system, due to cytoskeletal disorganization resulting in abnormal cell fragility. Nevertheless, it is still unclear how plakins bind to IFs, although some general rules are slowly emerging. We here describe in detail a recently developed protein-protein fluorescence binding assay, based on the production of recombinant proteins tagged with green fluorescent protein (GFP) and their use as fluid-phase fluorescent ligands on immobilized IF proteins. Using this method, we have been able to assess the ability of C-terminal regions of GFP-tagged plakin proteins to bind to distinct IF proteins and IF domains. This simple and sensitive technique, which is expected to facilitate further studies in this area, can also be potentially employed for any kind of protein-protein interaction studies.

The SH3p4/Sh3p8/SH3p13 protein family: Binding partners for synaptojanin and dynamin via a Grb2-like Src homology 3 domain

The GTPase dynamin I and the inositol 5-phosphatase synaptojanin are nerve terminal proteins implicated in synaptic vesicle recycling. Both proteins contain COOH-terminal proline-rich domains that can interact with a variety of Src homology 3 (SH3) domains. A major physiological binding partner for dynamin I and synaptojanin in the nervous system is amphiphysin I, an SH3 domain-containing protein also concentrated in nerve terminals. We have used the proline-rich tail of synaptojanin to screen a rat brain library by the two-hybrid method to identify additional interacting partners of synaptojanin. Three related proteins containing SH3 domains that are closely related to the SH3 domains of Grb2 were isolated: SH3p4, SH3p8, and SH3p13. Further biochemical studies demonstrated that the SH3p4/8/13 proteins bind to both synaptojanin and dynamin I. The SH3p4/8/13 transcripts are differentially expressed in tissues: SH3p4 mRNA was detected only in brain, SH3p13 mRNA was present in brain and testis, and the SH3p8 transcript was detected at similar levels in multiple tissues. Members of the SH3p4/8/13 protein family were found to be concentrated in nerve terminals, and pools of synaptojanin and dynamin I were coprecipitated from brain extracts with antibodies recognizing SH3p4/8/13. These findings underscore the important role of SH3-mediated interactions in synaptic vesicle recycling.

Evidence for protein X binding to a discontinuous epitope on the cellular prion protein during scrapie prion propagation

Studies on the transmission of human (Hu) prions to transgenic (Tg) mice suggested that another molecule provisionally designated protein X participates in the formation of nascent scrapie isoform of prion protein (PrPSc). We report the identification of the site at which protein X binds to the cellular isoform of PrP (PrPC) using scrapie-infected mouse (Mo) neuroblastoma cells transfected with chimeric Hu/MoPrP genes even though protein X has not yet been isolated. Substitution of a Hu residue at position 214 or 218 prevented PrPSc formation. The side chains of these residues protrude from the same surface of the C-terminal α-helix and form a discontinuous epitope with residues 167 and 171 in an adjacent loop. Substitution of a basic residue at positions 167, 171, or 218 also prevented PrPSc formation: at a mechanistic level, these mutant PrPs appear to act as “dominant negatives” by binding protein X and rendering it unavailable for prion propagation. Our findings seem to explain the protective effects of basic polymorphic residues in PrP of humans and sheep and suggest therapeutic and prophylactic approaches to prion diseases.

Sodium channel selectivity filter regulates antiarrhythmic drug binding

Local anesthetic antiarrhythmic drugs block Na+ channels and have important clinical uses. However, the molecular mechanism by which these drugs block the channel has not been established. The family of drugs is characterized by having an ionizable amino group and a hydrophobic tail. We hypothesized that the charged amino group of the drug may interact with charged residues in the channel’s selectivity filter. Mutation of the putative domain III selectivity filter residue of the adult rat skeletal muscle Na+ channel (μ1) K1237E increased resting lidocaine block, but no change was observed in block by neutral analogs of lidocaine. An intermediate effect on the lidocaine block resulted from K1237S and there was no effect from K1237R, implying an electrostatic effect of Lys. Mutation of the other selectivity residues, D400A (domain I), E755A (domain II), and A1529D (domain IV) allowed block by externally applied quaternary membrane-impermeant derivatives of lidocaine (QX314 and QX222) and accelerated recovery from block by internal QX314. Neo-saxitoxin and tetrodotoxin, which occlude the channel pore, reduced the amount of QX314 bound in D400A and A1529D, respectively. Block by outside QX314 in E755A was inhibited by mutation of residues in transmembrane segment S6 of domain IV that are thought to be part of an internal binding site. The results demonstrate that the Na+ channel selectivity filter is involved in interactions with the hydrophilic part of the drugs, and it normally limits extracellular access to and escape from their binding site just within the selectivity filter. Participation of the selectivity ring in antiarrhythmic drug binding and access locates this structure adjacent to the S6 segment.

Polarity of human replication protein A binding to DNA

Replication protein A (RPA), the nuclear single-stranded DNA binding protein is involved in DNA replication, nucleotide excision repair (NER) and homologous recombination. It is a stable heterotrimer consisting of subunits with molecular masses of 70, 32 and 14 kDa (p70, p32 and p14, respectively). Gapped DNA structures are common intermediates during DNA replication and NER. To analyze the interaction of RPA and its subunits with gapped DNA we designed structures containing 9 and 30 nucleotide gaps with a photoreactive arylazido group at the 3′-end of the upstream oligonucleotide or at the 5′-end of the downstream oligonucleotide. UV crosslinking and subsequent analysis showed that the p70 subunit mainly interacts with the 5′-end of DNA irrespective of DNA structure, while the subunit orientation towards the 3′-end of DNA in the gap structures strongly depends on the gap size. The results are compared with the data obtained previously with the primer–template systems containing 5′- or 3′-protruding DNA strands. Our results suggest a model of polar RPA binding to the gapped DNA.