Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular modelling”

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.

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

At least 793 records · Page 44

Delta excitations and shell-model information in heavy-ion, charge-exchange reactions

We calculate total cross sections for coherent pion production using localized plane-wave approximations for the shell-structure of valence nucleons that are excited to delta particles in the intermediate state in the (12C, 12B) and (12C, 12N) charge-exchange, heavy-ion reactions. We find comparable agreement to projectile downshift data for 12C(12C, 12B)12N. Then we improve the formalism by replacing the localized plane wave bound states with harmonic oscillator states which are imbedded in a multipole expansion approach and calculate pion differential cross sections to test for the sensitivity of the spectra to the single-particle mass parameter.

NASA Program Space Medicine↗

Evolutionary anomalies among the aminoacyl-tRNA synthetases

Unexpected relationships among the various aminoacyl-tRNA synthetases continue to be uncovered. The question arises - is this mainly the result of promiscuous exchange, or is the confusion really a reflection of the differential loss of past duplications? Phylogenetic analysis may yet provide the answer.

Non-NASA Center↗

Molecular Data for a Biochemical Model of DNA Radiation Damage: Electron Impact Ionization and Dissociative Ionization of DNA Bases and Sugar-Phosphate Backbone

As part of the database for building up a biochemical model of DNA radiation damage, electron impact ionization cross sections of sugar-phosphate backbone and DNA bases have been calculated using the improved binary-encounter dipole (iBED) model. It is found that the total ionization cross sections of C3'- and C5'-deoxyribose-phospate, two conformers of the sugar-phosphate backbone, are close to each other. Furthermore, the sum of the ionization cross sections of the separate deoxyribose and phosphate fragments is in close agreement with the C3'- and C5'-deoxyribose-phospate cross sections, differing by less than 10%. Of the four DNA bases, the ionization cross section of guanine is the largest, then in decreasing order, adenine, thymine, and cytosine. The order is in accordance with the known propensity of oxidation of the bases by ionizing radiation. Dissociative ionization (DI), a process that both ionizes and dissociates a molecule, is investigated for cytosine. The DI cross section for the formation of H and (cytosine-Hl)(+), with the cytosine ion losing H at the 1 position, is also reported. The threshold of this process is calculated to be 17.1 eV. Detailed analysis of ionization products such as in DI is important to trace the sequential steps in the biochemical process of DNA damage.

Dateo, Christopher E.↗

Energetics of the lattice: packing elements in crystals of four-stranded intercalated cytosine-rich DNA molecules

Condensation of single molecules from solution into crystals represents a transition between distinct energetic states. In solution, the atomic interactions within the molecule dominate. In the crystalline state, however, a set of additional interactions are formed between molecules in close contact in the lattice--these are the packing interactions. The crystal structures of d(CCCT), d(TAACCC), d(CCCAAT), and d(AACCCC) have in common a four-stranded intercalated cytosine segment, built by stacked layers of cytosine.cytosine+ (C.C+) base pairs coming from two parallel duplexes that intercalate into each other with opposite polarity. The intercalated cytosine segments in these structures are similar in their geometry, even though the sequences crystallized in different space groups. In the crystals, adenine and thymine residues of the sequences are used to build the three-dimensional crystal lattice by elaborately interacting with symmetry-related molecules. The packing elements observed provide novel insight about the copious ways in which nucleic acid molecules can interact with each other--for example, when folded in more complicated higher order structures, such as mRNA and chromatin.

Non-NASA Center↗

Structure of an anti-HIV-1 hammerhead ribozyme complex with a 17-mer DNA substrate analog of HIV-1 gag RNA and a mechanism for the cleavage reaction: 750 MHz NMR and computer experiments

The structure of an anti-HIV-1 ribozyme-DNA abortive substrate complex was investigated by 750 MHz NMR and computer modeling experiments. The ribozyme was a chimeric molecule with 30 residues-18 DNA nucleotides, and 12 RNA residues in the conserved core. The DNA substrate analog had 17 residues. The chimeric ribozyme and the DNA substrate formed a shortened ribozyme-abortive substrate complex of 47 nucleotides with two DNA stems (stems I and III) and a loop consisting of the conserved core residues. Circular dichroism spectra showed that the DNA stems assume A-family conformation at the NMR concentration and a temperature of 15 degrees C, contrary to the conventional wisdom that DNA duplexes in aqueous solution populate entirely in the B-form. It is proposed that the A-family RNA residues at the core expand the A-family initiated at the core into the DNA stems because of the large free energy requirement for the formation of A/B junctions. Assignments of the base H8/H6 protons and H1' of the 47 residues were made by a NOESY walk. In addition to the methyl groups of all T's, the imino resonances of stems I and III and AH2's were assigned from appropriate NOESY walks. The extracted NMR data along with available crystallographic data, were used to derive a structural model of the complex. Stems I and III of the final model displayed a remarkable similarity to the A form of DNA; in stem III, a GC base pair was found to be moving into the floor of the minor groove defined by flanking AT pairs; data suggest the formation of a buckled rhombic structure with the adjacent pair; in addition, the base pair at the interface of stem III and the loop region displayed deformed geometry. The loop with the catalytic core, and the immediate region of the stems displayed conformational multiplicity within the NMR time scale. A catalytic mechanism for ribozyme action based on the derived structure, and consistent with biochemical data in the literature, is proposed. The complex between the anti HIV-1 gag ribozyme and its abortive DNA substrate manifests in the detection of a continuous track of A.T base pairs; this suggests that the interaction between the ribozyme and its DNA substrate is stronger than the one observed in the case of the free ribozyme where the bases in stem I and stem III regions interact strongly with the ribozyme core region (Sarma, R. H., et al. FEBS Letters 375, 317-23, 1995). The complex formation provides certain guidelines in the design of suitable therapeutic ribozymes. If the residues in the ribozyme stem regions interact with the conserved core, it may either prevent or interfere with the formation of a catalytically active tertiary structure.

NASA Discipline Exobiology↗

Entropy and charge in molecular evolution--the case of phosphate

Biopoesis, the creation of life, implies molecular evolution from simple components, randomly distributed and in a dilute state, to form highly organized, concentrated systems capable of metabolism, replication and mutation. This chain of events must involve environmental processes that can locally lower entropy in several steps; by specific selection from an indiscriminate mixture, by concentration from dilute solution, and in the case of the mineral-induced processes, by particular effectiveness in ordering and selective reaction, directed toward formation of functional biomolecules. Numerous circumstances provide support for the notion that negatively charged molecules were functionally required and geochemically available for biopoesis. Sulfite ion may have been important in bisulfite complex formation with simple aldehydes, facilitating the initial concentration by sorption of aldehydes in positively charged surface active minerals. Borate ion may have played a similar, albeit less investigated role in forming charged sugar complexes. Among anionic species, oligophosphate ions and charged phosphate esters are likely to have been of even more wide ranging importance, reflected in the continued need for phosphate in a proposed RNA world, and extending its central role to evolved biochemistry. Phosphorylation is shown to result in selective concentration by surface sorption of compounds, otherwise too dilute to support condensation reactions. It provides protection against rapid hydrolysis of sugars and, by selective concentration, induces the oligomerization of aldehydes. As a manifestation of life arisen, phosphate already appears in an organic context in the oldest preserved sedimentary record.

NASA Discipline Exobiology↗

Modeling of a possible conformational change associated with the catalytic mechanism in the hammerhead ribozyme

Here we describe a possible model of the cleavage mechanism in the hammerhead ribozyme. In this model, the 2' hydroxyl of C17 is moved into an appropriate orientation for an in-line attack on the G1.1 phosphate through a change in its sugar pucker from C3' endo to C2' endo. This conformational change in the active site is caused by a change in the uridine turn placing the N2 and N3 atoms of G5 of the conserved core in hydrogen bonding geometry with the N3 and N2 atoms on the conserved G16.2 residue. The observed conformational change in the uridine turn suggests an explanation for the conservation of G5. In the crystal structure of H.M. Pley et al., Nature 372, 68-74 (1994), G5 is situated 5.3A away from G16.2. However, the uridine turn is sufficiently flexible to allow this conformational change with relatively modest changes in the backbone torsion angles (average change of 14.2 degrees). Two magnesium ions were modeled into the active site with positions analogous to those described in the functionally similar Klenow fragment 3'-5' exonuclease (L.S. Beese and T.A. Steitz, EMBO J. 10, 25-33 (1991)), the Group I intron (T.A. Steitz and J.A. Steitz, P.N.A.S. U.S.A. 90, 6498-6502 (1993); R.F. Setlik et al., J. Biomol. Str. Dyn. 10, 945-972 (1993)) and other phosphotransferases. Comparison of this model with one in which the uridine turn conformation was not changed showed that although the changes in the C17 sugar pucker could be modeled, insufficient space existed for the magnesium ions in the active site.

Non-NASA Center↗

Stable loop in the crystal structure of the intercalated four-stranded cytosine-rich metazoan telomere

In most metazoans, the telomeric cytosine-rich strand repeating sequence is d(TAACCC). The crystal structure of this sequence was solved to 1.9-A resolution. Four strands associate via the cytosine-containing parts to form a four-stranded intercalated structure held together by C.C+ hydrogen bonds. The base-paired strands are parallel to each other, and the two duplexes are intercalated into each other in opposite orientations. One TAA end forms a highly stabilized loop with the 5' thymine Hoogsteen-base-paired to the third adenine. The 5' end of this loop is in close proximity to the 3' end of one of the other intercalated cytosine strands. Instead of being entirely in a DNA duplex, this structure suggests the possibility of an alternative conformation for the cytosine-rich telomere strands.

NASA Discipline Exobiology↗

Alpha-helical hydrophobic polypeptides form proton-selective channels in lipid bilayers

Proton translocation is important in membrane-mediated processes such as ATP-dependent proton pumps, ATP synthesis, bacteriorhodopsin, and cytochrome oxidase function. The fundamental mechanism, however, is poorly understood. To test the theoretical possibility that bundles of hydrophobic alpha-helices could provide a low energy pathway for ion translocation through the lipid bilayer, polyamino acids were incorporated into extruded liposomes and planar lipid membranes, and proton translocation was measured. Liposomes with incorporated long-chain poly-L-alanine or poly-L-leucine were found to have proton permeability coefficients 5 to 7 times greater than control liposomes, whereas short-chain polyamino acids had relatively little effect. Potassium permeability was not increased markedly by any of the polyamino acids tested. Analytical thin layer chromatography measurements of lipid content and a fluorescamine assay for amino acids showed that there were approximately 135 polyleucine or 65 polyalanine molecules associated with each liposome. Fourier transform infrared spectroscopy indicated that a major fraction of the long-chain hydrophobic peptides existed in an alpha-helical conformation. Single-channel recording in both 0.1 N HCl and 0.1 M KCl was also used to determine whether proton-conducting channels formed in planar lipid membranes (phosphatidylcholine/phosphatidylethanolamine, 1:1). Poly-L-leucine and poly-L-alanine in HCl caused a 10- to 30-fold increase in frequency of conductive events compared to that seen in KCl or by the other polyamino acids in either solution. This finding correlates well with the liposome observations in which these two polyamino acids caused the largest increase in membrane proton permeability but had little effect on potassium permeability. Poly-L-leucine was considerably more conductive than poly-L-alanine due primarily to larger event amplitudes and, to a lesser extent, a higher event frequency. Poly-L-leucine caused two populations of conductive events, one in the 0.1-0.5 pA range, and one in the 1.0-5.0 pA range, whereas nearly all events caused by poly-L-alanine were in the 0.1-0.5 pA range at an applied voltage of +60 mV. The channel-like activity appeared to switch between conductive and nonconductive states, with most open-times in the range of 50-200 ms. We conclude that hydrophobic polyamino acids produce proton-conducting defects in lipid bilayers that may be used to model functional proton channels in biological membranes.

NASA Discipline Exobiology↗

Conformational influence of the ribose 2'-hydroxyl group: crystal structures of DNA-RNA chimeric duplexes

We have crystallized three double-helical DNA-RNA chimeric duplexes and determined their structures by X-ray crystallography at resolutions between 2 and 2.25 A. The two self-complementary duplexes [r(G)d(CGTATACGC)]2 and [d(GCGT)r(A)d(TACGC)]2, as well as the Okazaki fragment d(GGGTATACGC).r(GCG)d(TATACCC), were found to adopt A-type conformations. The crystal structures are non-isomorphous, and the crystallographic environments for the three chimeras are different. A number of intramolecular interactions of the ribose 2'-hydroxyl groups contribute to the stabilization of the A-conformation. Hydrogen bonds between 2'-hydroxyls and 5'-oxygens or phosphate oxygens, in addition to the previously observed hydrogen bonds to 1'-oxygens of adjacent riboses and deoxyriboses, are observed in the DNA-RNA chimeric duplexes. The crystalline chimeric duplexes do not show a transition between the DNA A- and B-conformations. CD spectra suggest that the Okazaki fragment assumes an A-conformation in solution as well. In this molecule the three RNA residues may therefore lock the complete decamer in the A-conformation. Crystals of an all-DNA strand with the same sequence as the self-complementary chimeras show a morphology which is different from those of the chimera crystals. Moreover, the oligonucleotide does not match any of the sequence characteristics of DNAs usually adopting the A-conformation in the crystalline state (e.g., octamers with short alternating stretches of purines and pyrimidines). In DNA-RNA chimeric duplexes, it is therefore possible that a single RNA residue can drive the conformational equilibrium toward the A-conformation.

NASA Discipline Number 52-20↗

A residue substitution near the beta-ionone ring of the retinal affects the M substates of bacteriorhodopsin

The switch in the bacteriorhodopsin photocycle, which reorients access of the retinal Schiff base from the extracellular to the cytoplasmic side, was suggested to be an M1----M2 reaction (Varo and Lanyi. 1991. Biochemistry. 30:5008-5015, 5016-5022). Thus, in this light-driven proton pump it is the interconversion of proposed M substates that gives direction to the transport. We find that in monomeric, although not purple membrane-lattice immobilized, D115N bacteriorhodopsin, the absorption maximum of M changes during the photocycle: in the time domain between its rise and decay it shifts 15 nm to the blue relative to the spectrum at earlier times. This large shift strongly supports the existence of two M substates. Since D115 is located near the beta-ionone ring of the retinal, the result raises questions about the possible involvement of the retinal chain or protein residues as far away as 10 A from the Schiff base in the mechanism of the switching reaction.

NASA Discipline Exobiology↗

Structural basis for stabilization of Z-DNA by cobalt hexaammine and magnesium cations

In the equilibrium between B-DNA and Z-DNA in poly(dC-dG), the [Co(NH3)6]3+ ion stabilizes the Z form 4 orders of magnitude more effectively than the Mg2+ ion. The structural basis of this difference is revealed in Z-DNA crystal structures of d(CpGpCpGpCpG) stabilized by either Na+/Mg2+ or Na+/Mg2+ plus [Co(NH3)6]3+. The crystals diffract X-rays to high resolution, and the structures were refined at 1.25 A. The [Co(NH3)6]3+ ion forms five hydrogen bonds onto the surface of Z-DNA, bonding to a guanine O6 and N7 as well as to a phosphate group in the ZII conformation. The Mg2+ ion binds through its hydration shell with up to three hydrogen bonds to guanine N7 and O6. Higher charge, specific fitting of more hydrogen bonds, and a more stable complex all contribute to the great effectiveness of [Co(NH3)6]3+ in stabilizing Z-DNA.

NASA Discipline Exobiology↗

Structural and energetic consequences of mutations in a solvated hydrophobic cavity

The structural and energetic consequences of modifications to the hydrophobic cavity of interleukin 1-beta (IL-1beta) are described. Previous reports demonstrated that the entirely hydrophobic cavity of IL-1beta contains positionally disordered water. To gain a better understanding of the nature of this cavity and the water therein, a number of mutant proteins were constructed by site-directed mutagenesis, designed to result in altered hydrophobicity of the cavity. These mutations involve the replacement of specific phenylalanine residues, which circumscribe the cavity, with tyrosine, tryptophan, leucine and isoleucine. Using differential scanning calorimetry to determine the relative stabilities of the wild-type and mutant proteins, we found all of the mutants to be destabilizing. X-ray crystallography was used to identify the structural consequences of the mutations. No clear correlation between the hydrophobicities of the specific side-chains introduced and the resulting stabilities was found.

Solvents/chemistry↗

The force exerted by the membrane potential during protein import into the mitochondrial matrix

The force exerted on a targeting sequence by the electrical potential across the inner mitochondrial membrane is calculated on the basis of continuum electrostatics. The force is found to vary from 3.0 pN to 2.2 pN (per unit elementary charge) as the radius of the inner membrane pore (assumed aqueous) is varied from 6.5 to 12 A, its measured range. In the present model, the decrease in force with increasing pore width arises from the shielding effect of water. Since the pore is not very much wider than the distance between water molecules, the full shielding effect of water may not be present; the extreme case of a purely membranous pore without water gives a force of 3.2 pN per unit charge, which should represent an upper limit. When applied to mitochondrial import experiments on the protein barnase, these results imply that forces between 11 +/- 2 pN and 13.5 +/- 2.5 pN catalyze the unfolding of barnase in those experiments. A comparison of these results with unfolding forces measured using atomic force microscopy is made.

Models, Molecular↗

Prediction of Mechanical Properties of Polymers With Various Force Fields

The effect of force field type on the predicted elastic properties of a polyimide is examined using a multiscale modeling technique. Molecular Dynamics simulations are used to predict the atomic structure and elastic properties of the polymer by subjecting a representative volume element of the material to bulk and shear finite deformations. The elastic properties of the polyimide are determined using three force fields: AMBER, OPLS-AA, and MM3. The predicted values of Young s modulus and shear modulus of the polyimide are compared with experimental values. The results indicate that the mechanical properties of the polyimide predicted with the OPLS-AA force field most closely matched those from experiment. The results also indicate that while the complexity of the force field does not have a significant effect on the accuracy of predicted properties, small differences in the force constants and the functional form of individual terms in the force fields determine the accuracy of the force field in predicting the elastic properties of the polyimide.

Odegard, Gregory M.↗