Engineering topics
Publications and source records attributed to Pohorille, Andrew.
This talk will provide an introduction to the goals of the Astrobiology program and outline some areas of interest to theoreticians and computer modelers. The contents of the talk will be I taken from the Astrobiology Roadmap, reports to the Astrobiology Institute and information, available from the Center for Computational Astrobiology web site (http://cca.arc.nasa.gov).
In the absence of extinct or extant record of protocells (the earliest ancestors of contemporary cells). the most direct way to test our understanding of the origin of cellular life is to construct laboratory models of protocells. Such efforts are currently underway in the NASA Astrobiology Program. They are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures and developing designs for molecules that perform proto-cellular functions. Many of these functions, such as import of nutrients, capture and storage of energy. and response to changes in the environment are carried out by proteins bound to membrane< We will discuss a series of large-scale, molecular-level computer simulations which demonstrate (a) how small proteins (peptides) organize themselves into ordered structures at water-membrane interfaces and insert into membranes, (b) how these peptides aggregate to form membrane-spanning structures (eg. channels), and (c) by what mechanisms such aggregates perform essential proto-cellular functions, such as proton transport of protons across cell walls, a key step in cellular bioenergetics. The simulations were performed using the molecular dynamics method, in which Newton's equations of motion for each item in the system are solved iteratively. The problems of interest required simulations on multi-nanosecond time scales, which corresponded to 10(exp 6)-10(exp 8) time steps.
On the path from inanimate to animate matter, a key step was the self-organization of molecules into protocells - the earliest ancestors of contemporary cells. Studies of the properties of protocells and the mechanisms by which they maintained themselves and reproduced are an important part of astrobiology. These studies also have the potential to greatly impact research in nanotechnology and computer science. Previous studies of protocells have focussed on self-replication. In these systems, Darwinian evolution occurs through a series of small alterations to functional molecules whose identities are stored. Protocells, however, may have been incapable of such storage. We hypothesize that under such conditions, the replication of functions and their interrelationships, rather than the precise identities of the functional molecules, is sufficient for survival and evolution. This process is called non-genomic evolution. Recent breakthroughs in experimental protein chemistry have opened the gates for experimental tests of non-genomic evolution. On the basis of these achievements, we have developed a stochastic model for examining the evolutionary potential of non-genomic systems. In this model, the formation and destruction (hydrolysis) of bonds joining amino acids in proteins occur through catalyzed, albeit possibly inefficient, pathways. Each protein can act as a substrate for polymerization or hydrolysis, or as a catalyst of these chemical reactions. When a protein is hydrolyzed to form two new proteins, or two proteins are joined into a single protein, the catalytic abilities of the product proteins are related to the catalytic abilities of the reactants. We will demonstrate that the catalytic capabilities of such a system can increase. Its evolutionary potential is dependent upon the competition between the formation of bond-forming and bond-cutting catalysts. The degree to which hydrolysis preferentially affects bonds in less efficient, and therefore less well-ordered, peptides is also critical to evolution of a non-genomic system. Based on these results, a new computational object called a "molnet" is defined. Like a neural network, it is formed of interconnected units that send "signals" to each other. Like molecules, neural networks have a specific function once their structure is defined. The difference between a molnet and traditional neural networks, is that input to molnets is not simply passed along and processed from input to output units, but rather it is utilized to form and break connections(bonds), and thus to form new structures. Molnets represent a powerful tool that can be used to understand the conditions under which chemical systems can form large molecules, such as proteins, and display ever more complex functions. This has direct applications, for example to the design of smart,synthetic fabrics. Additional information is contained in the original.
One of the main goals of astrobiology is to understand the origin of cellular life. The most direct approach to this problem 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 that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing 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 self-organize into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides form membrane-spanning structures (e.g. channels)? (3) By what mechanisms do such structures perform their functions? The simulations are performed using the molecular dynamics 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 in real space while long-range forces are evaluated in reciprocal space, using a particle-mesh algorithm which is of order O(NInN). With a time step of 2 femtoseconds, problems occurring on multi-nanosecond time scales (10(exp 6)-10(exp 8) time steps) are accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude, which requires algorithmic improvements and codes scalable to a large number of processors. Work in this direction is in progress. Two series of simulations are discussed. In one series, it is shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and fold into helical structures (see Figure). Once in the membrane, the peptides exhibit orientational flexibility with changing conditions, which may have provided a mechanism of transmitting signals between the protocell and its environment. 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 identify the mechanisms by which protons move through the gate.
We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.
Extraterrestrially delivered organics in the origin of cellular life. Various processes leading to the emergence of cellular life from organics delivered from space to earth or other planetary bodies in the solar system will be reviewed. The focus will be on: (1) self-assembly of amphiphilic material to vesicles and other structures, such as micelles and multilayers, and its role in creating environments suitable for chemical catalysis, (2) a possible role of extraterrestrial delivery of organics in the formation of the simplest bioenergetics (3) mechanisms leading from amino acids or their precursors to simple peptides and, subsequently, to the evolution of metabolism. These issues will be discussed from two opposite points of view: (1) Which molecules could have been particularly useful in the protobiological evolution; this may provide focus for searching for these molecules in interstellar media. (2) Assuming that a considerable part of the inventory of organic matter on the early earth was delivered extraterrestrially, what does relative abundance of different organics in space tell us about the scenario leading to the origin of life.
Two of Ernest Overton's lasting contributions to biology are the Meyer-Overton relationship between the potency of an anesthetic and its solubility in oil, and the Overton rule which relates the permeability of a membrane to the oil-water partition coefficient of the permeating molecule. A growing body of experimental evidence, however, cannot be reconciled with these theories. In particular, the molecular nature of membranes, unknown to Overton, needs to be included in any description of these phenomena. Computer simulations are ideally suited for providing atomic-level information about the behavior of small molecules in membranes. The authors discuss simulation studies relevant to Overton's ideas. Through simulations it was found that anesthetics tend to concentrate at interfaces and their anesthetic potency correlates better with solubility at the water-membrane interface than with solubility in oil. Simulation studies of membrane permeation revealed the anisotropic nature of the membranes, as evidenced, for example, by the highly nonuniform distribution of free volume in the bilayer. This, in turn, influences the diffusion rates of solutes, which increase with the depth in the membrane. Small solutes tend to move by hopping between voids in the bilayer, and this hopping motion may be responsible for the deviation from the Overton rule of the permeation rates of these molecules.
In the absence of extinct or extant 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: First, a very large population of candidate molecules is generated using a random synthetic approach. Next, the small numbers of molecules that can accomplish the desired task are selected. These 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. To date, we have obtained "a proof of concept" by evolving simple, novel proteins capable of selectively binding adenosine tri-phosphate (ATP). Our next goal is to create an enzyme that can phosphorylate amino acids and another to catalyze the formation of peptide bonds in the absence of nucleic acid templates. This latter reaction does not take place in contemporary cells. once developed, these enzymes will be encapsulated in liposomes so that they will function in a simulated cellular environment. To provide a continuous energy supply, usually needed to activate the substrates, an energy transduction complex which generates ATP from adenosine diphosphate, inorganic phosphate and light will be used. This system, consisting of two modern proteins, ATP synthase and bacteriorhodopsin, has already been built and shown to work efficiently. By coupling chemical synthesis to such a system, it will be possible to drive chemical reactions by light if only the substrates for these reactions are supplied.
Folding of the capped LQQLLQQLLQL peptide is investigated at the water-hexane interface by molecular dynamics simulations over 161.5 nanoseconds. Initially placed in the aqueous phase as a beta-strand, the peptide rapidly adsorbs to the interface, where it adopts an amphipathic conformation. The marginal presence of non-amphipathic structures throughout the complete trajectory indicate- that the corresponding conformations are strongly disfavored at the interface. It is further suggestive that folding in an interfacial environment proceeds through a pathway of successive amphipathic intermediates. The energetic and entropic penalties involved in the conformational changes along this pathway markedly increase the folding time-scales of LQQLLQQLLQL, explaining why the alpha-helix, the hypothesized lowest free energy structure for a sequence with a hydrophobic periodicity of 3.6, has not been reached yet. The formation of a type I beta-turn at the end of the simulation confirms the importance of such motifs as initiation sites allowing the peptide to coalesce towards a secondary structure.
A model for the evolution of biological systems in the absence of a nucleic acid-like genome is proposed and applied to model the earliest living organisms -- protocells composed of membrane encapsulated peptides. Assuming that the peptides can make and break bonds between amino acids, and bonds in non-functional peptides are more likely to be destroyed than in functional peptides, it is demonstrated that the catalytic capabilities of the system as a whole can increase. This increase is defined to be non-genomic evolution. The relationship between the proposed mechanism for evolution and recent experiments on self-replicating peptides is discussed.
Our objective is to help explain how the earliest ancestors of contemporary cells (protocells) performed their essential functions employing only the molecules available in the protobiological milieu. Our hypothesis is that vesicles, built of amphiphilic, membrane-forming materials, emerged early in protobiological evolution and served as precursors to protocells. We further assume that the cellular functions associated with contemporary membranes, such as capturing and, transducing of energy, signaling, or sequestering organic molecules and ions, evolved in these membrane environments. An alternative hypothesis is that these functions evolved in different environments and were incorporated into membrane-bound structures at some later stage of evolution. We focus on the application of the fundamental principles of physics and chemistry to determine how they apply to the formation of a primitive, functional cell. Rather than attempting to develop specific models for cellular functions and to identify the origin of the molecules which perform these functions, our goal is to define the structural and energetic conditions that any successful model must fulfill, therefore providing physico-chemical boundaries for these models. We do this by carrying out large-scale, molecular level computer simulations on systems of interest.
The 25 amino acids long, transmembrane fragment of the Influenza virus M2 protein forms a homotetrameric channel that transports protons across lipid bilayers. It has been postulated that high efficiency and selectivity of this process is due to gating by four histidine residues that occlude the channel lumen in the closed state. Two mechanisms of gating have been postulated. In one mechanism, the proton is "shuttled" through the gate by attaching to the delta nitrogen atom on the extracellular side of the imidazole ring, followed by the release of the proton attached to the epsilon nitrogen atom on the opposite side. In the second mechanism, the four histidines move away from each other due to electrostatic repulsion upon protonation, thus opening the gate sufficiently that a wire of water molecules can penetrate the gate. Then, protons are transported by "hopping" along the wire. In this paper, both mechanisms are evaluated in a series of molecular dynamics simulations by investigating stability of different protonation states of the channel that are involved in these mechanisms. For the shuttle mechanism, these are states with all epsilon protonated histidines, one biprotonated residue or one histidine protonated in the delta position. For the gate opening mechanism, this is the state in which all four histidines are biprotonated. In addition, a state with two biprotonated histidines is considered. For each system, composed of the protein channel embedded in phospholipid bilayer located between two water lamellae, a molecular dynamics trajectory of approximately 1.3 ns (after equilibration) was obtained. It is found that the states involved in the shuttle mechanism are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfers involved in the shuttle. In contract, the fully biprotonated state, implicated in the gate opening mechanism, is not stable and the channel looses its structural integrity. If only two histidines are biprotonated the channel deforms but remains intact with the gate mostly closed. In summary, the results of this study lend support to the shuttle mechanism but not to the gate opening mechanism of proton gating in M2.
The thermal conductivity was found to be only weakly dependent on the microstructure of the amorphous ice. In general, the amorphous ices were found to have thermal conductivities of the same order of magnitude as liquid water. This is in contradiction to recent experimental estimates of the thermal conductivity of amorphous ice, and it is suggested that the extremely low value obtained experimentally is due to larger-scale defects in the ice, such as cracks, but is not an intrinsic property of the bulk amorphous ice.
The main objective of this research was to estimate the thermal conductivity of cometry ices from computer simulations of model amorphous ices. This was divided into four specific tasks: (1) Generating samples of amorphous water ices at different microporosities; (2) Comparing the resulting molecular structures of the ices with experimental results, for those densities where data was available; (3) Calculating the thermal conductivities of liquid water and bulk amorphous ices and comparing these results with experimentally determined thermal conductivities; and (4) Investigating how the thermal conductivity of amorphous ice depends upon the microscopic porosity of the samples. The thermal conductivity was found to be only weakly dependent on the microstructure of the amorphous ice. In general, the amorphous ices were found to have thermal conductivities of the same order of magnitude as liquid water. This is in contradiction to recent experimental estimates of the thermal conductivity of amorphous ice, and it is suggested that the extremely low value obtained experimentally is due to larger-scale defects in the ice, such as cracks, but it is not an intrinsic property of the bulk amorphous ice.
The use of transmembrane proton gradients in energy transduction is an almost universal feature of life on earth. These proton gradients are established and maintained by specialized assemblies of proteins which actively pump protons across membranes. One broad class of proton pumps uses captured light energy to drive the proton pumping. Our goal is to elucidate the minimum structural requirements of a light-driven proton-pump. There are two basic components to a simple light-driven proton pump: a source of photo-generated protons and a "gate-keeper" which prevents these protons from reattaching themselves to their source. A wide variety of molecules in the membrane, even as simple as polycyclic aromatic hydrocarbons, are capable of releasing protons when illuminated. Our work is therefore focused on the design of the "gate-keeper." Our initial model involves a pair of proton acceptors, coupled to each other by a transient water bridge, and supported in the membrane by a small bundle of peptide helices. Upon illumination, the proton source transfers its proton to the:- first acceptor of the gate-keeper. While the reverse reaction is highly probable, all that is needed to ensure irreversibility is a nonvanishing probability that the proton will be transferred to the second acceptor across a transient water bridge. Back transfer of the proton to the first acceptor, and thence to the proton source, is impeded by the free energy required to move the proton uphill towards the. proton source and by the disruption of the transient water bridge. As a prototypical water-bridged proton transfer system, we are studying the transfer of a proton across a water bridge from a formic acid to a formate anion. With a pK(sub alpha), of 3.7. formic acid is a good model for the acidic amino acids glutamate and aspartate which are good candidates for gate-keeper proton acceptors. Simulations of proton transfer reactions in a membrane are complicated by the quantum mechanical nature of the process breaking and forming chemical bonds.
Even though most monomeric peptides are disordered in water they can adopt sequence-dependent, ordered structures, such as a-helices, at aqueous interfaces. This property is relevant to cellular signaling, membrane fusion, and the action of toxins and antibiotics. The mechanism of folding nonpolar peptides at the water-hexane interface was studied in the example of an 11-mer, of poly-L-leucine. Initially placed as a random coil on the water side of the interface, the peptide folded into an a-helix in 36 ns. Simultaneously, the peptide translocated into the hexane side of the interface. Folding was not sequential and involved a 3/10-helix as an intermediate. The folded peptide was either parallel to the interface or had its C-terminus exposed to water. An 11-mer, LQQLLQQLLQL, composed of leucine (L) and glutamine (G), was taken as a model amphiphilic peptide. It rapidly adopted an amphiphilic, disordered structure at the interface. Further folding proceeded through a series of amphiphilic intermediates.
The ability of small peptides to organize at aqueous interfaces was examined by performing a series of large-scale, molecular dynamics computer simulations of several peptides composed of two amino acids, nonpolar leucine (L) and polar glutamine (Q). The peptides differed in size and sequence of the amino acids. Studies on dipeptides LL, LQ, QL and QQ were extended to two heptamers, LQQLLQL and LQLQLQL, designed to maximize interfacial stability of an alpha-helix and a beta-strand, respectively, by exposing polar side chains to water and nonpolar side chains to a nonpolar phase. Finally, a transition of an undecamer, composed entirely of leucine residues, from a disordered structure in water to an alpha-helix in a nonpolar phase representing the interior of the membrane was investigated. Complete folding of a peptide in solution was accomplished for the first time in computer simulations. The simulations revealed several basic principles governing the sequence-dependent organization of peptides at interfaces. Short peptides tend to accumulate at interfaces and acquire ordered structures, providing that they have a proper sequence of polar and nonpolar amino acids. The dominant factor determining the interfacial structure of peptides is the hydrophobic effect, which is manifested at aqueous interfaces as a tendency for polar and nonpolar groups of the solute to segregate into the aqueous and nonpolar phases, respectively. If peptides consist of nonpolar residue's only, they become inserted into the nonpolar phase. As demonstrated by the example of the leucine undecamer, such peptides fold into an alpha-helix as they partition into the nonpolar medium. The folding proceeds through an intermediate, called 3-10-helix, which remains in equilibrium with the alpha-helix. Once in the nonpolar environment, the peptides can readily change their orientation with respect to the interface from parallel to perpendicular, especially in response to local electric fields. The ability of nonpolar peptides to modify both the structure and orientation with respect to the interface from parallel to perpendicular, especially in response to local electric fields. The ability of nonpolar peptides to modify both the structure and orientation with changing external conditions may have provided a simple mechanism of transmitting signals from the environment to the interior of a cell.
Some forms of nanotechnology appear to have enormous potential to improve aerospace and computer systems; computational nanotechnology, the design and simulation of programmable molecular machines, is crucial to progress. NASA Ames Research Center has begun a computational nanotechnology program including in-house work, external research grants, and grants of supercomputer time. Four goals have been established: (1) Simulate a hypothetical programmable molecular machine replicating itself and building other products. (2) Develop molecular manufacturing CAD (computer aided design) software and use it to design molecular manufacturing systems and products of aerospace interest, including computer components. (3) Characterize nanotechnologically accessible materials of aerospace interest. Such materials may have excellent strength and thermal properties. (4) Collaborate with experimentalists. Current in-house activities include: (1) Development of NanoDesign, software to design and simulate a nanotechnology based on functionalized fullerenes. Early work focuses on gears. (2) A design for high density atomically precise memory. (3) Design of nanotechnology systems based on biology. (4) Characterization of diamonoid mechanosynthetic pathways. (5) Studies of the laplacian of the electronic charge density to understand molecular structure and reactivity. (6) Studies of entropic effects during self-assembly. Characterization of properties of matter for clusters up to sizes exhibiting bulk properties. In addition, the NAS (NASA Advanced Supercomputing) supercomputer division sponsored a workshop on computational molecular nanotechnology on March 4-5, 1996 held at NASA Ames Research Center. Finally, collaborations with Bill Goddard at CalTech, Ralph Merkle at Xerox Parc, Don Brenner at NCSU (North Carolina State University), Tom McKendree at Hughes, and Todd Wipke at UCSC are underway.