Vladimír Baumruk

Vibrational Circular Dichroism of Proteins in H2O Solution

Spectroscopic studies have historically been important for estimation of the average secondary structure of proteins. Vibrational circular dichroism (VCD) is a hybrid of the more established electronic CD and infrared spectroscopies that combines the advantages of both (relatively high spectral resolution in IR region with high conformational sensitivity of CD methods). 1 For proteins the most valuable transition to study has proven to be the amide I band (assigned primarily to the C=O stretching of the amide group). Due to interference by water absorption band, all aqueous-phase VCD studies of the amide I have been done in D2O solution (termed as the amide I’). But the extent of deuteration is a poorly controlled parameter which leads to an ambiguity in interpretation of the data.

Anharmonic effects in IR, Raman, and Raman optical activity spectra of alanine and proline zwitterions

The difference spectroscopy of the Raman optical activity (ROA) provides extended information about molecular structure. However, interpretation of the spectra is based on complex and often inaccurate simulations. Previously, the authors attempted to make the calculations more robust by including the solvent and exploring the role of molecular flexibility for alanine and proline zwitterions. In the current study, they analyze the IR, Raman, and ROA spectra of these molecules with the emphasis on the force field modeling. Vibrational harmonic frequencies obtained with 25 ab initio methods are compared to experimental band positions. The role of anharmonic terms in the potential and intensity tensors is also systematically explored using the vibrational self-consistent field, vibrational configuration interaction (VCI), and degeneracy-corrected perturbation calculations. The harmonic approach appeared satisfactory for most of the lower-wavelength (200-1800 cm(-1)) vibrations. Modern generalized gradient approximation and hybrid density functionals, such as the common B3LYP method, provided a very good statistical agreement with the experiment. Although the inclusion of the anharmonic corrections still did not lead to complete agreement between the simulations and the experiment, occasional enhancements were achieved across the entire region of wave numbers. Not only the transitional frequencies of the C-H stretching modes were significantly improved but also Raman and ROA spectral profiles including N-H and C-H lower-frequency bending modes were more realistic after application of the VCI correction. A limited Boltzmann averaging for the lowest-frequency modes that could not be included directly in the anharmonic calculus provided a realistic inhomogeneous band broadening. The anharmonic parts of the intensity tensors (second dipole and polarizability derivatives) were found less important for the entire spectral profiles than the force field anharmonicities (third and fourth energy derivatives), except for a few weak combination bands which were dominated by the anharmonic tensor contributions.

Metalation of 5,10,15,20‐tetrakis(1‐methyl‐4‐pyridyl)porphyrin in silver colloids studied via time dependence of surface‐enhanced resonance Raman spectra

In SERS-active systems, when a free base porphyrin is adsorbed at a roughened metal surface, the metal ion may be incorporated into the porphyrin core. We have investigated the metalation kinetics of 5,10,15,20-tetrakis(1-methyl-4-pyridyl)porphyrin (TMPyP), reflected in the time evolution of its SERRS spectra. To establish proper metalation Raman markers, resonance Raman spectra of the free base and chemically Ag metalated TMPyP forms were measured under different conditions. The increase in intensity of the 395 cm-1 line, the decrease in intensity of the 329 cm-1 line, the reduction of the 1337+1360 cm-1 doublet to one strong line at 1340 cm-1 and the vanishing of the 1140 cm-1 medium line (the last exceptionally for 514.5 nm excitation) were found to be suitable markers for quantitative analysis. All other spectral changes related to metalation may coincide with effects caused by other variations of the porphyrin state. Time evolution of the porphyrin metalation was studied in various Ag colloid systems, with or without phosphate anions, citrate and/or Triton X-100. Time-dependent SERRS spectra of TMPyP were recorded within periods from minutes to several hours. Factor analysis of the SERRS spectra proved that the spectral changes were induced by only one phenomenon (metalation). The SERRS spectra were decomposed into the spectrum of the free base TMPyP and that of its Ag metalated form, and the metalation kinetics were determined by means of the time dependence of the metalated form portion in the SERRS spectrum. The results show a remarkable reliance on the state of the metal surface.

Comparison of Quantitative Conformer Analyses by Nuclear Magnetic Resonance and Raman Optical Activity Spectra for Model Dipeptides

Interpretation of the Raman optical activity (ROA) of peptides is difficult because of molecular flexibility and interaction with the solvent. Typically, simulations and experiments are compared in terms of a qualitative agreement between the spectra. However, on a series of the Pro-Gly, Gly-Pro, Pro-Ala, and Ala-Pro dipeptides more precise conformer ratios could be obtained with the aid of the density functional computations and numerical decomposition of the spectral shapes. All observed transitions were assigned, and the computed transition frequencies were scaled accordingly. Then the populations predicted by the optical spectroscopy agreed within a few percent with an analysis of the spin-spin coupling constants based on the Karplus equations, which was confirmed also by a comparison of calculated and experimental NMR couplings. The results are supported by molecular dynamics simulations and related to the previous conformational studies of similar molecules.

Interpretation of Raman and Raman Optical Activity Spectra of a Flexible Sugar Derivative, the Gluconic Acid Anion

Raman scattering and its polarized extension, Raman optical activity (ROA), are commonly used for monitoring of molecular conformational equilibria in solutions. This is complicated for saccharides due to extensive motions of the hydroxyl groups and other molecular parts. Standard interpretation procedures involving ab initio spectral simulations for a limited set of conformers are not adequate. In this study, a more general approach is proposed for the gluconic acid anion taken as a model compound, where quantum simulations of the spectra are directly coupled with molecular dynamics (MD) techniques. Such a multiscale approach reveals how the structural information is encoded in the broadened spectral lines. The spectra were simulated for solvent-solute clusters generated by MD. Conformational averaging was enabled by a limited library of conformers for which the spectral parameters could be calculated ab initio and moved on the MD geometries by Cartesian coordinate tensor transfer techniques. The B3LYP/CPCM/6-31+G** approximation was used as a default for computation of the source force fields and polarizability derivatives. The spectra thus obtained relatively faithfully reproduced most of the experimental features. The Amber and polarizable Amoeba MD force fields produced similar results; application of the latter, however, was limited by the long time necessary to achieve a converged conformational equilibrium. Both MD simulation and spectral averaging suggest that the hydroxyl groups as well as the backbone C-C bonds rotate relatively freely, with some restrictions in the vicinity of the carboxyl group. In spite of the averaging, spectral response of characteristic vibrational normal mode families, such as CH and OH bending, can clearly be identified in the spectra. The simulations thus confirm the experimental fact that flexible saccharides exhibit significant vibrational activity that reveals precious information about molecular structure and dynamics encoded in the Raman and ROA spectral shapes.

Side chain and flexibility contributions to the Raman optical activity spectra of a model cyclic hexapeptide

A model peptide, cyclo-(Phe-d-Pro-Gly-Arg-Gly-Asp), with a distinct folded structure containing short beta-hairpin and beta-sheet patterns was studied by Raman and Raman optical activity (ROA) spectroscopies. Unlike for previously analyzed vibrational circular dichroism of the same compound (Chirality 2008, 20, 1104), the Raman spectrum is dominated by side chain contributions and is more sensitive to their geometry fluctuations. The spectra and molecular motion were analyzed with the aid of the density functional theory simulations combined with molecular dynamics (MD). The side chain geometry fluctuations were found to significantly contribute to the broadening of the spectral bands, while dynamics of the backbone is rather restricted. According to our MD results, the side chains do not move freely but largely oscillate around preferred conformations. Averaging of computed spectra for many structures derived from the MD trajectories provided better spectral profiles than did a fixed geometry. The Raman and ROA scattering is dominated by the more polarizable phenylalanine and proline groups, as could be verified both by the computations and by comparison to experiments with a model Phe-d-Pro dipeptide. Computational analyses suggest that the ROA spectrum mostly senses local side chain conformation, whereas a vibrational coupling between different side chains contributes less. The coupling is mostly mediated by the peptide backbone and is restricted to specific vibrational region. The ROA spectroscopic technique thus provides important local structural information that needs, however, to be extracted by multiscale (QM/MM) simulation techniques.

Structure of the dimeric N-glycosylated form of fungal β-N-acetylhexosaminidase revealed by computer modeling, vibrational spectroscopy, and biochemical studies

BackgroundFungal β-N-acetylhexosaminidases catalyze the hydrolysis of chitobiose into its constituent monosaccharides. These enzymes are physiologically important during the life cycle of the fungus for the formation of septa, germ tubes and fruit-bodies. Crystal structures are known for two monomeric bacterial enzymes and the dimeric human lysosomal β-N-acetylhexosaminidase. The fungal β-N-acetylhexosaminidases are robust enzymes commonly used in chemoenzymatic syntheses of oligosaccharides. The enzyme from Aspergillus oryzae was purified and its sequence was determined.ResultsThe complete primary structure of the fungal β-N-acetylhexosaminidase from Aspergillus oryzae CCF1066 was used to construct molecular models of the catalytic subunit of the enzyme, the enzyme dimer, and the N-glycosylated dimer. Experimental data were obtained from infrared and Raman spectroscopy, and biochemical studies of the native and deglycosylated enzyme, and are in good agreement with the models. Enzyme deglycosylated under native conditions displays identical kinetic parameters but is significantly less stable in acidic conditions, consistent with model predictions. The molecular model of the deglycosylated enzyme was solvated and a molecular dynamics simulation was run over 20 ns. The molecular model is able to bind the natural substrate – chitobiose with a stable value of binding energy during the molecular dynamics simulation.ConclusionWhereas the intracellular bacterial β-N-acetylhexosaminidases are monomeric, the extracellular secreted enzymes of fungi and humans occur as dimers. Dimerization of the fungal β-N-acetylhexosaminidase appears to be a reversible process that is strictly pH dependent. Oligosaccharide moieties may also participate in the dimerization process that might represent a unique feature of the exclusively extracellular enzymes. Deglycosylation had only limited effect on enzyme activity, but it significantly affected enzyme stability in acidic conditions. Dimerization and N-glycosylation are the enzymes strategy for catalytic subunit stabilization. The disulfide bridge that connects Cys448 with Cys483 stabilizes a hinge region in a flexible loop close to the active site, which is an exclusive feature of the fungal enzymes, neither present in bacterial nor mammalian structures. This loop may play the role of a substrate binding site lid, anchored by a disulphide bridge that prevents the substrate binding site from being influenced by the flexible motion of the loop.

l-Alanyl-l-alanine Conformational Changes Induced by pH As Monitored by the Raman Optical Activity Spectra

Fine effects of the hydration, charge, and conformational structural changes in L-alanyl-L-alanine (Ala-Ala) dipeptide were studied with the aid of Raman and Raman optical activity (ROA) spectra. The spectra were recorded experimentally and analyzed by means of density functional computations. A (15)N and (13)C isotopically labeled analogue was synthesized and used to verify the vibrational mode assignment. Calculated shifts in vibrational frequencies for isotopically labeled molecule agreed well with the experiment. The assignment made it possible to scale computed vibrational frequencies and extract better structural information from the intensities. Solvent modeling with clusters obtained from molecular dynamics led to a qualitatively correct inhomogeneous broadening of Raman spectral lines but did not bring a convincing improvement of ROA signal when compared to a standard dielectric solvent correction. In comparison with the zwitterionic form, charged anionic and cationic dipeptides provided spectral variations that indicated different conformational behavior. Only minor backbone conformational change occurs in the cation, whereas the results indicate the presence of more anion conformers differing in the rotation of the NH(2) group and the backbone psi-angle. These findings are in agreement with previous electronic circular dichroism (ECD) and NMR studies. The results confirm the large potential of the ROA technique for the determination of final details in molecular structure and conformation.

Conformational properties of the Pro-Gly motif in the D-Ala-l-Pro-Gly-D-Ala model peptide explored by a statistical analysis of the NMR, Raman, and Raman optical activity spectra.

The Pro-Gly sequence in designed peptides and proteins is often used to mimic natural beta-hairpin turns. Shorter peptides containing this moiety, however, adopt multiple conformations, and their propensity to form the turn is not obvious. In this study, conformational flexibility of Pro-Gly was investigated with the aid of NMR, Raman scattering, and Raman optical activity (ROA). The spectra of a model tetrapeptide NH3+-D-Ala-L-Pro-Gly-D-Ala-CO2- were analyzed on the basis of statistical methods and density functional computations. Other peptide derivatives were adopted as supporting theoretical and experimental models. The results suggest that the loop conformation of the Pro-Gly core is not inherently stable in vacuum. On the other hand, in aqueous environment the propensity to form the beta-hairpin loops is an intrinsic property of the Pro-Gly sequence. It was observed also in a shorter Ac-D-Pro-Gly-NH-Me dipeptide. The attached alanine residues in the tetrapeptide stabilize the structure only partially. Thus an inclusion of the solvent in the calculations is important for correct description of peptide folding in the aqueous environment. The agreement of the optical spectra with the experiment was determined from overlaps between simulated and measured spectral curves. The comparison of the computed NMR, Raman, and ROA was hampered by experimental noise and limited accuracy of the computations. However, the statistical analysis of the spectroscopic data provided conformer distribution consistent with the computation of the relative energies. The combination of the NMR and Raman techniques with the quantum computations appeared very beneficial for the investigation of Pro-Gly conformational behavior, and can be recommended for future peptide folding studies.

METALATION OF POSITIVELY CHARGED WATER SOLUBLE MESOPORPHYRINS STUDIED VIA TIME-RESOLVED SERRS SPECTROSCOPY

Abstract Time-resolved SERRS spectra of 5,10,15,20-tetrakis[4-(trimethylammonio)phenyl]21 H ,23 H -porphine (TMAP) were recorded (using a multichannel Raman spectrometer) in various SERS-active Ag colloid/porphyrin systems. Data treatment based on a factor analysis was used to decompose all the SERRS spectra into two main components: SERRS spectrum of the free base TMAP and that of its Ag metalated form. The metalation kinetics obtained in this way was found to be highly dependent on the presence of phosphate anions, citrate and/or Triton X-100 in the colloidal system. The results are analogous to those previously obtained for 5,10,15,20-tetrakis(1-methyl-4-pyridyl)21 H ,23 H -porphine, a porphyrin with a substantially stronger tendency towards metalation.