Engineering topics

Pohorille, Andrzej

Publications and source records attributed to Pohorille, Andrzej.

7 records

NASA NTRS2018

Simple Ion Channels: From Structure to Electrophysiology and Back

A reliable way to establish whether our understanding of a channel is satisfactory is to reproduce its measured ionic conductance over a broad range of applied voltages in computer simulations. In molecular dynamics (MD), this can be done by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. Since this approach is computationally very expensive, we have developed a markedly more efficient alternative in which MD is combined with the electrodiffusion (ED) equation. In this approach, the assumptions of the ED equation can be rigorously tested, and the precision and consistency of the calculated conductance can be determined. We have demonstrated that the full current/voltage dependence and the underlying free energy profile for a simple channel can be reliably calculated from equilibrium or non-equilibrium MD simulations at a single voltage. To carry out MD simulations, a structural model of a channel has to be assumed, which is an important constraint, considering that high-resolution structures are available for only very few simple channels. If the comparison of calculated ionic conductance with electrophysiological data is satisfactory, it greatly increases our confidence that the structure and the function are described sufficiently accurately. We examined the validity of the ED for several channels embedded in phospholipid membranes - four naturally occurring channels: trichotoxin, alamethicin, p7 from hepatitis C virus (HCV) and Vpu from the HIV-1 virus, and a synthetic, hexameric channel, formed by a 21-residue peptide that contains only leucine and serine. All these channels mediate transport of potassium and chloride ions. It was found that the ED equation is satisfactory for these systems. In some of them experimental and calculated electrophysiological properties are in good agreement, whereas in others there are strong indications that the structural models are incorrect.

molecular dynamics (MD↗

NASA NTRS2017

Small Ion Channel Linking Molecular Simulations and Electrophysiology

Ion channels are pore-forming protein assemblies that mediate the transport of small ions across cell membranes. Otherwise, membrane bilayers would be almost impermeable to ions incapable to traverse the low dielectric constant, hydrophobic membrane core. Ion channels are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses, cardiac functions, muscle contraction and apoptosis. On the other extreme of biological complexity, viral ion channels (viroporins) influence many stages of the virus infection cycle either through regulating virus replication, such as entry, assembly and release or modulating the electrochemical balance in the subcellular compartments of host cells. Ion channels were crucial components of protocells. Their emergence facilitated adaptation of nascent life to different environmental conditions. The earliest ion channels must have been much simpler than most of their modern ancestors. Viral channels are among only a few naturally occurring models to study the structure, function and evolution of primordial channels. Experimental studies of these properties are difficult and often unreliable. In principle, computational methods, and molecular dynamics (MD) simulations in particular, can aid in providing information about both the structure and the function of ion channels. However, MD suffers from its own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. It is, therefore, essential to determine the reliability of MD simulations. We propose to do so on the basis of two criteria. One is channel stability on time scales that extend for several microseconds or longer. The other is the ability to reproduce the measured ionic conductance as a function of applied voltage. If both the stability and the calculated ionic conductance are satisfactory it will greatly increase our confidence that the structure and the function of a channel are described sufficiently accurately. To our knowledge, long time scale stability (approx.10 micro-sec) and the correct electrophysiology have been shown so far for only one channel - the synthetic LS3 hexamer). In this presentation, this approach will be discussed in application to two viral channels - Vpu, encoded by the HIV-1 genome and p7 of hepatitis C.

Pohorille, Andrzej↗

NASA NTRS2017

Exploring Connectivity in Sequence Space of Functional RNA

Emergence of replicable genetic molecules was one of the marking points in the origin of life, evolution of which can be conceptualized as a walk through the space of all possible sequences. A theoretical concept of fitness landscape helps to understand evolutionary processes through assigning a value of fitness to each genotype. Then, evolution of a phenotype is viewed as a series of consecutive, single-point mutations. Natural selection biases evolution toward peaks of high fitness and away from valleys of low fitness. whereas neutral drift occurs in the sequence space without direction as mutations are introduced at random. Large networks of neutral or near-neutral mutations on a fitness landscape, especially for sufficiently long genomes, are possible or even inevitable. Their detection in experiments, however, has been elusive. Although a few near-neutral evolutionary pathways have been found, recent experimental evidence indicates landscapes consist of largely isolated islands. The generality of these results, however, is not clear, as the genome length or the fraction of functional molecules in the genotypic space might have been insufficient for the emergence of large, neutral networks. Thorough investigation on the structure of the fitness landscape is essential to understand the mechanisms of evolution of early genomes. RNA molecules are commonly assumed to play the pivotal role in the origin of genetic systems. They are widely believed to be early, if not the earliest, genetic and catalytic molecules, with abundant biochemical activities as aptamers and ribozymes, i.e. RNA molecules capable, respectively, to bind small molecules or catalyze chemical reactions. Here, we present results of our recent studies on the structure of the sequence space of RNA ligase ribozymes selected through in vitro evolution. Several hundred thousands of sequences active to a different degree were obtained by way of deep sequencing. Analysis of these sequences revealed several large clusters defined such that every sequence in a cluster can be reached from any other sequence in the same cluster through a series of single point mutations. Sequences in a single cluster appear to adopt more than one secondary structure. The mechanism of refolding within a single cluster was examined. To shed light on possible evolutionary paths in the space of ribozymes, the connectivity between clusters was investigated. The effect of length of RNA molecules on the structure of the fitness landscape and possible evolutionary paths was examined by way of comparing functional sequences of 20 and 80 nucleobases in length. It was found that sequences of different lengths shared secondary structure motifs that were presumed responsible for catalytic activity, with increasing complexity and global structural rearrangements emerging in longer molecules.

Wei, Chenyu↗

NASA NTRS2017

Flexible Proteins at the Origin of Life

Almost all modern proteins possess well-defined, relatively rigid scaffolds that provide structural preorganization for desired functions. Such scaffolds require the sufficient length of a polypeptide chain and extensive evolutionary optimization. How ancestral proteins attained functionality, even though they were most likely markedly smaller than their contemporary descendants, remains a major, unresolved question in the origin of life. On the basis of evidence from experiments and computer simulations, we argue that at least some of the earliest water-soluble and membrane proteins were markedly more flexible than their modern counterparts. As an example, we consider a small, evolved in vitro ligase, based on a novel architecture that may be the archetype of primordial enzymes. The protein does not contain a hydrophobic core or conventional elements of the secondary structure characteristic of modern water-soluble proteins, but instead is built of a flexible, catalytic loop supported by a small hydrophilic core containing zinc atoms. It appears that disorder in the polypeptide chain imparts robustness to mutations in the protein core. Simple ion channels, likely the earliest membrane protein assemblies, could also be quite flexible, but still retain their functionality, again in contrast to their modern descendants. This is demonstrated in the example of antiamoebin, which can serve as a useful model of small peptides forming ancestral ion channels. Common features of the earliest, functional protein architectures discussed here include not only their flexibility, but also a low level of evolutionary optimization and heterogeneity in amino acid composition and, possibly, the type of peptide bonds in the protein backbone.

Flexible protein↗

NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗

NASA NTRS2014

Towards Co-evolution of Membrane Transport and Metabolism

Protocellular boundaries were inextricably connected to the metabolism they encapsulated: to be inheritable, early metabolism must have led to an increased rate of growth and division of vesicles and, similarly, transport through vesicle boundaries must have supported the evolution of metabolism. Even though explaining how this coupling emerged and evolved in the absence of the complex machinery of modern cells is one of the key issues in studies on the origin of life, little is known about the biochemical and biophysical processes that might have been involved. This gap in our knowledge is a major impediment in efforts to construct scenarios for the origin of life and laboratory models of protocells. A combination of experimental and computational studies carried out by us and our collaborators is aimed at helping to close this gap. Properties of membranes might have contributed to the selection of RNA as an early biopolymer. A kinetic mechanism was proposed (Sacerdote & Szostak, 2005) in which ribose was supplied more quickly than other aldopentoses to primordial cells for preferential incorporation of ribonucleotides into nucleic acids. This proposal is based on a finding that ribose permeates membranes an order of magnitude faster than its diastereomers, arabinose and xylose. Our computer simulations, which yield permeation rates in excellent agreement with experiment, and kinetic modeling explain this phenomenon in terms of inter- and intramolecular interactions involving exocyclic hydroxyl groups attached to carbon atoms of the pyranose ring (Wei and Pohorille, 2009). They also constrain scenarios for the formation of the earliest nucleic acids (Wei and Pohorille, 2013). In one scenario, sugars permeate protocellular walls and subsequently are used to synthesize nucleic acids inside protocells. As long as this process proceeds at the rate faster than 6x10(exp -3)/s, ribose derivatives will be available for synthesis easier than their diastereomers. If nucleosides or their activated derivatives are synthesized outside protocells and subsequently transported across protocellular membranes the kinetic mechanism does not apply because all diastereomers, which have their sugars in the furanose rather than pyranose form, permeate the membrane at approximately the same rate. Properties of membranes might have been also coupled to metabolism involving peptides. Recently, Adamala and Szostak (2013) have shown that a dipeptide inside fatty-acid vesicles catalyzes the formation of another dipeptide that binds to vesicle walls and, by doing so, promotes their growth at the expense of other vesicles. This coupling of metabolism, permeability of vesicles and their growth is the first demonstration of evolutionary advantage imparted by small, membrane-bound peptides. Building on this work we have calculated the rate at which different blocked amino acids are delivered to a protocell for synthesis of dipeptides. We have further shown that the dipeptides are located at the water-membrane interface rather than in the center of the bilayer. On these basis it is anticipated that other dipeptides containing aromatic, but not necessarily hydrophobic amino acids (e.g. tyrosine) could have the same catalytic effects. Insight from these studies allows for estimating the rate of vesicle growth and the rates of dipeptide synthesis required to keep the system in balance. These results, in combination with our earlier studies, lead to a general scenario for evolution from membrane-bound dipeptides to ion channels in the origin of life.

Metabolism↗

NASA NTRS2011

Origins of Protein Functions in Cells

In modern organisms proteins perform a majority of cellular functions, such as chemical catalysis, energy transduction and transport of material across cell walls. Although great strides have been made towards understanding protein evolution, a meaningful extrapolation from contemporary proteins to their earliest ancestors is virtually impossible. In an alternative approach, the origin of water-soluble proteins was probed through the synthesis and in vitro evolution of very large libraries of random amino acid sequences. In combination with computer modeling and simulations, these experiments allow us to address a number of fundamental questions about the origins of proteins. Can functionality emerge from random sequences of proteins? How did the initial repertoire of functional proteins diversify to facilitate new functions? Did this diversification proceed primarily through drawing novel functionalities from random sequences or through evolution of already existing proto-enzymes? Did protein evolution start from a pool of proteins defined by a frozen accident and other collections of proteins could start a different evolutionary pathway? Although we do not have definitive answers to these questions yet, important clues have been uncovered. In one example (Keefe and Szostak, 2001), novel ATP binding proteins were identified that appear to be unrelated in both sequence and structure to any known ATP binding proteins. One of these proteins was subsequently redesigned computationally to bind GTP through introducing several mutations that introduce targeted structural changes to the protein, improve its binding to guanine and prevent water from accessing the active center. This study facilitates further investigations of individual evolutionary steps that lead to a change of function in primordial proteins. In a second study (Seelig and Szostak, 2007), novel enzymes were generated that can join two pieces of RNA in a reaction for which no natural enzymes are known. Recently it was found that, as in the previous case, the proteins have a structure unknown among modern enzymes. In this case, in vitro evolution started from a small, non-enzymatic protein. A similar selection process initiated from a library of random polypeptides is in progress. These results not only allow for estimating the occurrence of function in random protein assemblies but also provide evidence for the possibility of alternative protein worlds. Extant proteins might simply represent a frozen accident in the world of possible proteins. Alternative collections of proteins, even with similar functions, could originate alternative evolutionary paths.

Seelig, Burchard↗