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 289 records · Page 16

An algorithm for converting a virtual-bond chain into a complete polypeptide backbone chain

A systematic analysis is presented of the algorithm for converting a virtual-bond chain, defined by the coordinates of the alpha-carbons of a given protein, into a complete polypeptide backbone. An alternative algorithm, based upon the same set of geometric parameters used in the Purisima-Scheraga algorithm but with a different "linkage map" of the algorithmic procedures, is proposed. The global virtual-bond chain geometric constraints are more easily separable from the loal peptide geometric and energetic constraints derived from, for example, the Ramachandran criterion, within the framework of this approach.

NASA Discipline Exobiology↗

Effects of carbonyl bond, metal cluster dissociation, and evaporation rates on predictions of nanotube production in high-pressure carbon monoxide

The high-pressure carbon monoxide (HiPco) process for producing single-wall carbon nanotubes (SWNTs) uses iron pentacarbonyl as the source of iron for catalyzing the Boudouard reaction. Attempts using nickel tetracarbonyl led to no production of SWNTs. This paper discusses simulations at a constant condition of 1300 K and 30 atm in which the chemical rate equations are solved for different reaction schemes. A lumped cluster model is developed to limit the number of species in the models, yet it includes fairly large clusters. Reaction rate coefficients in these schemes are based on bond energies of iron and nickel species and on estimates of chemical rates for formation of SWNTs. SWNT growth is measured by the conformation of CO2. It is shown that the production of CO2 is significantly greater for FeCO because of its lower bond energy as compared with that of NiCO. It is also shown that the dissociation and evaporation rates of atoms from small metal clusters have a significant effect on CO2 production. A high rate of evaporation leads to a smaller number of metal clusters available to catalyze the Boudouard reaction. This suggests that if CO reacts with metal clusters and removes atoms from them by forming MeCO, this has the effect of enhancing the evaporation rate and reducing SWNT production. The study also investigates some other reactions in the model that have a less dramatic influence.

Models, Chemical↗

Mineral induced phosphorylation of glycolate ion--a metaphor in chemical evolution

Bilateral surface-active minerals with excess positive charge concentrate glycolate and trimetaphosphate ion from l0(-3) m aqueous solution to half-saturation of the internal surface sites, and induce phosphorylation of glycolate ion in the mineral with trimetaphosphate, sorbed from l0(-2) m solution. By utilizing reactants from dilute solution at near-neutral pH, and eliminating the need for participating organic nitrogen compounds, the reaction comprises several elements considered necessary for geochemical realism in models for molecular evolution.

NASA Discipline Exobiology↗

Structural models of the MscL gating mechanism

Three-dimensional structural models of the mechanosensitive channel of large conductance, MscL, from the bacteria Mycobacterium tuberculosis and Escherichia coli were developed for closed, intermediate, and open conformations. The modeling began with the crystal structure of M. tuberculosis MscL, a homopentamer with two transmembrane alpha-helices, M1 and M2, per subunit. The first 12 N-terminal residues, not resolved in the crystal structure, were modeled as an amphipathic alpha-helix, called S1. A bundle of five parallel S1 helices are postulated to form a cytoplasmic gate. As membrane tension induces expansion, the tilts of M1 and M2 are postulated to increase as they move away from the axis of the pore. Substantial expansion is postulated to occur before the increased stress in the S1 to M1 linkers pulls the S1 bundle apart. During the opening transition, the S1 helices and C-terminus amphipathic alpha-helices, S3, are postulated to dock parallel to the membrane surface on the perimeter of the complex. The proposed gating mechanism reveals critical spatial relationships between the expandable transmembrane barrel formed by M1 and M2, the gate formed by S1 helices, and "strings" that link S1s to M1s. These models are consistent with numerous experimental results and modeling criteria.

NASA Discipline Cell Biology↗

A novel reagentless sensing system for measuring glucose based on the galactose/glucose-binding protein

The galactose/glucose-binding protein (GBP) is synthesized in the cytoplasm of Escherichia coli in a precursor form and exported into the periplasmic space upon cleavage of a 23-amino-acid leader sequence. GBP binds galactose and glucose in a highly specific manner. The ligand induces a hinge motion in GBP and the resultant protein conformational change constitutes the basis of the sensing system. The mglB gene, which codes for GBP, was isolated from the chromosome of E. coli using the polymerase chain reaction (PCR). Since wild-type GBP lacks cysteines in its structure, introducing this amino acid by site-directed mutagenesis ensures single-label attachment at specific sites with a sulfhydro-specific fluorescent probe. Site-directed mutagenesis by overlap extension PCR was performed to prepare three different mutants to introduce a single cysteine residue at positions 148, 152, and 182. Since these residues are not involved in ligand binding and since they are located at the edge of the binding cleft, they experience a significant change in environment upon binding of galactose or glucose. The sensing system strategy is based on the fluorescence changes of the probe as the protein undergoes a structural change on binding. In this work a reagentless sensing system has been rationally designed that can detect submicromolar concentrations of glucose. The calibration plots have a linear working range of three orders of magnitude. Although the system can sense galactose as well, this epimer is not a potential interfering substance since its concentration in blood is negligible. Copyright 2001 Academic Press.

Non-NASA Center↗

Autoinhibition of a calmodulin-dependent calcium pump involves a structure in the stalk that connects the transmembrane domain to the ATPase catalytic domain

The regulation of Ca(2+)-pumps is important for controlling [Ca(2+)] in the cytosol and organelles of all eukaryotes. Here, we report a genetic strategy to identify residues that function in autoinhibition of a novel calmodulin-activated Ca(2+)-pump with an N-terminal regulatory domain (isoform ACA2 from Arabidopsis). Mutant pumps with constitutive activity were identified by complementation of a yeast (K616) deficient in two Ca(2+)-pumps. Fifteen mutations were found that disrupted a segment of the N-terminal autoinhibitor located between Lys(23) and Arg(54). Three mutations (E167K, D219N, and E341K) were found associated with the stalk that connects the ATPase catalytic domain (head) and with the transmembrane domain. Enzyme assays indicated that the stalk mutations resulted in calmodulin-independent activity, with V(max), K(mATP), and K(mCa(2+)) similar to that of a pump in which the N-terminal autoinhibitor had been deleted. A highly conservative substitution at Asp(219) (D219E) still produced a deregulated pump, indicating that the autoinhibitory structure in the stalk is highly sensitive to perturbation. In plasma membrane H(+)-ATPases from yeast and plants, similarly positioned mutations resulted in hyperactive pumps. Together, these results suggest that a structural feature of the stalk is of general importance in regulating diverse P-type ATPases.

NASA Discipline Plant Biology↗

Calcium-binding proteins and development

The known roles for calcium-binding proteins in developmental signaling pathways are reviewed. Current information on the calcium-binding characteristics of three classes of cell-surface developmental signaling proteins (EGF-domain proteins, cadherins and integrins) is presented together with an overview of the intracellular pathways downstream of these surface receptors. The developmental roles delineated to date for the universal intracellular calcium sensor, calmodulin, and its targets, and for calcium-binding regulators of the cytoskeleton are also reviewed.

Review↗

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↗