Three-dimensional structure of proteolytic fragment 163-231 of bacterioopsin determined from nuclear magnetic resonance data in solution.

546 NOESY cross-peak volumes were measured in the two-dimensional NOESY spectrum of proteolytic fragment 163-231 of bacterioopsin in organic solution. These data and 42 detected hydrogen bonds were applied for determining the peptide spatial structure. The fold of the polypeptide chain was determined by local structure analysis, a distance geometry approach and systematic search for energetically allowed side-chain rotamers which are consistent with experimental NOESY cross-peak volumes. The effective rotational correlation time of 6 ns for the molecule was evaluated from optimization of the local structure to meet NOE data and from the dependence on mixing time of the NiH/Ci alpha H cross-peak volumes of the residues in alpha-helical conformation. The resulting structure has two well defined alpha-helical regions, 168-191 and 198-227, with root-mean-square deviation 44 pm and 69 pm, respectively, between the backbone atoms in 14 final energy refined conformations. The alpha-helices correspond to transmembrane segments F and G of bacteriorhodopsin. The segment F contains proline 186, which introduces a kink of about 25 degrees with a disruption of the hydrogen bond with the NH group of the following residue. The segments are connected by a flexible loop region 192-197. Torsion angles chi 1 are unequivocally defined for 62% of side chains in the alpha-helices but half of them differ from electron cryo-microscopy (ECM) model of bacteriorhodopsin, apparently because of the low resolution of ECM. Nevertheless, the F and G segments can be packed as in the ECM model and with side-chain conformations consistent with all NMR data in solution.     

Calculation of conformational ensembles from potentials of mean force. An approach to the knowledge-based prediction of local structures in globular proteins

We present a prototype of a new approach to the folding problem of polypeptide chains. This approach is based on the analysis of known protein structures. It derives the energy potentials for the atomic interactions of all amino acid residue pairs as a function of the distance between the involved atoms. These potentials are then used to calculate the energies of all conformations that exist in the data base with respect to a given sequence. Then, by using only the most stable conformations, clusters of the most probable conformations for the given sequence are obtained. To discuss the results properly we introduce a new classification of segments based on their conformational stability. Special care is taken to allow for sparse data sets. The use of the method is demonstrated in the discussion of the identical oligopeptide sequences found in different conformations in unrelated proteins. VNTFV, for example, adopts a beta-strand in ribonuclease but it is found in an alpha-helical conformation in erythrocruorin. In the case of VNTFV the ensemble obtained consists of a single cluster of beta-strand conformations, indicating that this may be the preferred conformation for the pentapeptide. When the flanking residues are included in the calculation the hepapeptide P-VNTFV-H (ribonuclease) again yields an ensemble of beta-strands. However, in the ensemble of D-VNTFV-A (erythrocruorin) the major cluster is of alpha-helical type. In the present study we concentrate on the local aspects of protein conformations. However, the theory presented is quite general and not restricted to oligopeptides. We indicate extensions of the approach to the calculation of global conformations of proteins as well as conceivable applications to a number of molecular systems.     

Conformation of a trimannoside bound to mannose-binding protein by nuclear magnetic resonance and molecular dynamics simulations

A model of the carbohydrate recognition domain of the serum form of mannose-binding protein (MBP) from rat complexed with methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside is presented. Allowed conformations for the bound sugar were derived from simulated annealing protocols incorporating distance restraints computed from transferred NOESY spectra. The resulting sugar conformations were then modeled into the MBP binding site, and these models of the complex were refined using molecular dynamics (MD) simulations in the presence of solvent water. These studies indicate that only one of the two major conformations of the alpha(1-->6) linkage found in solution is significantly populated in the bound state (omega = 60 degrees ), whereas the alpha(1-->3) linkage samples at least two states, similar to its behavior in free solution. The bound conformation allows direct hydrogen bonds to form between the sugar and K182 of MBP, in addition to other water-mediated hydrogen bonds. Estimates of binding constants of candidate complexes based on changes in solvent-accessible surface areas upon binding support the NMR and MD results. These estimates further suggest that the enthalpic gains of the additional sugar-MBP interactions in a trisaccharide as opposed to a monosaccharide are offset by entropic penalties, offering an explanation for previous binding data.     

Hydrogen bonding in helical polypeptides from molecular dynamics simulations and amide hydrogen exchange analysis: alamethicin and melittin in methanol.

Molecular dynamics simulations of ion channel peptides alamethicin and melittin, solvated in methanol at 27 degrees C, were run with either regular alpha-helical starting structures (alamethicin, 1 ns; melittin 500 ps either with or without chloride counterions), or with the x-ray crystal coordinates of alamethicin as a starting structure (1 ns). The hydrogen bond patterns and stabilities were characterized by analysis of the dynamics trajectories with specified hydrogen bond angle and distance criteria, and were compared with hydrogen bond patterns and stabilities previously determined from high-resolution NMR structural analysis and amide hydrogen exchange measurements in methanol. The two alamethicin simulations rapidly converged to a persistent hydrogen bond pattern with a high level of 3(10) hydrogen bonding involving the amide NH's of residues 3, 4, 9, 15, and 18. The 3(10) hydrogen bonds stabilizing amide NH's of residues C-terminal to P2 and P14 were previously proposed to explain their high amide exchange stabilities. The absence, or low levels of 3(10) hydrogen bonds at the N-terminus or for A15 NH, respectively, in the melittin simulations, is also consistent with interpretations from amide exchange analysis. Perturbation of helical hydrogen bonding in the residues before P14 (Aib10-P14, alamethicin; T11-P14, melittin) was characterized in both peptides by variable hydrogen bond patterns that included pi and gamma hydrogen bonds. The general agreement in hydrogen bond patterns determined in the simulations and from spectroscopic analysis indicates that with suitable conditions (including solvent composition and counterions where required), local hydrogen-bonded secondary structure in helical peptides may be predicted from dynamics simulations from alpha-helical starting structures. Each peptide, particularly alamethicin, underwent some large amplitude structural fluctuations in which several hydrogen bonds were cooperatively broken. The recovery of the persistent hydrogen bonding patterns after these fluctuations demonstrates the stability of intramolecular hydrogen-bonded secondary structure in methanol (consistent with spectroscopic observations), and is promising for simulations on extended timescales to characterize the nature of the backbone fluctuations that underlie amide exchange from isolated helical polypeptides.     

The NMR structure of cyclosporin A bound to cyclophilin in aqueous solution

Cyclosporin A bound to the presumed receptor protein cyclophilin was studied in aqueous solution at pH 6.0 by nuclear magnetic resonance spectroscopy using uniform 15N- or 13C-labeling of cyclosporin A and heteronuclear spectral editing techniques. Sequence-specific assignments were obtained for all but one of the cyclosporin A proton resonances. With an input of 108 intramolecular NOEs and four vicinal 3JHN alpha coupling constants, the three-dimensional structure of cyclosporin A bound to cyclophilin was calculated with the distance geometry program DISMAN, and the structures resulting from 181 converged calculations were energy refined with the program FANTOM. A group of 120 conformers was selected on the basis of the residual constraint violations and energy criteria to represent the solution structure. The average of the pairwise root-mean-square distances calculated for the backbone atoms of the 120 structures was 0.58 A. The structure represents a novel conformation of cyclosporin A, for which the backbone conformation is significantly different from the previously reported structures in single crystals and in chloroform solution. The structure has all peptide bonds in the trans form, contains no elements of regular secondary structure and no intramolecular hydrogen bonds, and exposes nearly all polar groups to its environment. The root-mean-square distance between the backbone atoms of the crystal structure of cyclosporin A and the mean of the 120 conformers representing the NMR structure of cyclosporin A bound to cyclophilin is 2.5 A.     

The conformation of echinomycin in DMSO solution

The conformation of the antibiotic echinomycin in DMSO solution has been determined from two-dimensional NMR and distance geometry calculation with distance constraints. Five converged conformations were calculated with NOE distance constraints followed by restraint energy minimization.     

High resolution nuclear magnetic resonance studies of the conformation and orientation of melittin bound to a lipid-water interface.

Previously, the size and stoichiometry of mixed micelles of perdeuterated dodecylphosphocholine and melittin were characterized and the 1H NMR spin systems of most amino acid residues of micelle-bound melittin identified. One- and two-dimensional 1H-1H Overhauser experiments have now been used to obtain qualitative information on intramolecular proton-proton distances. These data show that the N-terminal and the C-terminal segments of melittin form two spatially distinct, compact domains; using lipid spin labels these could be located near the micelle surface. For the C-terminal domain a detailed conformation was determined by using the distance contraints from the Overhauser studies as input for a distance geometry algorithm.     

Structure of myohemerythrin in the azidomet state at 1.7/1.3 A resolution

The molecular model of myohemerythrin, an oxygen-carrying protein from sipunculan worms, has been refined by stereochemically restrained least-squares minimization at 1.7/1.3 A resolution to a conventional R-value of 0.158. The estimated positional standard deviation is better than 0.15 A for most of the 979 protein atoms. The average isotropic displacement parameter, B, for the protein atoms is 23.1 A2. This high average B parameter appears to be due to the overall motion of the molecule, which correlates with the observed anisotropic diffraction. The side-chains of seven residues were modeled in two conformations, i.e. the side-chains were discretely disordered, and B parameters for several lysine and glutamate side-chains indicate that they are poorly localized. Of the residues in myohemerythrin, 66% are helical, with 62% occurring in four long alpha-helices with mean values for the backbone torsion angles of phi = -65 degrees, psi = -42 degrees, and for the hydrogen bonds distances of N ... O, 3.0 A and H ... O, 2.1 A, and angles of N ... O = C, 153 degrees, N-H ... O, 157 degrees, and H ... O = C, 147 degrees. For two-thirds of the alpha-helical residues, the torsional rotation of the C alpha-C beta bond, chi 1, is approximately -60 degrees, and for one-third chi 1 is approximately 180 degrees. Although most turns in myohemerythrin are well-categorized by previous classification, two do not fit in established patterns. Also included in the refined model are three sulfate ions, all partially occupied, and 157 water molecules, 40% of which are modeled fully occupied. Only one water molecule is internal to the protein, the remainder occur on the surface and are observed principally between symmetry-related molecules contributing, along with van der Waals' contacts, most of the interactions between molecules. There are eight intermolecular protein-protein hydrogen bonds, of which only four are between well-located atoms.     

Closed conformation of the active site loop of rabbit muscle triosephosphate isomerase in the absence of substrate: evidence of conformational heterogeneity

The active site loop of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, alternating between the two well defined "open" and "closed" conformations. Until now the closed conformation had only been observed in protein complexes with substrate analogues. Here, we present the first rabbit muscle apo TIM structure, refined to 1.5A resolution, in which the active site loop is either in the open or in the closed conformation in different subunits of the enzyme. In the closed conformation described here, the lid loop residues participate in stabilizing hydrogen bonds characteristic of holo TIM structures, whereas chemical interactions observed in the open loop conformation are similar to those found in the apo structures of TIM. In the closed conformation, a number of water molecules are observed at the projected ligand atom positions that are hydrogen bonded to the active site residues. Additives used during crystallization (DMSO and Tris molecules and magnesium atoms) were modeled in the electron density maps. However, no specific binding of these molecules is observed at, or close to, the active site and the lid loop. To further investigate this unusual closed conformation of the apo enzyme, two more rabbit muscle TIM structures, one in the same and another in a different crystal form, were determined. These structures present the open lid conformation only, indicating that the closed conformation cannot be explained by crystal contact effects. To rationalize why the active site loop is closed in the absence of ligand in one of the subunits, extensive comparison with previously solved TIM structures was carried out, supported by the bulk of available experimental information about enzyme kinetics and reaction mechanism of TIM. The observation of both open and closed lid conformations in TIM crystals might be related to a persistent conformational heterogeneity of this protein in solution.     

Refined 1.83 A structure of trypanosomal triosephosphate isomerase crystallized in the presence of 2.4 M-ammonium sulphate. A comparison with the structure of the trypanosomal triosephosphate isomerase-glycerol-3-phosphate complex

Triosephosphate isomerase (TIM) is a dimeric glycolytic enzyme. TIM from Trypanosoma brucei brucei has been crystallized at pH 7.0 in 2.4 M-ammonium sulphate. The well-diffracting crystals have one dimer per asymmetric unit. The structure has been refined at 1.83 A resolution with an R-factor of 18.3% for all data between 6 A and 1.83 A (37,568 reflections). The model consists of 3778 protein atoms and 297 solvent atoms. Subunit 1 is involved in considerably more crystal contacts than subunit 2. Correlated with these differences in crystal packing is the observation that only in the active site of subunit 2 is a sulphate ion bound. Furthermore, significant differences with respect to structure and flexibility are observed in three loops near the active site. In particular, there is a 7 A positional difference of the tip of the flexible loop (loop 6) when comparing subunit 1 and subunit 2. Also, the neighbouring loops (loop 5 and loop 7) have significantly different conformations and flexibility. In subunit 1, loop 6 is in an "open" conformation, in subunit 2, loop 6 is in an "almost closed" conformation. Only in the presence of a phosphate-containing ligand, such as glycerol-3-phosphate, does loop 6 take up the "closed" conformation. Loop 6 and loop 7 (and also to some extent loop 5) are rather flexible in the almost closed conformation, but well defined in the open and closed conformations. The closing of loop 6 (167 to 180), as observed in the almost closed conformation, slightly changes the main-chain conformation of the catalytic glutamate, Glu167, leading to a change of the chi 1 angle of this residue from approximately -60 degrees to approximately 60 degrees and the weakening of the hydrogen bonds between its polar side-chain atoms and Ser96. In the closed conformation, in the presence of glycerol-3-phosphate, the main-chain atoms of Glu167 remain in the same position as in the almost closed conformation, but the side-chain has rotated around the CA-CB bond changing chi 1 from approximately 60 degrees to approximately -60 degrees. In this new position the hydrogen bonding to Ser96 is completely lost and also a water-mediated salt bridge between OE2(Glu167) and NE(Arg99) is lost. Comparison of the two independently refined subunits, showed that the root-mean-square deviation for all 249 CA atoms is 0.9 A; for the CA atoms of the beta-strands this is only 0.2 A. The average B-factor for all subunit 1 and subunit 2 atoms is 20 A2 and 25 A2, respectively.(ABSTRACT TRUNCATED AT 400 WORDS)     

Stability of secondary structural elements in a solvent-free environment: the α helix

The stability of the α helix as an element of secondary structure is examined in the absence of solvation, in the gas phase. Mass-analyzed ion kinetic energy (MIKE) spectrometry was applied to measure intercharge repulsion and intercharge distance in multiply protonated melittin, a polypeptide known to possess a stable helical structure in a number of different environments. The experimental results, interpreted in combination with molecular mechanics calculations, suggest that triply charged melittin retains its secondary structure in the gas phase. The stability if the α-helical conformation of the polypeptide in the absence of solvent molecules reflects the fact that a network of intrinsic helical hydrogen bonds is energetically more favorable than unfolded conformations. © 1997 Wiley-Liss, Inc.     

Conformation in solution of porcine brain natriuretic peptide determined by combined use of nuclear magnetic resonance and distance geometry

The conformation in solution of porcine brain natriuretic peptide was determined by combined use of NMR spectroscopy and distance geometry. A set of 157 inter-proton-distance constraints was derived from the two-dimensional NOE spectra, and further a set of three hydrogen bond constraints was obtained from analysis of the temperature dependence of labile protons. The five structures with minimal violations were selected after performing distance-geometry calculations starting from 40 random initial conformations. The distance-geometry structures were further refined by the use of restrained energy minimization and restrained molecular dynamics. This structure shows a compact conformation with the carboxy-terminal region, Asn21-Tyr26, folded back to the disulfide-linked loop region, Cys4-Cys20. The characteristics of the conformation determined are as follows: conformations of the three segments interposed by glycine residues, which are Arg7-Ile12, Ser14-Leu18 and Cys20-Arg25, were well defined and the segments Arg7-Ile12 and Cys20-Arg25 are rather close to each other and nearly parallel. The biological significance of these local conformations is discussed on the basis of comparisons with those of atrial natriuretic peptide reported by Kobayashi et al.     

The embedding problem for predistance matrices

A fundamental problem in molecular biology is the determination of the conformation of macromolecules from NMR data. Several successful distance geometry programs have been developed for this purpose, for example DISGEO. A particularly difficult facet of these programs is the embedding problem, that is the problem of determining those conformations whose distances between atoms are nearest those measured by the NMR techniques. The embedding problem is the distance geometry equivalent of the multiple minima problem, which arises in energy minimization approaches to conformation determination. We show that the distance geometry approach has some nice geometry not associated with other methods that allows one to prove detailed results with regard to the location of local minima. We exploit this geometry to develop some algorithms which are faster and find more minima than the algorithms presently used. The authors were partially supported by National Science Foundation Grant CHE-8802341.     

Combined use of proton-proton Overhauser enhancements and a distance geometry algorithm for determination of polypeptide conformations. Application to micelle-bound glucagon

In a new approach for the determination of polypeptide conformation, experimental data on intramolecular distances between pairs of hydrogen atoms obtained from nuclear Overhauser enhancement studies are used as input for a distance geometry algorithm. The algorithm determines the limits of the conformation space occupied by the polypeptide chain. The experimental data are used in such a way that the real conformation should in all cases be within these limits. Two important features of the method are that the results do not depend critically on the accuracy of the distance measurements by nuclear Overhauser enhancement studies and that internal mobility of the polypeptide conformation is explicitly taken into consideration. The use of this new procedure is illustrated with a structural study of the region 19-27 of glucagon bound to perdeuterated dodecylphosphocholine micelles.     

Statistical and molecular dynamics studies of buried waters in globular proteins

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.     

High-resolution study of the three-dimensional structure of horse heart metmyoglobin

The three-dimensional structure of horse heart metmyoglobin has been refined to a final R-factor of 15.5% for all observed data in the 6.0 to 1.9 A resolution range. The final model consists of 1242 non-hydrogen protein atoms, 154 water molecules and one sulfate ion. This structure has nearly ideal bonding and bond angle geometry. A Luzzati plot of the variation in R-factor with resolution yields an estimated mean co-ordinate error of 0.18 A. An extensive analysis of the pattern of hydrogen bonds formed in horse heart metmyoglobin has been completed. Over 80% of the polypeptide chain is involved in eight helical segments, of which seven are composed mainly of alpha-helical (3.6(13))-type hydrogen bonds; the remaining helix is composed entirely of 3(10) hydrogen bonds. Altogether, of 102 hydrogen bonds between main-chain atoms only six are not involved in helical structures, and four of these six occur within beta-turns. The majority of water molecules in horse heart metmyoglobin are found in solvent networks that range in size from two to 35 members. The size of water molecule networks can be rationalized on the basis of three factors: the number of hydrogen bonds to the protein surface, the presence of charged side-chain atoms, and the ability to bridge to neighboring molecules in the crystal lattice. Bridging water networks form the dominant intermolecular interactions. The backbone conformation of horse heart metmyoglobin is very similar to sperm whale metmyoglobin, with significant differences in secondary structure occurring only near residues 119 and 120, where residues 120 to 123 in sperm whale form a distorted type I reverse turn and the horse heart protein has a type II turn at residues 119 to 122. Nearly all of the hydrogen bonds between main-chain atoms (occurring mainly in helical regions) are common to both proteins, and more than half of the hydrogen bonds involving side-chain atoms observed in horse heart are also found in sperm whale metmyoglobin. Unlike sperm whale metmyoglobin, the heme iron atom in horse heart metmyoglobin is not significantly displaced from the plane of the heme group.     

Conformation of di-n-propylglycine residues (Dpg) in peptides: crystal structures of a type I' beta-turn forming tetrapeptide and an alpha-helical tetradecapeptide.

The crystal structures of two oligopeptides containing di-n-propylglycine (Dpg) residues, Boc-Gly-Dpg-Gly-Leu-OMe (1) and Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (2) are presented. Peptide 1 adopts a type I'beta-turn conformation with Dpg(2)-Gly(3) at the corner positions. The 14-residue peptide 2 crystallizes with two molecules in the asymmetric unit, both of which adopt alpha-helical conformations stabilized by 11 successive 5 --> 1 hydrogen bonds. In addition, a single 4 --> 1 hydrogen bond is also observed at the N-terminus. All five Dpg residues adopt backbone torsion angles (phi, psi) in the helical region of conformational space. Evaluation of the available structural data on Dpg peptides confirm the correlation between backbone bond angle N-C(alpha)-C' (tau) and the observed backbone phi,psi values. For tau > 106 degrees, helices are observed, while fully extended structures are characterized by tau < 106 degrees. The mean tau values for extended and folded conformations for the Dpg residue are 103.6 degrees +/- 1.7 degrees and 109.9 degrees +/- 2.6 degrees, respectively.     

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.     

Synthesis,structure, and conformation of anti-tumor agents in the solid and solution states: hydroxyl derivatives of Ftorafur

The pyrimidine antimetabolite Ftorafur [FT; 5-fluoro-1-(tetrahydro-2-furyl)uracil] has shown significant antitumor activity in several adenocarcinomas with a spectrum of activity similar to, but less toxic than, 5-fluorouracil (5-FU). It is considered as a prodrug that acts as a depot form of 5-FU, and hence the two drugs exhibit a similar spectrum of chemotherapeutic activity. Ftorafur is metabolized in animals and humans when hydroxyl groups are introduced into the tetrahydrofuran moiety. These metabolites are also thought to be as active as ftorafur but less toxic than 5-FU. Hydroxyl derivatives: 2'-hydroxyftorafur (III), 3'-hydroxyftorafur (IV) and 2',3'-dihydroxyftorafur (II) were synthesized and X-ray and NMR studies of these hydroxyl derivatives were undertaken in our laboratories to study the structural and conformational features of Ftorafur and its metabolites in the solid and solution states. X-ray crystallographic investigations were carried out with data collected on a CAD-4 diffractometer. The structures were solved and refined using the SDP crystallographic package of Enraf-Nonius on PDP 11/34 and Microvax computers. All of the compounds studied had the base in the anti conformation. The glycosidic torsion angles varied from -20 to 60 degrees. There is an inverse correlation between the glycosyl bond distances and the chi angle. Molecules with a lower chi angle have a larger bond distance and vice versa. The sugar rings show a wide variation of conformations ranging from C2'-endo through C3'-endo to C4'-exo. The crystal structures are stabilized by hydrogen bonds involving the base nitrogen atom N3 and the hydroxyl oxygen atoms of the sugar rings as donors and the keto oxygens O2 and O4 of the base and the hydroxyl oxygen atoms O2' and O3' as acceptors. The NMR studies were carried out on Brüker 400 and 600 MHz instruments. Simulated proton spectra were obtained through Laocoon, and pseudorotational parameters were solved by Pseurot. Presence of syn or anti forms was demonstrated with the use of NOE experiments. The glycosyl conformations in solution vary more widely than in the solid state. The conformations of the sugar molecules are in agreement with the values obtained in the solid state. The studies of the structure and conformation in the solid and solution states give a model for the Ftorafur molecule that could be used in structure, function and biological activity correlation studies.     

Influence of N-terminal residue stereochemistry on the prolyl amide geometry and the conformation of 5-tert-butylproline type VI beta-turn mimics.

The effects of N-terminal amino acid stereochemistry on prolyl amide geometry and peptide turn conformation were investigated by coupling both L- and D-amino acids to (2S, 5R)-5-tert-butylproline and L-proline to generate, respectively, N-(acetyl)dipeptide N'-methylamides 1 and 2. Prolyl amide cis- and trans-isomers were, respectively, favored for peptides 1 and 2 as observed by proton NMR spectroscopy in water, DMSO and chloroform. The influence of solvent composition on amide proton chemical shift indicated an intramolecular hydrogen bond between the N'-methylamide proton and the acetamide carbonyl for the major conformer of dipeptides (S)-1, that became less favorable in (R)-1 and 2. The coupling constant (3J(NH,alpha)) values for the cis-isomer of (R)-1 indicated a phi2 dihedral angle value characteristic of a type VIb beta-turn conformation in solution. X-ray crystallographic analysis of N-acetyl-D-leucyl-5-tert-butylproline N'-methylamide (R)-lb showed the prolyl residue in a type VIb beta-turn geometry possessing an amide cis-isomer and psi3-dihedral angle having a value of 157 degrees, which precluded an intramolecular hydrogen bond. Intermolecular hydrogen bonding between the leucyl residues of two turn structures within the unit cell positioned the N-terminal residue in a geometry where their phi2 and psi2 dihedral angle values were not characteristic of an ideal type VIb turn. The circular dichroism spectra of tert-butylprolyl peptides (S)- and (R)-1b were found not to be influenced by changes in solvent composition from water to acetonitrile. The type B spectrum exhibited by (S)-1b has been previously assigned to a type VIa beta-turn conformation [Halab L, Lubell WD. J. Org. Chem. 1999; 64: 3312-3321]. The type C spectrum exhibited by the (R)-lb has previously been associated with type II' beta-turn and alpha-helical conformations in solution and appears now to be also characteristic for a type VIb geometry.