First-principles calculations of protein circular dichroism in the near ultraviolet

Electronic transitions in aromatic side chains are responsible for the characteristics of proteins in the near UV. We present the first systematic study of a large number of proteins focused on the accurate calculation of near-UV circular dichroism (CD) spectra. We report new parameter sets derived from ab initio calculations for benzene, phenol, and indole that describe the valence electronic transitions to the (1)L(b), (1)L(a), (1)B(b), and (1)B(a) states in the side chains of amino acids phenylalanine, tyrosine, and tryptophan. CD spectra were calculated, using the matrix method with the new side-chain parameters, for 30 proteins whose CD spectra and crystal structures have been made publicly available. The new parameter sets are fully self-consistent and yield near-UV spectra better than those obtained using previous parameter sets. The mean absolute errors for computed wild-type spectra in the near UV are reduced by a factor of approximately 2. A similiar reduction is found for the near-UV spectra (and difference spectra) of mutants involving aromatic amino acids. Empirical modifications to model vibronic coupling in the side-chain chromophore of phenylalanine offer no overall improvement. Protein CD calculations from first principles coupled with atomic-level modeling enhance the utility and interpretability of CD measurements in the near UV.     

Individual tyrosine side-chain contributions to circular dichroism of ribonuclease

Experimental and theoretical studies using site-directed mutants of ribonuclease A (RNase A) offer more extensive information on the tyrosine side-chain contributions to the circular dichroism (CD) of the enzyme. Bovine pancreatic RNase A has three exposed tyrosine residues (Tyr73, Tyr76, and Tyr115) and three buried tyrosine residues (Tyr25, Tyr92 and Tyr97). The difference CD spectra between the wild type and the mutants at pH 7.0 (Deltaepsilon(277,wt) - Deltaepsilon(277,mut)) show bands with more negative DeltaDeltaepsilon(277) values for Y73F and Y115F than those for Y25F and Y92F and bands with positive DeltaDeltaepsilon(277) values for Y76F and Y97F. The theoretical calculations are in good semiquantitative agreement for all the mutants. The pH difference spectrum (pH 11.3-7.0) for the wild type shows a negative band at 295 nm and an enhanced positive band at 245 nm. The three mutants at buried tyrosine sites and one mutant at an exposed tyrosine site (Y76F) exhibit pH-difference spectra that are similar to that of the wild type. In contrast, two mutants at exposed tyrosine sites (Y73F and Y115F) exhibit diminished 295-nm negative bands and, instead of positive bands at 245 nm, negative bands are observed. Our results indicate that Tyr73 and Tyr115, two of the exposed tyrosine residues, are the largest contributors to the 277- and 245-nm CD bands of RNaseA, but the buried tyrosine residues and the one remaining exposed residue also contribute to these bands. Disulfide contributions to the 277- and 240-nm bands and the peptide contribution to the 240-nm band are confirmed theoretically.     

The prediction of the circular dichroism of the benzene chromophore: TDDFT calculations and sector rules

The CD spectra of the series PhCH(Me)R, with R = Et (1), nPr (2), iPr (3), and tBu (4), are reported (1-3 for the first time) at room temperature in the 185-280 nm range and at 183 K. These purely hydrocarbon compounds represent the simplest chiral systems containing the phenyl chromophore and exhibit Cotton effects exclusively allied with the benzene transitions. The bands in 1La and 1Lb regions were checked against the available sector rules, with discordant outcomes. Time-dependent density-functional theory calculations, with various functionals and basis sets tested, correctly reproduced the prominent CD bands observed for 1-4.     

Characterization of Met95 mutants of a heme-regulated phosphodiesterase from Escherichia coli. Optical absorption, magnetic circular dichroism, circular dichroism, and redox potentials.

On the basis of amino acid sequences and crystal structures of similar enzymes, it is proposed that Met95 of the heme-regulated phosphodiesterase from Escherichia coli (Ec DOS) acts as a heme axial ligand. In accordance with this proposal, the Soret and visible optical absorption and magnetic circular dichroism spectra of the Fe(II) complexes of the Met95Ala and Met95Leu mutant proteins indicate that these complexes are five-coordinated high-spin, suggesting that Met95 is an axial ligand for the Fe(II) complex. However, the Fe(III) complexes of these mutants are six-coordinated low-spin, like the wild-type enzyme. The latter spectral findings are inconsistent with the proposal that the axial ligand to the Fe(III) heme is Met95. To determine the possibility of a redox-dependent ligand switch in Ec DOS, we further analyzed Soret CD spectra and redox potentials, which provide direct evidence on the environmental structure of the heme protein. CD spectra of Fe(III) Met95 mutants were all different from those of the wild-type protein, suggesting indirect coordination of Met95 to the Fe(III) wild-type heme. The redox potentials of the Met95Leu, Met95Ala and Met95His mutants were considerably lower than that of the wild-type enzyme (+70 mV) at -1, -26, and -122 mV vs. SHE, respectively. Thus, it is reasonable to speculate that water (or hydroxy anion) interacting with Met95, rather than Met95 itself, is the axial ligand to the Fe(III) heme.     

The calculated circular dichroism of polyproline II in the polarizability approximation

The circular dichroism (CD) spectrum of polyproline II (PPII) has heretofore been moderately well calculated from exciton theory only at the expense of assuming unreasonable chain conformations and accepting a conservative spectrum in the 180–250-nm region (which is not observed). We have incorporated far uv transitions in the polarizability approximation and, together with the π₂π* transition, have calculated the resulting correction to the exciton model. This has been accompanied by a modified assignment of the ππ* transition in PPII, and a simultaneous calculation of the absorption and CD spectra of the α-helix, β structure, PPI, and PPII. We obtain good agreement with the observed CD spectrum of PPII in the 180–250-nm region for acceptable chain conformations. In addition, we predict a negative CD into the far uv, in agreement with recent experimental observations. Our calculations also reproduce features of the far uv CD spectrum of the α-helix, and are in agreement with the CD spectra of the β chain and PPI. The calculated CD of the unordered polypeptide chain is not significantly influenced by far uv contributions, indicating that our previous calculation is valid for such a system. These results demonstrate the importance of incorporating far uv transitions in order to achieve an adequate theoretical explanation of the CD spectra of polypeptides.     

Effects of five-tryptophan mutations on structure, stability and function of Escherichia coli dihydrofolate reductase

To elucidate the roles of tryptophan residues in the structure, stability, and function of Escherichia coli dihydrofolate reductase (DHFR), its five tryptophan residues were replaced by site-directed mutagenesis with leucine, phenylalanine or valine (W22F, W22L, W30L, W47L, W74F, W74L, W133F, and W133V). Far-ultraviolet circular dichroism (CD) spectra of these mutants reveal that exciton coupling between Trp47 and Trp74 strongly affects the peptide CD of wild-type DHFR, and that Trp133 also contributes appreciably. No additivity was observed in the contributions of individual tryptophan residues to the fluorescence spectrum of wild-type DHFR, Trp74 having a dominant effect. These single-tryptophan mutations induce large changes in the free energy of urea unfolding, which showed values of 1.79-7.14 kcal/mol, compared with the value for wild-type DHFR of 6.08 kcal/mol. Analysis of CD and fluorescence spectra suggests that thermal unfolding involves an intermediate with the native-like secondary structure, the disrupted Trp47-Trp74 exciton coupling, and the solvent-exposed Trp30 and Trp47 side chains. All the mutants except W22L (13%) retain more than 50% of the enzyme activity of wild-type DHFR. These results demonstrate that the five tryptophan residues of DHFR play important roles in its structure and stability but do not crucially affect its enzymatic function.     

Crevice-forming mutants of bovine pancreatic trypsin inhibitor: stability changes and new hydrophobic surface.

Four mutants of bovine pancreatic trypsin inhibitor (BPTI) with replacements in the rigid core result in the creation of deep crevices on the surface of the protein. Other than crevices at the site of the mutation, few other differences are observed in the crystal structures of wild-type BPTI and the mutants F22A, Y23A, N43G, and F45A. These mutants are highly destabilized relative to wild type (WT). The differences between WT and mutants in the free energy change associated with cooperative folding/unfolding, delta delta G0 (WT-->mut), have been measured by calorimetry, and they are in good agreement with delta delta G0(WT-->mut) values from hydrogen exchange rates. For F22A the change in free energy difference is about 1.7 kcal/mol at 25 degrees C; for the other three mutants it is in the range of 5-7 kcal/mol at 25 degrees C. The experimental delta delta G0(WT-->mut) values of F22A, Y23A, and F45A are reasonably well accounted for as the sum of two terms: the difference in transfer free energy change, and a contribution from exposure to solvent of new surface (Eriksson, A.E., et al., 1992, Science 255, 178-183), if the recently corrected transfer free energies and surface hydrophobicities (De Young, L. & Dill, K., 1990, J. Phys. Chem. 94, 801-809; Sharp, K.A., et al., 1991a, Science 252, 106-109) are used and only nonpolar surface is taken into account. In N43G, three protein-protein hydrogen bonds are replaced by protein-water hydrogen bonds.(ABSTRACT TRUNCATED AT 250 WORDS)     

Tryptophan mutants of intestinal fatty acid-binding protein: ultraviolet absorption and circular dichroism studies.

An UV absorption and CD study of intestinal fatty acid-binding protein is presented. Since there are only two Trp residues in the molecule, two single-Trp mutants were prepared to deconvolute their signals. The individual contribution of the eight Phe and four Tyr residues was not established; however, Phe global contribution is relatively free of interferences from the other chromophores and was observed directly. CD spectra showed that Phe vibronic structure was unusually sharp and seems to monitor very specific details in the three-dimensional structure. The global signal from Tyr was assigned only approximately due to band broadening and overlapping. At the upper end of the CD spectrum, strong positive (1)L(b) Trp transitions from Trp 82 and strong negative (1)L(b) Trp transitions from Trp 6 were observed. (1)L(a) transitions were overall weak, positive for Trp 82 and negative for Trp 6, nearly cancelling each other out in the final spectrum. The above assignment is of practical and fundamental interest to monitor folding, binding, and molecular dynamics down to microdomain resolution. The assignment of Trp bands allowed comparison with previous data from CRABP1, another member of the IFABP family with 28% identical residues. It was found that structural homology extends beyond sequence and tertiary fold to include optical properties of equivalent Trp residues in the structure.     

Independent tyrosyl contributions to the CD of Ff gene 5 protein and the distinctive effects of Y41H and Y41F mutants on protein-protein cooperative interactions

The gene 5 protein (g5p) of the Ff virus contains five Tyr, individual mutants of which have now all been characterized by CD spectroscopy. The protein has a dominant tyrosyl 229-nm L(a) CD band that is shown to be approximately the sum of the five individual Tyr contributions. Tyr41 is particularly important in contributing to the high cooperativity with which the g5p binds to ssDNA, and Y41F and Y41H mutants are known to differ in dimer-dimer packing interactions in crystal structures. We compared the solution structures and binding properties of the Y41F and Y41H mutants using CD spectroscopy. Secondary structures of the mutants were similar by CD analyses and close to those derived from the crystal structures. However, there were significant differences in the binding properties of the two mutant proteins. The Y41H protein had an especially low binding affinity and perturbed the spectrum of poly[d(A)] in 2 mM Na(+) much less than did Y41F and the wild-type gene 5 proteins. Moreover, a change in the Tyr 229 nm band, assigned to the perturbation of Tyr34 at the dimer-dimer interface, was absent in titrations with the Y41H mutant under low salt conditions. In contrast, titrations with the Y41H mutant in 50 mM Na(+) exhibited typical CD changes of both the nucleic acid and the Tyr 229-nm band. Thus, protein-protein and g5p-ssDNA interactions appeared to be mutually influenced by ionic strength, indicative of correlated changes in the ssDNA binding and cooperativity loops of the protein or of indirect structural constraints.     

Circular dichroism of diglycosyl dichalcogenides in solution and solid state

Solution and solid-state CD spectra of nine peracetylated and deacetalyted diglycosyl disulfides were measured to study the relationship between the low-energy CD transitions (n1-->sigma*(S--S) and n2-->sigma*(S--S)) and helicity of the inherently chiral disulfide chromophore as perturbed by chiral carbohydrate moieties. The solid-state CD spectra were directly correlated with the reported X-ray structures of Ac4GlcSSGlcAc4 and Ac4GlcSSGalAc4, and the CD data revealed that all the disulfides have M helicity with C1--S--S--C1' dihedral angles -90 degrees < phi < 0 degrees both in solution and in the solid state. A TDDFT CD study was carried out on dimethyl diselenide which confirmed that the same quadrant rule is relevant between the signs of the low-energy CD transitions and the diselenide torsional angle as formulated previously for the disulfide chromophore. The CD spectra of Ac4GlcSeSeGlcAc4 measured in solution and in the solid state were correlated with its X-ray structure and reproduced well by TDDFT CD calculations performed on its tetra-O-acetyl derivative.     

Resolving near-ultraviolet circular dichroism spectra of single trp mutants in tear lipocalin

Near-ultraviolet circular dichroism (near-UV CD) spectra of tryptophan residues in proteins are complicated because the line shapes are derived from the overlap of both the 1L(a) and the 1L(b) electronic bands that vary independently. Contributing to this complexity, tryptophan near-UV CD spectra differ in the relative amplitude of the 0-0 vibronic band compared to the rest of the 1L(b) spectrum, an inherent feature that may result in poor fitting. To resolve this problem, a computer program that incorporated the separation of the 0-0 transition of 1L(b) component from the rest of the 1L(b) was written in LabVIEW and its amplitude was allowed to vary independently. This method showed dramatically improved fitting of 1L(a) and 1L(b) components in the near-UV CD tryptophan spectra in tear lipocalin mutants featuring low intensity of the 0-0 1L(b) component. Side chain dynamic characteristics (mobility and accessibility to the solvent) identified from different spectroscopic techniques were related to differences in Trp near-UV CD spectra. This method is broadly applicable to different types of Trp near-UV CD spectra.     

Stereochemical study of tolperisone, a muscle relaxant agent, by circular dichroism and ultraviolet spectroscopy

The stereochemistry of tolperisone, a chiral aryl-alkyl basic ketone was investigated by means of circular dichroism (CD) and ultraviolet (UV) spectroscopy. The unusually high optical activity of tolperisone hydrochloride in the n-->pi* region is interpreted by the presence of a chiral conformer in solution. For stereochemical reasons, the C = O group and the aromatic moiety lack coplanarity by forming an inherently dissymetric chromophore, of M helicity. Similar helicity prevails in the crystal phase, according to the solid-state CD spectrum of (-)-tolperisone HCl salt. The chirality rule proposed by Snatzke for nonplanar benzoyl chromophores predicts the absolute configuration of (-)-tolperisone hydrochloride to be R, in agreement with other alpha-methyl-beta-amino-ketones.     

Extreme heat- and pressure-resistant 7-kDa protein P2 from the archaeon Sulfolobus solfataricus is dramatically destabilized by a single-point amino acid substitution

This study reports the characterization of the recombinant 7-kDa protein P2 from Sulfolobus solfataricus and the mutants F31A and F31Y with respect to temperature and pressure stability. As observed in the NMR, FTIR, and CD spectra, wild-type protein and mutants showed substantially similar structures under ambient conditions. However, midpoint transition temperatures of the denaturation process were 361, 334, and 347 K for wild type, F31A, and F31Y mutants, respectively: thus, alanine substitution of phenylalanine destabilized the protein by as much as 27 K. Midpoint transition pressures for wild type and F31Y mutant could not be accurately determined because they lay either beyond (wild type) or close to (F31Y) 14 kbar, a pressure at which water undergoes a phase transition. However, a midpoint transition pressure of 4 kbar could be determined for the F31A mutant, implying a shift in transition of at least 10 kbar. The pressure-induced denaturation was fully reversible; in contrast, thermal denaturation of wild type and mutants was only partially reversible. To our knowledge, both the pressure resistance of protein P2 and the dramatic pressure and temperature destabilization of the F31A mutant are unprecedented. These properties may be largely accounted for by the role of an aromatic cluster where Phe31 is found at the core, because interactions among aromatics are believed to be almost pressure insensitive; furthermore, the alanine substitution of phenylalanine should create a cavity with increased compressibility and flexibility, which also involves an impaired pressure and temperature resistance. Proteins 29:381–390, 1997. © 1997 Wiley-Liss, Inc.     

The circular dichroism of aromatic polypeptides: Theoretical studies of poly-L-phenylalanine and some para-substituted derivatives

We have calculated rotational strengths and circular dichroism (CD) curves for sidechain and backbone transitions in poly-L -Phenylalanine (PLP), POLY-p-amino-L -phenylalanine (PPALP), poly-p-chloro-L -phenylalanine (PPCLP), poly-o-acetyl-L -tyrosine (POALT), and poly-p-nitro-L -phenylanine (PPNLP), using methods applied previously to poly-L -tyrosine (PLT). Comparison of the theoretical CD curves with available experimental data for PLP and PPALP indicate that these polypeptides form right-handed helices with side-chain conformations similar to that of PLT. For PPNLP, where experimental data are also available, no conformational assignment could be made, as none of the calculated curves gave good agreement with experiment. Possible reasons for this lack of agreement are discussed. For the other two polypeptides, PPCLP and POALT, although no experimental data are yet available, the calculated curves indicate that an unambiguous assignment should be possible. For the conformations (RA and LA) in which the side chains are packed more loosely, there are strong similarities in the calculated CD curves of a particular conformation, regardless of the para substituent. In the tighter R1 and L1 conformations, few generalizations can be drawn, each derivative having a distinctive pattern. In PLP, PPCLP, and POALT, where the side-chain L^a band is in the 200–210 nm region, the L1 conformation exhibits a negative nπ* rotational strength, opposite to that expected for a left-handed helix. One must therefore be cautious about assigning the helix sense of aromatic polypeptides on the basis of the sign of the nπ* CD band. Side-chain nπ* transitions present in POALT and PPNLP were found to have small rotational strength.     

PPP calculations of the circular dichroism spectra of some dinucleoside phosphates

Circular dichriosm (CD) spectra have been calculated for serveral dinucleoside phosphates using a variant of the Pariser-Parr-Pople π-electron molecular orbital method. This method does not require the prior knowledge of the experimental absorption spectra of transition moments of the bases forming the dinucleoside phosphates. Calculated spectra were obtained in good agreement with experimental spectra for four dinucleoside phosphates, ApA, UpU, GpA, and UpA, and reasonable agreement was obtained for ApG and ApU. The effect of changing conformation on the CD spectrum was studied for ApA, UpU, UpA, and ApU; the spectra of UpU, UpA, and ApU were sensitive to small change in conformation, whereas ApA was insensitive over the range of conformation studied. Further studies await detailed knowledge of the structure of dinucleoside phosphates in solution.     

Understanding the changes in the circular dichroism of light harvesting complex II upon varying its pigment composition and organization

In this work we modeled the circular dichroism (CD) spectrum of LHCII, the main light harvesting antenna of photosystem II of higher plants. Excitonic calculations are performed for a monomeric subunit, taken from the crystal structure of trimeric LHCII from spinach [Liu, Z. F., Yan, H. C., Wang, K. B., Kuang, T. Y., Zhang, J. P., Gui, L. L., An, X. M., and Chang, W. R. (2004) Nature 428, 287-292]. All of the major features of the CD spectrum above 450 nm are satisfactorily reproduced, and possible orientations of the Chl and carotenoid transition dipole moments are identified. The obtained modeling parameters are used to simulate the CD spectra of two complexes with altered pigment composition: a mutant lacking Chls a 611-612 and a complex lacking the carotenoid neoxanthin. By removing the relevant pigment(s) from the structure, we are able to reproduce their spectra, which implies that the alteration does not disturb the overall structure. The CD spectrum of trimeric LHCII shows a reversed relative intensity of the two negative bands around 470 and 490 nm as compared to monomeric LHCII. The simulations reproduce this reversal, indicating that it is mainly due to interactions between chromophores in different monomeric subunits, and the trimerization does not induce observable changes in the monomeric structure. Our simulated spectrum resembles one of two different trimeric CD spectra reported in literature. We argue that the differences in the experimental trimeric CD spectra are caused by changes in the strength of the monomer-monomer interactions due to the differences in detergents used for the purification of the complexes.     

Solution conformations of cyclosporins and magnesium-cyclosporin complexes determined by vibrational circular dichroism

Vibrational circular dichroism (VCD) spectroscopy was used to investigate the solution conformations of cyclosporins A, C, D, G, and H in CDCl(3), in the amide I and NH/OH-stretching regions, and their corresponding magnesium complexes in CD(3)CN, in the amide I region. VCD spectra are sensitive to the chiral arrangement of Cdbond;O and NH bonds in this cyclic undecapeptide. Calculations of molecular geometries, as well as IR and VCD intensities of model cyclosporin fragments that include the intramolecular hydrogen bonds of the crystal conformations of cyclosporins A and H (CsA and CsH), were carried out at the density functional theory (DFT; BPW91 functional/6-31G* basis set) level. The good agreement between IR and VCD spectra from experiment and DFT calculations provides evidence that the crystal conformation of CsA is dominant in CDCl(3) solution; CsH, however, assumes both an intramolecularly hydrogen-bonded crystal conformation and more open forms in solution. Comparisons of the experimental and calculated VCD spectra in the NH/OH-stretching region of the noncomplexed cyclosporins indicate that conformers with both free and hydrogen-bonded NH and OH groups are present in solution. Differences between the IR and VCD spectra for the metal-free and magnesium-complexed cyclosporins are indicative of strong interactions between cyclosporins and magnesium ions.     

Infrared magnetic circular dichroism of myoglobin derivatives.

By use of a newly constructed CD instrument, infrared magnetic circular dichroism (MCD) spectra were observed for various myoglobin derivatives. The ferric high spin myoglobin derivatives such as fluoride, water and hydroxide complexes, commonly exhibited the MCD spectra consisting of positive A terms. Therefore, the results reinforced the assignment that the infrared band is the charge transfer transition to the degenerate excited state (eg (dpi)). Since the fraction of A term estimated was approximately 80% for myoglobin fluoride and approximately 35% for myoglobin water, the effective symmetry for myoglobin fluoride is determined to be as close as D4h, while that for myoglobin water seems to have lower symmetry components. The ferric low spin derivatives such as myoglobin cyanide, myoglobin imidazole and myoglobin azide showed positive MCD spectra which are very similar to the electronic absorption spectra. These MCD spectra were assigned to the charge transfer transitions from porphyrin pi to iron d orbitals on the ground that they were observed only for the ferric low spin groups and insensitive to the axial ligands. The lack of temperature dependence in the MCD magnitude indicated that the MCD spectra are attributable to the Faraday B terms. Deoxymyoglobin, the ferrous high spin derivative, had fairly strong positive MCD around 760 nm with an anisotropy factor (delta epsilon/epsilon) of 1.4-10(-4). It shows some small MCD bands from 800 to 1800 nm. Among the ferrous low spin derivatives, carbonmonoxymyoglobin did not give any observable MCD in the infrared region while oxymyoglobin seemed to have significant MCD in the range from 700 to 1000 nm.