Application of isotopic labeling,and gas chromatography mass spectrometry,to understanding degradation products and pathways in the thermal-oxidative aging of Nylon 6.6

Nylon 6.6 containing ¹³C isotopic labels at specific positions along the macromolecular backbone has been subjected to extensive thermal-oxidative aging at 138 °C for time periods up to 243 days. In complementary experiments, unlabeled Nylon 6.6 was subjected to the same aging conditions under an atmosphere of ¹⁸O₂. Volatile organic degradation products were analyzed by cryofocusing gas chromatography mass spectrometry (cryo-GC/MS) to identify the isotopic labeling. The labeling results, combined with basic considerations of free radical reaction chemistry, provided insights to the origin of degradation species, with respect to the macromolecular structure. A number of inferences on chemical mechanisms were drawn, based on 1) the presence (or absence) of the isotopic labels in the various products, 2) the location of the isotope within the product molecule, and 3) the relative abundance of products as indicated by large differences in peak intensities in the gas chromatogram. The overall degradation results can be understood in terms of free radical pathways originating from initial attacks on three different positions along the nylon chain which include hydrogen abstraction from: the (CH₂) group adjacent to the nitrogen atom, at the (CH₂) adjacent the carbonyl group, and direct radical attack on the carbonyl. Understanding the pathways which lead to Nylon 6.6 degradation ultimately provides new insight into changes that can be leveraged to detect and reduce early aging and minimize problems associated with material degradation.     

Accurate double many‐body expansion potential energy surface by extrapolation to the complete basis set limit and dynamics calculations for ground state of NH2

An accurate single‐sheeted double many‐body expansion potential energy surface is reported for the title system. A switching function formalism has been used to warrant the correct behavior at the and dissociation channels involving nitrogen in the ground and first excited states. The topographical features of the novel global potential energy surface are examined in detail, and found to be in good agreement with those calculated directly from the raw ab initio energies, as well as previous calculations available in the literature. The novel surface can be using to treat well the Renner–Teller degeneracy of the and states of . Such a work can both be recommended for dynamics studies of the reaction and as building blocks for constructing the double many‐body expansion potential energy surface of larger nitrogen/hydrogen‐containing systems. In turn, a test theoretical study of the reaction has been carried out with the method of quantum wave packet on the new potential energy surface. Reaction probabilities, integral cross sections, and differential cross sections have been calculated. Threshold exists because of the energy barrier (68.5 meV) along the minimum energy path. On the curve of reaction probability for total angular momentum J = 0, there are two sharp peaks just above threshold. The value of integral cross section increases quickly from zero to maximum with the increase of collision energy, and then stays stable with small oscillations. The differential cross section result shows that the reaction is a typical forward and backward scatter in agreement with experimental measurement result. © 2013 Wiley Periodicals, Inc.     

Local Hartree–Fock orbitals using a three‐level optimization strategy for the energy

Using the three‐level energy optimization procedure combined with a refined version of the least‐change strategy for the orbitals—where an explicit localization is performed at the valence basis level—it is shown how to more efficiently determine a set of local Hartree–Fock orbitals. Further, a core–valence separation of the least‐change occupied orbital space is introduced. Numerical results comparing valence basis localized orbitals and canonical molecular orbitals as starting guesses for the full basis localization are presented. The results show that the localization of the occupied orbitals may be performed at a small computational cost if valence basis localized orbitals are used as a starting guess. For the unoccupied space, about half the number of iterations are required if valence localized orbitals are used as a starting guess compared to a canonical set of unoccupied Hartree–Fock orbitals. Different local minima may be obtained when different starting guesses are used. However, the different minima all correspond to orbitals with approximately the same locality. © 2013 Wiley Periodicals, Inc.     

Computational gibberellin‐binding channel discovery unraveling the unexpected perception mechanism of hormone signal by gibberellin receptor

Gibberellins (GAs) are phytohormones essential for many developmental processes in plants. In this work, fundamental mechanism of hormone perception by receptor GID1 has been studied by performing computational simulations, revealing a new GA‐binding channel of GID1 and a novel hormone perception mechanism involving only one conformational state of GID1. The novel hormone perception mechanism demonstrated here is remarkably different from the previously proposed/speculated mechanism [Murase et al., Nature 2008 , 456, 459] involving two conformational states (“OPEN” and “CLOSED”) of GID1. According to the new perception mechanism, GA acts as a “conformational stabilizer,” rather than the previously speculated “allosteric inducer,” to induce the recognition of protein DELLA by GID1. The novel mechanistic insights obtained in this study provide a new starting point for further studies on the detailed molecular mechanisms of GID1 interacting with DELLA and various hormones and for mechanism‐based rational design of novel, potent growth regulators that target crops and ornamental plants. © 2013 Wiley Periodicals, Inc.     

Multiscale geometric modeling of macromolecules II: Lagrangian representation

Geometric modeling of biomolecules plays an essential role in the conceptualization of biolmolecular structure, function, dynamics, and transport. Qualitatively, geometric modeling offers a basis for molecular visualization, which is crucial for the understanding of molecular structure and interactions. Quantitatively, geometric modeling bridges the gap between molecular information, such as that from X‐ray, NMR, and cryo‐electron microscopy, and theoretical/mathematical models, such as molecular dynamics, the Poisson–Boltzmann equation, and the Nernst–Planck equation. In this work, we present a family of variational multiscale geometric models for macromolecular systems. Our models are able to combine multiresolution geometric modeling with multiscale electrostatic modeling in a unified variational framework. We discuss a suite of techniques for molecular surface generation, molecular surface meshing, molecular volumetric meshing, and the estimation of Hadwiger's functionals. Emphasis is given to the multiresolution representations of biomolecules and the associated multiscale electrostatic analyses as well as multiresolution curvature characterizations. The resulting fine resolution representations of a biomolecular system enable the detailed analysis of solvent–solute interaction, and ion channel dynamics, whereas our coarse resolution representations highlight the compatibility of protein‐ligand bindings and possibility of protein–protein interactions. © 2013 Wiley Periodicals, Inc.     

Transiting the molecular potential energy surface along low energy pathways: The TRREAT algorithm

The Transition Rapidly exploring Random Eigenvector Assisted Tree (TRREAT) algorithm is introduced to perform searches along low curvature pathways on a potential energy surface (PES). The method combines local curvature information about the PES with an iterative Rapidly exploring Random Tree algorithm (LaValle, Computer Science Department, Iowa State University, 1998, TR98–11) that quickly searches high‐dimensional spaces for feasible pathways between local minima. Herein, the method is applied to identifying conformational changes of molecular systems using Cartesian coordinates while avoiding a priori definition of collective variables. We analyze the pathway identification problem for alanine dipeptide, cyclohexane and glycine using nonreactive and reactive forcefields. We show how TRREAT‐identified pathways can be used as valuable input guesses for double‐ended methods such as the Nudged Elastic Band when ascertaining transition state energies. This method can be utilized to improve/extend the reaction databases that lie at the core of automatic chemical reaction mechanism generator software currently developed to build kinetic models of chemical reactions. © 2013 Wiley Periodicals, Inc.     

FEW: A workflow tool for free energy calculations of ligand binding

In the later stages of drug design projects, accurately predicting relative binding affinities of chemically similar compounds to a biomolecular target is of utmost importance for making decisions based on the ranking of such compounds. So far, the extensive application of binding free energy approaches has been hampered by the complex and time‐consuming setup of such calculations. We introduce the free energy workflow (FEW) tool that facilitates setup and execution of binding free energy calculations with the AMBER suite for multiple ligands. FEW allows performing free energy calculations according to the implicit solvent molecular mechanics (MM‐PB(GB)SA), the linear interaction energy, and the thermodynamic integration approaches. We describe the tool's architecture and functionality and demonstrate in a show case study on Factor Xa inhibitors that the time needed for the preparation and analysis of free energy calculations is considerably reduced with FEW compared to a fully manual procedure. © 2013 Wiley Periodicals, Inc.     

Understanding the electronic energy transfer pathways in the trimeric and hexameric aggregation state of cyanobacteria phycocyanin within the framework of Förster theory

In the present study, the electronic energy transfer pathways in trimeric and hexameric aggregation state of cyanobacteria C‐phycocyanin (C‐PC) were investigated in term of the Förster theory. The corresponding excited states and transition dipole moments of phycocyanobilins (PCBs) located into C‐PC were examined by model chemistry in gas phase at time‐dependent density functional theory (TDDFT), configuration interaction‐singles (CIS), and Zerner's intermediate neglect of differential overlap (ZINDO) levels, respectively. Then, the long‐range pigment‐protein interactions were approximately taken into account by using polarizable continuum model (PCM) at TDDFT level to estimate the influence of protein environment on the preceding calculated physical quantities. The influence of the short‐range interaction caused by aspartate residue nearby PCBs was examined as well. Only when the protonation of PCBs and its long‐ and short‐range interactions were properly taken into account, the calculated energy transfer rates (1/K) in the framework of Förster model at TDDFT/B3LYP/6‐31+G* level were in good agreement with the experimental results of C‐PC monomer and trimer. Furthermore, the present calculated results suggested that the energy transfer pathway in C‐PC monomer is predominant from β‐155 to β‐84 (1/K = 13.4 ps), however, from α‐84 of one monomer to β‐84 (1/K = 0.3–0.4 ps) in a neighbor monomer in C‐PC trimer. In C‐PC hexamer, an additional energy flow was predicted to be from β‐155 (or α‐84) in top trimer to adjacent β‐155 (or α‐84) (1/K = 0.5–2.7 ps) in bottom trimer. © 2013 Wiley Periodicals, Inc.