Physical understanding of pore formation on supported lipid bilayer by bacterial toxins

Pore forming toxins are being classified in the protein community based on their ability of forming pores in living cell membranes. Some initial study has apparently pointed out the crystallographic pathway rather can be viewed as a structural as well as morphological changes of proteins in terms of self assembly before and during the pore formation process in surfactant medium. Being a water soluble compound, it changes its conformation and originates some pre-pore complex, which later partially goes inside the cell membrane causing a pore. The physical mechanism for this whole process is still unknown. In this study we have tried to understand these types of biological processes from physical point of view by using supported lipid bilayer as a model system.

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) C alpha, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

In addition to the biologically active monomer of the protein insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The insulin monomer contains two distinct surfaces, namely, the dimer forming surface (DFS) and the hexamer forming surface (HFS), that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behavior of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. In particular, we find that more structured water molecules, with higher residence times (similar to 300-500 ps), are present near HFS than those near DFS. A significant slowing down is observed in the decay of associated rotational auto time correlation functions of O-H bond vector of water in the vicinity of HFS. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly an outward H-atoms directed, clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inward H-atoms directed orientation owing to the flat curvature of hydrophobic surface and the interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces that reveal the significant differences between the two hydration layers.

PLIC: protein-ligand interaction clusters

Most of the biological processes are governed through specific protein-ligand interactions. Discerning different components that contribute toward a favorable protein-ligand interaction could contribute significantly toward better understanding protein function, rationalizing drug design and obtaining design principles for protein engineering. The Protein Data Bank (PDB) currently hosts the structure of similar to 68 000 protein-ligand complexes. Although several databases exist that classify proteins according to sequence and structure, a mere handful of them annotate and classify protein-ligand interactions and provide information on different attributes of molecular recognition. In this study, an exhaustive comparison of all the biologically relevant ligand-binding sites (84 846 sites) has been conducted using PocketMatch: a rapid, parallel, in-house algorithm. PocketMatch quantifies the similarity between binding sites based on structural descriptors and residue attributes. A similarity network was constructed using binding sites whose PocketMatch scores exceeded a high similarity threshold (0.80). The binding site similarity network was clustered into discrete sets of similar sites using the Markov clustering (MCL) algorithm. Furthermore, various computational tools have been used to study different attributes of interactions within the individual clusters. The attributes can be roughly divided into (i) binding site characteristics including pocket shape, nature of residues and interaction profiles with different kinds of atomic probes, (ii) atomic contacts consisting of various types of polar, hydrophobic and aromatic contacts along with binding site water molecules that could play crucial roles in protein-ligand interactions and (iii) binding energetics involved in interactions derived from scoring functions developed for docking. For each ligand-binding site in each protein in the PDB, site similarity information, clusters they belong to and description of site attributes are provided as a relational database-protein-ligand interaction clusters (PLIC).

PocketMatch (version 2.0): A parallel algorithm for the detection of structural similarities between protein ligand binding-sites

Knowledge of protein-ligand interactions is essential to understand several biological processes and important for applications ranging from understanding protein function to drug discovery and protein engineering. Here, we describe an algorithm for the comparison of three-dimensional ligand-binding sites in protein structures. A previously described algorithm, PocketMatch (version 1.0) is optimised, expanded, and MPI-enabled for parallel execution. PocketMatch (version 2.0) rapidly quantifies binding-site similarity based on structural descriptors such as residue nature and interatomic distances. Atomic-scale alignments may also be obtained from amino acid residue pairings generated. It allows an end-user to compute database-wide, all-to-all comparisons in a matter of hours. The use of our algorithm on a sample dataset, performance-analysis, and annotated source code is also included.

Sparsely populated residue conformations in protein structures: Revisiting ``experimental'' Ramachandran maps

The Ramachandran map clearly delineates the regions of accessible conformational (phi-) space for amino acid residues in proteins. Experimental distributions of phi, values in high-resolution protein structures, reveal sparsely populated zones within fully allowed regions and distinct clusters in apparently disallowed regions. Conformational space has been divided into 14 distinct bins. Residues adopting these relatively rare conformations are presented and amino acid propensities for these regions are estimated. Inspection of specific examples in a completely arid, fully allowed region in the top left quadrant establishes that side-chain and backbone interactions may provide the energetic compensation necessary for populating this region of phi- space. Asn, Asp, and His residues showed the highest propensities in this region. The two distinct clusters in the bottom right quadrant which are formally disallowed on strict steric considerations correspond to the gamma turn (C7 axial) conformation (Bin 12) and the i + 1 position of Type II turns (Bin 13). Of the 516 non-Gly residues in Bin 13, 384 occupied the i + 1 position of Type II turns. Further examination of these turn segments revealed a high propensity to occur at the N-terminus of helices and as a tight turn in hairpins. The strand-helix motif with the Type II turn as a connecting element was also found in as many as 57 examples. Proteins 2014; 82:1101-1112. (c) 2013 Wiley Periodicals, Inc.

Direct control of type IIA topoisomerase activity by a chromosomally encoded regulatory protein

Precise control of supercoiling homeostasis is critical to DNA-dependent processes such as gene expression, replication, and damage response. Topoisomerases are central regulators of DNA supercoiling commonly thought to act independently in the recognition and modulation of chromosome superstructure; however, recent evidence has indicated that cells tightly regulate topoisomerase activity to support chromosome dynamics, transcriptional response, and replicative events. How topoisomerase control is executed and linked to the internal status of a cell is poorly understood. To investigate these connections, we determined the structure of Escherichia coil gyrase, a type HA topoisomerase bound to YacG, a recently identified chromosomally encoded inhibitor protein. Phylogenetic analyses indicate that YacG is frequently associated with coenzyme A (CoA) production enzymes, linking the protein to metabolism and stress. The structure, along with supporting solution studies, shows that YacG represses gyrase by sterically occluding the principal DNA-binding site of the enzyme. Unexpectedly, YacG acts by both engaging two spatially segregated regions associated with small-molecule inhibitor interactions (fluoroquinolone antibiotics and the newly reported antagonist GSK299423) and remodeling the gyrase holo enzyme into an inactive, ATP-trapped configuration. This study establishes a new mechanism for the protein-based control of topoisomerases, an approach that may be used to alter supercoiling levels for responding to changes in cellular state.

Mycobacterium tuberculosis Cell Division Protein, FtsE, is an ATPase in Dimeric Form

FtsE is one of the earliest cell division proteins that assembles along with FtsX at the mid-cell site during cell division in Escherichia coli. Both these proteins are highly conserved across diverse bacterial genera and are predicted to constitute an ABC transporter type complex, in which FtsE is predicted to bind ATP and hydrolyse it, and FtsX is predicted to be an integral membrane protein. We had earlier reported that the MtFtsE of the human pathogen, Mycobacterium tuberculosis, binds ATP and interacts with MtFtsX on the cell membrane of M. tuberculosis and E. coli. In this study, we demonstrate that MtFtsE is an ATPase, the active form of which is a dimer, wherein the participating monomers are held together by non-covalent interactions, with the Cys84 of each monomer present at the dimer interface. Under oxidising environment, the dimer gets stabilised by the formation of Cys84-Cys84 disulphide bond. While the recombinant MtFtsE forms a dimer on the membrane of E. coli, the native MtFtsE seems to be in a different conformation in the M. tuberculosis membrane. Although disulphide bridges were not observed on the cytoplasmic side (reducing environment) of the membrane, the two participating monomers could be isolated as dimers held together by non-covalent interactions. Taken together, these findings show that MtFtsE is an ATPase in the non-covalent dimer form, with the Cys84 of each monomer present in the reduced form at the dimer interface, without participating in the dimerisation or the catalytic activity of the protein.

Jacalin-carbohydrate interactions: distortion of the ligand molecule as a determinant of affinity

Jacalin is among the most thoroughly studied lectins. Its carbohydrate-binding site has also been well characterized. It has been postulated that the lower affinity of beta-galactosides for jacalin compared with beta-galactosides is caused by steric interactions of the substituents in the former with the protein. This issue has been explored energetically and structurally using different appropriate carbohydrate complexes of jacalin. It turns out that the earlier postulation is not correct. The interactions of the substituent with the binding site remain essentially the same irrespective of the anomeric nature of the substitution. This is achieved through a distortion of the sugar ring in beta-galactosides. The difference in energy, and therefore in affinity, is caused by a distortion of the sugar ring in beta-galactosides. The elucidation of this unprecedented distortion of the ligand as a strategy for modulating affinity is of general interest. The crystal structures also provide a rationale for the relative affinities of the different carbohydrate ligands for jacalin.

Structural and functional analysis of two universal stress proteins YdaA and YnaF from Salmonella typhimurium: possible roles in microbial stress tolerance

In many organisms ``Universal Stress Proteins'' CUSPS) are induced in response to a variety of environmental stresses. Here we report the structures of two USPs, YnaF and YdaA from Salmonella typhimurium determined at 1.8 angstrom and 2.4 angstrom resolutions, respectively. YnaF consists of a single USP domain and forms a tetrameric organization stabilized by interactions mediated through chloride ions. YdaA is a larger protein consisting of two tandem USP domains. Two protomers of YdaA associate to form a structure similar to the YnaF tetramer. YdaA showed ATPase activity and an ATP binding motif G-2X-G-9X-G(S/T/N) was found in its C-terminal domain. The residues corresponding to this motif were not conserved in YnaF although YnaF could bind ATP. However, unlike YdaA, YnaF did not hydrolyse ATP in vitro. Disruption of interactions mediated through chloride ions by selected mutations converted YnaF into an ATPase. Residues that might be important for ATP hydrolysis could be identified by comparing the active sites of native and mutant structures. Only the C-terminal domain of YdaA appears to be involved in ATP hydrolysis. The structurally similar N-terminal domain was found to bind a zinc ion near the segment equivalent to the phosphate binding loop of the C-terminal domain. Mass spectrometric analysis showed that YdaA might bind a ligand of approximate molecular weight 800 daltons. Structural comparisons suggest that the ligand, probably related to an intermediate in lipid A biosynthesis, might bind at a site close to the zinc ion. Therefore, the N-terminal domain of YdaA binds zinc and might play a role in lipid metabolism. Thus, USPs appear to perform several distinct functions such as ATP hydrolysis, altering membrane properties and chloride sensing. (C) 2015 Elsevier Inc. All rights reserved.

A Gaussian network model study suggests that structural fluctuations are higher for inactive states than active states of protein kinases

We performed Gaussian network model based normal mode analysis of 3-dimensional structures of multiple active and inactive forms of protein kinases. In 14 different kinases, a more number of residues (1095) show higher structural fluctuations in inactive states than those in active states (525), suggesting that, in general, mobility of inactive states is higher than active states. This statistically significant difference is consistent with higher crystallographic B-factors and conformational energies for inactive than active states, suggesting lower stability of inactive forms. Only a small number of inactive conformations with the DFG motif in the ``in'' state were found to have fluctuation magnitudes comparable to the active conformation. Therefore our study reports for the first time, intrinsic higher structural fluctuation for almost all inactive conformations compared to the active forms. Regions with higher fluctuations in the inactive states are often localized to the aC-helix, aG-helix and activation loop which are involved in the regulation and/or in structural transitions between active and inactive states. Further analysis of 476 kinase structures involved in interactions with another domain/protein showed that many of the regions with higher inactive-state fluctuation correspond to contact interfaces. We also performed extensive GNM analysis of (i) insulin receptor kinase bound to another protein and (ii) holo and apo forms of active and inactive conformations followed by multi-factor analysis of variance. We conclude that binding of small molecules or other domains/proteins reduce the extent of fluctuation irrespective of active or inactive forms. Finally, we show that the perceived fluctuations serve as a useful input to predict the functional state of a kinase.

Calcium binding properties of calcium dependent protein kinase 1 (CaCDPK1) from Cicer arietinum

Calcium plays a crucial role as a secondary messenger in all aspects of plant growth, development and survival. Calcium dependent protein kinases (CDPKs) are the major calcium decoders, which couple the changes in calcium level to an appropriate physiological response. The mechanism by which calcium regulates CDPK protein is not well understood. In this study, we investigated the interactions of Ca2+ ions with the CDPK1 isoform of Cicer arietinum (CaCDPK1) using a combination of biophysical tools. CaCDPK1 has four different EF hands as predicted by protein sequence analysis. The fluorescence emission spectrum of CaCDPK1 showed quenching with a 5 nm red shift upon addition of calcium, indicating conformational changes in the tertiary structure. The plot of changes in intensity against calcium concentrations showed a biphasic curve with binding constants of 1.29 mu M and 120 mu M indicating two kinds of binding sites. Isothermal calorimetric (ITC) titration with CaCl2 also showed a biphasic curve with two binding constants of 0.027 mu M and 1.7 mu M. Circular dichroism (CD) spectra showed two prominent peaks at 208 and 222 nm indicating that CaCDPK1 is a alpha-helical rich protein. Calcium binding further increased the alpha-helical content of CaCDPK1 from 75 to 81%. Addition of calcium to CaCDPK1 also increased fluorescence of 8-anilinonaphthalene-1-sulfonic acid (ANS) indicating exposure of hydrophobic surfaces. Thus, on the whole this study provides evidence for calcium induced conformational changes, exposure of hydrophobic surfaces and heterogeneity of EF hands in CaCDPK1. (C) 2015 Elsevier GmbH. All rights reserved.

Interplay of electrostatics and lipid packing determines the binding of charged polymer coated nanoparticles to model membranes

Understanding of nanoparticle-membrane interactions is useful for various applications of nanoparticles like drug delivery and imaging. Here we report on the studies of interaction between hydrophilic charged polymer coated semiconductor quantum dot nanoparticles with model lipid membranes. Atomic force microscopy and X-ray reflectivity measurements suggest that cationic nanoparticles bind and penetrate bilayers of zwitterionic lipids. Penetration and binding depend on the extent of lipid packing and result in the disruption of the lipid bilayer accompanied by enhanced lipid diffusion. On the other hand, anionic nanoparticles show minimal membrane binding although, curiously, their interaction leads to reduction in lipid diffusivity. It is suggested that the enhanced binding of cationic QDs at higher lipid packing can be understood in terms of the effective surface potential of the bilayers which is tunable through membrane lipid packing. Our results bring forth the subtle interplay of membrane lipid packing and electrostatics which determine nanoparticle binding and penetration of model membranes with further implications for real cell membranes.

Structural insights into N-terminal to C-terminal interactions and implications for thermostability of a (beta/alpha)(8)-triosephosphate isomerase barrel enzyme

Although several factors have been suggested to contribute to thermostability, the stabilization strategies used by proteins are still enigmatic. Studies on a recombinant xylanase from Bacilllus sp. NG-27 (RBSX), which has the ubiquitous (beta/alpha)(8)-triosephosphate isomerase barrel fold, showed that just a single mutation, V1L, although not located in any secondary structural element, markedly enhanced the stability from 70 degrees C to 75 degrees C without loss of catalytic activity. Conversely, the V1A mutation at the same position decreased the stability of the enzyme from 70 degrees C to 68 degrees C. To gain structural insights into how a single extreme N-terminus mutation can markedly influence the thermostability of the enzyme, we determined the crystal structure of RBSX and the two mutants. On the basis of computational analysis of their crystal structures, including residue interaction networks, we established a link between N-terminal to C-terminal contacts and RBSX thermostability. Our study reveals that augmenting N-terminal to C-terminal noncovalent interactions is associated with enhancement of the stability of the enzyme. In addition, we discuss several lines of evidence supporting a connection between N-terminal to C-terminal noncovalent interactions and protein stability in different proteins. We propose that the strategy of mutations at the termini could be exploited with a view to modulate stability without compromising enzymatic activity, or in general, protein function in diverse folds where N and C termini are in close proximity. Database The coordinates of RBSX, V1A and V1L have been deposited in the PDB database under the accession numbers 4QCE, 4QCF, and 4QDM, respectively

First Structural View of a Peptide Interacting with the Nucleotide Binding Domain of Heat Shock Protein 90

The involvement of Hsp90 in progression of diseases like cancer, neurological disorders and several pathogen related conditions is well established. Hsp90, therefore, has emerged as an attractive drug target for many of these diseases. Several small molecule inhibitors of Hsp90, such as geldanamycin derivatives, that display antitumor activity, have been developed and are under clinical trials. However, none of these tested inhibitors or drugs are peptide-based compounds. Here we report the first crystal structure of a peptide bound at the ATP binding site of the N-terminal domain of Hsp90. The peptide makes several specific interactions with the binding site residues, which are comparable to those made by the nucleotide and geldanamycin. A modified peptide was designed based on these interactions. Inhibition of ATPase activity of Hsp90 was observed in the presence of the modified peptide. This study provides an alternative approach and a lead peptide molecule for the rational design of effective inhibitors of Hsp90 function.