Sequential Connectivities Research Articles

NMR identification and characterization of the flexible regions in the 160 kDa molten globule-like aggregate of barstar at low pH.

Barstar is known to form a molten globule-like A form below pH 4. This form exists as a soluble aggregate of 16 monomeric subunits, and appears to remain homogeneous in solution for at least two weeks. Here, structural characterization by NMR of the flexible regions in the A form of barstar has been carried out at pH 2.7 and 25 degrees C. Significantly, the A form appears to be a symmetrical aggregate. Using the recently described fast assignment strategy from HNN and HN(C)N spectra, along with the standard triple resonance and three-dimensional NMR experiments, the flexible segment of the aggregate has been identified to belong largely to the N-terminal end of the polypeptide chain; sequential connectivities were obtained for the first 20 residues (except two) from these experiments. This segment is free in each of the monomeric subunits, and does not form a part of the aggregated core of the A form. The secondary chemical shifts of these residues suggest propensity toward an extended structure. Their (3)J(HN,H)(alpha) coupling constants have values corresponding to those in a random coil structure. However, a few medium-range NOEs, some of them involving side chain atoms, are observed between some residues in this segment. The lowered temperature coefficients of the H(N) chemical shifts compared to random coil values indicate possibilities of some hydrogen bonding in this region. Analysis of the (15)N relaxation parameters and reduced spectral density functions, in particular the negative values of heteronuclear NOEs, indicates large-amplitude high-frequency motions in the N-terminal segments; the first three residues show more negative NOEs than the others. The (15)N transverse relaxation rates and the J(0) spectral density values for residues Ser12 and Ser69 are significantly larger than for the rest, indicating some microsecond to millisecond time scale conformational exchange contributions to the relaxation of these residues. Taken all together, the data suggest that the A form of barstar is an aggregate with a rigid core, but with the N-terminal 20 residues of each of the monomeric subunits, in a highly dynamic random coil conformation which shows transient local ordering of structure. The N-terminal segment, anchored to the aggregated core, exhibits free-flight motion.

Assignment Strategy for Proteins with Known Structure

In protein NMR the assignment of nuclear spin resonances is a prerequisite for all subsequent applications, such as studies of ligand binding, protein-DNA interactions, and dynamics. Resonance assignment is a time consuming step even when the 3D x-ray structure of the protein is available. A new strategy is presented to solve the “inverse” assignment problem, which is the determination of the NMR resonance assignment from a known 3D protein structure. The protocol employs NMR data in the form of residual dipolar couplings and chemical shifts, while it does not require any sequential NMR connectivity information. The assignment problem is mathematically formulated in terms of a weighted matching problem that can be computationally efficiently solved by a combinatorial optimization algorithm. The protocol is applied to ubiquitin using two or three residual dipolar couplings per amino acid measured in Pfl phage medium together with chemical shift information. The algorithm yields for more than 90% of the protein backbone resonances the correct assignment.

Protein structure elucidation from NMR proton densities.

The NMR-generated foc proton density affords a template to which the molecule has to be fitted to derive the structure. Here we present a computational protocol that achieves this goal. H(N) atoms are readily recognizable from (1)H/(2)H exchange or (1)H/(15)N heteronuclear single quantum correlation (HSQC) experiments. The primary structure is threaded through the unassigned foc by leapfrogging along peptidyl amide H(N)s and the connected H(alpha)s. Via a Bayesian approach, the probabilities of the sequential connectivity hypotheses are inferred from likelihoods of H(N)/H(N), H(N)/H(alpha), and H(alpha)/H(alpha) interatomic distances as well as (1)H NMR chemical shifts, both derived from public databases. Once the polypeptide sequence is identified, directionality becomes established, and the foc N and C termini are recognized. After a similar procedure, side chain H atoms are found, including discriminated cis/trans proline loci. The folded structure then is derived via a direct molecular dynamics embedding into mirror image-related representations of the foc and selected according to a lowest energy criterion. The method was applied to foc densities calculated for two protein domains, col 2 and kringle 2. The obtained structures are within 1.0-1.5 A (backbone heavy atoms) and 1.5-2.0 A (all heavy atoms) rms deviations from reported x-ray and/or NMR structures.

The Structural Basis of Compstatin Activity Examined by Structure-Function-based Design of Peptide Analogs and NMR

We have previously identified compstatin, a 13-residue cyclic peptide, that inhibits complement activation by binding to C3 and preventing C3 cleavage to C3a and C3b. The structure of compstatin consists of a disulfide bridge and a type I beta-turn located at opposite sides to each other. The disulfide bridge is part of a hydrophobic cluster, and the beta-turn is part of a polar surface. We present the design of compstatin analogs in which we have introduced a series of perturbations in key structural elements of their parent peptide, compstatin. We have examined the consistency of the structures of the designed analogs compared with compstatin using NMR, and we have used the resulting structural information to make structure-complement inhibitory activity correlations. We propose the following. 1) Even in the absence of the disulfide bridge, a linear analog has a propensity for structure formation consistent with a turn of a 3(10)-helix or a beta-turn. 2) The type I beta-turn is a necessary but not a sufficient condition for activity. 3) Our substitutions outside the type I beta-turn of compstatin have altered the turn population but not the turn structure. 4) Flexibility of the beta-turn is essential for activity. 5) The type I beta-turn introduces reversibility and sufficiently separates the two sides of the peptide, whereas the disulfide bridge prevents the termini from drifting apart, thus aiding in the formation of the hydrophobic cluster. 6) The hydrophobic cluster at the linked termini is involved in binding to C3 and activity but alone is not sufficient for activity. 7) beta-Turn residues Gln(5) (Asn(5))-Asp(6)-Trp(7)(Phe(7))-Gly(8) are specific for the turn formation, but only Gln(5)(Asn(5))-Asp(6)-Trp(7)-Gly(8) residues are specific for activity. 8) Trp(7) is likely to be involved in direct interaction with C3, possibly through the formation of a hydrogen bond. Finally we propose a binding model for the C3-compstatin complex.

Protein sequential resonance assignments by combinatorial enumeration using 13C alpha chemical shifts and their (i, i-1) sequential connectivities.

The need for the structural characterization of proteins on a genomic scale has brought with it demands for new technology to speed the structure determination process. In NMR, one bottleneck is the sequential assignment of backbone resonances. In this paper, we explore the computational complexity of the sequential assignment problem using only 13C alpha chemical shift data and C alpha (i, i - 1) sequential connectivity information, all of which can potentially be obtained from a single three-dimensional NMR spectrum. Although it is generally believed that there is too much ambiguity in such data to provide sufficient information for sequential assignment, we show that a straightforward combinatorial search algorithm can be used to find correct and unambiguous sequential assignments in a reasonable amount of CPU time for small proteins (approximately 80 residues or smaller) when there is little missing data. The deleterious effect of missing or spurious peaks and the dependence on match tolerances is also explored. This simple algorithm could be used as part of a semi-automated, interactive assignment procedure, e.g., to test partial manually determined solutions for uniqueness and to extend these solutions.

Intercalation of the (1R,2S,3R,4S)-N6-[1-(1,2,3,4-tetrahydro-2,3,4-trihydroxybenz[a]anthracenyl)]-2'-deoxyadenosyl adduct in the N-ras codon 61 sequence: DNA sequence effects.

The structure of the bay region (1R,2S,3R,4S)-N6-[1-(1,2,3,4-tetrahydro-2,3,4-trihydroxybenz[a]anthracenyl)]-2'-deoxyadenosyl adduct at X(7) of 5'-d(CGGACAXGAAG)-3'.5'-d(CTTCTTGTCCG)-3', incorporating codons 60, 61 (underlined), and 62 of the human N-ras protooncogene, was determined by NMR. This was the bay region benz[a]anthracene RSRS (61,3) adduct. The BA moiety intercalated above the 5'-face of the modified base pair. NOE connectivities between imino protons were disrupted at T16 and T17. Large chemical shifts at the lesion site were consistent with ring current shielding arising from the BA moiety. A large chemical shift dispersion was observed for the BA aromatic protons. An increased rise of 8.17 A was observed between base pairs A6 x T17 and X7 x T(16). The PAH moiety stacked with the purine ring of A6, the 5'-neighbor nucleotide. This resulted in buckling of the 5'-neighbor A6 x T17 base pair, evidenced by exchange broadening for the T17 imino resonance. It also interrupted sequential NOE connectivities between nucleotides C5 and A6. The A6 deoxyribose ring showed an increased percentage of the C3'-endo conformation. This differed from the bay region BA RSRS (61,2) adduct, in which the lesion was located at position X6 [Li, Z., Mao, H., Kim, H.-Y., Tamura, P. J., Harris, C. M., Harris, T. M., and Stone, M. P. (1999) Biochemistry 38, 2969-2981], but was similar to the benzo[a]pyrene BP SRSR (61,3) adduct [Zegar I. S., Chary, P., Jabil, R. J., Tamura, P. J., Johansen, T. N., Lloyd, R. S., Harris, C. M., Harris, T. M., and Stone, M. P. (1998) Biochemistry 37, 16516-16528]. The altered sugar pseudorotation at A6 appears to be common to both bay region BA RSRS (61,3) and BP SRSR (61,3) adducts. It could not be discerned if the C3'-endo conformation at A6 in the BA RSRS (61,3) adduct altered base pairing geometry at X7 x T16, as compared to the C2'-endo conformation. The structural studies suggest that the mutational spectrum of this adduct may be more complex than that of the BA RSRS (61,2) adduct.

Cation-dependent structural features of beta-casein-(1-25).

Complete sequence-specific, proton-resonance assignments have been determined for the calcium phosphate-stabilizing tryptic peptide beta-casein-(1-25) containing the phosphorylated sequence motif Ser(P)(17)-Ser(P)-Ser(P)-Glu-Glu(21). Spectra of the peptide have been recorded, in separate experiments, in the presence of excess ammonium ions, sodium ions and calcium ions, and of the dephosphorylated peptide in the presence of excess sodium ions. We observed significant changes to chemical shifts for backbone and side-chain resonances that were dependent upon the nature of the cation present. Medium-range nuclear Overhauser effect (nOe) enhancements, characteristic of small structured regions in the peptide, were observed and also found to be cation dependent. The secondary structure of the peptide was characterized by sequential and medium-range (i, i+2/3/4, which denotes an interaction between residue i and residue i+2, i+3 or i+4 in the peptide) nOe connectivities, and Halpha chemical shifts. Four structured regions were identified in the calcium-bound peptide: residues Arg(1) to Glu(4) were involved in a loop-type structure, and residues Val(8) to Glu(11), Ser(P)(17) to Glu(20) and Glu(21) to Thr(24) were implicated in beta-turn conformations. Comparison of the patterns of medium-range nOe connectivities in beta-casein-(1-25) with those in alpha(S1)-casein-(59-79) suggest that the two peptides have distinctly different conformations in the presence of calcium ions, despite having a high degree of sequential and functional similarity.

Cation-dependent structural features of β-casein-(1–25)

Complete sequence-specific, proton-resonance assignments have been determined for the calcium phosphate-stabilizing tryptic peptide β-casein-(1–25) containing the phosphorylated sequence motif Ser(P)17-Ser(P)-Ser(P)-Glu-Glu21. Spectra of the peptide have been recorded, in separate experiments, in the presence of excess ammonium ions, sodium ions and calcium ions, and of the dephosphorylated peptide in the presence of excess sodium ions. We observed significant changes to chemical shifts for backbone and side-chain resonances that were dependent upon the nature of the cation present. Medium-range nuclear Overhauser effect (nOe) enhancements, characteristic of small structured regions in the peptide, were observed and also found to be cation dependent. The secondary structure of the peptide was characterized by sequential and medium-range (i, i+2/3/4, which denotes an interaction between residue i and residue i+2, i+3 or i+4 in the peptide) nOe connectivities, and Hα chemical shifts. Four structured regions were identified in the calcium-bound peptide: residues Arg1 to Glu4 were involved in a loop-type structure, and residues Val8 to Glu11, Ser(P)17 to Glu20 and Glu21 to Thr24 were implicated in β-turn conformations. Comparison of the patterns of medium-range nOe connectivities in β-casein-(1–25) with those in αS1-casein-(59–79) suggest that the two peptides have distinctly different conformations in the presence of calcium ions, despite having a high degree of sequential and functional similarity.

NMR Identification of Heavy Metal-Binding Sites in a Synthetic Zinc Finger Peptide: Toxicological Implications for the Interactions of Xenobiotic Metals with Zinc Finger Proteins

Lead (Pb), mercury (Hg), and cadmium (Cd) are toxic and interfere with protein metal-binding sites. The Cys2/His2 zinc finger is a structural motif required for sequence-specific DNA binding and is present in zinc finger transcription factors (ZFP): Sp1, Egr-1, and TFIIIA. Neurotoxic studies have shown that heavy metals directly inhibit the DNA binding of ZFP and result in adverse cellular effects. Recently, we demonstrated the ability of heavy metals to alter the DNA binding of a synthetic Cys2/His2 finger peptide (Razmiafshari and Zawia, Toxicol. Appl. Pharmacol. 166, 1–12, 2000). To determine the precise site of interactions between heavy metals and this protein domain, Pb, Hg, Cd, and Ca were reconstituted with the synthetic apopeptide and studied by one- and two-dimensional NMR spectroscopy. In the presence of Zn, Cd, Hg, and Pb, but not Ca, distinct peptide NMR signal changes in the aliphatic region were observed and attributed to metal–cystiene interactions. However, chemical shifts indicative of metal–histidine binding were elicited by all the metals in the peptide's aromatic region. Chemical shift assignments and sequential connectivity were established in the presence and absence of Zn, Pb, and Ca through TOCSY and NOESY spectra. Cysteine and histidine residues showed a distinct change in their amide and beta resonances in the presence of Zn and Pb, suggesting the metal-ligand binding sites were near these residues. However, Ca led to no significant spectral changes in these regions, suggesting that it is not actively involved in the binding site. These studies reveal this structure as a mediator of metal-induced alterations in protein function.

Conformational analysis of the nonapeptide leuprorelin using NMR and molecular modeling

The nonapeptide Leuprorelin, one of the LHRH agonists, was studied by means of 2D nuclear magnetic resonance spectroscopy and molecular modeling. NOESY spectra in aqueous/deuterated methanol solution (50%H2O/CD3OD) at low temperature (268 K) revealed short-range nOe connectivities (i, i+1), characteristic of flexibility of the molecule. The HN–HN sequential connectivities observed provide evidence that the sequence has the propensity to form a bend involving residues 5 and 6 and the N-terminal segment. The α-proton chemical shifts compared to random coil and additional data from the amide proton temperature coefficients support this assumption. One long-range nOe cross peak between H2α–HNEth is indicative of proximity between C- and N-termini.

NMR Assignment and Secondary Structure Determination of an Octameric 110 kDa Protein Using TROSY in Triple Resonance Experiments

TROSY-type triple resonance experiments with the uniformly 2H,13C,15N-labeled 7,8-dihydroneopterin aldolase (DHNA) from Staphylococcus aureus, which is a symmetric homooctamer protein of molecular mass 110 kDa, showed 20-fold to 50-fold sensitivity gains when compared to the corresponding conventional triple resonance NMR experiments. On this basis, sequential connectivities could be established for nearly all pairs of neighboring residues in DHNA. TROSY-type nuclear Overhauser enhancement spectroscopy yielded additional data to close the remaining gaps in the sequential assignment, and provided supplementary information on the secondary structure. Complete sequence-specific assignments of the 121-residue polypeptide chain in this 110 kDa octamer could thus be obtained in aqueous solution at 20 °C, and the regular secondary structures in the solution conformation were found to coincide nearly identically with those in the crystal structure of the DHNA octamer.

Protein structure alignment using a genetic algorithm.

We have developed a novel, fully automatic method for aligning the three-dimensional structures of two proteins. The basic approach is to first align the proteins' secondary structure elements and then extend the alignment to include any equivalent residues found in loops or turns. The initial secondary structure element alignment is determined by a genetic algorithm. After refinement of the secondary structure element alignment, the protein backbones are superposed and a search is performed to identify any additional equivalent residues in a convergent process. Alignments are evaluated using intramolecular distance matrices. Alignments can be performed with or without sequential connectivity constraints. We have applied the method to proteins from several well-studied families: globins, immunoglobulins, serine proteases, dihydrofolate reductases, and DNA methyltransferases. Agreement with manually curated alignments is excellent. A web-based server and additional supporting information are available at http://engpub1.bu.edu/-josephs.

Sequence-specific NMR assignment of proteins by global fragment mapping with the program MAPPER.

A new program, MAPPER, for semiautomatic sequence-specific NMR assignment in proteins is introduced. The program uses an input of short fragments of sequentially neighboring residues, which have been assembled based on sequential NMR connectivities and for which either the 13C(alpha) and 13C(beta) chemical shifts or data on the amino acid type from other sources are known. MAPPER then performs an exhaustive search for self-consistent simultaneous mappings of all these fragments onto the protein sequence. Compared to using only the individual mappings of the spectroscopically connected fragments, the global mapping adds a powerful new constraint, which results in resolving many otherwise intractable ambiguities. In an initial application, virtually complete sequence-specific assignments were obtained for a 110 kDa homooctameric protein, 7,8-dihydroneopterin aldolase from Staphylococcus aureus.

A novel NMR experiment for the sequential assignment of proline residues and proline stretches in 13C/15N-labeled proteins.

A new pulse sequence is described for the sequential assignment of proline residues in 13C/15N-labeled proteins by correlating C delta and C alpha chemical shifts of proline residues with the H alpha chemical shift of the preceding residue. Notably, the experiment can provide the sequential connectivities in poly-proline stretches, which cannot be determined using standard triple resonance experiments. Excellent solvent suppression is achieved by coherence selection via a heteronuclear gradient echo. The new pulse sequence has been successfully applied to the 11 kDa HRDC domain.

Characterization of Zn(II)·Deglycobleomycin A2 and Interaction with d(CGCTAGCG)2: Direct Evidence for Minor Groove Binding of the Bithiazole Moiety

Deglycobleomycin (dgBLM) binds to and degrades the self-complementary oligonucleotide d(CGCTAGCG)2 in a sequence-selective fashion. To characterize the binding interaction, a 1:1 complex of Zn(II)·dgBLM A2 with the DNA octanucleotide has been examined using two-dimensional NMR experiments and restrained molecular dynamics calculations. Critical elements of the mode of DNA interaction within two structural domains of Zn(II)·dgBLM A2 were fundamentally different than those observed previously for Zn(II)·BLM A2 (Manderville, R. A.; Ellena, J. F.; Hecht, S. M. J. Am. Chem. Soc. 1995, 117, 7891). A minor groove mode of binding by the bithiazole moiety is supported by the present study. Only slight upfield shifting of the bithiazole (Bit) protons Bit H5 and Bit H5‘ is observed, and the sequential intrastrand NOE connectivities are retained upon Zn(II)·dgBLM A2 binding. The orientation of the drug molecule in the complex is based on the finding of 16 intermolecular Zn(II)·dgBLM A2−DNA NOEs. The cationic C-substi...

NMR investigation of the multidrug transporter EmrE, an integral membrane protein.

EmrE is an Escherichia coli multidrug transport protein that confers resistance to a wide range of toxicants by active transport across the bacterial cell membrane. The highly hydrophobic polytopic integral membrane protein has been purified and studied in its full-length form by high-resolution NMR spectroscopy in a mixture of chloroform/methanol/water (6:6:1, by vol.). Full activity is maintained after reconstitution of the protein into proteoliposomes from this solvent mixture. A series of heteronuclear (1H-15N) two- and three-dimensional experiments, as well as triple resonance experiments, were applied to the 110-residue protein and led to the assignment of the 1H, 15N and a large part of the 13C backbone resonances as well as many of the sidechain resonances. A preliminary analysis of the secondary structure, based on sequential NOE connectivities, deviation of chemical shifts from random coil values and 3J(NH-H alpha) coupling constants supports a model where the protein forms four alpha-helices between residues 4-26 (TM1), 32-53 (TM2), 58-76 (TM3) and 85-106 (TM4). For the residues of helices TM2 and TM3 a significant line broadening occurs due to slow conformational processes.

An efficient and economic site-specific deuteration strategy for NMR studies of homologous oligonucleotide repeat sequences.

We describe herein the use of a 2H-labeling strategy to achieve specific assignments of considerably overlapped cross peaks in the 1H-NMR spectrum of a DNA trinucleotide repeat sequence. Our strategy focuses on site-specific 2H-labeling of base moieties to simplify the NMR spectral regions which contain the major portion of the structural information. To achieve efficient preparation of 2H8- or 2H6-labeled DNA and RNA nucleosides and nucleotides, the existing synthetic and purification procedures were significantly improved. Our experiments demonstrate that pyrimidine H6 deuteration reactions may be carried out using non-deuterated base reagents with DMSO-d6 as a 2H donor. These reactions are simple and economic to perform and produce base deuterated nucleosides and nucleotides in high yield. The 2H-labeled residues have been incorporated into oligonucleotides with minor modifications of the existing reaction conditions. Using the homologous CGG repeat sequence, d(CGG)5, as an example, the effectiveness of the site-specific base deuteration strategy is demonstrated. In the otherwise extensively overlapped spectra of d(CGG)5, 2H-labeling has permitted unambiguous identification of a sequential connectivity at a central CG step and confirmation of several other NOE assignments. This information is critical for elucidation of the structure and the folding of the CGG repeat sequences and will contribute to the intensive effort to understand the mechanisms of triplet expansion, which has been implicated in the development of a number of hereditary neurodegenerative diseases. In addition to the two dimensional spectral simplification in a key spectral region using site-specific 2H8/2H6-labeling, the potential applications of the prescribed strategy in homonuclear three dimensional experiments are also discussed.

Structure and dynamics of the glucocorticoid receptor DNA-binding domain: comparison of wild type and a mutant with altered specificity.

Nuclear magnetic resonance was used to compare parameters reflecting solution structure and dynamics of the glucocorticoid receptor DNA-binding domain (GRDBD), which binds specifically to a GRE binding site on DNA, and a triple mutant (GRDBDEGA), which binds to an ERE site. The studies were prompted by an earlier observation that the cooperativity for dimeric DNA-binding is 10 times higher for the GRDBDEGA-ERE association than for the GRDBD-GRE association (Lundbäck et al., 1994). The higher binding cooperativity of the mutant was unexpected since the triple mutation (G458E, S459G, and V462A) is made in the recognition helix and distant from the dimerization surface which is formed by residues in the fragment A477-N491. Sequential and long-range NOE connectivities and measured 3JHNHalpha coupling constants indicate that the overall structures of the two proteins are very similar, possibly with a less well-defined structure of the fragment K486-N491 in GRDBDEGA. However, chemical shift changes, line broadening, and increased amide proton exchange rates are observed for several residues at, or close to, the dimerization surface of the mutant. These observations are interpreted as a lower stability and/or several slowly interconverting folded conformations of this region of GRDBDEGA. The effects are likely to be due to the loss of a hydrogen bond which links S459 to the dimerization region in GRDBD. Different mechanisms for the increased binding cooperativity of the mutant are discussed, and it is noted that the properties of the GRDBDEGA dimerization region are reminiscent of those reported for the estrogen receptor DBD, which also binds to an ERE site.

CATH--a hierarchic classification of protein domain structures.

Protein evolution gives rise to families of structurally related proteins, within which sequence identities can be extremely low. As a result, structure-based classifications can be effective at identifying unanticipated relationships in known structures and in optimal cases function can also be assigned. The ever increasing number of known protein structures is too large to classify all proteins manually, therefore, automatic methods are needed for fast evaluation of protein structures. We present a semi-automatic procedure for deriving a novel hierarchical classification of protein domain structures (CATH). The four main levels of our classification are protein class (C), architecture (A), topology (T) and homologous superfamily (H). Class is the simplest level, and it essentially describes the secondary structure composition of each domain. In contrast, architecture summarises the shape revealed by the orientations of the secondary structure units, such as barrels and sandwiches. At the topology level, sequential connectivity is considered, such that members of the same architecture might have quite different topologies. When structures belonging to the same T-level have suitably high similarities combined with similar functions, the proteins are assumed to be evolutionarily related and put into the same homologous superfamily. Analysis of the structural families generated by CATH reveals the prominent features of protein structure space. We find that nearly a third of the homologous superfamilies (H-levels) belong to ten major T-levels, which we call superfolds, and furthermore that nearly two-thirds of these H-levels cluster into nine simple architectures. A database of well-characterised protein structure families, such as CATH, will facilitate the assignment of structure-function/evolution relationships to both known and newly determined protein structures.

Interactions of deglycosylated cobalt(III)-pepleomycin (green form) with DNA based on NMR structural studies,.

Pepleomycin (PEP)1 is a metalloglycopeptide antitumor antibiotic that has improved pharmacological properties than does bleomycin (BLM). Both PEP and BLM bind to and degrade DNA in a sequence-selective manner. The binding interactions of HO2--Co(III)-CodPEP (CodPEP) with CGTACG have been studied by 2D NMR and molecular modeling. Inspection of the 2D-NMR data revealed 60 notable intermolecular NOEs between CodPEP and CGTACG which place the drug's metal binding domain and peptide linker in the minor groove of the DNA close to G8 and T9. On the basis of the NOEs, the drug's DNA binding domain is located close to the T9.A4 and A10.T3 base pairs. Intercalation of the bithiazole tail between these base pairs is indicated by the loss of DNA symmetry upon complexation with CodPEP, by a break in the sequential connectivity at the TpA steps, and by the upfield shift of the bithiazole H-H5 and H-H5' proton resonances. Intercalation of the bithiazole moiety unfolds the CodPEP molecule and exposes its hydroperoxide group to the DNA. The hydroperoxide group in the refined model of CodPEP-CGTACG is close to the C4' proton of T9, consistent with cleavage at this position. The NOE pattern between the pyrimidine ring of CodPEP and G8 of DNA suggests a specific pairing recognition via hydrogen bonds between these groups, thus establishing a 5'-GT-3' sequence preference. The structural elucidation of the free CodPEP and CoPEP [Caceres-Cortes et al. (1997) Eur. J. Biochem. 244, 818-828], and of the complex of CodPEP-CGTACG afford a plausible mechanism for the recognition and its subsequent cleaving of DNA by the drug. The process involves the unfolding of the compact CodPEP, recognition of a guanine base using the metal binding domain, threading of the bithiazole tail between base pairs, and finally positioning of the HO2- group close to the T or C found 3' to the specific G site.