SEARCH · Engineering Papers

Results for “Molecular simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1998

Mars Pathfider Rarefied Aerodynamics: Computations and Measurements

On July 4, 1997, after traveling close to 500 million km, the Pathfinder spacecraft successfully completed entry, descent, and landing at Mars. In the present paper, the focus is on the hypersonic rarefied portion of Pathfinder's atmospheric entry where the synergy of flight measurements, aerodynamic calculations, and atmospheric modeling tools are used to extract Pathfinder's attitude and the freestream density. Accuracy of the capsule aerodynamics directly impacts the inferred atmospheric properties extracted from deceleration measurements made by on-board accelerometers. The range of rarefaction considered in this study extends from the free molecular to continuum conditions and angles of attack from O to 30 deg. The aerodynamic computations are made with free-molecular and direct simulation Monte Carlo codes. The calculations show that Pathfinder is statically unstable for much of the transitional rarefied regime. Due to the relatively modest forces and the gyroscopic motion of the spacecraft, the angle of attack excursions were less than 5 deg as inferred from force measurements for the rarefied portion of entry and approached a nominal zero degree trim angle near hypersonic continuum conditions.

Moss, James N.↗

NASA NTRS1999

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↗

NASA NTRS2001

Thermodynamic Behavior of Nano-sized Gold Clusters on the (001) Surface

We have studied thermal expansion of the surface layers of the hexagonally reconstructed Au (001) surface using a classical Molecular Dynamics (MD) simulation technique with an Embedded Atomic Method (EAM) type many-body potential. We find that the top-most hexagonal layer contracts as temperature increases, whereas the second layer expands or contracts depending on the system size. The magnitude of expansion coefficient of the top layer is much larger than that of the other layers. The calculated thermal expansion coefficients of the top-most layer are about -4.93 x 10(exp -5)Angstroms/Kelvin for the (262 x 227)Angstrom cluster and -3.05 x 10(exp -5)Angstroms/Kelvin for (101 x 87)Angstrom cluster. The Fast Fourier Transform (FFT) image of the atomic density shows that there exists a rotated domain of the top-most hexagonal cluster with rotation angle close to 1 degree at temperature T less than 1000Kelvin. As the temperature increases this domain undergoes a surface orientational phase transition. These predictions are in good agreement with previous phenomenological theories and experimental studies.

Paik, Sun M.↗

NASA NTRS2001

Thermal Expansion and Diffusion Coefficients of Carbon Nanotube-Polymer Composites

Classical molecular dynamics (MD) simulations employing Brenner potential for intra-nanotube interactions and van der Waals forces for polymer-nanotube interface have been used to investigate thermal expansion and diffusion characteristics of carbon nanotube-polyethylene composites. Addition of carbon nanotubes to polymer matrix is found to significantly increase the glass transition temperature Tg, and thermal expansion and diffusion coefficients in the composite above Tg. The increase has been attributed to the temperature dependent increase of the excluded volume for the polymer chains, and the findings could have implications in the composite processing, coating and painting applications.

Wei, Chengyu↗

NASA NTRS1997

Peptide Folding and Translocation Across the Water-Membrane Interface

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.

Pohorille, Andrew↗

NASA NTRS2002

Analysis of Carbon Nanotube Pull-out from a Polymer Matrix

Molecular dynamics (MD) simulations of carbon nanotube (NT) pull-out from a polymer matrix are carried out. As the NT pull-out develops in the simulation, variations in the displacement and velocities of the NT are monitored. The existence of a carbon-ring-based period in NT sliding during pull-out is identified. Linear trends in the NT velocity-force relation are observed and used to estimate an effective viscosity coefficient for interfacial sliding at the NT/polymer interface. As a result, the entire process of NT pull-out is characterized by an interfacial friction model that is based on a critical pull-out force, and an analog of Newton's friction law used to describe the NT/polymer interfacial sliding.

Frankland, S. J. V.↗

NASA NTRS2002

Peptides at Membrane Surfaces and their Role in Prebiotic Evolution

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

Pohorille, Andrew↗

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 NTRS2003

Membrane-Based Functions in the Origin of Cellular Life

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

Wilson, Michael A.↗

NASA NTRS2003

Structural Basis for "Flip-Flop" Action of Human Pyruvate Dehydrogenase

The derivative of vitamin B1, thiamin pyrophosphate is a cofactor of pyruvate dehydrogenase, a component enzyme of the mitochondrial pyruvate dehydrogenase multienzyme complex that plays a major role in directing energy metabolism in the cell. This cofactor is used to cleave the C(sup alpha)-C(=O) bond of pyruvate followed by reductive acetyl transfer to lipoyl-dihydrolipoamide acetyltransferase. In alpha(sub 2)beta(sub 2)-tetrameric human pyruvate dehydrogenase, there are two cofactor binding sites, each of them being a center of independently conducted, although highly coordinated enzymatic reactions. The dynamic nonequivalence of two, otherwise chemically equivalent, catalytic sites can now be understood based on the recently determined crystal structure of the holo-form of human pyruvate dehydrogenase at 1.95A resolution. The structure of pyruvate dehydrogenase was determined using a combination of MAD phasing and molecular replacement followed by rounds of torsion-angles molecular-dynamics simulated-annealing refinement. The final pyruvate dehydrogenase structure included coordinates for all protein amino acids two cofactor molecules, two magnesium and two potassium ions, and 742 water molecules. The structure was refined to R = 0.202 and R(sub free) = 0.244. Our structural analysis of the enzyme folding and domain assembly identified a simple mechanism of this protein motion required for the conduct of catalytic action.

Ciszak, Ewa↗

NASA NTRS2003

Theory of Transport of Long Polymer Molecules through Carbon Nanotube Channels

A theory of transport of long chain polymer molecules through carbon nanotube (CNT) channels is developed using Fokker-Planck equation and direct molecular dynamics (MD) simulations. The mean transport or translocation time tau is found to depend on the chemical potential energy, entropy and diffusion coefficient. A power law dependence tau approx. N(sup 2)is found where N is number of monomers in a molecule. For 10(exp 5)-unit long polyethylene molecules, tau is estimated to be approx. 1micro-s. The diffusion coefficient of long polymer molecules inside CNTs, like that of short ones, are found to be few orders of magnitude larger than in ordinary silicate based zeolite systems.

Wei, Chenyu↗

NASA NTRS2003

A Dynamic Pathway for Stone-Wales Bond Rotation on Carbon Nanotubes through Diamond-Like Bonds

A new lower energy barrier with a two-step pathway of Stone-Wales (SW) ,ond rotation on carbon nanotubes (CNTs) is found through molecular dynamics (MD) simulations of CNTs under tension. The first step involves going over to a stable sp3-like metastable configuration with half rotated and partially tilted C-C bond. The second step involves going over to the fully rotated C-C bond with the formation of a SW defect in the nanotube. The energy barrier for this two-step dynamic pathway is significantly lower than the previously known static barrier for in-plane rotation of the C-C bond on a tensile strained (> 4%) CNT.

Wei, Chen-Yu↗

NASA NTRS2003

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

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

Pohorille, Andrew↗

NASA NTRS2003

Membrane peptides and their role in protobiological evolution

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

Models, Chemical↗

NASA NTRS2004

Modeling the interaction Between Ethylene Diamine and Water Films on the Surface of a Carbon Nanotube

It has been observed that a carbon nanotube (CNT) AFM tip coated with ethylene diamine (EDA) penetrates the liquid water-air interface more easily than an uncoated nanotube tip. The EDA coating remains intact through repeated cycles of dipping and removal. In order to understand the physical basis for this observation, we use ab initio quantum chemistry calculations to study the EDA-CNT-water interaction and to parameterize a force field describing this system. Molecular dynamics (MD) simulations are carried out for EDA-water mixtures and an EDA-coated carbon nanotube immmed in water. These simulations are similar to our earlier MD study that characterized the CNT-water interface. The attractive CNT-EDA and CNT-water interactions arise primarily from van der Waals forces, and the EDA-EDA, EDA-water and water-water interactions are mainly due to hydrogen bond formation. The binding energ of single EDA molecule to the nanotube is nearly three times larger than the corresponding value found for water (4.3 versus 1.5 kcal mol, respectively). The EDA molecules readily stick to and diffuse along the CNT surface. As a resulf mixing of the EDA and water films does not occur on the timescale of the MD simulations. The EDA film reduces the hydrophobicity of the nanotube surface and acts like a prototypical surfactant in stabilizing the suspension of carbon nanotubes in water. For this presentation, we use the MD simulations to determine how the presence of the carbon nanotube surface perturbs the properties of EDA-water mixtures.

Jaffe, Richard L.↗

NASA NTRS2006

Cassini-Huygens Aerodynamics with Comparison to Flight

An analysis has been performed on the aerodynamics of the Cassini spacecraft as it passed through the atmosphere of Titan with and without the Huygens probe attached. The free stream density of Titan s atmosphere was measured by two methods. However, these methods resulted in very different values of density, one result being 3-5 times higher than the other. In an attempt to understand the source of this discrepancy and verify the assumptions made by and the original engineering analysis performed by the Jet Propulsion Laboratory (JPL), free molecular and direct simulation Monte Carlo (DSMC) analyses were performed for two atmospheric passes. The drag coefficient was calculated using an area based on a Monte Carlo area determination program. Although the source of the discrepancy has not been determined, it has been confirmed that the original, simplified force-and-moment analysis performed by JPL produced results that were comparable to the high-fidelity DSMC analysis and that the source of the discrepancy lies elsewhere.

Liechty, Derek S.↗

NASA NTRS2007

Interactions of the C-terminal Domain of Human Ku70 with DNA Substrate: A Molecular Dynamics Study

NASA is developing a systems biology approach to improve the assessment of health risks associated with space radiation. The primary toxic and mutagenic lesion following radiation exposure is the DNA double strand break (DSB), thus a model incorporating proteins and pathways important in response and repair of this lesion is critical. One key protein heterodimer for systems models of radiation effects is the Ku(sub 70/80) complex. The Ku70/80 complex is important in the initial binding of DSB ends following DNA damage, and is a component of nonhomologous end joining repair, the primary pathway for DSB repair in mammalian cells. The C-terminal domain of Ku70 (Ku70c, residues 559-609), contains an helix-extended strand-helix motif and similar motifs have been found in other nucleic acid-binding proteins critical for DNA repair. However, the exact mechanism of damage recognition and substrate specificity for the Ku heterodimer remains unclear in part due to the absence of a high-resolution structure of the Ku70c/DNA complex. We performed a series of molecular dynamics (MD) simulations on a system with the subunit Ku70c and a 14 base pairs DNA duplex, whose starting structures are designed to be variable so as to mimic their different binding modes. By analyzing conformational changes and energetic properties of the complex during MD simulations, we found that interactions are preferred at DNA ends, and within the major groove, which is consistent with previous experimental investigations. In addition, the results indicate that cooperation of Ku70c with other subunits of Ku(sub 70/80) is necessary to explain the high affinity of binding as observed in experiments.

Hu, Shaowen↗

NASA NTRS2007

Interactions of the SAP Domain of Human Ku70 with DNA Substrate: A Molecular Dynamics Study

NASA is developing a systems biology approach to improve the assessment of health risks associated with space radiation. The primary toxic and mutagenic lesion following radiation exposure is the DNA double strand break (DSB), thus a model incorporating proteins and pathways important in response and repair of this lesion is critical. One key protein heterodimer for systems models of radiation effects is the Ku70/80 complex. The Ku70/80 complex is important in the initial binding of DSB ends following DNA damage, and is a component of nonhomologous end joining repair, the primary pathway for DSB repair in mammalian cells. The SAP domain of Ku70 (residues 556-609), contains an a helix-extended strand-helix motif and similar motifs have been found in other nucleic acid-binding proteins critical for DNA repair. However, the exact mechanism of damage recognition and substrate specificity for the Ku heterodimer remains unclear in part due to the absence of a high-resolution structure of the SAP/DNA complex. We performed a series of molecular dynamics (MD) simulations on a system with the SAP domain of Ku70 and a 10 base pairs DNA duplex. Large-scale conformational changes were observed and some putative binding modes were suggested based on energetic analysis. These modes are consistent with previous experimental investigations. In addition, the results indicate that cooperation of SAP with other domains of Ku70/80 is necessary to explain the high affinity of binding as observed in experiments.

Hu, Shaowen↗