The estrogen receptor α-selective agonist propyl pyrazole triol improves glucose tolerance in ob/ob mice: potential molecular mechanisms

The authors and journal apologise for an error in the above paper, which appeared in volume 199 part 2, pages 275­286. The error relates to Fig. 10, given on page 283.

The estrogen receptor {alpha}-selective agonist propyl pyrazole triol improves glucose tolerance in ob/ob mice; potential molecular mechanisms.

The aim of this study was to validate the role of estrogen receptor alpha (ERalpha) signaling in the regulation of glucose metabolism, and to compare the molecular events upon treatment with the ERalpha-selective agonist propyl pyrazole triol (PPT) or 17beta-estradiol (E(2)) in ob/ob mice. Female ob/ob mice were treated with PPT, E(2) or vehicle for 7 or 30 days. Intraperitoneal glucose and insulin tolerance tests were performed, and insulin secretion was determined from isolated islets. Glucose uptake was assayed in isolated skeletal muscle and adipocytes. Gene expression profiling in the liver was performed using Affymetrix microarrays, and the expression of selected genes was studied by real-time PCR analysis. PPT and E(2) treatment improved glucose tolerance and insulin sensitivity. Fasting blood glucose levels decreased after 30 days of PPT and E(2) treatment. However, PPT and E(2) had no effect on insulin secretion from isolated islets. Basal and insulin-stimulated glucose uptake in skeletal muscle and adipose tissue were similar in PPT and vehicle-treated ob/ob mice. Hepatic lipid content was decreased after E(2) treatment. In the liver, treatment with E(2) and PPT increased and decreased the respective expression levels of the transcription factor signal transducer and activator of transcription 3, and of glucose-6-phosphatase. In summary, our data demonstrate that PPT exerts anti-diabetic effects, and these effects are mediated via ERalpha.

Structural role of conserved Asn179 in the short-chain dehydrogenase/reductase scaffold.

Short-chain dehydrogenases/reductases (SDR) constitute a large family of enzymes found in all forms of life. Despite a low level of sequence identity, the three-dimensional structures determined display a nearly superimposable alpha/beta folding pattern. We identified a conserved asparagine residue located within strand betaF and analyzed its role in the short-chain dehydrogenase/reductase architecture. Mutagenetic replacement of Asn179 by Ala in bacterial 3beta/17beta-hydroxysteroid dehydrogenase yields a folded, but enzymatically inactive enzyme, which is significantly more resistant to denaturation by guanidinium hydrochloride. Crystallographic analysis of the wild-type enzyme at 1.2-A resolution reveals a hydrogen bonding network, including a buried and well-ordered water molecule connecting strands betaE to betaF, a common feature found in 16 of 21 known three-dimensional structures of the family. Based on these results, we hypothesize that in mammalian 11beta-hydroxysteroid dehydrogenase the essential Asn-linked glycosylation site, which corresponds to the conserved segment, displays similar structural features and has a central role to maintain the SDR scaffold.

How transcriptional activators bind target proteins.

The product of the proto-oncogene c-myc influences many cellular processes through the regulation of specific target genes. Through its transactivation domain (TAD), c-Myc protein interacts with several transcription factors, including TATA-binding protein (TBP). We present data that suggest that in contrast to some other transcriptional activators, an extended length of the c-Myc TAD is required for its binding to TBP. Our data also show that this interaction is a multistep process, in which a rapidly forming low affinity complex slowly converts to a more stable form. The initial complex formation results from ionic or polar interactions, whereas the slow conversion to a more stable form is hydrophobic in nature. Based on our results, we suggest two alternative models for activation domain/target protein interactions, which together provide a single universal paradigm for understanding activator-target factor interactions.

Structure, dynamics and electrostatics of the active site of glutaredoxin 3 from Escherichia coli: comparison with functionally related proteins.

The chemistry of active-site cysteine residues is central to the activity of thiol-disulfide oxidoreductases of the thioredoxin superfamily. In these reactions, a nucleophilic thiolate is required, but the associated pK(a) values differ vastly in the superfamily, from less than 4 in DsbA to greater than 7 in Trx. The factors that stabilize this thiolate are, however, not clearly established. The glutaredoxins (Grxs), which are members of this superfamily, contain a Cys-Pro-Tyr-Cys motif in their active site. In reduced Grxs, the pK(a) of the N-terminal active-site nucleophilic cysteine residue is lowered significantly, and the stabilization of the corresponding thiolate is expected to influence the redox potential of these enzymes. Here, we use a combination of long molecular dynamics (MD) simulations, pK(a) calculations, and experimental investigations to derive the structure and dynamics of the reduced active site from Escherichia coli Grx3, and investigate the factors that stabilize the thiolate. Several different MD simulations converged toward a consensus conformation for the active-site cysteine residues (Cys11 and Cys14), after a number of local conformational changes. Key features of the model were tested experimentally by measurement of NMR scalar coupling constants, and determination of pK(a) values of selected residues. The pK(a) values of the Grx3 active-site residues were calculated during the MD simulations, and support the underlying structural model. The structure of Grx3, in combination with the pK(a) calculations, indicate that the pK(a) of the N-terminal active-site cysteine residue in Grx3 is intermediate between that of its counterpart in DsbA and Trx. The pK(a) values in best agreement with experiment are obtained with a low (<4) protein dielectric constant. The calculated pK(a) values fluctuate significantly in response to protein dynamics, which underscores the importance of the details of the underlying structures when calculating pK(a) values. The thiolate of Cys11 is stabilized primarily by direct hydrogen bonding with the amide protons of Tyr13 and Cys14 and the thiol proton of Cys14, rather than by long-range interactions from charged groups or from a helix macrodipole. From the comparison of reduced Grx3 with other members of the thioredoxin superfamily, a unifying theme for the structural basis of thiol pK(a) differences in this superfamily begins to emerge.

NMR structure of oxidized glutaredoxin 3 from Escherichia coli.

A high precision NMR structure of oxidized glutaredoxin 3 [C65Y] from Escherichia coli has been determined. The conformation of the active site including the disulphide bridge is highly similar to those in glutaredoxins from pig liver and T4 phage. A comparison with the previously determined structure of glutaredoxin 3 [C14S, C65Y] in a complex with glutathione reveals conformational changes between the free and substrate-bound form which includes the sidechain of the conserved, active site tyrosine residue. In the oxidized form this tyrosine is solvent exposed, while it adopts a less exposed conformation, stabilized by hydrogen bonds, in the mixed disulfide with glutathione. The structures further suggest that the formation of a covalent linkage between glutathione and glutaredoxin 3 is necessary in order to induce these structural changes upon binding of the glutathione peptide. This could explain the observed low affinity of glutaredoxins for S-blocked glutathione analogues, in spite of the fact that glutaredoxins are highly specific reductants of glutathione mixed disulfides.

Direct NMR observation of the Cys-14 thiol proton of reduced Escherichia coli glutaredoxin-3 supports the presence of an active site thiol-thiolate hydrogen bond.

The active site of Escherichia coli glutaredoxin-3 (Grx3) consists of two redox active cysteine residues in the sequence -C11-P-Y-C14-H-. The 1H NMR resonance of the cysteine thiol proton of Cys-14 in reduced Grx3 is observed at 7.6 ppm. The large downfield shift and NOEs observed with this thiol proton resonance suggest the presence of a hydrogen bond with the Cys-11 thiolate, which is shown to have an abnormally low pKa value. A hydrogen bond would also agree with activity data of Grx3 active site mutants. Furthermore, the activity is reduced in a Grx3 H15V mutant, indicating electrostatic contributions to the stabilization of the Cys-11 thiolate.

NMR structure of Escherichia coli glutaredoxin 3-glutathione mixed disulfide complex: implications for the enzymatic mechanism.

Glutaredoxins (Grxs) catalyze reversible oxidation/reduction of protein disulfide groups and glutathione-containing mixed disulfide groups via an active site Grx-glutathione mixed disulfide (Grx-SG) intermediate. The NMR solution structure of the Escherichia coli Grx3 mixed disulfide with glutathione (Grx3-SG) was determined using a C14S mutant which traps this intermediate in the redox reaction. The structure contains a thioredoxin fold, with a well-defined binding site for glutathione which involves two intermolecular backbone-backbone hydrogen bonds forming an antiparallel intermolecular beta-bridge between the protein and glutathione. The solution structure of E. coli Grx3-SG also suggests a binding site for a second glutathione in the reduction of the Grx3-SG intermediate, which is consistent with the specificity of reduction observed in Grxs. Molecular details of the structure in relation to the stability of the intermediate and the activity of Grx3 as a reductant of glutathione mixed disulfide groups are discussed. A comparison of glutathione binding in Grx3-SG and ligand binding in other members of the thioredoxin superfamily is presented, which illustrates the highly conserved intermolecular interactions in this protein family.

Amyloid beta-peptide polymerization studied using fluorescence correlation spectroscopy.

BACKGROUND: The accumulation of fibrillar deposits of amyloid beta-peptide (Abeta) in brain parenchyma and cerebromeningeal blood vessels is a key step in the pathogenesis of Alzheimer's disease. In this report, polymerization of Abeta was studied using fluorescence correlation spectroscopy (FCS), a technique capable of detecting small molecules and large aggregates simultaneously in solution. RESULTS: The polymerization of Abeta dissolved in Tris-buffered saline, pH 7.4, occurred above a critical concentration of 50 microM and proceeded from monomers/dimers into two discrete populations of large aggregates, without any detectable amount of oligomers. The aggregation showed very high cooperativity and reached a maximum after 40 min, followed by an increase in the amount of monomers/dimers and a decrease in the size of the large aggregates. Electron micrographs of samples prepared at the time for maximum aggregation showed a mixture of an amorphous network and short diffuse fibrils, whereas only mature amyloid fibrils were detected after one day of incubation. The aggregation was reduced when Abeta was incubated in the presence of Abeta ligands, oligopeptides previously shown to inhibit fibril formation, and aggregates were partly dissociated after the addition of the ligands. CONCLUSIONS: The polymerization of Abeta is a highly cooperative process in which the formation of very large aggregates precedes the formation of fibrils. The entire process can be inhibited and, at least in early stages, partly reversed by Abeta ligands.

Thermal unfolding of small proteins with SH3 domain folding pattern.

The thermal unfolding of three SH3 domains of the Tec family of tyrosine kinases was studied by differential scanning calorimetry and CD spectroscopy. The unfolding transition of the three protein domains in the acidic pH region can be described as a reversible two-state process. For all three SH3 domains maximum stability was observed in the pH region 4.5 < pH < 7.0 where these domains unfold at temperatures of 353K (Btk), 342K (Itk), and 344K (Tec). At these temperatures an enthalpy change of 196 kJ/mol, 178 kJ/mol, and 169 kJ/mol was measured for Btk-, Itk-, and Tec-SH3 domains, respectively. The determined changes in heat capacity between the native and the denatured state are in an usual range expected for small proteins. Our analysis revealed that all SH3 domains studied are only weakly stabilized and have free energies of unfolding which do not exceed 12-16 kJ/mol but show quite high melting temperatures. Comparing unfolding free energies measured for eukaryotic SH3 domains with those of the topologically identical Sso7d protein from the hyperthermophile Sulfolobus solfataricus, the increased melting temperature of the thermostable protein is due to a broadening as well as a significant lifting of its stability curve. However, at their physiological temperatures, 310K for mesophilic SH3 domains and 350K for Sso7d, eukaryotic SH3 domains and Sso7d show very similar stabilities.

Conformation-dependent antibacterial activity of the naturally occurring human peptide LL-37.

The influence of ion composition, pH, and peptide concentration on the conformation and activity of the 37-residue human antibacterial peptide LL-37 has been studied. At micromolar concentration in water, LL-37 exhibits a circular dichroism spectrum consistent with a disordered structure. The addition of 15 mM HCO3-, SO42-, or CF3CO2- causes the peptide to adopt a helical structure, with approximately equal efficiency, while 160 mM Cl- is less efficient. A cooperative transition from disordered to helical structure is observed as the peptide concentration is increased, consistent with formation of an oligomer. The extent of alpha-helicity correlates with the antibacterial activity of LL-37 against both Gram-positive and Gram-negative bacteria. Two homologous peptides, FF-33 and SK-29, containing 4 and 8 residue deletions at the N terminus, respectively, require higher concentrations of anions for helix formation and are less active than LL-37 against Escherichia coli D21. Below pH 5, the helical content of LL-37 gradually decreases, and at pH 2 it is entirely disordered. In contrast, the helical structure is retained at pH over 13. The minimal inhibitory concentration of LL-37 against E. coli is 5 microM, and at 13-25 microM the peptide is cytotoxic against several eukaryotic cells. In solutions containing the ion compositions of plasma, intracellular fluid, or interstitial fluid, LL-37 is helical, and hence it could pose a danger to human cells upon release. However, in the presence of human serum, the antibacterial and the cytotoxic activities of LL-37 are inhibited.

A proteinaceous gene regulatory thermometer in Salmonella.

Novel utilization of the coiled-coil motif is presented that enables TlpA, an autoregulatory repressor protein in Salmonella, to sense temperature shifts directly and thereby to modulate the extent of transcription repression. Salmonella cells shifted to higher temperatures, such as those encountered at host entry, showed derepressed tlpA activity. tlpA::lacZ fusions indicated that the promoter itself is insensitive to thermal shifts and that transcription control was exerted by the autorepressor TlpA only. In vitro studies with highly purified TlpA showed concentration and temperature dependence for both fully folded conformation and function, indicating that the thermosensing in TlpA is based on monomer-to-coiled-coil equilibrium.

Active site directed mutagenesis of 3 beta/17 beta-hydroxysteroid dehydrogenase establishes differential effects on short-chain dehydrogenase/reductase reactions.

Mutagenetic replacements of conserved residues within the active site of the short-chain dehydrogenase/reductase (SDR) superfamily were studied using prokaryotic 3 beta/17 beta-hydroxysteroid dehydrogenase (3 beta/17 beta-HSD) from Comamonas testosteroni as a model system. The results provide novel data to establish Ser 138 as a member of a catalytically important "triad" of residues also involving Tyr151 and Lys155. A Ser-->Ala exchange at position 138 results in an almost complete (> 99.9%) loss of enzymatic activity, which is not observed with a Ser-->Thr replacement. This indicates that an essential factor for catalysis is the ability of side chain 138 to form hydrogen bond interactions. Mutations in the NAD(H) binding region, in strands beta A, beta D, and adjacent turns, reveal two additional residues, Thr12 and Asn87, which are important for correct binding of the coenzyme and with a differential effect on the reactions catalyzed. Thus, mutation of Thr12 to Ala results in a complete loss of the 3 beta-dehydrogenase activity, whereas the 3-oxoreductase activity remains unchanged. On the other hand, a T12S substitution yields a protein with unaltered catalytic constants for both reactions, revealing that a specific hydrogen bond is critical for the dehydrogenase activity. Our interpretation of the available crystal structure of 3 alpha/20 beta-HSD from Streptomyces hydrogenans suggests a hydrogen bond in that enzyme between the Thr12 side chain and the backbone NH of Asn87 rather than the coenzyme, indicating that this hydrogen bond to the beta D strand might determine a crucial difference between the reductive and the oxidative reaction types. Similarly, mutation of Asn87 to Ala results in an 80% reduction of kcat/Km in the dehydrogenase direction but also unchanged 3-oxoreductase properties. It appears that the binding of NAD+ to the protein is influenced by local structural changes involving strand beta D and turn beta A to alpha B.

Redox potentials of glutaredoxins and other thiol-disulfide oxidoreductases of the thioredoxin superfamily determined by direct protein-protein redox equilibria.

Glutaredoxins belong to the thioredoxin superfamily of structurally similar thiol-disulfide oxidoreductases catalyzing thiol-disulfide exchange reactions via reversible oxidation of two active-site cysteine residues separated by two amino acids (CX1X2C). Standard state redox potential (E degrees ') values for glutaredoxins are presently unknown, and use of glutathione/glutathione disulfide (GSH/GSSG) redox buffers for determining E degrees ' resulted in variable levels of GSH-mixed disulfides. To overcome this complication, we have used reverse-phase high performance liquid chromatography to separate and quantify the oxidized and reduced forms present in the thiol-disulfide exchange reaction at equilibrium after mixing one oxidized and one reduced protein. This allowed for direct and quantitative pair-wise comparisons of the reducing capacities of the proteins and mutant forms. Equilibrium constants from pair-wise reaction with thioredoxin or its P34H mutant, which have accurately determined E degrees ' values from their redox equilibrium with NADPH catalyzed by thioredoxin reductase, allowed for transformation into standard state values. Using this new procedure, the standard state redox potentials for the Escherichia coli glutaredoxins 1 and 3, which contain identical active site sequences CPYC, were found to be E degrees ' = -233 and -198 mV, respectively. These values were confirmed independently by using the thermodynamic linkage between the stability of the disulfide bond and the stability of the protein to denaturation. Comparison of calculated E degrees ' values from a number of proteins ranging from -270 mV for E. coli Trx to -124 mV for DsbA obtained using this method with those determined using glutathione redox buffers provides independent confirmation of the standard state redox potential of glutathione as -240 mV. Determining redox potentials through direct protein-protein equilibria is of general interest as it overcomes errors in determining redox potentials calculated from large equilibrium constants with the strongly reducing NADPH or by accumulating mixed disulfides with GSH.

Thermal unfolding of the DNA-binding protein Sso7d from the hyperthermophile Sulfolobus solfataricus.

Thermal unfolding of the small hyperthermophilic DNA-binding protein Sso7d was studied by circular dichroism spectroscopy and differential scanning calorimetry. The unfolding transition can be described by a reversible two state process. Maximum stability was observed in the region between pH 4.5 and 7.0 where Sso7d unfolds with a melting temperature between 370.8 to 371.9 K and an unfolding enthalpy between 62.9 and 65.4 kcal/mol. The heat capacity differences between the native and the heat denatured states obtained by differential scanning calorimetry (620 cal/(molK)) and circular dichroism spectroscopy (580 cal/(mol K)) resulted in comparable values. The thermodynamic reason for the high melting temperature of Sso7d is the shallow stability curve with a broad free energy maximum, corresponding to the relatively small heat capacity change which was obtained. The calculated stability curve shows that Sso7d has, despite of its high melting temperature, an only moderate intrinsic stability, which reaches its maximum (approximately 7 kcal/mol) at 282 K. Sso7d is particularly poorly stabilized (approximately 1 kcal/mol) at the maximum physiological growth temperature of Sulfolobus solfataricus. Sso7d has furthermore untypically low specific enthalpy (0.99 kcal/(mol residue)) and entropy (2.99 cal/(mol K)) values at convergence temperatures. No significant differences in thermal stability of the partially methylated Sso7d from Sulfolobus solfataricus and the cloned non-methylated form of the protein expressed in Escherichia coli were observed.

Solution structure of a naturally-occurring zinc-peptide complex demonstrates that the N-terminal zinc-binding module of the Lasp-1 LIM domain is an independent folding unit.

The three-dimensional solution structure of the 1:1 complex between the synthetic peptide ZF-1 and zinc was determined by 1H NMR spectroscopy. The peptide, initially isolated from pig intestines, is identical in sequence to the 30 N-terminal amino acid residues of the human protein Lasp-1 belonging to the LIM domain protein family. The final set of 20 energy-refined NMR conformers has an average rmsd relative to the mean structure of 0.55 A for the backbone atoms of residues 3-30. Calculations without zinc atom constraints unambiguously identified Cys 5, Cys 8, His 26, and Cys 29 as the zinc-coordinating residues. LIM domains consist of two sequential zinc-binding modules and the NMR structure of the ZF-1-zinc complex is the first example of a structure of an isolated module. Comparison with the known structures of the N-terminal zinc-binding modules of both the second LIM domain of chicken CRP and rat CRIP with which ZF-1 shares 50% and 43% sequence identity, respectively, supports the notion that the zinc-binding modules of the LIM domain have a conserved structural motif and identifies local regions of structural diversity. The similarities include conserved zinc-coordinating residues, a rubredoxin knuckle involving Cys 5 and Cys 8, and the coordination of the zinc ion by histidine N delta in contrast to the more usual coordination by N epsilon observed for other zinc-finger domains. The present structure determination of the ZF-1-zinc complex establishes the N-terminal half of a LIM domain as an independent folding unit. The structural similarities of N- and C-terminal zinc-binding modules of the LIM domains, despite limited sequence identity, lead to the proposal of a single zinc-binding motif in LIM domains. The coordinates are available from the Brookhaven protein data bank, entry 1ZFO.

DNA binding exerted by a bacterial gene regulator with an extensive coiled-coil domain.

Although quite common in the eukaryotic cell, bacterial proteins with an extensive coiled-coil domain are still relatively rare. One of the few thus far documented examples, TlpA from Salmonella typhimurium, is characterized by a remarkably long (250 amino acids) alpha-helical coiled-coil domain. Herein, we demonstrate that TlpA is a novel sequence-specific DNA-binding protein. Several tlpA deletion mutants have been constructed, and their corresponding protein products were purified and tested for DNA binding. Two of the mutant proteins were shown to be deficient in DNA binding. Both mutants were analyzed by circular dichroism and electron microscopy, supporting the notion that mutant proteins wre shown to be deficient in DNA binding. Both mutants were analyzed by circular dichroism and electron microscopy, supporting the notion that mutant proteins were largely intact despite lacking the amino acid residues necessary for DNA binding. In vivo studies with transcriptional tlpA-lacZ fusions demonstrated that TlpA acts as a repressor. Using the repressor phenotype as a readout, the chain exchange previously described in vitro could also be confirmed in vivo. We believe the coiled-coil domain acts not only as a dimerization interface but could also serve a role as a flexible modulator of the protein-DNA interaction.

Glutaredoxin-3 from Escherichia coli. Amino acid sequence, 1H AND 15N NMR assignments, and structural analysis.

The primary and secondary structure of glutaredoxin-3 (Grx3), a glutathione-disulfide oxidoreductase from Escherichia coli, has been determined. The amino acid sequence of Grx3 consists of 82 residues and contains a redox-active motif, Cys-Pro-Tyr-Cys, typical of the glutaredoxin family. Sequence comparison reveals a homology (33% identity) to that of glutaredoxin-1 (Grx1) from E. coli as well as to other members of the thioredoxin superfamily. In addition to the active site cysteine residues, Grx3 contains one additional cysteine (Cys65) corresponding to one of the two non-active site (or structural) cysteine residues present in mammalian glutaredoxins. The sequence-specific 1H and 15N nuclear magnetic resonance assignments of reduced Grx3 have been obtained. From a combined analysis of chemical shifts, 3JHNalpha coupling constants, sequential and medium range NOEs, and amide proton exchange rates, the secondary structure of reduced Grx3 was determined and found to be very similar to that inferred from amino acid sequence comparison to homologous proteins. The consequences of the proposed structural similarity to Grx1 are that Grx3, while possessing a largely intact GSH binding cleft, would have a very different spatial distribution of charged residues, most notably surrounding the active site cysteine residues and occurring in the proposed hydrophobic protein-protein interaction area. These differences may contribute to the observed very low Kcat of Grx3 as a reductant of insulin disulfides or as a hydrogen donor for ribonucleotide reductase. Thus, despite an identical active site disulfide motif and a similar secondary structure and tertiary fold, Grx3 and Grx1 display large functional differences in in vitro protein disulfide oxido-reduction reactions.

Conformational sampling by NMR solution structures calculated with the program DIANA evaluated by comparison with long-time molecular dynamics calculations in explicit water.

The NMR solution structure of bovine pancreatic trypsin inhibitor (BPTI) obtained by distance geometry calculations with the program DIANA is compared with groups of conformers generated by molecular dynamics (MD) simulations in explicit water at ambient temperature and pressure. The MD simulations started from a single conformer and were free or restrained either by the experimental NOE distance restraints or by time-averaged restraints; the groups of conformers were collected either in 10 ps intervals during 200 ps periods of simulation, or in 50 ps intervals during a 1 ns period of simulation. Overall, these comparisons show that the standard protein structure determination protocol with the program DIANA provides a picture of the protein structure that is in agreement with MD simulations using "realistic" potential functions over a nanosecond timescale. For well-constrained molecular regions there is a trend in the free MD simulation of duration 1 ns that the sampling of the conformation space is slightly increased relative to the DIANA calculations. In contrast, for surface-exposed side-chains that are less extensively constrained by the NMR data, the DIANA conformers tend to sample larger regions of conformational space than conformers selected from any of the MD trajectories. Additional insights into the behavior of surface side-chains come from comparison of the MD runs of 200 ps or 1 ns duration. In this time range the sampling of conformation space by the protein surface depends strongly on the length of the simulation, which indicates that significant side-chain transitions occur on the nanosecond timescale and that much longer simulations will be needed to obtain statistically significant data on side-chain dynamics.

Structure and internal dynamics of the bovine pancreatic trypsin inhibitor in aqueous solution from long-time molecular dynamics simulations.

Structural and dynamic properties of bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution are investigated using two molecular dynamics (MD) simulations: one of 1.4 ns length and one of 0.8 ns length in which atom-atom distance bounds derived from NMR spectroscopy are included in the potential energy function to make the trajectory satisfy these experimental data more closely. The simulated properties of BPTI are compared with crystal and solution structures of BPTI, and found to be in agreement with the available experimental data. The best agreement with experiment was obtained when atom-atom distance restraints were applied in a time-averaged manner in the simulation. The polypeptide segments found to be most flexible in the MD simulations coincide closely with those showing differences between the crystal and solution structures of BPTI.