Charge Transfer Transitions in Neutral and Ionic Polypeptides: A Theoretical Study

An Improved Diabatization Scheme for Computing the Electronic Circular Dichroism of Proteins

We advance the quality of first-principles calculations of protein electronic circular dichroism (CD) through an amelioration of a key deficiency of a previous procedure that involved diabatization of electronic states on the amide chromophore (to obtain interamide couplings) in a β-strand conformation of a diamide. This yields substantially improved calculated far-ultraviolet (far-UV) electronic circular dichroism (CD) spectra for β-sheet conformations. The interamide couplings from the diabatization procedure for 13 secondary structural elements (13 diamide structures) are applied to compute the CD spectra for seven example proteins: myoglobin (α helix), jacalin (β strand), concanavalin A (β type I), elastase (β type II), papain (α + β), 310-helix bundle (310-helix) and snow flea antifreeze protein (polyproline). In all cases, except concanavalin A and papain, the CD spectra computed using the interamide couplings from the diabatization procedure yield improved agreement with experiment with respect to previous first-principles calculations.

UV-VUV Photofragmentation Spectroscopy of Isolated Neutral Fragile Macromolecules: A Proof-of-Principle Based on a Deprotonated Vancomycin-Peptide Noncovalent Complex

The gas phase offers the possibility to analyze organic molecules by ultraviolet-vacuum ultraviolet (UV-VUV) spectroscopy without any solvent effect or limitation in terms of spectral range due to absorption by the solvent. Up to now, the size and chemical composition of neutral molecular systems under study have been limited by the use of vaporization methods based on thermal heating. Soft sources of gas-phase thermolabile molecular systems such as electrospray or matrix-assisted laser desorption ionization are appealing alternatives to heating-based techniques, but they lead to the production of ions. In such cases, UV-VUV action spectroscopy is then the method of choice to study the electronic structure and corresponding photodynamics of these gas-phase molecular ions. However, previous investigations have shown that the UV-VUV action spectrum of a given molecular ion depends on the charge state, which in many cases might be a caveat. Here, by means of synchrotron radiation coupled to mass spectrometry and through the test case of the glycopeptide antibiotic vancomycin noncovalently bound to a deprotonated small peptide, we show that the UV-VUV photofragmentation spectrum of neutral thermally fragile organic molecules can be obtained via charge-tagging action spectroscopy.

UV-Visible Lysine-Glutamate Dimer Excitations in Protein Charge Transfer Spectra: TDDFT Descriptions Using an Optimally Tuned CAM-B3LYP Functional

Recent reports of distinctive UV-vis absorption profiles for monomeric proteins rich in charged amino acids that span 250-800 nm have opened up a new label-free optical spectral window for probing biomolecular structure and interactions. Combined experimental-computational studies have revealed that such broad absorption profiles of these proteins arise from photoexcited charge transfer (CT) transitions in spatially proximal charged amino acids such as lysine (Lys) and glutamate (Glu). Here, using time-dependent density functional theory (TDDFT) with an optimally tuned CAM-B3LYP functional, we refine the computed UV-vis spectra for Lys-Glu dimers within protein folds and quantify the percentage CT character of the constituent transitions. The optimally tuned functionals are derived through a careful analysis of the CAM-B3LYP parameter space for Lys-Glu dimers as a function of amino-acid conformation and side chain separation. Our studies reveal that the tuned Lys-Glu dimer spectrum spans 150-650 nm and possesses 5 specific types of CT excitations with diverse and large spatial charge separation length scales of 2-10 Å. These include inter-/intra-residue peptide backbone to peptide backbone (BB-CT) excitations spanning 160-210 nm, inter-/intra-residue peptide backbone to side chain (BS-CT) excitations spanning 160-260 nm, and side chain to side chain (SS-CT) excitations, which show the broadest absorption range spanning 260-650 nm.

Optical backbone-sidechain charge transfer transitions in proteins sensitive to secondary structure and modifications

We report the UV-vis absorption originating from protein backbone-sidechain charge transfer transitions in charged amino acids with anionic and cationic sidechains.

A model for ultra-fast charge transport in membrane proteins

Isolated proteins have recently been observed to transport charge and reactivity over very long distances with extraordinary rates and near perfect efficiencies in spite of their site. This is not the case if the peptide is in water, where the efficiency of charge hopping to the next site is reduced to approximately 2%. Here, water is not an ideal solvent for charge transport. The issue at hand is how to explain such enormous charge transfer quenching in water compared to another typical medium, namely lipid. We performed molecular dynamics simulations to computationally substantiate the novel long-distance charge transfer yield of the polypeptides in lipids. This is characterized by the charge transfer persistent-distance decay constant and not by the rate, which is seldom, if ever, measured and hence not directly addressed here. This model can encompass an extremely wide range of yields over very long distances in peptides in various media. The calculations here demonstrate the good charge transport efficiency in lipids in contrast to the poor efficiency in water. The protein charge transport also exhibits a very strong anisotropic effect in lipids. The peptide secondary structure effect of charge transfer in membranes is analyzed in contrast to that in water. These results suggest that this model can be useful for the prediction of charge transfer efficiency in various environments of interest and indicate that the charge transfer is highly efficient in membrane proteins.

Action spectroscopy of a protonated peptide in the ultraviolet range

Action spectroscopy of substance P, a model undecapeptide, has been probed from 5.2 eV to 20 eV.

Contribution of synchrotron radiation to photoactivation studies of biomolecular ions in the gas phase

Photon activation of ions in the visible and ultraviolet range attracts a growing interest, partly for its promising applications in tandem mass spectrometry. However, this task is not trivial, as it requires notably high brilliance photon sources. Hence, most of the work in that field has been performed using lasers. Synchrotron radiation is a source continuously tunable over a wide photon energy range and which possesses the necessary characteristics for ion activation. This review focuses on the array of applications of synchrotron radiation in photon activation of ions ranging from near UV to soft X-rays.

First principles study of photostability within hydrogen-bonded amino acids

The photochemistry and photophysics of a two-glycine minimal model is studied at the CASPT2//CASSCF level of theory. Different photoinduced processes are discussed, on the basis of the calculated minimum energy paths and the characterization of the electronic state crossings. Two main processes could provide UV-photostability to the hydrogen-bonded peptide system: (i) forward-backward photoinduced electron/proton transfer involving the H in the hydrogen bond, (ii) singlet-singlet energy transfer between two amino acids, providing ultrafast population of the low-energy n,π* state.

Modeling the infrared and circular dichroism spectroscopy of a bridged cyclic diamide

Density functional theory based molecular dynamics simulations are used to study the structure, infrared (IR) spectroscopy, circular dichroism (CD) spectroscopy, and coupling between the amide I vibrations of a bridged cyclic diamide in the gas phase and in aqueous solution. IR spectra computed via the dipole moment time correlation function show a large red-shift of 30 cm(-1) in the amide I vibration in solution compared to the gas phase, and are in good agreement with experiment. Conformationally averaged CD spectra computed using the CIS(D) method are highly sensitive to the structures used, and structures sampled in the aqueous phase simulation are required to obtain qualitatively correct CD spectra. Analysis of the coupling between the amide I modes shows that in the aqueous phase there is an increased localization of the vibrations on the individual peptide groups and a reduction in the mode coupling parameter compared to the gas phase. Overall, the results illustrate the significance of incorporating molecular dynamics in the simulation of IR and CD spectra.

UV resonance raman investigation of electronic transitions in alpha-helical and polyproline II-like conformations

UV resonance Raman (UVRR) excitation profiles and Raman depolarization ratios were measured for a 21-residue predominantly alanine peptide, AAAAA(AAARA) 3A (AP), excited between 194 and 218 nm. Excitation within the pi-->pi* electronic transitions of the amide group results in UVRR spectra dominated by amide vibrations. The Raman cross sections and excitation profiles provide information about the nature of the electronic transitions of the alpha-helix and polyproline II (PPII)-like peptide conformations. AP is known to be predominantly alpha-helical at low temperatures and to take on a PPII helix-like conformation at high temperatures. The PPII-like and alpha-helix conformations show distinctly different Raman excitation profiles. The PPII-like conformation cross sections are approximately twice those of the alpha-helix. This is due to hypochromism that results from excitonic interactions between the NV 1 transition of one amide group with higher energy electronic transitions of other amide groups, which decreases the alpha-helical NV 1 (pi-->pi*) oscillator strengths. Excitation profiles of the alpha-helix and PPII-like conformations indicate that the highest signal-to-noise Raman spectra of alpha-helix and PPII-like conformations are obtained at excitation wavelengths of 194 and 198 nm, respectively. We also see evidence of at least two electronic transitions underlying the Raman excitation profiles of both the alpha-helical and the PPII-like conformations. In addition to the well-known approximately 190 nm pi-->pi* transitions, the Raman excitation profiles and Raman depolarization ratio measurements show features between 205-207 nm, which in the alpha-helix likely results from the parallel excitonic component. The PPII-like helix appears to also undergo excitonic splitting of its pi-->pi* transition which leads to a 207 nm feature.

Characterization of dry globular proteins and protein fibrils by synchrotron radiation vacuum UV circular dichroism

Circular dichroism using synchrotron radiation (SRCD) can extend the spectral range down to approximately 130 nm for dry proteins, potentially providing new structural information. Using a selection of dried model proteins, including alpha-helical, beta-sheet, and mixed-structure proteins, we observe a low-wavelength band in the range 130-160 nm, whose intensity and peak position is sensitive to the secondary structure of the protein and may also reflect changes in super-secondary structure. This band has previously been observed for peptides but not for globular proteins, and is compatible with previously published theoretical calculations related to pi-orbital transitions. We also show that drying does not lead to large changes in the secondary structure and does not induce orientational artifacts. In combination with principal component analysis, our SRCD data allow us to distinguish between two different types of protein fibrils, highlighting that bona fide fibrils formed by lysozyme are structurally more similar to the nonclassical fibrillar aggregates formed by the SerADan peptide than with the amyloid formed by alpha-synuclein. Thus, despite the lack of direct structural conclusions, a comprehensive SRCD-based database of dried protein spectra may provide a useful method to differentiate between various types of supersecondary structure and aggregated protein species.

Charge-transfer transitions in protein circular dichroism calculations

Charge-transfer transitions in proteins play a key role in many biophysical processes, from the behavior of redox proteins to photochemical reactions. We present ab initio calculations on a model dipeptide and more approximate calculations of the electronic excited states of proteins which, taken together, provide the most definitive assignment and characterization of charge-transfer transitions in proteins to date. We have calculated from first principles the electronic circular dichroism (CD) spectra of 31 proteins on the basis of their structures. Compared to previous studies, we achieve more accurate calculated CD spectra between 170 and 190 nm, owing mainly to the importance in alpha-helices of a charge-transfer transition from the lone pair on one peptide group to the pi* orbital on the next peptide group.

Circular and linear dichroism of proteins

Circular dichroism (CD) is an important technique in the structural characterisation of proteins, and especially for secondary structure determination. The CD of proteins can be calculated from first principles using the so-called matrix method, with an accuracy which is almost quantitative for helical proteins. Thus, for proteins of unknown structure, CD calculations and experimental data can be used in conjunction to aid structure analysis. Linear dichroism (LD) can be calculated using analogous methodology and has been used to establish the relative orientations of subunits in proteins and protein orientation in an environment such as a membrane. However, simple analysis of LD data is not possible, due to overlapping transitions. So coupling the calculations and experiment is an important strategy. In this paper, the use of LD for the determination of protein orientation and how these data can be interpreted with the aid of calculations, are discussed. We review methods for the calculation of CD spectra, focusing on semiempirical and ab initio parameter sets used in the matrix method. Lastly, a new web interface for online CD and LD calculation is presented.

Theoretical investigation of the photoinitiated folding of HP-36

A computational model was developed to examine the phototriggered folding of a caged protein, a protein modified with an organic photolabile cross-linker. Molecular dynamics simulations of the modified 36-residue fragment of subdomain B of chicken villin head piece with a photolabile linker were performed, starting from both the caged and the uncaged structures. Construction of a free-energy landscape, based on principal components as well as on radius of gyration versus root-mean-square deviation, and circular dichroism calculations were employed to characterize folding behavior and structures. The folded structures observed in the molecular dynamics trajectories were found to be similar to that of the wild-type protein, in agreement with the published experimental results. The free-energy landscapes of the modified and wild-type proteins have similar topology, suggesting common thermodynamic/kinetic behavior. The existence of small differences in the free-energy surface of the modified protein from that of the native protein, however, indicates subtle differences in the folding behavior.

Environment-Controlled Interchromophore Charge Transfer Transitions in Dipeptides Probed by UV Absorption and Electronic Circular Dichroism Spectroscopy

Charge transfer (CT) transitions between the C-terminal carboxylate and peptide group have been investigated for alanyl-X and X-alanine dipeptides by far-UV absorption and electronic circular dichroism (ECD) spectroscopy (where X represents different amino acid residues). The spectra used in the present study were obtained by subtracting the spectrum of the cationic species from that of the corresponding zwitterionic peptide spectrum. These spectra displayed three bands, e.g., band I between 44 and 50 kK (kK = 10(3) cm(-1)), band II at 53 kK, and band III above 55 kK, which were, respectively, assigned to a n(COO-) --> pi* CT transition, a pi(COO-) --> pi* CT transition, and a carboxylate pi --> pi* (NV1) transition, respectively By comparison of the intensity, bandwidth, and wavenumber position of band I of some of the investigated dipeptides, we found that positive charges on the N-terminal side chain (for X = K), and to a minor extent also the N-terminal proton, reduce its intensity. This can be understood in terms of attractive Coulomb interactions that stabilize the ground state over the charge transfer state. For alanylphenylalanine, we assigned band I to a n(COO-) --> pi* CT transition into the aromatic side chain, indicating that aromatic side chains interact electronically with the backbone. We also performed ECD measurements at different pH values (pH 1-6) for a selected subset of XA and AX peptides. By subtraction of the pH 1 spectrum from that observed at pH 6, the ECD spectrum of the CT transition was obtained. A titration curve of their spectra reveals a substantial dependence on the protonation state of the aspartic acid side chain of AD, which is absent in DA and AE. This most likely reflects a conformational transition of the C-terminus into a less extended state, though the involvement of a side chain --> peptide CT transition cannot be completely ruled out.

Dissociation Kinetics of Peptide Ions

The dissociation of peptide ions has been found to have ultrafast components that in many ways are uniquely different from typical unimolecular kinetics. As such, some peptide reactions provide new channels, which do not conform to statistical models of reaction kinetics. When the dissociation rates are in the 100 fs range, they are in a time scale where statistical methods do not yet apply, although molecules that have not yet dissociated will later in time undergo statistical redistribution of their excess energy, which, however, may not lead to noticeable reactivity within the experimental time frames for large peptides and hence are simply dissipative. This work is meant to reconcile the long time statistical results of Lifshitz et al. (2003) with the work of Schlag et al. (1995/6) that suggests an alternate parallel and much faster time scale for dissociation. It is argued that the two sets of results and interpretations augment one another and in fact open up a most interesting new field of peptide kinetics in addition to the unimolecular behavior, which becomes de facto arrested by the shear size of the molecule being unable to find a transition state on any reasonable time scale.

Recent developments in the electronic spectroscopy of amides and alpha-helical polypeptides

Recent experimental and theoretical advances in understanding the electronic excited states of simple amides are reviewed. Polarized reflection spectroscopy of single crystals of N-acetylglycine shows that the direction of the first pipi* (NV(1)) transition dipole moment of a secondary amide differs by approximately 15 degrees from that of a primary amide. Ab initio calculations on simple amides support this conclusion. Ab initio studies of di- and tri-amides demonstrate that several inter-amide charge-transfer (CT) transitions occur in the 150-175-nm region, between the NV(1) and NV(2) transitions. When the correct dipole transition moment direction for peptides is used in calculations of the circular dichroism of the alpha-helix, the results are much improved over those from earlier calculations that used the direction for primary amides. Studies that consider the mixing of the NV(1) transition with CT transitions are reviewed. These indicate that such mixing is likely to have a significant effect on the absorption and CD spectra of the alpha-helix and other types of peptide conformation. Nevertheless, the independent systems model gives a reasonable first approximation to the absorption and CD spectra of the alpha-helix.