Engineering topics

Pohorille, Andrew

Publications and source records attributed to Pohorille, Andrew.

At least 37 records · Page 2

NASA NTRS2011

Molecular Dynamics Simulation of the Antiamoebin Ion Channel: Linking Structure and Conductance

Molecular dynamics simulations were carried out in order to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistant with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive.The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be non-conducting. The conductance of the hexamer was estimated to be 115+/-34 pS and 74+/-20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K(+) and Cl(-) with their first solvation shells intact. The free energy barrier encountered by K(+) is only 2.2 kcal/mol whereas Cl(-) encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics. keywords: ion channels, peptaibols, channel conductance, molecular dynamics

Wilson, Michael A.↗

NASA NTRS2010

The Earliest Ion Channels in Protocellular Membranes

Cellular membranes with their hydrophobic interior are virtually impermeable to ions. Bulk of ion transport through them is enabled through ion channels. Ion channels of contemporary cells are complex protein molecules which span the membrane creating a cylindrical pore filled with water. Protocells, which are widely regarded as precursors to modern cells, had similarly impermeable membranes, but the set of proteins in their disposal was much simpler and more limited. We have been, therefore, exploring an idea that the first ion channels in protocellular membranes were formed by much smaller peptide molecules that could spontaneously selfassemble into short-lived cylindrical bundles in a membrane. Earlier studies have shown that a group of peptides known as peptaibols is capable of forming ion channels in lipid bilayers when they are exposed to an electric field. Peptaibols are small, non-genetically encoded peptides produced by some fungi as a part of their system of defense against bacteria. They are usually only 14-20 residues long, which is just enough to span the membrane. Their sequence is characterized by the presence of non-standard amino acids which, interestingly, are also expected to have existed on the early earth. In particular, the presence of 2-aminoisobutyric acid (AIB) gives peptaibols strong helix forming propensities. Association of the helices inside membranes leads to the formation of cylindrical bundles, typically containing 4 to 10 monomers. Although peptaibols are excellent candidates for models of the earliest ion channels their structures, which are stabilized only by van der Waals forces and occasional hydrogen bonds between neighboring helices, are not very stable. Although it might properly reflect protobiological reality, it is also a major obstacle in studying channel behavior. For this reason we focused on two members of the peptaibol family, trichotoxin and antiamoebin, which are characterized by a single conductance level. This indicates that their structures are unique and stable. In addition, it is also believed that the trichotoxin channel displays some selectivity between potassium and chloride ions. This makes trichotoxin and antiamoebin ideal models of the earliest ion channels that could provide insight into the origins of ion conductance and selectivity. In the absence of crystal structure of the trichotoxin and antiamoebin channels, we propose their molecular models based on experimentally determined number of monomers forming the bundles. We use molecular dynamics simulations to validate the models in terms of their conductance and selectivity. On the basis of our simulations we show that the emergence of channels built of small, alpha-helical peptides was protobiologically plausible and did not require highly specific amino acid sequences, which is a convenient evolutionary trait. Despite their simple structure, such channels could possess properties that, at the first sight, appear to require markedly larger complexity. To this end, we will discuss how the amino acid sequence and structure of primitive channels give rise to the phenomena of ionic conductance and selectivity across the earliest cell walls, which were essential functions for the emergence and early evolution of protocells. Furthermore, we will argue that even though architectures of membrane proteins are not nearly as diverse as those of water-soluble proteins, they are sufficiently flexible to adapt readily to the functional demands arising during evolution.

Mijajlovic, Milan↗

NASA NTRS2006

Chance of Necessity: Modeling Origins of Life

The fundamental nature of processes that led to the emergence of life has been a subject of long-standing debate. One view holds that the origin of life is an event governed by chance, and the result of so many random events is unpredictable. This view was eloquently expressed by Jacques Monod in his book Chance or Necessity. In an alternative view, the origin of life is considered a deterministic event. Its details need not be deterministic in every respect, but the overall behavior is predictable. A corollary to the deterministic view is that the emergence of life must have been determined primarily by universal chemistry and biochemistry rather than by subtle details of environmental conditions. In my lecture I will explore two different paradigms for the emergence of life and discuss their implications for predictability and universality of life-forming processes. The dominant approach is that the origin of life was guided by information stored in nucleic acids (the RNA World hypothesis). In this view, selection of improved combinations of nucleic acids obtained through random mutations drove evolution of biological systems from their conception. An alternative hypothesis states that the formation of protocellular metabolism was driven by non-genomic processes. Even though these processes were highly stochastic the outcome was largely deterministic, strongly constrained by laws of chemistry. I will argue that self-replication of macromolecules was not required at the early stages of evolution; the reproduction of cellular functions alone was sufficient for self-maintenance of protocells. In fact, the precise transfer of information between successive generations of the earliest protocells was unnecessary and could have impeded the discovery of cellular metabolism. I will also show that such concepts as speciation and fitness to the environment, developed in the context of genomic evolution also hold in the absence of a genome.

Pohorille, Andrew↗

NASA NTRS2006

A Bioinformatics Facility for NASA

Building on an existing prototype, we have fielded a facility with bioinformatics technologies that will help NASA meet its unique requirements for biological research. This facility consists of a cluster of computers capable of performing computationally intensive tasks, software tools, databases and knowledge management systems. Novel computational technologies for analyzing and integrating new biological data and already existing knowledge have been developed. With continued development and support, the facility will fulfill strategic NASA s bioinformatics needs in astrobiology and space exploration. . As a demonstration of these capabilities, we will present a detailed analysis of how spaceflight factors impact gene expression in the liver and kidney for mice flown aboard shuttle flight STS-108. We have found that many genes involved in signal transduction, cell cycle, and development respond to changes in microgravity, but that most metabolic pathways appear unchanged.

Schweighofer, Karl↗

NASA NTRS2006

Water as a matrix for life

"Follow the water" is the canonical strategy in searching for life in the universe. Conventionally, discussion of this topic is focused on how solvent supports organic chemistry sufficiently rich to seed life. Perhaps more importantly, solvent must promote self-organization of organic matter into functional structures capable of responding to environmental changes. This process is based on non-covalent interactions. They are constantly formed and broken in response to internal and external stimuli. This requires that their strength must be properly tuned. If they were too weak, the system would exhibit undesired, uncontrolled response to natural fluctuations of physical and chemical parameters. If they were too strong kinetics of biological processes would be slow and energetics costly. Non-covalent interactions are strongly mediated by the solvent. Specifically, high dielectric solvents for life are needed for solubility of polar species and flexibility of biological structures stabilized by electrostatic interactions. Water exhibits a remarkable trait that it promotes solvophobic interactions between non-polar species, which are responsible for self-organization phenomena such as the formation of cellular boundary structures, and protein folding and aggregation. Unusual temperature dependence of hydrophobic interactions - they often become stronger as temperature increases - is a consequence of the temperature insensitivity of properties of the liquid water. This contributes to the existence of robust life over a wide temperature range. Water is not the only liquid with favorable properties for supporting life. Other pure liquids or their mixtures that have high dielectric constants and simultaneously support some level of self-organization will be discussed.

Pohorille, Andrew↗

NASA NTRS2006

Mathematical Modeling of the Origins of Life

The emergence of early metabolism - a network of catalyzed chemical reactions that supported self-maintenance, growth, reproduction and evolution of the ancestors of contemporary cells (protocells) was a critical, but still very poorly understood step on the path from inanimate to animate matter. Here, it is proposed and tested through mathematical modeling of biochemically plausible systems that the emergence of metabolism and its initial evolution towards higher complexity preceded the emergence of a genome. Even though the formation of protocellular metabolism was driven by non-genomic, highly stochastic processes the outcome was largely deterministic, strongly constrained by laws of chemistry. It is shown that such concepts as speciation and fitness to the environment, developed in the context of genomic evolution, also held in the absence of a genome.

Pohorille, Andrew↗

NASA NTRS2006

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

NASA NTRS2006

The Origin and Early Evolution of Membrane Proteins

The origin and early evolution of membrane proteins, and in particular ion channels, are considered from the point of view that the transmembrane segments of membrane proteins are structurally quite simple and do not require specific sequences to fold. We argue that the transport of solute species, especially ions, required an early evolution of efficient transport mechanisms, and that the emergence of simple ion channels was protobiologically plausible. We also argue that, despite their simple structure, such channels could possess properties that, at the first sight, appear to require markedly larger complexity. These properties can be subtly modulated by local modifications to the sequence rather than global changes in molecular architecture. In order to address the evolution and development of ion channels, we focus on identifying those protein domains that are commonly associated with ion channel proteins and are conserved throughout the three main domains of life (Eukarya, Prokarya, and Archaea). We discuss the potassium-sodium-calcium superfamily of voltage-gated ion channels, mechanosensitive channels, porins, and ABC-transporters and argue that these families of membrane channels have sufficiently universal architectures that they can readily adapt to the diverse functional demands arising during evolution.

Pohorille, Andrew↗

NASA NTRS2005

Optimal Sampling of a Reaction Coordinate in Molecular Dynamics

Estimating how free energy changes with the state of a system is a central goal in applications of statistical mechanics to problems of chemical or biological interest. From these free energy changes it is possible, for example, to establish which states of the system are stable, what are their probabilities and how the equilibria between these states are influenced by external conditions. Free energies are also of great utility in determining kinetics of transitions between different states. A variety of methods have been developed to compute free energies of condensed phase systems. Here, I will focus on one class of methods - those that allow for calculating free energy changes along one or several generalized coordinates in the system, often called reaction coordinates or order parameters . Considering that in almost all cases of practical interest a significant computational effort is required to determine free energy changes along such coordinates it is hardly surprising that efficiencies of different methods are of great concern. In most cases, the main difficulty is associated with its shape along the reaction coordinate. If the free energy changes markedly along this coordinate Boltzmann sampling of its different values becomes highly non-uniform. This, in turn, may have considerable, detrimental effect on the performance of many methods for calculating free energies.

Pohorille, Andrew↗

NASA NTRS2005

The Origin and Early Evolution of Membrane Proteins

Membrane proteins mediate functions that are essential to all cells. These functions include transport of ions, nutrients and waste products across cell walls, capture of energy and its transduction into the form usable in chemical reactions, transmission of environmental signals to the interior of the cell, cellular growth and cell volume regulation. In the absence of membrane proteins, ancestors of cell (protocells), would have had only very limited capabilities to communicate with their environment. Thus, it is not surprising that membrane proteins are quite common even in simplest prokaryotic cells. Considering that contemporary membrane channels are large and complex, both structurally and functionally, a question arises how their presumably much simpler ancestors could have emerged, perform functions and diversify in early protobiological evolution. Remarkably, despite their overall complexity, structural motifs in membrane proteins are quite simple, with a-helices being most common. This suggests that these proteins might have evolved from simple building blocks. To explain how these blocks could have organized into functional structures, we performed large-scale, accurate computer simulations of folding peptides at a water-membrane interface, their insertion into the membrane, self-assembly into higher-order structures and function. The results of these simulations, combined with analysis of structural and functional experimental data led to the first integrated view of the origin and early evolution of membrane proteins.

Pohorille, Andrew↗

NASA NTRS2005

Water as a matrix for life

Life is based on non-covalent interactions. They might be either specific (enzyme-substrate interactions, selective ion transport) or nonspecific (lipid-lipid and lipid-protein interactions needed for membrane integrity, fusion and division). Their strength needs to be properly tuned, and this is mediated by the solvent. If interactions are too weak, there might be undesired response to natural fluctuations of physical and chemical parameters. If they are too strong it could impede kinetics and energetics of cellular processes. Thus, the solvent must allow for balancing these interactions. Physical and chemical properties of solvent provide strong constraints for life. Water exhibits a remarkable trait that it promotes both solvophobic and solvophilic interactions. Solvophobic interactions; related to high dielectric constant of the solvent) are necessary for self-organization of matter whereas solvophilic interactions are needed to ensure solubility of polar species. Water offers a large temperature domain of stable liquid and the characteristics hydrophobic effects are a consequence of the temperature in sensitivity of essential properties of its liquid state. Water, however, is not the only liquid with these favorable properties. I will compare in detail properties of water and other pure liquids or their mixtures that have a high dielectric constant and simultaneously support self-organization. I will also discuss properties of water that are unfavorable to life (e.g. its chemical activity against polymerization reactions) and close with summarizing what are alternatives to water as a matrix of life in space.

Pohorille, Andrew↗

NASA NTRS2005

Hydrophobic effect at aqueous interfaces

Conceptual basis for hydrophobic effects in bulk water and at aqueous interfaces have similar conceptual basis but often manifests itself differently. Using a wide range of computer simulations as the basis, I will review different forms of hydrophobic effects at a variety of interfaces starting from simple liquid-vapor and water-oil interfaces and progressing to water-membrane interfaces. I will start with discussing how water is organized at different interfaces, stressing both similarities and differences. The main thread is that, as in the bulk liquid, hydrophobic effects have profound influence on conformational equilibria and organization of both small molecules and macromolecules, but the result of this influence is quite different. Specifically, it will be shown that many small, but not necessarily amphiphilic molecules tend to accumulate at the interface and, and this tendency will be explained. Furthermore, I will show that many short peptides that are disordered in water spontaneously fold into well-defined structures in the interfacial environment. Biological implications of this self-organizing effect will be discussed.

Pohorille, Andrew↗

NASA NTRS2003

Structure, Function, Self-Assembly and Origin of Simple Membrane Proteins

Integral membrane proteins perform such essential cellular functions as transport of ions, nutrients and waste products across cell walls, transduction of environmental signals, regulation of cell fusion, recognition of other cells, energy capture and its conversion into high-energy compounds. In fact, 30-40% of genes in modem organisms codes for membrane proteins. Although contemporary membrane proteins or their functional assemblies can be quite complex, their transmembrane fragments are usually remarkably simple. The most common structural motif for these fragments is a bundle of alpha-helices, but occasionally it could be a beta-barrel. In a series of molecular dynamics computer simulations we investigated self-organizing properties of simple membrane proteins based on these structural motifs. Specifically, we studied folding and insertion into membranes of short, nonpolar or amphiphatic peptides. We also investigated glycophorin A, a peptide that forms sequence-specific dimers, and a transmembrane aggregate of four identical alpha-helices that forms an efficient and selective voltage-gated proton channel was investigated. Many peptides are attracted to water-membrane interfaces. Once at the interface, nonpolar peptides spontaneously fold to a-helices. Whenever the sequence permits, peptides that contain both polar and nonpolar amino also adopt helical structures, in which polar and nonpolar amino acid side chains are immersed in water and membrane, respectively. Specific identity of side chains is less important. Helical peptides at the interface could insert into the membrane and adopt a transmembrane conformation. However, insertion of a single helix is unfavorable because polar groups in the peptide become completely dehydrated upon insertion. The unfavorable free energy of insertion can be regained by spontaneous association of peptides in the membrane. The first step in this process is the formation of dimers, although the most common are aggregates of 4-7 helices. The helices could arrange themselves such that they formed pores capable of transporting ions and small molecules across membranes. Stability of transmembrane aggregates of simple proteins is often only marginal and, therefore, it can be regulated by environmental signals or small sequence modifications in the region of interhelical interactions. A key step in the earliest evolution of membrane proteins was the emergence of selectivity for specific substrates. Many channels could become selective if one or only a few properly chosen amino acids are properly placed along the channel, acting as filters or gates. This is a convenient evolutionary solution because it does not require imposing conditions on the whole sequence.

Pohorille, Andrew↗

NASA NTRS2003

Non-Genomic Origins of Proteins and Metabolism

It is proposed that evolution of inanimate matter to cells endowed with a nucleic acid- based coding of genetic information was preceded by an evolutionary phase, in which peptides not coded by nucleic acids were able to self-organize into networks capable of evolution towards increasing metabolic complexity. Recent findings that truly different, simple peptides (Keefe and Szostak, 2001) can perform the same function (such as ATP binding) provide experimental support for this mechanism of early protobiological evolution. The central concept underlying this mechanism is that the reproduction of cellular functions alone was sufficient for self-maintenance of protocells, and that self- replication of macromolecules was not required at this stage of evolution. The precise transfer of information between successive generations of the earliest protocells was unnecessary and, possibly, undesirable. The key requirement in the initial stage of protocellular evolution was an ability to rapidly explore a large number of protein sequences in order to discover a set of molecules capable of supporting self- maintenance and growth of protocells. Undoubtedly, the essential protocellular functions were carried out by molecules not nearly as efficient or as specific as contemporary proteins. Many, potentially unrelated sequences could have performed each of these functions at an evolutionarily acceptable level. As evolution progressed, however proteins must have performed their functions with increasing efficiency and specificity. This, in turn, put additional constraints on protein sequences and the fraction of proteins capable of performing their functions at the required level decreased. At some point, the likelihood of generating a sufficiently efficient set of proteins through a non-coded synthesis was so small that further evolution was not possible without storing information about the sequences of these proteins. Beyond this point, further evolution required coupling between proteins and informational polymers that is characteristic to all known forms of life. The emergence of such coupling must be postulated in any scenario of the origin of life, no matter whether it starts with RNA or proteins. To examine the evolutionary potential of non-genomic systems, a simple, computationally tractable model, which is still capable of capturing the essential features of the real system, has been studied computationally. Both constructive and destructive processes have been introduced into the model in a stochastic manner. Instead of assuming random reaction sets, only a suite of protobiologically plausible reactions has been considered. Peptides have been explicitly considered as protoenzymes and their catalytic efficiencies have been assigned on the basis of biochemical principles and experimental estimates. Simulations have been carried out using a novel approach (The Next Reaction Method) that is appropriate even for very low concentrations of reactants. Studies have focused on global autocatalytic processes and their diversity.

Pohorille, Andrew↗

NASA NTRS2003

Membrane peptides and their role in protobiological evolution

How simple membrane peptides performed such essential protocellular functions as transport of ions and organic matter across membranes separating the interior of the cell from the environment, capture and utilization of energy, and transduction of environmental signals, is a key question in protobiological evolution. On the basis of detailed, molecular-level computer simulations we explain how these peptides fold at water-membrane interfaces, insert into membranes, self-assemble into higher-order structures and acquire functions. We have investigated the interfacial behavior and folding of several peptides built of leucine and glutamine residues and have demonstrated that many of them tend to adopt ordered structures. Further, we have studied the insertion of an alpha-helical peptide containing leucine (L) and serine (S) of the form (LSLLLSL)3 into a model membrane. The transmembrane state is metastable, and approximately 15 kcal mol(-1) is required to insert the peptide into the membrane. Investigations of dimers formed by (LSLLLSL)3 and glycophorin A demonstrate how the favorable free energy of helix association can offset the unfavorable free energy of insertion, leading to self-assembly of peptide helices in the membrane. An example of a self-assembled structure is the tetrameric transmembrane pore of the influenza virus M2 protein, which is an efficient and selective voltage-gated proton channel. Our simulations explain the gating mechanism and provide guidelines how to re-engineer the channel to act as a simple proton pump. In general, emergence of integral membrane proteins appears to be quite feasible and may be easier to envision than the emergence of water-soluble proteins.

Models, Chemical↗

NASA NTRS2003

Research in Computational Astrobiology

We report on several projects in the field of computational astrobiology, which is devoted to advancing our understanding of the origin, evolution and distribution of life in the Universe using theoretical and computational tools. Research projects included modifying existing computer simulation codes to use efficient, multiple time step algorithms, statistical methods for analysis of astrophysical data via optimal partitioning methods, electronic structure calculations on water-nuclei acid complexes, incorporation of structural information into genomic sequence analysis methods and calculations of shock-induced formation of polycylic aromatic hydrocarbon compounds.

Chaban, Galina↗

NASA NTRS2002

Incorporating Biological Knowledge into Evaluation of Casual Regulatory Hypothesis

Biological data can be scarce and costly to obtain. The small number of samples available typically limits statistical power and makes reliable inference of causal relations extremely difficult. However, we argue that statistical power can be increased substantially by incorporating prior knowledge and data from diverse sources. We present a Bayesian framework that combines information from different sources and we show empirically that this lets one make correct causal inferences with small sample sizes that otherwise would be impossible.

Chrisman, Lonnie↗

NASA NTRS2002

Peptides at Membrane Surfaces and their Role in Prebiotic Evolution

Protocells had to transport ions and organic matter across membranes separating the interior of the cell from the environment, capture and utilize energy and transduce environmental signals. In a series of detailed, molecular-level computer simulations we show how these peptides in contact with membranes can acquire ordered structures and functions. We have investigated the stability of a simple alpha-helical peptide containing Leucine (L) and Serine (S) of the form (LSLLLSL)3 in a model membrane system. The parallel in-plane state is the most stable configuration. The transmembrane state is metastable, and about 15 kcal/mol is required to insert the peptide into the membrane. We investigated dimes of both (LSLLLSL)3 and glycophorin A, and show how the free energy of helix association can, at least partially, offset the free energy of insertion. We have also investigated the transmembrane pore of the influenza M2 protein. This aggregate of four identical alpha-helices, each built of 25 amino acids, forms an efficient and selective voltage-gated proton channel. Our simulations explain the gating mechanism, which can involve strands of hydrogen-bonded water through the pore or proton transfer through tautomerization of protein residues. The channel can be re-engineered to act as a simple proton pump.

Pohorille, Andrew↗