Engineering PapersSearch

Engineering topics

New, Michael H.

Publications and source records attributed to New, Michael H..

Effects of artificial defoliation of pines on the structure and physiology of the soil fungal community of a mixed pine-spruce forest

Loss of photosynthetic area can affect soil microbial communities by altering the availability of fixed carbon. We used denaturing gradient gel electrophoresis (DGGE) and Biolog filamentous-fungus plates to determine the effects of artificial defoliation of pines in a mixed pine-spruce forest on the composition of the fungal community in a forest soil. As measured by DGGE, two fungal species were affected significantly by the defoliation of pines (P < 0.001); the frequency of members of the ectomycorrhizal fungus genus Cenococcum decreased significantly, while the frequency of organisms of an unidentified soil fungus increased. The decrease in the amount of Cenococcum organisms may have occurred because of the formation of extensive hyphal networks by species of this genus, which require more of the carbon fixed by their host, or because this fungus is dependent upon quantitative differences in spruce root exudates. The defoliation of pines did not affect the overall composition of the soil fungal community or fungal-species richness (number of species per core). Biolog filamentous-fungus plate assays indicated a significant increase (P < 0.001) in the number of carbon substrates utilized by the soil fungi and the rate at which these substrates were used, which could indicate an increase in fungal-species richness. Thus, either small changes in the soil fungal community give rise to significant increases in physiological capabilities or PCR bias limits the reliability of the DGGE results. These data indicate that combined genetic and physiological assessments of the soil fungal community are needed to accurately assess the effect of disturbance on indigenous microbial systems.

Pinus/growth & development/physiology

Effects of litter addition on ectomycorrhizal associates of a lodgepole pine (Pinus contorta) stand in Yellowstone National Park

Increasing soil nutrients through litter manipulation, pollution, or fertilization can adversely affect ectomycorrhizal (EM) communities by inhibiting fungal growth. In this study, we used molecular genetic methods to determine the effects of litter addition on the EM community of a Pinus contorta stand in Yellowstone National Park that regenerated after a stand-replacing fire. Two controls were used; in unmodified control plots nothing was added to the soil, and in perlite plots perlite, a chemically neutral substance, was added to maintain soil moisture and temperature at levels similar to those under litter. We found that (i) species richness did not change significantly following perlite addition (2.6 +/- 0.3 species/core in control plots, compared with 2.3 +/- 0.3 species/core in perlite plots) but decreased significantly (P < 0.05) following litter addition (1.8 +/- 0.3 species/core); (ii) EM infection was not affected by the addition of perlite but increased significantly (P < 0.001) in response to litter addition, and the increase occurred only in the upper soil layer, directly adjacent to the added litter; and (iii) Suillus granulatus, Wilcoxina mikolae, and agaricoid DD were the dominant organisms in controls, but the levels of W. mikolae and agaricoid DD decreased significantly in response to both perlite and litter addition. The relative levels of S. granulatus and a fourth fungus, Cortinariaceae species 2, increased significantly (P < 0.01 and P < 0.05, respectively) following litter addition. Thus, litter addition resulted in some negative effects that may be attributable to moisture-temperature relationships rather than to the increased nutrients associated with litter. Some species respond positively to litter addition, indicating that there are differences in their physiologies. Hence, changes in the EM community induced by litter accumulation also may affect ecosystem function.

Fungi/classification/growth & development

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

Protein Kinase Classification with 2866 Hidden Markov Models and One Support Vector Machine

The main application considered in this paper is predicting true kinases from randomly permuted kinases that share the same length and amino acid distributions as the true kinases. Numerous methods already exist for this classification task, such as HMMs, motif-matchers, and sequence comparison algorithms. We build on some of these efforts by creating a vector from the output of thousands of structurally based HMMs, created offline with Pfam-A seed alignments using SAM-T99, which then must be combined into an overall classification for the protein. Then we use a Support Vector Machine for classifying this large ensemble Pfam-Vector, with a polynomial and chisquared kernel. In particular, the chi-squared kernel SVM performs better than the HMMs and better than the BLAST pairwise comparisons, when predicting true from false kinases in some respects, but no one algorithm is best for all purposes or in all instances so we consider the particular strengths and weaknesses of each.

Weber, Ryan

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.

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

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.

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

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.

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

Evolutionary Cell Computing: From Protocells to Self-Organized Computing

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.

Colombano, Silvano

Molecular Simulations in Astrobiology

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.

Pohorille, Andrew

Computer Simulations of Small Molecules in Membranes: Insights from Computer Simulations into the Interactions of Small Molecules with Lipid Bilayers

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.

Pohorille, Andrew

An Inherited Efficiencies Model of Non-Genomic Evolution

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.

New, Michael H.

Membrane-Based Functions in the Origin of Cellular Life

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.

Chipot, Christophe

Computer-aided design of a proton pump

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.

New, Michael H.