A parallel multistate framework for atomistic non-equilibrium reaction dynamics of solutes in strongly interacting organic solvents

We describe a parallel linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM. Forces are obtained using the Hellman-Feynmann relationship, giving continuous gradients, and excellent energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to CCSD(T)-F12 electronic structure theory, we built a 64-state MS-EVB model designed to study the F + CD3CN -> DF + CD2CN reaction in CD3CN solvent. This approach allows us to build a reactive potential energy surface (PES) whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We use our PES to run MD simulations, and examine a range of transient observables which follow in the wake of reaction, including transient spectra of the DF vibrational band, time dependent profiles of vibrationally excited DF in CD3CN solvent, and relaxation rates for energy flow from DF into the solvent, all of which agree well with experimental observations. Immediately following deuterium abstraction, the nascent DF is in a non-equilibrium regime in two different respects: (1) it is highly excited, with ~23 kcal mol-1 localized in the stretch; and (2) not yet Hydrogen bonded to the CD3CN solvent, its microsolvation environment is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of its microsolvation environment results in a red shift. These two competing effects result in a post-reaction relaxation profile distinct from that observed when DF vibration excitation occurs within an equilibrium microsolvation environment. The parallel software framework presented in this paper should be more broadly applicable to a range of complex reactive systems.Comment: 58 pages and 29 Figure

Mechanics of lipid bilayer junctions affecting the size of a connecting lipid nanotube

In this study we report a physical analysis of the membrane mechanics affecting the size of the highly curved region of a lipid nanotube (LNT) that is either connected between a lipid bilayer vesicle and the tip of a glass microinjection pipette (tube-only) or between a lipid bilayer vesicle and a vesicle that is attached to the tip of a glass microinjection pipette (two-vesicle). For the tube-only configuration (TOC), a micropipette is used to pull a LNT into the interior of a surface-immobilized vesicle, where the length of the tube L is determined by the distance of the micropipette to the vesicle wall. For the two-vesicle configuration (TVC), a small vesicle is inflated at the tip of the micropipette tip and the length of the tube L is in this case determined by the distance between the two interconnected vesicles. An electrochemical method monitoring diffusion of electroactive molecules through the nanotube has been used to determine the radius of the nanotube R as a function of nanotube length L for the two configurations. The data show that the LNT connected in the TVC constricts to a smaller radius in comparison to the tube-only mode and that tube radius shrinks at shorter tube lengths. To explain these electrochemical data, we developed a theoretical model taking into account the free energy of the membrane regions of the vesicles, the LNT and the high curvature junctions. In particular, this model allows us to estimate the surface tension coefficients from R(L) measurements

Recombination, Solvation and Reaction of CN Radicals Following Ultraviolet Photolysis of ICN in Organic Solvents

The fates of CN radicals produced by ultraviolet (UV) photolysis of ICN in various organic solvents have been examined by transient electronic and vibrational absorption spectroscopy (TEAS and TVAS). Near-UV and visible bands in the TEAS measurement enable direct observation of the CN radicals and their complexes with the solvent molecules. Complementary TVAS measurements probe the products of CN–radical reactions. Geminate recombination to form ICN and INC is a minor pathway on the 150 fs −1300 ps time scales of our experiments in the chosen organic solvents; nonetheless, large infrared transition dipole moments permit direct observation of INC that is vibrationally excited in the CN stretching mode. The time constants for INC vibrational cooling range from 30 ps in tetrahydrofuran (THF) to 1400 ps in more weakly interacting solvents such as chloroform. The major channel for CN removal in the organic solvents is reaction with solvent molecules, as revealed by depletion of solvent absorption bands and growth of product bands in the TVA spectra. HCN is a reaction product of hydrogen atom abstraction in most of the photoexcited solutions, and forms with vibrational excitation in both the C–H and CN stretching modes. The vibrational cooling rate of the CN stretch in HCN depends on the solvent, and follows the same trend as the cooling rate of the CN stretch in INC. However, in acetonitrile solution an additional reaction pathway produces C₃H₃N₂^• radicals, which release HCN on a much longer time scale

Contrasting ring-opening propensities in UV-excited α-pyrone and coumarin

Ring-opening quantum yields following UV-photoexcitation of coumarin and α-pyrone are influenced by the dynamics through, rather than just the geometries of, conical intersections.</p

Probing the ultrafast energy dissipation mechanism of the sunscreen oxybenzone after UVA irradiation

Oxybenzone is a common constituent of many commercially available sunscreens providing photoprotection from ultraviolet light incident on the skin. Femtosecond transient electronic and vibrational absorption spectroscopies have been used to investigate the non-radiative relaxation pathways of oxybenzone in cyclohexane and methanol after excitation in the UVA region. The present data suggest that the photoprotective properties of oxybenzone can be understood in terms of an initial ultrafast excited state enol -> keto tautomerization, followed by efficient internal conversion and subsequent vibrational relaxation to the ground state (enol) tautomer

Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury

<p>Abstract</p> <p>Background</p> <p>The sequencing of the <it>D.melanogaster </it>genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of <it>Drosophila </it>protein-coding genes contain one or more alternative exons. A recent transcription map of the <it>Drosophila </it>embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of <it>Drosophila </it>transcriptome.</p> <p>Results</p> <p>Bioinformatic analysis of 1,303 <it>Drosophila </it>ORESTES clusters identified 68 sequences derived from unannotated regions in the current <it>Drosophila </it>genome version (4.3). Of these, a set of 38 was analysed by polyA⁺northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The <it>SP212 </it>gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this <it>locus </it>is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury.</p> <p>Conclusion</p> <p>Using the ORESTES methodology we identified 17 novel exons from low abundance <it>Drosophila </it>transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.</p

Predisposition to Cancer Caused by Genetic and Functional Defects of Mammalian Atad5

ATAD5, the human ortholog of yeast Elg1, plays a role in PCNA deubiquitination. Since PCNA modification is important to regulate DNA damage bypass, ATAD5 may be important for suppression of genomic instability in mammals in vivo. To test this hypothesis, we generated heterozygous (Atad5+/m) mice that were haploinsuffficient for Atad5. Atad5+/m mice displayed high levels of genomic instability in vivo, and Atad5+/m mouse embryonic fibroblasts (MEFs) exhibited molecular defects in PCNA deubiquitination in response to DNA damage, as well as DNA damage hypersensitivity and high levels of genomic instability, apoptosis, and aneuploidy. Importantly, 90% of haploinsufficient Atad5+/m mice developed tumors, including sarcomas, carcinomas, and adenocarcinomas, between 11 and 20 months of age. High levels of genomic alterations were evident in tumors that arose in the Atad5+/m mice. Consistent with a role for Atad5 in suppressing tumorigenesis, we also identified somatic mutations of ATAD5 in 4.6% of sporadic human endometrial tumors, including two nonsense mutations that resulted in loss of proper ATAD5 function. Taken together, our findings indicate that loss-of-function mutations in mammalian Atad5 are sufficient to cause genomic instability and tumorigenesis