A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.     

Angular variances for internal bond rotations of side chains in GXG-based tripeptides derived from (13)C-NMR relaxation measurements: Implications to protein folding

The study of backbone and side-chain internal motions in proteins and peptides is crucial to having a better understanding of protein/peptide "structure" and to characterizing unfolded and partially folded states of proteins and peptides. To achieve this, however, requires establishing a baseline for internal motions and motional restrictions for all residues in the fully, solvent-exposed "unfolded state." GXG-based tripeptides are the simpliest peptides where residue X is fully solvent exposed in the context of an actual peptide. In this study, a series of GXG-based tripeptides has been synthesized with X being varied to include all twenty common amino acid residues. Proton-coupled and -decoupled (13)C-nmr relaxation measurements have been performed on these twenty tripeptides and various motional models (Lipari-Szabo model free approach, rotational anisotropic diffusion, rotational fluctuations within a potential well, rotational jump model) have been used to analyze relaxation data for derivation of angular variances and motional correlation times for backbone and side-chain chi(1) and chi(2) bonds and methyl group rotations. At 298 K, backbone motional correlation times range from about 50 to 85 ps, whereas side-chain motional correlation times show a much broader spread from about 18 to 80 ps. Angular variances for backbone phi,psi bond rotations range from 11 degrees to 23 degrees and those for side chains vary from 5 degrees to 24 degrees for chi(1) bond rotations and from 5 degrees to 27 degrees for chi(2) bond rotations. Even in these peptide models of the "unfolded state," side-chain angular variances can be as restricted as those for backbone and beta-branched (valine, threonine, and isoleucine) and aromatic side chains display the most restricted motions probably due to steric hinderence with backbone atoms. Comparison with motional data on residues in partially folded, beta-sheet-forming peptides indicates that side-chain motions of at least hydrophobic residues are less restricted in the partially folded state, suggesting that an increase in side-chain conformational entropy may help drive early-stage protein folding. Copyright 1999 John Wiley & Sons, Inc.     

Backrub-like backbone simulation recapitulates natural protein conformational variability and improves mutant side-chain prediction

Incorporation of effective backbone sampling into protein simulation and design is an important step in increasing the accuracy of computational protein modeling. Recent analysis of high-resolution crystal structures has suggested a new model, termed backrub, to describe localized, hinge-like alternative backbone and side-chain conformations observed in the crystal lattice. The model involves internal backbone rotations about axes between C-alpha atoms. Based on this observation, we have implemented a backrub-inspired sampling method in the Rosetta structure prediction and design program. We evaluate this model of backbone flexibility using three different tests. First, we show that Rosetta backrub simulations recapitulate the correlation between backbone and side-chain conformations in the high-resolution crystal structures upon which the model was based. As a second test of backrub sampling, we show that backbone flexibility improves the accuracy of predicting point-mutant side-chain conformations over fixed backbone rotameric sampling alone. Finally, we show that backrub sampling of triosephosphate isomerase loop 6 can capture the millisecond/microsecond oscillation between the open and closed states observed in solution. Our results suggest that backrub sampling captures a sizable fraction of localized conformational changes that occur in natural proteins. Application of this simple model of backbone motions may significantly improve both protein design and atomistic simulations of localized protein flexibility.     

Side-chain flexibility in protein-ligand binding: the minimal rotation hypothesis

The goal of this work is to learn from nature about the magnitudes of side-chain motions that occur when proteins bind small organic molecules, and model these motions to improve the prediction of protein-ligand complexes. Following analysis of protein side-chain motions upon ligand binding in 63 complexes, we tested the ability of the docking tool SLIDE to model these motions without being restricted to rotameric transitions or deciding which side chains should be considered as flexible. The model tested is that side-chain conformational changes involving more atoms or larger rotations are likely to be more costly and less prevalent than small motions due to energy barriers between rotamers and the potential of large motions to cause new steric clashes. Accordingly, SLIDE adjusts the protein and ligand side groups as little as necessary to achieve steric complementarity. We tested the hypothesis that small motions are sufficient to achieve good dockings using 63 ligands and the apo structures of 20 different proteins and compared SLIDE side-chain rotations to those experimentally observed. None of these proteins undergoes major main-chain conformational change upon ligand binding, ensuring that side-chain flexibility modeling is not required to compensate for main-chain motions. Although more frugal in the number of side-chain rotations performed, this model substantially mimics the experimentally observed motions. Most side chains do not shift to a new rotamer, and small motions are both necessary and sufficient to predict the correct binding orientation and most protein-ligand interactions for the 20 proteins analyzed.     

Motional dynamics of residues in a beta-hairpin peptide measured by 13C-NMR relaxation.

Structurally characterizing partially folded peptides is problematic given the nature of their transient conformational states. 13C-NMR relaxation data can provide information on the geometry of bond rotations, motional restrictions, and correlated bond rotations of the backbone and side chains and, therefore, is one approach that is useful to assess the presence of folded structure within a conformational ensemble. A peptide 12mer, R1GITVNG7KTYGR12, has been shown to partially fold in a relatively stable beta-hairpin conformation centered at NG. Here, five residues, G2, V5, G7, Y10, G11, were selectively 13C-enriched, and 13C-NMR relaxation experiments were performed to obtain auto- and cross-correlation motional order parameters, correlation times, bond rotation angular variances, and bond rotational correlation coefficients. Our results indicate that, of the three glycines, G7 within the hairpin beta-turn displays the most correlated phi(t),psi(t) rotations with its axis of rotation bisecting the angle defined by the H-C-H bonds. These positively correlated bond rotations give rise to "twisting" type motions of the HCH group. V5 and Y10 phi,psi bond rotations are also positively correlated, with their CbetaCalphaH groups undergoing similar "twisting" type motions. Motions of near-terminal residues G2 and G11 are less restricted and less correlated and are best described as wobbling-in-a-cone. V5 and Y10 side-chain motions, aside from being highly restricted, were found to be correlated with phi,psi bond rotations. At 303 K, where the hairpin is considered "unfolded," the peptide exists in a transient, collapsed state because backbone and side-chain motions of V5, G7, and Y10 remain relatively restricted, unlike their counterparts in GXG-based tripeptides. These results provide unique information toward understanding conformational variability in the unfolded state of proteins, which is necessary to solve the protein folding problem.     

Side-chain control of beta-peptide secondary structures.

As one of the most important families of non-natural polymers with the propensity to form well-defined secondary structures, the beta-peptides are attracting increasing attention. The compounds incorporating beta-amino acid residues have found various applications in medicinal chemistry and biochemistry. The conformational pool of beta-peptides comprises several periodic folded conformations, which can be classified as helices, and nonpolar and polar strands. The latter two are prone to form pleated sheets. The numerous studies that have been performed on the side-chain dependence of the stability of the folded structures allow certain conclusions concerning the principles of design of the beta-peptide foldamers. The folding propensity is influenced by local torsional, side-chain to backbone and long-range side-chain interactions. Although beta-peptide foldamers are sensitive to solvent, the systematic choice of the side-chain pattern and spatiality allows the design of the desired specific secondary structure. The application of beta-peptide foldamers may open up new directions in the synthesis of highly organized artificial tertiary structures with biochemical functions.     

A survey of amino acid side-chain interactions in 21 proteins

Based on the atomic co-ordinate data for 21 representative proteins, the frequencies of long-range interactions between side-chain groups and 15 different types of side-chain atoms have been determined. The observed frequencies are compared to the results expected for random association in order to define a scale of relative affinities. Thirty-five residue-atom pairs exhibit frequencies of interaction that differ by at least 50% from the expected values. The amino acids tend to fall into three classes: non-polar, neutral and polar amino acids. The data are regrouped in a different way to determine the average affinity of each amino acid side-chain group for all other types of side-chain groups. Fourteen side-chain pairs have at least 50% fewer interactions than expected, while 21 side-chain pairs have at least 50% more interactions than expected. Unusual patterns of association are discussed and compared with current ideas about the organization of protein structure.     

Accurate and efficient description of protein vibrational dynamics: comparing molecular dynamics and Gaussian models

Current all-atom potential based molecular dynamics (MD) allows the identification of a protein's functional motions on a wide-range of timescales, up to few tens of nanoseconds. However, functional, large-scale motions of proteins may occur on a timescale currently not accessible by all-atom potential based MD. To avoid the massive computational effort required by this approach, several simplified schemes have been introduced. One of the most satisfactory is the Gaussian network approach based on the energy expansion in terms of the deviation of the protein backbone from its native configuration. Here, we consider an extension of this model that captures in a more realistic way the distribution of native interactions due to the introduction of effective side-chain centroids. Since their location is entirely determined by the protein backbone, the model is amenable to the same exact and computationally efficient treatment as previous simpler models. The ability of the model to describe the correlated motion of protein residues in thermodynamic equilibrium is established through a series of successful comparisons with an extensive (14 ns) MD simulation based on the AMBER potential of HIV-1 protease in complex with a peptide substrate. Thus, the model presented here emerges as a powerful tool to provide preliminary, fast yet accurate characterizations of protein near-native motion.     

Identification of domains and domain interface residues in multidomain proteins from graph spectral method

We present a novel method for the identification of structural domains and domain interface residues in proteins by graph spectral method. This method converts the three-dimensional structure of the protein into a graph by using atomic coordinates from the PDB file. Domain definitions are obtained by constructing either a protein backbone graph or a protein side-chain graph. The graph is constructed based on the interactions between amino acid residues in the three-dimensional structure of the proteins. The spectral parameters of such a graph contain information regarding the domains and subdomains in the protein structure. This is based on the fact that the interactions among amino acids are higher within a domain than across domains. This is evident in the spectra of the protein backbone and the side-chain graphs, thus differentiating the structural domains from one another. Further, residues that occur at the interface of two domains can also be easily identified from the spectra. This method is simple, elegant, and robust. Moreover, a single numeric computation yields both the domain definitions and the interface residues.     

Structure and internal mobility of proteins: a molecular dynamics study of hen egg white lysozyme

The structure and internal motions of the protein hen egg white lysozyme are studied by analysis of simulation and experimental data. A molecular dynamics simulation and an energy minimization of the protein in vacuum have been made and the results compared with high-resolution structures and temperature factors of hen egg white lysozyme in two different crystal forms and of the homologous protein human lysozyme. The structures obtained from molecular dynamics and energy minimization have root-mean-square deviations for backbone atoms of 2.3 Å and 1.1–1.3 Å, respectively, relative to the crystal structures; the different crystal structures have root-mean-square deviations of 0.73–0.81 Å for the backbone atoms. In comparing the backbone dihedral angles, the difference between the dynamics and the crystal structure on which it is based is the same as that between any two crystal structures. The internal fluctuations of atomic positions calculated from the molecular dynamics trajectory agree well with the temperature factors from the three structures. Simulation and crystal results both show that there are large motions for residues involved in exposed turns of the backbone chain, relatively smaller motions for residues involved in the middle of helices or β-sheet structures, and relatively small motions of residues near disulfide bridges. Also, both the simulation and crystal data show that side-chain atoms have larger fluctuations than main-chain atoms. Moreover, the regions that have large deviations among the x-ray crystal structures, which indicates flexibility, are found to have large fluctuations in the simulation.     

Application of cross-correlated NMR spin relaxation to the zinc-finger protein CRP2(LIM2): evidence for collective motions in LIM domains

The solution structure of quail CRP2(LIM2) was significantly improved by using an increased number of NOE constraints obtained from a 13C,15N-labeled protein sample and by applying a recently developed triple-resonance cross-correlated relaxation experiment for the determination of the backbone dihedral angle psi. Additionally, the relative orientation of the 15N(i)-1HN(i) dipole and the 13CO(i) CSA tensor, which is related to both backbone angles phi and psi, was probed by nitrogen-carbonyl multiple-quantum relaxation and used as an additional constraint for the refinement of the local geometry of the metal-coordination sites in CRP2(LIM2). The backbone dynamics of residues located in the folded part of CRP2(LIM2) have been characterized by proton-detected 13C'(i-1)-15N(i) and 15N(i)-1HN(i) multiple-quantum relaxation, respectively. We show that regions having cross-correlated time modulation of backbone isotropic chemical shifts on the millisecond to microsecond time scale correlate with residues that are structurally altered in the mutant protein CRP2(LIM2)R122A (disruption of the CCHC zinc-finger stabilizing side-chain hydrogen bond) and that these residues are part of an extended hydrogen-bonding network connecting the two zinc-binding sites. This indicates the presence of long-range collective motions in the two zinc-binding subdomains. The conformational plasticity of the LIM domain may be of functional relevance for this important protein recognition motif.     

Structure-based analysis of protein-RNA interactions using the program ENTANGLE

Until recently, drawing general conclusions about RNA recognition by proteins has been hindered by the paucity of high-resolution structures. We have analyzed 45 PDB entries of protein-RNA complexes to explore the underlying chemical principles governing both specific and non-sequence specific binding. To facilitate the analysis, we have constructed a database of interactions using ENTANGLE, a JAVA-based program that uses available structural models in their PDB format and searches for appropriate hydrogen bonding, stacking, electrostatic, hydrophobic and van der Waals interactions. The resulting database of interactions reveals correlations that suggest the basis for the discrimination of RNA from DNA and for base-specific recognition. The data illustrate both major and minor interaction strategies employed by families of proteins such as tRNA synthetases, ribosomal proteins, or RNA recognition motifs with their RNA targets. Perhaps most surprisingly, specific RNA recognition appears to be mediated largely by interactions of amide and carbonyl groups in the protein backbone with the edge of the RNA base. In cases where a base accepts a proton, the dominant amino acid donor is arginine, whereas in cases where the base donates a proton, the predominant acceptor is the backbone carbonyl group, not a side-chain group. This is in marked contrast to DNA-protein interactions, which are governed predominantly by amino acid side-chain interactions with functional groups that are presented in the accessible major groove. RNA recognition often proceeds through loops, bulges, kinks and other irregular structures that permit use of all the RNA functional groups and this is seen throughout the protein-RNA interaction database.     

Measurements of protein sequence-structure correlations

Correlations between protein structures and amino acid sequences are widely used for protein structure prediction. For example, secondary structure predictors generally use correlations between a secondary structure sequence and corresponding primary structure sequence, whereas threading algorithms and similar tertiary structure predictors typically incorporate interresidue contact potentials. To investigate the relative importance of these sequence-structure interactions, we measured the mutual information among the primary structure, secondary structure and side-chain surface exposure, both for adjacent residues along the amino acid sequence and for tertiary structure contacts between residues distantly separated along the backbone. We found that local interactions along the amino acid chain are far more important than non-local contacts and that correlations between proximate amino acids are essentially uninformative. This suggests that knowledge-based contact potentials may be less important for structure predication than is generally believed.     

Four-dimensional heteronuclear correlation experiments for chemical shift assignment of solid proteins

Chemical shift assignment is the first step in all established protocols for structure determination of uniformly labeled proteins by NMR. The explosive growth in recent years of magic-angle spinning (MAS) solid-state NMR (SSNMR) applications is largely attributable to improved methods for backbone and side-chain chemical shift correlation spectroscopy. However, the techniques developed so far have been applied primarily to proteins in the size range of 5–10 kDa, despite the fact that SSNMR has no inherent molecular weight limits. Rather, the degeneracy inherent to many 2D and 3D SSNMR spectra of larger proteins has prevented complete unambiguous chemical shift assignment. Here we demonstrate the implementation of 4D backbone chemical shift correlation experiments for assignment of solid proteins. The experiments greatly reduce spectral degeneracy at a modest cost in sensitivity, which is accurately described by theory. We consider several possible implementations and investigate the CANCOCX pulse sequence in detail. This experiment involves three cross polarization steps, from H to CA[i], CA[i] to N[i], and N[i] to C′[i−1], followed by a final homonuclear mixing period. With short homonuclear mixing times (<20 ms), backbone correlations are observed with high sensitivity; with longer mixing times (>200 ms), long-range correlations are revealed. For example, a single 4D experiment with 225 ms homonuclear mixing time reveals ∼200 uniquely resolved medium and long-range correlations in the 56-residue protein GB1. In addition to experimental demonstrations in the 56-residue protein GB1, we present a theoretical analysis of anticipated improvements in resolution for much larger proteins and compare these results in detail with the experiments, finding good agreement between experiment and theory under conditions of stable instrumental performance.     

The response of internal dynamics to hydrophobic core mutations in the SH3 domain from the Fyn tyrosine kinase

We have used (15)N- and (2)H-NMR spin relaxation experiments to study the response of backbone and side-chain dynamics when a leucine or valine is substituted for a completely buried phenylalanine residue in the SH3 domain from the Fyn tyrosine kinase. Several residues show differences in the time scales and temperature dependences of internal motions when data for the three proteins are compared. Changes were also observed in the magnitude of dynamics, with the valine, and to a lesser extent leucine mutant, showing enhanced flexibility compared to the wild-type (WT) protein. The motions of many of the same amide and methyl groups are affected by both mutations, identifying a set of loci where dynamics are sensitive to interactions involving the targeted side chain. These results show that contacts within the hydrophobic core affect many aspects of internal mobility throughout the Fyn SH3 domain.     

Side-chain dynamics of the SAP SH2 domain correlate with a binding hot spot and a region with conformational plasticity

X-linked lymphoproliferative disease is caused by mutations in the protein SAP, which consists almost entirely of a single SH2 domain. SAP interacts with the Tyr281 site of the T<-->B cell signaling protein SLAM via its SH2 domain. Interestingly, binding is not dependent on phosphorylation but does involve interactions with residues N-terminal to the Tyr. We have used 15N and 2H NMR relaxation experiments to investigate the motional properties of the SAP SH2 domain backbone amides and side-chain methyl groups in the free protein and complexes with phosphorylated and non-phosphorylated peptides derived from the Tyr281 site of SLAM. The most mobile methyl groups are in side-chains with large RMSD values between the three crystal structures of SAP, suggesting that fast time-scale dynamics in side-chains is associated with conformational plasticity. The backbone amides of two residues which interact with the C-terminal part of the peptides experience fast time-scale motions in the free SH2 domain that are quenched upon binding of either the phosphorylated or non-phosphorylated peptide. Of most importance, the mobility of methyl groups in and around the binding site for residues in the N-terminus of the peptide is significantly restricted in the complexes, underscoring the dominance of this interaction with SAP and demonstrating a correlation between changes in rapid side-chain motion upon binding with local binding energy.     

Efficient side-chain and backbone assignment in large proteins: Application to tGCN5

In determining the structure of large proteins by NMR, it would be desirable to obtain complete backbone, side-chain, and NOE assignments efficiently, with a minimum number of experiments and samples. Although new strategies have made backbone assignment highly efficient, side-chain assignment has remained more difficult. Faced with the task of assigning side-chains in a protein with poor relaxation properties, the Tetrahymena histone acetyltransferase tGCN5, we have developed an assignment strategy that would provide complete side-chain assignments in cases where fast ¹³C transverse relaxation causes HCCH-TOCSY experiments to fail. Using the strategy presented here, the majority of aliphatic side-chain proton and carbon resonances can be efficiently obtained using optimized H(CC-CO)NH-TOCSY and (H)C(C-CO)NH-TOCSY experiments on a partially deuterated protein sample. Assignments can be completed readily using additional information from a ¹³ C-dispersed NOESY-HSQC spectrum. Combination of these experiments with H(CC)NH-TOCSY and (H)C(C)NH-TOCSY may provide complete backbone and side-chain assignments for large proteins using only one or two samples.     

Conformational properties of trans Ac-Asn-Pro-Tyr-NHMe and trans Ac-Tyr-Pro-Asn-NHMe in dimethylsulfoxide and in water determined by multinuclear n.m.r. spectroscopy

Vicinal coupling constants between various nuclei provide backbone and side-chain conformational information for a series of asparagine- and tyrosine-containing peptides in DMSO and in H2O. By enriching Tyr of Ac-Asn-Pro-Tyr-NHMe with 15N, it has been possible to distinguish between the resonances of the two side-chain beta protons of Tyr. Analysis of the coupling constants in terms of the distributions of side-chain conformations in these peptides indicates that the addition of Asn to the Pro-Tyr sequence leads to a less random conformational distribution. When compared to the side-chain rotamer distribution of Ac-Asn-NHMe and Ac-Tyr-NHMe, particular Asn and Tyr side-chain conformations of Ac-Asn-Pro-Tyr-NHMe are stabilized in dimethylsulfoxide solution. The interaction(s) which stabilize a unique Tyr side-chain conformation of Ac-Asn-Pro-Tyr-NHMe in dimethylsulfoxide are not present in Ac-Ala-Pro-Tyr-NHMe and are unaffected by the addition of Val-Pro to the C-terminus of Asn-Pro-Tyr. In water, a preferential stabilization of one Asn side-chain conformation of Ac-Asn-Pro-Tyr-NHMe is also observed, while the Tyr side-chain rotamer distribution is similar to that of Ac-Tyr-NHMe. An interaction between the Asn side chain and the Pro-Tyr-NHMe backbone was previously shown to stabilize a beta-bend conformation at Pro-Tyr in water. Data are also presented for Ac-Tyr-Pro-Asn-NHMe, for which local interactions do not stabilize particular backbone conformations in dimethylsulfoxide or in water. The conformations of the peptides studied here are relatively insensitive to temperatures between 27 degrees and 62 degrees, both in dimethylsulfoxide and in water. The sequences Asn-Pro-Tyr and Tyr-Pro-Asn occur in ribonuclease A, and these tripeptides serve as models for the interactions involved in the folding of this protein.     

Dynamics of fd coat protein in the bacteriophage

The dynamics of the coat protein in fd bacteriophage are described with solid-state 15N and 2H NMR experiments. The virus particles and the coat protein subunits are immobile on the time scales of the 15N chemical shift anisotropy (10(3) Hz) and 2H quadrupole (10(6) Hz) interactions. Previously we have shown that the Trp-26 side chain is immobile, that the two Tyr and three Phe side chains undergo only rapid twofold jump motions about their C beta-C gamma bond axis [Gall, C. M., Cross, T. A., DiVerdi, J. A., & Opella, S. J. (1982) Proc. Natl. Acad. Sci. U.S.A. 79, 101-105], and that most of the backbone peptide linkages are highly constrained but do undergo rapid small amplitude motions [Cross, T. A., & Opella, S. J. (1982) J. Mol. Biol. 159, 543-549] in the coat protein subunits in the virus particles. In this paper, we demonstrate that the four N-terminal residues of the coat protein subunits are highly mobile, since both backbone and side-chain sites of these residues undergo large amplitude motions that are rapid on the time scales of the solid-state NMR experiments. In addition, the dynamics of the methyl-containing aliphatic residues Ala, Leu, Val, Thr, and Met are analyzed. Large amplitude jump motions are observed in nearly all of these side chains even though, with the exception of the N-terminal residue Ala-1, their backbone peptide linkages are highly constrained. The established information about the dynamics of the structural form of fd coat protein in the virus particle is summarized qualitatively.(ABSTRACT TRUNCATED AT 250 WORDS)     

Aromatic side-chain contributions to the far ultraviolet circular dichroism of peptides and proteins

The rotational strength of the L_a transition in phenylalanine and tyrosine side chains has been calculated for dipeptides with various backbone and side-chain conformations. Similar calculations have also been performed for tripeptides in the β-turn conformation with aromatic residues at the corners of the turn. The interaction of the aromatic ring with neighboring peptides generates rotational strengths in the L_a transition of the order of 0.1 Debye-Bohr magneton. When the preferred backbone and side-chain conformations are considered, it is found that the most probable conformations have positive L_a bonds. This result accounts for the observation that the N-acyl amino acid amides of L -Tyr and L Phe have positive L_a bands. It also suggests that, although other interactions may affect the numerical value and even the sign, there will be a significant positive contribution to the rotational strength of aromatic residues in globular proteins from nearest-neighbor interactions. Calculations on proteins of known conformation at the nearest-neighbor level confirm the tendency toward positive L_a contributions for Phe and Tyr residues. This contribution can be of the order of 10% of the observed CD even in proteins with rather strong amide contributions. In some proteins, such as the gene 5 protein from bacteriophage fd and many snake-venom toxins, side-chain contributions from Tyr and Trp residues manifest themselves as positive CD bands in the 225–250-nm region. The magnitude of the nearest-neighbor contributions and the trend toward positive contributions are consistent with the observation of such CD bands in globular proteins. No special stacking interaction among aromatic side chains needs to be invoked.