One protein�many functions

The concept of one enzyme-one activity had influenced biochemistry for over half a century. Over 1000 enzymes are now described. Many of them are highly 'specific'. Some of them are crystallized and their three-dimensional structures determined. They range from 12 to 1000 kDa in molecular weight and possess 124 to several hundreds of amino acids. They occur as single polypeptides or multiple-subunit proteins. The active sites are assembled on these by appropriate tertiary folding of the polypeptide chain, or by interaction of the constituent subunits. The substrate is held by the side-chains of a few amino acids at the active site on the surface, occupying a tiny fraction of the total area. What is the bulk of the protein behind the active site doing? Do all proteins have only one function each? Why not a protein have more than one active site on its large surface? Will we discover more than one activity for some proteins? These newer possibilities are emerging and are finding experimental support. Some proteins purified to homogeneity using assay methods for different activities are now recognized to have the same molecular weight and a high degree of homology of amino acid sequence. Obviously they are identical. They represent the phenomenon of one protein-many functions.

Using Correlated Parameters for Improved Ranking of Protein-Protein Docking Decoys

A successful protein-protein docking study culminates in identification of decoys at top ranks with near-native quaternary structures. However, this task remains enigmatic because no generalized scoring functions exist that effectively infer decoys according to the similarity to near-native quaternary structures. Difficulties arise because of the highly irregular nature of the protein surface and the significant variation of the nonbonding and solvation energies based on the chemical composition of the protein-protein interface. In this work, we describe a novel method combining an interface-size filter, a regression model for geometric compatibility (based on two correlated surface and packing parameters), and normalized interaction energy (calculated from correlated nonbonded and solvation energies), to effectively rank decoys from a set of 10,000 decoys. Tests on 30 unbound binary protein-protein complexes show that in 16 cases we can identify at least one decoy in top three ranks having <= 10 angstrom backbone root mean square deviation from true binding geometry. Comparisons with other state-of-art methods confirm the improved ranking power of our method without the use of any experiment-guided restraints, evolutionary information, statistical propensities, or modified interaction energy equations. Tests on 118 less-difficult bound binary protein-protein complexes with <= 35% sequence redundancy at the interface showed that in 77% cases, at least 1 in 10,000 decoys were identified with <= 5 angstrom backbone root mean square deviation from true geometry at first rank. The work will promote the use of new concepts where correlations among parameters provide more robust scoring models. It will facilitate studies involving molecular interactions, including modeling of large macromolecular assemblies and protein structure prediction. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 32: 787-796, 2011.

Ricin-membrane interaction: membrane penetration depth by fluorescence quenching and resonance energy transfer

The entry of the plant toxin ricin and its A- and B-subunits in model membranes in the presence as well as absence of monosialoganglioside (GM(1)) has been studied. Dioleoylphosphatidylcholine and 5-, 10-, and 12-doxyl- or 9,10-dibromophosphatidylcholines serve as quenchers of intrinsic tryptophan fluorescence of the proteins. The parallax method of Chattopadhyay and London [(1987) Biochemistry 26, 39-45] has been employed to measure the average membrane penetration depth of tryptophans of ricin and its B-chain and the actual depth of the sole Trp 211 in the A-chain. The results indicate that both of the chains as well as intact ricin penetrate the membrane deeply and the C-terminal end of the A-chain is well inside the bilayer, especially at pH 4.5. An extrinsic probe N-(iodoacetyl)-N'-(5-sulfo-1-naphthyl) ethylenediamine (I-AEDANS) has been attached to Cys 259 of the A-chain, and the kinetics of penetration has been followed by monitoring the increase in AEDANS fluorescence at 480 nm. The insertion follows first-order kinetics, and the rate constant is higher at a lower pH. The energy transfer distance analysis between Trp 211 and AEDANS points out that the conformation of the A-chain changes as it inserts into the membrane. CD studies indicate that the helicity of the proteins increases after penetration, which implies that some of the unordered structure in the native protein is converted to the ordered form during this process. Hydrophobic forces seem to be responsible for stabilizing a particular protein conformation inside the membrane.

Testis-specific histone (H1t) is not phosphorylated and has a weak interaction with chromatin

Soluble chromatin was prepared from rat testes after a brief micrococcal nuclease digestion. After adsorption onto hydroxylapatite at low ionic strength, the histone Hl subtypes were eluted with a shallow salt gradient of 0.3 M NaCl to 0.7 M NaCl. Histone Hlt was eluted at 0.4 MNaCl, while histones H1a and Hlc were eluted at 0.43 M NaCl and 0.45 M respectively. The extreme divergence of the amino acid sequence of the C-terminal half of histone Hlt, the major DNA binding domain of histone Hl, from that of the somatic consensus sequence may contribute to the weaker interaction of histone Hlt with the rat testis chromatin. Further, histone Hlt was not phosphorylated in vivo in contrast to histone Hla and Hlc, as is evident from the observation that histone Hlt lacks the SPKK motif recognized by the CDC-2kinase or the RR/KXS motif recognized by protein kinase A.

Thermodynamics of lectin-sugar interaction: Binding of sugars to winged bean (Psophocarpus tetragonolobus) basic agglutinin (WBAI)

Combining site of WBAI is extended and encompasses all the residues of blood group A-reactive trisaccharide [GalNAcalpha3Galbeta4Glc]. Though both of the fucose residues of A-pentasaccharide [GalNAcalpha(Fucalpha2)3Galbeta(Fucalpha3)4Glc] do not directly interact, with the combining site they thermodynamically favour the interaction of GalNAcalpha3Galbeta4Glc part of the molecule by imposing a sterically favourable orientation of the binding epitope viz. GalNAcalpha3Galbeta4Glc of the saccharide. Binding of sugars is driven by enthalpy and is devoid of heat capacity changes. This together with enthalpy-entropy compensation observed for these processes underscore the importance of water reorganization as being one of the principal determinant of protein-sugar interactions.

Characterization of the DNA-binding domain of beta protein, a component of phage lambda Red-pathway by UV catalyzed cross-linking

beta protein, a key component of Red-pathway of phage lambda is necessary for its growth and general genetic recombination in recombination-deficient mutants of Escherichia coli. To facilitate studies on structure-function relationships, we overexpressed beta protein and purified it to homogeneity. A chemical cross-linking reagent, glutaraldehyde, was used to stabilize the physical association of beta protein in solution. A 67-kDa band, corresponding to homodimer, was identified after separation by SDS-polyacrylamide gel electrophoresis. Stoichiometric measurements indicated a site-size of 1 monomer of beta protein/5 nucleotide residues. Electrophoretic gel mobility shift assays suggested that beta protein formed stable nucleoprotein complexes with 36-mer, but not with 27- or 17-mer DNA. Interestingly, the interaction of beta protein with DNA and the stability of nucleoprotein complexes was dependent on the presence of MgCl2, and the binding was abolished by 250 mM NaCl. The K-d of beta protein binding to 36-mer DNA was on the order of 1.8 x 10(-6) M. Photochemical cross-linking of native beta protein or its fragments, generated by chymotrypsin, to 36-mer DNA was performed to identify its DNA-binding domain. Characterization of the cross-linked peptide disclosed that amino acids required for DNA-binding specificity resided within a 20-kDa peptide at the N-terminal end. These findings provide a basis for further understanding oi the structure and function of beta protein.

Interaction of Bacillus polymyxa with some oxide minerals with reference to mineral beneficiation and environmental control

Interaction of Bacillus polymyxa with calcite, hematite, corundum and quartz resulted in significant surface chemical changes not only of the cells but also in the minerals. Both the cell surfaces as well as quartz particles were rendered more hydrophobic after mutual interaction, whilst the rest of the minerals exhibited enhanced hydrophilicity after interaction with the bacteria. The bacteria were also observed to be capable of dissolving calcite, hematite and corundum and biosorbing the dissolved metal ions to varying extents. An excess of polysaccharides could be observed on biotreated calcite, hematite and corundum while the predominance of a protein-based metabolic product was evident on quartz surfaces. The utility of bioprocessing in the beneficiation of the above minerals through bioflotation and bioflocculation is demonstrated. (C) 1997 Elsevier Science Ltd.

Interaction of A beta peptide (1-40) with amino acid aluminium complexes: relevance to Alzheimer's disease

A beta (39-43 aminoacid residues) is the principal peptide component of amyloid deposits in Alzheimer's disease (AD). A beta peptide is derived from the amyloid precursor protein (APP) in which mutations give rise to many forms of familial AD. Aluminium is reported to play a key role in inducing conformational change in the synthetic beta-amyloid peptide (1-40)from alpha-helix to beta-pleated sheet, leading to aggregation and fibrillar formation. We have studied the interaction of amino acid-Al complexes such as D-Asp-Al and L-Glu-Al with A beta(1-40) in TFE/buffer (70% TFE and 30% H2O v/v pH 6.7) mixture using CD spectroscopy. The interaction of either of these amino acid complexes with A beta(1-40) results in loss of alpha-helical content and the peptide is more unstructured compared to free Al3+ in the solution. Our data strongly support the idea, that the Al3+ in the form of aminoacid-Al complexes is more effective in inducing random coil conformation in the A beta peptide than the free Al3+ present in the solution.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

PRUNE and PROBE-two modular web services for protein-protein docking

The protein-protein docking programs typically perform four major tasks: (i) generation of docking poses, (ii) selecting a subset of poses, (iii) their structural refinement and (iv) scoring, ranking for the final assessment of the true quaternary structure. Although the tasks can be integrated or performed in a serial order, they are by nature modular, allowing an opportunity to substitute one algorithm with another. We have implemented two modular web services, (i) PRUNE: to select a subset of docking poses generated during sampling search (http://pallab.serc.iisc.ernet.in/prune) and (ii) PROBE: to refine, score and rank them (http://pallab.serc.iisc.ernet.in/probe). The former uses a new interface area based edge-scoring function to eliminate > 95% of the poses generated during docking search. In contrast to other multi-parameter-based screening functions, this single parameter based elimination reduces the computational time significantly, in addition to increasing the chances of selecting native-like models in the top rank list. The PROBE server performs ranking of pruned poses, after structure refinement and scoring using a regression model for geometric compatibility, and normalized interaction energy. While web-service similar to PROBE is infrequent, no web-service akin to PRUNE has been described before. Both the servers are publicly accessible and free for use.

The repertoire of protein kinases encoded in the draft version of the human genome: atypical variations and uncommon domain combinations

Background: Phosphorylation by protein kinases is central to cellular signal transduction. Abnormal functioning of kinases has been implicated in developmental disorders and malignancies. Their activity is regulated by second messengers and by the binding of associated domains, which are also influential in translocating the catalytic component to their substrate sites, in mediating interaction with other proteins and carrying out their biological roles. Results: Using sensitive profile-search methods and manual analysis, the human genome has been surveyed for protein kinases. A set of 448 sequences, which show significant similarity to protein kinases and contain the critical residues essential for kinase function, have been selected for an analysis of domain combinations after classifying the kinase domains into subfamilies. The unusual domain combinations in particular kinases suggest their involvement in ubiquitination pathways and alternative modes of regulation for mitogen-activated protein kinase kinases (MAPKKs) and cyclin-dependent kinase (CDK)-like kinases. Previously unexplored kinases have been implicated in osteoblast differentiation and embryonic development on the basis of homology with kinases of known functions from other organisms. Kinases potentially unique to vertebrates are involved in highly evolved processes such as apoptosis, protein translation and tyrosine kinase signaling. In addition to coevolution with the kinase domain, duplication and recruitment of non-catalytic domains is apparent in signaling domains such as the PH, DAG-PE, SH2 and SH3 domains. Conclusions: Expansion of the functional repertoire and possible existence of alternative modes of regulation of certain kinases is suggested by their uncommon domain combinations. Experimental verification of the predicted implications of these kinases could enhance our understanding of their biological roles.

The restriction mapping problem revisited

In computational molecular biology, the aim of restriction mapping is to locate the restriction sites of a given enzyme on a DNA molecule. Double digest and partial digest are two well-studied techniques for restriction mapping. While double digest is NP-complete, there is no known polynomial-time algorithm for partial digest. Another disadvantage of the above techniques is that there can be multiple solutions for reconstruction. In this paper, we study a simple technique called labeled partial digest for restriction mapping. We give a fast polynomial time (O(n(2) log n) worst-case) algorithm for finding all the n sites of a DNA molecule using this technique. An important advantage of the algorithm is the unique reconstruction of the DNA molecule from the digest. The technique is also robust in handling errors in fragment lengths which arises in the laboratory. We give a robust O(n(4)) worst-case algorithm that can provably tolerate an absolute error of O(Delta/n) (where Delta is the minimum inter-site distance), while giving a unique reconstruction. We test our theoretical results by simulating the performance of the algorithm on a real DNA molecule. Motivated by the similarity to the labeled partial digest problem, we address a related problem of interest-the de novo peptide sequencing problem (ACM-SIAM Symposium on Discrete Algorithms (SODA), 2000, pp. 389-398), which arises in the reconstruction of the peptide sequence of a protein molecule. We give a simple and efficient algorithm for the problem without using dynamic programming. The algorithm runs in time O(k log k), where k is the number of ions and is an improvement over the algorithm in Chen et al. (C) 2002 Elsevier Science (USA). All rights reserved.

Amino acid interaction preferences in helical membrane proteins

Membrane proteins are involved in a number of important biological functions. Yet, they are poorly understood from the structure and folding point of view. The external environment being drastically different from that of globular proteins, the intra-protein interactions in membrane proteins are also expected to be different. Hence, statistical potentials representing the features of inter-residue interactions based exclusively on the structures of membrane proteins are much needed. Currently, a reasonable number of structures are available, making it possible to undertake such an analysis on membrane proteins. In this study we have examined the inter-residue interaction propensities of amino acids in the membrane spanning regions of the alpha-helical membrane (HM) proteins. Recently we have shown that valuable information can be obtained on globular proteins by the evaluation of the pair-wise interactions of amino acids by classifying them into different structural environments, based on factors such as the secondary structure or the number of contacts that a residue can make. Here we have explored the possible ways of classifying the intra-protein environment of HM proteins and have developed scoring functions based on different classification schemes. On evaluation of different schemes, we find that the scheme which classifies amino acids to different intra-contact environment is the most promising one. Based on this classification scheme, we also redefine the hydrophobicity scale of amino acids in HM proteins.

PocketAlign A Novel Algorithm for Aligning Binding Sites in Protein Structures

A fundamental task in bioinformatics involves a transfer of knowledge from one protein molecule onto another by way of recognizing similarities. Such similarities are obtained at different levels, that of sequence, whole fold, or important substructures. Comparison of binding sites is important to understand functional similarities among the proteins and also to understand drug cross-reactivities. Current methods in literature have their own merits and demerits, warranting exploration of newer concepts and algorithms, especially for large-scale comparisons and for obtaining accurate residue-wise mappings. Here, we report the development of a new algorithm, PocketAlign, for obtaining structural superpositions of binding sites. The software is available as a web-service at http://proline.physicslisc.emetin/pocketalign/. The algorithm encodes shape descriptors in the form of geometric perspectives, supplemented by chemical group classification. The shape descriptor considers several perspectives with each residue as the focus and captures relative distribution of residues around it in a given site. Residue-wise pairings are computed by comparing the set of perspectives of the first site with that of the second, followed by a greedy approach that incrementally combines residue pairings into a mapping. The mappings in different frames are then evaluated by different metrics encoding the extent of alignment of individual geometric perspectives. Different initial seed alignments are computed, each subsequently extended by detecting consequential atomic alignments in a three-dimensional grid, and the best 500 stored in a database. Alignments are then ranked, and the top scoring alignments reported, which are then streamed into Pymol for visualization and analyses. The method is validated for accuracy and sensitivity and benchmarked against existing methods. An advantage of PocketAlign, as compared to some of the existing tools available for binding site comparison in literature, is that it explores different schemes for identifying an alignment thus has a better potential to capture similarities in ligand recognition abilities. PocketAlign, by finding a detailed alignment of a pair of sites, provides insights as to why two sites are similar and which set of residues and atoms contribute to the similarity.

Interplay between NS3 protease and human La protein regulates translation-replication switch of Hepatitis C virus

HCV NS3 protein plays a central role in viral polyprotein processing and RNA replication. We demonstrate that the NS3 protease (NS3(pro)) domain alone can specifically bind to HCV-IRES RNA, predominantly in the SLIV region. The cleavage activity of the NS3 protease domain is reduced upon HCV-RNA binding. More importantly, NS3(pro) binding to the SLIV hinders the interaction of La protein, a cellular IRES-trans acting factor required for HCV IRES-mediated translation, resulting in inhibition of HCV-IRES activity. Although overexpression of both NS3(pro) as well as the full length NS3 protein decreased the level of HCV IRES mediated translation, replication of HCV replicon RNA was enhanced significantly. These observations suggest that the NS3(pro) binding to HCV IRES reduces translation in favor of RNA replication. The competition between the host factor (La) and the viral protein (NS3) for binding to HCV IRES might regulate the molecular switch from translation to replication of HCV.