1H and 15N NMR assignment and solution structure of the SH3 domain of spectrin: Comparison of unrefined and refined structure sets with the crystal structure

The assignment of the ¹H and ¹⁵Nnuclear magnetic resonance spectra of the Src-homology region 3 domain ofchicken brain -spectrin has been obtained. A set of solutionstructures has been determined from distance and dihedral angle restraints,which provide a reasonable representation of the protein structure insolution, as evaluated by a principal component analysis of the globalpairwise root-mean-square deviation (rmsd) in a large set of structuresconsisting of the refined and unrefined solution structures and the crystalstructure. The solution structure is well defined, with a lower degree ofconvergence between the structures in the loop regions than in the secondarystructure elements. The average pairwise rmsd between the 15 refinedsolution structures is 0.71 ± 0.13 Å for the backbone atoms and1.43 ± 0.14 Å for all heavy atoms. The solution structure isbasically the same as the crystal structure. The average rmsd between the 15refined solution structures and the crystal structure is 0.76 Å forthe backbone atoms and 1.45 ± 0.09 Å for all heavy atoms. Thereare, however, small differences probably caused by intermolecular contactsin the crystal structure.     

Prediction of protein secondary structure from circular dichroism using theoretically derived spectra

Circular dichroism (CD) is a spectroscopic technique commonly used to investigate the structure of proteins. Major secondary structure types, alpha‐helices and beta‐strands, produce distinctive CD spectra. Thus, by comparing the CD spectrum of a protein of interest to a reference set consisting of CD spectra of proteins of known structure, predictive methods can estimate the secondary structure of the protein. Currently available methods, including K2D2, use such experimental CD reference sets, which are very small in size when compared to the number of tertiary structures available in the Protein Data Bank (PDB). Conversely, given a PDB structure, it is possible to predict a theoretical CD spectrum from it. The methodological framework for this calculation was established long ago but only recently a convenient implementation called DichroCalc has been developed. In this study, we set to determine whether theoretically derived spectra could be used as reference set for accurate CD based predictions of secondary structure. We used DichroCalc to calculate the theoretical CD spectra of a nonredundant set of structures representing most proteins in the PDB, and applied a straightforward approach for predicting protein secondary structure content using these theoretical CD spectra as reference set. We show that this method improves the predictions, particularly for the wavelength interval between 200 and 240 nm and for beta‐strand content. We have implemented this method, called K2D3, in a publicly accessible web server at http://www. ogic.ca/projects/k2d3 . Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Numbers and exchangeability with water of oxygen-17 atoms coupled to molybdenum (V) in different reduced forms of xanthine oxidase

The effect of using [17O]water (24-50% enriched) as solvent on the Mo(V) electron paramagnetic resonance spectra of different reduced forms of xanthine oxidase has been investigated. All the Mo(V) signals are affected. Procedures are described, based on the use of difference spectral techniques, that facilitate interpretation of such spectra. The number of coupled oxygen atoms may be determined by estimation of the fraction of the spectrum that remains unchanged by the isotope at a known enrichment. For a species having two coupled oxygen atoms, the use of two different isotope enrichments permits elimination from the difference spectra of the contribution of the two singly substituted species. From the application of these methods, it is concluded that not only the strength of the hyperfine coupling of oxygen ligands of molybdenum but also their number and their exchangeability with the solvent vary from one reduced form of the enzyme to another. The inhibited species from active xanthine oxidase has been studied in the most detail. It has two weakly coupled oxygen atoms [A(17O)av = 0.1-0.2 mT] that do not exchange with the solvent. A cyclic structure is proposed for this species in which two oxygen ligands of molybdenum are bonded to the carbon of the formaldehyde or other alcohol or aldehyde molecule that reacted in producing the signal. Structures of the other signal-giving species from active xanthine oxidase (Very Rapid and Rapid types 1 and 2) are discussed, as is corresponding information on species from the desulfo enzyme and from sulfite oxidase.     

A holistic approach for protein secondary structure estimation from infrared spectra in H(2)O solutions

We present an improved technique for estimating protein secondary structure content from amide I and amide III band infrared spectra. This technique combines the superposition of reference spectra of pure secondary structure elements with simultaneous aromatic side chain, water vapor, and solvent background subtraction. Previous attempts to generate structural reference spectra from a basis set of reference protein spectra have had limited success because of inaccuracies arising from sequential background subtractions and spectral normalization, arbitrary spectral band truncation, and attempted resolution of spectroscopically degenerate structure classes. We eliminated these inaccuracies by defining a single mathematical function for protein spectra, permitting all subtractions, normalizations, and amide band deconvolution steps to be performed simultaneously using a single optimization algorithm. This approach circumvents many of the problems associated with the sequential nature of previous methods, especially with regard to removing the subjectivity involved in each processing step. A key element of this technique was the calculation of reference spectra for ordered helix, unordered helix, sheet, turns, and unordered structures from a basis set of spectra of well-characterized proteins. Structural reference spectra were generated in the amide I and amide III bands, both of which have been shown to be sensitive to protein secondary structure content. We accurately account for overlaps between amide and nonamide regions and allow different structure types to have different extinction coefficients. The agreement between our structure estimates, for proteins both inside and outside the basis set, and the corresponding determinations from X-ray crystallography is good.     

Direct evidence of nitrogen coupling in the copper(II) complex of bovine serum albumin by S-band electron spin resonance technique

ESR spectra of the tight binding Cu(II) complex of bovine serum albumin (BSA) has been studied using S-band. At physiological pH, only one form of copper binding to BSA was detected from the ESR spectra. From previous X-band ESR spectra, nitrogen superhyperfine splittings were observable in the g perpendicular region; however, the resolution of the g parallel region was not sufficient to confirm the exact donor atoms of the complex. Using low-frequency ESR (2-4 GHz) at 77 K, we have resolved the nitrogen superhyperfine structure in the g parallel region. A computer simulation method has been developed for distinguishing between three and four nitrogen donor atoms. The Hyde-Froncisz theory of g and A strain broadening has been modified to use a field-swept calculation for the line shape. The observed intensity pattern and the computer simulation of such spectra positively confirm the structure of Cu(II) ion coordinated to four in-plane nitrogen atoms in frozen aqueous solutions of Cu(II)-BSA complexes at physiological pH. This is the first time that this binding site has been confirmed on the protein instead of a protein fragment or model compound. This work is another example of the usefulness of the S-band ESR technique for characterizing the metal-protein interactions when random variation in g factors cause line broadening in conventional X-band ESR spectra.     

Tertiary templates for proteins. Use of packing criteria in the enumeration of allowed sequences for different structural classes

We assume that each class of protein has a core structure that is defined by internal residues, and that the external, solvent-contacting residues contribute to the stability of the structure, are of primary importance to function, but do not determine the architecture of the core portions of the polypeptide chain. An algorithm has been developed to supply a list of permitted sequences of internal residues compatible with a known core structure. This list is referred to as the tertiary template for that structure. In general the positions in the template are not sequentially adjacent and are distributed throughout the polypeptide chain. The template is derived using the fixed positions for the main-chain and beta-carbon atoms in the test structure and selected stereochemical rules. The focus of this paper is on the use of two packing criteria: avoidance of steric overlap and complete filling of available space. The program also notes potential polar group interactions and disulfide bonds as well as possible burial of formal charges. Central to the algorithm is the side-chain rotamer library. In an update of earlier studies by others, we show that 17 of the 20 amino acids (omitting Met, Lys and Arg) can be represented adequately by 67 side-chain rotamers. A list of chi angles and their standard deviations is given. The newer, high-resolution, refined structures in the Brookhaven Protein Data Bank show similar mean chi values, but have much smaller deviations than those of earlier studies. This suggests that a rotamer library may be a better structural approximation than was previously thought. In using packing constraints, it has been found essential to include all hydrogen atoms specifically. The "unified atom" representation is not adequate. The permitted rotamer sequences are severely restricted by the main-chain plus beta-carbon atoms of the test structure. Further restriction is introduced if the full set of atoms of the external residues are held fixed, the full-chain model. The space-filling requirement has a major role in restricting the template lists. The preliminary tests reported here make it appear likely that templates prepared from the currently known core structures will be able to discriminate between these structures. The templates should thus be useful in deciding whether a sequence of unknown tertiary structure fits any of the known core classes and, if a fit is found, how the sequence should be aligned in three dimensions to fit the core of that class.(ABSTRACT TRUNCATED AT 400 WORDS)     

Application of molecular dynamics with interproton distance restraints to three-dimensional protein structure determination. A model study of crambin

The applicability of restrained molecular dynamics for the determination of three-dimensional protein structures on the basis of short interproton distances (less than 4 A) that can be realistically determined from nuclear magnetic resonance measurements in solution is assessed. The model system used is the 1.2 A resolution crystal structure of the 46 residue protein crambin, from which a set of 240 approximate distance restraints, divided into three ranges (2.5 +/- 0.5, 3.0+0.5(-1.0) and 4 +/- 1 A), is derived. This interproton distance set comprises 159 short-range ([i-j] less than or equal to 5) and 56 ([i-j] greater than 5) long-range inter-residue distances and 25 intra-residue distances. Restrained molecular dynamics are carried out using a number of different protocols starting from two initial structures: a completely extended beta-strand; and an extended structure with two alpha-helices in the same positions as in the crystal structure (residues 7 to 19, and 23 to 30) and all other residues in the form of extended beta-strands. The root-mean-square (r.m.s.) atomic differences between these two initial structures and the crystal structure are 43 A and 23 A, respectively. It is shown that, provided protocols are used that permit the secondary structure elements to form at least partially prior to folding into a tertiary structure, convergence to the correct final structure, both globally and locally, is achieved. The r.m.s. atomic differences between the converged restrained dynamics structures and the crystal structure range from 1.5 to 2.2 A for the backbone atoms and from 2.0 to 2.8 A for all atoms. The r.m.s. atomic difference between the X-ray structure and the structure obtained by first averaging the co-ordinates of the converged restrained dynamics structures is even smaller: 1.0 A for the backbone atoms and 1.6 A for all atoms. These results provide a measure with which to judge future experimental results on proteins whose crystal structures are unknown. In addition, from an examination of the dynamics trajectories, it is shown that the convergence pathways followed by the various simulations are different.     

Towards developing a protein infrared spectra databank (PISD) for proteomics research

Fourier transform infrared (FTIR) spectroscopy is an attractive tool for proteomics research as it can be used to rapidly characterize protein secondary structure in aqueous solution. Recently, a number of secondary structure prediction methods based on reference sets of FTIR spectra from proteins with known structure from X-ray crystallography have been suggested. These prediction methods, often referred to as pattern recognition based approaches, demonstrated good prediction accuracy using some error measure, e.g., the standard error of prediction (SEP). However, to avoid possible adverse effects from differences in recording, the analysis has been mostly based on reference sets of FTIR spectra from proteins recorded in one laboratory only. As a result, these studies were based on reference sets of FTIR spectra from a limited number of proteins. Pattern recognition based approaches, however, rely on reference sets of FTIR spectra from as many proteins as possible representing all possible band shape variation to be related to the diversity of protein structural classes. Hence, if we want to build reliable pattern recognition based systems to support proteomics research, which are capable of making good predictions from spectral data of any unknown protein, one common goal should be to build a comprehensive protein infrared spectra databank (PISD) containing FTIR spectra of proteins of known structure. We have started the process of developing a comprehensive PISD composed of spectra recorded in different laboratories. As part of this work, here we investigate possible effects on prediction accuracy achieved by a neural network analysis when using reference sets composed of FTIR spectra from different laboratories. Surprisingly low magnitude of difference in SEPs throughout all our experiments suggests that FTIR spectra recorded in different laboratories may be safely combined into one reference set with only minor deterioration of prediction accuracy in the worst case.     

The pH-dependent conformational states of kyotorphin: a constant-pH molecular dynamics study

An extensive conformational study of the analgesic dipeptide kyotorphin (L-Tyr-L-Arg) at different pH values was performed using a constant-pH molecular dynamics method. This dipeptide showed a remarkable pH-dependent conformational variety. The protonation of the N-terminal amine was identified as a key element in the transition between the more extended and the more packed conformational states, as monitored by the dihedral angle defined by the atoms 1Cbeta-1Calpha-2Calpha-2Cbeta. The principal-component analysis of kyotorphin identified two major conformational populations (the extended trans and the packed cis) together with conformations that occur exclusively at extreme pH values. Other, less stable conformations were also identified, which help us to understand the transitions between the predominant populations. The fitting of kyotorphin's conformational space to the structure of morphine resulted in a set of conformers that were able to fulfill most of the constraints for the mu-receptor. These results suggest that there may be strong similarities between the kyotorphin receptor and the structural family of opioid receptors.     

Determination of the three-dimensional solution structure of barnase using nuclear magnetic resonance spectroscopy

The solution conformation of the ribonuclease barnase has been determined by using 1H nuclear magnetic resonance (NMR) spectroscopy. The 20 structures were calculated by using 853 interproton distance restraints obtained from analyses of two-dimensional nuclear Overhauser spectra, 72 phi and 53 chi 1 torsion angle restraints, and 17 hydrogen-bond distance restraints. The calculated structures contain two alpha-helices (residues 6-18 and 26-34) and a five-stranded antiparallel beta-sheet (residues 50-55, 70-75, 85-91, 94-101, and 105-108). The core of the protein is formed by the packing of one of the alpha-helices (residues 6-18) onto the beta-sheet. The average RMS deviation between the calculated structures and the mean structure is 1.11 A for the backbone atoms and 1.75 A for all atoms. The protein is least well-defined in the N-terminal region and in three large loops. When these regions are excluded, the average RMS deviation between the calculated structures and the mean structure for residues 5-34, 50-56, 71-76, 85-109 is 0.62 A for the backbone atoms and 1.0 A for all atoms. The NMR-derived structure has been compared with the crystal structure of barnase [Mauguen et al. (1982) Nature (London) 297, 162-164].     

Bacteriochlorophyll a-protein interactions in a complex from Prosthecochloris aestuarii. A Resonance Raman study

Raman spectra of bacteriochlorophyll a (BChl) bound to the soluble protein complex from Prosthecochloris aestuarii have been obtained at low temperature, using the resonance effect on their Q_x for Soret electronic bands. These spectra show that the acetyl carbonyls of at least four of the seven molecules bound to the monomer subunit of the complex and the ketone carbonyls of at least five of them are oriented close to the mean plane of the conjugated part of the dihydrophorbin macrocycle. Up to three bacteriochlorophyll molecules may have their ketone carbonyls free from hydrogen-bonding and up to two may have their acetyl carbonyls similarly free. Several of the binding sites of the remaining conjugated carbonyls are probably the same as those binding the conjugated carbonyls of bacteriochlorophyll (and of bacteriopheophytin) in reaction centers and in antenna structures of purple bacteria and as those binding chlorophyll in the antenna of higher plants and algae. The present resonance Raman spectra confirm that the magnesium atoms of most of the seven bacteriochlorophylls are pentacoordinated. They also show that polarisation effects from their local environments induce changes in the groundstate structures of the dihydrophorbin skeletons of the complexed molecules with respect to those of isolated, monomeric bacteriochlorophyll. These changes are quasi-identical for the seven molecules. These environmental effects predominate over any structural change brought about by intermolecular bonding of the conjugated carbonyls or of the magnesium atoms. The dihydrophorbin rings of the seven molecules thus appear to be immersed in a nearly homogeneous medium of low permittivity, although specific van der Waals interactions may polarise the free carbonyls to quite different extents. The possible implications of these observations on the interpretation of the electronic spectrum of the set of complexed bacteriochlorophylls are discussed.     

Electron spin resonance study of the copper(II) complexes of human and dog serum albumins and some peptide analogs

Electron spin resonance spectra of the first Cu(II) complexes of human serum albumin, dog serum albumin, l-aspartyl-l-histidine N-methylamide and glycyl-glycyl-l-histidine N-methylamide have been studied using isotopically pure ⁶⁵Cu in its chloride form. At 77° K, the esr spectra of Cu(II) complex of human serum albumin exhibited only one form of esr signal between pH 6.5 and 11. No intermediate forms were detected. The presence of an equally spaced nine-line superhyperfine structure with spacing ～15 G indicated considerable covalent bonding between Cu(II) and four nitrogen atoms derived from the protein. The esr spectrum form of Cu(II) bound to human serum albumin detected at neutral pH would be consistent with the participation of four nitrogens from the α-NH₂ group, two peptide groups, and the imidazole group of a histidine residue. In contrast, the esr spectra of Cu(II)-dog serum albumin complex showed a transition from a low pH form to a high pH form as the pH was increased to 9.5. These spectral changes were found to be reversible upon lowering the pH. Ligand superhyperfine splittings in the low pH form of the esr signal of Cu(II)-dog albumin were not resolved. The distinct pH dependence of the esr signals observed in human and dog serum albumin complexes could be correlated to their respective optical spectra changes as a function of pH. At room temperature and in the pH range between 6 and 11, the esr spectra of Cu(II) complexes of l-aspartyl-l-alanyl-l-histidine N-methylamide and glycyl-glycyl-l-histidine N-methylamide exhibited a well-resolved nine-line superhyperfine structure indicating metal coordination with four equivalent nitrogen atoms of peptide.     

Determination of the complete three-dimensional structure of the trypsin inhibitor from squash seeds in aqueous solution by nuclear magnetic resonance and a combination of distance geometry and dynamical simulated annealing

The complete three-dimensional structure of the trypsin inhibitor from seeds of the squash Cucurbita maxima in aqueous solution was determined on the basis of 324 interproton distance constraints, 80 non-nuclear Overhauser effect distances, and 22 hydrogen-bonding constraints, supplemented by 27 phi backbone angle constraints derived from nuclear magnetic resonance measurements. The nuclear magnetic resonance input data were converted to the distance constraints in a semiquantitative manner after a sequence specific assignment of 1H spectra was obtained using two-dimensional nuclear magnetic resonance techniques. Stereospecific assignments were obtained for 17 of the 48 prochiral centers of the squash trypsin inhibitor using the floating chirality assignment introduced at the dynamical simulated annealing stage of the calculations. A total of 34 structures calculated by a hybrid distance geometry-dynamical simulated annealing method exhibit well-defined positions for both backbone and side-chain atoms. The average atomic root-mean-square difference between the individual structures and the minimized mean structure is 0.35(+/- 0.08) A for the backbone atoms and 0.89(+/- 0.17) A for all heavy atoms. The precision of the structure determination is discussed and correlated to the experimental input data.     

Peakmatch: a simple and robust method for peak list matching

Peak lists are commonly used in NMR as input data for various software tools such as automatic assignment and structure calculation programs. Inconsistencies of chemical shift referencing among different peak lists or between peak and chemical shift lists can cause severe problems during peak assignment. Here we present a simple and robust tool to achieve self-consistency of the chemical shift referencing among a set of peak lists. The Peakmatch algorithm matches a set of peak lists to a specified reference peak list, neither of which have to be assigned. The chemical shift referencing offset between two peak lists is determined by optimizing an assignment-free match score function using either a complete grid search or downhill simplex optimization. It is shown that peak lists from many different types of spectra can be matched reliably as long as they contain at least two corresponding dimensions. Using a simulated peak list, the Peakmatch algorithm can also be used to obtain the optimal agreement between a chemical shift list and experimental peak lists. Combining these features makes Peakmatch a useful tool that can be applied routinely before automatic assignment or structure calculation in order to obtain an optimized input data set.     

The structural evolution of cobra venom cytotoxins

Summary In order to analyze the evolutionary behavior of the cobra venom cytotoxins, their probable tertiary structure was predicted using computer graphics. The 41 amino acid sequences known show that the major evolutionary changes have taken place in two particularly exposed areas of the molecular surface. In each area, neighboring residue positions seem to have evolved interdependently, but there is no obvious interdependence between the two areas. Indeed, the relative evolution of these two areas prompts a subdivision of the sequence set into four groups. According to the known cytotoxin circular dichroism spectra, one of these four groups could be characterized by a difference in molecular secondary structure. Sine the two variable areas have functional associations, it is suggested that their evolution may be governed by a target with several similar binding sites.     

An algorithm for constraint-based structural template matching: application to 3D templates with statistical analysis

MOTIVATION: Structural templates consisting of a few atoms in a specific geometric conformation provide a powerful tool for studying the relationship between protein structure and function. Current methods for template searching constrain template syntax and semantics by their design. Hence there is a need for a more flexible core algorithm upon which to build more sophisticated tools. Statistical analysis of structural similarity is still in its infancy when compared with its analogue in sequence alignment. In the context of template matching, there is an urgent need for normalization of scores so that results from templates with differing sensitivity may be compared directly. RESULTS: We introduce Jess, a fast and flexible algorithm for searching protein structures for small groups of atoms under arbitrary constraints on geometry and chemistry. We apply the algorithm to a set of manually derived enzyme active site templates, and derive an empirical measure for estimating the relative significance of hits encountered using differing templates.     

Empirical solvation models can be used to differentiate native from near-native conformations of bovine pancreatic trypsin inhibitor.

Several hydration models for peptides and proteins based on solvent accessible surface area have been proposed previously. We have evaluated some of these models as well as four new ones in the context of near-native conformations of a protein. In addition, we propose an empirical site-site distance-dependent correction that can be used in conjunction with any of these models. The set of near-native structures consisted of 39 conformations of bovine pancreatic trypsin inhibitor (BPTI) each of which was a local minimum of an empirical energy function (ECEPP) in the absence of solvent. Root-mean-square (rms) deviations from the crystallographically determined structure were in the following ranges: 1.06-1.94 A for all heavy atoms, 0.77-1.36 A for all backbone heavy atoms, 0.68-1.33 A for all alpha-carbon atoms, and 1.41-2.72 A for all side-chain heavy atoms. We have found that there is considerable variation among the solvent models when evaluated in terms of concordance between the solvation free energy and the rms deviations from the crystallographically determined conformation. The solvation model for which the best concordance (0.939) with the rms deviations of the C alpha atoms was found was derived from NMR coupling constants of peptides in water combined with an exponential site-site distance dependence of the potential of mean force. Our results indicate that solvation free energy parameters derived from nonpeptide free energies of hydration may not be transferrable to peptides. Parameters derived from peptide and protein data may be more applicable to conformational analysis of proteins. A general approach to derive parameters for free energy of hydration from ensemble-averaged properties of peptides in solution is described.     

A study of the conformations of valinomycin in solution phase

Vibrational absorption and vibrational circular dichroism (VCD) spectra of valinomycin are measured, in different solvents, in the ester and amide carbonyl stretching regions. The influence of cations, namely Li(+), Na(+), K(+), and Cs(+), in methanol-d(4) solvent is also investigated. Ab initio quantum mechanical calculations using density functional theory and 6-31G* basis set are used to predict the absorption and VCD spectra. A bracelet-type structure for valinomycin that reproduces the experimental absorption and VCD spectra in inert solvents is identified. For the structure of valinomycin in polar solvents, a propeller-type structure was optimized, but further investigations are required to confirm this structure. A symmetric octahedral environment for the ester carbonyl groups in the valinomycin-K(+) complex is supported by the experimental VCD spectra. The results obtained in the present study demonstrate that even for large macrocyclic peptides, such as valinomycin, VCD can be used as an independent structural tool for the study of conformations in solution.     

Residue-level global and local ensemble-ensemble comparisons of protein domains

Many methods of protein structure generation such as NMR-based solution structure determination and template-based modeling do not produce a single model, but an ensemble of models consistent with the available information. Current strategies for comparing ensembles lose information because they use only a single representative structure. Here, we describe the ENSEMBLATOR and its novel strategy to directly compare two ensembles containing the same atoms to identify significant global and local backbone differences between them on per-atom and per-residue levels, respectively. The ENSEMBLATOR has four components: eePREP (ee for ensemble-ensemble), which selects atoms common to all models; eeCORE, which identifies atoms belonging to a cutoff-distance dependent common core; eeGLOBAL, which globally superimposes all models using the defined core atoms and calculates for each atom the two intraensemble variations, the interensemble variation, and the closest approach of members of the two ensembles; and eeLOCAL, which performs a local overlay of each dipeptide and, using a novel measure of local backbone similarity, reports the same four variations as eeGLOBAL. The combination of eeGLOBAL and eeLOCAL analyses identifies the most significant differences between ensembles. We illustrate the ENSEMBLATOR''s capabilities by showing how using it to analyze NMR ensembles and to compare NMR ensembles with crystal structures provides novel insights compared to published studies. One of these studies leads us to suggest that a “consistency check” of NMR-derived ensembles may be a useful analysis step for NMR-based structure determinations in general. The ENSEMBLATOR 1.0 is available as a first generation tool to carry out ensemble-ensemble comparisons.     

Nitrogen chemical structure in DNA and related molecules by X-ray absorption spectroscopy.

The electronic environment of nitrogen in nucleic acid bases, nucleotides, polynucleotides and DNA has been studied, for the first time using X-Ray Absorption Near-Edge Spectroscopy (XANES). Generally, the spectra of these complex molecules consist of low energy bands corresponding to 1s-->pi* transitions and high energy bands corresponding to 1s-->sigma* transition, as illustrated using several nitrogen model compounds. The 1s-->pi* transitions show particular sensitivity to the chemical environment of the nitrogen. Oxygen substitution on ring carbon atoms generally results in a significant blue shift of the lowest 1s-->pi* bands while halogen substitution results in a small blue shift. These observations illustrate the significance of the disturbance of the aromatic ring system produced by exocyclic carbonyl groups. Direct substitution on the nitrogen frequently results in significant spectral perturbations. Differences between the spectra of the polynucleotides and the sums of spectra of the individual nucleotides point to the effects of hydrogen-bonding in complementary double-helix structures. The XANES spectrum of a DNA sample with a known ratio of the polynucleotides is equivalent to the weighted sum of the spectra of individual polynucleotides, indicating that the difference in base stacking interactions produces negligible spectral effects. The variability of nitrogen K-edge spectra in these samples and in protein may be useful for chemically specific imaging using X-ray microscopes.