Water-Mediated Electronic Structure of Oligopeptides Probed by Their UV Circular Dichroism, Absorption Spectra, and Time-Dependent DFT Calculations

Ordered Aggregates of Fmoc-Diphenylalanine at Alkaline pH as a Precursor of Fibril Formation and Peptide Gelation

The ultrashort peptide N-fluorenylmethoxycarbonyl-phenylalanyl-phenylalanine (FmocFF) has been largely investigated due to its ability to self-assemble into fibrils (100 nm-μm scale) that can form a sample-spanning gel network. The initiation of the gelation process requires either a solvent switch (water added to dimethyl sulfoxide) or a pH-switch (alkaline to neutral) protocol, both of which ensure the solubility of the peptide as a necessary step preceding gelation. While the respective gel phases are well understood in structural and material characteristics terms the pregelation conditions are known to a lesser extent. The question we asked is to what extent the gel-forming fibrils are already partially formed, i.e., oligomers or protofibrils. Focusing on the pregelation conditions for the pH-switch method, we investigated the self-assembly of soluble FmocFF aggregates in alkaline pH by UV circular dichroism, IR, vibrational circular dichroism, and 1H NMR spectroscopy for different peptide concentrations and more systematically as a function of temperature. The temperature dependence of the UVCD spectra of FmocFF in H2O and D2O revealed a complicated isotope effect that affects the peptide backbone and fluorene conformations in peptide aggregates differently. Moreover, we found that the melting of formed aggregates depends on peptide concentration in a nonmonotonic way. At 20 mM the UVCD data revealed the population of at least two different thermodynamic intermediate states, which seem to differ in terms of the relative arrangement of the fluorene moiety. The IR spectrum of this sample at room temperature indicates an antiparallel β-sheet arrangement, as suggested earlier in the literature. However, we show that this interpretation can only be valid if one invokes a nondispersive redshift of the two amide I' bands in a locally crystalline environment. The respective vibrational circular dichroism spectrum of the amide I' region is consistent with a left-handed helically twisted structure of the formed aggregates. A comparison of our data with spectra of the aqueous gel phase suggests that fibrils in the latter resemble the ones at alkaline pH probed by our experiments.

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.

A Joint Experimental and Theoretical Study on the Structural and Spectroscopic Properties of the Piv-Pro-d-Ser-NHMe Peptide

In this paper, we investigate the secondary structure of the Piv-Pro-d-Ser-NHMe peptide by means of nuclear magnetic resonance (NMR) and electronic circular dichroism (ECD) experiments, in conjunction with theoretical simulations based on molecular dynamics and time-dependent density functional theory calculations including polarizable embedding to account for solvent effects. The various experimental and theoretical protocols are assessed and validated, and are shown to provide a consistent description of the turn structure adopted by this peptide in solution. In addition, a simple fitting procedure is proposed to make the simulated and experimental ECD almost perfectly match. This full methodology is finally tested on another small peptide, enlightening its efficiency and robustness.

Theoretical-Computational Modeling of CD Spectra of Aqueous Monosaccharides by Means of Molecular Dynamics Simulations and Perturbed Matrix Method

The electronic circular dichroism (ECD) spectra of aqueous d-glucose and d-galactose were modeled using a theoretical-computational approach combining molecular dynamics (MD) simulations and perturbed matrix method (PMM) calculations, hereafter termed MD-PMM. The experimental spectra were reproduced with a satisfactory accuracy, confirming the good performances of MD-PMM in modeling different spectral features in complex atomic-molecular systems, as already reported in previous studies. The underlying strategy of the method was to perform a preliminary long timescale MD simulation of the chromophore followed by the extraction of the relevant conformations through essential dynamics analysis. On this (limited) number of relevant conformations, the ECD spectrum was calculated via the PMM approach. This study showed that MD-PMM was able to reproduce the essential features of the ECD spectrum (i.e., the position, the intensity, and the shape of the bands) of d-glucose and d-galactose while avoiding the rather computationally expensive aspects, which were demonstrated to be important for the final outcome, such as (i) the use of a large number of chromophore conformations; (ii) the inclusion of quantum vibronic coupling; and (iii) the inclusion of explicit solvent molecules interacting with the chromophore atoms within the chromophore itself (e.g., via hydrogen bonds).

Electronic circular dichroism from real-time propagation in state space

In this paper, we propose to compute the electronic circular dichroism (ECD) spectra of chiral molecules using a real-time propagation of the time-dependent Schrödinger equation (TDSE) in the space of electronic field-free eigenstates, by coupling TDSE with a given treatment of the electronic structure of the target. The time-dependent induced magnetic moment is used to compute the ECD spectrum from an explicit electric perturbation. The full matrix representing the transition magnetic moment in the space of electronic states is generated from that among pairs of molecular orbitals. In the present work, we show the ECD spectra of methyloxirane, of several conformers of L-alanine, and of the Λ-Co(acac)₃ complex, computed from a singly excited ansatz of time-dependent density functional theory eigenstates. The time-domain ECD spectra properly reproduce the frequency-domain ones obtained in the linear-response regime and quantitatively agree with the available experimental data. Moreover, the time-domain approach to ECD allows us to naturally go beyond the ground-state rotationally averaged ECD spectrum, which is the standard outcome of the linear-response theory, e.g., by computing the ECD spectra from electronic excited states.

A computational approach for modeling electronic circular dichroism of solvated chromophores

The present study consists in a novel computational protocol to model the UV-circular dichroism spectra of solvated species. It makes use of quantum-chemical calculations on a series of conformations of a flexible chromophore or on a series of chromophore/solvent clusters extracted from molecular dynamic simulations. The protocol is described and applied to the aqueous cationic tripeptide GAG⁺ and to the aqueous neutral decapeptide (GVGVP)₂ . The protocol has proven able to: (i) properly consider the conformational motion of solute in the given environment; (ii) give the actual statistical weight of each conformational state; (iii) provide a reliable quantum mechanical method able to reproduce the spectral features. Temperature effects on conformations and spectral properties are properly taken into account. The role of explicit solvent on the conformational analysis and the spectra calculation is discussed. The comparison of the calculated circular dichroism spectra with experimental ones recorded at different temperatures represents a strict validation test of the method.

Forbidden Secondary Structures Found in Gel-Forming Fibrils of Glycylphenylalanylglycine

The zwitterionic l-tripeptide glycylphenylalanylglycine self-assembles into very long crystalline fibrils in an aqueous solution, which causes the formation of an exceptionally strong gel phase (G' ∼ 5 × 10⁶ Pa). The Rietveld refinement analysis of its powder X-ray diffraction (PXRD) pattern reveals a unit cell with four peptides forming a P2₁2₁2₁ space group and adopting an inverse polyproline II conformation, that is, a right-handed helical structure that occupies the "forbidden" region of the Ramachandran plot. This unusual structure is stabilized by a plethora of intermolecular interactions facilitated by the large number of different functional groups of the unblocked tripeptide. Comparisons of simulated and experimental Fourier transform infrared and vibrational circular dichroism (VCD) amide I' profiles corroborate the PXRD structure. Our experimental setup reduces the sample to a quasi-two-dimensional network of fibrils. We exploited the influence of this reduced dimensionality on the amide I VCD to identify the main fibril axis. We demonstrate that PXRD, vibrational spectroscopy, and amide I simulations provide a powerful toolset for secondary structure and fibril axis determination.

Novel Isoindolinone-Based Analogs of the Natural Cyclic Peptide Fenestin A: Synthesis and Antitumor Activity

Small- and medium-sized cyclopeptides have been found to have extensive bioactivities and have drawn much attention from medicinal chemists. In the work described in this paper, various cyclic peptide analogs of Fenestin A were synthesized by intramolecular photoinduced electron-transfer cyclization reactions to study the influence of slight structural changes on the bioactivity of small cyclopeptides. The incorporation of thiazole and rigid isoindolinone fragments was found to improve the bioactivity of the cyclopeptide. Detailed in vitro studies of the apoptosis mechanism, mitochondrial membrane potential, cell cycle, intracellular Ca²⁺ concentration, and lactate dehydrogenase activity following treatment with a cyclopeptide showed that the cyclopeptide could induce apoptosis of tumor cells and lead to cycle arrest in the G2/M phase. The research also suggested that the photoinduced reaction could be applied to construct cyclic peptides stereoselectively, and the introduction of rigid fragments could enhance the biological activity of cyclopeptide drugs.

Far-UV circular dichroism signatures indicate fluorophore labeling induced conformational changes of penetratin

Fluorescent labeling is a broadly utilized approach to assess in vitro and in vivo behavior of biologically active, especially cell-penetrating and antimicrobial peptides. In this communication, far-UV circular dichroism (CD) spectra of penetratin (PEN) fluorophore conjugates reported previously have been re-evaluated. Compared to the intrinsically disordered native peptide, rhodamine B and carboxyfluorescein derivatives of free and membrane-bound PEN exhibit extrinsic CD features. Potential sources of these signals displayed above 220 nm are discussed suggesting the contributions of both intra- and intermolecular chiral exciton coupling mechanisms. Careful evaluation of the CD spectra of fluorophore-labeled peptides is a valuable tool for early detection of labeling-provoked structural alterations which in turn may modify the membrane binding and cellular uptake compared to the unconjugated form.

Remarkably high solvatochromism in the circular dichroism spectra of the polyproline-II conformation: limitations or new opportunities?

Circular dichroism is a conventional method for studying the secondary structures of peptides and proteins and their transitions. While certain circular dichroism features are characteristic of α-helices and β-strands, the third most abundant secondary structure, the polyproline-II helix, does not exhibit a strictly conserved spectroscopic appearance. Due to its extended nature, the polyproline-II helix is highly accessible to the surrounding solvent; thus, the environment has a critical influence on the lineshape of the circular dichroism spectra of this structure. To showcase possible effects due to the medium, in this work, we report an experimental spectroscopic study of polyproline-II-forming oligomeric peptides in various environments: solvents, detergent micelles, and liposomes. Strikingly, the examination of an oligomeric peptide in a solvent series showed a remarkable 7 nm solvatochromic shift in the main negative band starting with hexafluoropropan-2-ol and moving to hexane. Furthermore, a previously predicted positive band below 200 nm was discovered in the spectra in nonpolar environments. In isotropic liposomes, the expected transition to the transmembrane state correlated with the appearance of a positive band at 228 nm. Our results demonstrate that changes in solvation should be taken into consideration when assessing the circular dichroism spectra of peptides expected to adopt the polyproline-II conformation. Although this precaution may complicate spectral analysis, characterization of solvent-induced spectral changes can generate new opportunities for testing the location of peptides in complex systems such as micelles or lipid bilayers.

Short peptides as predictors for the structure of polyarginine sequences in disordered proteins

Intrinsically disordered proteins and intrinsically disordered regions are frequently enriched in charged amino acids. Intrinsically disordered regions are regularly involved in important biological processes in which one or more charged residues is the driving force behind a protein-biomolecule interaction. Several lines of experimental and computational evidence suggest that polypeptides and proteins that carry high net charges have a high preference for extended conformations with average end-to-end distances exceeding expectations for self-avoiding random coils. Here, we show that charged arginine residues even in short glycine-capped model peptides (GRRG and GRRRG) significantly affect the conformational propensities of each other when compared with the intrinsic propensities of a mostly unperturbed arginine in the tripeptide GRG. A conformational analysis based on experimentally determined J-coupling constants from heteronuclear NMR spectroscopy and amide I' band profiles from vibrational spectroscopy reveals that nearest-neighbor interactions stabilize extended β-strand conformations at the expense of polyproline II and turn conformations. The results from molecular dynamics simulations with a CHARMM36m force field and TIP3P water reproduce our results only to a limited extent. The use of the Ramachandran distribution of the central residue of GRRRG in a calculation of end-to-end distances of polyarginines of different length yielded the expected power law behavior. The scaling coefficient of 0.66 suggests that such peptides would be more extended than predicted by a self-avoiding random walk. Our findings thus support in principle theoretical predictions.