Browsing by Subject "protein folding"
Now showing 1 - 10 of 10
Results Per Page
Sort Options
Item A motion planning approach to protein folding(Texas A&M University, 2004-09-30) Song, GuangProtein folding is considered to be one of the grand challenge problems in biology. Protein folding refers to how a protein's amino acid sequence, under certain physiological conditions, folds into a stable close-packed three-dimensional structure known as the native state. There are two major problems in protein folding. One, usually called protein structure prediction, is to predict the structure of the protein's native state given only the amino acid sequence. Another important and strongly related problem, often called protein folding, is to study how the amino acid sequence dynamically transitions from an unstructured state to the native state. In this dissertation, we concentrate on the second problem. There are several approaches that have been applied to the protein folding problem, including molecular dynamics, Monte Carlo methods, statistical mechanical models, and lattice models. However, most of these approaches suffer from either overly-detailed simulations, requiring impractical computation times, or overly-simplified models, resulting in unrealistic solutions. In this work, we present a novel motion planning based framework for studying protein folding. We describe how it can be used to approximately map a protein's energy landscape, and then discuss how to find approximate folding pathways and kinetics on this approximate energy landscape. In particular, our technique can produce potential energy landscapes, free energy landscapes, and many folding pathways all from a single roadmap. The roadmap can be computed in a few hours on a desktop PC using a coarse potential energy function. In addition, our motion planning based approach is the first simulation method that enables the study of protein folding kinetics at a level of detail that is appropriate (i.e., not too detailed or too coarse) for capturing possible 2-state and 3-state folding kinetics that may coexist in one protein. Indeed, the unique ability of our method to produce large sets of unrelated folding pathways may potentially provide crucial insight into some aspects of folding kinetics that are not available to other theoretical techniques.Item Computational and experimental investigations of forces in protein folding(Texas A&M University, 2005-02-17) Schell, David AndrewProperly 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.Item Experimental and Computational Studies on Protein Folding, Misfolding and Stability(2010-07-14) Wei, YunProteins need fold to perform their biological function. Thus, understanding how proteins fold could be the key to understanding life. In the first study, the stability and structure of several !-hairpin peptide variants derived from the C-terminus of the B1 domain of protein G (PGB1) were investigated by a number of experimental and computational techniques. Our analysis shows that the structure and stability of this hairpin can be greatly affected by one or a few simple mutations. For example, removing an unfavorable charge near the N-terminus of the peptide (Glu42 to Gln or Thr) or optimization of the N-terminal charge-charge interactions (Gly41 to Lys) both stabilize the peptide, even in water. Furthermore, a simple replacement of a charged residue in the turn (Asp47 to Ala) changes the !-turn conformation. Our results indicate that the structure and stability of this !?hairpin peptide can be modulated in numerous ways and thus contributes towards a more complete understanding of this important model !-hairpin as well as to the folding and stability of larger peptides and proteins. The second study revealed that PGB1 and its variants can form amyloid fibrils in vitro under certain conditions and these fibrils resemble those from other proteins that have been implicated in diseases. To gain a further understanding of molecular mechanism of PGB1 amyloid formation, we designed a set of variants with mutations that change the local secondary structure propensity in PGB1, but have similar global conformational stability. The kinetics of amyloid formation of all these variants have been studied and compared. Our results show that different locations of even a single mutation can have a dramatic effect on PGB1 amyloid formation, which is in sharp contrast with a previous report. Our results also suggest that the "-helix in PGB1 plays an important role in the amyloid formation process of PGB1. In the final study, we investigate the forces that contribute to protein stability in a very general manner. Based on what we have learned about the major forces that contribute to the stability of globular proteins, protein stability should increase as the size of the protein increases. This is not observed: the conformational stability of globular proteins is independent of protein size. In an effort to understand why large proteins are not more stable than small proteins, twenty single-domain globular proteins ranging in size from 35 to 470 residues have been analyzed. Our study shows that nature buries more charged groups and more non-hydrogen-bonded polar groups to destabilize large proteins.Item Intelligent Motion Planning and Analysis with Probabilistic Roadmap Methods for the Study of Complex and High-Dimensional Motions(2011-02-22) Tapia, LydiaAt 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.Item Investigation of Peptide Folding by Nuclear Magnetic Resonance Spectroscopy(2012-07-16) Hwang, SoYounUnderstanding structure and folding of a protein is the key to understanding its biological function and potential role in diseases. Despite the importance of protein folding, a molecular level understanding of this process is still lacking. Solution-state nuclear magnetic resonance (NMR) is a powerful technique to investigate protein structure, dynamics, and folding mechanisms, since it provides residue-specific information. One of the major contributions that govern protein structure appears to be the interaction with the solvent. The importance of these interactions is particularly apparent in membrane proteins, which exist in an amphiphilic environment. Here, individual peptide fragments taken from the disulfide bond forming protein B (DsbB) were investigated in various solvents. The alpha-helical structures that were obtained, suggest that DsbB follows the two-stage model for folding. However, side chains of polar residues showed different conformations compared to the X-ray structure of fulllength protein, implying that polar side-chains may re-orient upon helix packing in order to form the necessary tertiary interactions that stabilize the global fold of DsbB. Model peptides in general represent attractive systems for the investigation of non-covalent interactions important for protein folding, including those with the solvent. NMR structures of the water soluble peptide, BBA5, were obtained in the presence an organic co-solvent, methanol. These structures indicate that the addition of methanol stabilizes an alpha-helix segment, but disrupts a hydrophobic cluster forming a beta-hairpin. Since dynamic effects reduce the ability for experimental observation of individual, bound solvent molecules, results were compared with molecular dynamics simulations. This comparison indicates that the observed effects of NMR structures are due to preferred binding of methanol and reduction of peptide-water hydrogen bonding. NMR structures, such as those determined here, represent a distribution of conformations under equilibrium. The dynamic process of protein unfolding can nevertheless be accessed through denaturation. A method was developed to probe thermal denaturation by measuring the temperature dependence of NOE intensity. Applied to a model peptide, trpzip4, it was confirmed that the beta-hairpin structure of this peptide is stabilized by the hydrophobic cluster formed by tryptophan residues. Together, the peptides investigated here illustrate the important roles that solvent-peptide interactions and side chain-side chain hydrophobic interactions play in forming stable secondary and tertiary structures.Item Probing the denatured state ensemble with fluorescence(Texas A&M University, 2004-09-30) Alston, Roy WillisTo understand protein stability and the mechanism of protein folding, it is essential that we gain a better understanding of the ensemble of conformations that make up the denatured state of a protein. The primary goal of the research described here was to see what we might learn about the denatured state using fluorescence. To this end, tryptophan was introduced at five sites in Ribonuclease Sa (RNase Sa): D1W, Y52W, Y55W, T76W, and Y81W. The fluorescent properties of the denatured states of these five proteins were studied and compared to the fluorescent properties of eight model compounds: N-acetyl-tryptophan-amide (NATA), N-acetyl-Ala-Trp-Ala-amide (AWA), N-acetyl-Ala-Ala-Trp-Ala-Ala-amide (AAWAA), and five pentapeptides based on the sequence around the original tryptophan substitutions in RNase Sa. Regardless of the denaturant, ?max for the proteins and model compounds differed very little, 349.3 ? 1.2 nm. However, significant differences were observed in the fluorescence intensity at ?max (IF), suggesting that IF is more sensitive to the immediate environment than ?max. The differences in IF are due in part to quenching by neighboring side chains. More importantly, IF was always significantly greater in the protein than in its corresponding pentapeptide, indicating that the protein exerts an effect on the tryptophan, which cannot be mimicked by the pentapeptide models. Acrylamide and iodide quenching experiments were also performed on the model compounds and proteins. Significant differences in the Stern-Volmer quenching constant (KSV) were also observed between the proteins and between the proteins and their corresponding pentapeptides. Importantly, the KSV for the protein was always less than in its corresponding pentapeptide. These data along with the IF data show that non-local structure in the unfolded state influences tryptophan fluorescence and accessibility. In summary, these and our other studies show that fluorescence can be used to gain a better understanding of the denatured states of proteins.Item Rational design of organophosphorus hydrolase for the degradation and detection of neurotoxic pesticides and chemical warfare agents(Texas A&M University, 2007-09-17) Reeves, Tony ElvernIt is critical to consider the balance between the catalytic capabilities of an enzyme and the inherent structural stability of the protein when developing enzymes for specific applications. Rational site directed mutagenesis has been used to explore the role of residues 254 and 257 in the global stability and catalytic specificities of organophosphorus hydrolase (OPH, EC 3.1.8.1). Substitution of residues H254 and H257, which are located near the active site, had a marked effect on both the global stability and substrate specificity of the enzyme. For example, the for the double mutation CoTG????2+ H254R H257L (RL) enzyme variant was 19.6 kcal/mol, 5.7 kcal/mol less than that of the wild type enzyme. At the same time, the altered enzyme was catalytically more effective against VX and VR (Russian VX), as compared to the wild type enzyme. Limited proteolysis verified the importance of residues 254 and 257 for functional stability, evidenced by enhanced resistance to irreversible unfolding associated with thermal denaturation. It has been possible to construct third generation OPH variants, which are more stable than the wild type enzyme, with a 10 ????C increase in the apparent melting temperature (TM app), yet retained desirable catalytic properties. It appeared that aromatic stacking and cation-???? interactions involving near active site residues not only affected activity but significantly contributed to the chemical and thermal stability of OPH. Rational design was used to develop an enzyme with an optimized orientation on a catalytically active biosensor surface. In these studies, lysine side chains located on the surface of OPH were used to create attachment sites to a surface plasmon resonance sensor resulting in an ensemble of enzyme orientations. Some of these orientations could be functionally restrictive if the active site is oriented toward the sensor surface. Substitution of a lysine near the active site resulted in 20% more activity with 53% less enzyme immobilized, thus increasing the specific activity of the decorated surface 2.5 fold.Item Rigidity Analysis for Modeling Protein Motion(2010-07-14) Thomas, Shawna L.Protein structure and motion plays an essential role in nearly all forms of life. Understanding both protein folding and protein conformational change can bring deeper insight to many biochemical processes and even into some devastating diseases thought to be the result of protein misfolding. Experimental methods are currently unable to capture detailed, large-scale motions. Traditional computational approaches (e.g., molecular dynamics and Monte Carlo simulations) are too expensive to simulate time periods long enough for anything but small peptide fragments. This research aims to model such molecular movement using a motion framework originally developed for robotic applications called the Probabilistic Roadmap Method. The Probabilistic Roadmap Method builds a graph, or roadmap, to model the connectivity of the movable object?s valid motion space. We previously applied this methodology to study protein folding and obtained promising results for several small proteins. Here, we extend our existing protein folding framework to handle larger proteins and to study a broader range of motion problems. We present a methodology for incrementally constructing roadmaps until they satisfy a set of evaluation criteria. We show the generality of this scheme by providing evaluation criteria for two types of motion problems: protein folding and protein transitions. Incremental Map Generation eliminates the burden of selecting a sampling density which in practice is highly sensitive to the protein under study and difficult to select. We also generalize the roadmap construction process to be biased towards multiple conformations of interest thereby allowing it to model transitions, i.e., motions between multiple known conformations, instead of just folding to a single known conformation. We provide evidence that this generalized motion framework models large-scale conformational change more realistically than competing methods. We use rigidity theory to increase the efficiency of roadmap construction by introducing a new sampling scheme and new distance metrics. It is only with these rigidity-based techniques that we were able to detect subtle folding differences between a set of structurally similar proteins. We also use it to study several problems related to protein motion including distinguishing secondary structure formation order, modeling hydrogen exchange, and folding core identification. We compare our results to both experimental data and other computational methods.Item Techniques for modeling and analyzing RNA and protein folding energy landscapes(2009-05-15) Tang, XinyuRNA 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.Item The folding kinetics of ribonuclease Sa and a charge-reversal variant(Texas A&M University, 2005-02-17) Trefethen, Jared M.The primary objective was to study the kinetics of folding of RNase Sa. Wild-type RNase Sa does not contain tryptophan. A tryptophan was substituted at residue 81 (WT*) to allow fluorescence spectroscopy to be used to monitor folding. This tryptophan mutation did not change the stability. An analysis of the folding kinetics of RNase Sa showed two folding phases, indicating the presence of an intermediate and consistent with the following mechanism: D ? I ? N. Both refolding limbs of the chevron plot (abcissa = final conc. of denaturant and ordinate = kinetic rate) had non-zero slopes suggesting that proline isomerization was not rate-limiting. The conformational stability of a charge-reversed variant, WT*(D17R), of a surface exposed residue on RNase Sa has been studied by equilibrium techniques. This mutant with a single amino acid charge reversal of a surface exposed residue resulted in decreased stability. Calculations using Coulomb?s Law suggested that favorable electrostatic interactions in the denatured state were the cause for the decreased stability for the charge-reversed variant. Folding and unfolding kinetic studies were designed and conducted to study the charge-reversal effect. Unfolding kinetics showed a 10-fold increase in the unfolding rate constant for WT*(D17R) over WT* and no difference in the rate of refolding. Kinetics experiments were also conducted at pH 3 where protonation of Asp17 (charge reversal site) would be expected to negate the observed kinetic effect. At pH 3 the kinetics of unfolding of WT* RNase Sa and the WT*(D17R) mutant were more similar. These kinetic results indicate that a single-site charge reversal lowered the free energy of the denatured state as suspected. Additionally, the results showed that the transition state was stabilized as well. These results show that a specific Coulombic interaction lowered the free energy in the denatured and transition state of the charge-reversal mutant, more than in WT*. To our knowledge, this is the first demonstration that a favorable electrostatic interaction in the denatured state ensemble has been shown to influence the unfolding kinetics of a protein.