A computational model of the 5-HT3 receptor extracellular domain: search for ligand binding sites

A three-dimensional model of the 5-HT₃ receptor extarcellular domain has been derived on the basis of the nicotinic acetylcholine receptor model recently published by Tsigelny et al. Maximum complementarity between the position and characteristics of mutated residues putatively involved in ligand interaction and the pharmacophoric elements derived by the indirect approach applied on several series of 5-HT₃ ligands have been exploited to gain insights into the ligand binding modalities and to speculate on the mechanistic role of the structural components. The analysis of the three-dimensional model allows one to distinguish among amino acids that exert key roles in ligand interactions, subunit architecture, receptor assembly and receptor dynamics. For some of these, alternative roles with respect to the ones hypothesized by experimentalists are assigned. Different binding modalities for agonists and antagonists are highlighted, and residues which probably play a role in the transduction of binding into a change in conformational state of the receptor are suggested. Received: 27 July 2000 / Accepted: 15 September 2000 / Published online: 21 December 2000     

Improved mapping of protein binding sites

Computational mapping methods place molecular probes – small molecules or functional groups – on a protein surface in order to identify the most favorable binding positions by calculating an interaction potential. Mapping is an important step in a number of flexible docking and drug design algorithms. We have developed improved algorithms for mapping protein surfaces using small organic molecules as molecular probes. The calculations reproduce the binding of eight organic solvents to lysozyme as observed by NMR, as well as the binding of four solvents to thermolysin, in good agreement with x-ray data. Application to protein tyrosine phosphatase 1B shows that the information provided by the mapping can be very useful for drug design. We also studied why the organic solvents bind in the active site of proteins, in spite of the availability of alternative pockets that can very tightly accommodate some of the probes. A possible explanation is that the binding in the relatively large active site retains a number of rotational states, and hence leads to smaller entropy loss than the binding elsewhere else. Indeed, the mapping reveals that the clusters of the ligand molecules in the protein's active site contain different rotational-translational conformers, which represent different local minima of the free energy surface. In order to study the transitions between different conformers, reaction path and molecular dynamics calculations were performed. Results show that most of the rotational states are separated by low free energy barriers at the experimental temperature, and hence the entropy of binding in the active site is expected to be high.     

Exploiting three kinds of interface propensities to identify protein binding sites

Predicting the binding sites between two interacting proteins provides important clues to the function of a protein. In this study, we present a building block of proteins called order profiles to use the evolutionary information of the protein sequence frequency profiles and apply this building block to produce a class of propensities called order profile interface propensities. For comparisons, we revisit the usage of residue interface propensities and binary profile interface propensities for protein binding site prediction. Each kind of propensities combined with sequence profiles and accessible surface areas are inputted into SVM. When tested on four types of complexes (hetero-permanent complexes, hetero-transient complexes, homo-permanent complexes and homo-transient complexes), experimental results show that the order profile interface propensities are better than residue interface propensities and binary profile interface propensities. Therefore, order profile is a suitable profile-level building block of the protein sequences and can be widely used in many tasks of computational biology, such as the sequence alignment, the prediction of domain boundary, the designation of knowledge-based potentials and the protein remote homology detection.     

Comprehensive study of proteins that interact with microcystin-LR

We carried out a comprehensive study of proteins that exhibit specific interactions with a naturally occurring toxin, microcystin (MC)-LR, in order to gain insight into the unknown underlying mechanism of MC virulence. This audacious study employed a simple affinity test that used MC-LR immobilized on an original ethylene oxide based monolithic solid phase (Moli-gel), and swine liver lysate. Some of the proteins that interacted with MC-LR on this original affinity resin were separated by SDS-PAGE, measured by nano-LC/MS/MS after trypsin digestion, and identified using a Mascot database search. Protein sequence analyses revealed that glutathione S-transferase (GST) was one of the candidate target proteins for MC-LR. This protein was confirmed as a target protein for MC-LR based on the results of for the inhibition of an enzymatic reaction by Dhb-MC-LR. Moreover, L-3-hydroxyacyl coenzyme A dehydrogenase (HDHA) was shown to be one of the proteins that specifically interacts with MC-LR. Our results demonstrated that our analytical systems based on an original affinity resin and nano-LC/MS/MS were effective for target protein research.     

MCSS-based predictions of RNA binding sites

The diversity of RNA tertiary structures provides the basis for specific recognition by proteins or small molecules. To investigate the structural basis and the energetics which control RNA-ligand interactions, favorable RNA binding sites are identified using the MCSS method, which has been employed previously only for protein receptors. Two different RNAs for which the structures have been determined by NMR spectroscopy were examined: two structures of the TAR RNA which contains an arginine binding site, and the structure of the 16S rRNA which contains an aminoglycoside binding site (paromomycin). In accord with the MCSS methodology, the functional groups representing the entire ligand or only part of it (one residue in the case of the aminoglycosides) are first replicated and distributed with random positions and orientations around the target and then energy minimized in the force field of the target RNA. The Coulombic term and the dielectric constant of the force field are adjusted to approximate the effects of solvent-screening and counterions. Optimal force field parameters are determined to reproduce the binding mode of arginine to the TAR RNA. The more favorable binding sites for each residue of the aminoglycoside ligands are then calculated and compared with the binding sites observed experimentally. The predictability of the method is evaluated and refinements are proposed to improve its accuracy. Received: 24 April 1998 / Accepted: 4 August 1998 / Published online: 7 December 1998     

Ligand-binding prediction in the resistance-nodulation-cell division (RND) proteins

The resistance-nodulation-cell division (RND) protein family is a ubiquitous group of proteins primarily present in bacteria. These proteins, involved in the transport of multiple drugs across the cell envelope in bacteria, exhibit broad substrate specificity and act like efflux pumps. In this work, a protein belonging to the RND protein family, AcrB of Escherichia coli was used as a working model to predict in silico the compounds transported by 47 RND proteins. From AcrB we extracted and clustered 14 amino acids directly involved in substrate interactions. This clustering provides enough information to identify 16 groups that correlates with the ligand they extrude, such as proteins expelling aromatic hydrocarbons (SrpB cluster) or proteins expelling heavy metals (CnrA cluster). The relationship between conserved, cluster-specific and variable residues indicates that although the ligand-binding domain is conserved in structure, it has enough flexibility to recognize specifically a diversity of molecules.     

Statistical relationships among docking scores for different protein binding sites

This report describes the existence of statistical relationships among scores computed with the DOCK program for a library of small molecules and a panel of protein binding sites. Multivariate relationships are observed in docking scores computed for a constant set of ligands in different binding sites of proteins that are dissimilar in structure and function. The structural basis for the correlations found among scores is analyzed in terms of size, shape and charge characteristics of the binding sites considered. Interestingly, these results parallel a growing body of evidence demonstrating the promiscuous behavior of small molecules in their interactions with macromolecules that could have an impact in future efforts in drug design.     

Enhancement of the StreptoTag method for isolation of endogenously expressed proteins with complex RNA binding targets

StreptoTag is a novel affinity chromatography-based method for the isolation of high- and low-affinity RNA binding proteins. Originally it was shown possible to isolate recombinant protein from yeast or bacterial extracts using small, specific, well-characterised RNA binding targets. Here we show that using an enhanced aptamer it is not only possible to efficiently immobilise large, highly structured RNA binding targets onto the streptomycin columns but also that the StreptoTag method can be used for the isolation and purification of endogenously expressed regulatory proteins, with relatively low abundance, from eukaryotic extracts. As an example for this we uncover the identity of a karyophilic cellular protein which specifically binds to an area within the large, highly folded structure that characterises the mRNA from the unique 3' region (U3) of the mouse mammary tumour virus (MMTV) long terminal repeat (LTR). Hence, this method is now suitable for the quick and efficient isolation and identification of novel RNA binding proteins such as regulatory factors.     

Isothermal titration microcalorimetric method for studying the combined ligand binding with application to the binding of ethylurea and (N,N)dimethylurea on urease

A simple, novel method was introduced for determining equilibrium constants and enthalpies of binding of two different competitive ligands on a macromolecule by isothermal titration microcalorimetry technique. This method was applied to the simultaneous binding of ethylurea (I) and (N,N)dimethylurea (X), on jack-bean urease at pH 7.0 (tris-base; 30 mM) at 27°C. The dissociation equilibrium constants measured by this method were markedly consistent with inhibition constants obtained from assay of enzyme activities in the presence of I and X.     

Design and synthesis of affinity ligands and relation of their structure with adsorption of proteins

Affinity chromatography has played an increasingly important role both in the pharmaceutical industry and academic research. In the present study, we report our preliminary investigation on the relationship between the affinity ligand structure and its adsorption to multi-protein samples. The structure of the ligands, including the size of the ring (cyclic group) and the length of the chain (linear group), has a great impact on the adsorption of ligands to proteins. Meanwhile, the functional groups that the ligands carry are also closely related to the adsorption of ligands to proteins. This research provides good guidance for the design and synthesis of affinity materials in affinity chromatography. It is also useful to other protein-ligand interaction-related research.     

Identification of nuclear proteins that interact with platinum-modified DNA by photoaffinity labeling

Interactions between cellular proteins and cisplatin-modified DNA are important in determining the anticancer activity of the drug. To develop a general approach for identifying proteins that mediate cellular responses to cisplatin, photoreactive cisplatin analogues having a tethered benzophenone moiety were prepared and used to form the major 1,2-intrastrand platinum-DNA cross-links. Upon irradiation of the platinated DNA dissolved in a HeLa nuclear extract, the appended photolabile benzophenone group generates a highly reactive species that binds irreversibly to cellular proteins that interact with the probe. Several DNA-protein cross-linked adducts were identified that may function in the cellular processing of cisplatin-DNA adducts. Of these, PARP-1 had not previously been demonstrated directly to contact Pt-DNA cross-links in human cells.     

Ion exchange chromatography of monoclonal antibodies: Effect of resin ligand density on dynamic binding capacity

Dynamic binding capacity (DBC) of a monoclonal antibody on agarose based strong cation exchange resins is determined as a function of resin ligand density, apparent pore size of the base matrix, and protein charge. The maximum DBC is found to be unaffected by resin ligand density, apparent pore size, or protein charge within the tested range. The critical conductivity (conductivity at maximum DBC) is seen to vary with ligand density. It is hypothesized that the maximum DBC is determined by the effective size of the proteins and the proximity to which they can approach one another. Once a certain minimum resin ligand density is supplied, additional ligand is not beneficial in terms of resin capacity. Additional ligand can provide flexibility in designing ion exchange resins for a particular application as the critical conductivity could be matched to the feedstock conductivity and it may also affect the selectivity.     

Experimental design and measurement uncertainty in ligand binding studies by affinity capillary electrophoresis

In all life sciences ligand binding assays (LBAs) play a crucial role. Unfortunately these are very error prone. One part of this uncertainty results from the unavoidable random measurement uncertainty, another part can be attributed to the experimental design. To investigate the latter, uncertainty propagation was evaluated as a function of the given experimental design. A design space including the normalized maximum response range (nMRR), the data point position (DPP), the data point range (DPR) and the number of data points (NoDP) was defined. Based on ten measured ms ACE source data sets 20 specific parameter sets were selected by Design of Experiments. Monte Carlo simulations using 100 000 repeats for every parameter set were employed. The resulting measurement uncertainty propagation factors (measurement uncertainty multiplier: MUM) were used to describe the whole design space by polynomial regression. The resulting 5‐dimensional response surface was investigated to evaluate the design parameter's influence and to find the minimal uncertainty propagation. It could be shown, that the nMRR is of highest importance, followed by DPP and DPR. Interestingly, the NoDP is less relevant. However, the interactions of the four parameters need to be carefully considered during design optimization. Using at least five data points which cover over 40% of the upper part of the binding hyperbola (DPP > 0.57) the MUM will be minimized (MUM approximately 1.5) when the nMRR is appropriate. It is possible to reduce the measurement uncertainty propagation more than one order of magnitude.     

A microcalorimetry and binding study on interaction of dodecyl trimethylammonium bromide with wigeon hemoglobin

The thermodynamic parameters for the binding of dodecyl trimethylammonium bromide (DTAB) with wigeon hemoglobin (Hb) in aqueous solution at various pH and 27 °C have been measured by equilibrium dialysis and titration microcalorimetry techniques. The Scatchard plots represent unusual features at neutral and alkaline pH and specific binding at acidic pH. This leads us to analyze the binding data by fitting the data to the Hill equation for multiclasses of binding sites. The best fit was obtained with the equation for one class at acidic pH and two classes at neutral and alkaline pH. The thermodynamic analysis of the binding process shows that the strength of binding at neutral pH is more than these at other pH values. This can be related to the more accessible hydrophobic surface area of wigeon hemoglobin at this pH. The endothermic enthalpy data which was measured by microcalorimetry confirms the binding data analysis and represents the more regular and stable structure of wigeon hemoglobin at neutral pH.