A network representation of protein structures: Implications for protein stability

This study views each protein structure as a network of noncovalent connections between amino acid side chains. Each amino acid in a protein structure is a node, and the strength of the noncovalent interactions between two amino acids is evaluated for edge determination. The protein structure graphs (PSGs) for 232 proteins have been constructed as a function of the cutoff of the amino acid interaction strength at a few carefully chosen values. Analysis of such PSGs constructed on the basis of edge weights has shown the following: 1), The PSGs exhibit a complex topological network behavior, which is dependent on the interaction cutoff chosen for PSG construction. 2), A transition is observed at a critical interaction cutoff, in all the proteins, as monitored by the size of the largest cluster (giant component) in the graph. Amazingly, this transition occurs within a narrow range of interaction cutoff for all the proteins, irrespective of the size or the fold topology. And 3), the amino acid preferences to be highly connected (hub frequency) have been evaluated as a function of the interaction cutoff. We observe that the aromatic residues along with arginine, histidine, and methionine act as strong hubs at high interaction cutoffs, whereas the hydrophobic leucine and isoleucine residues get added to these hubs at low interaction cutoffs, forming weak hubs. The hubs identified are found to play a role in bringing together different secondary structural elements in the tertiary structure of the proteins. They are also found to contribute to the additional stability of the thermophilic proteins when compared to their mesophilic counterparts and hence could be crucial for the folding and stability of the unique three-dimensional structure of proteins. Based on these results, we also predict a few residues in the thermophilic and mesophilic proteins that can be mutated to alter their thermal stability.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Polypyrimidine tract-binding protein interacts with coxsackievirus B3 RNA and influences its translation

We have investigated the possible role of trans-acting factors interacting with the untranslated regions (UTRs) of coxsackievirus B3 (CVB3) RNA. We show here that polypyrimidine tract-binding protein (PTB) binds specifically to both 5' and 3' UTRs, but with different affinity. We have demonstrated that PTB is a bona fide internal ribosome entry site (IRES) trans-acting factor (ITAF) for CVB3 RNA by characterizing the effect of partial silencing of FIB ex vivo in He La cells. Furthermore, IRES activity in BSC-1 cells, which are reported to have a very low level of endogenous FIB, was found to be significantly lower than that in He La cells. Additionally, we have mapped the putative contact points of PTB on the 5' and 3' UTRs by an RNA toe-printing assay. We have shown that the 3' UTR is able to stimulate CVB3 IRES-mediated translation. Interestingly, a deletion of 15 nt at the 5' end or 14 rut at the 3' end of the CVB3 3' UTR reduced the 3' UTR-mediated enhancement of IRES activity ex vivo significantly, and a reduced interaction was shown with PTB. It appears that the FIB protein might help in circularization of the CVB3 RNA by bridging the ends necessary for efficient translation of the viral RNA.

Molecular Basis for the Role of Staphylococcus aureus Penicillin Binding Protein 4 in Antimicrobial Resistance

Penicillin binding proteins (PBPs) are membrane-associated proteins that catalyze the final step of murein biosynthesis. These proteins function as either transpeptidases or carboxypeptidases and in a few cases demonstrate transglycosylase activity. Both transpeptidase and carboxypeptidase activities of PBPs occur at the D-Ala-D-Ala terminus of a murein precursor containing a disaccharide pentapeptide comprising N-acetyl-glucosamine and N-acetyl-muramic acid-L-Ala-D-Glu-L-Lys-D-Ala-D-Ala. beta-Lactam antibiotics inhibit these enzymes by competing with the pentapeptide precursor for binding to the active site of the enzyme. Here we describe the crystal structure, biochemical characteristics, and expression profile of PBP4, a low-molecular-mass PBP from Staphylococcus aureus strain COL. The crystal structures of PBP4-antibiotic complexes reported here were determined by molecular replacement, using the atomic coordinates deposited by the New York Structural Genomics Consortium. While the pbp4 gene is not essential for the viability of S. aureus, the knockout phenotype of this gene is characterized by a marked reduction in cross-linked muropeptide and increased vancomycin resistance. Unlike other PBPs, we note that expression of PBP4 was not substantially altered under different experimental conditions, nor did it change across representative hospital- or community-associated strains of S. aureus that were examined. In vitro data on purified recombinant S. aureus PBP4 suggest that it is a beta-lactamase and is not trapped as an acyl intermediate with beta-lactam antibiotics. Put together, the expression analysis and biochemical features of PBP4 provide a framework for understanding the function of this protein in S. aureus and its role in antimicrobial resistance.

Mycobacterium tuberculosis nucleoid-associated DNA-binding protein H-NS binds with high-affinity to the Holliday junction and inhibits strand exchange promoted by RecA protein

A number of studies have shown that the structure and composition of bacterial nucleoid influences many a processes related to DNA metabolism. The nucleoid-associated proteins modulate not only the DNA conformation but also regulate the DNA metabolic processes such as replication, recombination, repair and transcription. Understanding of how these processes occur in the context of Mycobacterium tuberculosis nucleoid is of considerable medical importance because the nucleoid structure may be constantly remodeled in response to environmental signals and/or growth conditions. Many studies have concluded that Escherichia coli H-NS binds to DNA in a sequence-independent manner, with a preference for A-/T-rich tracts in curved DNA; however, recent studies have identified the existence of medium- and low-affinity binding sites in the vicinity of the curved DNA. Here, we show that the M. tuberculosis H-NS protein binds in a more structure-specific manner to DNA replication and repair intermediates, but displays lower affinity for double-stranded DNA with relatively higher GC content. Notably, M. tuberculosis H-NS was able to bind Holliday junction (HJ), the central recombination intermediate, with substantially higher affinity and inhibited the three-strand exchange promoted by its cognate RecA. Likewise, E. coli H-NS was able to bind the HJ and suppress DNA strand exchange promoted by E. coli RecA, although much less efficiently compared to M. tuberculosis H-NS. Our results provide new insights into a previously unrecognized function of H-NS protein, with implications for blocking the genome integration of horizontally transferred genes by homologous and/or homeologous recombination.

HARMONY: a server for the assessment of protein structures

Protein structure validation is an important step in computational modeling and structure determination. Stereochemical assessment of protein structures examine internal parameters such as bond lengths and Ramachandran (phi, psi) angles. Gross structure prediction methods such as inverse folding procedure and structure determination especially at low resolution can sometimes give rise to models that are incorrect due to assignment of misfolds or mistracing of electron density maps. Such errors are not reflected as strain in internal parameters. HARMONY is a procedure that examines the compatibility between the sequence and the structure of a protein by assigning scores to individual residues and their amino acid exchange patterns after considering their local environments. Local environments are described by the backbone conformation, solvent accessibility and hydrogen bonding patterns. We are now providing HARMONY through a web server such that users can submit their protein structure files and, if required, the alignment of homologous sequences. Scores are mapped on the structure for subsequent examination that is useful to also recognize regions of possible local errors in protein structures. HARMONY server is located at http://caps.ncbs.res.in/harmony/

Computer modelling approach to study the modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein

The modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein (ABP) have been studied by energy minimization using the low resolution (2.4 A) X-ray data of the protein. These studies suggest that these sugars preferentially bind in the alpha-form to ABP, unlike L-arabinose where both alpha- and beta-anomers bind almost equally. The best modes of binding of alpha- and beta-anomers of D-galactose and D-fucose differ slightly in the nature of the possible hydrogen bonds with the protein. The residues Arg 151 and Asn 232 of ABP from bidentate hydrogen bonds with both L-arabinose and D-galactose, but not with D-fucose or D-glucose. However in the case of L-arabinose, Arg 151 forms hydrogen bonds with the hydroxyl group at the C-4 atom and the ring oxygen, whereas in case of D-galactose it forms bonds with the hydroxyl groups at the C-4 and C-6 atoms of the pyranose ring. The calculated conformational energies also predict that D-galactose is a better inhibitor than D-fucose and D-glucose, in agreement with kinetic studies. The weak inhibitor D-glucose binds preferentially to one domain of ABP leading to the formation of a weaker complex. Thus these studies provide information about the most probable binding modes of these sugars and also provide a theoretical explanation for the observed differences in their binding affinities.

Computer simulation of protein—carbohydrate complexes: application to arabinose-binding protein and pea lectin

The CCEM method (Contact Criteria and Energy Minimisation) has been developed and applied to study protein-carbohydrate interactions. The method uses available X-ray data even on the native protein at low resolution (above 2.4 Å) to generate realistic models of a variety of proteins with various ligands.The two examples discussed in this paper are arabinose-binding protein (ABP) and pea lectin. The X-ray crystal structure data reported on ABP-β-l-arabinose complex at 2.8, 2.4 and 1.7 Å resolution differ drastically in predicting the nature of the interactions between the protein and ligand. It is shown that, using the data at 2.4 Å resolution, the CCEM method generates complexes which are as good as the higher (1.7 Å) resolution data. The CCEM method predicts some of the important hydrogen bonds between the ligand and the protein which are missing in the interpretation of the X-ray data at 2.4 Å resolution. The theoretically predicted hydrogen bonds are in good agreement with those reported at 1.7 Å resolution. Pea lectin has been solved only in the native form at 3 Å resolution. Application of the CCEM method also enables us to generate complexes of pea lectin with methyl-α-d-glucopyranoside and methyl-2,3-dimethyl-α-d-glucopyranoside which explain well the available experimental data in solution.

Vitamin A and thyroxine carrier proteins in chicken plasma Steady-state control of the plasma level of free retinol-binding protein and free thyroxine

1. 1. The binding parameters of prealbumin-2 with retinol-binding protein and thyroxine (T4) revealed the existence of distinct and multiple sites for both retinol-binding protein and T4. 2. 2. From the analysis of binding parameters of retinol-binding protein with prealbumin-2 it is clear that under steady-state conditions about 99% of the holo-retinol-binding protein remains bound to prealbumin-2. 3. 3. Equilibrium dialysis studies on binding properties of thyroid hormones with prealbumin-2 revealed that it has a single high affinity site and three low affinity sites. 4. 4. The occurrence of three carrier proteins for thyroid hormones, thyroxine-binding globulin, prealbumin-2 and albumin has been demonstrated. However, the chicken thyroxine-binding globulin differs from human thyroxine-binding globulin by being relatively less acidic and occuring at a two-fold lower concentration. But the thyroid hormone binding parameters are comparable. 5. 5. Highly sensitive methods were developed for determination of T4 binding capacities of the various proteins and plasma level of total T4 by fractionation of carrier proteins and further quantitatively employing in electrophoresis and equilibrium dialysis. 6. 6. The thyroxine-binding proteins were found to be two types, one (viz., thyroxine-binding globulin) of great affinity but of low binding capacity, which mainly acts as reservoir of T4, and another (viz.,prealbumin-2) of low affinity but of high binding capacity, which can participate predominantly in the control of the free T4 pool.

Microsomal redox systems in brown adipose tissue: high lipid peroxidation, low cholesterol biosynthesis and no detectable cytochrome P-450

The presence of redox systems in microsomes of brown adipose tissue (BAT) in cold exposed rats was investigated and compared with liver. BAT microsomes showed high activity of lipid peroxidation measured both by the formation of malondialdehyde (MDA) and by oxygen uptake. NADH and NADPH dependent cytochrome c reductase activity were present in both BAT and liver microsomes. Aminopyrine demethylase and aniline hydroxylase activities, the characteristic detoxification enzymes in liver microsomes could not be detected in BAT microsomes. BAT minces showed very poor incorporation of [1-14C]acetate and [2-14C]-mevalonate in unsaponifiable lipid fraction compared to liver. Biosynthesis of cholesterol and ubiquinone, but not fatty acids, and the activity of 3-hydroxy-3-methyl glutaryl CoA reductase appear to be very low in BAT. Examination of difference spectra showed the presence of only cytochrome b 5 in BAT microsomes. In addition to the inability to detect the enzyme activities dependent on cytochrome P-450, a protein with the characteristic spectrum, molecular size in SDS-PAGE and interaction with antibodies in double diffusion test, also could not be detected in BAT microsomes. The high activity of lipid peroxidation in microsomes, being associated with large oxygen uptake and oxidation of NADPH, will also contribute to the energy dissipation as heat in BAT, considered important in thermogenesis.

Structure of Sesbania Mosaic Virus at 4·7 Å Resolution and Partial Sequence of the Coat Protein

Sesbania mosaic virus (SMV) is a plant virus infecting Sesbania grandiflora plants in Andhra Pradesh, India. Amino acid sequence of the tryptic peptides of SMV coat protein were determined using a gas phase sequenator. These sequences showed identical amino acids at 69% of the positions when aligned with the corresponding residues of southern bean mosaic virus (SBMV).Crystals diffracting to better than 3 Å resolution were obtained by precipitating the virus with ammonium sulphate. The crystals belonged to rhombohedral space group R3 with α = 291·4 Å and α = 61·9°. Three-dimensional X-ray diffraction data on these crystals were collected to a resolution of 4·7 Å, using a Siemens-Nicolet area detector system. Self-rotation function studies revealed the icosahedral symmetry of the virus particles, as well as their precise orientation in the unit cell. Cross-rotation function and modelling studies with SBMV showed that it is a valid starting model for SMV structure determination. Low resolution phases computed using a polyalanine model of SBMV were subjected to refinement and extension by real-space electron density averaging and solvent flattening. The final electron density map revealed a polypeptide fold similar to SBMV. The single disulphide bridge of SBMV coat protein is retained in SMV. Four icosahedrally independent cation binding sites have been tentatively identified. Three of these sites, related by a quasi threefold axis, are also found in SBMV. The fourth site is situated on the quasi threefold axis. Aspartic acid residues, which replace Ile218 of SBMV from the quasi threefold-related subunits are suitable ligands to the cation at this site

The variants of the protein toxins abrin and ricin A useful guide to understanding the processing events in the toxin transport

Kinetic data on inhibition of protein synthesis in thymocyte by three abrins and ricin have been obtained. The intrinsic efficiencies of A chains of four toxins to inactivate ribosomes, as analyzed by k1-versus-concentration plots were abrin II, III > ricin > abrin I. The lag times were 90, 66, 75 and 105 min at a 0.0744 nM concentration of each of abrin I, II, III and ricin, respectively. To account for the observed differences in the dose-dependent lag time, functional and structural variables of toxins such as binding efficiency of B chains to receptors and low-pH-induced structural alterations have been analyzed. The association constants obtained by stopped flow studies showed that abrin-I (4.13 × 105 M−1 s−1) association with putative receptor (4-methylumbelliferyl-α-D-galactoside) is nearly two times more often than abrin III (2.6 × 105 M−1 s−1) at 20°C. Equillibrium binding constants of abrin I and II to thymocyte at 37°C were 2.26 × 107 M−1 and 2.8 × 107 M−1 respectively. pH-induced structural alterations as studied by a parallel enhancement in 8-anilino-L-naphthalene sulfonate fluorescence revealed a high degree of qualitative similarity. These results taken with a nearly identical concentration-independent lag time (minimum lag of 41–42 min) indicated that the binding efficiencies and internalization efficiencies of these toxins are the same and that the observed difference in the dose-dependent lag time is causally related to the proposed processing event. The rates of reduction of inter-subunit disulfide bond, an obligatory step in the intoxication process, have been measured and compared under a variety of conditions. Intersubunit disulfide reduction of abrin I is fourfold faster than that of abrin II at pH 7.2. The rate of disulfide reduction in abrin I could be decreased 1 I-fold by adding lactose, compared to that without lactose. The observed differences in the efficiencies of A chains, the dose-dependent lag period, the modulating effect of lactose on the rates of disulfide reduction and similarity in binding properties make the variants a valuable tool to probe the processing events in toxin transport in detail.

Characterization of lysozyme crystals with unusually low solvent content

Studies on the low-humidity (88%) forms of tetragonal and monoclinic lysozyme, resulting from water-mediated transformations, have provided a wealth of information on the variability in protein hydration, its structural consequences and the water structure associated with proteins, in addition to facilitating the delineation of the rigid and the flexible regions in the protein molecule and the invariant features in its hydration shell. Surprisingly, monoclinic lysozyme continues to diffract even when the environmental humidity is drastically reduced, thus permitting the structural study of the enzyme at different levels of hydration. As part of a study in this direction, three very low humidity forms, two of them occuring at a nominal relative humidity of 38% and the other at 5% relative humidity, have been characterized. These have unprecedented low solvent contents of 16.9, 17.6 and 9.4%, respectively, as determined by the Matthews method.

Heat shock protein 70 as a biomarker of heat stress in a simulated hot cockpit

Background: Fighter pilots are frequently exposed to high temperatures during high-speed low-level flight. Heat strain can result in temporary impairment of cognitive functions and when severe, loss of consciousness and consequent loss of life and equipment. Induction of stress proteins is a highly conserved stress response mechanism from bacteria to humans. induced stress protein levels are known to be cytoprotective and have been correlated with stress tolerance. Although many studies on the heat shock response mechanisms have been performed in cell culture and animal model systems, there is very limited information on stress protein induction in human subjects. Hypothesis: Heat shock proteins (Hsp), especially Hsp70, may be induced in human subjects exposed to high temperatures in a hot cockpit designed to simulate heat stress experienced in low flying sorties. Methods: Six healthy volunteers were subjected to heat stress at 55degreesC in a high temperature cockpit simulator for a period of 1 h at 30% humidity. Physiological parameters such as oral and skin temperatures, heart rate, and sweat rate were monitored regularly during this time. The level of Hsp70 in leukocytes was examined before and after the heat exposure in each subject. Conclusions: Hsp70 was found to be significantly induced in all the six subjects exposed to heat stress. The level of induced Hsp70 appears to correlate with other strain indicators such as accumulative circulatory strain and Craig's modified index. The usefulness of Hsp70 as a molecular marker of heat stress in humans is discussed.