Understanding how proteins fold and bind is interesting since these processes are central to most biological activity. Protein folding and protein-protein interaction are by themselves very complex but using a good and robust system to study them could ease some of the hurdles. In this thesis I have tried to answer some of the fundamental questions of protein folding and binding. I chose to work with PDZ domains, which are protein domains consisting of 90-100 amino acids. They are found in more than 400 human proteins and function mostly as protein-protein interaction units. These proteins are very stable, easy to express and purify and their folding reaction is reversible under most laboratory conditions. I have characterized the interaction of PSD-95 PDZ3 domain with its putative ligand under different experimental conditions and found out that its binding kinetics is sensitive to salt and pH.  I also demonstrated that the two conserved residues R318 and H372 in PDZ3 are responsible for the salt and pH effect, respectively, on the binding reaction. Moreover, I determined that for PSD 95 PDZ3 coupling of distal residues to peptide binding was better described by a distance relationship and there was a very weak evidence of an allosteric network. Further, I showed that another PDZ domain, SAP97 PDZ2 undergoes conformational change upon ligand binding. Also, I characterized the binding mechanism of a dimeirc ligand/PDZ1-2 tandem interaction and showed that despite its apparent complexity the binding reaction is best described by a square scheme. Additionally, I determined that for the SAP 97 PDZ/HPV E6 interaction that all three PDZ domains each bind one molecule of the E6 protein and that a set of residues in the PDZ2 of SAP 97 could operate in an unexpected long-range manner during E6 interaction. Finally, I showed that perhaps all members in the PDZ family could fold via a three state folding mechanism. I characterized the folding mechanism of five different PDZ domains having similar overall fold but different primary structure and the results indicate that all five fold via an intermediate with two transition states. Transition state one is rate limiting at low denaturant concentration and vice versa for transition state two. Comparing and characterizing the structures of the transition states of two PDZ domains using phi value analysis indicated that their early transition states are less similar as compared to their late transition states.

Proteins are important building blocks of all living organisms. They are composed of a defined sequence of different amino acids, and fold into a specific three-dimensional, ordered structure. The three-dimensional structure largely determines the function of the protein, but protein function always requires motion. Small movements within the protein structure govern the functional properties, and this thesis aims to better understand these discrete protein movements. The motions within the protein structure are governed by thermodynamics, which therefore is useful to predict protein interactions. Nuclear magnetic resonance (NMR) is a powerful tool to study proteins at atomic resolution. Therefore, NMR is the primary method used within this thesis, along with other biophysical techniques such as Fluorescence spectroscopy, Circular Dichroism spectroscopy and in silico modeling. In paper I, NMR in combination with molecular engineering is used to show that the folding of the catalytical subdomains of the enzyme Adenylate kinase does not affect the core of the protein, and thus takes a first step to linking folding, thermodynamic stability and catalysis. In paper II, the structure of the primary allergen from Brazil nut, Ber e 1, is presented along with biophysical measurements that help explain the allergenic potential of the protein. Paper III describes the need for a specific Brazil nut lipid fraction needed to induce an allergenic response. NMR and fluorescence spectroscopy is used to show that there is a direct interaction between Ber e 1 and one or several components in the lipid fraction.

Proteins present at the synapse need to be multitasking in order to perform all vital functions in this limited space. In this thesis I have analyzed the function and evolution of such proteins, focusing on the PDZ domain and the paralemmin family. The PDZ domains bind to a wide variety of interaction partners. The affinity for each partner is regulated by residues at the binding site, but also through intradomain allostery. How this intradomain allostery is transferred to the binding site is not established. I here show that side chain interactions can explain all transfer of intradomain allostery in three analyzed PDZ domains. A circularly permuted PDZ domain has an identical set of amino acids as the original protein and a very similar structure with only a few perturbed side chains. By using the circular permutant I show that a slight alteration in the position of a side chain leads to a corresponding change in allosteric signal. I further study the folding of several PDZ domains and show that they all fold via a conserved folding mechanism, supporting the notion that the final structure has a part in deciding folding mechanism. The folding mechanism of the circularly permuted PDZ domain is conserved compared to the original protein illustrating how circular permutations can be tolerated through evolution. The multifunctionality of paralemmins probably lies in their highly flexible structures. I have studied the evolution of the paralemmins and found that the four mammalian paralemmins arose in the two whole-genome duplications that occurred early in the vertebrate evolution. The fact that all four paralemmins have survived evolution since the gene duplications suggests that they have important functions, possibly in the development of the nervous system. Synaptic proteins are crucial for many biological processes, and their misfolding implicated in many diseases. The results presented here shed light on the mechanisms of action of the synaptic proteins and will help us to understand how they generate disease.

Misfolding of membrane proteins as a result of mutations that disrupt their functions in substrate transport across the membrane or signal transduction is the cause of many significant human diseases. Yet, we still have a limited understanding of the direct consequences of these mutations on folding and function - a necessary step toward the rational design of corrective therapeutics. This thesis addresses the gap in understanding the residue-specific implications for folding through a series of experiments that utilize the cystic fibrosis transmembrane conductance regulator (CFTR) as a model in various contexts. We first examined the thermodynamic implications of mutations in the soluble nucleotide binding domain 1 (NBD1) of CFTR. We found that mutations can have a significant effect on thermodynamic stability that is masked in non-physiological conditions. Our studies were then focussed on a membrane-embedded hairpin CFTR fragment comprised of transmembrane segments 3 (TM3) and 4 (TM4) to evaluate the direct effects of mutations on folding in a systematic manner. It was found that the translocon-mediated membrane insertion of helices closely parallels a basic hydrophobic-aqueous partitioning event. This study was then extended to determine residue-specific effects on helix-helix association. We found that this process is not solely dependent on hydropathy, but there is a context dependence of these results with regard to residue position within the helix. Overall, these findings constitute a key step in relating mutation-derived effects on membrane protein folding to the underlying basis of human disease such as cystic fibrosis.

Recombinant human macrophage colony stimulating factor beta (rhm-CSFβ) is a glycoprotein that stimulates the proliferation, differentiation and survival of cells belonging to the monocyte-macrophage lineage. It contains nine inter-subunit and intra-subunit disulfide bonds and represents an excellent model system for studying disulfide bond formation during protein folding because the assembly of its monomeric subunits and the maturation of its biological activity depend on the progressive formation of the correct disulfide structure during in vitro folding. Knowledge obtained from these studies can be potentially useful in understanding the roles of disulfide bond formation during protein folding in general. rhm-CSF8 was modified by partial reduction of disulfide bonds, yielding CN¹⁵⁷'¹⁵⁹-modified rhm-CSFβ. The modification did not affect the biological activity, stability, or the overall conformation of the protein. However, the C-terminal regions near the modification sites were shown to exhibit faster deuterium exchange behavior as a result of the chemical modification, indicating that the C-terminal regions became more flexible. Folding kinetics of rhm-CSFβ and CN¹⁵⁷'¹⁵⁹-modified rhm-CSFβ were shown to be essentially the same, suggesting that the modification did not affect the folding kinetics of the oxidized rhm-CSFβ. The denatured and reduced rhm-CSFβ was refolded with the aid of a chemical oxidant. The data indicated that the in vitro folding rhm-CSFβ proceeded via multiple pathways involving monomeric and dimeric intermediates. Disulfide bond shuffling catalyzed by GSH/GSSG represented an important isomerization step in folding. A dimeric intermediate, D-SS8-cam2, was isolated and identified as a kinetic trap, perhaps requiring significant structural arrangement to convert to the native protein. The heterogeneous folding mixture detected by both disulfide bond quenching and H/D pulsed labeling indicate that rhm -CSFβ folding is a diffusion like process as described by the folding funnel model. / Graduation date: 2003

We show that a single-copy tempering method is useful in protein-folding simulations of large scale and high accuracy (explicit solvent, atomic representation, and physics-based potential). The method uses a runtime estimate of the average potential energy from an integral identity to guide a random walk in the continuous temperature space. It was used for folding three mini-proteins, trpzip2 (PDB ID: 1LE1), trp-cage (1L2Y), and villin headpiece (1VII) within atomic accuracy. Further, using a modification of the method with a dihedral bias potential added on the roof temperature, we were able to fold four larger helical proteins: α3D (2A3D), α3W (1LQ7), Fap1-NRα (2KUB) and S-836 (2JUA). We also discuss how to optimally use simulation data through an integral identity. With the help of a general mean force formula, the identity makes better use of data collected in a molecular dynamics simulation and is more accurate and precise than the common histogram approach.

<p>Understanding the interconversion between the thermodynamically distinguishable states present in a protein folding pathway provides not only the kinetics and energetics of protein folding but also insights into the functional roles of these states in biological systems. The protein component of bacterial RNase P holoenzyme from Bacillus subtilis (P protein) was used as a model system to elucidate the general folding/unfolding of an intrinsically unstructured protein (IUP) both in the absence and presence of ligands.</p><p>P protein was previously characterized as an intrinsically unstructured protein, and it is predominantly unfolded in the absence of ligands. Addition of small anions can induce the protein to fold. Therefore, the folding and binding are tightly coupled. Trimethylamine-N oxide (TMAO), an osmolyte that stabilizes the unliganded folded form of the protein, enabled us to study the folding process of P protein in the absence of ligand. Transient stopped-flow kinetic time courses at various final TMAO concentrations showed multiphase kinetics. Equilibrium "cotitration" experiments were performed using both TMAO and urea to obtain a TMAO-urea titration surface of P protein. Both kinetic and equilibrium studies show evidence of an intermediate state in the P protein folding process. The intermediate state is significantly populated and the folding rate constants involved in the reaction are slow relative to similar size proteins. </p><p>NMR spectroscopy was used to characterize the structural properties of the folding intermediate of P protein. The results indicate that the N-terminal (residues 2-19) and C-terminal regions (residues 91-116, 118 is the last residue) are mostly unfolded. 1H-15N HSQC NMR spectra were collected at various pH values. The results suggest that His 22 may play a major role in the energetics of the equilibria between the unfolded, intermediate, and native states of P protein.</p><p>Ligand-induced folding kinetics were also investigated to elucidate the overall coupled folding and binding mechanism of P protein and the holoenzyme assembly process. Stopped flow fluorescence experiments were performed at various final ligand concentrations and the data were analyzed using a minimal complexity model that included three conformational states (unfolded, intermediate and folded) in each of three possible liganding states (0, 1 and 2 ligands). The kinetic and equilibrium model parameters that best fit the data were used to calculate the flux through each of the six possible folding/binding pathways. This novel flux-based analysis allows evaluation of the relative importance of pathways in which folding precedes binding or vice versa. The results indicate that the coupled folding and binding mechanism of P protein is strongly dependent on ligand concentration. This conclusion can be generalized to other protein systems for which ligand binding is coupled to conformational changes.</p> / Dissertation

RNA and protein molecules undergo a dynamic folding process that is important to their function. Computational methods are critical for studying this folding pro- cess because it is difficult to observe experimentally. In this work, we introduce new computational techniques to study RNA and protein energy landscapes, includ- ing a method to approximate an RNA energy landscape with a coarse graph (map) and new tools for analyzing graph-based approximations of RNA and protein energy landscapes. These analysis techniques can be used to study RNA and protein fold- ing kinetics such as population kinetics, folding rates, and the folding of particular subsequences. In particular, a map-based Master Equation (MME) method can be used to analyze the population kinetics of the maps, while another map analysis tool, map-based Monte Carlo (MMC) simulation, can extract stochastic folding pathways from the map. To validate the results, I compared our methods with other computational meth- ods and with experimental studies of RNA and protein. I first compared our MMC and MME methods for RNA with other computational methods working on the com- plete energy landscape and show that the approximate map captures the major fea- tures of a much larger (e.g., by orders of magnitude) complete energy landscape. Moreover, I show that the methods scale well to large molecules, e.g., RNA with 200+ nucleotides. Then, I correlate the computational results with experimental findings. I present comparisons with two experimental cases to show how I can pre- dict kinetics-based functional rates of ColE1 RNAII and MS2 phage RNA and their mutants using our MME and MMC tools respectively. I also show that the MME and MMC tools can be applied to map-based approximations of protein energy energy landscapes and present kinetics analysis results for several proteins.

At first glance, robots and proteins have little in common. Robots are commonly thought of as tools that perform tasks such as vacuuming the floor, while proteins play essential roles in many biochemical processes. However, the functionality of both robots and proteins is highly dependent on their motions. In order to study motions in these two divergent domains, the same underlying algorithmic framework can be applied. This method is derived from probabilistic roadmap methods (PRMs) originally developed for robotic motion planning. It builds a graph, or roadmap, where configurations are represented as vertices and transitions between configurations are edges. The contribution of this work is a set of intelligent methods applied to PRMs. These methods facilitate both the modeling and analysis of motions, and have enabled the study of complex and high-dimensional problems in both robotic and molecular domains. In order to efficiently study biologically relevant molecular folding behaviors we have developed new techniques based on Monte Carlo solution, master equation calculation, and non-linear dimensionality reduction to run simulations and analysis on the roadmap. The first method, Map-based master equation calculation (MME), extracts global properties of the folding landscape such as global folding rates. On the other hand, another method, Map-based Monte Carlo solution (MMC), can be used to extract microscopic features of the folding process. Also, the application of dimensionality reduction returns a lower-dimensional representation that still retains the principal features while facilitating both modeling and analysis of motion landscapes. A key contribution of our methods is the flexibility to study larger and more complex structures, e.g., 372 residue Alpha-1 antitrypsin and 200 nucleotide ColE1 RNAII. We also applied intelligent roadmap-based techniques to the area of robotic motion. These methods take advantage of unsupervised learning methods at all stages of the planning process and produces solutions in complex spaces with little cost and less manual intervention compared to other adaptive methods. Our results show that our methods have low overhead and that they out-perform two existing adaptive methods in all complex cases studied.

Properly folded proteins are necessary for all living organisms. Incorrectly folded proteins can lead to a variety of diseases such as Alzheimer’s Disease or Bovine Spongiform Encephalitis (Mad Cow Disease). Understanding the forces involved in protein folding is essential to the understanding and treatment of protein misfolding diseases. When proteins fold, a significant amount of surface area is buried in the protein interior. It has long been known that burial of hydrophobic surface area was important to the stability of the folded structure. However, the impact of burying polar surface area is not well understood. Theoretical results suggest that burying polar groups decreases the stability, but experimental evidence supports the belief that polar group burial increases the stability. Studies of tyrosine to phenylalanine mutations have shown the removal of the tyrosine OH group generally decreases stability. Through computational investigations into the effect of buried tyrosine on protein stability, favorable van der Waals interactions are shown to correlate with the change in stability caused by replacing the tyrosine with phenylalanine to remove the polar OH group. Two large-scale studies on nearly 1000 high-resolution x-ray structures are presented. The first investigates the electrostatic and van der Waals interactions, analyzing the energetics of burying various atom groups in the protein interior. The second large-scale study analyzes the packing differences in the interior of the protein and shows that hydrogen bonding increases packing, decreasing the volume of a hydrogen bonded backbone by about 1.5 Å3 per hydrogen bond. Finally, a structural comparison between RNase Sa and a variant in which five lysines replaced five acidic groups to reverse the net charge is presented. It is shown that these mutations have a marginal impact on the structure, with only small changes in some loop regions.