Matching Items (11)
Filtering by
- All Subjects: Biophysics
- Creators: Shumway, John
Description
Molecular dynamics (MD) simulations provide a particularly useful approach to understanding conformational change in biomolecular systems. MD simulations provide an atomistic, physics-based description of the motions accessible to biomolecular systems on the pico- to micro-second timescale, yielding important insight into the free energy of the system, the dynamical stability of contacts and the role of correlated motions in directing the motions of the system. In this thesis, I use molecular dynamics simulations to provide molecular mechanisms that rationalize structural, thermodynamic, and mutation data on the interactions between the lac repressor headpiece and its O1 operator DNA as well as the ERK2 protein kinase. I performed molecular dynamics simulations of the lac repressor headpiece - O1 operator complex at the natural angle as well as at under- and overbent angles to assess the factors that determine the natural DNA bending angle. I find both energetic and entropic factors contribute to recognition of the natural angle. At the natural angle the energy of the system is minimized by optimization of protein-DNA contacts and the entropy of the system is maximized by release of water from the protein-DNA interface and decorrelation of protein motions. To identify the mechanism by which mutations lead to auto-activation of ERK2, I performed a series of molecular dynamics simulations of ERK1/2 in various stages of activation as well as the constitutively active Q103A, I84A, L73P and R65S ERK2 mutants. My simulations indicate the importance of domain closure for auto-activation and activity regulation. My results enable me to predict two loss-of-function mutants of ERK2, G83A and Q64C, that have been confirmed in experiments by collaborators. One of the powerful capabilities of MD simulations in biochemistry is the ability to find low free energy pathways that connect and explain disparate structural data on biomolecular systems. An extention of the targeted molecular dynamics technique using constraints on internal coordinates will be presented and evaluated. The method gives good results for the alanine dipeptide, but breaks down when applied to study conformational changes in GroEL and adenylate kinase.
ContributorsBarr, Daniel Alan (Author) / van der Vaart, Arjan (Thesis advisor) / Matyushov, Dmitry (Committee member) / Wolf, George (Committee member) / Shumway, John (Committee member) / Arizona State University (Publisher)
Created2011
Description
This thesis describes several experiments based on carbon nanotube nanofludic devices and field-effect transistors. The first experiment detected ion and molecule translocation through one single-walled carbon nanotube (SWCNT) that spans a barrier between two fluid reservoirs. The electrical ionic current is measured. Translocation of small single stranded DNA oligomers is marked by large transient increases in current through the tube and confirmed by a PCR (polymerase chain reaction) analysis. Carbon nanotubes simplify the construction of nanopores, permit new types of electrical measurement, and open new avenues for control of DNA translocation. The second experiment constructed devices in which the interior of a single-walled carbon nanotube field-effect transistor (CNT-FET) acts as a nanofluidic channel that connects two fluid reservoirs, permitting measurement of the electronic properties of the SWCNT as it is wetted by an analyte. Wetting of the inside of the SWCNT by water turns the transistor on, while wetting of the outside has little effect. This finding may provide a new method to investigate water behavior at nanoscale. This also opens a new avenue for building sensors in which the SWCNT functions as an electronic detector. This thesis also presents some experiments that related to nanofabrication, such as construction of FET with tin sulfide (SnS) quantum ribbon. This work demonstrates the application of solution processed IV-VI semiconductor nanostructures in nanoscale devices.
ContributorsCao, Zhai (Author) / Lindsay, Stuart (Thesis advisor) / Vaiana, Sara (Committee member) / Ros, Robert (Committee member) / Marzke, Robert (Committee member) / Shumway, John (Committee member) / Arizona State University (Publisher)
Created2011
Description
Nanofluidic devices in which one single-walled carbon nanotube (SWCNT) spans a barrier between two fluid reservoirs were constructed, enabling direct electrical measurement of the transport of ions and molecules. Ion current through these devices is about 2 orders of magnitude larger than that predicted from the bulk resistivity of the electrolyte. Electroosmosis drives excess current, carried by cations, and is found to be the origin of giant ionic current through SWCNT as shown by building an ionic field-effect transistor with a gate electrode embedded in the fluid barrier. Wetting of inside of the semi-conducting SWCNT by water showed the change of its electronic property, turning the electronic SWCNT field-effect transistor to "on" state. These findings provide a new method to investigate and control the ion and molecule behavior at nanoscale.
ContributorsPang, Pei (Author) / Lindsay, Stuart (Thesis advisor) / Ros, Robert (Committee member) / Shumway, John (Committee member) / Tao, Nongjian (Committee member) / Menéndez, Jose (Committee member) / Arizona State University (Publisher)
Created2011
Description
Proteins are a fundamental unit in biology. Although proteins have been extensively studied, there is still much to investigate. The mechanism by which proteins fold into their native state, how evolution shapes structural dynamics, and the dynamic mechanisms of many diseases are not well understood. In this thesis, protein folding is explored using a multi-scale modeling method including (i) geometric constraint based simulations that efficiently search for native like topologies and (ii) reservoir replica exchange molecular dynamics, which identify the low free energy structures and refines these structures toward the native conformation. A test set of eight proteins and three ancestral steroid receptor proteins are folded to 2.7Å all-atom RMSD from their experimental crystal structures. Protein evolution and disease associated mutations (DAMs) are most commonly studied by in silico multiple sequence alignment methods. Here, however, the structural dynamics are incorporated to give insight into the evolution of three ancestral proteins and the mechanism of several diseases in human ferritin protein. The differences in conformational dynamics of these evolutionary related, functionally diverged ancestral steroid receptor proteins are investigated by obtaining the most collective motion through essential dynamics. Strikingly, this analysis shows that evolutionary diverged proteins of the same family do not share the same dynamic subspace. Rather, those sharing the same function are simultaneously clustered together and distant from those functionally diverged homologs. This dynamics analysis also identifies 77% of mutations (functional and permissive) necessary to evolve new function. In silico methods for prediction of DAMs rely on differences in evolution rate due to purifying selection and therefore the accuracy of DAM prediction decreases at fast and slow evolvable sites. Here, we investigate structural dynamics through computing the contribution of each residue to the biologically relevant fluctuations and from this define a metric: the dynamic stability index (DSI). Using DSI we study the mechanism for three diseases observed in the human ferritin protein. The T30I and R40G DAMs show a loss of dynamic stability at the C-terminus helix and nearby regulatory loop, agreeing with experimental results implicating the same regulatory loop as a cause in cataracts syndrome.
ContributorsGlembo, Tyler J (Author) / Ozkan, Sefika B (Thesis advisor) / Thorpe, Michael F (Committee member) / Ros, Robert (Committee member) / Kumar, Sudhir (Committee member) / Shumway, John (Committee member) / Arizona State University (Publisher)
Created2011
Interactions driving the collapse of islet amyloid polypeptide: implications for amyloid aggregation
Description
Human islet amyloid polypeptide (hIAPP), also known as amylin, is a 37-residue intrinsically disordered hormone involved in glucose regulation and gastric emptying. The aggregation of hIAPP into amyloid fibrils is believed to play a causal role in type 2 diabetes. To date, not much is known about the monomeric state of hIAPP or how it undergoes an irreversible transformation from disordered peptide to insoluble aggregate. IAPP contains a highly conserved disulfide bond that restricts hIAPP(1-8) into a short ring-like structure: N_loop. Removal or chemical reduction of N_loop not only prevents cell response upon binding to the CGRP receptor, but also alters the mass per length distribution of hIAPP fibers and the kinetics of fibril formation. The mechanism by which N_loop affects hIAPP aggregation is not yet understood, but is important for rationalizing kinetics and developing potential inhibitors. By measuring end-to-end contact formation rates, Vaiana et al. showed that N_loop induces collapsed states in IAPP monomers, implying attractive interactions between N_loop and other regions of the disordered polypeptide chain . We show that in addition to being involved in intra-protein interactions, the N_loop is involved in inter-protein interactions, which lead to the formation of extremely long and stable β-turn fibers. These non-amyloid fibers are present in the 10 μM concentration range, under the same solution conditions in which hIAPP forms amyloid fibers. We discuss the effect of peptide cyclization on both intra- and inter-protein interactions, and its possible implications for aggregation. Our findings indicate a potential role of N_loop-N_loop interactions in hIAPP aggregation, which has not previously been explored. Though our findings suggest that N_loop plays an important role in the pathway of amyloid formation, other naturally occurring IAPP variants that contain this structural feature are incapable of forming amyloids. For example, hIAPP readily forms amyloid fibrils in vitro, whereas the rat variant (rIAPP), differing by six amino acids, does not. In addition to being highly soluble, rIAPP is an effective inhibitor of hIAPP fibril formation . Both of these properties have been attributed to rIAPP's three proline residues: A25P, S28P and S29P. Single proline mutants of hIAPP have also been shown to kinetically inhibit hIAPP fibril formation. Because of their intrinsic dihedral angle preferences, prolines are expected to affect conformational ensembles of intrinsically disordered proteins. The specific effect of proline substitutions on IAPP structure and dynamics has not yet been explored, as the detection of such properties is experimentally challenging due to the low molecular weight, fast reconfiguration times, and very low solubility of IAPP peptides. High-resolution techniques able to measure tertiary contact formations are needed to address this issue. We employ a nanosecond laser spectroscopy technique to measure end-to-end contact formation rates in IAPP mutants. We explore the proline substitutions in IAPP and quantify their effects in terms of intrinsic chain stiffness. We find that the three proline mutations found in rIAPP increase chain stiffness. Interestingly, we also find that residue R18 plays an important role in rIAPP's unique chain stiffness and, together with the proline residues, is a determinant for its non-amyloidogenic properties. We discuss the implications of our findings on the role of prolines in IDPs.
ContributorsCope, Stephanie M (Author) / Vaiana, Sara M (Thesis advisor) / Ghirlanda, Giovanna (Committee member) / Ros, Robert (Committee member) / Lindsay, Stuart M (Committee member) / Ozkan, Sefika B (Committee member) / Arizona State University (Publisher)
Created2013
Description
Single molecule DNA Sequencing technology has been a hot research topic in the recent decades because it holds the promise to sequence a human genome in a fast and affordable way, which will eventually make personalized medicine possible. Single molecule differentiation and DNA translocation control are the two main challenges in all single molecule DNA sequencing methods. In this thesis, I will first introduce DNA sequencing technology development and its application, and then explain the performance and limitation of prior art in detail. Following that, I will show a single molecule DNA base differentiation result obtained in recognition tunneling experiments. Furthermore, I will explain the assembly of a nanofluidic platform for single strand DNA translocation, which holds the promised to be integrated into a single molecule DNA sequencing instrument for DNA translocation control. Taken together, my dissertation research demonstrated the potential of using recognition tunneling techniques to serve as a general readout system for single molecule DNA sequencing application.
ContributorsLiu, Hao (Author) / Lindsay, Stuart M (Committee member) / Yan, Hao (Committee member) / Levitus, Marcia (Committee member) / Arizona State University (Publisher)
Created2013
Description
Single molecules in a tunnel junction can now be interrogated reliably using chemically-functionalized electrodes. Monitoring stochastic bonding fluctuations between a ligand bound to one electrode and its target bound to a second electrode ("tethered molecule-pair" configuration) gives insight into the nature of the intermolecular bonding at a single molecule-pair level, and defines the requirements for reproducible tunneling data. Importantly, at large tunnel gaps, there exists a regime for many molecules in which the tunneling is influenced more by the chemical identity of the molecules than by variability in the molecule-metal contact. Functionalizing a pair of electrodes with recognition reagents (the "free analyte" configuration) can generate a distinct tunneling signal when an analyte molecule is trapped in the gap. This opens up a new interface between chemistry and electronics with immediate implications for rapid sequencing of single DNA molecules.
ContributorsChang, Shuai (Author) / Lindsay, Stuart (Thesis advisor) / Ros, Robert (Committee member) / Zhang, Peiming (Committee member) / Tao, Nongjian (Committee member) / Shumway, John (Committee member) / Arizona State University (Publisher)
Created2012
Description
The photosynthetic reaction center is a type of pigment-protein complex found widely in photosynthetic bacteria, algae and higher plants. Its function is to convert the energy of sunlight into a chemical form that can be used to support other life processes. The high efficiency and structural simplicity make the bacterial reaction center a paradigm for studying electron transfer in biomolecules. This thesis starts with a comparison of the primary electron transfer process in the reaction centers from the Rhodobacter shperoides bacterium and those from its thermophilic homolog, Chloroflexus aurantiacus. Different temperature dependences in the primary electron transfer were found in these two type of reaction centers. Analyses of the structural differences between these two proteins suggested that the excess surface charged amino acids as well as a larger solvent exposure area in the Chloroflexus aurantiacus reaction center could explain the different temperature depenence. The conclusion from this work is that the electrostatic interaction potentially has a major effect on the electron transfer. Inspired by these results, a single point mutant was designed for Rhodobacter shperoides reaction centers by placing an ionizable amino acid in the protein interior to perturb the dielectrics. The ionizable group in the mutation site largely deprotonated in the ground state judging from the cofactor absorption spectra as a function of pH. By contrast, a fast charge recombination assoicated with protein dielectric relaxation was observed in this mutant, suggesting the possibility that dynamic protonation/deprotonation may be taking place during the electron transfer. The fast protein dielectric relaxation occuring in this mutant complicates the electron transfer pathway and reduces the yield of electron transfer to QA. Considering the importance of the protein dielectric environment, efforts have been made in quantifying variations of the internal field during charge separation. An analysis protocol based on the Stark effect of reaction center cofactor spectra during charge separation has been developed to characterize the charge-separated radical field acting on probe chromophores. The field change, monitored by the dynamic Stark shift, correlates with, but is not identical to, the electron transfer kinetics. The dynamic Stark shift results have lead to a dynamic model for the time-dependent dielectric that is complementary to the static dielectric asymmetry observed in past steady state experiments. Taken together, the work in this thesis emphasizes the importance of protein electrostatics and its dielectric response to electron transfer.
ContributorsGuo, Zhi (Author) / Woodbury, Neal W (Thesis advisor) / Lindsay, Stuart M (Committee member) / Ross, Robert (Committee member) / Ozkan, Banu S (Committee member) / Moore, Thomas A. (Committee member) / Arizona State University (Publisher)
Created2012
Description
The goal of this theoretical study of infrared spectra was to ascertain to what degree molecules may be identified from their IR spectra and which spectral regions are best suited for this purpose. The frequencies considered range from the lowest frequency molecular vibrations in the far-IR, terahertz region (below ~3 THz or 100 cm-1) up to the highest frequency vibrations (~120 THz or 4000 cm-1). An emphasis was placed on the IR spectra of chemical and biological threat molecules in the interest of detection and prevention. To calculate IR spectra, the technique of normal mode analysis was applied to organic molecules ranging in size from 8 to 11,352 atoms. The IR intensities of the vibrational modes were calculated in terms of the derivative of the molecular dipole moment with respect to each normal coordinate. Three sets of molecules were studied: the organophosphorus G- and V-type nerve agents and chemically related simulants (15 molecules ranging in size from 11 to 40 atoms); 21 other small molecules ranging in size from 8 to 24 atoms; and 13 proteins ranging in size from 304 to 11,352 atoms. Spectra for the first two sets of molecules were calculated using quantum chemistry software, the last two sets using force fields. The "middle" set used both methods, allowing for comparison between them and with experimental spectra from the NIST/EPA Gas-Phase Infrared Library. The calculated spectra of proteins, for which only force field calculations are practical, reproduced the experimentally observed amide I and II bands, but they were shifted by approximately +40 cm-1 relative to experiment. Considering the entire spectrum of protein vibrations, the most promising frequency range for differentiating between proteins was approximately 600-1300 cm-1 where water has low absorption and the proteins show some differences.
ContributorsMott, Adam J (Author) / Rez, Peter (Thesis advisor) / Ozkan, Banu (Committee member) / Shumway, John (Committee member) / Thorpe, Michael (Committee member) / Vaiana, Sara (Committee member) / Arizona State University (Publisher)
Created2012
Description
In eukaryotes, DNA is packed in a highly condensed and hierarchically organized structure called chromatin, in which DNA tightly wraps around the histone octamer consisting of one histone 3-histone 4 (H3-H4) tetramer and two histone 2A- histone 2B (H2A-H2B) dimers with 147 base pairs in an almost two left handed turns. Almost all DNA dependent cellular processes, such as DNA duplication, transcription, DNA repair and recombination, take place in the chromatin form. Based on the critical importance of appropriate chromatin condensation, this thesis focused on the folding behavior of the nucleosome array reconstituted using different templates with various controllable factors such as histone tail modification, linker DNA length, and DNA binding proteins. Firstly, the folding behaviors of wild type (WT) and nucleosome arrays reconstituted with acetylation on the histone H4 at lysine 16 (H4K16 (Ac)) were studied. In contrast to the sedimentation result, atomic force microscopy (AFM) measurements revealed no apparent difference in the compact nucleosome arrays between WT and H4K16 (Ac) and WT. Instead, an optimal loading of nucleosome along the template was found necessary for the Mg2+ induced nucleosome array compaction. This finding leads to the further study on the role of linker DNA in the nucleosome compaction. A method of constructing DNA templates with varied linker DNA lengths was developed, and uniformly and randomly spaced nucleosome arrays with average linker DNA lengths of 30 bp and 60 bp were constructed. After comprehensive analyses of the nucleosome arrays' structure in mica surface, the lengths of the linker DNA were found playing an important role in controlling the structural geometries of nucleosome arrays in both their extended and compact forms. In addition, higher concentration of the DNA binding domain of the telomere repeat factor 2 (TRF2) was found to stimulate the compaction of the telomeric nucleosome array. Finally, AFM was successfully applied to investigate the nucleosome positioning behaviors on the Mouse Mammary Tumor Virus (MMTV) promoter region, and two highly positioned region corresponded to nucleosome A and B were identified by this method.
ContributorsFu, Qiang (Author) / Lindsay, Stuart M (Thesis advisor) / Yan, Hao (Committee member) / Ghirlanda, Giovanna (Committee member) / Arizona State University (Publisher)
Created2010