Charge-transfer transitions in protein circular dichroism spectra

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.

DichroWeb, a website for calculating protein secondary structure from circular dichroism spectroscopic data

Circular dichroism (CD) spectroscopy is a widely‐used method for characterizing the secondary structures of proteins. The well‐established and highly used analysis website, DichroWeb (located at: http://dichroweb.cryst.bbk.ac.uk/html/home.shtml) enables the facile quantitative determination of helix, sheet, and other secondary structure contents of proteins based on their CD spectra. DichroWeb includes a range of reference datasets and algorithms, plus graphical and quantitative methods for determining the quality of the analyses produced. This article describes the current website content, usage and accessibility, as well as the many upgraded features now present in this highly popular tool that was originally created nearly two decades ago.

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.

Accurate secondary structure prediction and fold recognition for circular dichroism spectroscopy

Circular dichroism (CD) spectroscopy is a widely used technique for the study of protein structure. Numerous algorithms have been developed for the estimation of the secondary structure composition from the CD spectra. These methods often fail to provide acceptable results on α/β-mixed or β-structure-rich proteins. The problem arises from the spectral diversity of β-structures, which has hitherto been considered as an intrinsic limitation of the technique. The predictions are less reliable for proteins of unusual β-structures such as membrane proteins, protein aggregates, and amyloid fibrils. Here, we show that the parallel/antiparallel orientation and the twisting of the β-sheets account for the observed spectral diversity. We have developed a method called β-structure selection (BeStSel) for the secondary structure estimation that takes into account the twist of β-structures. This method can reliably distinguish parallel and antiparallel β-sheets and accurately estimates the secondary structure for a broad range of proteins. Moreover, the secondary structure components applied by the method are characteristic to the protein fold, and thus the fold can be predicted to the level of topology in the CATH classification from a single CD spectrum. By constructing a web server, we offer a general tool for a quick and reliable structure analysis using conventional CD or synchrotron radiation CD (SRCD) spectroscopy for the protein science research community. The method is especially useful when X-ray or NMR techniques fail. Using BeStSel on data collected by SRCD spectroscopy, we investigated the structure of amyloid fibrils of various disease-related proteins and peptides.

Protein characterisation by synchrotron radiation circular dichroism spectroscopy

Circular dichroism (CD) spectroscopy is a well-established technique for the study of proteins. Synchrotron radiation circular dichroism (SRCD) spectroscopy extends the utility of conventional CD spectroscopy (i.e. using laboratory-based instruments) because the high light flux from a synchrotron enables collection of data to lower wavelengths, detection of spectra with higher signal-to-noise levels and measurements in the presence of strongly absorbing non-chiral components such as salts, buffers, lipids and detergents. This review describes developments in instrumentation, methodologies and bioinformatics that have enabled new applications of the SRCD technique for the study of proteins. It includes examples of the use of SRCD spectroscopy for providing static and dynamic structural information on molecules, including determinations of secondary structures of intact proteins and domains, assessment of protein stability, detection of conformational changes associated with ligand and drug binding, monitoring of environmental effects, examination of the processes of protein folding and membrane insertion, comparisons of mutant and modified proteins, identification of intermolecular interactions and complex formation, determination of the dispositions of proteins in membranes, identification of natively disordered proteins and their binding partners and examination of the carbohydrate components of glycoproteins. It also discusses how SRCD can be used in conjunction with macromolecular crystallography and other biophysical techniques to provide a more complete picture of protein structures and functions, including how proteins interact with other macromolecules and ligands. This review also includes a discussion of potential new applications in structural and functional genomics using SRCD spectroscopy and future instrumentation and bioinformatics developments that will enable such studies. Finally, the appendix describes a number of computational/bioinformatics resources for secondary structure analyses that take advantage of the improved data quality available from SRCD. In summary, this review discusses how SRCD can be used for a wide range of structural and functional studies of proteins.

Influence of the solvent on the self-assembly of a modified amyloid beta peptide fragment. II. NMR and computer simulation investigation

The conformation of a model peptide AAKLVFF based on a fragment of the amyloid beta peptide Abeta16-20, KLVFF, is investigated in methanol and water via solution NMR experiments and molecular dynamics computer simulations. In previous work, we have shown that AAKLVFF forms peptide nanotubes in methanol and twisted fibrils in water. Chemical shift measurements were used to investigate the solubility of the peptide as a function of concentration in methanol and water. This enabled the determination of critical aggregation concentrations. The solubility was lower in water. In dilute solution, diffusion coefficients revealed the presence of intermediate aggregates in concentrated solution, coexisting with NMR-silent larger aggregates, presumed to be beta-sheets. In water, diffusion coefficients did not change appreciably with concentration, indicating the presence mainly of monomers, coexisting with larger aggregates in more concentrated solution. Concentration-dependent chemical shift measurements indicated a folded conformation for the monomers/intermediate aggregates in dilute methanol, with unfolding at higher concentration. In water, an antiparallel arrangement of strands was indicated by certain ROESY peak correlations. The temperature-dependent solubility of AAKLVFF in methanol was well described by a van't Hoff analysis, providing a solubilization enthalpy and entropy. This pointed to the importance of solvophobic interactions in the self-assembly process. Molecular dynamics simulations constrained by NOE values from NMR suggested disordered reverse turn structures for the monomer, with an antiparallel twisted conformation for dimers. To model the beta-sheet structures formed at higher concentration, possible model arrangements of strands into beta-sheets with parallel and antiparallel configurations and different stacking sequences were used as the basis for MD simulations; two particular arrangements of antiparallel beta-sheets were found to be stable, one being linear and twisted and the other twisted in two directions. These structures were used to simulate circular dichroism spectra. The roles of aromatic stacking interactions and charge transfer effects were also examined. Simulated spectra were found to be similar to those observed experimentally (in water or methanol) which show a maximum at 215 or 218 nm due to pi-pi* interactions, when allowance is made for a 15-18 nm red-shift that may be due to light scattering effects.

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.

Interpretation of synchrotron radiation circular dichroism spectra of anionic, cationic, and zwitterionic dialanine forms

Electronic absorption and synchrotron radiation circular dichroism (SRCD) spectra of the anionic, cationic, and zwitterionic forms of L-alanyl-L-alanine (AA) in aqueous solutions were measured and interpreted by molecular dynamics (MD) and ab initio computations. Time-dependent density functional theory (TD DFT) was applied to predict the electronic excited states. The modeling enabled the assessment of the role of molecular conformation, charge, and interaction with the polar environment in the formation of the spectral shapes. Particularly, inclusion of explicit solvent molecules in the computations appeared to be imperative because of the participation of water orbitals in the amide electronic structure. Implicit dielectric continuum solvent models gave inferior results for clusters, especially at low-energy transitions. Because of the dispersion of transition energies, tens of water/AA clusters had to be averaged in order to obtain reasonable spectral shapes with a more realistic inhomogeneous broadening. The modeling explained most of the observed differences, as the anionic and zwitterionic SRCD spectra were similar and significantly different from the cationic spectrum. The greatest deviation between the experimental and theoretical curves observed for the lowest-energy negative anion signal can be explained by the limited precision of the TD DFT method, but also by the complex dynamics of the amine group. The results also indicate that differences in the experimental spectral shapes do not directly correlate with the peptide main-chain conformation. Future peptide and protein conformational studies based on circular dichroic spectroscopy can be reliable only if such effects of molecular dynamics, solvent structure, and polar solvent-solute interactions are taken into account.

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.

A reference database for circular dichroism spectroscopy covering fold and secondary structure space

MOTIVATION

Circular Dichroism (CD) spectroscopy is a long-established technique for studying protein secondary structures in solution. Empirical analyses of CD data rely on the availability of reference datasets comprised of far-UV CD spectra of proteins whose crystal structures have been determined. This article reports on the creation of a new reference dataset which effectively covers both secondary structure and fold space, and uses the higher information content available in synchrotron radiation circular dichroism (SRCD) spectra to more accurately predict secondary structure than has been possible with existing reference datasets. It also examines the effects of wavelength range, structural redundancy and different means of categorizing secondary structures on the accuracy of the analyses. In addition, it describes a novel use of hierarchical cluster analyses to identify protein relatedness based on spectral properties alone. The databases are shown to be applicable in both conventional CD and SRCD spectroscopic analyses of proteins. Hence, by combining new bioinformatics and biophysical methods, a database has been produced that should have wide applicability as a tool for structural molecular biology.

Geometry and Solvent Dependence of the Electronic Spectra of the Amide Group and Consequences for Peptide Circular Dichroism

The influence of geometry variations and solvent environment of N-methylacetamide on its energies and absorption intensities was systematically analyzed with the aid of the time-dependent density functional theory (TD DFT). Selective and often complicated reactions of individual electronic levels on the perturbations were found important for the resultant spectral profile. For example, the n-pi band position varied by tens of nanometers due to the C=O bond length oscillations, while it was rather unaffected by surrounding water. On the contrary, pi-pi type transition energies and intensities were broadly dispersed by the aqueous environment but exhibited a modest coordinate dependence. A simple electrostatic model used previously for absorption in the IR region (J. Chem. Phys. 2005, 122, 144501) explained these changes only partially. Additionally, electronic transfer between the solute and the solvent had to be considered for faithful modeling of the ultraviolet light absorption. The inclusion of the environment and dynamics in the modeling then provided more accurate positions, intensities, and realistic inhomogeneous widths of spectral lines. These factors were found important for absorption and circular dichroism spectra of larger peptides and proteins. This was demonstrated with a combined DFT/coupled oscillator model providing principal features observed in electronic circular dichroism spectra of standard peptide conformations.

First-principles calculations of protein circular dichroism in the far-ultraviolet and beyond

Understanding the relationship between the amino acid sequence of a protein and its unique, compact three-dimensional structure is one of the grand challenges in molecular biophysics. One exciting approach to the protein-folding problem is fast time-resolved spectroscopy in the ultra-violet (UV). Time-resolved electronic circular dichroism (CD) spectroscopy offers resolution on a nanosecond (or faster) timescale, but does not provide the spatial resolution of techniques like X-ray crystallography or NMR. There is a need to underpin fast timescale spectroscopic studies of protein folding with a stronger theoretical foundation. We review some recent studies in this regard and briefly highlight how modern quantum chemical models of aromatic groups have improved the accuracy of calculations of protein CD spectra near-UV. On the other side of the far-UV, we describe calculations indicating that charge-transfer transitions are likely to be responsible for bands observed in the vacuum UV in protein CD.

Synchrotron radiation circular dichroism spectroscopy of proteins and applications in structural and functional genomics

The technique of Synchrotron Radiation Circular Dichroism (SRCD) spectroscopy and its advantages over conventional circular dichroism spectroscopy are described in this tutorial review, as well as recent applications of the technique in structural and functional genomics. Circular dichroism (CD) spectroscopy is a well-established method in biological chemistry and structural biology, but its utility can be limited by the low flux of the light source in the far ultraviolet and vacuum ultraviolet wavelength regions in conventional CD instruments. The development of synchrotron radiation circular dichroism (SRCD), using the intense light of a synchrotron beam, has greatly expanded the utility of the method, especially as a tool for both structural and functional genomics. These applications take advantage of the enhanced features of SRCD relative to conventional CD: the ability to measure lower wavelength data containing more electronic transitions and hence more structural information, the higher signal-to-noise hence requiring smaller samples, the higher intensity enabling measurements in absorbing buffers and in the presence of lipids and detergents, and the ability to do faster measurements enabling high throughput and time-resolved spectroscopy.This article discusses recent developments in SRCD instrumentation, software, sample preparation and methods of analyses, with particular emphasis on their applications to the study of proteins. These advances have led to new applications in structural genomics (SG), including the potential for fold recognition as a means of target selection and the examination of membrane proteins, a class of proteins usually excluded from SG programmes. Other SG uses include detection of macromolecular interactions as a screen for complex formation, and examination of glycoproteins and sugar components. In functional genomics (FG) new applications include screening for ligand binding as a means of identifying function, and examination of structural differences in mutant proteins as a means of gaining insight into function.

Study of conformational properties of a biologically active peptide of fibronectin by circular dichroism, NMR and molecular dynamics simulation

Circular dichroism (CD), and NMR spectra have been recorded and molecular dynamics (MD) simulations have been performed in water and water-trifluoroethanol (TFE) mixed solvent for a synthetic biologically active 13-amino-acid fragment of human fibronectin and two related peptides. The CD results are interpreted on the basis of statistical analyses of MD trajectories and of ensuing calculations of CD spectra based on Schellman's matrix method. It is observed that the peptide conformation is quite variable in water and loses its mobility with the addition of TFE. (1)H-NOE data were found to be consistent with the most abundant calculated conformation.