Effect of hydrostatic pressure on the solvent in crystals of hen egg-white lysozyme

The mass density of protein crystals can be measured in Ficoll gradients as a function of hydrostatic pressure. Carbon tetrachloride-toluene mixtures provide convenient density markers, and the compressibility of these standards is reported. Measurements on tetragonal crystals of hen egg-white lysozyme yielded densities at room temperature of 1.2367(+/- 0.0010) g cm-3 at 1 atm and 1.2586(+/- 0.0017) g cm-3 at 1000 atm (1 atm = 101,325 Pa). When combined with the unit cell dimensions at these two pressures these values lead to an estimated compression (fractional change in volume) of the crystal solvent at 1000 atm of 0.0369(+/- 0.0054). This value is comparable to that of a 0.7 M solution of NaCl. From an approximate estimate of the Donnan effect for the crystal in the 1.4 M-NaCl mother liquor, the crystal solvent contains 0.8 M-Na+ and 2.5 M-Cl-. It is concluded that the compressibility of solvent in lysozyme crystals is, within experimental error, the same as bulk solvent and does not exhibit the dramatically altered compressibility expected of an ice or glass-like solid. The crystallographically observable water sites, 151 at 1 atm and 163 at 1000 atm, showed a tendency to increase the number of hydrogen bonds made to other water sites at the expense of hydrogen bonds made to protein. The explanation for this phenomenon is presently unknown. Water sites that occur in both structures tend to have comparable temperature factors and show some tendency to follow the pressure-induced changes in protein atom positions. The compression expected for the water molecules themselves is too small to be observable at the resolution of the X-ray data collected in this study.     

Conformational deformation in deoxymyoglobin by hydrostatic pressure

Pressure effect on the equilibrium conformation in sperm whale deoxymyoglobin and its volume fluctuation are studied by the normal mode analysis and strain tensor analysis. The pressure-induced deformation of interhelix regions are found to be remarkably more compressed than the other parts of the molecule. The intrahelix compressibility is shown to be relatively small. We also calculate the compressibility of the three hydrophobic clusters, located at the bottom, distal, and proximal side of the heme. Its value is found to decrease in the indicated order. The average compressibility of these hydrophobic clusters is less than the average interhelix compressibility, even though there are large cavities in these clusters. In order to study overall deformation, we define a linear compressibility and calculate it for all pairs of C^αatoms. The pressure-induced deformation is observed to be heterogeneous also in this analysis. The calculated root-mean-square displacement of the constituent atoms in the equilibrium conformation at 1,000 atm from those at 0 atm is 0.12 Å, which is much smaller in magnitude than the average value of the atomic fluctuations at room temperature. In the water solvent, the volume excluded by the protein molecule in the equilibrium conformation is reduced by 0.9%, when the pressure is raised from 0 to 1,000 atm. The calculated magnitude of the root-mean-square volume fluctuation is 0.3% of the excluded volume at room temperature. The square of the volume fluctuation is given as a sum of contributions from individual normal modes. Contributions from low frequency normal modes are found to dominate over those from higher frequency normal modes. The estimated value of the isothermal compressibility of deoxymyoglobin is 9.37 × 10^?12 cm²/dyn. © 1993 Wiley-Liss, Inc.     

Molecular dynamics simulation of solvated protein at high pressure.

We have completed a molecular dynamics simulation of protein (bovine pancreatic trypsin inhibitor, BPTI) in solution at high pressure (10 kbar). The structural and energetic effects of the application of high pressure to solvated protein are analyzed by comparing the results of the high-pressure simulation with a corresponding simulation at low pressure. The volume of the simulation cell containing one protein molecule plus 2943 water molecules decreases by 24.7% at high pressure. This corresponds to a compressibility for the protein solution of beta = 1.8 x 10(-2) kbar-1. The compressibility of the protein is estimated to be about one-tenth that of bulk water, while the protein hydration layer water is found to have a greater compressibility as compared to the bulk, especially for water associated with hydrophobic groups. The radius of gyration of BPTI decreases by 2% and there is a one third decrease in the protein backbone atomic fluctuations at high pressure. We have analyzed pressure effects on the hydration energy of the protein. The total hydration energy is slightly (4%) more favorable at high pressure even though the surface accessibility of the protein has decreased by a corresponding amount. Large pressure-induced changes in the structure of the hydration shell are observed. Overall, the solvation shell waters appear more ordered at high pressure; the pressure-induced ordering is greatest for nonpolar surface groups. We do not observe evidence of pressure-induced unfolding of the protein over the 100-ps duration of the high-pressure simulation. This is consistent with the results of high-pressure optical experiments on BPTI.(ABSTRACT TRUNCATED AT 250 WORDS)     

Determination of the three-dimensional solution structure of barnase using nuclear magnetic resonance spectroscopy

The solution conformation of the ribonuclease barnase has been determined by using 1H nuclear magnetic resonance (NMR) spectroscopy. The 20 structures were calculated by using 853 interproton distance restraints obtained from analyses of two-dimensional nuclear Overhauser spectra, 72 phi and 53 chi 1 torsion angle restraints, and 17 hydrogen-bond distance restraints. The calculated structures contain two alpha-helices (residues 6-18 and 26-34) and a five-stranded antiparallel beta-sheet (residues 50-55, 70-75, 85-91, 94-101, and 105-108). The core of the protein is formed by the packing of one of the alpha-helices (residues 6-18) onto the beta-sheet. The average RMS deviation between the calculated structures and the mean structure is 1.11 A for the backbone atoms and 1.75 A for all atoms. The protein is least well-defined in the N-terminal region and in three large loops. When these regions are excluded, the average RMS deviation between the calculated structures and the mean structure for residues 5-34, 50-56, 71-76, 85-109 is 0.62 A for the backbone atoms and 1.0 A for all atoms. The NMR-derived structure has been compared with the crystal structure of barnase [Mauguen et al. (1982) Nature (London) 297, 162-164].     

Pressure response of protein backbone structure. Pressure-induced amide 15N chemical shifts in BPTI.

The effect of pressure on amide 15N chemical shifts was studied in uniformly 15N-labeled basic pancreatic trypsin inhibitor (BPTI) in 90%1H2O/10%2H2O, pH 4.6, by 1H-15N heteronuclear correlation spectroscopy between 1 and 2,000 bar. Most 15N signals were low field shifted linearly and reversibly with pressure (0.468 +/- 0.285 ppm/2 kbar), indicating that the entire polypeptide backbone structure is sensitive to pressure. A significant variation of shifts among different amide groups (0-1.5 ppm/2 kbar) indicates a heterogeneous response throughout within the three-dimensional structure of the protein. A tendency toward low field shifts is correlated with a decrease in hydrogen bond distance on the order of 0.03 A/2 kbar for the bond between the amide nitrogen atom and the oxygen atom of either carbonyl or water. The variation of 15N shifts is considered to reflect site-specific changes in phi, psi angles. For beta-sheet residues, a decrease in psi angles by 1-2 degrees/2 kbar is estimated. On average, shifts are larger for helical and loop regions (0.553 +/- 0.343 and 0.519 +/- 0.261 ppm/2 kbar, respectively) than for beta-sheet (0.295 +/- 0.195 ppm/2 kbar), suggesting that the pressure-induced structural changes (local compressibilities) are larger in helical and loop regions than in beta-sheet. Because compressibility is correlated with volume fluctuation, the result is taken to indicate that the volume fluctuation is larger in helical and loop regions than in beta-sheet. An important aspect of the volume fluctuation inferred from pressure shifts is that they include motions in slower time ranges (less than milliseconds) in which many biological processes may take place.     

Pressure dependence of human fibrinogen correlated to the conformational alpha-helix to beta-sheet transition: an Fourier transform infrared study microspectroscopic study

We used Fourier transform infrared (FTIR) microspectroscopy to investigate pressure-induced conformational changes in secondary structure of fibrinogen (FBG). Solid state FBG was compressed on a KBr pellet (1KBr method) or between two KBr pellets (2KBr method). The peak positions of the original and second-derivative ir spectra of compressed FBG samples prepared by the 1KBr method were similar to FBG sample without pressure. When FBG was prepared by the 2KBr method and pressure was increased up to 400 kg/cm(2), peaks at 1625 (intermolecular beta-sheet) and 1611 (beta-sheet aggregates structure and/or the side-chain absorption of the tyrosine residues) cm(-1) were enhanced. The peaks near 1661 (beta-sheet) and 1652 (alpha-helix) cm(-1) also exhibited a marked change with pressure. A linear correlation was found between the peak intensity ratio of 1611/1652 cm(-1) (r = 0.9879) or 1625/1652 cm(-1) (r = 0.9752) and applied pressure. The curve-fitted compositional changes in secondary structure of FBG also indicate that the composition of the alpha-helix structure (1657-1659 cm(-1)) was gradually reduced with the increase in compression pressure, but the composition of the beta-sheet structure (1681, 1629, and 1609 cm(-1)) gradually increased. This indicates that pressure-induced conformational changes in FBG include not only transformations from alpha-helix to beta-sheet structure, but also unfolding and denaturation of FBG and the formation of aggregates.     

Determination of the spatial structure of insectotoxin 15A from Buthus erpeus by (1)H-NMR spectroscopy data]

The solution structure of insectotoxin 15A (35 residues) from scorpion Buthus eupeus was determined on the basis of 386 interproton distance restraints 12 hydrogen-bonding restraints and 113 dihedral angle restraints derived from 1H NMR experiments. A group of 20 structures was calculated with the distance geometry program DIANA followed by the restrained energy minimization with the program CHARMM. The atomic RMS distribution about the mean coordinate position is 0.64 +/- 0.11 A for the backbone atoms and 1.35 +/- 0.20 A for all atoms. The structure contains an alpha-helix (residues 10-20) and a three-stranded antiparallel beta-sheet (residues 2-5, 24-28 and 29-33). A pairing of the eight cysteine residues of insectotoxin 15A was established basing on NMR data. Three disulfide bridges (residues 2-19, 16-31 and 20-33) connect the alpha-helix with the beta-sheet, and the fourth one (5-26) joins beta-strands together. The spatial fold of secondary structure elements (the alpha-helix and the beta-sheet) of the insectotoxin 15A is very similar to those of the other short and long scorpion toxins in spite of a low (about 20%) sequence homology.     

The three-dimensional structure of guanine-specific ribonuclease F1 in solution determined by NMR spectroscopy and distance geometry.

Two-dimensional 1H-NMR studies have been performed on ribonuclease F1 (RNase F1), which contains 106 amino acid residues. Sequence-specific resonance assignments were accomplished for the backbone protons of 99 amino acid residues and for most of their side-chain protons. The three-dimensional structures were constructed on the basis of 820 interproton-distance restraints derived from NOE, 64 distance restraints for 32 hydrogen bonds and 33 phi torsion-angle restraints. A total of 40 structures were obtained by distance geometry and simulated-annealing calculations. The average root-mean-square deviation (residues 1-106) between the 40 converged structures and the mean structure obtained by averaging their coordinates was 0.116 +/- 0.018 nm for the backbone atoms and 0.182 +/- 0.015 nm for all atoms including the hydrogen atoms. RNase F1 was determined to be an alpha/beta-type protein. A well-defined structure constitutes the core region, which consists of a small N-terminal beta-sheet (beta 1, beta 2) and a central five-stranded beta-sheet (beta 3-beta 7) packed on a long helix. The structure of RNase F1 has been compared with that of RNase T1, which was determined by X-ray crystallography. Both belong to the same family of microbial ribonucleases. The polypeptide backbone fold of RNase F1 is basically identical to that of RNase T1. The conformation-dependent chemical shifts of the C alpha protons are well conserved between RNase F1 and RNase T1. The residues implicated in catalysis are all located on the central beta-sheet in a geometry similar to that of RNase T1.     

Determination of the three-dimensional solution structure of the antihypertensive and antiviral protein BDS-I from the sea anemone Anemonia sulcata: a study using nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing

The three-dimensional solution structure of the antihypertensive and antiviral protein BDS-I from the sea anemone Anemonia sulcata has been determined on the basis of 489 interproton and 24 hydrogen-bonding distance restraints supplemented by 23 phi backbone and 21 chi 1 side-chain torsion angle restraints derived from nuclear magnetic resonance (NMR) measurements. A total of 42 structures is calculated by a hybrid metric matrix distance geometry-dynamical simulated annealing approach. Both the backbone and side-chain atom positions are well defined. The average atomic rms difference between the 42 individual SA structures and the mean structure obtained by averaging their coordinates is 0.67 +/- 0.12 A for the backbone atoms and 0.90 +/- 0.17 A for all atoms. The core of the protein is formed by a triple-stranded antiparallel beta-sheet composed of residues 14-16 (strand 1), 30-34 (strand 2), and 37-41 (strand 3) with an additional mini-antiparallel beta-sheet at the N-terminus (residues 6-9). The first and second strands of the triple-stranded antiparallel beta-sheet are connected by a long exposed loop (residues 17-30). A number of side-chain interactions are discussed in light of the structure.     

Solution structure of the catalytic domain of gammadelta resolvase. Implications for the mechanism of catalysis

The site-specific DNA recombinase, gammadelta resolvase, from Escherichia coli catalyzes recombination of res site-containing plasmid DNA to two catenated circular DNA products. The catalytic domain (residues 1-105), lacking a C-terminal dimerization interface, has been constructed and the NMR solution structure of the monomer determined. The RMSD of the NMR conformers for residues 2-92 excluding residues 37-45 and 64-73 is 0.41 A for backbone atoms and 0.88 A for all heavy atoms. The NMR solution structure of the monomeric catalytic domain (residues 1-105) was found to be formed by a four-stranded parallel beta-sheet surrounded by three helices. The catalytic domain (residues 1-105), deficient in the C-terminal dimerization domain, was monomeric at high salt concentration, but displayed unexpected dimerization at lower ionic strength. The unique solution dimerization interface at low ionic strength was mapped by NMR. With respect to previous crystal structures of the dimeric catalytic domain (residues 1-140), differences in the average conformation of active-site residues were found at loop 1 containing the catalytic S10 nucleophile, the beta1 strand containing R8, and at loop 3 containing D67, R68 and R71, which are required for catalysis. The active-site loops display high-frequency and conformational backbone dynamics and are less well defined than the secondary structures. In the solution structure, the D67 side-chain is proximal to the S10 side-chain making the D67 carboxylate group a candidate for activation of S10 through general base catalysis. Four conserved Arg residues can function in the activation of the phosphodiester for nucleophilic attack by the S10 hydroxyl group. A mechanism for covalent catalysis by this class of recombinases is proposed that may be related to dimer interface dissociation.     

Plasminogen activator inhibitor 1. Structure of the native serpin, comparison to its other conformers and implications for serpin inactivation

The crystal structure of a constitutively active multiple site mutant of plasminogen activator inhibitor 1 (PAI-1) was determined and refined at a resolution of 2.7 A.The present structure comprises a dimer of two crystallographically independent PAI-1 molecules that pack by association of the residues P6 to P3 of the reactive centre loop of one molecule (A) with the edge of the main beta-sheet A of the other molecule (B).Thus, the reactive centre loop is ordered for molecule A by crystal packing forces, while for molecule B it is unconstrained by crystal packing contacts and is disordered.The overall structure of active PAI-1 is similar to the structures of other active inhibitory serpins exhibiting as the major secondary structural feature a five-stranded beta-sheet A and an intact proteinase-binding loop protruding from the one end of the elongated molecule. No preinsertion of the reactive centre loop is observed in this structure.A comparison of the present structure with the previously determined crystal structures of PAI-1 in its alternative conformations reveals that, upon cleavage of an intact form of PAI-1 or formation of latent PAI-1, the well-characterised rearrangements of the serpin secondary structural elements are accompanied by dramatic and partly unexpected conformational changes of helical and loop structures proximal to beta-sheet A.The present structure explains the stabilising effects of the mutated residues, reveals the structural cause for the observed spectroscopic differences between active and latent PAI-1, and provides new insights into possible mechanisms of stabilisation by its natural binding partner, vitronectin.     

Porcine pancreatic phospholipase A2: sequence-specific 1H and 15N NMR assignments and secondary structure

The solution structure of porcine pancreatic phospholipase A2 (124 residues, 14 kDa) has been studied by two-dimensional homonuclear 1H and two- and three-dimensional heteronuclear 15N-1H nuclear magnetic resonance spectroscopy. Backbone assignments were made for 117 of the 124 amino acids. Short-range nuclear Overhauser effect (NOE) data show three alpha-helices from residues 1-13, 40-58, and 90-109, an antiparallel beta-sheet for residues 74-85, and a small antiparallel beta-sheet between residues 25-26 and 115-116. A 15N-1H heteronuclear multiple-quantum correlation experiment was used to monitor amide proton exchange over a period of 22 h. In total, 61 amide protons showed slow or intermediate exchange, 46 of which are located in the three large helices. Helix 90-109 was found to be considerably more stable than the other helices. For the beta-sheets, four hydrogen bonds could be identified. The secondary structure of porcine PLA in solution, as deduced from NMR, is basically the same as the structure of porcine PLA in the crystalline state. Differences were found in the following regions, however. Residues 1-6 in the first alpha-helix are less structured in solution than in the crystal structure. Whereas in the crystal structure residues 24-29 are involved both in a beta-sheet with residues 115-117 and in a hairpin turn, the expected hydrogen bonds between residues 24-117 and 25-29 do not show slow exchange behavior. This and the absence of several expected NOEs imply that this region has a less well defined structure in solution. Finally, the hydrogen bond between residues 78-81, which is part of a beta-sheet, does not show slow exchange behavior.     

The effects of pressure and cholesterol on rotational motions of perylene in lipid bilayers

Using steady-state fluorescence polarization measurements, an isothermal pressure-induced phase transition was observed in dimyristoyl-L-alpha-phosphatidylcholine multilamellar vesicles containing perylene. The temperature-to-pressure equivalence, dT/dP, estimated from the phase transition pressure, P1/2, is about 22 K X kbar-1, which is comparable to values determined from diphenylhexatriene polarization (Chong, P.L.-G. and Weber, G. (1983) Biochemistry 22, 5544-5550). In addition, we have employed a new method, introduced in this paper, to calculate the rate of in-plane rotation (Rip) and the rate of out-of-plane rotation (Rop) of perylene in lipid bilayers. The effects of pressure and cholesterol on the rotational rates of perylene in two lipid bilayer systems have been examined. They are 1-palmitoyl-2-oleoyl-L-alpha-phosphatidylcholine (POPC) multilamellar vesicles (MLV) and 50 mol% cholesterol in POPC (MLV). Rop is smaller than Rip due to the fact that the out-of-plane rotation requires a larger volume change than the in-plane rotation. Cholesterol seems not to affect Rop significantly, but pressure causes a decrease in Rop by about a factor of three. In contrast, the effects of pressure and cholesterol on Rip are less straightforward. At 1 atm cholesterol increases Rip by a factor of about two. Similarly, in the absence of cholesterol 1.5 kbar pressure essentially triples Rip. However, if both cholesterol is added and pressure is applied, Rip decreases sharply. The possible interactions between cholesterol and perylene are discussed.     

Long-term molecular dynamics simulation of copper plastocyanin in water

A long molecular dynamics simulation (1.1 ns) of fully hydrated plastocyanin has been performed and analysed to relate protein dynamics to structural elements and functional properties. The solvated structure is described in detail by the analysis of H-bond network. During all the simulation, the crystal H-bond network is maintained in the beta-sheet regions, while several H-bonds are broken or formed on the external surface of the protein. To evaluate whether such changes could be due to conformational rearrangements or to solvent competition, we have examined the average number of H-bonds between protein atoms and water molecules, and the root mean square deviations from crystal structure as a function of protein residues. Protein mobility and flexibility have been examined by positional and dihedral angle rms fluctuations. Finally, cross-correlation maps have revealed the existence of correlated motions among residues connected by hydrogen bonds.     

Three-dimensional solution structure of beta cryptogein, a beta elicitin secreted by a phytopathogenic fungus Phytophthora cryptogea.

Cryptogein belongs to a new family of 10-kDa proteins called elicitins. Elicitins are necrotic and signaling proteins secreted by Phytophthora spp. responsible for the incompatible reaction and systemic hypersensitive-like necroses of diverse plant species leading to resistance against fungal or bacterial plant pathogens. The solution structure of beta cryptogein from Phytophthora cryptogea fungus was determined by using multidimensional heteronuclear nuclear magnetic resonance spectroscopy. A set of 18 structures was calculated using 1360 NOE-derived distance restraints and 40 dihedral angle restraints obtained from 3JHNH alpha couplings. The RMS deviation from the mean structure is 0.87 +/- 0.14 A for backbone atoms and 1.34 +/- 0.14 A for all the non-hydrogen atoms of residues 2 to 98. The structure of beta cryptogein reveals a novel protein fold, with five helices and a double-stranded beta-sheet facing an omega-loop. One edge of the beta-sheet and the adjacent face of the omega-loop form a hydrophobic cavity. This cavity made of highly conserved residues represents a plausible binding site. Residue 13, which has been identified from directed mutagenesis and natural sequence comparison studies as a key amino acid involved in the differential control of necrosis, is surface exposed and could contribute to the binding to a ligand or a receptor. The solution structure is close to the X-ray structure, with slight differences lightly due to the crystal packing.     

The three-dimensional structure of a glutathione S-transferase from the mu gene class. Structural analysis of the binary complex of isoenzyme 3-3 and glutathione at 2.2-A resolution.

The crystal structure of a mu class glutathione S-transferase (EC 2.5.1.18) from rat liver (isoenzyme 3-3) in complex with the physiological substrate glutathione (GSH) has been solved at 2.2-A resolution by multiple isomorphous replacement methods. The enzyme crystallized in the monoclinic space group C2 with unit cell dimensions of a = 87.98 A, b = 69.41 A, c = 81.34 A, and beta = 106.07 degrees. Oligonucleotide-directed site-specific mutagenesis played an important role in the solution of the structure in that the cysteine mutants C86S, C114S, and C173S were used to help locate the positions of mercuric ion sites in nonisomorphous derivatives with ethylmercuric phosphate and to align the sequence with the model derived from MIR phases. A complete model for the protein was not obtained until part of the solvent structure was interpreted. The dimer in the asymmetric unit refined to a crystallographic R = 0.171 for 19,298 data and I > or = 1.5 sigma (I). The final model consists of 4150 atoms, including all non-hydrogen atoms of 434 amino acid residues, two GSH molecules, and oxygen atoms of 474 water molecules. The dimeric enzyme is globular in shape with dimensions of 53 x 62 x 56 A. Crystal contacts are primarily responsible for conformational differences between the two subunits which are related by a noncrystallographic 2-fold axis. The structure of the type 3 subunit can be divided into two domains separated by a short linker, a smaller alpha/beta domain (domain I, residues 1-82), and a larger alpha domain (domain II, residues 90-217). Domain I contains four beta-strands which form a central mixed beta-sheet and three alpha-helices which are arranged in a beta alpha beta alpha beta beta alpha motif. Domain II is composed of five alpha-helices. Domain I can be considered the glutathione binding domain, while domain II seems to be primarily responsible for xenobiotic substrate binding. The active site is located in a deep (19-A) cavity which is composed of three relatively mobile structural elements: the long loop (residues 33-42) of domain I, the alpha 4/alpha 5 helix-turn-helix segment, and the C-terminal tail. GSH is bound at the active site in an extended conformation at one end of the beta-sheet of domain I with its backbone facing the cavity and the sulfur pointing toward the subunit to which it is bound.(ABSTRACT TRUNCATED AT 400 WORDS)     

Crystallographic dissection of the thermal motion of protein-sugar complex

The crystal structure of turkey egg lysozyme (TEL) complexed with di-N-acetylchitobiose (NAG2) was refined at 1.19 A resolution by the full-matrix least-squares method with anisotropic temperature factors, and its thermal motion was evaluated by the TLS method. The average ESDs of atomic parameters of nonhydrogen atoms were 0.030 A for coordinates and 0.025 A(2) for anisotropic temperature factors. The active site cleft of TEL binds the alpha-anomer of NAG2 in a nonproductive binding mode with its pyranose rings parallel to a beta-sheet. The TEL structure was compared with the re-refined 1.12 A structure of native TEL. The RMS difference for equivalent Calpha atoms was 0.103 A and a relatively large difference was observed in the region of residues 104-125 rather than in the beta-sheet region where NAG2 was bound. In contrast, the temperature factor of the beta-sheet region was significantly decreased by the NAG2 binding. The TLS model that describes the rigid body motion in translation, libration, and screw motion was adopted for the evaluation of the molecular motion of TEL and NAG2, and the TLS parameters were determined by the least-squares fit to U(ij). The contribution of the external motion of TEL was estimated to be 55.8% of the observed temperature factor for the native structure and 45.9% for the NAG2 complex. The internal motion of TEL represented with atomic thermal ellipsoids was very similar between the native and complex structures except the NAG2 binding region. In the structure of NAG2, the rigid body motion dominates the thermal motion. The center of rotation of NAG2, 4.45A far from the center of gravity, is on the nitrogen atom of the acetylamino group that is hydrogen bonded to the main-chain peptide groups of Asn49 and Ala107. The rigid body motion of NAG2 indicates that the acetylamino group is most strongly bound to the active site, and the recognition of this group is a crucial step of the substrate binding.     

Refined structure of porcine pepsinogen at 1.8 A resolution

The molecular structure of porcine pepsinogen at 1.8 A resolution has been determined by a combination of molecular replacement and multiple isomorphous phasing techniques. The resulting structure was refined by restrained-parameter least-squares methods. The final R factor [formula: see text] is 0.164 for 32,264 reflections with I greater than or equal to sigma (I) in the resolution range of 8.0 to 1.8 A. The model consists of 2785 protein atoms in 370 residues, a phosphoryl group on Ser68 and 238 ordered water molecules. The resulting molecular stereochemistry is consistent with a well-refined crystal structure with co-ordinate accuracy in the range of 0.10 to 0.15 A for the well-ordered regions of the molecule (B less than 15 A2). For the enzyme portion of the zymogen, the root-mean-square difference in C alpha atom co-ordinates with the refined porcine pepsin structure is 0.90 A (284 common atoms) and with the C alpha atoms of penicillopepsin it is 1.63 A (275 common atoms). The additional 44 N-terminal amino acids of the prosegment (Leu1p to Leu44p, using the letter p after the residue number to distinguish the residues of the prosegment) adopt a relatively compact structure consisting of a long beta-strand followed by two approximately orthogonal alpha-helices and a short 3(10)-helix. Intimate contacts, both electrostatic and hydrophobic interactions, are made with residues in the pepsin active site. The N-terminal beta-strand, Leu1p to Leu6p, forms part of the six-stranded beta-sheet common to the aspartic proteinases. In the zymogen the first 13 residues of pepsin, Ile1 to Glu13, adopt a completely different conformation from that of the mature enzyme. The C alpha atom of Ile1 must move approximately 44 A in going from its position in the inactive zymogen to its observed position in active pepsin. Electrostatic interactions of Lys36pN and hydrogen-bonding interactions of Tyr37pOH, and Tyr90H with the two catalytic aspartate groups, Asp32 and Asp215, prevent substrate access to the active site of the zymogen. We have made a detailed comparison of the mammalian pepsinogen fold with the fungal aspartic proteinase fold of penicillopepsin, used for the molecular replacement solution. A structurally derived alignment of the two sequences is presented.     

Solution structure of a sweet protein single-chain monellin determined by nuclear magnetic resonance and dynamical simulated annealing calculations

Single-chain monellin (SCM), which is an engineered 94-residue polypeptide, has proven to be as sweet as native two-chain monellin. SCM is more stable than the native monellin for both heat and acidic environments. Data from gel filtration HPLC and NMR indicate that the SCM exists as a monomer in aqueous solution. The solution structure of SCM has been determined by nuclear magnetic resonance (NMR) spectroscopy and dynamical simulated annealing calculations. A stable alpha-helix spanning residues Phe11-Ile26 and an antiparallel beta-sheet formed by residues 2-5, 36-38, 41-47, 54-64, 69-75, and 83-88 have been identified. The sheet was well defined by backbone-backbone NOEs, and the corresponding beta-strands were further confirmed by hydrogen bond networks based on amide hydrogen exchange data. Strands beta2 and beta3 are connected by a small bulge comprising residues Ile38-Cys41. A total of 993 distance and 56 dihedral angle restraints were used for simulated annealing calculations. The final simulated annealing structures (k) converged well with a root-mean-square deviation (rmsd) between backbone atoms of 0.49 A for secondary structural regions and 0.70 A for backbone atoms excluding two loop regions. The average restraint energy-minimized (REM) structure exhibited root-mean-square deviations of 1.19 A for backbone atoms and 0.85 A for backbone atoms excluding two loop regions with respect to 20 k structures. The solution structure of SCM revealed that the long alpha-helix was folded into the concave side of a six-stranded antiparallel beta-sheet. The side chains of Tyr63 and Asp66 which are common to all sweet peptides showed an opposite orientation relative to H1 helix, and they were all solvent-exposed. Residues at the proposed dimeric interface in the X-ray structure were observed to be mostly solvent-exposed and demonstrated high degrees of flexibility.     

Solution structure of PAFP-S: a new knottin-type antifungal peptide from the seeds of Phytolacca americana

The three-dimensional solution structure of PAFP-S, an antifungal peptide extracted from the seeds of Phytolacca americana, was determined using 1H NMR spectroscopy. This cationic peptide contains 38 amino acid residues. Its structure was determined from 302 distance restraints and 36 dihedral restraints derived from NOEs and coupling constants. The peptide has six cysteines involved in three disulfide bonds. The previously unassigned parings have now been determined from NMR data. The solution structure of PAFP-S is presented as a set of 20 structures using ab initio dynamic simulated annealing, with an average RMS deviation of 1.68 A for the backbone heavy atoms and 2.19 A for all heavy atoms, respectively. For the well-defined triple-stranded beta-sheet involving residues 8-10, 23-27, and 32-36, the corresponding values were 0.39 and 1.25 A. The global fold involves a cystine-knotted three-stranded antiparallel beta-sheet (residues 8-10, 23-27, 32-36), a flexible loop (residues 14-19), and four beta-reverse turns (residues 4-8, 11-14, 19-22, 28-32). This structure features all the characteristics of the knottin fold. It is the first structural model of an antifungal peptide that adopts a knottin-type structure. PAFP-S has an extended hydrophobic surface comprised of residues Tyr23, Phe25, Ile27, Tyr32, and Val34. The side chains of these residues are well-defined in the NMR structure. Several hydrophilic and positively charged residues (Arg9, Arg38, and Lys36) surround the hydrophobic surface, giving PAFP-S an amphiphilic character which would be the main structural basis of its biological function.