3D domain swapping provides a minor alternative refolding pathway for ribonuclease A

The C-terminal β-hairpin of RNase A contains a turn with a cis Asn113-Pro114 peptide bond. Pioneering pulsed HX experiments have shown that the C-terminal β-hairpin forms early during refolding. This is puzzling since the Asn113-Pro114 bond is predominately trans at this stage and this conformation destabilizes the native monomer. RNase A, when refolded at high concentration, forms a series of 3D domain-swapped oligomers. In the oligomers formed by C-terminal β-strand swapping, Asn113-Pro114 is trans and permits the formation of a new intersubunit β-sheet. We hypothesize that oligomeric species with trans Asn113-Pro114 may form during refolding. Such species could account for the HX results while comfortably accommodating Asn113-Pro114 in the trans conformation. Here, we test this hypothesis by employing chromatographic methods to detect oligomers forming in refolding conditions and find significant amounts of dimer. We propose that a 3D domain-swapped dimeric intermediate provides a minor alternative pathway for RNase A refolding.     

Proline isomerization in bovine pancreatic ribonuclease A. 2. Folding conditions

The kinetics of cis-trans isomerization of individual X-Pro peptide groups is used to study the backbone dynamics of bovine pancreatic ribonuclease A (RNase A). We previously developed and validated a fluorescence method for monitoring the cis-trans isomerization of the Tyr92-Pro93 and Asn113-Pro114 peptide groups of RNase A under unfolding conditions [Juminaga, D., Wedemeyer, W. J., and Scheraga, H. A. (1998) Biochemistry 37, 11614-11620]. The essence of this method is to introduce a fluorescent residue (Tyr or Trp) in a position adjacent to the isomerizing proline (if one is not already present) and to eliminate the fluorescence of other such residues adjacent to prolines by mutating them to phenylalanine. Here, we extend this method to observe the cis-trans isomerization of these peptide groups under folding conditions using two site-directed mutants (Y92F and Y115F) of RNase A. Both isomerizations decelerate with increasing concentrations of GdnHCl, with nearly identical m values (1.11 and 1.19 M(-1), respectively) and extrapolated zero-GdnHCl time constants (42 and 32 s, respectively); by contrast, under unfolding conditions, the cis-trans isomerizations of both Pro93 and Pro114 are independent of GdnHCl concentration. Remarkably, the isomerization rates under folding conditions at GdnHCl concentrations above 1 M are significantly slower than those measured under unfolding conditions. The temperature dependence of the Pro114 isomerization under folding conditions is also unusual; whereas Pro93 exhibits an activation energy typical of proline isomerization (19.4 kcal/mol), Pro114 exhibits a sharply reduced activation energy of 5.7 kcal/mol. A structurally plausible model accounts for these results and, in particular, shows that folding conditions strongly accelerate the cis-trans isomerization of both peptide groups to their native cis conformation, suggesting the presence of flickering local structure in their beta-hairpins.     

Solution NMR evidence for a cis Tyr-Ala peptide group in the structure of [Pro93Ala] bovine pancreatic ribonuclease A

Proline peptide group isomerization can result in kinetic barriers in protein folding. In particular, the cis proline peptide conformation at Tyr92-Pro93 of bovine pancreatic ribonuclease A (RNase A) has been proposed to be crucial for chain folding initiation. Mutation of this proline-93 to alanine results in an RNase A molecule, P93A, that exhibits unfolding/refolding kinetics consistent with a cis Tyr92-Ala93 peptide group conformation in the folded structure (Dodge RW, Scheraga HA, 1996, Biochemistry 35:1548-1559). Here, we describe the analysis of backbone proton resonance assignments for P93A together with nuclear Overhauser effect data that provide spectroscopic evidence for a type VI beta-bend conformation with a cis Tyr92-Ala93 peptide group in the folded structure. This is in contrast to the reported X-ray crystal structure of [Pro93Gly]-RNase A (Schultz LW, Hargraves SR, Klink TA, Raines RT, 1998, Protein Sci 7:1620-1625), in which Tyr92-Gly93 forms a type-II beta-bend with a trans peptide group conformation. While a glycine residue at position 93 accommodates a type-II bend (with a positive value of phi93), RNase A molecules with either proline or alanine residues at this position appear to require a cis peptide group with a type-VI beta-bend for proper folding. These results support the view that a cis Pro93 conformation is crucial for proper folding of wild-type RNase A.     

Energy difference associated with proline isomerization in ribonuclease A

We examined energy differences caused by the cis-trans transformation of every proline residues in the native structure of RNAase A. The results show that the cis form of Pro-93 and Pro-114 gave the lowest conformational energy, i.e., conformations after conversion of one of the proline residues to the trans form had a little higher energy; the transformation of Pro-42 and Pro-117 to the cis form, on the other hand, caused a much higher energy increase.     

Isolation of four novel tachykinins from frog (Rana catesbeiana) brain and intestine

In a survey for unknown bioactive peptides in frog (Rana catesbeiana) brain and intestine, we isolated four novel peptides that exhibit potent stimulant effects on smooth muscle preparation of guinea pig ileum. By microsequencing and synthesis, these peptides were identified as Lys- Pro- Ser- Pro- Asp- Arg- Phe- Tyr- Gly- Leu- Met- NH2 (ranatachykinin A), Tyr- Lys- Ser- Asp- Ser- Phe- Tyr- Gly- Leu- Met- NH2 (ranatachykinin B), His- Asn- Pro- Ala- Ser- Phe- Ile- Gly- Leu- Met- NH2 (ranatachykinin C) and Lys- Pro- Ans- Pro- Glu- Arg- Phe- Tyr- Ala- Pro- Met- NH2 (ranatachykinin D). Ranatachykinin (RTK) A, B and C conserve the C- terminal sequence, Phe- X- Gly- Leu- Met- NH2, which is common to known members of the tachykinin family. On the other hand, RTK-D has a striking feature in its C-terminal sequence, Phe- Tyr- Ala- Pro- Met- NH2, which has never been found in other known tachykinins, and may constitute a new subclass in the tachykinin family.     

A Hinge Region cis-Proline in Ribonuclease A Acts as a Conformational Gatekeeper for C-Terminal Domain Swapping

Domain swapping, the process in which a structural unit is exchanged between monomers to create a dimer containing two subunits of the monomeric fold, is believed to be an important mechanism for oligomerization and the formation of amyloid fibrils. Structural studies have implicated proline as an important residue for domain swapping due to its increased frequency in hinge regions preceding swapped arms. We hypothesized that proline's unique ability to populate both cis and trans peptide bond conformations may allow proline to act as a conformational gatekeeper, regulating interconversion between monomer and domain-swapped dimer forms. The hinge region of RNase A contains a proline at residue 114 that adopts a cis conformation in the monomer and extends to a trans conformation in the C-terminal domain-swapped dimer. Substitution of P114 with residues that strongly prefer a trans peptide bond (Ala, Gly) results in significant population of the C-terminal domain-swapped dimer under near-physiological conditions (pH 8.0, 37 °C). This is in stark contrast to dimerization of wild-type RNase A, which requires incubation under extreme conditions such as lyophilization from acetic acid or elevated temperature. In addition, we observe similar results when cis-P114 is mutated to glycine in a homologous RNase, human pancreatic RNase 1. Our results suggest that isomerization at P114 may facilitate population of a partially unfolded intermediate or alternative structure competent for domain swapping and provide support for a more general role for proline isomerization as a conformational gatekeeper in domain swapping and oligomerization.     

The aromatic residues of bovine pancreatic ribonuclease studied by 1H nuclear magnetic resonance.

1. The aromatic proton resonances in the 360-MHz 1H nuclear magnetic resonance (NMR) spectrum of bovine pancreatic ribonuclease were divided into histidine, tyrosine and phenylalanine resonances by means of pH titrations and double resonance experiments. 2. Photochemically induced dynamic nuclear polarization spectra showed that one histidine (His-119) and two tyrosines are accessibly to photo-excited flavin. This permitted the identification of the C-4 proton resonance of His-119. 3. The resonances of the ring protons of Tyr-25, Tyr-76 and Tyr-115 and the C-4 proton of His-12 were identified by comparison with subtilisin-modified and nitrated ribonucleases. Other resonances were assigned tentatively to Tyr-73, Tyr-92 and Phe-46. 4. On addition of active-site inhibitors, all phenylalanine resonances broadened or disappeared. The resonance that was most affected was assigned tentatively to Phe-120. 5. Four of the six tyrosines of bovine RNase, identified as Tyr-76, Tyr-115 and, tentatively, Tyr-73 and Tyr-92, are titratable above pH 9. The rings of Tyr-73 and Tyr-115 are rapidly rotating or flipping by 180 degrees about their C beta--C gamma bond and are accessible to flavin in photochemically induced dynamic nuclear polarization experiments. Tyr-25 is involved in a pH-dependent conformational transition, together with Asp-14 and His-48. A scheme for this transition is proposed. 6. Binding of active-site inhibitors to bovine RNase only influences the active site and its immediate surroundings. These conformational changes are probably not connected with the pH-dependent transition in the region of Asp-14, Tyr-25 and His-48. 7. In NMR spectra of RNase A at elevated temperatures, no local unfolding below the temperature of the thermal denaturation was observed. NMR spectra of thermally unfolded RNase A indicated that the deviations from a random coil are small and might be caused by interactions between neighbouring residues.     

Assignment of asparagine-44 side-chain primary amide 1H NMR resonances and the peptide amide N1H resonance of glycine-37 in basic pancreatic trypsin inhibitor

New assignments of three previously undetected amide proton NMR resonance lines in bovine pancreatic trypsin inhibitor are reported. These are the peptide amide proton of Gly-37 and the primary amide protons of Asn-44. Specific assignments of Asn-44 and Asn-43 HE and HZ resonances are also reported. The Gly-37 NH and Asn-44 HZ resonances are shifted upfield to 4.3 and 3.4 ppm, respectively, by the ring current of the Tyr-35 aromatic group, while Asn-44 HE resonates at 7.8 ppm. The abnormal chemical shifts of Asn-44 HZ and Gly-37 NH indicate that both NH's interact with the pi-electron cloud of the Tyr-35 ring. This is consistent with their location in the crystal structure. The resonances are resolved by differential labeling techniques and are studied by combined use of NOE and exchange difference spectroscopy.     

Complete assignment of the 500 MHz 1H-NMR spectra of bleomycin A2 in H2O and D2O solution by means of two-dimensional NMR spectroscopy

1H-NMR spectra of bleomycin A2 recorded at 500 MHz in D2O and H2O at 24 degrees C and 3 degrees C were investigated. Resonances of the individual spin systems were identified by using two-dimensional correlation spectroscopy (COSY), two-dimensional spin echo correlated spectroscopy (SECSY) and by the application of two-dimensional Nuclear Overhauser Effect spectroscopy (NOESY). Employment of these techniques allowed the assignment of 113 exchangeable and 59 non-exchangeable protons in the 1H NMR spectrum of bleomycin A2. By means of 2D NOE spectroscopy also interresidual connectivities could be observed. Comparison of the NOESY spectra at 3 degrees C and 24 degrees C suggest that at low temperatures the central party of the bleomycin A2 molecule tends to adopt an extended conformation.     

Role of two proline-containing turns in the folding of porcine ribonuclease

Unfolded ribonuclease (RNase) from porcine pancreas consists of a mixture of fast and slow-refolding species. The equilibrium distribution of these species differs strongly from other homologous RNases, because an additional proline residue is present at position 115 of the porcine protein. The major slow-folding species of porcine RNase contains incorrect proline isomers at Pro93 and at Pro114-Pro115. Both positions are presumably part of beta-turn structures in the native protein, as deduced from the structure of the homologous bovine RNase A. The folding kinetics of these molecules depend strongly on the conditions used. Under unfavorable conditions (near the unfolding transition), refolding is virtually blocked by the presence of the incorrect proline peptide bonds and partially folded intermediates with incorrect isomers could not be detected. As a consequence, folding is very slow under such conditions and the re-isomerization of Pro114-Pro115 is the first and rate-limiting step of folding. Under strongly native conditions (such as in the presence of ammonium sulfate), refolding is much faster. A largely folded intermediate accumulates with the turns around Pro93 and Pro114-Pro115 still in the non-native conformation. These results suggest that incorrect proline isomers strongly influence protein folding and that, under favorable conditions, the polypeptide chain can fold with two beta-turns locked into a non-native conformation. We conclude, therefore, that early formation of correct turn structure is not necessarily required for protein folding. However, the presence of incorrect turns, locked-in by non-native proline isomers, strongly decreases the rate of refolding. Alternative pathways of folding exist. The choice of pathway depends on the number and distribution of incorrect proline isomers and on the folding conditions.     

NMR structure of the alpha-hemoglobin stabilizing protein: insights into conformational heterogeneity and binding

The structure of alpha-hemoglobin stabilizing protein (AHSP), a molecular chaperone for free alpha-hemoglobin, has been determined using NMR spectroscopy. The protein native state shows conformational heterogeneity attributable to the isomerization of the peptide bond preceding a conserved proline residue. The two equally populated cis and trans forms both adopt an elongated antiparallel three alpha-helix bundle fold but display major differences in the loop between the first two helices and at the C terminus of helix 3. Proline to alanine single point mutation of the residue Pro-30 prevents the cis/trans isomerization. The structure of the P30A mutant is similar to the structure of the trans form of AHSP in the loop 1 region. Both the wild-type AHSP and the P30A mutant bind to alpha-hemoglobin, and the wild-type conformational heterogeneity is quenched upon complex formation, suggesting that just one conformation is the active form. Changes in chemical shift observed upon complex formation identify a binding interface comprising the C terminus of helix 1, the loop 1, and the N terminus of helix 2, with the exposed residues Phe-47 and Tyr-51 being attractive targets for molecular recognition. The characteristics of this interface suggest that AHSP binds at the intradimer alpha1beta1 interface in tetrameric HbA.     

Structure of the toxic domain of the Escherichia coli heat-stable enterotoxin ST I

Active fragments of the heat-stable enterotoxin ST I of Escherichia coli were chemically synthesized with the sequence Cys-Cys-Glu-Leu-Cys-Cys-Asn-Pro-Ala-Cys-Thr-Gly-Cys-(Tyr) and studied by proton (1H NMR) and carbon-13 (13C NMR) nuclear magnetic resonance spectroscopy as a function of pH and temperature. All of the nonexchangeable protons in the 1H NMR spectrum were assigned. Although all amide protons were present at temperatures below 25 degrees C and and pH values below 6, some of the resonances are broad and could not be assigned. The temperature dependence of these broad resonances indicates a change in conformation that is localized in the N-terminus. Other amide protons disappear at higher temperatures owing to chemical exchange with the solvent. Sufficient resonance assignments can be made at high and low temperatures to permit structural conclusions to be made. The chemical shifts of the alpha-carbon protons indicate the presence of substantial structure, which was further defined with the observed pattern of nuclear Overhauser enhancements (NOEs), coupling constants, and exchange rates. The NMR data identify a turn from Ala-14 to Cys-18. A second likely turn is centered around the proline residue. An interresidue NOE between the alpha-carbon protons of Asn-12 and Gly-17 indicates that the molecule folds back on itself. The NMR information is sufficient to define the structure of the C-terminal region of ST I. Manual model building then indicated that one arrangement of the three disulfides is particularly compatible with the NMR data and van der Waals constraints. A model incorporating the disulfide arrangement proposed by Houghten and his co-workers [Houghten, R.A., Ostresh, J.M., & Klipstein, F.A. (1984) Eur. J. Biochem. 145, 157-162] and the NMR constraints was derived with the programs PROTO [Frayman, F. (1985) Ph.D. Thesis, Northwestern University] and NOEMOT [Lane, A.N., Lefévre, J.-F., & Jardetsky, O. (1986) Biochim. Biophys. Acta 867, 45-56].     

Cis proline mutants of ribonuclease A. I. Thermal stability.

A chemically synthesized gene for ribonuclease A has been expressed in Escherichia coli using a T7 expression system (Studier, F.W., Rosenberg, A.H., Dunn, J.J., & Dubendorff, J.W., 1990, Methods Enzymol. 185, 60-89). The expressed protein, which contains an additional N-terminal methionine residue, has physical and catalytic properties close to those of bovine ribonuclease A. The expressed protein accumulates in inclusion bodies and has scrambled disulfide bonds; the native disulfide bonds are regenerated during purification. Site-directed mutations have been made at each of the two cis proline residues, 93 and 114, and a double mutant has been made. In contrast to results reported for replacement of trans proline residues, replacement of either cis proline is strongly destabilizing. Thermal unfolding experiments on four single mutants give delta Tm approximately equal to 10 degrees C and delta delta G0 (apparent) = 2-3 kcal/mol. The reason is that either the substituted amino acid goes in cis, and cis<==>trans isomerization after unfolding pulls the unfolding equilibrium toward the unfolded state, or else there is a conformational change, which by itself is destabilizing relative to the wild-type conformation, that allows the substituted amino acid to form a trans peptide bond.     

Involvement of prolines-114 and -117 in the slow refolding phase of ribonuclease A as determined by isomer-specific proteolysis

Using the method of isomer-specific proteolysis (ISP), the cis-trans nature of the peptide bonds involving prolines-114 and -117 in ribonuclease (RNase) has been investigated. These studies involve the pretreatment of RNase first with either a short pepsin pulse or a short mercaptoethanol pulse to irreversibly unfold the protein and then with a short chymotrypsin pulse to quickly cleave the Tyr115-Val116 bond so that the chain is suitably trimmed for the subsequent stereospecific cleavage either by aminopeptidase P, to investigate proline-117, or by a proline-specific endopeptidase, to investigate proline-114. The most reasonable interpretation of our results suggests that proline-117 is essentially 100% trans in both the native and unfolded states, so it apparently makes no direct contribution to the slow refolding kinetics of RNase. It is also determined that proline-114 is 100% cis in native RNase and ca. 95% cis in reversibly unfolded RNase so only 5% of the unfolded RNase can be rate limited by trans to cis isomerization of proline-114 during refolding. Careful spectroscopic studies of refolding show that the smallest and slowest of the refolding phases, the ct phase, has the proper amplitude (5%), relaxation time (400 s at 10 degrees C), and activation energy (17 kcal) for a phase that is rate limited by the trans to cis isomerization of proline-114. Measurements of the kinetics of binding of cytidine 2'-monophosphate during refolding further show that RNase does not become active until proline-114 has isomerized to the native cis configuration. It is concluded that none of the three prolines thus far examined (i.e., prolines-93, -114, and -117) by the ISP method is involved in the formation of a fully active, nativelike intermediate which has "incorrect" proline isomers. The specific structural process which is responsible for the largest of the three slow refolding phases, the XY phase, is still undetermined. Although ISP results on proline-42 are not yet available, it seems possible that this slow phase may be rate limited by a process other than proline isomerization. In unrelated studies, results from chymotrypsin hydrolyses of several short peptides containing the sequence -X-Y-Pro- show that cleavage of an active X-Y bond is very slow when it is immediately adjacent on the amino side of a proline peptide bond. Thus, chymotrypsin cleavage may not be generally useful as the analytical step in isomer-specific proteolysis.     

Two-dimensional NMR and structure determination of salmon calcitonin in methanol

The structure of the 32-residue peptide salmon calcitonin (sCT) in 90% MeOH-10% H2O has been investigated by two-dimensional NMR techniques and molecular modeling. Sequential assignments for nearly all of the 32 spin systems have been obtained, and results indicate that the heptaresidue loop formed by the disulfide bond between Cys-1 and Cys-7 is followed by an alpha-helical segment from Val-8 through Tyr-22. A region of conformational heterogeneity is observed for residues 20-25, resulting from the slow isomerism of the cis and trans forms of Pro-23. The C-terminal segment is found to exist in an extended conformation.     

Hydrogen bonding modulates binding of exogenous ligands in a myoglobin proximal cavity mutant.

In the sperm whale myoglobin mutant H93G, the proximal histidine is replaced by glycine, leaving a cavity in which exogenous imidazole can bind and ligate the heme iron (Barrick, D. (1994) Biochemistry 33, 6545-6554). Structural studies of this mutant suggest that serine 92 may play an important role in imidazole binding by serving as a hydrogen bond acceptor. Serine 92 is highly conserved in myoglobins, forming a well-characterized weak hydrogen bond with the proximal histidine in the native protein. We have probed the importance of this hydrogen bond through studies of the double mutants S92A/H93G and S92T/H93G incorporating exogenous imidazole and methylimidazoles. (1)H NMR spectra reveal that loss of the hydrogen bond in S92A/H93G does not affect the conformation of the bound imidazole. However, the binding constants for imidazoles to the ferrous nitrosyl complex of S92A/H93G are much weaker than in H93G. These results are discussed in terms of hydrogen bonding and steric packing within the proximal cavity. The results also highlight the importance of the trans diatomic ligand in altering the binding and sensitivity to perturbation of the ligand in the proximal cavity.     

Two-dimensional NMR studies of [Pro-10] atrial natriuretic factor [7-23

The conformational properties of an analog of atrial naturiuretic factor, [Pro-10] ANF(7-23), were examined in H₂O, H₂O/DMSO-d₆ (2/1), and DMSO-d₆ using two-dimensional nmr techniques. The sequence differs from the native peptide by the absence of the exocyclic N- and C-terminal residues, and the substitution of a proline for a glycine at position 10—a modification expected to reduce the conformational flexibility of this analog. The backbone proton nmr resonances were assigned from two-dimensional correlated spectroscopy (2D-COSY), relayed COSY, and 2D nuclear Overhàuser enhancement (NOE) experiments, and the solution conformation was evaluated from vicinal spin–spin coupling constants and NOE data. Despite the substitution of a proline in the sequence, [Pro-10] ANF(7-23) exhibits a considerable amount of flexibility in all of the solvents employed.     

Crystal structures of native and inactivated cis-3-chloroacrylic acid dehalogenase. Structural basis for substrate specificity and inactivation by (R)-oxirane-2-carboxylate

The bacterial degradation pathways for the nematocide 1,3-dichloropropene rely on hydrolytic dehalogenation reactions catalyzed by cis- and trans-3-chloroacrylic acid dehalogenases (cis-CaaD and CaaD, respectively). X-ray crystal structures of native cis-CaaD and cis-CaaD inactivated by (R)-oxirane-2-carboxylate were elucidated. They locate four known catalytic residues (Pro-1, Arg-70, Arg-73, and Glu-114) and two previously unknown, potential catalytic residues (His-28 and Tyr-103'). The Y103F and H28A mutants of these latter two residues displayed reductions in cis-CaaD activity confirming their importance in catalysis. The structure of the inactivated enzyme shows covalent modification of the Pro-1 nitrogen atom by (R)-2-hydroxypropanoate at the C3 position. The interactions in the complex implicate Arg-70 or a water molecule bound to Arg-70 as the proton donor for the epoxide ring-opening reaction and Arg-73 and His-28 as primary binding contacts for the carboxylate group. This proposed binding mode places the (R)-enantiomer, but not the (S)-enantiomer, in position to covalently modify Pro-1. The absence of His-28 (or an equivalent) in CaaD could account for the fact that CaaD is not inactivated by either enantiomer. The cis-CaaD structures support a mechanism in which Glu-114 and Tyr-103' activate a water molecule for addition to C3 of the substrate and His-28, Arg-70, and Arg-73 interact with the C1 carboxylate group to assist in substrate binding and polarization. Pro-1 provides a proton at C2. The involvement of His-28 and Tyr-103' distinguishes the cis-CaaD mechanism from the otherwise parallel CaaD mechanism. The two mechanisms probably evolved independently as the result of an early gene duplication of a common ancestor.     

Structural analysis of d(GCAATTGC)2 and its complex with berenil by nuclear magnetic resonance spectroscopy

The structures of d(GCAATTGC)2 and its complex with berenil in solution were analyzed by two-dimensional 1H NMR spectroscopy. Intra- and internucleotide nuclear Overhauser effect (NOE) connectivities demonstrate that the octanucleotide duplex is primarily in the B conformation. Binding with berenil stabilizes the duplex with respect to thermal denaturation by about 10 degrees C, based on the appearance of the imino proton signals. The berenil-d(GCAATTGC)2 system is in fast exchange on the NMR time scale. The two-dimensional NMR data reveal that berenil binds in the minor groove of d(GCAATTGC)2. The aromatic drug protons are placed within 5 A of the H2 proton of both adenines, the H1', H5', and H5" of both thymidines, and the H4', H5', and H5" of the internal guanosine. The amidine protons on berenil are also close to the H2 proton of both adenines. The duplex retains an overall B conformation in the complex with berenil. At 18 degrees C, NOE contacts at longer mixing times indicate the presence of end-to-end association both in the duplex alone and also in its complex with berenil. These intermolecular contacts either vanished or diminished substantially at 45 degrees C. Two molecular models are proposed for the berenil-(GCAATTGC)2 complex; one has hydrogen bonds between the berenil amidine protons and the carbonyl oxygen, O2, of the external thymines, and the other has hydrogen bonds between the drug amidine protons and the purine nitrogen, N3, of the internal adenines. Quantitative analysis of the NOE data favors the second model.