The conformation of deamino-oxytocin: X-ray analysis of the 'dry' and 'wet' forms

Two crystal structures of (1 beta-mercaptopropionic acid) deamino-oxytocin are reported. The 'dry form' in space group C2 has cell dimensions a = 27.08 +/- 0.03, b = 9.06 +/- 0.01, c = 22.98 +/- 0.02 A, beta = 102.06 +/- 0.03 with one deamino-oxytocin and six water molecules per asymmetric unit. The 'wet form' in space group P2(1) has cell dimensions a = 27.27 +/- 0.02, b = 9.04 +/- 0.01, c = 23.04 +/- 0.02 A, beta = 102.24 +/- 0.02, with two deamino-oxytocin and 13 water molecules per asymmetric unit. A local twofold parallel to the monoclinic axis gives a pseudo C2 packing. Initial phases of the 'dry form' were calculated by the heavy-atom method from the isomorphous and anomalous difference Pattersons and anomalous difference Fouier synthesis. The structure was refined by using restrained least-squares at 1.2 A resolution to a crystallographic R = 0.10. The molecular replacement method yielded the P2(1) structure that was refined with geometric restraints to R less than 0.09, by using all data to 1.09 A resolution. Deamino-oxytocin consists of a cyclic tocin ring formed by six amino acids, closed by a disulphide bridge, S1-S6, and held by two trans-annular hydrogen bonds N2-O5 and N5-O2 with a type II turn at residues 3 and 4. A flexible tripeptide tail has a loosely hydrogen-bonded type I beta-turn between N9 and O6. The sulphur of cysteine at position 1 is disordered in all the molecules leading to alternative hands of disulphide. The conformational flexibility of Ile 3, Asn 5, Pro 7 side chains and the disulphide bridge is consistent with previous models of oxytocin in which flexibility is necessary for biological activity.     

NMR study of the solution conformation of actinomycin D.

The solution conformation of actinomycin D, the Gram-positive antibiotic and DNA-binding drug, has been determined by 1H-NMR in deuterated dimethyl sulfoxide. The structure determination is based on the experimental data set of NOE restraints. Four structures were obtained from the distance geometry and restrained molecular dynamics calculation. The resultant structures satisfy the experimental restraints very well. These structures are found to be compatible with the X-ray crystal structures.     

Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: An alternative method for generating a high-resolution solution structure

Summary A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTA-TAATG)·d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of <0.3 Å between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD<0.5 Å) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.Abbreviations MC Monte Carlo - rMC restrained Monte Carlo - MD molecular dynamics - rMD restrained molecular dynamics - DG distance geometry - EM energy minimization - 2D NOE two-dimensional nuclear Overhauser effect - DQF-COSY double-quantum-filtered correlation spectroscopy - RMSD root-mean-square deviation To whom correspondence should be addressed.     

NMR conformational studies on a synthetic peptide reproducing the [1-20] processing domain of the pro-ocytocin-neurophysin precursor

The combined use of several nuclear magnetic resonance and restrained molecular dynamics techniques allowed the formulation of a molecular model for the preferred solution conformation of a synthetic peptide reproducing the [1-20] processing domain of the pro-ocytocin-neurophysin precursor. In the model, the conformation of the 20-membered tocin ring, with the two Cys¹ and Cys⁶ residues bridged by a disulphide bond, is very close to that observed for isolated ocytocin in the solid state; in addition, a type II β-turn is postulated for the 7-10 segment of the acyclic tail containing the Lys¹¹-Arg¹² processing site, and connecting ocytocin to the neurophysin domain, while the C-terminal 13-20 segment of the molecule is believed to assume a helical structure. This particular structural organization could be important in participating as the favorable conformation for optimal substrate-enzyme active site recognition and processing by specific endoproteases. © 1996 John Wiley & Sons, Inc.     

Refinement of the solution structure of the ribonucleotide 5'r(GCAUGC)2: combined use of nuclear magnetic resonance and restrained molecular dynamics

The solution structure of the self-complementary hexamer 5'r(GCAUGC)2 is investigated by means of nuclear magnetic resonance spectroscopy and restrained molecular dynamics. The proton resonances are assigned in a sequential manner, and a set of 110 approximate interproton distance restraints are derived from the two-dimensional nuclear Overhauser enhancement spectra. These distances are used as the basis of a structure refinement by restrained molecular dynamics in which the experimental restraints are incorporated into the total energy function of the system in the form of effective potentials. Eight restrained molecular dynamics simulations are carried out, four starting from a structure with regular A-type geometry and four from one with regular B-type geometry. The atomic root mean square (rms) difference between the initial structures is 3.2 A. In the case of all eight simulations, convergence is achieved both globally and locally to a set of very similar A-type structures with an average atomic rms difference between them of 0.8 +/- 0.2 A. Further, the atomic rms differences between the restrained dynamics structures obtained by starting out from the same initial structures but with different random number seeds for the assignment of the initial velocities are the same as those between the restrained dynamics structures starting out from the two different initial structures. These results suggest that the restrained dynamics structures represent good approximations of the solution structure. The converged structures exhibit clear sequence-dependent variation in some of the helical parameters, in particular helix twist, roll, slide, and propellor twist.(ABSTRACT TRUNCATED AT 250 WORDS)     

Determination of the conformation of d(GGAAATTTCC)2 in solution by use of 1H NMR and restrained molecular dynamics

The conformation of the putative bent DNA d(GGAAATTTCC)2 in solution was studied by use of 1H NMR and restrained molecular dynamics. Most of the resonances were assigned sequentially. A total of 182 interproton distance restraints were determined from two-dimensional nuclear Overhauser effect spectra with short mixing times. Torsion angle restraints for each sugar moiety were determined by qualitative analysis of a two-dimensional correlated spectrum. Restrained molecular dynamics was carried out with the interproton distances and torsion angles incorporated into the total energy function of the system in the form of effective potential terms. As initial conformations for restrained molecular dynamics, classical A-DNA and B-DNA were adopted. The root mean square deviation (rmsd) between these two conformations is 5.5 A. The conformations obtained by use of restrained molecular dynamics are very similar to each other, the rmsd being 0.8 A. On the other hand, the conformations obtained by use of molecular dynamics without experimental restraints or restrained energy minimization depended heavily on the initial conformations, and convergence to a similar conformation was not attained. The conformation obtained by use of restrained molecular dynamics exhibits a few remarkable features. The second G residue takes on the BII conformation [Fratini, A. V., Kopka, M. L., Drew, H. R., & Dickerson, R. E. (1982) J. Biol. Chem. 257, 14686-14707] rather than the standard BI conformation. There is discontinuity of the sugar puckering between the eighth T and ninth C. The minor groove of the oligo(dA) tract is rather compressed. As a result, d(GGAAATTTCC)2 is bent.     

Backbone folding of the polypeptide cardiac stimulant anthopleurin-A determined by nuclear magnetic resonance, distance geometry and molecular dynamics

The solution conformation of the cardiac stimulatory sea anemone polypeptide anthopleurin-A has been characterised using distance geometry and restrained molecular dynamics calculations. A set of 253 approximate interproton distance restraints and 14 peptide backbone torsion angle restraints derived from two-dimensional 1H-NMR spectra at 500 MHz were used as input for these calculations. 13 structures generated by either metric matrix or variable target function distance geometry calculations were refined using energy minimisation and restrained molecular dynamics. The resulting structures contain a region of twisted antiparellel beta-sheet to which two separate regions of unordered chain are linked by three disulphide bonds. Two loops, one including Pro-41 and the other encompassing residues 10-18, are poorly defined by the NOE data.     

Determination of three-dimensional solution structure of waglerin I,a toxin from Trimeresurus wagleri,using 2D-NMR and molecular dynamics simulation

The solution conformation of a synthetic snake venom toxin waglerin I, has been determined by using proton nuclear magnetic resonance spectroscopy. By y a combination of various two-dimensional NMR techniques, the ¹H-NMR spectrum of waglerin I was completely assigned. A set of 247 interproton distance restraints was derived from nuclear Overhauser enhancement (NOE) measurements. These NOE constraints, in addition to the 2 dihedral angle restraints (from coupling constant measurements) and 7 ω torsion angle restraints for prolines, formed the basis of three-dimensional structure determined by molecular dynamics techniques. The 19 structures that were obtained satisfy the experimental restraints, and display small deviation from idealized covalent geometry. Analysis of converged structures indicates that the toxin has no special secondary structure. In the solution structure of waglerin I, the central ring region is well defined but the N- and C-termini possesses more disorder.     

Structure refinement of oligonucleotides by molecular dynamics with nuclear Overhauser effect interproton distance restraints: application to 5' d(C-G-T-A-C-G)2

The solution structure of the self-complementary DNA hexamer 5' d(C-G-T-A-C-G)2 is refined by restrained molecular dynamics in which 192 interproton distances, determined from pre-steady-state nuclear Overhauser enhancement measurements, are incorporated into the total energy of the system in the form of effective potentials. First the method is tested by applying an idealized set of distance restraints taken from classical B-DNA to a simulation starting off from A-DNA and vice versa. It is shown that in both cases the expected transition between A- and B-DNA occurs. Second, a set of restrained molecular dynamics calculations is carried out starting from both A- and B-DNA with the experimental interproton distances for 5' d(C-G-T-A-C-G)2 as restraints. Convergence to the same B-type structure is achieved with the interproton distances equal to the measured values within experimental error. The root-mean-square atomic difference between the two average restrained dynamics structures (less than 1 A) is approximately the same as the root-mean-square fluctuations of the atoms.     

Hydrogen--deuterium exchange kinetics of the amide protons of oxytocin studied by nuclear magnetic resonance

The contribution of intramolecular hydrogen bonding to the solution structure of oxytocin was evaluated by study of amide hydrogen exchange rates in D2O by Fourier transform 1H NMR spectroscopy. Resolution enhancement filtering was employed in the determination of individual pseudo-first-order rate constants. Apparent barriers to exchange of 0.5 and 0.6 kcal mol-1 were measured for Asn5 and Cys6 peptide NH, respectively. The slowing is best explained by steric hindrance to solvent access in the case of Asn5, while for the Cys6 participation in a weak intramolecular hydrogen bond is possible. Fourfold acceleration of base-catalyzed exchange was observed for Tyr2 NH; it is proposed that this is the result of electronic effects induced by hydrogen bonding of Cys1 C=0, either to Cys6 NH or to the N-terminal amino group. Exchange proceeds near the random coil limit for each of the remaining residues. Comparison with exchange data for the model tripeptide N-acetyl-L-prolyl-L-leucylglycinamide demonstrates no evidence of noncovalent association of the tocin ring with the tripeptide tail of the hormone.     

Refinement of the solution structure of the RNA-DNA hybrid 5'-[r(GCA)d(TGC)]2. Combined use of nuclear magnetic resonance and restrained molecular dynamics

The solution conformation of the self-complementary RNA-DNA hybrid hexamer 5'-[r(GCA)d(TGC)]2 is investigated by NMR spectroscopy and restrained molecular dynamics. The 1H-NMR spectrum is assigned in a sequential manner using two-dimensional homonuclear Hartmann-Hahn and nuclear Overhauser enhancement spectroscopy. From the latter a set of 178 approximate interproton distance restraints are determined and used as the basis of a structure refinement by restrained molecular dynamics. Eight independent calculations are carried out, four from a classical A-type geometry and four from a classical B-type one. Convergence is achieved to very similar A-type structures with an average atomic root mean square difference between them of 1.0 +/- 0.2 A. The converged structures exhibit variations in helical parameters similar to those found previously for the analogue RNA hexamer 5'-r(GCAUGC)2 [(1988) Biochemistry 27, 1735-1743].     

Molecular dynamics re-refinement of two different small RNA loop structures using the original NMR data suggest a common structure

Restrained molecular dynamics simulations are a robust, though perhaps underused, tool for the end-stage refinement of biomolecular structures. We demonstrate their utility-using modern simulation protocols, optimized force fields, and inclusion of explicit solvent and mobile counterions-by re-investigating the solution structures of two RNA hairpins that had previously been refined using conventional techniques. The structures, both domain 5 group II intron ribozymes from yeast ai5γ and Pylaiella littoralis, share a nearly identical primary sequence yet the published 3D structures appear quite different. Relatively long restrained MD simulations using the original NMR restraint data identified the presence of a small set of violated distance restraints in one structure and a possibly incorrect trapped bulge nucleotide conformation in the other structure. The removal of problematic distance restraints and the addition of a heating step yielded representative ensembles with very similar 3D structures and much lower pairwise RMSD values. Analysis of ion density during the restrained simulations helped to explain chemical shift perturbation data published previously. These results suggest that restrained MD simulations, with proper caution, can be used to "update" older structures or aid in the refinement of new structures that lack sufficient experimental data to produce a high quality result. Notable cautions include the need for sufficient sampling, awareness of potential force field bias (such as small angle deviations with the current AMBER force fields), and a proper balance between the various restraint weights.     

The solution structure of human transforming growth factor alpha

The solution structure of transforming growth factor alpha has been determined by a combination of high-resolution 1H-nuclear magnetic resonance and distance geometry and restrained molecular dynamics. The 382 restraints derived from the NMR experiments were used to calculate many distance geometry structures, which were then refined by restrained molecular mechanics. Five of these structures were further refined using a variety of methods. Comparison of independently measured parameters, such as calculated hydrogen bonding patterns and experimental amide exchange rates, have been used to evaluate the accuracy of the structures. Also, possible mechanisms to explain the pH-dependent conformational interconversion observed are suggested. Finally comparisons between this work and others on this topic have been made.     

Three-dimensional solution structure of an insulin dimer. A study of the B9(Asp) mutant of human insulin using nuclear magnetic resonance, distance geometry and restrained molecular dynamics.

The solution structure of the B9(Asp) mutant of human insulin has been determined by two-dimensional 1H nuclear magnetic resonance spectroscopy. Thirty structures were calculated by distance geometry from 451 interproton distance restraints based on intra-residue, sequential and long-range nuclear Overhauser enhancement data, 17 restraints on phi torsional angles obtained from 3JH alpha HN coupling constants, and the restraints from 17 hydrogen bonds, and the three disulphide bridges. The distance geometry structures were optimized using restrained molecular dynamics (RMD) and energy minimization. The average root-mean-square deviation for the best 20 RMD refined structures is 2.26 A for the backbone and 3.14 A for all atoms if the less well-defined N and C-terminal residues are excluded. The helical regions are better defined, with root-mean-square deviation values of 1.11 A for the backbone and 2.03 A for all atoms. The data analysis and the calculations show that B9(Asp) insulin, in water solution at the applied pH (1.8 to 1.9), is a well-defined dimer with no detectable difference between the two monomers. The association of the two monomers in the solution dimer is relatively loose as compared with the crystal dimer. The overall secondary and tertiary structures of the monomers in the 2Zn crystal hexamer is found to be preserved. The conformation-averaged NMR structures obtained for the monomer is close to the structure of molecule 1 in the hexamer of the 2Zn insulin crystal. However, minor, but significant deviations from this structure, as well as from the structure of monomeric insulin in solution, exist and are ascribed to the absence of the hexamer and crystal packing forces, and to the presence of monomer-monomer interactions, respectively. Thus, the monomer in the solution dimer shows a conformation similar to that of the crystal monomer in molecular regions close to the monomer-monomer interface, whereas it assumes a conformation similar to that of the solution structure of monomeric insulin in other regions, suggesting that B9(Asp) insulin adopts a monomer-like conformation when this is not inconsistent with the monomer-monomer arrangement in the dimer.     

Refinement of the solution structure of the DNA hexamer 5'd(GCATGC)2: combined use of nuclear magnetic resonance and restrained molecular dynamics

The solution structure of the self-complementary DNA hexamer 5'd(GCATGC)2 comprising the specific target site for the restriction endonuclease Sph 1 is investigated by using nuclear magnetic resonance spectroscopy and restrained molecular dynamics. All the nonexchangeable proton resonances are assigned sequentially, and from time-dependent nuclear Overhauser enhancement measurements a set of 158 approximate interproton distances are determined. These distances are used as the basis of a structure refinement using restrained molecular dynamics in which the interproton distances are incorporated into the total energy function of the system in the form of an effective potential term. Two restrained molecular dynamics simulations are carried out, starting from classical B- and A-DNA [atomic root mean square (rms) difference 3.3 A]. In both cases convergence is achieved to essentially identical structures satisfying the experimental restraints and having a root mean square difference of only 0.3 A between them, which is within the rms fluctuations of the atoms about their average positions. These results suggest that the restrained molecular dynamics structures represent reasonable approximations of the solution structure. The converged structures are of the B type and exhibit clear sequence-dependent variations of helical parameters, some of which follow Calladine's rules and can be attributed to the relief of interstrand purine-purine clash at adjacent base pairs. In addition, the converged restrained dynamics structures appear bent with a radius of curvature of approximately 20 A. This bending appears to be due almost entirely to the large positive base roll angles, particularly at the Pyr-Pur steps. Further, the global and local helix axes are not coincident, and the global helix axis represents a superhelical axis which the bent DNA, when extended into an "infinite" helix by repeated translation and rotation, wraps around.     

Tertiary structure of two-electron reduced Megasphaera elsdenii flavodoxin and some implications, as determined by two-dimensional 1H-NMR and restrained molecular dynamics

The tertiary structure of the non-crystallizable two-electron-reduced Megasphaera elsdenii flavodoxin (15 kDa, 137 amino acid residues) has been determined using nuclear Overhauser enhancement restraints extracted from two-dimensional 1H-NMR spectra. A tertiary structure satisfying the experimental restraints very well (maximum NOE violation of 66 pm) was obtained with use of restrained molecular dynamics, using 509 distance restraints (including one non-NOE) on a starting structure modeled from the crystal structure of one-electron-reduced Clostridium MP flavodoxin. The protein consists of a central parallel beta-sheet surrounded on both sides by two alpha-helices. The flavin is positioned at the periphery of the molecule. The tertiary structure of the protein is highly defined with the exception of the flavin. The latter is expected to result from performing the restrained molecular dynamics simulation without water molecules and without proper charges on the flavin. The flavin, including the phosphate, the ribityl side chain and the isoalloxazine ring, is solvent accessible under the experimental conditions used and evidenced by a two-dimensional amide exchange experiment. This accessibility is expected to be important in the redox potential regulation of the semiquinone/hydroquinone couple of the protein. The amide exchange against deuterons and several typical line shapes in the two-dimensional NMR spectra are consistent with the structure generated. The structure is discussed in detail.     

Solid-State NMR-Restrained Ensemble Dynamics of a Membrane Protein in Explicit Membranes

Solid-state NMR has been used to determine the structures of membrane proteins in native-like lipid bilayer environments. Most structure calculations based on solid-state NMR observables are performed using simulated annealing with restrained molecular dynamics and an energy function, where all nonbonded interactions are represented by a single, purely repulsive term with no contributions from van der Waals attractive, electrostatic, or solvation energy. To our knowledge, this is the first application of an ensemble dynamics technique performed in explicit membranes that uses experimental solid-state NMR observables to obtain the refined structure of a membrane protein together with information about its dynamics and its interactions with lipids. Using the membrane-bound form of the fd coat protein as a model membrane protein and its experimental solid-state NMR data, we performed restrained ensemble dynamics simulations with different ensemble sizes in explicit membranes. For comparison, a molecular dynamics simulation of fd coat protein was also performed without any restraints. The average orientation of each protein helix is similar to a structure determined by traditional single-conformer approaches. However, their variations are limited in the resulting ensemble of structures with one or two replicas, as they are under the strong influence of solid-state NMR restraints. Although highly consistent with all solid-state NMR observables, the ensembles of more than two replicas show larger orientational variations similar to those observed in the molecular dynamics simulation without restraints. In particular, in these explicit membrane simulations, Lys⁴⁰, residing at the C-terminal side of the transmembrane helix, is observed to cause local membrane curvature. Therefore, compared to traditional single-conformer approaches in implicit environments, solid-state NMR restrained ensemble simulations in explicit membranes readily characterize not only protein dynamics but also protein-lipid interactions in detail.     

Determination of the solution structures of domains II and III of protein G from Streptococcus by 1H nuclear magnetic resonance.

We have used 1H nuclear magnetic resonance spectroscopy to determine the solution structures of two small (61 and 64 residue) immunoglobulin G (IgG)-binding domains from protein G, a cell-surface protein from Streptococcus strain G148. The two domains differ in sequence by four amino acid substitutions, and differ in their affinity for some subclasses of IgG. The structure of domain II was determined using a total of 478 distance restraints, 31 phi and 9 chi 1 dihedral angle restraints; that of domain III was determined using a total of 445 distance restraints, 31 phi and 9 chi 1 dihedral angle restraints. A protocol which involved distance geometry, simulated annealing and restrained molecular dynamics was used to determine ensembles of 40 structures consistent with these restraints. The structures are found to consist of an alpha-helix packed against a four-stranded antiparallel-parallel-antiparallel beta-sheet. The structures of the two domains are compared to each other and to the reported structure of a similar domain from a protein G from a different strain of Streptococcus. We conclude that the difference in affinity of domains II and III for IgG is due to local changes in amino acid side-chains, rather than a more extensive change in conformation, suggesting that one or more of the residues which differ between them are directly involved in interaction with IgG.     

Solution structure of neurotoxin I from the sea anemone Stichodactyla helianthus. A nuclear magnetic resonance, distance geometry, and restrained molecular dynamics study

The three-dimensional structure of the sea anemone polypeptide Stichodactyla helianthus neurotoxin I in aqueous solution has been determined using distance geometry and restrained molecular dynamics simulations based on NMR data acquired at 500 MHz. A set of 470 nuclear Overhauser enhancement values was measured, of which 216 were used as distance restraints in the structure determination along with 15 dihedral angles derived from coupling constants. After restrained molecular dynamics refinement, the eight structures that best fit the input data form a closely related family. They describe a structure that consists of a core of twisted, four-stranded, antiparallel beta-sheet encompassing residues 1-3, 19-24, 29-34, and 40-47, joined by three loops, two of which are well defined by the NMR data. The third loop, encompassing residues 7-16, is poorly defined by the data and is assumed to undergo conformational averaging in solution. Pairwise root mean square displacement values for the backbone heavy atoms of the eight best structures are 1.3 +/- 0.2A when the poorly defined loop is excluded and 3.6 +/- 1.0A for all backbone atoms. Refinement using restrained molecular dynamics improved the quality of the structures generated by distance geometry calculations with respect to the number of nuclear Overhauser enhancements violated, the size of the total distance violations and the total potential energies of the structures. The family of structures for S. heliathus neurotoxin I is compared with structures of related sea anemone proteins that also bind to the voltage-gated sodium channel.     

Solution structure of a natural CPPC active site variant, the reduced form of thioredoxin h1 from poplar

Assignment of heteronuclear and homonuclear multidimensional NMR spectra permits determination of the first three-dimensional solution structure of a higher-plant thioredoxin h. The collection of 1906 distance restraints, 137 TALOS-derived dihedral restraints, and 66 hydrogen bonds was used in the restrained molecular dynamics protocol to calculate the structure of the reduced form of thioredoxin h1 from poplar with an atomic rmsd of 0.60 +/- 0.12 A. This enzyme exhibits an unusual active site with the sequence WCPPC and original properties in terms of stability and specificity. Compared to other known thioredoxin structures, thioredoxin h1 from poplar adopts the classical "Trx fold". Its atypical active site possesses a conformation similar to that of other common thioredoxins but appears to be more rigid. Moreover, the hydrogen bond network, stabilizing the in-core beta-sheet, is tighter than in Chlamydomonas reinhardtii, explaining the difference in thermostability.