Arabidopsis thaliana J-class heat shock proteins: cellular stress sensors

Plants are sessile organisms that have evolved a variety of mechanisms to maintain their cellular homeostasis under stressful environmental conditions. Survival of plants under abiotic stress conditions requires specialized group of heat shock protein machinery, belonging to Hsp70:J-protein family. These heat shock proteins are most ubiquitous types of chaperone machineries involved in diverse cellular processes including protein folding, translocation across cell membranes, and protein degradation. They play a crucial role in maintaining the protein homeostasis by reestablishing functional native conformations under environmental stress conditions, thus providing protection to the cell. J-proteins are co-chaperones of Hsp70 machine, which play a critical role by stimulating Hsp70s ATPase activity, thereby stabilizing its interaction with client proteins. Using genome-wide analysis of Arabidopsis thaliana, here we have outlined identification and systematic classification of J-protein co-chaperones which are key regulators of Hsp70s function. In comparison with Saccharomyces cerevisiae model system, a comprehensive domain structural organization, cellular localization, and functional diversity of A. thaliana J-proteins have also been summarized. Electronic supplementary material The online version of this article (doi:10.1007/s10142-009-0132-0) contains supplementary material, which is available to authorized users.

Hyperexpression of chicken riboflavin carrier protein: Antibodies to the recombinant protein curtail pregnancy in rodents

Chicken riboflavin carrier protein (RCP) is a phosphoglycoprotein present in the egg white and yolk of egg-laying animals and in the sera of laying hens and of estrogenized chicks. The RCP cDNA, encoding a protein of predictedMr27,000, has been cloned into a T7 polymerase-driven vector, and high-level expression was observed on induction with IPTG inEscherichia coli.The protein was largely localized in inclusion bodies when expressed at 37°C but was present in the cytosolic fraction when induced at 22°C. At 37°C, two major bands were detected in whole-cell lysates of the strain expressing the protein. N-terminal sequence analysis indicated that the two proteins represented translated products with and without the pelB leader sequence encoded in the pET20b vector, but both included an additional 10 amino acids generated during cloning procedures. The inclusion body obtained at 37°C, on extraction with detergent, led to preferential solubilization of the protein without the pelB signal sequence. The solubilized recombinant RCP was recognized by polyclonal antisera to native RCP but radioimmunoassay revealed quantitative differences in the epitopes exhibited by the recombinant protein. Thus, sequence-specific monoclonal antibodies to chicken RCP also cross-reacted with the recombinant protein with almost equal efficiency, but antibodies which recognize conformation-dependent epitopes showed relatively reduced cross-reactivity with the recombinant protein. Polyclonal antibodies to recombinant RCP were able to recognize both the native and the denatured RCP. Administration of recombinant RCP antisera to pregnant mice led to embryonic resorption leading to early pregnancy termination. These findings reveal that the recombinant protein will be useful for investigations related to the mechanism of pregnancy termination on immunoneutralization of RCP in mammals, as well as in unraveling folding properties of RCP in terms of its ligand binding and antigenetic determinants exposed at its surface.

Identification and thermodynamic characterization of molten globule states of periplasmic binding proteins

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of Delta C-p of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.

Amino acid interaction preferences in proteins

Understanding the key factors that influence the interaction preferences of amino acids in the folding of proteins have remained a challenge. Here we present a knowledge-based approach for determining the effective interactions between amino acids based on amino acid type, their secondary structure, and the contact based environment that they find themselves in the native state structure as measured by their number of neighbors. We find that the optimal information is approximately encoded in a 60 x 60 matrix describing the 20 types of amino acids in three distinct secondary structures (helix, beta strand, and loop). We carry out a clustering scheme to understand the similarity between these interactions and to elucidate a nonredundant set. We demonstrate that the inferred energy parameters can be used for assessing the fit of a given sequence into a putative native state structure.

High-level expression of biologically active glycoprotein hormones in Pichia pastoris strains—selection of strain GS115, and not X-33, for the production of biologically active N-glycosylated 15N-labeled phCG

The methylotrophic yeast Pichia pastoris is widely used for the production of recombinant glycoproteins. With the aim to generate biologically active 15N-labeled glycohormones for conformational studies focused on the unravelling of the NMR structures in solution, the P. pastoris strains GS115 and X-33 were explored for the expression of human chorionic gonadotropin (phCG) and human follicle-stimulating hormone (phFSH). In agreement with recent investigations on the N-glycosylation of phCG, produced in P. pastoris GS115, using ammonia/glycerol-methanol as nitrogen/carbon sources, the N-glycosylation pattern of phCG, synthesized using NH4Cl/glucose–glycerol–methanol, comprised neutral and charged, phosphorylated high-mannose-type N-glycans (Man8–15GlcNAc2). However, the changed culturing protocol led to much higher amounts of glycoprotein material, which is of importance for an economical realistic approach of the aimed NMR research. In the context of these studies, attention was also paid to the site specific N-glycosylation in phCG produced in P. pastoris GS115. In contrast to the rather simple N-glycosylation pattern of phCG expressed in the GS115 strain, phCG and phFSH expressed in the X-33 strain revealed, besides neutral high-mannose-type N-glycans, also high concentrations of neutral hypermannose-type N-glycans (Manup-to-30GlcNAc2). The latter finding made the X-33 strain not very suitable for generating 15N-labeled material. Therefore, 15N-phCG was expressed in the GS115 strain using the new optimized protocol. The 15N-enrichment was evaluated by 15N-HSQC NMR spectroscopy and GLC-EI/MS. Circular dichroism studies indicated that 15N-phCG/GS115 had the same folding as urinary hCG. Furthermore, 15N-phCG/GS115 was found to be similar to the unlabeled protein in every respect as judged by radioimmunoassay, radioreceptor assays, and in vitro bioassays.

Stereochemistry of peptides containing 1-aminocycloheptane-1-carboxylic acid (Ac7c)

The crystal structures of four peptides incorporating 1-aminocycloheptane-1-carboxylic acid (Ac7c) are described. Boc-Aib-Ac7c-NHMe and Boc-Pro-Ac7c-Ala-OMe adopt beta-turn conformations stabilized by an intramolecular 4----1 hydrogen bond, the former folding into a type-I/III beta-turn and the latter into a type-II beta-turn. In the dipeptide esters, Boc-Aib-Ac7c-OMe and Boc-Pro-Ac7c-OMe, the Ac7c and Aib residues adopt helical conformations, while the Pro residue remains semi-extended in both the molecules of Boc-Pro-Ac7c-OMe found in the asymmetric unit. The cycloheptane ring of Ac7c residues adopts a twist-chair conformation in all the peptides studied. 1H-NMR studies in CDCl3 and (CD3)2SO and IR studies in CDCl3 suggest that Boc-Aib-Ac7c-NHMe and Boc-Pro-Ac7c-Ala-OMe maintain the beta-turn conformations in solution.

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.

Stereochemistry of peptides containing 1-aminocycloheptane-1-carboxylic acid (Ac-c)

The crystal structures of four peptides incorporating 1-aminocycloheptane-1-carboxylic acid (Ac7c) are described. Boc-Aib-Ac7c-NHMe and Boc-Pro-Ac7c-Ala-OMe adopt beta-turn conformations stabilized by an intramolecular 4----1 hydrogen bond, the former folding into a type-I/III beta-turn and the latter into a type-II beta-turn. In the dipeptide esters, Boc-Aib-Ac7c-OMe and Boc-Pro-Ac7c-OMe, the Ac7c and Aib residues adopt helical conformations, while the Pro residue remains semi-extended in both the molecules of Boc-Pro-Ac7c-OMe found in the asymmetric unit. The cycloheptane ring of Ac7c residues adopts a twist-chair conformation in all the peptides studied. 1H-NMR studies in CDCl3 and (CD3)2SO and IR studies in CDCl3 suggest that Boc-Aib-Ac7c-NHMe and Boc-Pro-Ac7c-Ala-OMe maintain the beta-turn conformations in solution.

Non-Protein Amino Acids in the Design of Secondary Structure Scaffolds

The use of stereochemically constrained amino acids permits the design of short peptides as models for protein secondary structures. Amino acid residues that are restrained to a limited range of backbone torsion angles (ϕ-ψ) may be used as folding nuclei in the design of helices and β-hairpins. α-Amino-isobutyric acid (Aib) and related Cαα dialkylated residues are strong promoters of helix formation, as exemplified by a large body of experimentally determined structures of helical peptides. DPro-Xxx sequences strongly favor type II’ turn conformations, which serve to nucleate registered β-hairpin formation. Appropriately positioned DPro-Xxx segments may be used to nucleate the formation of multistranded antiparallel β-sheet structures. Mixed (α/β) secondary structures can be generated by linking rigid modules of helices and β-hairpins. The approach of using stereochemically constrained residues promotes folding by limiting the local structural space at specific residues. Several aspects of secondary structure design are outlined in this chapter, along with commonly used methods of spectroscopic characterization.

Helical Conformations of Hexapeptides Containing N-Terminus Diproline Segments

The role of N-terminus diproline segments in facilitating helical folding in short peptides has been investigated in a set of model hexapeptides of the type Piv-Xxx-Yyy-Aib-Leu-Aib-Phe-OMe (Piv, pivaloyl). Nine sequences have been investigated with the following N-terminus dipeptide segments: (D)Pro-Ala (4) and Pro-Psi Pro (5, Psi, pseudoproline), Ala-Ala (6), Ala-Pro (7), Pro-Ala (8), Aib-Ala (9), Ala-Aib (10). The analog sequences Piv-Pro-Pro-Ala-Leu-Aib-Phe-OMe (2) and Piv-Pro-Pro-Ala-Aib-Ala-Aib-OMe (3) have also been studied. Solid state conformations have been determined by X-ray crystallography for peptides 4, 6, and 8 and compared with the previously determined crystal structure of peptide 1 (Boc-Pro-Pro-Aib-Leu-Aib-Val-OMe); (Rai et al., JACS 2006, 128, 7916-7928). Peptides 1 and 6 adopt almost identical helical conformations with unfolding of the helix at the N-terminus Pro (1) residue. Peptide 4 reveals the anticipated (D)Pro-Ala type II' beta-turn, followed by a stretch of 3(10)-helix. Peptide 8 adopts a folded conformation stabilized by four successive 4 -> 1 intramolecular hydrogen bonds. Ala (2) adopts an alpha(L) conformation, resulting in a type II beta-turn conformation followed by a stretch of 3(10)-helix. Conformational properties in solution were probed using solvent perturbation of NH chemical shifts which permit delineation of hydrogen bonded NH groups and nuclear Overhauser effects (NOEs) between backbone protons, which are diagnostic of local residue conformations. The results suggest that continuous helical conformations are indeed significantly populated for peptides 2 and 3. Comparison of the results for peptides 1 and 2, suggest that there is a significant influence of the residue that follows diproline segments in influencing backbone folding. (C) 2010 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 94: 360-370, 2010.

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/

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Thermodynamic and structural studies of cavity formation in proteins suggest that loss of packing interactions rather than the hydrophobic effect dominates the observed energetics

The hydrophobic effect is widely believed to be an important determinant of protein stability. However, it is difficult to obtain unambiguous experimental estimates of the contribution of the hydrophobic driving force to the overall free energy of folding. Thermodynamic and structural studies of large to small substitutions in proteins are the most direct method of measuring this contribution. We have substituted the buried residue Phe8 in RNase S with alanine, methionine, and norleucine, Binding thermodynamics and structures were characterized by titration calorimetry and crystallography, respectively. The crystal structures of the RNase S F8A, F8M, and F8Nle mutants indicate that the protein tolerates the changes without any main chain adjustments, The correlation of structural and thermodynamic parameters associated with large to small substitutions was analyzed for nine mutants of RNase S as well as 32 additional cavity-containing mutants of T4 lysozyme, human lysozyme, and barnase. Such substitutions were typically found to result in negligible changes in Delta C-p and positive values of both Delta Delta H degrees and aas of folding. Enthalpic effects were dominant, and the sign of Delta Delta S is the opposite of that expected from the hydrophobic effect. Values of Delta Delta G degrees and Delta Delta H degrees correlated better with changes in packing parameters such as residue depth or occluded surface than with the change in accessible surface area upon folding. These results suggest that the loss of packing interactions rather than the hydrophobic effect is a dominant contributor to the observed energetics for large to small substitutions. Hence, estimates of the magnitude of the hydrophobic driving force derived from earlier mutational studies are likely to be significantly in excess of the actual value.

Interplay between multiple length and time scales in complex chemical systems

Processes in complex chemical systems, such as macromolecules, electrolytes, interfaces, micelles and enzymes, can span several orders of magnitude in length and time scales. The length and time scales of processes occurring over this broad time and space window are frequently coupled to give rise to the control necessary to ensure specificity and the uniqueness of the chemical phenomena. A combination of experimental, theoretical and computational techniques that can address a multiplicity of length and time scales is required in order to understand and predict structure and dynamics in such complex systems. This review highlights recent experimental developments that allow one to probe structure and dynamics at increasingly smaller length and time scales. The key theoretical approaches and computational strategies for integrating information across time-scales are discussed. The application of these ideas to understand phenomena in various areas, ranging from materials science to biology, is illustrated in the context of current developments in the areas of liquids and solvation, protein folding and aggregation and phase transitions, nucleation and self-assembly.

Extended colored Petri net: An efficient tool for analyzing concurrent systems

This paper proposes a novel and simple definition of general colored Petri nets. This definition is coherent with that of (uncolored) Petri nets, preserves the reflexivity of the original net and is extended to represent inhibitors. Also suggested are systematic and formal merging rules to obtain a well-formed structure of the extended colored Petri net by folding a given uncolored net. Finally, we present a technique to compute colored invariants by selecting colored RP-subnets. On the average, the proposed technique performs better than the existing ones. The analysis procedure is explained through an illustrative example of a three-level interrupt-priority-handler scheme.