Search NASA⌕ Search

SEARCH · Search NASA

Results for “binding sites”

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 199 records · Page 11

Glycogen synthase kinase 3 phosphorylates kinesin light chains and negatively regulates kinesin-based motility

Membrane-bounded organelles (MBOs) are delivered to different domains in neurons by fast axonal transport. The importance of kinesin for fast antero grade transport is well established, but mechanisms for regulating kinesin-based motility are largely unknown. In this report, we provide biochemical and in vivo evidence that kinesin light chains (KLCs) interact with and are in vivo substrates for glycogen synthase kinase 3 (GSK3). Active GSK3 inhibited anterograde, but not retrograde, transport in squid axoplasm and reduced the amount of kinesin bound to MBOs. Kinesin microtubule binding and microtubule-stimulated ATPase activities were unaffected by GSK3 phosphorylation of KLCs. Active GSK3 was also localized preferentially to regions known to be sites of membrane delivery. These data suggest that GSK3 can regulate fast anterograde axonal transport and targeting of cargos to specific subcellular domains in neurons.

NASA Discipline Neuroscience↗

Synthesis of adriamycin-coupled polyglutaraldehyde microspheres and evaluation of their cytostatic activity

Adriamycin was coupled to polyglutaraldehyde microspheres having an average diameter of 4500 A. The coupled microspheres remained stable during incubation with cells. Full cytostatic activity was observed when the coupled adriamycin was tested with murine or human leukemia and murine sarcoma cell lines. A 10-fold increase in sensitivity was obtained with drug-resistant human leukemia cell lines. Repeated use of the coupled microspheres in the cytostatic assays did not decrease their activity, indicating that these complexes can be recycled. The results suggest that coupled adriamycin sufficiently perturbs the plasma membrane to lead to cytostatic activity. It is proposed that this mode of drug delivery provides multiple and repetitious sites for drug-cell interactions. In addition, the drug-polymer complexes may overcome those forms of resistance that are the result of decreased drug binding at the cell surface.

Tokes, Z. A.↗

Specialization of the DNA-Cleaving Activity of a Group I Ribozyme Through In Vitro Evolution

In an earlier study, an in vitro evolution procedure was applied to a large population of variants of the Tetrahymena group 1 ribozyme to obtain individuals with a 10(exp 5)-fold improved ability to cleave a target single-stranded DNA substrate under simulated physiological conditions. The evolved ribozymes also showed a twofold improvement, compared to the wild-type, in their ability to cleave a single-stranded RNA substrate. Here, we report continuation of the in vitro evolution process using a new selection strategy to achieve both enhanced DNA and diminished RNA-cleavage activity. Our strategy combines a positive selection for DNA cleavage with a negative selection against RNA binding. After 36 "generations" of in vitro evolution, the evolved population showed an approx. 100-fold increase in the ratio of DNA to RNA-cleavage activity. Site-directed mutagenesis experiment confirmed the selective advantage of two covarying mutations within the catalytic core of ribozyme that are largely responsible for this modified behavior. The population of ribozymes has now undergone a total of 63 successive generations of evolution, resulting in an average 28 mutations relative to the wild-type that are responsible for the altered phenotype.

Tsang, Joyce↗

Purification method for recombinant proteins based on a fusion between the target protein and the C-terminus of calmodulin

Calmodulin (CaM) was used as an affinity tail to facilitate the purification of the green fluorescent protein (GFP), which was used as a model target protein. The protein GFP was fused to the C-terminus of CaM, and a factor Xa cleavage site was introduced between the two proteins. A CaM-GFP fusion protein was expressed in E. coli and purified on a phenothiazine-derivatized silica column. CaM binds to the phenothiazine on the column in a Ca(2+)-dependent fashion and it was, therefore, used as an affinity tail for the purification of GFP. The fusion protein bound to the affinity column was then subjected to a proteolytic digestion with factor Xa. Pure GFP was eluted with a Ca(2+)-containing buffer, while CaM was eluted later with a buffer containing the Ca(2+)-chelating agent EGTA. The purity of the isolated GFP was verified by SDS-PAGE, and the fluorescence properties of the purified GFP were characterized.

Non-NASA Center↗

Cellular localization of the Ca2+ binding TCH3 protein of Arabidopsis

TCH3 is an Arabidopsis touch (TCH) gene isolated as a result of its strong and rapid upregulation in response to mechanical stimuli, such as touch and wind. TCH3 encodes an unusual calcium ion-binding protein that is closely related to calmodulin but has the potential to bind six calcium ions. Here it is shown that TCH3 shows a restricted pattern of accumulation during Arabidopsis vegetative development. These data provide insight into the endogenous signals that may regulate TCH3 expression and the sites of TCH3 action. TCH3 is abundant in the shoot apical meristem, vascular tissue, the root columella and pericycle cells that give rise to lateral roots. In addition, TCH3 accumulation in cells of developing shoots and roots closely correlates with the process of cellular expansion. Following wind stimulation, TCH3 becomes more abundant in specific regions including the branchpoints of leaf primordia and stipules, pith parenchyma, and the vascular tissue. The consequences of TCH3 upregulation by wind are therefore spatially restricted and TCH3 may function at these sites to modify cell or tissue characteristics following mechanical stimulation. Because TCH3 accumulates specifically in cells and tissues that are thought to be under the influence of auxin, auxin levels may regulate TCH3 expression during development. TCH3 is upregulated in response to low levels of exogenous indole-3-acetic acid (IAA), but not by inactive auxin-related compounds. These results suggest that TCH3 protein may play roles in mediating physiological responses to auxin and mechanical environmental stimuli.

NASA Discipline Plant Biology↗

Molecular cloning and characterization of a tumor-associated, growth-related, and time-keeping hydroquinone (NADH) oxidase (tNOX) of the HeLa cell surface

NOX proteins are growth-related cell surface proteins that catalyze both hydroquinone or NADH oxidation and protein disulfide interchange and exhibit prion-like properties. The two enzymatic activities alternate to generate a regular period length of about 24 min. Here we report the expression, cloning, and characterization of a tumor-associated NADH oxidase (tNOX). The cDNA sequence of 1830 bp is located on gene Xq25-26 with an open reading frame encoding 610 amino acids. The activities of the bacterially expressed tNOX oscillate with a period length of 22 min as is characteristic of tNOX activities in situ. The activities are inhibited completely by capsaicin, which represents a defining characteristic of tNOX activity. Functional motifs identified by site-directed mutagenesis within the C-terminal portion of the tNOX protein corresponding to the processed plasma membrane-associated form include quinone (capsaicin), copper and adenine nucleotide binding domains, and two cysteines essential for catalytic activity. Four of the six cysteine to alanine replacements retained enzymatic activity, but the period lengths of the oscillations were increased. A single protein with two alternating enzymatic activities indicative of a time-keeping function is unprecedented in the biochemical literature.

NASA Discipline Cell Biology↗

Differential Processing of Low and High LET Radiation Induced DNA Damage: Investigation of Switch from ATM to ATR Signaling

The members of the phosphatidylinositol kinase-like kinase family of proteins namely ataxia-telangiectasia mutated (ATM) and ATM- and Rad3-related (ATR) are directly responsible for the maintenance of genomic integrity by mounting DDR through signaling and facilitating the recruitment of repair factors at the sites of DNA damage along with coordinating the deployment of cell cycle checkpoints to permit repair by phosphorylating Checkpoint kinase Chk1, Chk2 and p53. High LET radiation from GCR (Galactic Cosmic Rays) consisting mainly of protons and high energy and charged (HZE) particles from SPE (Solar Particle Event) pose a major health risk for astronauts on their space flight missions. The determination of these risks and the design of potential safeguards require sound knowledge of the biological consequences of lesion induction and the capability of the cells to counter them. We here strive to determine the coordination of ATM and ATR kinases at the break sites directly affecting checkpoint signaling and DNA repair and whether differential processing of breaks induced by low and high LET radiation leads to possible augmentation of swap of these damage sensors at the sites of DNA damage. Exposure of cells to IR triggers rapid autophosphorylation of serine-1981 that causes dimer dissociation and initiates monomer formation of ATM. ATM kinase activity depends on the disruption of the dimer, which allows access and phosphorylation of downstream ATM substrates like Chk2. Evidence suggests that ATM is activated by the alterations in higher-order chromatin structure although direct binding of ATM to DSB ends may be a crucial step in its activation. On the other hand, in case of ATR, RPA (replication protein A)-coated ssDNA (single-stranded DNA) generated as a result of stalled DNA replication or during processing of chromosomal lesions is crucial for the localization of ATR to sites of DNA damage in association with ATR-interacting protein (ATRIP). Although the majority of RPA-coated ssDNA is generally present only during DNA replication, ATR activation in G1 and G2-phase might still require formation of RPA-coated ssDNA, probably initiated by the MRN-CtIP complex and then extended by the Exo1- or BLM-dependent mechanisms at the sites of DSBs. Evidence accumulates that activation of ATM and ATR are oppositely regulated by the length of single stranded overhangs generated at the break sites by processes mentioned above and these stretches of single stranded overhangs hold the clue for ATM to ATR switch at broken DNA ends. We irradiated 82-6hTERT human fibroblast cells with low LET gamma-rays and high LET Fe and Si particles. Preliminary results with cells exposed to 1Gy gamma-rays show that the kinetics of pChk2-pT68 foci formation is comparable to that of gamma-H2AX although they appear to recede quicker. The number and intensity of observed foci reaches a maximum at 30 min and 60 min post IR for Chk2-pT68 and gamma-H2AX foci respectively and all Chk2-pT68 foci colocalize with gamma-H2AX foci. The kinetics of Chk1-pS345 and ATRIP are being determined. Results of Chk2-pT68 foci kinetics was also corroborated by western blot experiments, although phosphorylation was detected as early as 10 min and started receding 30 min post IR with 2Gy of gamma-rays. On the other hand, level of ATR-pS428 reached its maximum between 60 and 120 min and was maintained until the last measured time point of 4 hours post IR as determined by western blotting. Experiments performed with high LET Fe and Si particles will be reported.

Saha, Janapriya↗

Ionene modified small polymeric beads

Linear ionene polyquaternary cationic polymeric segments are bonded by means of the Menshutkin reaction (quaternization) to biocompatible, extremely small, porous particles containing halide or tertiary amine sites which are centers for attachment of the segments. The modified beads in the form of emulsions or suspensions offer a large, positively-charged surface area capable of irreversibly binding polyanions such as heparin, DNA, RNA or bile acids to remove them from solution or of reversibly binding monoanions such as penicillin, pesticides, sex attractants and the like for slow release from the suspension.

Rembaum, Alan↗

Effect of salinity on the preservation of biomarkers in hypersaline microbial mat kerogens

Kerogen is the residue that remains after the minerals and organic compounds have been extracted from a rock using acid and organic solvents, respectively [1]. It is not only the largest organic carbon sink in the geological record but the best-preserved fraction due to its recalcitrance and ability to retain syngenetic source information during deposition [2]. Due to the complex, macromolecular, cross-linked polymeric structure of kerogen, it is difficult to analyse intact, often requiring breakdown into smaller subunits through thermal (pyrolysis) or chemical treatment (chemolysis) [3]. Saline and hypersaline environments host diverse microbial communities, which can trap and bind organic matrices such as cellular components (e.g., extracellular polymeric substances, cell wall, peptidoglycan, sheath material, etc.) with sediments, forming layered, laminated structures called microbial mats. Here we present depth profiles of two hypersaline cores through a microbial mat with differing total salinities. We looked at the composition of organic compounds released from kerogen during high temperature, hydrogen-assisted, pyrolysis to understand the effect of salinity on the composition and timing of their incorporation [4]. Field Site: The Exportadora del Sal S. A. (ESSA) saltern of Guerrero Negro (Baja California Sur, Mexico) lies within a typical sabhka environment adjacent to Laguna Ojo de Liebre. Note the spatial extent of the 13 major ponds of the ESSA system. The GN ponds thrive under well-monitored and stable conditions with steady state accretion and degradation of mat layers (around 0.5–1.4 cm per year). Organic mats typically attain up to 10 cm in thickness, depending on the location, which represents about 60 years of growth since these ponds were established. Sedimentation rate is rapid in these non-lithifying and non-mineralising mats, 1–1.5 cm per year. This environment is ideal for investigating in situ lipid preservation and diagenesis without the influence of constantly changing physical parameters. Methods: P4n5 and P5AB cores were collected as ~9 cm and 6.5 cm fresh cores, respectively, in June 2001 and September 2010, respectively from the ESSA saltern. The P4n5 core (9-9.1% salinity) was taken from Pond 4 near 5 and sub-divided into 10 layers. The P5AB core (11.2% salinity) was taken from Pond 5A near 5B and sub-divided into 8 layers. Lyophilised microbial mat powders were solvent-extracted using a modified Bligh-Dyer method and analysed by gas chromatography-mass spectrometry (GC-MS). Additionally, a subset of layers—layer 4 and layer 7—which represented active biomass and sediment-processed, respectively, were subjected to mild-acid methanolysis. The respective residues from conventional solvent extraction and acid methanolysis were loaded with a molybdenum sulfide catalyst and placed into a stainless-steel reactor. Samples were run on the hydropyrolysis set up (heated to 500 °C with 13-15 MPa of H2 pressure), extracted, and analysed by GC-MS. Results: Molecular profiles from the two ponds were significantly different and reflect the contributions of different photosynthetic and respiratory microbial communities to preserved organic matter. For example, the higher salinity P5AB core showed evidence of greater archaeal inputs (similarly observed in lab culture experiments [5]) as well as potential kerogen-bound carotenoids/carotenoid rearranged products or fragments. A hydrophobic emulsion formed during the acid methanolysis processing of layer 4 from both P4n5 and P5AB and layer 7 of P5AB. Hydropyrolysis of this hydrophobic residues released a greater abundance of polycyclic lipids relative to the pre-extracted control. Implications and Future Work: Catalytic hydropyrolysis of kerogen can rapidly generate abundant saturated pyrolysate products from the bound biomarker pool without altering the structures or stereochemistries of the products [6]. Chemical processing of solvent-extractable residues indicated that cellular matrices such as extracellular polymeric substances (EPS) may play a role in the sequestration of polycyclic lipid biomarkers such as steranes and hopanes. Their higher relative abundance compared to the pre-extracted control in the higher salinity layers and cores preliminarily aligns with this hypothesis, although further experiments are required to confirm this. It has been well-documented that EPS plays an important role in mineralisation, specifically carbonate formation, in microbial ecosystems. EPS can enhance calcium carbonate precipitation by providing diffusion-limited sites that create alkalinity gradients in response to microbial processes [7]. It has been experimentally demonstrated that salinity influences the total amount of EPS (both loosely- and tightly-bound) which increases with increasing salinity [8]. The higher salinity at Guerrero Negro switches the microbial population from filamentous Microcoleus to being dominated by Phormidium, Oscillatoria, and unicellular cyanobacteria. We demonstrated that lipid binding into kerogen via strong covalent linkages occurs at the very earliest stages of sedimentary diagenesis. We are currently investigating the influence of salinity on organic preservation through experimental and modelling approaches. Understanding how key biomarkers transform into preserved organic matter (i.e., from precursor biolipids to bound geolipids) in brine ecosystems will aid the search for organic biosignatures on other planetary bodies, especially Icy Moons and modern Mars.

hypersaline↗

B-DNA to Z-DNA structural transitions in the SV40 enhancer: stabilization of Z-DNA in negatively supercoiled DNA minicircles

During replication and transcription, the SV40 control region is subjected to significant levels of DNA unwinding. There are three, alternating purine-pyrimidine tracts within this region that can adopt the Z-DNA conformation in response to negative superhelix density: a single copy of ACACACAT and two copies of ATGCATGC. Since the control region is essential for both efficient transcription and replication, B-DNA to Z-DNA transitions in these vital sequence tracts may have significant biological consequences. We have synthesized DNA minicircles to detect B-DNA to Z-DNA transitions in the SV40 enhancer, and to determine the negative superhelix density required to stabilize the Z-DNA. A variety of DNA sequences, including the entire SV40 enhancer and the two segments of the enhancer with alternating purine-pyrimidine tracts, were incorporated into topologically relaxed minicircles. Negative supercoils were generated, and the resulting topoisomers were resolved by electrophoresis. Using an anti-Z-DNA Fab and an electrophoretic mobility shift assay, Z-DNA was detected in the enhancer-containing minicircles at a superhelix density of -0.05. Fab saturation binding experiments demonstrated that three, independent Z-DNA tracts were stabilized in the supercoiled minicircles. Two other minicircles, each with one of the two alternating purine-pyrimidine tracts, also contained single Z-DNA sites. These results confirm the identities of the Z-DNA-forming sequences within the control region. Moreover, the B-DNA to Z-DNA transitions were detected at superhelix densities observed during normal replication and transcription processes in the SV40 life cycle.

Non-NASA Center↗

Origin of the catalytic activity of bovine seminal ribonuclease against double-stranded RNA

Bovine seminal ribonuclease (RNase) binds, melts, and (in the case of RNA) catalyzes the hydrolysis of double-stranded nucleic acid 30-fold better under physiological conditions than its pancreatic homologue, the well-known RNase A. Reported here are site-directed mutagenesis experiments that identify the sequence determinants of this enhanced catalytic activity. These experiments have been guided in part by experimental reconstructions of ancestral RNases from extinct organisms that were intermediates in the evolution of the RNase superfamily. It is shown that the enhanced interactions between bovine seminal RNase and double-stranded nucleic acid do not arise from the increased number of basic residues carried by the seminal enzyme. Rather, a combination of a dimeric structure and the introduction of two glycine residues at positions 38 and 111 on the periphery of the active site confers the full catalytic activity of bovine seminal RNase against duplex RNA. A structural model is presented to explain these data, the use of evolutionary reconstructions to guide protein engineering experiments is discussed, and a new variant of RNase A, A(Q28L K31C S32C D38G E111G), which contains all of the elements identified in these experiments as being important for duplex activity, is prepared. This is the most powerful catalyst within this subfamily yet observed, some 46-fold more active against duplex RNA than RNase A.

Non-NASA Center↗

Changes in root cap pH are required for the gravity response of the Arabidopsis root

Although the columella cells of the root cap have been identified as the site of gravity perception, the cellular events that mediate gravity signaling remain poorly understood. To determine if cytoplasmic and/or wall pH mediates the initial stages of root gravitropism, we combined a novel cell wall pH sensor (a cellulose binding domain peptide-Oregon green conjugate) and a cytoplasmic pH sensor (plants expressing pH-sensitive green fluorescent protein) to monitor pH dynamics throughout the graviresponding Arabidopsis root. The root cap apoplast acidified from pH 5.5 to 4.5 within 2 min of gravistimulation. Concomitantly, cytoplasmic pH increased in columella cells from 7.2 to 7.6 but was unchanged elsewhere in the root. These changes in cap pH preceded detectable tropic growth or growth-related pH changes in the elongation zone cell wall by 10 min. Altering the gravity-related columella cytoplasmic pH shift with caged protons delayed the gravitropic response. Together, these results suggest that alterations in root cap pH likely are involved in the initial events that mediate root gravity perception or signal transduction.

NASA Discipline Plant Biology↗

Crystal structures of pinoresinol-lariciresinol and phenylcoumaran benzylic ether reductases and their relationship to isoflavone reductases

Despite the importance of plant lignans and isoflavonoids in human health protection (e.g. for both treatment and prevention of onset of various cancers) as well as in plant biology (e.g. in defense functions and in heartwood development), systematic studies on the enzymes involved in their biosynthesis have only recently begun. In this investigation, three NADPH-dependent aromatic alcohol reductases were comprehensively studied, namely pinoresinol-lariciresinol reductase (PLR), phenylcoumaran benzylic ether reductase (PCBER), and isoflavone reductase (IFR), which are involved in central steps to the various important bioactive lignans and isoflavonoids. Of particular interest was in determining how differing regio- and enantiospecificities are achieved with the different enzymes, despite each apparently going through similar enone intermediates. Initially, the three-dimensional x-ray crystal structures of both PLR_Tp1 and PCBER_Pt1 were solved and refined to 2.5 and 2.2 A resolutions, respectively. Not only do they share high gene sequence similarity, but their structures are similar, having a continuous alpha/beta NADPH-binding domain and a smaller substrate-binding domain. IFR (whose crystal structure is not yet obtained) was also compared (modeled) with PLR and PCBER and was deduced to have the same overall basic structure. The basis for the distinct enantio-specific and regio-specific reactions of PCBER, PLR, and IFR, as well as the reaction mechanism and participating residues involved (as identified by site-directed mutagenesis), are discussed.

NASA Discipline Plant Biology↗

The Cytoskeleton and Force Response Mechanisms

The long term aim of this project was to define the mechanisms by which cells sense and respond to the physical forces experienced at 1g and missing in microgravity. Identification and characterization of the elements of the cells force response mechanism could provide pathways and molecules to serve as targets for pharmacological intervention to mitigate the pathologic effects of microgravity. Mechanical forces experienced by the organism can be transmitted to cells through molecules that allow cells to bind to the extracellular matrix and through other types of molecules which bind cells to each other. These molecules are coupled in large complexes of proteins to structural elements such as the actin cytoskeleton that give the cell the ability to sense, resist and respond to force. Application of small forces to tissue culture cells causes local elevation of intracellular calcium through stretch activated ion channels, increased tyrosine phosphorylation and a restructuring of the actin cytoskeleton. Using collagen coated iron oxide beads and strong magnets, we can apply different levels of force to cells in culture. We have found that force application causes the cells to polymerize actin at the site of mechanical deformation and unexpectedly, to depolymerize actin across the rest of the cell. Observations of GFP- actin expressing cells demonstrate that actin accumulates at the site of deformation within the first five minutes of force application and is maintained for many tens of minutes after force is removed. Consistent with the reinforcement of the cytoskeletal structures underlying the integrin-bead interaction, force also alters the motion of bound magnetic beads. This effect is seen following the removal of the magnetic field, and is only partially ablated by actin disruption with cytochalsin B. While actin is polymerizing locally at the site of force application, force also stimulates a global reduction in actin filament content within the cells. We have examined the roles of several actin filament disassembly factors in the global reduction of cellular actin filaments. The calcium regulated actin filament severing protein gelsolin is not necessary for the increased actin turnover, as cells derived from gelsolin null and wildtype mice still show a reduction in total actin filament content. Instead, our work suggests that the actin binding protein cofilin may be important for these changes in actin dynamics. Cofilin binds to and enhances the disassembly of actin filaments. Using immunological methods, we observe transient changes in the phosphorylation state of cofilin upon force application that suggests that cofilin may mediate actin filament turnover. Early after force application, cofilin is transiently dephosphorylated, activating its actin disassembly activity. Subsequently, we find a hyper-phosphorylation of cofilin, rendering it inactive. This reduction in cofilin activity may explain the stability of the force induced actin structuttes. In testing this hypothesis, we aimed to generate cells that express the constituitively active kinase (LIM-kinase) that phosphorylates cofilin. lnidial attempts in the cell lines used for the our previous studies proved unsuccessful. While we prepare this work for pubication, we are continuing to study other cell lines and tissue sources to determine whether they show a reduction in F-actin content after force application.

Allen, Philip Goodwin↗

Characterization of the growth and auxin physiology of roots of the tomato mutant, diageotropica

Roots of the tomato (Lycopersicon esculentum, Mill.) mutant (diageotropica (dgt) exhibit an altered phenotype. These roots are agravitropic and lack lateral roots. Relative to wild-type (VFN8) roots, dgt roots are less sensitive to growth inhibition by exogenously applied IAA and auxin transport inhibitors (phytotropins), and the roots exhibit a reduction in maximal growth inhibition in response to ethylene. However, IAA transport through roots, binding of the phytotropin, tritiated naphthylphthalamic acid ([3H]NPA), to root microsomal membranes, NPA-sensitive IAA uptake by root segments, and uptake of [3H]NPA into root segments are all similar in mutant and wild-type roots. We speculate that the reduced sensitivity of dgt root growth to auxin-transport inhibitors and ethylene is an indirect result of the reduction in sensitivity to auxin in this single gene, recessive mutant. We conclude that dgt roots, like dgt shoots, exhibit abnormalities indicating they have a defect associated with or affecting a primary site of auxin perception or action.

NASA Discipline Plant Biology↗

REACTER 2.0: Quantum-Informed Reaction Constraints and Automated Interaction Typing

REACTER is a heuristic method for modeling chemical reactions in classical molecular dynamics simulations, implemented in LAMMPS as fix bond/react. The authors recently extended LAMMPS to support alphanumeric labels for atom types, bond types etc., which enables the pre- and post-reaction templates required by the REACTER protocol to be portable between different simulations and greatly simplifies the task of creating simulation-ready reaction templates. To further increase the generality of reaction templates, support for wildcard characters within atom types has been added, along with the automatic assignment of interaction types for new bonds, angles, etc. based on the involved atom types. In some cases, this feature can express a class of reactions with one pair of reaction templates, where previously dozens may have been required. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, as well as a custom constraint that leverages LAMMPS’ powerful built-in variable framework. Other new features include variable support for various inputs (e.g., to allow reaction rates or cutoffs to be dependent on overall conversion), on-the-fly update of molecule IDs, and the ability to create new atoms positioned with respect to the reaction site. The new features are applied to modeling polymeric, thermosetting and composite materials, and advanced applications of the new reaction constraints are demonstrated. For example, REACTER is shown to accurately reproduce mechanically-induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a constraint on the total potential energy of the involved atoms.

polymer simulations↗

Studies of H2O on beta-AgI surfaces - An effective pair potential model

The adsorption of a water molecule on surfaces of beta-AgI, the hexagonal crystal believed to be primarily responsible for the ice-nucleating properties of AgI, is studied on the basis of an effective pair potential model. The water molecule is represented by a rigid point charge ST-2 model and the AgI substrate by an array of point atoms, and maximal binding energy surfaces and optimal H2O configurations are generated for the water molecule adsorbed on the rigid and unrelaxed basal and prism faces. Modeling of H2O adsorption above a two-layer ledge, an iodine vacancy and an H2O molecule trapped in the vacancy indicates that H2O adsorption is favored at interstitial sites where no substrate atoms lie directly below. The prism face is found to attract the water molecule more strongly and provide larger energy barriers to surface diffusion than basal face sites, with the ideal basal faces providing hexagonal patterns of adsorption sites for the H2O with preferred dipole moments aligned.

Hale, B. N.↗

Shielding, the bulk chemical potential, and cohesion in alloys

It is shown that the bulk chemical potential in alloys is intimately related to the spatial dependence of the shielding cloud that results when the electronic charge rearranges itself as one atom type is replaced by another at a given site. Such a relationship fixes the relative energy scale between the alloy and its pure constituents, important in determining the stability of alloys. A correct treatment of shielding is thus essential to quantitative calculations of alloy stability. A model calculation of the bulk chemical potential and cohesion of alloys in the tight-binding approximation is presented as a numerical example. In the course of this investigation a general invariant of an integral over the shielding cloud is derived.

Stern, E. A.↗