Evolutionary trace analysis of eukaryotic DNA topoisomerase I superfamily: Identification of novel antitumor drug binding site

The studies of novel inhibitors of DNA topoisomerase I (Topo I) have already become very promising in cancer chemotherapy. Identifying the new drug-binding residues is playing an important role in the design and optimization of Topo I inhibitors. The designed compounds may have novel scaffolds, thus will be helpful to overcome the toxicities of current camptothecin (CPT) drugs and may provide a solution to cross resistance with these drugs. Multiple sequence alignments were performed on eukaryotic DNA topoisomerase I superfamily and thus the evolutionary tree was constructed. The Evolutionary Trace method was applied to identify functionally important residues of human Topo I. It has been demonstrated that class-specific hydrophobic residues Ala351, Met428, Pro431 are located around the 7,9-position of CPT, indicating suitable substitution of hydrophobic group on CPT will increase antitumor activity. The conservative residue Lys436 in the superfamily is of particular interest and new CPT derivatives designed based on this residue may greatly increase water solubility of such drugs. It has also been demonstrated that the residues Asn352 and Arg364 were conservative in the superfamily, whose mutation will render CPT resistance. As our molecular docking studies demonstrated they did not make any direct interaction with CPT, they are important drug-binding site residues for future design of novel non-camptothecin lead compounds. This work provided a strong basis for the design and synthesis of novel highly potent CPT derivatives and virtual screening for novel lead compounds.     

Evolutionary trace analysis of eukaryotic DNA topoisomerase I superfamily: Identification of novel antitumor drug binding site

During protein evolutionary processes, protein fam-ily members undergo extensive random mutations and a long period of natural selections, and thus induce the functional evolution and the emergence of subfamily. The evolutionary variation events were recorded in the sequences of protein family members. Therefore, identification of functionally important residues can be achieved by studying residue conservation in protein sequence families. Generally, the residues conserved across the family of…     

Evolutionary trace analysis of the alpha-D-phosphohexomutase superfamily

The alpha-D-phosphohexomutase superfamily is composed of four related enzymes that catalyze a reversible, intramolecular phosphoryl transfer on their sugar substrates. The enzymes in this superfamily play important and diverse roles in carbohydrate metabolism in organisms from bacteria to humans. Recent structural and mechanistic studies of one member of this superfamily, phosphomannomutase/phosphoglucomutase (PMM/PGM) from Pseudomonas aeruginosa, have provided new insights into enzyme mechanism and substrate recognition. Here we use sequence-sequence and sequence-structure comparisons via evolutionary trace analysis to examine 71 members of the alpha-D-phosphohexomutase superfamily. These analyses show that key residues in the active site, including many of those involved in substrate contacts in the P. aeruginosa PMM/PGM complexes, are conserved throughout the enzyme family. Several important regions show class-specific differences in sequence that appear to be correlated with differences in substrate specificity exhibited by subgroups of the family. In addition, we describe the translocation of a 20-residue segment containing the catalytic phosphoserine of phosphoacetylglucosamine mutase, which uniquely identifies members of this subgroup.     

The ligand binding site of the synaptosomal choline transporter: A provisional model based on inhibition studies

A topographic model of the ligand binding site of the choline transporter was deduced from inhibition studies with the help of CPK molecular models. It is posited that there are two identical or closely similar hydrophilic anionic sites separated from each other by an hinged, essentially planar but conformationally flexible cationic hydrophobic domain. Subsequently to attachment of external choline to either one of the anionic sites, both sites cooperate in enveloping the ligand by a Venus fly-trap mechanism. This leads to rapid configurational changes by which the closed-liganded form of the transporter opens up to the interior to release the bound choline. Intracellular K⁺, a ligand for the choline-binding site, is proposed to be countertransported by a reversal of the above mechanism.Special issue dedicated to Dr. Morris H. Aprison.     

Localization of ligand binding site in proteins identified in silico

Knowledge-based models for protein folding assume that the early-stage structural form of a polypeptide is determined by the backbone conformation, followed by hydrophobic collapse. Side chain–side chain interactions, mostly of hydrophobic character, lead to the formation of the hydrophobic core, which seems to stabilize the structure of the protein in its natural environment. The fuzzy-oil-drop model is employed to represent the idealized hydrophobicity distribution in the protein molecule. Comparing it with the one empirically observed in the protein molecule reveals that they are not in agreement. It is shown in this study that the irregularity of hydrophobic distributions is aim-oriented. The character and strength of these irregularities in the organization of the hydrophobic core point to the specificity of a particular protein’s structure/function. When the location of these irregularities is determined versus the idealized fuzzy-oil-drop, function-related areas in the protein molecule can be identified. The presented model can also be used to identify ways in which protein–protein complexes can possibly be created. Active sites can be predicted for any protein structure according to the presented model with the free prediction server at . The implication based on the model presented in this work suggests the necessity of active presence of ligand during the protein folding process simulation. Figure Fuzzy-oil-drop model applied to identify the ligation site in lysozyme complexed with N-acetylglucosamine (PDB ID:1LMQ) in form of hydrophobicity deficiency (ΔH) profile and three-dimensional distribution of on protein surface     

Sequence and structure continuity of evolutionary importance improves protein functional site discovery and annotation

Protein functional sites control most biological processes and are important targets for drug design and protein engineering. To characterize them, the evolutionary trace (ET) ranks the relative importance of residues according to their evolutionary variations. Generally, top‐ranked residues cluster spatially to define evolutionary hotspots that predict functional sites in structures. Here, various functions that measure the physical continuity of ET ranks among neighboring residues in the structure, or in the sequence, are shown to inform sequence selection and to improve functional site resolution. This is shown first, in 110 proteins, for which the overlap between top‐ranked residues and actual functional sites rose by 8% in significance. Then, on a structural proteomic scale, optimized ET led to better 3D structure‐function motifs (3D templates) and, in turn, to enzyme function prediction by the Evolutionary Trace Annotation (ETA) method with better sensitivity of (40% to 53%) and positive predictive value (93% to 94%). This suggests that the similarity of evolutionary importance among neighboring residues in the sequence and in the structure is a universal feature of protein evolution. In practice, this yields a tool for optimizing sequence selections for comparative analysis and, via ET, for better predictions of functional site and function. This should prove useful for the efficient mutational redesign of protein function and for pharmaceutical targeting.     

Binding site on human immunoglobulin G for the affinity ligand HWRGWV

Affinity ligand HWRGWV has demonstrated the ability to isolate human immunoglobulin G (hIgG) from mammalian cell culture media. The ligand specifically binds hIgG through its Fc portion. This work shows that deglycosylation of hIgG has no influence on its binding to the HWRGWV ligand and the ligand does not compete with Protein A or Protein G in binding hIgG. It is suggested by the mass spectrometry (MS) data and docking simulation that HWRGWV binds to the pFc portion of hIgG and interacts with the amino acids in the loop Ser383–Asn389 (SNGQPEN) located in the C_H3 domain. Subsequent modeling has suggested a possible three‐dimensional minimized solution structure for the interaction of hIgG and the HWRGWV ligand. The results support the fact that a peptide as small as a hexamer can have specific interactions with large proteins such as hIgG. Copyright © 2009 John Wiley & Sons, Ltd.     

Reconciliation of classical and reacted‐site probability approaches to allowance for ligand multivalence in binding studies

The objective of this investigation is to engender greater confidence in the validity of binding equations derived for multivalent ligands on the basis of reacted‐site probability theory. To that end, a demonstration of the theoretical interconnection between expressions derived by the classical stepwise equilibria and reacted‐site probability approaches for univalent ligands is followed by the use of the traditional stepwise procedure to derive binding equations for bivalent and trivalent ligands. As well as demonstrating the unwieldy nature of the classical binding equation for multivalent ligand systems, that exercise has allowed numerical simulation to be used to illustrate the equivalence of binding curves generated by the two approaches. The advantages of employing a redefined binding function for multivalent ligands are also confirmed by subjecting the simulated results to a published analytical procedure that has long been overlooked. Copyright © 2013 John Wiley & Sons, Ltd.     

Identification of amino acid residues responsible for the alpha5 subunit binding selectivity of L-655,708, a benzodiazepine binding site ligand at the GABA(A) receptor

L-655,708 is a ligand for the benzodiazepine site of the gamma-aminobutyric acid type A (GABA(A)) receptor that exhibits a 100-fold higher affinity for alpha5-containing receptors compared with alpha1-containing receptors. Molecular biology approaches have been used to determine which residues in the alpha5 subunit are responsible for this selectivity. Two amino acids have been identified, alpha5Thr208 and alpha5Ile215, each of which individually confer approximately 10-fold binding selectivity for the ligand and which together account for the 100-fold higher affinity of this ligand at alpha5-containing receptors. L-655,708 is a partial inverse agonist at the GABA(A) receptor which exhibited no functional selectivity between alpha1- and alpha5-containing receptors and showed no change in efficacy at receptors containing alpha1 subunits where amino acids at both of the sites had been altered to their alpha5 counterparts (alpha1Ser205-Thr,Val212-Ile). In addition to determining the binding selectivity of L-655,708, these amino acid residues also influence the binding affinities of a number of other benzodiazepine (BZ) site ligands. They are thus important elements of the BZ site of the GABA(A) receptor, and further delineate a region just N-terminal to the first transmembrane domain of the receptor alpha subunit that contributes to this binding site.     

Multiple diverse ligands binding at a single protein site: A matter of pre-existing populations

Here, we comment on the steadily increasing body of data showing that proteins with specificity actually bind ligands of diverse shapes, sizes, and composition. Such a phenomenon is not surprising when one considers that binding is a dynamic process with populations in equilibrium and that the shape of the binding site is strongly influenced by the molecular partner. It derives implicitly from the concept of populations. All proteins, specific and nonspecific, exist in ensembles of substates. If the library of ligands in solution is large enough, favorably matching ligands with altered shapes and sizes can be expected to bind, with a redistribution of the protein populations. Point mutations at spatially distant sites may exert large conformational rearrangements and hinge effects, consistent with mutations away from the binding site leading to population shifts and (cross-)drug resistance. A similar effect is observed in protein superfamilies, in which different sequences with similar topologies display similar large-scale dynamic motions. The hinges are frequently at analogous sites, yet with different substrate specificity. Similar topologies yield similar conformational isomers, although with different distributions of population times, owing to the change in the conditions, that is, the change in the sequences. In turn, different distributions relate to binding of different sizes and shapes. Hence, the binding site shape and size are defined by the ligand. They are not independent entities of fixed proportions and cannot be analyzed independently of the binding partner. Such a proposition derives from viewing proteins as dynamic distributions, presenting to the incoming ligands a range of binding site shapes. It illustrates how presumably specific binding molecules can bind multiple ligands. In terms of drug design, the ability of a single receptor to recognize many dissimilar ligands shows the need to consider more diverse molecules. It provides a rationale for higher affinity inhibitors that are not derived from substrates at their transition states and indicates flexible docking schemes.     

The specificity and kinetic properties of trypsin ethylated at the binding site

Treatment of trypsin with triethyloxonium tetrafluoroborate at pH 8, 25 °C, results in abolition of binding to the enzyme of specific cationic substrates and inhibitors. The binding constant of soybean trypsin inhibitor to ethylated trypsin is 10000-fold smaller than to intact trypsin. However, the intrinsic ability of trypsin to recognize and react with nonspecific neutral substrates and inhibitors is not lost, and in several cases even considerably enhanced. Thus ethylated trypsin (Tret) resembles chymotrypsin in its behavior. Trypsin-like enzymes are also affected in a similar manner.     

Characterization of electrostatic binding sites of extracellular polymers by linear programming analysis of titration data

Electrostatic binding sites of extracellular polymeric substances (EPS) were characterized from titration data using linear programming analysis. Test results for three synthetic solutions of given solutes comprising amino, carboxyl, and phenolic groups indicated that this method was able to identify the electrostatic binding sites. For the six sites with pK(a) between 3 and 10, the estimated pK(a) deviated 0.11 +/- 0.09 from the theoretical values, and the estimated concentrations deviated 3.0% +/- 0.9% from the actual concentrations. Two EPS samples were then extracted from a hydrogen-producing sludge (HPS) and a sulfate-reducing biofilm (SRB). Analysis of charge excess data in titration from pH 3 to 11 indicated that the EPS of HPS comprised of five electrostatic binding sites with pK(a) ranging from 3 to 11. The pK(a) values of these binding sites and the possible corresponding functional groups were pK(a) 4.8 (carboxyl), pK(a) 6.0 (carboxyl/phosphoric), pK(a) 7.0 (phosphoric), pK(a) 9.8 (amine/phenolic), and pK(a) 11.0 (hydroxyl). EPS of the SRB comprised five of similar binding sites (with corresponding pK(a) values of 4.4, 6.0, 7.4, 9.4, and 11.0), plus one extra site at pK(a) 8.2, which was likely corresponding to the sulfhydryl group. The total electrostatic binding site concentration of EPS extracted from HPS were 10.88 mmol/g-EPS, of which the highest concentration was from the site of pK(a) 11.0. The corresponding values for the EPS extracted from SRB were 16.44 mmol/g-EPS and pK(a) 4.4. The total concentrations of electrostatic binding sites found in this study were 20- to 30-fold of those reported for bacterial cell surface, implying that EPS might be more crucial in biosorption of metals than bacterial cell surface in wastewater treatment and in bioremediation.     

Insilco analysis of functionally important residues in folate receptors

Lack of crystal structure data of folate binding proteins has left so many questions unanswered (for example, important residues in active site, binding domain, important amino acid residues involved in interactions between ligand and receptor). With sequence alignment and PROSITE motif identification, we attempted to answer evolutionarily significant residues that are of functional importance for ligand binding and that form catalytic sites. We have analyzed 46 different FRs and FBP sequences of various organisms obtained from Genbank. Multiple sequence alignment identified 44 highly conserved identical amino acid residues with 10 cysteine residues and 12 motifs including ECSPNLGPW (which might help in the structural stability of FR).     

Distal heme pocket regulation of ligand binding and stability in soybean leghemoglobin

Leghemoglobins facilitate diffusion of oxygen through root tissue to a bacterial terminal oxidase in much the same way that myoglobin transports oxygen from blood to muscle cell mitochondria. Leghemoglobin serves an additional role as an oxygen scavenger to prevent inhibition of nitrogen fixation. For this purpose, the oxygen affinity of soybean leghemoglobin is 20-fold greater than myoglobin, resulting from an 8-fold faster association rate constant combined with a 3-fold slower dissociation rate constant. Although the biochemical mechanism used by myoglobin to bind oxygen has been described in elegant detail, an explanation for the difference in affinity between these two structurally similar proteins is not obvious. The present work demonstrates that, despite their similar structures, leghemoglobin uses methods different from myoglobin to regulate ligand affinity. Oxygen and carbon monoxide binding to a comprehensive set of leghemoglobin distal heme pocket mutant proteins in comparison to their myoglobin counterparts has revealed some of these mechanisms. The "distal histidine" provides a crucial hydrogen bond to stabilize oxygen in myoglobin but has little effect on bound oxygen in leghemoglobin and is retained mainly for reasons of protein stability and prevention of heme loss. Furthermore, soybean leghemoglobin uses an unusual combination of HisE7 and TyrB10 to sustain a weak stabilizing interaction with bound oxygen. Thus, the leghemoglobin distal heme pocket provides a much lower barrier to oxygen association than occurs in myoglobin and oxygen dissociation is regulated from the proximal heme pocket.     

Identification of a novel ligand binding site in phosphoserine phosphatase from the hyperthermophilic archaeon Thermococcus onnurineus

Phosphoserine phosphatase (PSP) catalyzes the final and irreversible step of L‐serine synthesis by hydrolyzing phosphoserine to produce L ‐serine and inorganic phosphate. Developing a therapeutic drug that interferes with serine production is of great interest to regulate the pathogenicity of some bacteria and control D ‐serine levels in neurological diseases. We determined the crystal structure of PSP from the hyperthermophilic archaeon Thermococcus onnurineus at 1.8 Å resolution, revealing an NDSB ligand bound to a novel site that is located in a fissure between the catalytic domain and the CAP module. The structure shows a half‐open conformation of the CAP 1 module with a unique protruding loop of residues 150–155 that possesses a helical conformation in other structures of homologous PSPs. Activity assays indicate that the enzyme exhibits marginal PSP activity at low temperature but a sharp increase in the k_cat/K_M value, approximately 22 fold, when the temperature is increased. Structural and biochemical analyses suggest that the protruding loop in the active site might be an essential component for the regulation of the activity of PSP from hyperthermophilic T. onnurineus. Identification of this novel binding site distantly located from the catalytic site may be exploited for the development of effective therapeutic allosteric inhibitors against PSP activity. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Location and nature of the residues important for ligand recognition in G-protein coupled receptors

The overall structure of the biogenic amine subclass of the G-protein-coupled receptors, and of their ligand binding sites, is discussed with the aim of highlighting the major structural features of these receptors that are responsible for ligand recognition. A comparison is made between biogenic amine receptors, peptide receptors of the rhodopsin class, and the secretin receptors which all have peptide ligands. The question of where the peptide ligands bind, whether at extracellular sites or within the transmembrane helix bundle, is discussed. The suitability of the rhodopsin crystal structure as a template for construction of homology models is discussed and it is concluded that there are many reasons why a caution should be issued against using it uncritically.