Globular Protein Ultraviolet Circular Dichroic Spectra. Calculation from Crystal Structures via the Dipole Interaction Model

Introducing DInaMo: A Package for Calculating Protein Circular Dichroism Using Classical Electromagnetic Theory

The dipole interaction model is a classical electromagnetic theory for calculating circular dichroism (CD) resulting from the π-π* transitions of amides. The theoretical model, pioneered by J. Applequist, is assembled into a package, DInaMo, written in Fortran allowing for treatment of proteins. DInaMo reads Protein Data Bank formatted files of structures generated by molecular mechanics or reconstructed secondary structures. Crystal structures cannot be used directly with DInaMo; they either need to be rebuilt with idealized bond angles and lengths, or they need to be energy minimized to adjust bond lengths and bond angles because it is common for crystal structure geometries to have slightly short bond lengths, and DInaMo is sensitive to this. DInaMo reduces all the amide chromophores to points with anisotropic polarizability and all nonchromophoric aliphatic atoms including hydrogens to points with isotropic polarizability; all other atoms are ignored. By determining the interactions among the chromophoric and nonchromophoric parts of the molecule using empirically derived polarizabilities, the rotational and dipole strengths are determined leading to the calculation of CD. Furthermore, ignoring hydrogens bound to methyl groups is initially explored and proves to be a good approximation. Theoretical calculations on 24 proteins agree with experiment showing bands with similar morphology and maxima.

Advantages of electronic circular dichroism detection for the stereochemical analysis and characterization of drugs and natural products by liquid chromatography

The need for analytical methods for the determination of the enantiomeric excess of chiral compounds increased significantly in the last decades, and enantioselective separation techniques resulted particularly efficient to this purpose. Moreover, when detection systems based on chiroptical properties (optical rotation or circular dichroism) are employed in high-performance liquid chromatography (HPLC), the stereochemistry of a chiral analyte can be fully determined. Indeed, the coupling of HPLC with chiroptical detection systems allows the simultaneous assessment of the absolute configuration of stereoisomers and the evaluation of the enantiomeric/diastereomeric excess of samples. These features are particularly important in the study of drugs and natural products provided with biological activity, because the assignment of their absolute stereochemistry is essential to establish reliable structure-activity relationships. The following review aims to discuss the analytical advantages arising from the employment of electronic circular dichroism (ECD) detection systems in stereochemical analysis by HPLC upon chiral and non-chiral stationary phases and their use for the stereochemical characterization of chiral drugs and natural compounds. The different methods for the correlation between absolute stereochemistry and chiroptical properties are critically discussed. Relevant HPLC applications of ECD detection systems are then reported, and their analytical advantages are highlighted. For instance, the importance of the concentration-independent anisotropy factor (g-factor; g=Δɛ/ɛ) for the determination of the stereoisomeric composition of samples upon non-chiral stationary phases is underlined, since its sensitivity makes ECD detection very well suited for the enantioselective analysis of large libraries of chiral compounds in relatively short times.

Effects of selenium on the structure and function of recombinant human S-adenosyl-L-methionine dependent arsenic (+3 oxidation state) methyltransferase in E. coli

The effects of Se(IV) on the structure and function of recombinant human arsenic (+3 oxidation state) methyltransferase (AS3MT) purified from the cytoplasm of Escherichia coli were studied. The coding region of human AS3MT complementary DNA was amplified from total RNA extracted from HepG2 cell by reverse transcription PCR. Soluble and active human AS3MT was expressed in the E. coli with a Trx fusion tag under a lower induction temperature of 25 degrees C. Spectra (UV-vis, circular dichroism, and fluorescence) were first used to probe the interaction of Se(IV) and recombinant human AS3MT and the structure-function relationship of the enzyme. The recombinant human AS3MT had a secondary structure of 29.0% alpha-helix, 23.9% beta-pleated sheet, 17.9% beta-turn, and 29.2% random coil. When Se(IV) was added, the content of the alpha-helix did not change, but that of the beta-pleated sheet increased remarkably in the conformation of recombinant human AS3MT. Se(IV) inhibited the enzymatic methylation of inorganic As(III) in a concentration-dependent manner. The IC(50) value for Se(IV) was 2.38 muM. Double-reciprocal (1/V vs. 1/[inorganic As(III)]) plots showed Se(IV) to be a noncompetitive inhibitor of the methylation of inorganic As(III) by recombinant human AS3MT with a K (i) value of 2.61 muM. We hypothesized that Se(IV) interacts with the sulfhydryl group of cysteine(s) in the structural residues rather than the cysteines of the active site (Cys156 and Cys206). When Se(IV) was combined with cysteine(s) in the structural residues, the conformation of recombinant human AS3MT changed and the enzymatic activity decreased. Considering the quenching of tryptophan fluorescence, Cys72 and/or Cys226 are deduced to be primary targets for Se(IV).

Two Dimensional Electronic Correlation Spectroscopy of the npi* and pipi* Protein Backbone Transitions: A Simulation Study

The two dimensional (2D) photon echo spectrum of the amide ultraviolet (UV) bands of proteins are simulated. Two effective exciton Hamiltonian parameter sets developed by Woody and Hirst, which predict similar CD spectra, may be distinguished by their very different 2DUV spectra. These differences are enhanced in specific configurations of pulse polarizations which provide chirality-induced signals.

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.

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.

Theoretical UV Circular Dichroism of Cyclo(l-Proline-l-Proline)

MP2, DFT, and molecular mechanics (AMBER, CVFF, and CFF91) geometry optimizations were performed on the cyclic dipeptide cyclo(L-Pro-L-Pro) starting from crystal structure data. Three stable conformations were identified as energy minima by all methods, but assignment of relative energy varied between the methods. The pi-pi transition feature of the UV circular dichroic (CD) spectrum was predicted for each minimized structure using the classical physics method of the dipole interaction model. The model was sensitive to the different conformations. The UV-CD predictions were compared individually and as a Boltzmann-weighted composite with published experimental CD spectra [Bowman, R. L.; Kellerman, M.; Johnson, W. C., Jr. Biopolymers 1983, 22, 1045]. For all structures, the original parameters of Applequist [Applequist, J. J. Chem. Phys. 1979, 71, 4324] with a bandwidth of 3000 cm(-1) most accurately replicated experiment, except for the CFF91 structures, which matched experiment best with a bandwidth of 4000 cm(-1). The inclusion of solvent by a continuum model did not significantly alter the minimized geometries obtained by molecular or quantum mechanics, but it did have an effect on the relative predicted energies of CFF91 and B3LYP conformations. The overall effect of solvent inclusion was negligible when Boltzmann-weighted spectra were considered. Gas-phase CFF91 structures were also reasonably good for prediction of CD spectra, and when water was included via a continuum model for energy calculations, the weighting scheme resembled that of the higher-level weightings. The CD calculated using the MP2/6-311G structures and energies for weighting were most descriptive of the 180 nm negative band in the experimental CD, but red-shifted the location of the 205 nm band. DFT structures were comparably, though not identically, as descriptive of the first pi-pi band, and did a better job with placement of the second (positive) pi-pi band. DFT calculations were less sensitive to basis set effect than the MP2 calculations, with 6-31G results in close agreement with 6-311G. The results suggest that it is possible to use geometries obtained from a variety of different methods (molecular mechanical or quantum mechanical) with the classical physics dipole interaction model to qualitatively reproduce the UV CD of model amides.

Theoretical UV Circular Dichroism of Aliphatic Cyclic Dipeptides

Four cyclic dipeptides (piperazine-2,5-diones), cyclo(L-Pro-Gly), cyclo(L-Pro-L-Leu), cyclo(L-Ala-L-Ala), and cyclo(L-Pro-L-Ala), were modeled from crystal structure data. Conformations resulting from energy minimization using molecular mechanics were compared with traditional ab initio and density functional theory geometric optimizations for each dipeptide. In all computational cases, the gas phase was assumed. The pi-pi transition feature of the UV circular dichroic (CD) spectra was predicted for each peptide structure via the classical dipole interaction model. The dipole interaction model predicted CD spectra that qualitatively agreed with experiment when MP2 or DFT geometries were used. By coupling MP2 or DFT geometric optimizations with the classical physics method of the dipole interaction model, significantly better CD spectra were calculated than those using geometries obtained by molecular mechanics. Thus, one can couple quantum mechanical geometries with a classical physics model for calculation of circular dichroism.

Heme-heme interactions in tetramers and dimers of hemoglobin subunits: DeVoe theory calculations

Detectable exciton couplets arising from heme-heme interactions in the hemoglobin (Hb) tetramers of HbO(2) and deoxyHb were predicted by DeVoe theory. This prediction was supported by the observation of an exciton couplet in the CD difference spectrum between the Hb tetramer and the alphabeta dimer of HbCO. In this paper, DeVoe theory is used to calculate the heme-heme interactions in the CO complex of the Hb tetramer (alpha(2)beta(2)) and dimer (alphabeta), the systems studied by Goldbeck et al. The couplet strength of the resulting theoretical CD difference spectrum agrees well with experiment, thus confirming that heme-heme interactions contribute significantly to the CD of HbCO. Given that the heme-heme distances in HbCO are 25 A and more, it is highly likely that heme-heme interactions also contribute significantly to the CD of other multi-heme proteins, e.g., cytochrome c(3), cytochrome oxidase, cytochrome bc(1), etc., where the hemes are in closer proximity.

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.

Calculations of protein circular dichroism from first principles

Understanding the relationship between the amino acid sequence of a protein and its unique, compact 3D structure is one of the grand challenges in molecular biophysics. One particularly exciting approach is time-resolved electronic circular dichroism (CD) spectroscopy, which offers resolution on a nanosecond (or faster) time scale, although it does not provide the spatial resolution of techniques like X-ray crystallography or NMR. The thrust of our work is to underpin fast time scale spectroscopic studies of protein folding with a stronger theoretical foundation. Ultimately, we seek to use molecular dynamics simulations to study the influence of conformational dynamics and conformational transitions on the electronic CD spectra of proteins. We discuss how improved quantum chemical models of individual chromophores, including aromatic sidechains, can be incorporated into calculations of the electronic structure of proteins and their CD.

Dipole interaction model predicted pi-pi* circular dichroism of cyclo(L-Pro)3 using structures created by semi-empirical, ab initio, and molecular mechanics methods

Cyclo(l-Pro)3 (CP3) is a synthetic peptide created to model cis and torsionally strained peptide bonds that also exhibits a strong distinctive UV circular dichroic (CD) spectrum. Circular dichroic spectra were computed for the amide pi-pi* transition using the dipole interaction model for various conformations of the peptide. Conformations of CP3 were created initially from crystal data, and followed by energy minimizations via molecular mechanics using the cvff force field; the effects of additional geometric optimizations by semi-empirical and ab initio quantum mechanics were investigated. The CD spectra for each conformation were calculated using a variety of different parameters, and each result was compared with the published experimental spectrum [Deber, C.M., Scatturin, A., Vaidya, V.M. & Blout, E.R. (1970) Small cyclic proline peptides: UV absorption and CD. In: Peptides: Chemistry and Biochemistry, Proceedings of the First American Peptide Symposium (Weinstein, B., ed.), Marcel Dekker, New York pp. 163-173]. Herein, two distinct conformations, a C3 symmetric and an asymmetric form, gave CD predictions that separately did not resemble the experimental spectrum. Energy differences were predicted at various theoretical levels, including MP2 and density functional theory. When the predicted CD spectra for each conformation were multiplied by Boltzmann weighting factors created using heats of formation determined by the AM1 optimizations, the weighted composite CD spectrum created did resemble experiment for the pi-pi* transition indicating that both conformations may exist simultaneously in solution.

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.

An improved model for calculating the optical rotation of simple saccharides

A calculational model for the optical rotation (OR) at the sodium D-line of simple saccharides has been developed that eliminates deficiencies of a previous model. Conformational conclusions based on the earlier model are not affected, as established by a recalculation of the OR phi,psi-map of methyl 3-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside. The model relocates the strong long-wavelength sigma-sigma* circular dichroism (CD) component, which is mainly responsible for the NaD OR from 160 to below 130 nm, where it is now known to occur. That correction allows future modeling of CD bands of different origins that appear in the 150-190 nm region. In order to demonstrate the utility of the revised model, it was applied to calculating the OR of methyl 3-O-(alpha-L-rhamnopyranosyl)-alpha-L-rhamnopyranoside. The results provide experimental support for conformational conclusions derived from a molecular mechanics study of that molecule.