Engineering topics

Pohorille, Andrew

Publications and source records attributed to Pohorille, Andrew.

At least 55 records · Page 3

NASA NTRS2002

The M2 Proton Channel of Influenza Virus: How Does It Work?

The transport of protons across membranes is an essential process for both bioenergetics of modem cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATR), synthesized from adenosine diphosphate. ATR, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this study was: how the same process can be accomplished with the aid of similar, but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC), which is a good model of the biological membranes focusing cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M2 protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M2 protein is 97 amino acids in length, but a fragment 25 amino acids long, which contains a transmembrane domain of 19 amino acids flanked by 3 amino acids on each side, is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules, and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional, important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel, in which the histidine residues were in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel looses its structural integrity and one helix moves away from the remaining three. This result indicates that such a mechanism of proton transport is unlikely. The simulations revealed that translocation along a network of water molecules in the channel and tautomerization of the histidine residues in the M2 proteins in the most likely mechanism of proton transport. The results not only explain how a remarkably simple protein system can efficiently aid in the formation of proton gradients across cell walls, but also suggest how this system can be genetically re-engineered to become a directional, reversible proton pump. Such a pump can provide energy to laboratory-built models of simple cellular systems. If they were successfully constructed it would greatly advance our understanding of the beginnings of life and find important applications in medicine and pharmacology.

Pohorille, Andrew↗

NASA NTRS2002

Artificial cells: prospects for biotechnology

A variety of techniques can now be used to alter the genome of a cell. Although these techniques are very powerful, they have limitations related to cost and efficiency of scale. Artificial cells designed for specific applications combine properties of biological systems such as nanoscale efficiency, self-organization and adaptability at relatively low cost. Individual components needed for such structures have already been developed, and now the main challenge is to integrate them in functional microscopic compartments. It will then become possible to design and construct communities of artificial cells that can perform different tasks related to therapeutic and diagnostic applications.

Review↗

NASA NTRS2002

Who or What? Self-Replication and Function-Reproduction in the Origin of Life

In this presentation, we will present results on the fundamental properties of two classes of replicating systems: autocatalytic replicators that reproduce exact copies of a template molecule, and function reproducers that maintain a set of essential functions without replicating the identities of the functional moieties. We will describe the stability and behavior in-the-large of autocatalytic replicators. Most importantly, we have found no sharp distinction between an autocatalytic and a non-autocatalytic domain. We will also present a new derivation of von Kiedrowski's square-root rate law. Function - reproducers are proposed as an important component of protocells and we will present theoretical results on a simple model system that incorporates known peptide biophysics. For a wide range of parameters, we have shown that this type of system can improve its overall performance, even in the absence of any method for information storage. This type of system improvement is defined to be non-genomic evolution.

New, Michael H.↗

NASA NTRS2002

Research in Computational Astrobiology

We present results from several projects in the new 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. We have developed a procedure for calculating long-range effects in molecular dynamics using a plane wave expansion of the electrostatic potential. This method is expected to be highly efficient for simulating biological systems on massively parallel supercomputers. We have perform genomics analysis on a family of actin binding proteins. We have performed quantum mechanical calculations on carbon nanotubes and nucleic acids, which simulations will allow us to investigate possible sources of organic material on the early earth. Finally, we have developed a model of protobiological chemistry using neural networks.

Chaban, Galina↗

NASA NTRS2002

Protein Solvation in Membranes and at Water-Membrane Interfaces

Different salvation properties of water and membranes mediate a host of biologically important processes, such as folding, insertion into a lipid bilayer, associations and functions of membrane proteins. These processes will be discussed in several examples involving synthetic and natural peptides. In particular, a mechanism by which a helical peptide becomes inserted into a model membrane will be described. Further, the molecular mechanism of recognition and association of protein helical segments in membranes will be discussed. These processes are crucial for proper functioning of a cell. A membrane-spanning domain of glycophorin A, which exists as a helical dimer, serves as the model system. For this system, the free energy of dissociation of the helices is being determined for both the wild type and a mutant, in which dimerization is disrupted.

Pohorille, Andrew↗

NASA NTRS2002

Peptides at Membrane Surfaces and their Role in the Origin of Life

All ancestors of contemporary cells (protocells) had to transport ions and organic matter across membranous walls, capture and utilize energy and transduce environmental signals. In modern organisms, all these functions are preformed by membrane proteins. We make the parsimonious assumption that in the protobiological milieu the same functions were carried out by their simple analogs - peptides. This, however, required that simple peptides could self-organize into ordered, functional structures. In a series of detailed, molecular-level computer simulations we demonstrated how this is possible. One example is the peptide (LSLLLSL)3 which forms a trameric bundle capable of transporting protons across membranes. Another example is 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 in this channel. The channel can be re-engineered into a simple proton pump.

Pohorille, Andrew↗

NASA NTRS2001

Artificial Cells: Prospects for Biotechnology

A variety of techniques can now be used to alter the genome of a cell. Although these techniques are very powerful, they also have limitations related to cost and efficiency of scale. Artificial cells designed for specific applications combine properties of biological systems such as nano-scale efficiency, self-organization and adaptability at relatively low cost. Individual components needed for such structures have already been developed, and now the main challenge is to integrate them in functional microscopic compartments. It will then become possible to design and construct communities of artificial cells that can perform different tasks related to therapeutic and diagnostic applications.

Pohorille, Andrew↗

NASA NTRS2001

Calculating Free Energies Using Average Force

A new, general formula that connects the derivatives of the free energy along the selected, generalized coordinates of the system with the instantaneous force acting on these coordinates is derived. The instantaneous force is defined as the force acting on the coordinate of interest so that when it is subtracted from the equations of motion the acceleration along this coordinate is zero. The formula applies to simulations in which the selected coordinates are either unconstrained or constrained to fixed values. It is shown that in the latter case the formula reduces to the expression previously derived by den Otter and Briels. If simulations are carried out without constraining the coordinates of interest, the formula leads to a new method for calculating the free energy changes along these coordinates. This method is tested in two examples - rotation around the C-C bond of 1,2-dichloroethane immersed in water and transfer of fluoromethane across the water-hexane interface. The calculated free energies are compared with those obtained by two commonly used methods. One of them relies on determining the probability density function of finding the system at different values of the selected coordinate and the other requires calculating the average force at discrete locations along this coordinate in a series of constrained simulations. The free energies calculated by these three methods are in excellent agreement. The relative advantages of each method are discussed.

Darve, Eric↗

NASA NTRS2001

Models of Protocellular Structure, Function and Evolution

In the absence of any record of protocells, the most direct way to test our understanding of the origin of cellular life is to construct laboratory models that capture important features of protocellular systems. Such efforts are currently underway in a collaborative project between NASA-Ames, Harvard Medical School and University of California. They are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures. The centerpiece of this project is a method for the in vitro evolution of protein enzymes toward arbitrary catalytic targets. A similar approach has already been developed for nucleic acids in which a small number of functional molecules are selected from a large, random population of candidates. The selected molecules are next vastly multiplied using the polymerase chain reaction. A mutagenic approach, in which the sequences of selected molecules are randomly altered, can yield further improvements in performance or alterations of specificities. Unfortunately, the catalytic potential of nucleic acids is rather limited. Proteins are more catalytically capable but cannot be directly amplified. In the new technique, this problem is circumvented by covalently linking each protein of the initial, diverse, pool to the RNA sequence that codes for it. Then, selection is performed on the proteins, but the nucleic acids are replicated. Additional information is contained in the original extended abstract.

New, Michael H.↗

NASA NTRS2001

Molecular Dynamics Studies of Structure and Functions of Water-Membrane Interfaces

A large number of essential cellular processes occur at the interfaces between water and membranes. The selectivity and dynamics of these processes are largely determined by the structural and electrical properties of the water-membrane interface. We investigate these properties by the molecular dynamics method. Over the time scales of the simulations, the membrane undergoes fluctuations described by the capillary wave model. These fluctuations produce occasional thinning defects in the membrane which provide effective pathways for passive transport of ions and small molecules across the membrane. Ions moving through the membrane markedly disrupt its structure and allow for significant water penetration into the membrane interior. Selectivity of transport, with respect to ionic charge, is determined by the interfacial electrostatic potential. Many small molecules. of potential significance in catalysis, bioenergetics and pharmacology, are shown to bind to the interface. The energetics and dynamics of this process will be discussed.

Pohorille, Andrew↗

NASA NTRS2001

Better, Cheaper, Faster Molecular Dynamics

Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.

Pohorille, Andrew↗

NASA NTRS2001

Molecular Simulations in Astrobiology

One of the main goals of astrobiology is to understand the origin of cellular life. In the absence of any record of the earliest ancestors of contemporary cells, protocells, the most direct way to test our understanding of their characteristics is to construct laboratory models of protocells. Such efforts, currently underway in the NASA Astrobiology Program, are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures and developing designs of molecules that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing and storing energy, and responding to changes in the environment, are carried out by proteins bound to membranes. We use computer simulations to address the following, questions about these proteins: (1) How do small proteins (peptides) organize themselves into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides aggregate to form membrane-spannin(y structures (e.g., channels)? (3) By what mechanisms do such aggregates perform their functions? The simulations are performed using the molecular dynamics (MD) method. In this method, Newton's equations of motion for each atom in the system are solved iteratively. At each time step, the forces exerted on each atom by the remaining atoms are evaluated by dividing them into two parts. Short-range forces are calculated directly in real space while long-range forces are evaluated in reciprocal space, usually using a particle-mesh algorithm which is of order O(NlnN). Currently, a time step of 2 femtoseconds is typically used, thereby making studies of problems occurring on multi-nanosecond time scales (10(exp 6) - 10(exp 8) time steps) accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude. Such an extension requires both algorithmic improvements and codes scalable to a large number of parallel processors. Work in this direction is in progress. Two specific series of simulations that demonstrate how peptides self-organize and function in membranes are discussed. In one series of simulations, it was shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and, simultaneously, fold into two different helical structures, which remain in equilibrium. Once in the membrane, the peptides can readily change their orientation, especially in response to local electric fields. This structural and orientational flexibility of peptides with changing conditions may have provided a mechanism of transmitting signals between the environment and the interior of the protocell. In another series of simulations, the mechanism by which a simple protein channel efficiently mediates proton transport across membranes was investigated. This process is a key step in cellular bioenergetics. In the channel under study, proton transport is gated by four histidines that occlude the channel pore. The simulations demonstrate that protons move through the gate by a "shuttle" mechanism, wherein one histidine is protonated on the extracellular side and, subsequently, the proton bound on the opposite side is released.

Pohorille, Andrew↗

NASA NTRS2001

Proton Pumps: Mechanism of Action and Applications

Recent progress in understanding molecular structures and mechanisms of action of proton pumps has paved the way to their novel applications in biotechnology. Proton pumps, in particular bacteriorhodopsin and ATP synthases, are capable of continuous, renewable conversion of light to chemical, mechanical or electrical energy, which can be used in macro- or nano-scale devices. The capability of protein systems incorporated into liposomes to generate ATP, which can be further used to drive chemical reactions, and to act as molecular motors has been already demonstrated. Other possible applications of such biochemical devices include targeted drug delivery and biocatalytic re actors. All these devices might prove superior to their inorganic alternatives.

Lanyi, Janos K.↗

NASA NTRS2001

Simulation of Peptides at Aqueous Interfaces

Behavior of peptides at water-membrane interfaces is of great interest in studies on cellular transport and signaling, membrane fusion, and the action of toxins and antibiotics. Many peptides, which exist in water only as random coils, can form sequence-dependent, ordered structures at aqueous interfaces, incorporate into membranes and self-assembly into functional units, such as simple ion channels. Multi -nanosecond molecular dynamics simulations have been carried out to study the mechanism and energetics of interfacial folding of both non-polar and amphiphilic peptides, their insertion into membranes and association into higher-order structures. The simulations indicate that peptides fold non-sequentially, often through a series of amphiphilic intermediates. They further incorporate into the membrane in a preferred direction as folded monomers, and only then aggregate into dimers and, possibly, further into "dimers of dimers".

Pohorille, Andrew↗

NASA NTRS2001

Models of Protocellular Structure, Function and Evolution

In the absence of any record of protocells, the most direct way to test our understanding, of the origin of cellular life is to construct laboratory models that capture important features of protocellular systems. Such efforts are currently underway in a collaborative project between NASA-Ames, Harvard Medical School and University of California. They are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures. The centerpiece of this project is a method for the in vitro evolution of protein enzymes toward arbitrary catalytic targets. A similar approach has already been developed for nucleic acids in which a small number of functional molecules are selected from a large, random population of candidates. The selected molecules are next vastly multiplied using the polymerase chain reaction.

New, Michael H.↗

NASA NTRS2001

Proton Transport

The transport of protons across membranes is an essential process for both bioenergetics of modern cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATP), synthesized from adenosine diphosphate. ATP, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this Study was: how the same process can be accomplished with the aid of similar but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC) which is a good model of the biological membranes forming cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M$_2$ protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M$_2$ protein is 97 amino acids in length, but a fragment 25 amino acids long. which contains a transmembrane domain of 19 amino acids flanked by three amino acids on each side. is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules. and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small ions, such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel with the histidine residues in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel loses its structural integrity and one helix moves away from the remaining three. Additional information is contained within the original extended abstract.

Pohorille, Andrew↗

NASA NTRS2000

Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm

One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous mixing. A similar difficulty may be encountered in any method relying on global modifications of phase space.

Darve, Eric↗

NASA NTRS2000

Models of protocellular structures, functions and evolution

The central step in the origin of life was the emergence of organized structures from organic molecules available on the early earth. These predecessors to modern cells, called 'proto-cells,' were simple, membrane bounded structures able to maintain themselves, grow, divide, and evolve. Since there is no fossil record of these earliest of life forms, it is a scientific challenge to discover plausible mechanisms for how these entities formed and functioned. To meet this challenge, it is essential to create laboratory models of protocells that capture the main attributes associated with living systems, while remaining consistent with known, or inferred, protobiological conditions. This report provides an overview of a project which has focused on protocellular metabolism and the coupling of metabolism to energy transduction. We have assumed that the emergence of systems endowed with genomes and capable of Darwinian evolution was preceded by a pre-genomic phase, in which protocells functioned and evolved using mostly proteins, without self-replicating nucleic acids such as RNA.

Pohorille, Andrew↗