Photochemical Formation of Diazenecarbaldehyde (HNNCHO) and Diazenecarbothialdehyde (HNNCHS) in Low-Temperature Matrices.

The photochemical decomposition of 1,2,4-oxadiazole-3,5-diamine and 1,2,4-thiadiazole-3,5-diamine was investigated in low-temperature Ar and Kr matrixes at different wavelengths. The analysis of matrix-isolation infrared (MI-IR) spectra aided by high-level quantum chemical computations showed not only that these photochemical reactions yield [NH2, C, N, X] (X = O, S) isomers but also that the bands of a novel, formerly unobserved species were observed. The comparison of computed IR spectra of potential products with the observed spectra suggests that these species are the diazenecarbaldehyde (HNNCHO) and diazenecarbothialdehyde (HNNCHS). Neither of the reactive HNNCHO and HNNCHS molecules was observed experimentally before. Both molecules are identified in the matrix as a complex with the other photoproduct, NH2CN. Comparison of the present experiments with former photochemical experiments on 1,2,5-oxadiazole-3,4-diamine and 1,2,5-thiadiazole-3,4-diamine and the analysis of the rate of formation of the different photoproducts indicate that HNNCHO and HNNCHS are formed in a different reaction path than H2NNCX and H2NC(NX) (X = O, S), and not by photoisomerization from these latter products.

Photochemical generation of H2NCNX, H2NNCX, H2NC(NX) (X = O, S) in low-temperature matrices.

The [NH2, C, N, O] and the [NH2, C, N, S] molecular systems were investigated by computational and matrix-isolation spectroscopic methods. The determination of the equilibrium structures and relative energies by CCSD(T) method was followed by the computation of the harmonic and anharmonic vibrational wavenumbers, infrared intensities, relative Raman activities, and UV excitation energies. These computed data were used to assist the identification of products obtained by UV laser photolysis of 3,4-diaminofurazan and 3,4-diaminothiadiazole in low-temperature Ar and Kr matrices. It is shown that two open-chain H2NNCX and H2NCNX and one cyclic H2NC(NX) (X = O, S) isomers are generated in the case of both systems. Except for H2NNCO and H2NCNS, the present study reports the first generation and spectroscopic identification of these compounds.

Huntingtin exon 1 fibrils feature an interdigitated ß-hairpin-based polyglutamine core.

Polyglutamine expansion within the exon1 of huntingtin leads to protein misfolding, aggregation, and cytotoxicity in Huntington's disease. This incurable neurodegenerative disease is the most prevalent member of a family of CAG repeat expansion disorders. Although mature exon1 fibrils are viable candidates for the toxic species, their molecular structure and how they form have remained poorly understood. Using advanced magic angle spinning solid-state NMR, we directly probe the structure of the rigid core that is at the heart of huntingtin exon1 fibrils and other polyglutamine aggregates, via measurements of long-range intramolecular and intermolecular contacts, backbone and side-chain torsion angles, relaxation measurements, and calculations of chemical shifts. These experiments reveal the presence of ß-hairpin-containing ß-sheets that are connected through interdigitating extended side chains. Despite dramatic differences in aggregation behavior, huntingtin exon1 fibrils and other polyglutamine-based aggregates contain identical ß-strand-based cores. Prior structural models, derived from X-ray fiber diffraction and computational analyses, are shown to be inconsistent with the solid-state NMR results. Internally, the polyglutamine amyloid fibrils are coassembled from differently structured monomers, which we describe as a type of "intrinsic" polymorphism. A stochastic polyglutamine-specific aggregation mechanism is introduced to explain this phenomenon. We show that the aggregation of mutant huntingtin exon1 proceeds via an intramolecular collapse of the expanded polyglutamine domain and discuss the implications of this observation for our understanding of its misfolding and aggregation mechanisms.

VCD Robustness of the Amide-I and Amide-II Vibrational Modes of Small Peptide Models.

The rotational strengths and the robustness values of amide-I and amide-II vibrational modes of For(AA)n NHMe (where AA is Val, Asn, Asp, or Cys, n = 1-5 for Val and Asn; n = 1 for Asp and Cys) model peptides with α-helix and ß-sheet backbone conformations were computed by density functional methods. The robustness results verify empirical rules drawn from experiments and from computed rotational strengths linking amide-I and amide-II patterns in the vibrational circular dichroism (VCD) spectra of peptides with their backbone structures. For peptides with at least three residues (n ≥ 3) these characteristic patterns from coupled amide vibrational modes have robust signatures. For shorter peptide models many vibrational modes are nonrobust, and the robust modes can be dependent on the residues or on their side chain conformations in addition to backbone conformations. These robust VCD bands, however, provide information for the detailed structural analysis of these smaller systems.

Near-infrared laser-induced structural changes of glycine·water complexes in an Ar matrix.

The structures of glycine·H2O complexes have been reinvestigated in low-temperature inert matrices. To go beyond the former matrix-isolation IR studies, NIR laser irradiation was used to change the relative abundances of the different complexes in the matrix. It is shown that the irradiation of the first overtone of the OH stretching mode of glycine as well as of the first overtone of the OH stretching mode of the water molecule in the complex can induce structural changes. Comparison of the experimental IR spectra with the IR spectra computed for different structures resulted in more reliable assignments of spectral patterns and identification of more structures than in former studies.

Exploring the conformational space of cysteine by matrix isolation spectroscopy combined with near-infrared laser induced conformational change.

Six conformers of α-cysteine were identified by matrix isolation IR spectroscopy combined with NIR laser irradiation. Five of these conformers are identical with the five out of six conformers that have recently been identified by microwave spectroscopy. The sixth conformer observed in the present study is a short-lived conformer, which decays by H-atom tunneling; its half-life in a 12 K N2 matrix is (1.1 ± 0.5) × 10(3) s. This study proves that matrix isolation IR spectroscopy combined with NIR laser irradiation is a suitable method to identify conformers of a complex system for which computations predict several dozens of conformers, and that the reliability of this method for conformational assignment is comparable to that of microwave spectroscopy.

Near-infrared laser induced conformational change of alanine in low-temperature matrixes and the tunneling lifetime of its conformer VI.

The near- and mid-IR spectra of α-alanine isolated in low-temperature Ar, Kr, and N2 matrixes were measured. Production of the short-lived conformer VI at the expense of the predominant conformer I was observed upon short irradiation with NIR laser light at the first O-H stretching overtone band of conformer I. Conformer VI decays by H-atom tunneling at 12 K with half-lives of 5.7 ± 1 s, 2.8 ± 1 s in Ar (two different sites), 7.0 ± 1 s in Kr, and 2.8 × 10(3) ± 1.2 × 10(3) s in N2. Upon prolonged irradiation, conformer I slowly transformed into conformer IIa. On the basis of these irradiation experiments, the unambiguous vibrational assignments of conformers I, IIa, and VI are given. In contrast to similar experiments for glycine, the irradiation experiments did not lead to the formation of conformer IIIb. This is explained by a very low IIIb → I barrier height computed for alanine, which results in a very fast depletion of conformer IIIb even in low-temperature matrixes.

Tunneling lifetime of the ttc/VIp conformer of glycine in low-temperature matrices.

Conformer ttc/VIp of glycine and glycine-N,N,O-d(3) has been prepared in low-temperature Ar, Kr, Xe, and N(2) matrices by near-infrared (NIR) laser irradiation of the first OH stretching overtone of conformer ttt/Ip. Glycine (and glycine-N,N,O-d(3)) ttc/VIp was found to convert back to ttt/Ip in the dark by hydrogen-atom tunneling. The observed half-lives of ttc/VIp in Ar, Kr, and Xe matrices at 12 K were 4.4 ± 1 s (50.0 ± 1 h), 4.0 ± 1 s (48.0 ± 1 h), and 2.8 ± 1 s (99.3 ± 2 h), respectively. In correspondence with the observation for the cis-to-trans conversion of formic and acetic acid, the tunneling half-life of glycine ttc/VIp in a N(2) matrix is more than 3 orders of magnitude longer (6.69 × 10(3) and 1.38 × 10(4) s for two different sites) than in noble gas matrices due to complex formation with the host molecules. The present results are important to understand the lack of experimental observation of some computationally predicted conformers of glycine and other amino acids.

Reliability of computed signs and intensities for vibrational circular dichroism spectra.

We present in detail a novel measure that improves the reliability of the assignment procedure for vibrational circular dichroism (VCD) spectra extending the useful robustness concept introduced by Nicu and Baerends. This measure enables spectroscopists to single out bands with unreliable VCD intensities that can be disregarded during analysis and determination of absolute configuration. The previously proposed robustness criterion is shown to be gauge dependent and less reliable than the one proposed here.

Effects of strong and weak hydrogen bond formation on VCD spectra: a case study of 2-chloropropionic acid.

The vibrational circular dichroism (VCD) spectrum of S-(-) and R-(+)-2-chloropropionic acid is thoroughly analyzed. Besides the VCD spectrum of the monomer, the dimers (stabilized by strong hydrogen bonds) and the 2-chloropropionic acid-CHCl(3) complexes (stabilized by a weak hydrogen bond) are studied both experimentally (in solution and in low-temperature Ar matrix) and by quantum chemical computations. It is shown that dimer formation drastically changes, and even weak complex formation can also substantially affect the overall shape of the VCD spectrum. The present and previous results can be generalized for the practice of absolute configuration determination of carboxylic acids by VCD spectroscopy. For these measurements, if bulky groups do not block dimer formation, comparison of the computed spectra of the dimers with the experimental spectra recorded in relatively concentrated (∼0.1 mol dm(-3)) solutions is suggested. Our study also shows that due to the stabilization of monomers and/or the formation of weak complexes, the VCD spectrum recorded in CHCl(3) is more complex and, like in the present case, can have a lower intensity than that of the spectrum recorded in CCl(4). Therefore, if solubility allows, CCl(4) is a much preferred solvent over CHCl(3).

Gas-phase and Ar-matrix SQM scaling factors for various DFT functionals with basis sets including polarization and diffuse functions.

Scaling factors for Pulay's scaled quantum mechanical (SQM) scheme have been determined for four different widely used DFT functionals (PBE, B3LYP, B3PW91, and M06-2X) and for two basis sets (6-31++G** and aug-cc-pVTZ) by fitting computed results to 347 fundamental experimental vibrational frequencies of 33 molecules. Measurements in the gas phase and in solid argon matrices were used independently in the fitting procedure in order to provide a simple method of estimating matrix shifts. The accuracy of the new scaling factors is demonstrated on test molecules including hydrogen-bonded systems and molecules containing chlorine and sulfur atoms.

Is ß-homo-proline a pseudo-γ-turn forming element of ß-peptides? An IR and VCD spectroscopic study on Ac-ß-HPro-NHMe in cryogenic matrices and solutions.

In order to test the pseudo-γ-turn forming capability of ß-homo-proline (ß(3)-HPro) 2-[(2S)-1-acetylpyrrolidin-2-yl]-N-methylacetamide (Ac-ß(3)-HPro-NHMe) was synthesized and its potential energy landscape was investigated by infrared (IR) and vibrational circular dichroism (VCD) spectroscopy combined with density functional calculations. Based upon a comparison between experimental and computed spectra three different pseudo-γ-turn-like trans conformers and a cis conformer were identified in low-temperature Ar and Kr matrices. The computations in agreement with the observations reveal that, in contrast to its α-Pro analogue, the room-temperature abundance of the cis conformer is significant, falling above 10% in the isolated phase. Furthermore, solution-phase vibrational spectra and computations show that the cis conformer is predominant in polar solvents. This result indicates that ß(3)-HPro is significantly less apt to form pseudo-γ-turns when compared to the γ-turn forming tendency of α-proline. The present study also shows that the interpretation of solution-phase VCD spectra of flexible molecules should be done with extra caution.

Chiral and achiral fundamental conformational building units of beta-peptides: a matrix isolation conformational study on Ac-beta-HGly-NHMe and Ac-beta-HAla-NHMe.

The infrared spectra of two model beta-peptides, N-acetyl-3-aminopropionic acid-N'-methylamide (Ac-beta-HGly-NHMe) and N-acetyl-3-aminobutanoic acid-N'-methylamide (Ac-beta-HAla-NHMe), have been recorded in low-temperature Ar and Kr matrixes. The spectra were assigned by the help of electronic structure calculations. The analysis of spectra, in line with the theoretical predictions, revealed that both biocompatible peptide building blocks have a single dominant backbone conformer. Besides this prevalent conformer, which has a six-membered H-bonded pseudoring, conformers with eight-membered H-bonded rings are also observed but in a significantly smaller amount. The calculated conformer distribution is consistent with the experimental findings. The present work along with other recent results supports the concept that the backbone conformation of longer biopolymers, such as alpha- and beta-peptides, can be deciphered using the conformers of their structural building blocks. In this respect, our conformational study on the simplest models for beta-peptides both by IR spectroscopic experiments and quantum chemical studies has significance for the better design and understanding of the backbone conformations of larger beta-peptides with biomedical potential. The present conformational analysis suggests that although beta-peptides, having an "extra" backbone torsion and hence more conformational freedom, should be more flexible than alpha-peptides, fewer backbone conformers are viable based on their relative energies. Thus, from a larger conformational arsenal, only a lower number of backbone conformers can emerge, which possibly had a fundamental effect on their applicability during prebiotic evolution.

Quantum chemistry in parallel with PQS.

This article describes the capabilities and performance of the latest release (version 4.0) of the Parallel Quantum Solutions (PQS) ab initio program package. The program was first released in 1998 and evolved from the TEXAS program package developed by Pulay and coworkers in the late 1970s. PQS was designed from the start to run on Linux-based clusters (which at the time were just becoming popular) with all major functionality being (a) fully parallel; and (b) capable of carrying out calculations on large-by ab initio standards-molecules, our initial aim being at least 100 atoms and 1000 basis functions with only modest memory requirements. With modern hardware and recent algorithmic developments, full accuracy, high-level calculations (DFT, MP2, CI, and Coupled-Cluster) can be performed on systems with up to several thousand basis functions on small (4-32 node) Linux clusters. We have also developed a graphical user interface with a model builder, job input preparation, parallel job submission, and post-job visualization and display.

A matrix isolation study on Ac-Gly-NHMe and Ac-l-Ala-NHMe, the simplest chiral and achiral building blocks of peptides and proteins.

The infrared absorption (IR) spectra of acetyl-N-methyl-glycine and acetyl-N-methyl-alanine have been recorded in dichloromethane and dimethyl sulfoxide-d(6) solution, as well as in Ar and Kr matrices. The spectra were assigned with the help of quantum chemical calculations. Based on the assignments of the matrix-isolation IR spectra, in line with theoretical predictions, two different hydrogen bonded conformers were identified, furthermore a third conformer is likely to be present, which cannot be unambiguously identified. In dichloromethane two conformers could be observed, while in dimethyl sulfoxide a single conformer could be identified. Vibrational circular dichroism (VCD) spectra of acetyl-N-methyl-l-alanine have also been recorded in solutions and matrices. These matrix-isolation VCD spectra not only support the assignments of the matrix-isolation IR spectra, but also demonstrate that these spectra can be interpreted much easier with the help of quantum chemical calculations than the VCD spectra recorded in solutions. It is also shown that the rotatory strength of some vibrational transitions changes rapidly as a function of the backbone torsional coordinates; hence the appearance of some regions in the VCD spectra is extremely dependent on any perturbations, e.g. weak intermolecular interactions.

Secondary structures of peptides and proteins via NMR chemical-shielding anisotropy (CSA) parameters.

Complete nuclear magnetic resonance (NMR) chemical-shielding tensors, sigma, have been computed at different levels of density-functional theory (DFT), within the gauge-including atomic orbital (GIAO) formalism, for the atoms of the peptide model For-L-Ala-NH2 as a function of the backbone dihedral angles phi and psi by employing a dense grid of 10 degrees. A complete set of rigorously orthogonal symmetric tensor invariants, {sigma iso, rho, tau}, is introduced, where sigma iso is the usual isotropic chemical shielding, while the newly introduced rho and tau parameters describe the magnitude and the orientation/shape of the chemical-shielding anisotropy (CSA), respectively. The set {sigma iso, rho, tau} is unaffected by unitary transformations of the symmetric part of the shielding tensor. The mathematically and physically motivated {rho, tau} anisotropy pair is easily connected to more traditional shielding anisotropy measures, like span (Omega) and skew (kappa). The effectiveness of the different partitions of the CSA information in predicting conformations of peptides and proteins has been tested throughout the Ramachandran space by generating theoretical NMR anisotropy surfaces for our For-L-Ala-NH2 model. The CSA surfaces, including Omega(phi, psi), kappa(phi, psi), rho(phi, psi), and tau(phi, psi) are highly structured. Individually, none of these surfaces is able to distinguish unequivocally between the alpha-helix and beta-strand secondary structural types of proteins. However, two- and three-dimensional correlated plots, including Omega versus kappa, rho versus tau, and sigma iso versus rho versus tau, especially for 13Calpha, have considerable promise in distinguishing among all four of the major secondary structural elements.

Conformers and photochemistry of propyl nitrites: a matrix isolation study.

The infrared spectra of both constitutional isomers (n and i) of propyl nitrite have been recorded in an Ar matrix. Conformational analysis and assignments of the vibrational transitions have been carried out on the basis of quantum chemical calculations. Assignment of spectral lines to different conformers was also aided experimentally, by utilizing the different rate of photodecomposition of the conformers, as well as by employing conformational cooling using a supersonic jet as the inlet source for matrix deposition. The rate of photodecomposition is primarily determined by the steric alignment of the nitrite group, whereas jet cooling affects mainly the conformation of the alkyl tail. On the basis of these experimental observations and computational predictions two to three conformers of isopropyl nitrite and eight conformers of n-propyl nitrite were identified. After broadband ultraviolet-visible (UV-vis) photolysis of isopropyl nitrite in the matrix, HNO, acetone, HNO.acetone complex, acetaldehyde, and nitrosomethane were identified as the main products. Furthermore, in a small amount, NO and possibly the isopropoxy radical were also present in the matrix. Photolysis of n-propyl nitrite yielded HNO, propanal, and their 1:1 complex as the main products together with a small amount of NO and cis-1-nitrosopropanol.