Kwang Im Oh

Nitrile and thiocyanate IR probes: Quantum chemistry calculation studies and multivariate least-square fitting analysis

Hydration effects on the C[Triple Bond]N stretching mode frequencies of MeCN and MeSCN are investigated by carrying out ab initio calculations for a number of MeCN-water and MeSCN-water complexes with varying number of water molecules. It is found that the CN frequency shift induced by the hydrogen-bonding interactions with water molecules originate from two different ways to form hydrogen bonds with the nitrogen atom of the CN group. Considering the MeCN- and MeSCN-water cluster calculation results as databases, we first examined the validity of vibrational Stark effect relationship between the CN frequency and the electric field component parallel to the CN bond and found no strong correlation between the two. However, taking into account of additional electric field vector components is a simple way to generalize the vibrational Stark theory for the nitrile chromophore. Also, the electrostatic potential calculation method has been proposed and examined in detail. It turned out that the interactions of water molecules with nitrogen atoms lone pair orbital and with nitrile pi orbitals can be well described by the electrostatic potential calculation method. The present computational results will be of use to quantitatively simulate various linear and nonlinear vibrational spectra of nitrile compounds in solutions.

Nitrile and thiocyanate IR probes: Molecular dynamics simulation studies

Nitrile- and thiocyanate-derivatized amino acids have been found to be useful IR probes for investigating their local electrostatic environments in proteins. To shed light on the CN stretch frequency shift and spectral lineshape change induced by interactions with hydrogen-bonding solvent molecules, we carried out both classical and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations for MeCN and MeSCN in water. These QM/MM and conventional force field MD simulation results were found to be inconsistent with the experimental results as well as with the high-level ab initio calculation results of MeCN-water and MeSCN-water potential energies. Thus, a new set of atomic partial charges of MeCN and MeSCN is obtained. By using the MD simulation trajectories and the electrostatic potential model recently developed, the CN and SCN stretching mode frequency trajectories were obtained and used to simulate the IR spectra. The C[Triple Bond]N frequency blueshifts of MeCN and MeSCN in water are estimated to be 9.0 and 1.9 cm(-1), respectively, in comparison with those of gas phase values. These values are found to be in reasonable agreement with the experimentally measured IR spectra of MeCN, MeSCN, beta-cyano-L-alanine, and cyanylated cysteine in water and other polar solvents.

β-Azidoalanine as an IR Probe : Application to Amyloid Aβ(16-22) Aggregation

Beta-azidoalanine dipeptide 1 was synthesized, and its azido stretching vibration in H2O and dimethyl sulfoxide (DMSO) was studied by using Fourier transform (FT) IR spectroscopy. The dipole strength of the azido stretch mode is found to be about 19 and 5 times larger than those of the CN and SCN stretch modes, respectively, which have been used as local environmental IR sensors. The azido stretch band in H2O is blue-shifted by about 14 cm(-1) in comparison to that in DMSO, indicative of its sensitivity to the electrostatic environment. To test the utility of beta-azidoalanine as an IR probe of the local electrostatic environment in proteins, azidopeptide 4 was prepared by its incorporation into Abeta(16-22) peptide of the Alzheimers disease amyloid beta-protein at position Ala21. The amide I IR spectrum of 4 in D2O suggests that the azidopeptide thus modified forms in-register beta-sheets in aggregates as observed for normal Abeta(16-22). The azido peak frequency of 4 in aggregates is almost identical to that in DMSO, indicating that the azido group is not exposed to water but to the hydrophobic environment. We believe that beta-azidoalanine will be used as an effective IR probe for providing site-specific information about the local electrostatic environments of proteins.

Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: comparisons with IR and vibrational circular dichroism spectra.

We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed alpha-helix and extended beta and P(II) conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the P(II) conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is alpha(R). Polarizable CHARMM force field gives almost exclusively P(II) conformation and other force fields predict that both alpha-helical and extended (beta and P(II)) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the P(II) conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.

Circular dichroism eigenspectra of polyproline II and β-strand conformers of trialanine in water: Singular value decomposition analysis†

Despite that a number of experimental and theoretical investigations have been carried out to determine the structure of trialanine in water, the reported populations of polyproline II (PPII) and β-strand conformers vary and were found to be dependent on which spectroscopic method was used. Such discrepancies are due to limitations of different spectroscopic methods used. Here, the temperature- and pH-dependent circular dichroism (CD) and NMR experiments have been carried out to develop a self-consistent singular value decomposition procedure. The temperature-dependent CD spectra indicate the presence of two conformers, but due to the two peptide bonds in a trialanine, one should take into consideration of four different conformers to fully interpret the NMR results. From the pH-dependent NMR coupling constant measurements, the conformation of zwitterionic trialanine is little different from that of cationic one. The strong pH dependency of CD spectrum is likely due to charge transfer transitions between carboxylate and nearby peptide groups or internal field effects not to pH-dependent conformational change. To simultaneously analyze the temperature-dependent CD and NMR data, a self-consistent procedure was used to newly determine the reference NMR coupling constants required to estimate one of the peptide dihedral angles. From the estimated enthalpy and entropy changes associated with the transition from enthalpically favorable PPII conformer to entropically favorable β-strand conformer, the relative populations of the four possible conformers of trialanine were determined and compared with the previous experimental findings. We anticipate that the present experimental results and interpretation procedure would be of use in determining the solution structures of small oligopeptides in the future.

A comprehensive library of blocked dipeptides reveals intrinsic backbone conformational propensities of unfolded proteins

Despite prolonged scientific efforts to elucidate the intrinsic peptide backbone preferences of amino‐acids based on understanding of intermolecular forces, many open questions remain, particularly concerning neighboring peptide interaction effects on the backbone conformational distribution of short peptides and unfolded proteins. Here, we show that spectroscopic studies of a complete library of 400 dipeptides reveal that, irrespective of side‐chain properties, the backbone conformation distribution is narrow and they adopt polyproline II and β‐strand, indicating the importance of backbone peptide solvation and electronic effects. By directly comparing the dipeptide circular dichroism and NMR results with those of unfolded proteins, the comprehensive dipeptides form a complete set of structural motifs of unfolded proteins. We thus anticipate that the present dipeptide library with spectroscopic data can serve as a useful database for understanding the nature of unfolded protein structures and for further refinements of molecular mechanical parameters. Proteins 2011;

Azido Gauche Effect on the Backbone Conformation of β-Azidoalanine Peptides

To study the azido gauche effect on the backbone conformation of β-azidoalanine (Aza) dipeptide (AAD, Ac-Aza-NHMe) and tripeptide (AAT, Ac-Aza-Aza-NH(2)), we used spectroscopic methods in combination with quantum chemistry calculations and molecular dynamics (MD) simulations. From the (1)H NMR coupling constants and (1)H,(1)H NOESY experimental data, we found that AAD in water mainly adopts a seven-membered cyclic (C(7)) rather than polyproline II (P(II)) backbone conformation and prefers the gauche- (g(-)) side-chain conformer. From the amide I IR absorption and circular dichroism (CD) spectra, the backbone conformation of AAD in water is found to deviate from P(II) but is rather close to C(7). Thus, the backbone conformation of AAD differs from that of alanine dipeptide (AD, Ac-Ala-NHMe), which is mainly P(II) in water. The underlying origin of the backbone conformational difference between AAD and AD in water was elucidated by quantum chemistry calculations with density functional theory (DFT). It was found that the C(7)/g(-) conformer is the lowest energy structure of an isolated AAD. Here, the β-azido group forms intramolecular electrostatic interactions with two neighboring peptide bonds, which are facilitated by the azido gauche effect. Thus, the β-azido group appears to be responsible for directing the peptide backbone conformation toward the C(7) structure. The quantum mechanical/molecular mechanical (QM/MM) MD simulations show that AAD in water adopts neither P(II) nor right-handed α-helix (α(R)) and prefers the g(-) conformer. Thus, the intramolecular electrostatic interactions between the β-azido group and two nearby peptide bonds are also found even in the aqueous solution structure of AAD. Consequently, the β-azido group appears to be an effective C(7)-conformation-directing element, which may also be useful for tuning the structures of other amino acids and polypeptides.

Conformational distributions of denatured and unstructured proteins are similar to those of 20 × 20 blocked dipeptides

Understanding intrinsic conformational preferences of amino-acids in unfolded proteins is important for elucidating the underlying principles of their stability and re-folding on biological timescales. Here, to investigate the neighbor interaction effects on the conformational propensities of amino-acids, we carried out 1H NMR experiments for a comprehensive set of blocked dipeptides and measured the scalar coupling constants between alpha protons and amide protons as well as their chemical shifts. Detailed inspection of these NMR properties shows that, irrespective of amino-acid side-chain properties, the distributions of the measured coupling constants and chemical shifts of the dipeptides are comparatively narrow, indicating small variances of their conformation distributions. They are further compared with those of blocked amino-acids (Ac–X–NHMe), oligopeptides (Ac–GGXGG–NH2), and native (lysozyme), denatured (lysozyme and outer membrane protein X from Escherichia coli), unstructured (Domain 2 of the protein 5A of Hepatitis C virus), and intrinsically disordered (hNlg3cyt: intracellular domain of human NL3) proteins. These comparative investigations suggest that the conformational preferences and local solvation environments of the blocked dipeptides are quite similar to not only those of other short oligopeptides but also those of denatured and natively unfolded proteins.

Neighboring Residue Effects in Terminally Blocked Dipeptides: Implications for Residual Secondary Structures in Intrinsically Unfolded/Disordered Proteins

For nuclear magnetic resonance (NMR)-based protein structure determinations, the random coil chemical shifts are very important because the secondary and tertiary protein structure predictions become possible by examining deviations of measured chemical shifts from those reference chemical shift values. In addition, neighboring residue effects on chemical shifts and J-coupling constants are crucial in understanding the nature of conformational propensities exhibited by unfolded or intrinsically disordered proteins. We recently reported the 1D NMR results for a complete set of terminally blocked dipeptides (Oh KI, Jung YS, Hwang GS, Cho M. J Biomol NMR 2012;53:25-41), but the NMR resonance assignments were not possible so that the average chemical shifts and J-coupling constants were only considered. In the present work, to thoroughly investigate the neighboring residue effects and random coil chemical shifts we extend the previous studies with 2D NMR, and measured all the (3) J(HNHα) values and H(α) and H(N) chemical shifts of the same set of terminally blocked dipeptides that are free from structural effects like secondary structure, hydrogen-bond, long-range backbone, and side-chain interactions. In particular, the preceding and following residue effects on amino-acid backbone conformational propensities are revealed and directly compared with previous works on either short peptides or empirical chemical shift database.