Electrostatic Contributions to the Thermodynamics of Ribonuclease P Protein Folding

Electrostatic interactions are of fundamental importance in determining the structure and stability of macromolecules. For example, charge-charge interactions modulate the folding and binding of proteins and influence protein solubility. Electrostatic interactions are highly variable and can be both favorable and unfavorable. The ability to quantify these interactions is challenging but vital to understanding the detailed balance and major roles that they have in different proteins and biological processes. Measuring pKa values of ionizable groups provides a sensitive method for experimentally probing the electrostatic properties of a protein.

pKa values report the free energy of site-specific proton binding and provide a direct means of studying protein folding and pH-dependent stability. Using a combination of NMR, circular dichroism, and fluorescence spectroscopy along with singular value decomposition, we investigated the contributions of electrostatic interactions to the thermodynamic stability and folding of the protein subunit of Bacillus subtilis ribonuclease P, P protein. Taken together, the results suggest that unfavorable electrostatics alone do not account for the fact that P protein is intrinsically unfolded in the absence of ligand because the pKa differences observed between the folded and unfolded state are small. Presumably, multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

Analysis of Molecular Dynamics Simulations of Protein Folding

Microsecond long Molecular Dynamics (MD) trajectories of biomolecular processes are now possible due to advances in computer technology. Soon, trajectories long enough to probe dynamics over many milliseconds will become available. Since these timescales match the physiological timescales over which many small proteins fold, all atom MD simulations of protein folding are now becoming popular. To distill features of such large folding trajectories, we must develop methods that can both compress trajectory data to enable visualization, and that can yield themselves to further analysis, such as the finding of collective coordinates and reduction of the dynamics. Conventionally, clustering has been the most popular MD trajectory analysis technique, followed by principal component analysis (PCA). Simple clustering used in MD trajectory analysis suffers from various serious drawbacks, namely, (i) it is not data driven, (ii) it is unstable to noise and change in cutoff parameters, and (iii) since it does not take into account interrelationships amongst data points, the separation of data into clusters can often be artificial. Usually, partitions generated by clustering techniques are validated visually, but such validation is not possible for MD trajectories of protein folding, as the underlying structural transitions are not well understood. Rigorous cluster validation techniques may be adapted, but it is more crucial to reduce the dimensions in which MD trajectories reside, while still preserving their salient features. PCA has often been used for dimension reduction and while it is computationally inexpensive, being a linear method, it does not achieve good data compression. In this thesis, I propose a different method, a nonmetric multidimensional scaling (nMDS) technique, which achieves superior data compression by virtue of being nonlinear, and also provides a clear insight into the structural processes underlying MD trajectories. I illustrate the capabilities of nMDS by analyzing three complete villin headpiece folding and six norleucine mutant (NLE) folding trajectories simulated by Freddolino and Schulten [1]. Using these trajectories, I make comparisons between nMDS, PCA and clustering to demonstrate the superiority of nMDS. The three villin headpiece trajectories showed great structural heterogeneity. Apart from a few trivial features like early formation of secondary structure, no commonalities between trajectories were found. There were no units of residues or atoms found moving in concert across the trajectories. A flipping transition, corresponding to the flipping of helix 1 relative to the plane formed by helices 2 and 3 was observed towards the end of the folding process in all trajectories, when nearly all native contacts had been formed. However, the transition occurred through a different series of steps in all trajectories, indicating that it may not be a common transition in villin folding. The trajectories showed competition between local structure formation/hydrophobic collapse and global structure formation in all trajectories. Our analysis on the NLE trajectories confirms the notion that a tight hydrophobic core inhibits correct 3-D rearrangement. Only one of the six NLE trajectories folded, and it showed no flipping transition. All the other trajectories get trapped in hydrophobically collapsed states. The NLE residues were found to be buried deeply into the core, compared to the corresponding lysines in the villin headpiece, thereby making the core tighter and harder to undo for 3-D rearrangement. Our results suggest that the NLE may not be a fast folder as experiments suggest. The tightness of the hydrophobic core may be a very important factor in the folding of larger proteins. It is likely that chaperones like GroEL act to undo the tight hydrophobic core of proteins, after most secondary structure elements have been formed, so that global rearrangement is easier. I conclude by presenting facts about chaperone-protein complexes and propose further directions for the study of protein folding.

STRUCTURE AND FUNCTION OF THE GROUP III CHAPERONINS, A UNIQUE CLADE OF PROTEIN FOLDING NANOMACHINES

The survival and descent of cells is universally dependent on maintaining their proteins in a properly folded condition. It is widely accepted that the information for the folding of the nascent polypeptide chain into a native protein is encrypted in the amino acid sequence, and the Nobel Laureate Christian Anfinsen was the first to demonstrate that a protein could spontaneously refold after complete unfolding. However, it became clear that the observed folding rates for many proteins were much slower than rates estimated in vivo. This led to the recognition of required protein-protein interactions that promote proper folding. A unique group of proteins, the molecular chaperones, are responsible for maintaining protein homeostasis during normal growth as well as stress conditions. Chaperonins (CPNs) are ubiquitous and essential chaperones. They form ATP-dependent, hollow complexes that encapsulate polypeptides in two back-to-back stacked multisubunit rings, facilitating protein folding through highly cooperative allosteric articulation. CPNs are usually classified into Group I and Group II. Here, I report the characterization of a novel CPN belonging to a third Group, recently discovered in bacteria. Group III CPNs have close phylogenetic association to the Group II CPNs found in Archaea and Eukarya, and may be a relic of the Last Common Ancestor of the CPN family. The gene encoding the Group III CPN from Carboxydothermus hydrogenoformans and Candidatus Desulforudis audaxviator was cloned in E. coli and overexpressed in order to both characterize the protein and to demonstrate its ability to function as an ATPase chaperone. The opening and closing cycle of the Chy chaperonin was examined via site-directed mutations affecting the ATP binding site at R155. To relate the mutational analysis to the structure of the CPN, the crystal structure of both the AMP-PNP (an ATP analogue) and ADP bound forms were obtained in collaboration with Sun-Shin Cha in Seoul, South Korea. The ADP and ATP binding site substitutions resulted in frozen forms of the structures in open and closed conformations. From this, mutants were designed to validate hypotheses regarding key ATP interacting sites as well as important stabilizing interactions, and to observe the physical properties of the resulting complexes by calorimetry.

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

EPR spectroscopic investigation of membrane protein structure and folding on light harvesting complex LHCIIb

Structure and folding of membrane proteins are important issues in molecular and cell biology. In this work new approaches are developed to characterize the structure of folded, unfolded and partially folded membrane proteins. These approaches combine site-directed spin labeling and pulse EPR techniques. The major plant light harvesting complex LHCIIb was used as a model system. Measurements of longitudinal and transversal relaxation times of electron spins and of hyperfine couplings to neighboring nuclei by electron spin echo envelope modulation(ESEEM) provide complementary information about the local environment of a single spin label. By double electron electron resonance (DEER) distances in the nanometer range between two spin labels can be determined. The results are analyzed in terms of relative water accessibilities of different sites in LHCIIb and its geometry. They reveal conformational changes as a function of micelle composition. This arsenal of methods is used to study protein folding during the LHCIIb self assembly and a spatially and temporally resolved folding model is proposed. The approaches developed here are potentially applicable for studying structure and folding of any protein or other self-assembling structure if site-directed spin labeling is feasible and the time scale of folding is accessible to freeze-quench techniques.

Folding funnels and frustration in off-lattice minimalist protein landscapes

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a β-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.

Lattice model for rapidly folding protein-like heteropolymers.

Protein folding is a relatively fast process considering the astronomical number of conformations in which a protein could find itself. Within the framework of a lattice model, we show that one can design rapidly folding sequences by assigning the strongest attractive couplings to the contacts present in a target native state. Our protein design can be extended to situations with both attractive and repulsive contacts. Frustration is minimized by ensuring that all the native contacts are again strongly attractive. Strikingly, this ensures the inevitability of folding and accelerates the folding process by an order of magnitude. The evolutionary implications of our findings are discussed.

A biochemical analysis into the co-translational folding of a G protein-coupled receptor; GPR35

The folding and targeting of membrane proteins poses a major challenge to the cell, as they must remain insertion competent while their highly hydrophobic transmembrane (TM) domains are transferred from the ribosome, through the aqueous cytosol and into the lipid bilayer. The biogenesis of a mature membrane protein takes place through the insertion and integration into the lipid bilayer. A number of TM proteins have been shown to gain some degree of secondary structure within the ribosome tunnel and to retain this conformation throughout maturation. Although studies into the folding and targeting of a number of membrane proteins have been carried out to date, there is little information on one of the largest class of eukaryotic membrane proteins; the G-protein-coupled receptors (GPCRs). This project studies the early folding events of the human ortholog of GPR35. To analyse the structure of the 1st TM domain, intermediates were generated and assessed by the biochemical method of pegylation (PEG-MAL). A structurally-similar microbial opsin (Bacterioopsin) was also used to investigate the differences in the early protein folding within eukaryotic and prokaryotic translation systems. Results showed that neither the 1st TM domain of GPR35 nor Bacterioopsin were capable of compacting in the ribosome tunnel before their N-terminus reached the ribosome exit point. The results for this assay remained consistent whether the proteins were translated in a eukaryotic or prokaryotic translation system. To examine the communication mechanism between the ribosome, the nascent chain and the protein targeting pathway, crosslinking experiments were carried out using the homobifunctional lysine cross-linker BS3. Specifically, the data generated here show that the nascent chain of GPR35 reaches the ribosomal protein uL23 in an extended conformation and interacts with the SRP protein as it exits the ribosome tunnel. This confirms the role of SRP in the co-translational targeting of GPR35. Using these methods insights into the early folding of GPCRs has been obtained. Further experiments using site-directed mutagenesis to reduce hydrophobicity in the 1st TM domain of GPR35, highlighted the mechanisms by which GPCRs are targeted to the endoplasmic reticulum. Confirming that hydrophobicity within the signal anchor sequence is essential of SRP-dependent targeting. Following the successful interaction of the nascent GPR35 and SRP, GPR35 is successfully targeted to ER membranes, shown here as dog pancreas microsomes (DPMs). Glycosylation of the GPR35 N-terminus was used to determine nascent chain structure as it is inserted into the ER membrane. These glycosylation experiments confirm that TM1 has obtained its compacted state whilst residing in the translocon. Finally, a site-specific cross-linking approach using the homobifunctional cysteine cross-linker, BMH, was used to study the lateral integration of GPR35 into the ER. Cross-linking of GPR35 TM1 and TM2 could be detected adjacent to a protein of ~45kDa, believed to be Sec61α. The loss of this adduct, as the nascent chain extends, showed the lateral movement of GPR35 TM1 from the translocon was dependent on the subsequent synthesis of TM2.

Distinct conformational properties determined by implicit and explicit representation of protein-solvent interactions. An analytical and computer simulation study

In the protein folding problem, solvent-mediated forces are commonly represented by intra-chain pairwise contact energy. Although this approximation has proven to be useful in several circumstances, it is limited in some other aspects of the problem. Here we show that it is possible to achieve two models to represent the chain-solvent system. one of them with implicit and other with explicit solvent, such that both reproduce the same thermodynamic results. Firstly, lattice models treated by analytical methods, were used to show that the implicit and explicitly representation of solvent effects can be energetically equivalent only if local solvent properties are time and spatially invariant. Following, applying the same reasoning Used for the lattice models, two inter-consistent Monte Carlo off-lattice models for implicit and explicit solvent are constructed, being that now in the latter the solvent properties are allowed to fluctuate. Then, it is shown that the chain configurational evolution as well as the globule equilibrium conformation are significantly distinct for implicit and explicit solvent systems. Actually, strongly contrasting with the implicit solvent version, the explicit solvent model predicts: (i) a malleable globule, in agreement with the estimated large protein-volume fluctuations; (ii) thermal conformational stability, resembling the conformational hear resistance of globular proteins, in which radii of gyration are practically insensitive to thermal effects over a relatively wide range of temperatures; and (iii) smaller radii of gyration at higher temperatures, indicating that the chain conformational entropy in the unfolded state is significantly smaller than that estimated from random coil configurations. Finally, we comment on the meaning of these results with respect to the understanding of the folding process. (C) 2009 Elsevier B.V. All rights reserved.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Immunoelectron microscopy provides evidence for the presence of mitochondrial heat shock 10-kDa protein (chaperonin 10) in red blood cells and a variety of secretory granules

Hsp10 (10-kDa heat shock protein, also known as chaperonin 10 or Cpn10) is a co-chaperone for Hsp60 in the protein folding process. This protein has also been shown to be identical to the early pregnancy factor, which is an immunosuppressive growth factor found in maternal serum. In this study we have used immunogold electron microscopy to study the subcellular localization of Hsp10 in rat tissues sections embedded in LR Gold resin employing polyclonal antibodies raised against different regions of human Hsp10. In all rat tissues examined including liver, heart, pancreas, kidney, anterior pituitary, salivary gland, thyroid, and adrenal gland, antibodies to Hsp10 showed strong labeling of mitochondria. However, in a number of tissues, in addition to the mitochondrial labeling, strong and highly specific labeling with the Hsp10 antibodies was also observed in several extramitochondrial compartments. These sites included zymogen granules in pancreatic acinar cells, growth hormone granules in anterior pituitary, and secretory granules in PP pancreatic islet cells. Additionally, the mature red blood cells which lack mitochondria, also showed strong reactivity with the Hsp10 antibodies. The observed labeling with the Hsp10 antibodies, both within mitochondria as well as in other compartments/cells, was abolished upon omission of the primary antibodies or upon preadsorption of the primary antibodies with the purified recombinant human Hsp10. These results provide evidence that similar to a number of other recently described mitochondrial proteins (viz., Hsp60, tumor necrosis factor receptor-associated protein- 1, P32 (gC1q-R) protein, and cytochrome c), Hsp10 is also found at a variety of specific extramitochondrial sites in normal rat tissue. These results raise important questions as to how these mitochondrial proteins are translocated to other compartments and their possible function(s) at these sites. The presence of these proteins at extramitochondrial sites in normal tissues has important implications concerning the role of mitochondria in apoptosis and genetic diseases.

Integrating protein structural information

Dissertação apresentada para obtenção de Grau de Doutor em Bioquímica,Bioquímica Estrutural, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

Role of the host cell's unfolded protein response in arenavirus infection.

Arenaviruses are enveloped RNA viruses with a nonlytic life cycle that cause acute and persistent infections. Here, we investigated the role of the host cell's unfolded protein response (UPR) in infection of the prototypic arenavirus lymphocytic choriomeningitis virus (LCMV). In mammalian cells, the endoplasmic reticulum (ER) chaperone protein GRP78/BiP functions as the principal sensor for the induction of the UPR and interacts with three mediators: kinase/endonuclease inositol-requiring protein 1 (IRE1), PKR-like ER kinase (PERK), and activating transcription factor 6 (ATF6). Acute infection with LCMV resulted in a selective induction of the ATF6-regulated branch of the UPR, whereas pathways controlled by PERK and IRE1 were neither activated nor blocked. Expression of individual LCMV proteins revealed that the viral glycoprotein precursor (GPC), but not that of other viral proteins, was responsible for the induction of ATF6. Rapid downregulation of the viral GPC during transition from acute to persistent LCMV infection restored basal levels of UPR signaling. To address a possible role of ATF6 signaling in LCMV infection, we used cells deficient in site 2 protease (S2P), a metalloprotease required for the activation of ATF6. Cells deficient in S2P showed significantly lower levels of production of infectious virus during acute but not persistent infection, indicating a requirement for ATF6-mediated signaling for optimal virus multiplication. In summary, acute LCMV infection seems to selectively induce the ATF6-regulated branch of the UPR that is likely beneficial for virus replication and cell viability, but it avoids induction of PERK and IRE1, whose activation may be detrimental for virus and the host cell.

Oxidative folding of synthetic polypeptides S-protected as tert-butylthio derivatives.

A new method for oxidative folding of synthetic polypeptides assembled by stepwise solid phase synthesis is introduced. Folding is obtained in excellent yields by reacting S-tert-butylthiolated polypeptides with a 100-fold molar excess of cysteine at 37 degrees C in a slightly alkaline buffer containing chaotropic salts, and in the presence of air-oxygen. This novel protocol has been applied to the folding of S-tert-butylthiolated human thymus and activation-regulated chemokine (hu-TARC) derivatives as well as to larger segments of Plasmodium falciparum and Plasmodium berghei circumsporozoite proteins. Folded P. falciparum polypeptides have been used as substrates of endoproteinase Glu-C (Glu-C) and endoproteinase Asp-N (Asp-N) in an attempt to identify their disulfide connectivities. Particular practical advantages of the present method are (i) easy purification and storage of the S-protected peptide derivatives, (ii) elimination of the risk of cysteine alkylation during the acidolytic cleavage deprotection and resin cleavage steps, (iii) possibility to precisely evaluate the extent of folding and disulfide bond formation by mass spectrometry, and (iv) facile recovery of the final folded product.

Dicarboxylic amino acids and glycine-betaine regulate chaperone-mediated protein-disaggregation under stress.

Active protein-disaggregation by a chaperone network composed of ClpB and DnaK + DnaJ + GrpE is essential for the recovery of stress-induced protein aggregates in vitro and in Escherichia coli cells. K-glutamate and glycine-betaine (betaine) naturally accumulate in salt-stressed cells. In addition to providing thermo-protection to native proteins, we found that these osmolytes can strongly and specifically activate ClpB, resulting in an increased efficiency of chaperone-mediated protein disaggregation. Moreover, factors that inhibited the chaperone network by impairing the stability of the ClpB oligomer, such as natural polyamines, dilution, or high salt, were efficiently counteracted by K-glutamate or betaine. The combined protective, counter-negative and net activatory effects of K-glutamate and betaine, allowed protein disaggregation and refolding under heat-shock temperatures that otherwise cause protein aggregation in vitro and in the cell. Mesophilic organisms may thus benefit from a thermotolerant osmolyte-activated chaperone mechanism that can actively rescue protein aggregates, correctly refold and maintain them in a native state under heat-shock conditions.