Insight into the early stages of thermal unfolding of peanut agglutinin by molecular dynamics simulations

Peanut agglutinin is a homotetrameric nonglycosylated protein. The protein has a unique open quaternary structure. Molecular dynamics simulations have been employed follow the atomistic details of its unfolding at different temperatures. The early events of the deoligomerization of the protein have been elucidated in the present study. Simulation trajectories of the monomer as well as those of the tetramer have been compared and the tetramer is found to be substantially more stable than its monomeric counterpart. The tetramer shows retention of most of its.. secondary structure but considerable loss of the tertiary structure at high temperature. e generation of a This observation impies the molten globule-like intermediate in the later stages of deoligomerization. The quaternary structure of the protein has weakened to a large extent, but none of the subunits are separated. In addition, the importance of the metal-binding to the stability of the protein structure has also been investigated. Binding of the metal ions not only enhances the local stability of the metal-ion binding loop, but also imparts a global stability to the overall structure. The dynamics of different interfaces vary significantly as probed through interface clusters. The differences are substantially enhanced at higher temperatures. The dynamics and the stability of the interfaces have been captured mainly by cluster analysis, which has provided detailed information on the thermal deoligomerization of the protein.

Studies on the conformation of amino acids XII. Energy calculations on prolyl residue

The favoured conformations of the prolyl residue have been obtained by calculating their potential energies arising from bond-angle strain, torsion-angle strain, non-bonded and electrostatic interatomic energies. In addition to the five membered ring, the peptide unit at the amino end (with ω = 180°) and the C′ atom at the carboxyl end have been taken into account. It is found that there are two local minima in the configurational space of the parameters defining the conformation, as is actually observed-one (denoted by B) with Cγ displaced on the same side as C′, which is lower in energy than the other (denoted by A) with Cγ displaced on the opposite side of C′. The other four atoms Cδ, N, Cα, Cβ are nearly in a plane. The conformations of minimum energy (for both A and B) have bond angles very close to the mean observed values while the torsion angles are well within the range observed in various structures for each type. Taking into account the fact that the influence of neighbouring molecules in a crystal structure may make the conformation of a molecule different from the minimal one, the ranges of the conformational parameters for which the energy is within 0.6 kcal/mole above the minimum value (called the "most probable range") and within 1.2 kcal/mole (called the "probable range") have been determined. The ranges thus obtained, agree well with observation, and most of the observed data lie within the most probable ranges, although differing appreciably from the conformation of minimum energy. The study has been extended, in a limited way, to the conformation of the ring in the amino acid proline. Since the nitrogen is tetrahedral in this (as contrasted with being planar in the prolyl residue), it is found that any one of the five atoms can be out of plane (either way), with the other four lying nearly in a plane. These correspond to low energy conformations (up to 1.2 kcal/mole above the minimum). One such example, in which the Cα atom is out of plane is known for dl-proline · HCl. It is also shown that in these calculations energies due to bond length distortions can be neglected to a good degree of approximation, provided the 'best' values of the bond lengths for the particular compound are used in the theoretical calculations.

A novel approach for large-scale polypeptide folding based on elastic networks using continuous optimization

We present a new computationally efficient method for large-scale polypeptide folding using coarse-grained elastic networks and gradient-based continuous optimization techniques. The folding is governed by minimization of energy based on Miyazawa–Jernigan contact potentials. Using this method we are able to substantially reduce the computation time on ordinary desktop computers for simulation of polypeptide folding starting from a fully unfolded state. We compare our results with available native state structures from Protein Data Bank (PDB) for a few de-novo proteins and two natural proteins, Ubiquitin and Lysozyme. Based on our simulations we are able to draw the energy landscape for a small de-novo protein, Chignolin. We also use two well known protein structure prediction software, MODELLER and GROMACS to compare our results. In the end, we show how a modification of normal elastic network model can lead to higher accuracy and lower time required for simulation.

The nature of inhibition of 3-hydroxy-3-methylglutaryl CoA reductase by garlic-derived diallyl disulfide

A concentration dependent inhibition of 3-hydroxy-3-methylglutaryl CoA (HMG CoA) reductase was found on preincubation of microsomal preparations with diallyl disulfide, a component of garlic oil. This inhibited state was only partially reversed even with high concentrations of DTT. Glutathione, a naturally occurring reducing thiol agent, was ineffective. The substrate, HMG CoA, but not NADPH, was able to give partial protection for the DTT-dependent, but not glutathione-dependent activity. The garlic-derived diallyl disulfide is the most effective among the sulfides tested for inhibition of HMG CoA reductase. Formation of protein internal disulfides, inaccessible for reduction by thiol agents, but not of protein dimer, is likely to be the cause of this inactivation.

Immunochemical analysis of prolamins in gluten-free foods

People with coeliac disease have to maintain a gluten-free diet, which means excluding wheat, barley and rye prolamin proteins from their diet. Immunochemical methods are used to analyse the harmful proteins and to control the purity of gluten-free foods. In this thesis, the behaviour of prolamins in immunological gluten assays and with different prolamin-specific antibodies was examined. The immunoassays were also used to detect residual rye prolamins in sourdough systems after enzymatic hydrolysis and wheat prolamins after deamidation. The aim was to characterize the ability of the gluten analysis assays to quantify different prolamins in varying matrices in order to improve the accuracy of the assays. Prolamin groups of cereals consist of a complex mixture of proteins that vary in their size and amino acid sequences. Two common characteristics distinguish prolamins from other cereal proteins. Firstly, they are soluble in aqueous alcohols, and secondly, most of the prolamins are mainly formed from repetitive amino acid sequences containing high amounts of proline and glutamine. The diversity among prolamin proteins sets high requirements for their quantification. In the present study, prolamin contents were evaluated using enzyme-linked immunosorbent assays based on ω- and R5 antibodies. In addition, assays based on A1 and G12 antibodies were used to examine the effect of deamidation on prolamin proteins. The prolamin compositions and the cross-reactivity of antibodies with prolamin groups were evaluated with electrophoretic separation and Western blotting. The results of this thesis research demonstrate that the currently used gluten analysis methods are not able to accurately quantify barley prolamins, especially when hydrolysed or mixed in oats. However, more precise results can be obtained when the standard more closely matches the sample proteins, as demonstrated with barley prolamin standards. The study also revealed that all of the harmful prolamins, i.e. wheat, barley and rye prolamins, are most efficiently extracted with 40% 1-propanol containing 1% dithiothreitol at 50 °C. The extractability of barley and rye prolamins was considerably higher with 40% 1-propanol than with 60% ethanol, which is typically used for prolamin extraction. The prolamin levels of rye were lowered by 99.5% from the original levels when an enzyme-active rye-malt sourdough system was used for prolamin degradation. Such extensive degradation of rye prolamins suggest the use of sourdough as a part of gluten-free baking. Deamidation increases the diversity of prolamins and improves their solubility and ability to form structures such as emulsions and foams. Deamidation changes the protein structure, which has consequences for antibody recognition in gluten analysis. According to the resuts of the present work, the analysis methods were not able to quantify wheat gluten after deamidation except at very high concentrations. Consequently, deamidated gluten peptides can exist in food products and remain undetected, and thus cause a risk for people with gluten intolerance. The results of this thesis demonstrate that current gluten analysis methods cannot accurately quantify prolamins in all food matrices. New information on the prolamins of rye and barley in addition to wheat prolamins is also provided in this thesis, which is essential for improving gluten analysis methods so that they can more accurately quantify prolamins from harmful cereals.

Inverse relationship of the dehydrogenase and ADP-ribosylation activities in sodium-nitroprusside-treated glyceraldehyde-3-phosphate dehydrogenase is coincidental

Incubation of glyceraldehyde-3-phosphate dehydrogenase (GAPD) with sodium nitroprusside (SNP) decreased its activity in concentration- and time-dependent fashion in the presence of a thiol compounds, with DTT being more effective than GSH. Both forward and backward reactions were effected. Coinciding with this, HgCl2-sensitive labelling of the protein by [32P]NAD+ also increased, indicating the stimulation of ADP-ribosylation. Treatment with SNP of GAPD samples from rabbit muscle, sheep brain and yeast inactivated the dehydrogenase activity of the three, but only the mammalian proteins showed ADP-ribosylation activity. The SNP-modified protein of rabbit muscle GAPD, freed from the reagent by Sephadex filtration showed a concentration-dependent restoration of the dehydrogenase activity on preincubation with DTT and GSH. Such thiol-treated preparations also gave increased ADP-ribosylation activity with DTT, and to a lesser extent with GSH. The SNP-modified protein was unable to catalyze this activity with the native yeast enzyme and native and heat-inactivated muscle enzyme. It was possible to generate the ADP-ribosylation activity in muscle GAPD, by an NO-independent mechanism, on dialysis in Tris buffer under aerobic conditions , and on incubating with NADPH, but not NADH, in muscle and brain, but not yeast, enzymes. The results suggest that the inverse relationship of the dehydrogenase and ADP-ribosylation activities is coincidental but not correlated

Competing peroxidase and oxidase reactions in scopoletin-dependent H2O2-initiated oxidation of NADH by horseradish peroxidase

Addition of NADH inhibited the peroxidative loss of scopoletin in presence of horseradish peroxidase and H2O2 and decreased the ratio of scopoletin (consumed):H2O2 (added). Concomitantly NADH was oxidized and oxygen was consumed with a stoichiometry of NADH: O-2 of 2:1. On step-wise addition of a small concentration of H2O2 a high rate of NADH oxidation was obtained for a progressively decreasing time period followed by termination of the reaction with NADH:H2O2 ratio decreasing from about 40 to 10. The rate of NADH oxidation increased linearly with increase in scopoletin concentration. Other phenolic compounds including p-coumarate also supported this reaction to a variable degree. A 418-nm absorbing compound;d accumulated during oxidation of NADH. The effectiveness of a small concentration of H2O2 in supporting NADH oxidation increased in presence of SOD and decreased in presence of cytochrome c, but the reaction terminated even in their presence. The results indicate that the peroxidase is not continuously generating H2O2 during scopolerin-mediated NADH oxidation and that both peroxidase and oxidase reactions occur simultaneously competing for an active form of the enzyme.

Interactions of L-serine at the active site of serine hydroxymethyltransferases: induction of thermal stability

Serine hydroxymethyltransferase (SHMT), EC 2.1.2.1, exhibits broad substrate and reaction specificity. In addition to cleaving many 3-hydroxyamino acids to glycine and an aldehyde, the enzyme also catalyzed the decarboxylation, transamination and racemization of several substrate analogues of amino acids. To elucidate the mechanism of interaction of substrates, especially L-serine with the enzyme, a comparative study of interaction of L-serine with the enzyme from sheep liver and Escherichia coli, was carried out. The heat stability of both the enzymes was enhanced in the presence of serine, although to different extents. Thermal denaturation monitored by spectral changes indicated an alteration in the apparent T, of sheep liver and E. coli SHMTs from 55 +/- 1 degrees C to 72 +/- 3 degrees C at 40 mM serine and from 67 +/- 1 degrees C to 72 +/- 1 degrees C at 20 mM serine, respectively. Using stopped flow spectrophotometry k values of (49 +/- 5)(.)10(-3) s(-1) and (69 +/- 7).10(-3) s(-1) for sheep liver and E. coli enzymes were determined at 50 mM serine. The binding of serine monitored by intrinsic fluorescence and sedimentation velocity measurements indicated that there was no generalized change in the structure of both proteins. However, visible CD measurements indicated a change in the asymmetric environment of pyridoxal 5'-phosphate at the active site upon binding of serine to both the enzymes. The formation of an external aldimine was accompanied by a change in the secondary structure of the enzymes monitored by far UV-CD spectra. Titration microcalorimetric studies in the presence of serine (8 mM) also demonstrated a single class of binding and the conformational changes accompanying the binding of serine to the enzyme resulted in a more compact structure leading to increased thermal stability of the enzyme.

Identification of the active-site peptide of 2,3-dihydroxybenzoic acid decarboxylase from Aspergillus oryzae

The non-oxidative decarboxylation of aromatic acids is a poorly understood reaction. The transformation of 2,3-dihydroxybenzoic acid to catechol in the fungal metabolism of indole is a prototype of such a reaction. 2,3-Dihydroxybenzoic acid decarboxylase (EC 4.1.1.46) which catalyzes this reaction was purified to homogeneity from anthranilate induced cultures of Aspergillus oryzae using affinity chromatography. The enzyme did not require cofactors like NAD(+), PLP, TPP or metal ions for its activity. There was no spectral evidence for the presence of enzyme bound cofactors. The preparation, which was adjudged homogeneous by the criteria of SDS-PAGE, sedimentation analysis and N-terminal analysis, was characterized for its physicochemical and kinetic parameters. The enzyme was inactivated by group-specific modifiers like diethyl pyrocarbonate (DEPC) and N-ethylmaleimide (NEM). The kinetics of inactivation by DEPC suggested the presence of a single class of essential histidine residues, the second order rate constant of inactivation for which was 12.5 M(-1) min(-1). A single class of cysteine residues was modified by NEM with a second order rate constant of 33 M(-1) min(-1). Substrate analogues protected the enzyme against inactivation by both DEPC and NEM, suggesting the Location of the essential histidine and cysteine to be at the active site of the enzyme. The incorporation of radiolabelled NEM in a differential labelling experiment was 0.73 mol per mol subunit confirming the presence of a single essential cysteine per active-site. Differentially labelled enzyme was enzymatically cleaved and the peptide bearing the label was purified and sequenced. The active-site peptide LLGLAETCK and the N-terminal sequence MLGKIALEEAFALPRFEEKT did not bear any similarity to sequences reported in the Swiss-Prot Protein Sequence Databank, a reflection probably of the unique primary structure of this novel enzyme. The sequences reported in this study will appear in the Swiss-Prot Protein Sequence Databank under the accession number P80402.

DEPTH: a web server to compute depth and predict small-molecule binding cavities in proteins

Depth measures the extent of atom/residue burial within a protein. It correlates with properties such as protein stability, hydrogen exchange rate, protein-protein interaction hot spots, post-translational modification sites and sequence variability. Our server, DEPTH, accurately computes depth and solvent-accessible surface area (SASA) values. We show that depth can be used to predict small molecule ligand binding cavities in proteins. Often, some of the residues lining a ligand binding cavity are both deep and solvent exposed. Using the depth-SASA pair values for a residue, its likelihood to form part of a small molecule binding cavity is estimated. The parameters of the method were calibrated over a training set of 900 high-resolution X-ray crystal structures of single-domain proteins bound to small molecules (molecular weight < 1.5 KDa). The prediction accuracy of DEPTH is comparable to that of other geometry-based prediction methods including LIGSITE, SURFNET and Pocket-Finder (all with Matthew's correlation coefficient of similar to 0.4) over a testing set of 225 single and multi-chain protein structures. Users have the option of tuning several parameters to detect cavities of different sizes, for example, geometrically flat binding sites. The input to the server is a protein 3D structure in PDB format. The users have the option of tuning the values of four parameters associated with the computation of residue depth and the prediction of binding cavities. The computed depths, SASA and binding cavity predictions are displayed in 2D plots and mapped onto 3D representations of the protein structure using Jmol. Links are provided to download the outputs. Our server is useful for all structural analysis based on residue depth and SASA, such as guiding site-directed mutagenesis experiments and small molecule docking exercises, in the context of protein functional annotation and drug discovery.

Allosteric Communication in Cysteinyl tRNA Synthetase A NETWORK OF DIRECT AND INDIRECT READOUT

Protein structure networks are constructed for the identification of long-range signaling pathways in cysteinyl tRNA synthetase (CysRS). Molecular dynamics simulation trajectory of CysRS-ligand complexes were used to determine conformational ensembles in order to gain insight into the allosteric signaling paths. Communication paths between the anticodon binding region and the aminoacylation region have been identified. Extensive interaction between the helix bundle domain and the anticodon binding domain, resulting in structural rigidity in the presence of tRNA, has been detected. Based on the predicted model, six residues along the communication paths have been examined by mutations (single and double) and shown to mediate a coordinated coupling between anticodon recognition and activation of amino acid at the active site. This study on CysRS clearly shows that specific key residues, which are involved in communication between distal sites in allosteric proteins but may be elusive in direct structure analysis, can be identified from dynamics of protein structure networks.

Asparagine and glutamine differ in their propensities to form specific side chain-backbone hydrogen bonded motifs in proteins

Short range side chain-backbone hydrogen bonded motifs involving Asn and Gln residues have been identified from a data set of 1370 protein crystal structures (resolution = 1.5 angstrom). Hydrogen bonds involving residues i - 5 to i + 5 have been considered. Out of 12,901 Asn residues, 3403 residues (26.4%) participate in such interactions, while out of 10,934 Gln residues, 1780 Gln residues (16.3%) are involved in these motifs. Hydrogen bonded ring sizes (Cn, where n is the number of atoms involved), directionality and internal torsion angles are used to classify motifs. The occurrence of the various motifs in the contexts of protein structure is illustrated. Distinct differences are established between the nature of motifs formed by Asn and Gln residues. For Asn, the most highly populated motifs are the C10 (COdi .NHi + 2), C13 (COdi .NHi + 3) and C17 (NdHi .COi - 4) structures. In contrast, Gln predominantly forms C16 (COei .NHi - 3), C12 (NeHi .COi - 2), C15 (NeHi .COi - 3) and C18 (NeHi .COi - 4) motifs, with only the C18motif being analogous to the Asn C17structure. Specific conformational types are established for the Asn containing motifs, which mimic backbone beta-turns and a-turns. Histidine residues are shown to serve as a mimic for Asn residues in side chain-backbone hydrogen bonded ring motifs. Illustrative examples from protein structures are considered. Proteins 2012; (c) 2011 Wiley Periodicals, Inc.

Dimethyl sulfoxide induced structural transformations and non-monotonic concentration dependence of conformational fluctuation around active site of lysozyme

Experimental studies have observed significant changes in both structure and function of lysozyme (and other proteins) on addition of a small amount of dimethyl sulfoxide (DMSO) in aqueous solution. Our atomistic molecular dynamic simulations of lysozyme in water-DMSO reveal the following sequence of changes on increasing DMSO concentration. (i) At the initial stage (around 5% DMSO concentration) protein's conformational flexibility gets markedly suppressed. From study of radial distribution functions, we attribute this to the preferential solvation of exposed protein hydrophobic residues by the methyl groups of DMSO. (ii) In the next stage (10-15% DMSO concentration range), lysozome partially unfolds accompanied by an increase both in fluctuation and in exposed protein surface area. (iii) Between 15-20% concentration ranges, both conformational fluctuation and solvent accessible protein surface area suddenly decrease again indicating the formation of an intermediate collapse state. These results are in good agreement with near-UV circular dichroism (CD) and fluorescence studies. We explain this apparently surprising behavior in terms of a structural transformation which involves clustering among the methyl groups of DMSO. (iv) Beyond 20% concentration of DMSO, the protein starts its final sojourn towards the unfolding state with further increase in conformational fluctuation and loss in native contacts. Most importantly, analysis of contact map and fluctuation near the active site reveal that both partial unfolding and conformational fluctuations are centered mostly on the hydrophobic core of active site of lysozyme. Our results could offer a general explanation and universal picture of the anomalous behavior of protein structure-function observed in the presence of cosolvents (DMSO, ethanol, tertiary butyl alcohol, dioxane) at their low concentrations. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3694268]

Local Structural Differences in Homologous Proteins: Specificities in Different SCOP Classes

The constant increase in the number of solved protein structures is of great help in understanding the basic principles behind protein folding and evolution. 3-D structural knowledge is valuable in designing and developing methods for comparison, modelling and prediction of protein structures. These approaches for structure analysis can be directly implicated in studying protein function and for drug design. The backbone of a protein structure favours certain local conformations which include alpha-helices, beta-strands and turns. Libraries of limited number of local conformations (Structural Alphabets) were developed in the past to obtain a useful categorization of backbone conformation. Protein Block (PB) is one such Structural Alphabet that gave a reasonable structure approximation of 0.42 angstrom. In this study, we use PB description of local structures to analyse conformations that are preferred sites for structural variations and insertions, among group of related folds. This knowledge can be utilized in improving tools for structure comparison that work by analysing local structure similarities. Conformational differences between homologous proteins are known to occur often in the regions comprising turns and loops. Interestingly, these differences are found to have specific preferences depending upon the structural classes of proteins. Such class-specific preferences are mainly seen in the all-beta class with changes involving short helical conformations and hairpin turns. A test carried out on a benchmark dataset also indicates that the use of knowledge on the class specific variations can improve the performance of a PB based structure comparison approach. The preference for the indel sites also seem to be confined to a few backbone conformations involving beta-turns and helix C-caps. These are mainly associated with short loops joining the regular secondary structures that mediate a reversal in the chain direction. Rare beta-turns of type I' and II' are also identified as preferred sites for insertions.

PocketAnnotate: towards site-based function annotation

A computational pipeline PocketAnnotate for functional annotation of proteins at the level of binding sites has been proposed in this study. The pipeline integrates three in-house algorithms for site-based function annotation: PocketDepth, for prediction of binding sites in protein structures; PocketMatch, for rapid comparison of binding sites and PocketAlign, to obtain detailed alignment between pair of binding sites. A novel scheme has been developed to rapidly generate a database of non-redundant binding sites. For a given input protein structure, putative ligand-binding sites are identified, matched in real time against the database and the query substructure aligned with the promising hits, to obtain a set of possible ligands that the given protein could bind to. The input can be either whole protein structures or merely the substructures corresponding to possible binding sites. Structure-based function annotation at the level of binding sites thus achieved could prove very useful for cases where no obvious functional inference can be obtained based purely on sequence or fold-level analyses. An attempt has also been made to analyse proteins of no known function from Protein Data Bank. PocketAnnotate would be a valuable tool for the scientific community and contribute towards structure-based functional inference. The web server can be freely accessed at http://proline.biochem.iisc.ernet.in/pocketannotate/.