Structural studies of protein unfolding

Methods for macromolecular structure determination (NMR and crystallography) are now being used to get structural information on partially folded and unfolded states of proteins. These techniques, in combination with proton hydrogen exchange studies are powerful tools to extract information on non-native states of proteins. This review discusses progress In this area of protein folding.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

Network properties of protein-decoy structures

Convergence of the vast sequence space of proteins into a highly restricted fold/conformational space suggests a simple yet unique underlying mechanism of protein folding that has been the subject of much debate in the last several decades. One of the major challenges related to the understanding of protein folding or in silico protein structure prediction is the discrimination of non-native structures/decoys from the native structure. Applications of knowledge-based potentials to attain this goal have been extensively reported in the literature. Also, scoring functions based on accessible surface area and amino acid neighbourhood considerations were used in discriminating the decoys from native structures. In this article, we have explored the potential of protein structure network (PSN) parameters to validate the native proteins against a large number of decoy structures generated by diverse methods. We are guided by two principles: (a) the PSNs capture the local properties from a global perspective and (b) inclusion of non-covalent interactions, at all-atom level, including the side-chain atoms, in the network construction accommodates the sequence dependent features. Several network parameters such as the size of the largest cluster, community size, clustering coefficient are evaluated and scored on the basis of the rank of the native structures and the Z-scores. The network analysis of decoy structures highlights the importance of the global properties contributing to the uniqueness of native structures. The analysis also exhibits that the network parameters can be used as metrics to identify the native structures and filter out non-native structures/decoys in a large number of data-sets; thus also has a potential to be used in the protein `structure prediction' problem.

Identification of an exported heat shock protein 70 in Plasmodium falciparum

Host cell remodelling is a hallmark of malaria pathogenesis. It involves protein folding, unfolding and trafficking events and thus participation of chaperones such as Hsp70s and Hsp40s is well speculated. Until recently, only Hsp40s were thought to be the sole representative of the parasite chaperones in the exportome. However, based on the re-annotated Plasmodium falciparum genome sequence, a putative candidate for exported Hsp70 has been reported, which otherwise was known to be a pseudogene. We raised a specific antiserum against a C-terminal peptide uniquely present in PfHsp70-x. Immunoblotting and immunofluorescence-based approaches in combination with sub-cellular fractionation by saponin and streptolysin-O have been taken to determine the expression and localization of PfHsp70-x in infected erythrocyte. The re-annotated sequence of PfHsp70-x reveals it to be a functional protein with an endoplasmic reticulum signal peptide. It gets maximally expressed at the schizont stage of intra-erythrocytic life cycle. Majority of the protein localizes to the parasitophorous vacuole and some of it gets exported to the erythrocyte compartment where it associates with Maurer's clefts. The identification of an exported parasite Hsp70 chaperone presents us with the fact that the parasite has evolved customized chaperones which might be playing crucial roles in aspects of trafficking and host cell remodelling.

Protein Structure and Function: Looking through the Network of Side-Chain Interactions

Network theory has become an excellent method of choice through which biological data are smoothly integrated to gain insights into complex biological problems. Understanding protein structure, folding, and function has been an important problem, which is being extensively investigated by the network approach. Since the sequence uniquely determines the structure, this review focuses on the networks of non-covalently connected amino acid side chains in proteins. Questions in structural biology are addressed within the framework of such a formalism. While general applications are mentioned in this review, challenging problems which have demanded the attention of scientific community for a long time, such as allostery and protein folding, are considered in greater detail. Our aim has been to explore these important problems through the eyes of networks. Various methods of constructing protein structure networks (PSN) are consolidated. They include the methods based on geometry, edges weighted by different schemes, and also bipartite network of protein-nucleic acid complexes. A number of network metrics that elegantly capture the general features as well as specific features related to phenomena, such as allostery and protein model validation, are described. Additionally, an integration of network theory with ensembles of equilibrium structures of a single protein or that of a large number of structures from the data bank has been presented to perceive complex phenomena from network perspective. Finally, we discuss briefly the capabilities, limitations, and the scope for further explorations of protein structure networks.

Cloning and expression of glucose regulated protein 78 (GRP78) in Fenneropenaeus chinensis

GRP78 (78 kDa glucose-regulated protein), also known as BiP (immunoglobulin heavy-chain-binding protein), is an essential regulator of endoplasmic reticulum (ER) homeostasis because of its multiple functions in protein folding, ER calcium binding, and controlling of the activation of transmembrane ER stress sensors. In this report, we cloned the full length cDNA of GRP78 (FcGRP78) from Chinese shrimp Fenneropenaeus chinensis. This cDNA revealed a 2,325 bp with 1,968 bp open reading frame encoding 655 amino acids. This is the first reported GRP78 gene in Crustacea. The deduced amino acid sequence of FcGRP78 shared high identity with previously reported insect GRP78s: 86, 87 and 85% identity with GRP78s of Drosophila melanogaster, Aedes aegypti and Bombyx mori, respectively. Northern blot analysis shows that FcGRP78 is ubiquitously expressed in tissues of shrimp. Heat shock at 35A degrees C significantly enhanced the expression of FcGRP78 at the first hour, reached the maximum at 4 h post heat shock, dropped after that and resumed to the normal level until 48 h of post recovery at 25A degrees C. Additionally, differential expression of FcGRP78 was detected in haemocytes, hepatopancreas and lymphoid organ when shrimp were challenged by white spot syndrome virus (WSSV). We inferred that FcGRP78 may play important roles in chaperoning, protein folding and immune function of shrimp.

Protein Docking by the Underestimation of Free Energy Funnels in the Space of Encounter Complexes

Similarly to protein folding, the association of two proteins is driven by a free energy funnel, determined by favorable interactions in some neighborhood of the native state. We describe a docking method based on stochastic global minimization of funnel-shaped energy functions in the space of rigid body motions (SE(3)) while accounting for flexibility of the interface side chains. The method, called semi-definite programming-based underestimation (SDU), employs a general quadratic function to underestimate a set of local energy minima and uses the resulting underestimator to bias further sampling. While SDU effectively minimizes functions with funnel-shaped basins, its application to docking in the rotational and translational space SE(3) is not straightforward due to the geometry of that space. We introduce a strategy that uses separate independent variables for side-chain optimization, center-to-center distance of the two proteins, and five angular descriptors of the relative orientations of the molecules. The removal of the center-to-center distance turns out to vastly improve the efficiency of the search, because the five-dimensional space now exhibits a well-behaved energy surface suitable for underestimation. This algorithm explores the free energy surface spanned by encounter complexes that correspond to local free energy minima and shows similarity to the model of macromolecular association that proceeds through a series of collisions. Results for standard protein docking benchmarks establish that in this space the free energy landscape is a funnel in a reasonably broad neighborhood of the native state and that the SDU strategy can generate docking predictions with less than 5 � ligand interface Ca root-mean-square deviation while achieving an approximately 20-fold efficiency gain compared to Monte Carlo methods.

Unfolded protein response genes regulated by CED-1 are required for Caenorhabditis elegans innate immunity.

The endoplasmic reticulum stress response, also known as the unfolded protein response (UPR), has been implicated in the normal physiology of immune defense and in several disorders, including diabetes, cancer, and neurodegenerative disease. Here, we show that the apoptotic receptor CED-1 and a network of PQN/ABU proteins involved in a noncanonical UPR response are required for proper defense to pathogen infection in Caenorhabditis elegans. A full-genome microarray analysis indicates that CED-1 functions to activate the expression of pqn/abu genes. We also show that ced-1 and pqn/abu genes are required for the survival of C. elegans exposed to live Salmonella enterica, and that overexpression of pqn/abu genes confers protection against pathogen-mediated killing. The results indicate that unfolded protein response genes, regulated in a CED-1-dependent manner, are involved in the C. elegans immune response to live bacteria.

Short-lived alpha-helical intermediates in the folding of beta-sheet proteins.

Several lines of evidence point strongly toward the importance of highly alpha-helical intermediates in the folding of all globular proteins, regardless of their native structure. However, experimental refolding studies demonstrate no observable alpha-helical intermediate during refolding of some beta-sheet proteins and have dampened enthusiasm for this model of protein folding. In this study, beta-sheet proteins were hypothesized to have potential to form amphiphilic helices at a period of <3.6 residues/turn that matches or exceeds the potential at 3.6 residues/turn. Hypothetically, such potential is the basis for an effective and unidirectional mechanism by which highly alpha-helical intermediates might be rapidly disassembled during folding and potentially accounts for the difficulty in detecting highly alpha-helical intermediates during the folding of some proteins. The presence of this potential was confirmed, indicating that a model entailing ubiquitous formation of alpha-helical intermediates during the folding of globular proteins predicts previously unrecognized features of primary structure. Further, the folding of fatty acid binding protein, a predominantly beta-sheet protein that exhibits no apparent highly alpha-helical intermediate during folding, was dramatically accelerated by 2,2,2-trifluoroethanol, a solvent that stabilizes alpha-helical structure. This observation suggests that formation of an alpha-helix can be a rate-limiting step during folding of a predominantly beta-sheet protein and further supports the role of highly alpha-helical intermediates in the folding of all globular proteins.

Thermodynamic analysis of a molecular chaperone binding to unfolded protein substrates.

Molecular chaperones are a highly diverse group of proteins that recognize and bind unfolded proteins to facilitate protein folding and prevent nonspecific protein aggregation. The mechanisms by which chaperones bind their protein substrates have been studied for decades. However, there are few reports about the affinity of molecular chaperones for their unfolded protein substrates. Thus, little is known about the relative binding affinities of different chaperones and about the relative binding affinities of chaperones for different unfolded protein substrates. Here we describe the application of SUPREX (stability of unpurified proteins from rates of H-D exchange), an H-D exchange and MALDI-based technique, in studying the binding interaction between the molecular chaperone Hsp33 and four different unfolded protein substrates, including citrate synthase, lactate dehydrogenase, malate dehydrogenase, and aldolase. The results of our studies suggest that the cooperativity of the Hsp33 folding-unfolding reaction increases upon binding with denatured protein substrates. This is consistent with the burial of significant hydrophobic surface area in Hsp33 when it interacts with its substrate proteins. The SUPREX-derived K(d) values for Hsp33 complexes with four different substrates were all found to be within the range of 3-300 nM.

De novo design and molecular assembly of a transmembrane diporphyrin-binding protein complex.

The de novo design of membrane proteins remains difficult despite recent advances in understanding the factors that drive membrane protein folding and association. We have designed a membrane protein PRIME (PoRphyrins In MEmbrane) that positions two non-natural iron diphenylporphyrins (Fe(III)DPP's) sufficiently close to provide a multicentered pathway for transmembrane electron transfer. Computational methods previously used for the design of multiporphyrin water-soluble helical proteins were extended to this membrane target. Four helices were arranged in a D(2)-symmetrical bundle to bind two Fe(II/III) diphenylporphyrins in a bis-His geometry further stabilized by second-shell hydrogen bonds. UV-vis absorbance, CD spectroscopy, analytical ultracentrifugation, redox potentiometry, and EPR demonstrate that PRIME binds the cofactor with high affinity and specificity in the expected geometry.

Conformational kinetics reveals affinities of protein conformational states.

Most biological reactions rely on interplay between binding and changes in both macromolecular structure and dynamics. Practical understanding of this interplay requires detection of critical intermediates and determination of their binding and conformational characteristics. However, many of these species are only transiently present and they have often been overlooked in mechanistic studies of reactions that couple binding to conformational change. We monitored the kinetics of ligand-induced conformational changes in a small protein using six different ligands. We analyzed the kinetic data to simultaneously determine both binding affinities for the conformational states and the rate constants of conformational change. The approach we used is sufficiently robust to determine the affinities of three conformational states and detect even modest differences in the protein's affinities for relatively similar ligands. Ligand binding favors higher-affinity conformational states by increasing forward conformational rate constants and/or decreasing reverse conformational rate constants. The amounts by which forward rate constants increase and reverse rate constants decrease are proportional to the ratio of affinities of the conformational states. We also show that both the affinity ratio and another parameter, which quantifies the changes in conformational rate constants upon ligand binding, are strong determinants of the mechanism (conformational selection and/or induced fit) of molecular recognition. Our results highlight the utility of analyzing the kinetics of conformational changes to determine affinities that cannot be determined from equilibrium experiments. Most importantly, they demonstrate an inextricable link between conformational dynamics and the binding affinities of conformational states.

Functional analysis of the C-terminal domain of the WbaP protein that mediates initiation of O antigen synthesis in Salmonella enterica

WbaP catalyzes the transfer of galactose-1-phosphate onto undecaprenyl phosphate (Und-P). The enzyme belongs to a large family of bacterial membrane proteins required for initiation of the synthesis of O antigen lipopolysaccharide and polysaccharide capsules. Previous work in our laboratory demonstrated that the last transmembrane helix and C-terminal tail region of WbaP (WbaP(CT)) are sufficient for enzymatic activity. Here, we demonstrate the cytoplasmic location of the WbaP C-terminal tail and show that WbaPCT domain N-terminally fused to thioredoxin (TrxA-WbaP(CT)) exhibits improved protein folding and enhanced transferase activity. Alanine replacement of highly conserved charged or polar amino acids identified seven critical residues for enzyme activity in vivo and in vitro. Four of these residues are located in regions predicted to be a-helical. These regions and their secondary structure predictions are conserved in distinct WbaP family members, suggesting they may contribute to form a conserved catalytic center.

First report of the use of a saxitoxin-protein conjugate to develop a DNA aptamer to a small molecule toxin

Saxitoxin (STX) is a low molecular weight neurotoxin mainly produced by certain marine dinoflagellates that, along with its family of similarly related paralytic shellfish toxins, may cause the potentially fatal intoxication known as paralytic shellfish poisoning. Illness and fatality rates are low due to the effective monitoring programs that determine when toxins exceed the established regulatory action level and effectuate shellfish harvesting closures accordingly. Such monitoring programs rely on the ability to rapidly screen large volumes of samples. Many of the screening assays currently available employ antibodies or live animals. This research focused on developing an analytical recognition element that would eliminate the challenges associated with the limited availability of antibodies and the use of animals. Here we report the discovery of a DNA aptamer that targets STX. Concentration-dependent and selective binding of the aptamer to STX was determined using a surface plasmon resonance sensor. Not only does this work represent the first reported aptamer to STX, but also the first aptamer to any marine biotoxin. A novel strategy of using a toxin-protein conjugate for DNA aptamer selection was successfully implemented to overcome the challenges associated with aptamer selection to small molecules. Taking advantage of such an approach could lead to increased diversity and accessibility of aptamers to low molecular weight toxins, which could then be incorporated as analytical recognition elements in diagnostic assays for foodborne toxin detection. The selected STX aptamer sequence is provided here, making it available to any investigator for use in assay development for the detection of STX.