The NMR solution structure of intestinal fatty acid-binding protein complexed with palmitate: application of a novel distance geometry algorithm

The three-dimensional solution structure of rat intestinal fatty acid-binding protein (I-FABP) complexed with palmitate has been determined using multidimensional triple-resonance NMR methods. The structure is based on 3889 conformational restraints derived mostly from 3-D 13C- and 15N-resolved nuclear Overhauser (NOESY) experiments. The 3-D NOESY data for this 15.4 kDa complex contained an average of nine possible interpretations per cross-peak. To circumvent this ambiguity, an eight-stage iterative procedure was employed to gradually interpret and introduce unambiguous distance restraints during subsequent rounds of structure calculations. The first stage of this procedure relied critically upon an initial structural model based on the consensus 1H/13C chemical shift-derived secondary structure and a set of symmetry-checked restraints derived from the 3-D 13C-resolved NOESY spectrum. The structures were calculated using DISTGEOM, a program that implements a novel distance geometry algorithm with pairwise Gaussian metrization. A central feature of this algorithm is the use of an iteratively optimized Gaussian distribution for the selection of trial distances, which overcomes the tendency of metrization to produce crushed structures. In addition, this algorithm randomly selects pairwise elements of the distance matrix, which results in an improved sampling of conformational space for a given computational effort. The final family of 20 distance geometry/simulated annealing structures exhibited an average pairwise C(alpha) root-mean-square deviation of 0.98 A, and their stereochemical quality, as assessed by PROCHECK, was comparable to that of 2.5 A X-ray crystal structures. The NMR structure was compared with the X-ray crystal structure of the same ligand/protein complex and was found to be essentially identical within the precision of the results. The NMR structure was also compared with that of the palmitate complex with bovine heart FABP, which shares 30% sequence identity with rat I-FABP. The overall folds were the same, but differences were noted with respect to the presence or absence of apparent conformational heterogeneity and the location and conformation of the bound fatty acid.     

Structure of hen lysozyme in solution

The structure of the 129-residue protein hen lysozyme has been determined in solution by two-dimensional 1H nuclear magnetic resonance methods. 1158 NOE distance restraints, and 68 phi and 24 chi 1 dihedral angle restraints were employed in conjunction with distance geometry and simulated annealing procedures. The overall C alpha root-mean-square deviation from the average for 16 calculated structures is 1.8(+/- 0.2) A, but excluding 14 residues in exposed disordered regions, this value reduces to 1.3(+/- 0.2) A. Regions of secondary structure, and the four alpha-helices in particular, are well defined (C alpha root-mean-square deviation 0.8(+/- 0.3) A for helices). The main-chain fold is closely similar to structures of the protein in the crystalline state. Furthermore, many of the internal side-chains are found in well-defined conformational states in the solution structures, and these correspond well with the conformational states found in the crystal. The general high level of definition of mainchain and many internal side-chains in the solution structures is reinforced by the results of an analysis of coupling constants and ring current shifts. Many side-chains on the surface, however, are highly disordered amongst the set of solution structures. In certain cases this disorder has been shown to be dynamic in origin by the examination of 3J alpha beta coupling constants.     

Alternative measures of pesticide use

tRNA (m5U54)-methyltransferase (RUMT) catalyzes the S-adenosylmethionine-dependent methylation of uridine-54 in the T psi C-loop of all transfer RNAs in E. coli to form the 54-ribosylthymine residue. However, in all tRNA structures, residue 54 is completely buried and the question arises as to how RUMT gains access to the methylation site. A 17-mer RNA hairpin consisting of nucleotides 49-65 of the T psi-loop is a substrate for RUMT. Homonuclear NMR methods in conjunction with restrained molecular dynamics (MD) methods were used to determine the solution structure of the 17-mer T-arm fragment. The loop of the hairpin exhibits enhanced flexibility which renders the conventional NMR average structure less useful compared to the more commonly found situation where a molecule exists in predominantly one major conformation. However, when resorting to softer refinement methods such as MD with time-averaged restraints, the conflicting restraints in the loop can be satisfied much better. The dynamic structure of the T-arm is represented as an ensemble of 10 time-clusters. In all of these, U54 is completely exposed. The flexibility of the T psi-loop in solution in conjunction with extensive binding studies of RUMT with the T psi C-loop and tRNA suggest that the specificity of the RUMT/ tRNA recognition is associated with tRNA tertiary structure elements. For the methylation, RUMT would simply have to break the tertiary interactions between the D- and T-loops, leading to a melting of the T-arm structure and making U54 available for methylation.     

Solution structure of recombinant human interleukin-6

Interleukin-6 (IL-6) is a 185 amino acid cytokine which exerts multiple biological effects in vivo and whose dysregulation underlies several disease processes. The solution structure of recombinant human interleukin-6 has now been determined using heteronuclear three and four-dimensional NMR spectroscopy. The structure of the molecule was determined using 3044 distance and torsion restraints derived by NMR spectroscopy to generate an ensemble of 32 structures using a combined distance geometry/simulated annealing protocol. The protein contains five alpha-helices interspersed with variable-length loops; four of these helices constitute a classical four-helix bundle with the fifth helix located in the CD loop. There were no distance violations greater than 0.3 A in any of the final 32 structures and the ensemble has an average-to-the-mean backbone root-mean-square deviation of 0.50 A for the core four-helix bundle. Although the amino-terminal 19 amino acids are disordered in solution, the remainder of the molecule has a well defined structure that shares many features displayed by other long-chain four-helix bundle cytokines. The high-resolution NMR structure of hIL-6 is used to rationalize available mutagenesis data in terms of a heteromeric receptor complex.     

Solution structure of the superactive monomeric des-[Phe(B25)] human insulin mutant: elucidation of the structural basis for the monomerization of des-[Phe(B25)] insulin and the dimerization of native insulin

The three-dimensional solution structure of des-[Phe(B25)] human insulin has been determined by nuclear magnetic resonance spectroscopy and restrained molecular dynamics calculations. Thirty-five structures were calculated by distance geometry from 581 nuclear Overhauser enhancement-derived distance constraints, ten phi torsional angle restraints, the restraints from 16 helical hydrogen bonds, and three disulfide bridges. The distance geometry structures were optimized using simulated annealing and restrained energy minimization. The average root-mean-square (r.m.s.) deviation for the best 20 refined structures is 1.07 angstroms for the backbone and 1.92 angstroms for all atoms if the less well-defined N and C-terminal residues are excluded. The helical regions are more well defined, with r.m.s. deviations of 0.64 angstroms for the backbone and 1.51 angstroms for all atoms. It is found that the des-[Phe(B25)] insulin is a monomer under the applied conditions (4.6 to 4.7 mM, pH 3.0, 310 K), that the overall secondary and tertiary structures of the monomers in the 2Zn crystal hexamer of native insulin are preserved, and that the conformation-averaged NMR solution structure is close to the structure of molecule 1 in the hexamer. The structure reveals that the lost ability of des-[Phe(B25)] insulin to self-associate is caused by a conformational change of the C-terminal region of the B-chain, which results in an intra-molecular hydrophobic interaction between Pro(B28) and the hydrophobic region Leu(B11)-Leu(B15) of the B-chain alpha-helix. This interaction interferes with the inter-molecular hydrophobic interactions responsible for the dimerization of native insulin, depriving the mutant of the ability to dimerize. Further, the structure displays a series of features that may explain the high potency of the mutant on the basis of the current model for the insulin-receptor interaction. These features are: a change in conformation of the C-terminal region of the B-chain, the absence of strong hydrogen bonds between this region and the rest of the molecule, and a relatively easy accessibility to the Val(A3) residue.     

NMR study suggests a major role for Arg111 in maintaining the structure and dynamical properties of type II human cellular retinoic acid binding protein

The solution structure of a site-directed mutant of type-II human cellular retinoic acid binding protein (CRABPII) with Arg111 replaced by methionine (R111M) has been determined by NMR spectroscopy. The sequential assignments of the 1H and 15N resonances of apo-R111M were established by multinuclear multidimensional NMR. The solution structure was calculated from 2302 distance restraints and 77 phi dihedral restraints derived from the NMR data. The root-mean-square deviation of the ensemble of 28 refined conformers that represent the structure from the mean coordinate set derived from them was 0.54 +/- 0.26 and 0.98 +/- 0.23 A for the backbone atoms and all heavy atoms, respectively. The solution structure of apo-R111M is similar to that of wild-type apo-CRABPII. However, there are significant conformational differences between the two proteins, localized mainly to three segments (Leu19-Ala36, Glu73-Cys81, and Leu99-Pro105) clustered around the ligand entrance more than 17 A away from the point mutation. In apo-R111M, all the three segments move toward the center of the ligand entrance so that the opening of the ligand-binding pocket in apo-R111M is much smaller than that in wild-type apo-CRABPII. Furthermore, the ligand-binding pocket of apo-R111M, especially the ligand entrance, is much less flexible than that of apo-CRABPII. Surprisingly, apo-R111M is more similar to holo-CRABPII than to apo-CRABPII in both structure and dynamical properties. The conformational and dynamical changes caused by the mutation are similar to those induced by binding of RA, although the magnitudes of the changes caused by the mutation are smaller than those induced by binding of RA. The results suggest that Arg111 plays a critical role in determining the structure and dynamical properties of CRABPII.     

Solution structure and dynamics of linked cell attachment modules of mouse fibronectin containing the RGD and synergy regions: comparison with the human fibronectin crystal structure

We report the three-dimensional solution structure of the mouse fibronectin cell attachment domain consisting of the linked ninth and tenth type III modules, mFnFn3(9,10). Because the tenth module contains the RGD cell attachment sequence while the ninth contains the synergy region, mFnFn3(9,10) has the cell attachment activity of intact fibronectin. Essentially complete signal assignments and approximately 1800 distance and angle restraints were derived from multidimensional heteronuclear NMR spectra. These restraints were used with a hybrid distance geometry/simulated annealing protocol to generate an ensemble of 20 NMR structures having no distance or angle violations greater than 0.3 A or 3 degrees. Although the beta-sheet core domains of the individual modules are well-ordered structures, having backbone atom rmsd values from the mean structure of 0.51(+/-0.12) and 0.40(+/-0.07) A, respectively, the rmsd of the core atom coordinates increases to 3.63(+/-1.41) A when the core domains of both modules are used to align the coordinates. The latter result is a consequence of the fact that the relative orientation of the two modules is not highly constrained by the NMR restraints. Hence, while structures of the beta-sheet core domains of the NMR structures are very similar to the core domains of the crystal structure of hFnFn3(9,10), the ensemble of NMR structures suggests that the two modules form a less extended and more flexible structure than the fully extended rod-like crystal structure. The radius of gyration, Rg, of mFnFn3(9,10) derived from small-angle neutron scattering measurements, 20.5(+/-0.5) A, agrees with the average Rg calculated for the NMR structures, 20.4 A, and is ca 1 A less than the value of Rg calculated for the X-ray structure. The values of the rotational anisotropy, D ||/D perpendicular, derived from an analysis of 15N relaxation data, range from 1.7 to 2.1, and are significantly less than the anisotropy of 2.67 predicted by hydrodynamic modeling of the crystal coordinates. In contrast, hydrodynamic modeling of the NMR coordinates yields anisotropies in the range of 1.9 to 2.7 (average 2.4(+/-0.2)), with NMR structures bent by more than 20 degrees relative the crystal structure having calculated anisotropies in best agreement with experiment. In addition, the relaxation parameters indicate that several loops in mFnFn3(9,10), including the RGD loop, are flexible on the nanosecond to picosecond time-scale. Taken together, our results suggest that, in solution, the limited set of interactions between the mFnFn3(9,10) modules position the RGD and synergy regions to interact specifically with cell surface integrins, and at the same time permit sufficient flexibility that allows mFnFn3(9,10) to adjust for some variation in integrin structure or environment.     

The three-dimensional NMR solution structure of alpha-cobratoxin at pH 7.5 and conformational differences with the NMR solution structure at pH 3.2

The 3D solution structure of alpha-cobratoxin, a neurotoxin purified from the Naja naja siamensis snake venom, has been determined by Nuclear Magnetic Resonance spectroscopy, in conjunction with distance geometry and restrained molecular dynamics, at pH 7.5. A total of 490 distance restraints were obtained from NOE intensities and 25 phi dihedral angle restraints deduced from J-coupling data. The generated structures are well defined with root mean square deviations from a geometrical mean structure of 0.107 +/- 0.036 nm for the backbone atoms and 0.128 +/- 0.073 nm for the side-chain atoms (considering residues 1 to 66 minus 26 to 35). A comparison between the generated structures at pH 7.5 and the mean NMR solution structure at pH 3.2 revealed that the 3D structure of alpha-cobratoxin is more compact at neutral pH. This major difference is mainly due to the pH-dependent conformational variations of three residues His18, Thr44 and Thr59.     

Solution structure and internal motion of a bioactive peptide derived from nerve growth factor

The conformation and internal dynamics of a bioactive cyclic peptide, N-acetyl-YCTDEKQCY, derived from the C-D loop of beta-nerve growth factor (beta-NGF) were analyzed by solution NMR spectroscopy. NMR experimental data were used to calculate an ensemble of peptide structures. All of the structures had a beta-turn at residues Asp4-Gln7 but could be divided into two families according the presence or absence of a hydrogen bond at Gln7. Comparison of the calculated structures with the corresponding C-D loops from the x-ray structures of the NGF revealed striking similarity. The orientation of Glu5, Lys6, and Gln7 side chains in the NGF mimetic was very similar to the C-D loop of NGF. These residues are known to participate in interactions with the TrkA receptor. Relaxation measurements of the peptidomimetic alpha-carbons at 13C natural abundance and calculated dynamic parameters suggest that the loop region of peptide is well structured but that residues Thr3, Asp4, Glu5, and Lys6 undergo slow conformational exchange. These results suggest that conformational similarity and possibly peptide dynamics are responsible for the bioactivity of the peptide.     

Structure and stability of an immunoglobulin superfamily domain from twitchin, a muscle protein of the nematode Caenorhabditis elegans

The NMR solution structure of an immunoglobulin superfamily module of twitchin (Ig 18') has been determined and the kinetic and equilibrium folding behaviour characterised. Thirty molecular coordinates were calculated using a hybrid distance geometry-simulated annealing protocol based on 1207 distance and 48 dihedral restraints. The atomic rms distributions about the mean coordinate for the ensemble of structures is 0.55( +/- 0.09) A for backbone atoms and 1.10( +/- 0.08) A for all heavy atoms. The protein has a topology very similar to that of telokin and the titin Ig domains and thus it falls into the I set of the immunoglobulin superfamily. The close agreement between the predicted and observed structures of Ig 18' demonstrates clearly that the I set profile can be applied in the structure prediction of immunoglobulin-like domains of diverse modular proteins. Folding studies reveal that the protein has relatively low thermodynamic stability, deltaG(H2O)U-F = 4.0 kcal mol(-1) at physiological pH. Unfolding studies suggest that the protein has considerable kinetic stability, the half life of the unfolding is greater than 40 minutes in the absence of denaturant.     

Solution structure of Lqh-8/6, a toxin-like peptide from a scorpion venom--structural heterogeneity induced by proline cis/trans isomerization

Lqh-8/6 is a minor fraction isolated from the venom of the scorpion Leiurus quinquestriatus hebraeus. Here we describe the purification, amino acid sequencing and solution structure determination by NMR and molecular modeling of this peptide. Lqh-8/6 is a small polypeptide (38 residues) which contains 8 half-cystines and is highly similar to another venom component, chlorotoxin. Standard homonuclear methods were used to sequentially assign the proton NMR spectra and to collect spatial restraints for structure determination. Two populations, identified early in the assignment step, are in slow interconversion on the NMR timescale. The two conformers were shown to originate from a cis/trans peptidyl-prolyl isomerization. Using a distance geometry program and simulated annealing protocol under the NMR restraints we obtained 10 final structures for the major conformation (trans isomer). None of the structures showed NOE violations larger than 0.05 nm, and the rmsd value relative to the mean structure (considering the main chain atoms in well-defined secondary structure) is 0.07 nm. The three-dimensional structure contains a short alpha-helix strapped on a small antiparallel beta-strand and an N-terminal extended fragment. The sequence/structure and structure/function relationships of the new scorpion toxin-like peptide are discussed in the context of the present structure determination. This toxin shows a stable, highly populated cis conformer of a peptidyl-prolyl peptide bond.     

Investigation of the solution structure of chymotrypsin inhibitor 2 using molecular dynamics: comparison to x-ray crystallographic and NMR data

The native solution structure and dynamics of chymotrypsin inhibitor 2 (CI2) have been studied using a long (5.3 ns) molecular dynamics (MD) simulation without any imposed restraints. The majority of the experimentally observed spin-spin coupling constants, short- and long-range nuclear Overhauser effect (NOE) cross peaks and the amide hydrogen exchange behavior were reproduced by the MD simulation. This good correspondence suggests that the major structural features of the protein during the simulation are representative of the true protein structure in solution. Two water molecules formed hydrogen bond bridges between beta2 and beta3, in agreement with X-ray crystallographic data and a recent reassessment of the solution structure using time-averaged NMR restraints during MD refinement. The active-site loop of the protein displayed the greatest structural changes and the highest mobility. When this loop region was excluded, the average Calpha r.m.s. deviation of the simulated solution structures from the crystal structure was approximately 1.5 Angstrom from 0.5 to 5.3 ns. There is structural heterogeneity in particular regions of the NMR-derived solution structures, which could be a result of imprecision or true internal motion. A study of the distribution of mobility through the protein allows us to distinguish between these two alternatives. In particular, deviations in the active-site loop appear to be a result of heightened mobility, which is also supported by good correspondence between calculated and experimental S2 N-H order parameters. On the other hand, other ill-defined regions of the NMR-derived structures are well defined in the simulation and are probably the result of a lack of structural restraints (i.e. NOEs), as opposed to reflecting the true mobility.     

Characterization of the three-dimensional solution structure of human profilin: 1H, 13C, and 15N NMR assignments and global folding pattern

Human profilin is a 15-kDa protein that plays a major role in the signaling pathway leading to cytoskeletal rearrangement. Essentially complete assignment of the 1H, 13C, and 15N resonances of human profilin have been made by analysis of multidimensional, double- and triple-resonance nuclear magnetic resonance (NMR) experiments. The deviation of the 13C alpha and 13C beta chemical shifts from their respective random coil values were analyzed and correlate well with the secondary structure determined from the NMR data. Twenty structures of human profilin were refined in the program X-PLOR using a total of 1186 experimentally derived conformational restraints. The structures converged to a root mean squared distance deviation of 1.5 A for the backbone atoms. The resultant conformational ensemble indicates that human profilin is an alpha/beta protein comprised of a seven-stranded, antiparallel beta-sheet and three helices. The secondary structure elements for human profilin are quite similar to those found in Acanthamoeba profilin I [Archer, S. J., Vinson, V. K., Pollard, T. D., & Torchia, D. A. (1993), Biochemistry 32, 6680-6687], suggesting that the three-dimensional structure of Acanthamoeba profilin I should be analogous to that determined here for human profilin. The structure determination of human profilin has facilitated the sequence alignment of lower eukaryotic and human profilins and provides a framework upon which the various functionalities of profilin can be explored. At least one element of the actin-binding region of human profilin is an alpha-helix. Two mechanisms by which phosphatidylinositol 4,5-bisphosphate can interfere with actin-binding by human profilin are proposed.     

Intrinsic nature of the three-dimensional structure of proteins as determined by distance geometry with good sampling properties

A protocol for distance geometry calculation is shown to have excellent sampling properties in the determination of three-dimensional structures of proteins from nuclear magnetic resonance (NMR) data. This protocol uses a simulated annealing optimization employing mass-weighted molecular dynamics in four-dimensional space (Havel, T.F. (1991) Prog. Biophys. Mol. Biol., 56, 43-78). It attains an extremely large radius of convergence, allowing a random coil conformation to be used as the initial estimate for the succeeding optimization process. Computations are performed with four systems of simulated distance data as tests of the protocol, using an unconstrained L-alanine 30mer and three different types of proteins, bovine pancreatic trypsin inhibitor, the alpha-amylase inhibitor Tendamistat, and the N-terminal domain of the 434-repressor. The test of the unconstrained polypeptide confirms that the sampled conformational space is that of the statistical random coil. In the larger and more complicated systems of the three proteins, the protocol gives complete convergence of the optimization without any trace of initial structure dependence. As a result of an exhaustive conformational sampling by the protocol, the intrinsic nature of the structures generated with distance restraints derived from NMR data has been revealed. When the sampled structures are compared with the corresponding X-ray structures, we find that the averages of the sampled structures always show a certain pattern of discrepancy from the X-ray structure. This discrepancy is due to the short distance nature of the distance restraints, and correlates with the characteristic shape of the protein molecule.     

Comparative conformational analyses of mu-selective dermorphin and delta-selective deltorphin-II in aqueous solution by 1H-NMR spectroscopy

Two-dimensional 1H-NMR methods have been used to obtain complete proton resonance assignments and possible solution conformations of dermorphin (H-Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH2) and deltorphin-II (H-Tyr-D-Ala-Phe-Glu-Val-Val-Gly-NH2), naturally occurring mu- and delta-selective opioids, respectively, in order to examine the conformational characteristics that are closely related to the selectivities towards mu/delta-opioid receptors. With the use of the proton-proton distances derived from ROESY measurements in aqueous solution, 50 possible 3D structures are generated by means of distance geometry calculations. The conformers which satisfy the distance constraints and the torsion angles estimated from JNHC alpha H vicinal coupling constants within the allowable range are then subjected to molecular dynamics simulations for 10 ps after equilibration. Although dermorphin and deltorphin-II are both in equilibrium among many flexible conformers, some conformational differences are observed between these peptides: many conformers of dermorphin show a structure rounded at the N-terminal Tyr-D-Ala-Phe-Gly-Tyr and C-terminal Gly-Tyr-Pro-Ser-NH2 moieties, which are almost at right angles to each other, while those of deltorphin-II are characterized by a 'hook'-shaped backbone structure in which the nearly extended conformation of the Val-Val-Gly-NH2 sequence is located under the folded conformation of the N-terminal Tyr-D-Ala-Phe-Glu sequence. The possible relationship between these conformational characteristics and the mu/delta-opioid receptor selectivities is discussed.     

Solution structure of nodularin. An inhibitor of serine/threonine-specific protein phosphatases

The three-dimensional solution structure of nodularin was studied by NMR and molecular dynamics simulations. The conformation in water was determined from the distance and dihedral data by distance geometry and refined by iterative relaxation matrix analysis. The cyclic backbone adopts a well defined conformation but the remote parts of the side chains of arginine as well as the amino acid derivative Adda have a large spatial dispersion. For the unusual amino acids the partial charges were calculated and nodularin was subjected to molecular dynamic simulations in water. A good agreement was found between experimental and computational data with hydrogen bonds, solvent accessibility, molecular motion, and conformational exchange. The three-dimensional structure resembles very closely that of microcystin-LR in the chemically equivalent segment. Therefore, it is expected that the binding of both microcystins and nodularins to serine/threonine-specific protein phosphatases is similar on an atomic level.     

Dynamical studies of peptide motifs in the Plasmodium falciparum circumsporozoite surface protein by restrained and unrestrained MD simulations

The immunodominant region on the circumsporozoite surface (CS) protein of the malaria parasite Plasmodium falciparum contains 37 repeated copies of a asparagine-alanine-asparagine-proline (NANP) motif NMR studies of linear synthetic peptides containing one, two or three repeat units provided evidence for nascent type I beta-turns within the NPNA cadence in aqueous solution. The beta-turns could be stabilised upon substituting proline for alpha-methylproline (p(Me)) in the dodecamer (NP(Me)NA)3, without loss of the ability to elicit antibodies cross-reactive with P. falciparum sporozoites. In this work, four 4 ns MD simulations of the dodecapeptide Acetyl-(NP(Me)NA)3, in water, using NOE distance restraints, using 3J-coupling constant restraints, using both these restraints and without restraints, were carried out to determine the conformations of this peptide in aqueous solution. An unrestrained MD simulation of the unmethylated Ac-(NPNA)3 peptide in water was also carried out to investigate the effect of the additional methyl groups on the structure and dynamics of the peptide. The application of NOE distance restraints and 3J-coupling constant restraints leads to contradictory results, probably due to different averaging time scales inherent to the measurement of these data, which exceed the 100 ps averaging applied in the simulations. The additional methyl groups lead to more compact structures, which display enhanced local fluctuations. The central tetrapeptide adopts a type I beta-turn, while the outer motifs display more conformational variability. The three motifs in the methylated dodecamer peptide, however, adopt frequently in the distance restrained MD simulation a compact structure such that the outer motifs appear to form a hydrophobic core by stacking of their two proline rings. This arrangement also suggests how a peptide containing multiple tandemly linked copies of a stable beta-turn NPNA motif might adopt a folded stem-like structure, which conceivably may be of biological relevance in the native CS protein.     

Structure of a DNA duplex containing a single 2'-O-methyl-beta-D-araT: combined use of NMR, restrained molecular dynamics, and full relaxation matrix refinement

Two-dimensional 1H NMR spectroscopy was used to determine the solution structure of the double-stranded DNA oligonucleotide d(5'-CGCATATAGCC-3'): d(5'-GGCTAXATGCG-3'), where X is 1-(2-O-methyl-beta-D-arabinofuranosyl)thymine. The structure determination was based on a total relaxation matrix analysis of NOESY cross-peak intensities using the MARDIGRAS program. The improved RANDMARDI procedure was used during the calculations to include the experimental "noise" in the NOESY spectra. The NOE-derived distance restraints were applied in restrained molecular dynamics calculations. Twenty final structures each were generated for the modified DNA duplex from both A-form and B-form DNA starting structures. The root-mean-square deviation of the coordinates for the 40 structures was 0.82 A. The duplex adopts a normal B-DNA-type helix, and the spectra as well as the structure show that the modified nucleotide X adopts a C2'-endo (S) sugar conformation. There are no significant changes in the helix originating from the modified nucleotide. The CH3O group on X is directed toward the major groove, and there seems to be free space for further modifications at this position.     

Internal motions and hydration of sucrose in a diluted water solution

This paper represents a synthesis of our most recent work on the hydration, internal and overall dynamics of sucrose in a diluted water solution. The studies were carried out as a total ensemble of 1.2 nanosecond condensed phase molecular dynamics trajectories. In this study the focus is on a 500 ps trajectory starting with the solute in the crystalline conformation. The presence of water was found to significantly alter the accessible conformational space of the solute. All potential intra-molecular hydrogen bonds were found to be exchanged to surrounding water molecules and the simulations suggest that the sucrose conformation is stabilized by the dynamic presence of two interring bridging water molecules: O-2g...Ow...O-3f and O-2g...Ow...O-lf. The overall shift in conformation of the solute induced by the presence of water was found to improve the theoretical models of experimental traits. It is demonstrated that the hydration structure and the internal and overall motions of sucrose compare extremely well with NMR data such as glycosidic heteronuclear coupling constants and the molecular tumbling time, with X-ray data of two partially hydrated sucrose structures in a protein complex and with translational diffusion coefficients and hydration numbers established from experimental studies.     

The solution structure of the complex of Lactobacillus casei dihydrofolate reductase with methotrexate

We have determined the three-dimensional solution structure of the complex of Lactobacillus casei dihydrofolate reductase (18.3 kDa, 162 amino acid residues) formed with the anticancer drug methotrexate using 2531 distance, 361 dihedral angle and 48 hydrogen bond restraints obtained from analysis of multidimensional NMR spectra. Simulated annealing calculations produced a family of 21 structures fully consistent with the constraints. The structure has four alpha-helices and eight beta-strands with two other regions, comprising residues 11 to 14 and 126 to 127, also interacting with each other in a beta-sheet manner. The methotrexate binding site is very well defined and the structure around its glutamate moiety was improved by including restraints reflecting the previously determined specific interactions between the glutamate alpha-carboxylate group with Arg57 and the gamma-carboxylate group with His28. The overall fold of the binary complex in solution is very similar to that observed in the X-ray studies of the ternary complex of L. casei dihydrofolate reductase formed with methotrexate and NADPH (the structures of the binary and ternary complexes have a root-mean-square difference over the backbone atoms of 0.97 A). Thus no major conformational change takes place when NADPH binds to the binary complex. In the binary complex, the loop comprising residues 9 to 23 which forms part of the active site has been shown to be in the "closed" conformation as defined by M. R. Sawaya & J. Kraut, who considered the corresponding loops in crystal structures of complexes of dihydrofolate reductases from several organisms. Thus the absence of the NADPH does not result in the "occluded" form of the loop as seen in crystal studies of some other dihydrofolate reductases in the absence of coenzyme. Some regions of the structure in the binary complex which form interaction sites for NADPH are less well defined than other regions. However, in general terms, the NADPH binding site appears to be essentially pre-formed in the binary complex. This may contribute to the tighter binding of coenzyme in the presence of methotrexate.