Heat and SDS insensitive NDK dimers are largely stabilised by hydrophobic interaction to form functional hexamer in Mycobacterium smegmatis

The primary structure and function of nucleoside diphosphate kinase (NDK), a substrate non-specific enzyme involved in the maintenance of nucleotide pools is also implicated to play pivotal roles in many other cellular processes. NDK is conserved from bacteria to human and forms a homotetramer or hexamer to exhibit its biological activity. However, the nature of the functional oligomeric form of the enzyme differs among different organisms. The functional form of NDKs from many bacterial systems, including that of the human pathogen, Mycobacterium tuberculosis (MtuNDK), is a hexamer, although some bacterial NDKs are tetrameric in nature. The present study addresses the oligomeric property of MsmNDK and how a dimer, the basic subunit of a functional hexamer, is stabilized by hydrogen bonds and hydrophobic interactions. Homology modeling was generated using the three-dimensional structure of MtuNDK as a template; the residues interacting at the monomer-monomer interface of MsmNDK were mapped. Using recombinant enzymes of wild type, catalytically inactive mutant, and monomer-monomer interactive mutants of MsmNDK, the stability of the dimer was verified under heat, SDS, low pH, and methanol. The predicted residues (Gln17, Ser24 and Glu27) were engaged in dimer formation, however the mutated proteins retained the ATPase and GTPase activity even after introducing single (MsmNDK- Q17A, MsmNDK-E27A, and MsmNDK-E27Q) and double (MsmNDK-E27A/Q17A) mutation. However, the monomer monomer interaction could be abolished using methanol, indicating the stabilization of the monomer-monomer interaction by hydrophobic interaction.

Probing p300/CBP associated Factor (PCAF)-dependent pathways with a small molecule inhibitor

PCAF (KAT2B) belongs to the GNAT family of lysine acetyltransferases (KAT) and specifically acetylates the histone H3K9 residue and several nonhistone proteins. PCAF is also a transcriptional coactivator. Due to the lack of a PCAF KAT-specific small molecule inhibitor, the exclusive role of the acetyltransferase activity of PCAF is not well understood. Here, we report that a natural compound of the hydroxybenzoquinone class, embelin, specifically inhibits H3Lys9 acetylation in mice and inhibits recombinant PCAF-mediated acetylation with near complete specificity in vitro. Furthermore, using embelin, we have identified the gene networks that are regulated by PCAF during muscle differentiation, further highlighting the broader regulatory functions of PCAF in muscle differentiation in addition to the regulation via MyoD acetylation.

Depth: a web server to compute depth, cavity sizes, detect potential small-molecule ligand-binding cavities and predict the pK(a) of ionizable residues in proteins

Residue depth accurately measures burial and parameterizes local protein environment. Depth is the distance of any atom/residue to the closest bulk water. We consider the non-bulk waters to occupy cavities, whose volumes are determined using a Voronoi procedure. Our estimation of cavity sizes is statistically superior to estimates made by CASTp and VOIDOO, and on par with McVol over a data set of 40 cavities. Our calculated cavity volumes correlated best with the experimentally determined destabilization of 34 mutants from five proteins. Some of the cavities identified are capable of binding small molecule ligands. In this study, we have enhanced our depth-based predictions of binding sites by including evolutionary information. We have demonstrated that on a database (LigASite) of similar to 200 proteins, we perform on par with ConCavity and better than MetaPocket 2.0. Our predictions, while less sensitive, are more specific and precise. Finally, we use depth (and other features) to predict pK(a)s of GLU, ASP, LYS and HIS residues. Our results produce an average error of just <1 pH unit over 60 predictions. Our simple empirical method is statistically on par with two and superior to three other methods while inferior to only one. The DEPTH server (http://mspc.bii.a-star.edu.sg/depth/) is an ideal tool for rapid yet accurate structural analyses of protein structures.

Interaction of Silver Nanoparticles with Serum Proteins Affects Their Antimicrobial Activity In Vivo

The emergence of multidrug-resistant bacteria is a global threat for human society. There exist recorded data that silver was used as an antimicrobial agent by the ancient Greeks and Romans during the 8th century. Silver nanoparticles (AgNPs) are of potential interest because of their effective antibacterial and antiviral activities, with minimal cytotoxic effects on the cells. However, very few reports have shown the usage of AgNPs for antibacterial therapy in vivo. In this study, we deciphered the importance of the chosen methods for synthesis and capping of AgNPs for their improved activity in vivo. The interaction of AgNPs with serum albumin has a significant effect on their antibacterial activity. It was observed that uncapped AgNPs exhibited no antibacterial activity in the presence of serum proteins, due to the interaction with bovine serum albumin (BSA), which was confirmed by UV-Vis spectroscopy. However, capped AgNPs with citrate or poly(vinylpyrrolidone)] exhibited antibacterial properties due to minimized interactions with serum proteins. The damage in the bacterial membrane was assessed by flow cytometry, which also showed that only capped AgNPs exhibited antibacterial properties, even in the presence of BSA. In order to understand the in vivo relevance of the antibacterial activities of different AgNPs, a murine salmonellosis model was used. It was conclusively proved that AgNPs capped with citrate or PVP exhibited significant antibacterial activities in vivo against Salmonella infection compared to uncapped AgNPs. These results clearly demonstrate the importance of capping agents and the synthesis method for AgNPs in their use as antimicrobial agents for therapeutic purposes.

Understanding the role of domain-domain linkers in the spatial orientation of domains in multi-domain proteins

Inter-domain linkers (IDLs)' bridge flanking domains and support inter-domain communication in multi-domain proteins. Their sequence and conformational preferences enable them to carry out varied functions. They also provide sufficient flexibility to facilitate domain motions and, in conjunction with the interacting interfaces, they also regulate the inter-domain geometry (IDG). In spite of the basic intuitive understanding of the inter-domain orientations with respect to linker conformations and interfaces, we still do not entirely understand the precise relationship among the three. We show that IDG is evolutionarily well conserved and is constrained by the domain-domain interface interactions. The IDLs modulate the interactions by varying their lengths, conformations and local structure, thereby affecting the overall IDG. Results of our analysis provide guidelines in modelling of multi-domain proteins from the tertiary structures of constituent domain components.

Form and Function of Proteins. Historical Background and the Indian Effort in Macromolecular Crystallography

Global efforts in macromolecular crystallography started in the thirties of the last century. However, definitive results began to emerge only in the late fifties and the early sixties. India has a long tradition in crystallography. The country had a head start in theoretical and computational structural biology, thanks to the efforts of G.N. Ramachandran and his colleagues in the fifties and the sixties. However, macromolecular crystallography got off the ground in India only in the eighties, particularly after the Bangalore group received adequate support from the Department of Science and Technology under their Thrust Area Programme. The Bangalore centre was also identified as a national nucleus for the development of the area in the country. Since then work in the area has spread widely and is being carried out by several groups, mainly led by scientists trained at Bangalore or their descendents, in about thirty institutions in India. In addition to the Department of Science and Technology, the effort is now supported by other agencies like the Department of Biotechnology and the Council of Scientific and Industrial Research. The problems addressed by macromolecular crystallographers in India encompass almost all aspects of modern biology. Indian efforts in macromolecular crystallography have also become an important component of the international efforts in the area.

Viscoelastic Behavior of Human Lamin A Proteins in the Context of Dilated Cardiomyopathy

Lamins are intermediate filament proteins of type V constituting a nuclear lamina or filamentous meshwork which lines the nucleoplasmic side of the inner nuclear membrane. This protein mesh provides a supporting scaffold for the nuclear envelope and tethers interphase chromosome to the nuclear periphery. Mutations of mainly A-type lamins are found to be causative for at least 11 human diseases collectively termed as laminopathies majority of which are characterised by aberrant nuclei with altered structural rigidity, deformability and poor mechanotransduction behaviour. But the investigation of viscoelastic behavior of lamin A continues to elude the field. In order to address this problem, we hereby present the very first report on viscoelastic properties of wild type human lamin A and some of its mutants linked with Dilated cardiomyopathy (DCM) using quantitative rheological measurements. We observed a dramatic strain-softening effect on lamin A network as an outcome of the strain amplitude sweep measurements which could arise from the large compliance of the quasi-cross-links in the network or that of the lamin A rods. In addition, the drastic stiffening of the differential elastic moduli on superposition of rotational and oscillatory shear stress reflect the increase in the stiffness of the laterally associated lamin A rods. These findings present a preliminary insight into distinct biomechanical properties of wild type lamin A protein and its mutants which in turn revealed interesting differences.

Common recognition principles across diverse sequence and structural families of sialic acid binding proteins

Sialic acids form a large family of 9-carbon monosaccharides and are integral components of glycoconjugates. They are known to bind to a wide range of receptors belonging to diverse sequence families and fold classes and are key mediators in a plethora of cellular processes. Thus, it is of great interest to understand the features that give rise to such a recognition capability. Structural analyses using a non-redundant data set of known sialic acid binding proteins was carried out, which included exhaustive binding site comparisons and site alignments using in-house algorithms, followed by clustering and tree computation, which has led to derivation of sialic acid recognition principles. Although the proteins in the data set belong to several sequence and structure families, their binding sites could be grouped into only six types. Structural comparison of the binding sites indicates that all sites contain one or more different combinations of key structural features over a common scaffold. The six binding site types thus serve as structural motifs for recognizing sialic acid. Scanning the motifs against a non-redundant set of binding sites from PDB indicated the motifs to be specific for sialic acid recognition. Knowledge of determinants obtained from this study will be useful for detecting function in unknown proteins. As an example analysis, a genome-wide scan for the motifs in structures of Mycobacterium tuberculosis proteome identified 17 hits that contain combinations of the features, suggesting a possible function of sialic acid binding by these proteins.

Flocculation behaviour of hematite-kaolinite suspensions in presence of extracellular bacterial proteins and polysaccharides

Cells of Bacillus subtilis exhibited higher affinity towards hematite than to kaolinite. Bacterial cells were grown and adapted in the presence of hematite and kaolinite. Higher amounts of mineral-specific proteinaceous compounds were secreted in the presence of kaolinite while hematite-grown cells produced higher amounts of exopolysaccharides. Extracellular proteins (EP) exhibited higher adsorption density on kaolinite which was rendered more hydrophobic. Hematite surfaces were rendered more hydrophilic due to increased adsorption of extracellular polysaccharides (ECP). Significant surface chemical changes were produced due to interaction between minerals and extracellular proteins and polysaccharides. Iron oxides such as hematite could be effectively removed from kaolinite clays using selective bioflocculation of hematite after interaction with EP and ECP extracted from mineral-grown cells. (C) 2013 Elsevier B.V. All rights reserved.

The relationship between classification of multi-domain proteins using an alignment-free approach and their functions: a case study with immunoglobulins

Establishing functional relationships between multi-domain protein sequences is a non-trivial task. Traditionally, delineating functional assignment and relationships of proteins requires domain assignments as a prerequisite. This process is sensitive to alignment quality and domain definitions. In multi-domain proteins due to multiple reasons, the quality of alignments is poor. We report the correspondence between the classification of proteins represented as full-length gene products and their functions. Our approach differs fundamentally from traditional methods in not performing the classification at the level of domains. Our method is based on an alignment free local matching scores (LMS) computation at the amino-acid sequence level followed by hierarchical clustering. As there are no gold standards for full-length protein sequence classification, we resorted to Gene Ontology and domain-architecture based similarity measures to assess our classification. The final clusters obtained using LMS show high functional and domain architectural similarities. Comparison of the current method with alignment based approaches at both domain and full-length protein showed superiority of the LMS scores. Using this method we have recreated objective relationships among different protein kinase sub-families and also classified immunoglobulin containing proteins where sub-family definitions do not exist currently. This method can be applied to any set of protein sequences and hence will be instrumental in analysis of large numbers of full-length protein sequences.

Correlation between local structural dynamics of proteins inferred from NMR ensembles and evolutionary dynamics of homologues of known structure

Conformational changes in proteins are extremely important for their biochemical functions. Correlation between inherent conformational variations in a protein and conformational differences in its homologues of known structure is still unclear. In this study, we have used a structural alphabet called Protein Blocks (PBs). PBs are used to perform abstraction of protein 3-D structures into a 1-D strings of 16 alphabets (a-p) based on dihedral angles of overlapping pentapeptides. We have analyzed the variations in local conformations in terms of PBs represented in the ensembles of 801 protein structures determined using NMR spectroscopy. In the analysis of concatenated data over all the residues in all the NMR ensembles, we observe that the overall nature of inherent local structural variations in NMR ensembles is similar to the nature of local structural differences in homologous proteins with a high correlation coefficient of .94. High correlation at the alignment positions corresponding to helical and beta-sheet regions is only expected. However, the correlation coefficient by considering only the loop regions is also quite high (.91). Surprisingly, segregated position-wise analysis shows that this high correlation does not hold true to loop regions at the structurally equivalent positions in NMR ensembles and their homologues of known structure. This suggests that the general nature of local structural changes is unique; however most of the local structural variations in loop regions of NMR ensembles do not correlate to their local structural differences at structurally equivalent positions in homologues.

Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety-six of these alignment positions correspond to conserved Gly residue with (phi, ) values allowed for non-glycyl residues. Reasons for this observation were investigated by in-silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (phi, ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (phi, ) values which are disallowed for Ala. In-silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.

PROPENSITIES OF AMINO ACID RESIDUES IN PROTEINS FOR DIFFERENT REGIONS OF THE RAMACHANDRAN MAP

Using a dataset of 1164 crystal structures of largely non-homologous proteins defined at a resolution of 1.5 angstrom or better, we have investigated the (phi,psi) preferences of 20 residue types by considering the residues which occur in loops. Propensities of residue types to occur in the loops with (phi,psi) values in the aa region of the Ramachandran map has a poor correlation coefficient of 0.48 to the Chou-Fasman propensities of the residue types to occur in the a-helical segments. However the correlation coefficient between propensities of residues in loops to adopt beta conformations and those in beta-sheet is much higher (0.95). These observations suggest that a-helix formation is well influenced by the local amino acid sequence while intrinsic preference of residue types for beta-sheet plays a major role in the formation of beta-sheet. The main chain polar groups of residues in loops, that can affect the (phi,psi) values, can be involved in intra-molecular hydrogen bonding. Therefore we investigated further by considering subset of residues in loops with low (0 to 2) number of intra-molecular hydrogen bonds per residue involving main chain polar atoms. For this subset, the correlation coefficients between propensities for alpha-helix and alpha(R) region and between beta-sheet and beta-region are 0.26 and 0.64 respectively. This reiterates higher intrinsic tendency of beta-region favouring residues to adopt beta-sheet than alpha(R) region favouring residues to adopt alpha-helical structure.

Accurate prediction of interfacial residues in two-domain proteins using evolutionary information: Implications for three-dimensional modeling

With the preponderance of multidomain proteins in eukaryotic genomes, it is essential to recognize the constituent domains and their functions. Often function involves communications across the domain interfaces, and the knowledge of the interacting sites is essential to our understanding of the structure-function relationship. Using evolutionary information extracted from homologous domains in at least two diverse domain architectures (single and multidomain), we predict the interface residues corresponding to domains from the two-domain proteins. We also use information from the three-dimensional structures of individual domains of two-domain proteins to train naive Bayes classifier model to predict the interfacial residues. Our predictions are highly accurate (approximate to 85%) and specific (approximate to 95%) to the domain-domain interfaces. This method is specific to multidomain proteins which contain domains in at least more than one protein architectural context. Using predicted residues to constrain domain-domain interaction, rigid-body docking was able to provide us with accurate full-length protein structures with correct orientation of domains. We believe that these results can be of considerable interest toward rational protein and interaction design, apart from providing us with valuable information on the nature of interactions. Proteins 2014; 82:1219-1234. (c) 2013 Wiley Periodicals, Inc.