Solution structure of apamin determined by nuclear magnetic resonance and distance geometry

The solution structure of the bee venom neurotoxin apamin has been determined with a distance geometry program using distance constraints derived from NMR. Twenty embedded structures were generated and refined by using the program DSPACE. After error minimization using both conjugate gradient and dynamics algorithms, six structures had very low residual error. Comparisons of these show that the backbone of the peptide is quite well-defined with the largest rms difference between backbone atoms in these structures of 1.34 A. The side chains have far fewer constraints and show greater variability in their positions. The structure derived here is generally consistent with the qualitative model previously described, with most differences occurring in the loop between the beta-turn (residues 2-5) and the C-terminal alpha-helix (residues 9-17). Comparisons are made with previously derived models from NMR data and other methods.     

NMR-pseudoenergy approach to the solution structure of acyl carrier protein

A method for protein structure determination from two-dimensional NMR cross-relaxation data is presented and explored by using short amino acid segments from acyl carrier protein as a test case. The method is based on a molecular mechanics program and incorporates NMR distance constraints in the form of a pseudoenergy term that accurately reflects the distance-dependent precision of NMR cross-relaxation data. When it is used in an indiscriminant fashion, the method has a tendency to produce structures representing local energy minima near starting structures, rather than structures representing a global energy minimum. However, stepwise inclusion of energy terms, beginning with a function heavily weighted by backbone distance constraints, appears to simplify the potential energy surface to a point where convergence to a common backbone structure from a variety of starting structures is possible. In the case of the segment from residues 3 to 15 in acyl carrier protein, a nearly perfect alpha-helix is produced starting with a linear chain, an alpha-helical chain, or a chain having residues with alternating linear and alpha-helical backbone torsional angles. In the case of the segment from residues 26 to 36 a structure having a right-handed loop is produced.     

Three-dimensional structure of the neurotoxin ATX Ia from Anemonia sulcata in aqueous solution determined by nuclear magnetic resonance spectroscopy

With the aid of 1H nuclear magnetic resonance (NMR) spectroscopy, the three-dimensional structure in aqueous solution was determined for ATX Ia, which is a 46 residue polypeptide neurotoxin of the sea anemone Anemonia sulcata. The input for the structure calculations consisted of 263 distance constraints from nuclear Overhauser effects (NOE) and 76 vicinal coupling constants. For the structure calculation several new or ammended programs were used in a revised strategy consisting of five successive computational steps. First, the program HABAS was used for a complete search of all backbone and chi 1 conformations that are compatible with the intraresidual and sequential NMR constraints. Second, using the program DISMAN, we extended this approach to pentapeptides by extensive sampling of all conformations that are consistent with the local and medium-range NMR constraints. Both steps resulted in the definition of additional dihedral angle constraints and in stereospecific assignments for a number of beta-methylene groups. In the next two steps DISMAN was used to obtain a group of eight conformers that contain no significant residual violations of the NMR constraints or van der Waals contacts. Finally, these structures were subjected to restrained energy refinement with a modified version of the molecular mechanics module of AMBER, which in addition to the energy force field includes potentials for the NOE distance constraints and the dihedral angle constraints. The average of the pairwise minimal RMS distances between the resulting refined conformers calculated for the well defined molecular core, which contains the backbone atoms of 35 residues and 20 interior side chains, is 1.5 +/- 0.3 A. This core is formed by a four-stranded beta-sheet connected by two well-defined loops, and there is an additional flexible loop consisting of the eleven residues 8-18. The core of the protein is stabilized by three disulfide bridges, which are surrounded by hydrophobic residues and shielded on one side by hydrophilic residues.     

The NMR solution structure of the pheromone Er-11 from the ciliated protozoan Euplotes raikovi.

The NMR solution structure of the pheromone Er-11, a 39-residue protein from the ciliated protozoan Euplotes raikovi, was calculated with the distance geometry program DIANA from 449 NOE upper distance constraints and 97 dihedral angle constraints, and the program OPAL was employed for structure refinement by molecular mechanics energy minimization in a water bath. For a group of 20 conformers used to characterize the solution structure, the average of the pairwise RMS deviations from the mean structure calculated for the backbone heavy atoms N, C alpha, and C' of residues 2-38 was 0.30 A. The molecular architecture is dominated by an up-down-up bundle of three short helices with residues 2-9, 12-19, and 22-32, which is closely similar to the previously determined structures of the homologous pheromones Er-1, Er-2, and Er-10. This finding provides structural evidence for the capability shown by these pheromones to compete with each other in binding reactions to their cell-surface receptors.     

Solution structure of murine epidermal growth factor determined by NMR spectroscopy and refined by energy minimization with restraints.

The solution structure of murine epidermal growth factor (mEGF) at pH 3.1 and a temperature of 28 degrees C has been determined from NMR data, using distance geometry calculations and restrained energy minimization. The structure determination is based on 730 conformational constraints derived from NMR data, including 644 NOE-derived upper bound distance constraints, constraints on the ranges of 32 dihedral angles based on measurements of vicinal coupling constants, and 54 upper and lower bound constraints associated with nine hydrogen bonds and the three disulfide bonds. The distance geometry interpretation of the NMR data is based on previously published sequence-specific 1H resonance assignments [Montelione et al. (1988) Biochemistry 27, 2235-2243], supplemented here with individual assignments for some side-chain amide, methylene, and isopropyl methyl protons. The molecular architecture of mEGF is the same as that described previously [Montelione et al. (1987) Proc. Natl. Acad. Sci. U.S.A. 84, 5226-5230], but the structure is overall more precisely determined by a more extensive set of NMR constraints. Analysis of proton NMR line widths, amide proton exchange rates, and side-chain 3J(H alpha-H beta) coupling constants provides evidence for internal motion in several regions of the mEGF molecule. Because mEGF is one member of a large family of homologous growth factors and protein domains for which X-ray crystal structures are not yet available, the atomic coordinates resulting from the present structure refinement (which we have deposited in the Brookhaven Protein Data Bank) are important data for understanding the structures of EGF-like proteins and for further detailed comparisons of these structures with mEGF.     

Solution structures of alpha-conotoxin G1 determined by two-dimensional NMR spectroscopy

Two-dimensional NMR data have been used to generate solution structures of alpha-conotoxin G1, a potent peptide antagonist of the acetylcholine receptor. Structural information was obtained in the form of proton-proton internuclear distance constraints, and initial structures were produced with a distance geometry algorithm. Energetically more favorable structures were generated by using the distance geometry structures as input for a constrained energy minimization program. The results of both of these calculations indicate that the overall backbone conformation of the molecule is well-defined by the NMR data whereas the side-chain conformations are generally less well-defined. The main structural features derived from the NMR data were the presence of tight turns centered on residues Pro5 and Arg9. The solution structures are compared with previous proposed models of conotoxin G1, and the NMR data are interpreted in conjunction with chemical modification studies and structural properties of other antagonists of the acetylcholine receptor to gain insight into structure-activity relationships in these peptide toxins.     

Improved strategies for the determination of protein structures from NMR data: the solution structure of acyl carrier protein

The hybrid method that combines the early stages of a distance geometry program with simulated annealing in the presence of NMR constraints was optimized to obtain structures fully consistent with the observed NMR data. This was achieved by using more restrictive bounds of the NOE constraints than those usually used in the literature and by grouping the NOEs into classes dependent on the quality of the experimental NOE data. The 'floating' stereospecific assignment introduced at the simulated annealing stage of the calculations further improved the definition of the local conformation. An improved sampling and convergence property of the hybrid method was obtained by means of fitting the substructure obtained from the distance geometry program to different conformations. Compared to the standard hybrid methods, this procedure gave superior structures for a 77 amino acid protein, acyl carrier protein from Escherichia coli.     

Three-dimensional solution structure of the HIV-1 protease complexed with DMP323, a novel cyclic urea-type inhibitor, determined by nuclear magnetic resonance spectroscopy.

The three-dimensional solution structure of the HIV-1 protease homodimer, MW 22.2 kDa, complexed to a potent, cyclic urea-based inhibitor, DMP323, is reported. This is the first solution structure of an HIV protease/inhibitor complex that has been elucidated. Multidimensional heteronuclear NMR spectra were used to assemble more than 4,200 distance and angle constraints. Using the constraints, together with a hybrid distance geometry/simulated annealing protocol, an ensemble of 28 NMR structures was calculated having no distance or angle violations greater than 0.3 A or 5 degrees, respectively. Neglecting residues in disordered loops, the RMS deviation (RMSD) for backbone atoms in the family of structures was 0.60 A relative to the average structure. The individual NMR structures had excellent covalent geometry and stereochemistry, as did the restrained minimized average structure. The latter structure is similar to the 1.8-A X-ray structure of the protease/DMP323 complex (Chang CH et al., 1995, Protein Science, submitted); the pairwise backbone RMSD calculated for the two structures is 1.22 A. As expected, the mismatch between the structures is greatest in the loops that are disordered and/or flexible. The flexibility of residues 37-42 and 50-51 may be important in facilitating substrate binding and product release, because these residues make up the respective hinges and tips of the protease flaps. Flexibility of residues 4-8 may play a role in protease regulation by facilitating autolysis.     

Three-dimensional structure of acyl carrier protein in solution determined by nuclear magnetic resonance and the combined use of dynamical simulated annealing and distance geometry

The solution conformation of acyl carrier protein from Escherichia coli (77 residues) has been determined on the basis of 423 interproton-distance restraints and 32 hydrogen-bonding restraints derived from NMR measurements. A total of nine structures were computed using a hybrid approach combining metric matrix distance geometry and dynamic simulated annealing. The polypeptide fold is well defined with an average backbone atomic root-mean-square difference of 0.20 +/- 0.03 nm between the final nine converged structures and the mean structure obtained by averaging their coordinates. The principal structural motif is composed of three helices: 1 (residues 3-12), 2 (residues 37-47) and 4 (residues 65-75) which line a hydrophobic cavity. Helices 2 and 4 are approximately parallel to each other and anti-parallel at an angle of approximately equal to 150 degrees to helix 1. The smaller helix 3 (residues 56-63) is at an angle of approximately equal to 100 degrees to helix 4.     

Three-dimensional structure of ectatomin from Ectatomma tuberculatum ant venom

Summary Two-dimensional ¹H NMR techniques were used to determine the spatial structure of ectatomin, a toxin from the venom of the ant Ectatomma tuberculatum. Nearly complete proton resonance assignments for two chains of ectatomin (37 and 34 amino acid residues, respectively) were obtained using 2D TOCSY, DQF-COSY and NOESY experiments. The cross-peak volumes in NOESY spectra were used to define the local structure of the protein and generate accurate proton-proton distance constraints employing the MARDIGRAS program. Disulfide bonds were located by analyzing the global fold of ectatomin, calculated with the distance geometry program DIANA. These data, combined with data on the rate of exchange of amide protons with deuterium, were used to obtain a final set of 20 structures by DIANA. These structures were refined by unrestrained energy minimization using the CHARMm program. The resulting rms deviations over 20 structures (excluding the mobile N- and C-termini of each chain) are 0.75 ? for backbone heavy atoms, and 1.25 ? for all heavy atoms. The conformations of the two chains are similar. Each chain consists of two α-helices and a hinge region of four residues; this forms a hairpin structure which is stabilized by disulfide bridges. The hinge regions of the two chains are connected together by a third disulfide bridge. Thus, ectatomin forms a four-α-helical bundle structure.     

Design, synthesis and solution structure of a renin inhibitor. Structural constraints from NOE, and homonuclear and heteronuclear coupling constants combined with distance geometry calculations.

A macrocyclic renin inhibitor was designed using molecular modeling and a model of human renin. The synthesized molecular displayed poor binding affinity. To investigate the reasons for the observed inactivity, the structure of the compound has been studied by NMR spectroscopy and distance geometry. Structural constraints for distance geometry calculations were derived from nuclear Overhauser effects and homonuclear and heteronuclear three bond coupling constants. Homonuclear coupling constants were measured directly from the resolution-enhanced proton spectra and heteronuclear coupling constants were measured from the natural abundance 15N- and 13C-edited TOCSY experiments. One phi angle was determined uniquely by this method and two were reduced to two possible values each. By using a statistical analysis of 400 structures generated with distance geometry, two families of structures were found to be consistent with the NMR data. The solution structures so derived were different from the originally designed structure, including an internal hydrogen bond. This provides a possible explanation for the lack of effectiveness of this compound.     

Determining the three-dimensional fold of a protein from approximate constraints

We propose a new approach for calculating the three-dimensional (3D) structure of a protein from distance and dihedral angle constraints derived from experimental data. We suggest that such constraints can be obtained from experiments such as tritium planigraphy, chemical or enzymatic cleavage of the polypeptide chain, paramagnetic perturbation of nuclear magnetic resonance (NMR) spectra, measurement of hydrogen-exchange rates, mutational studies, mass spectrometry, and electron paramagnetic resonance. These can be supplemented with constraints from theoretical prediction of secondary structures and of buried/exposed residues. We report here distance geometry calculations to generate the structures of a test protein Staphylococcal nuclease (STN), and the HIV-1 rev protein (REV) of unknown structure. From the available 3D atomic coordinates of STN, we set up simulated data sets consisting of varying number and quality of constraints, and used our group's Self Correcting Distance Geometry (SECODG) program DIAMOD to generate structures. We could generate the correct tertiary fold from qualitative (approximate) as well as precise distance constraints. The root mean square deviations of backbone atoms from the native structure were in the range of 2.0 A to 8.3 A, depending on the number of constraints used. We could also generate the correct fold starting from a subset of atoms that are on the surface and those that are buried. When we used data sets containing a small fraction of incorrect distance constraints, the SECODG technique was able to detect and correct them. In the case of REV, we used a combination of constraints obtained from mutagenic data and structure predictions. DIAMOD generated helix-loop-helix models, which, after four self-correcting cycles, populated one family exclusively. The features of the energy-minimized model are consistent with the available data on REV-RNA interaction. Our method could thus be an attractive alternative for calculating protein 3D structures, especially in cases where the traditional methods of X-ray crystallography and multidimensional NMR spectroscopy have been unsuccessful.     

A refined three-dimensional solution structure of a carboxy terminal fragment of apolipoprotein CII

The three-dimensional structure of a synthetic fragment of human apolipoprotein CII (apo-CII) in 35%, 1,1,1,3,3,3-hexafluoro-2-propanol (HFP) has been determined on the basis of distance and intensity constraints derived from two-dimensional proton nuclear magnetic resonance measurements. The NOE crosspeak build-up rates were converted to distance constraints which were used in the distance geometry program DIANA. A set of one hundred structures were generated and of these ten structures were used in molecular dynamics simulations using the program XPLOR. This program enabled a direct minimization between the difference of the two-dimensional NOE intensities and those calculated from the full relaxation matrix. In this way spin diffusion is fully taken into account, which can be seen from the considerable improvement of the R-factor after the relaxation matrix refinement. These calculations show that this fragment, which corresponds to the carboxy terminal 30 amino acids of intact apo-CII and which retains its ability to activate lipoprotein lipase, is essentially flexible, but has three defined secondary structural elements. The most significant one is an -helix between residues 67 and 74. The following three residues adopt a turn-like structure. Another turn of -helix is seen between residues 56 and 59. The effect of the solvent system on the secondary structure was studied by circular dichroism spectroscopy. The results show that the mixed aqueous 35% HFP solvent induces secondary structure of a very similar nature to the one induced by sodium dodecyl sulphate.Abbreviations Apo-CII Apolipoprotein CII - CD Circular Dichroism - DOPC 1,2-dioleoyl-sn-glycero-3-phosphocholine - DOPG 1,2-dioleoyl-sn-glycero-3-phosphoglycerol - HAc Acetic Acid - HFP 1,1,1,3,3,3-hexafluoro-2-propanol - ISPA Isolated Spin Pair Approximation - NMR Nuclear Magnetic Resonance - NOE Nuclear Overhauser Enhancement - NOESY Nuclear Overhauser Enhancement Spectroscopy - RMSD Root Mean Square Deviation - SDS Sodium Dodecyl Sulfate     

Solution structure of Fe(II) cytochrome c551 from Pseudomonas aeruginosa as determined by two-dimensional 1H NMR.

The solution structure of Fe(II) cytochrome c551 from Pseudomonas aeruginosa based on 2D 1H NMR data is reported. Two sets of structure calculations were completed with a combination of simulated annealing and distance geometry calculations: one set of 20 structures included the heme-peptide covalent linkages, and one set of 10 structures excluded them. The main-chain atoms were well constrained within the two structural ensembles (1.30 and 1.35 A average RMSD, respectively) except for two regions spanning residues 30-40 and 60-70. The results were essentially the same when global fold comparisons were made between the ensembles with an average RMSD of 1.33 A. In total, 556 constraints were used, including 479 NOEs, 53 volume constraints, and 24 other distances. This report represents the first solution structure determination of a heme protein by 2D 1H NMR and should provide a basis for the application of these techniques to other proteins containing large prosthetic groups or cofactors.     

Conformational analysis of human calcitonin in solution.

The solution conformation of human calcitonin in a mixture of 60% water and 40% trifluoroethanol has been determined by the combined use of 1H NMR spectroscopy and distance geometry calculations with a distributed computing technique. 1H NMR spectroscopy provided 195 distance constraints and 13 hydrogen bond constraints. The 20 best converged structures exhibit atomic rmsd of 0.43 A for the backbone atoms from the averaged coordinate position in the region of Asn3-Phe22. The conformation is characterized by a nearly amphiphilic alpha-helix domain that extends from Leu4 in the cyclic region to His20. There are no significant differences observed among the overall structures of a series of calcitonins obtained from ultimobranchial bodies, including those that possess 20- to 50-fold greater activity. Three aromatic amino acid residues, Tyr12, Phe16 and Phe19, form a hydrophobic surface of human calcitonin. Bulky side chains on the surface could interfere with the ligand-receptor interaction thereby causing its low activity, relative to those of other species.     

Solution structure of amphibian tachykinin Uperolein bound to DPC micelles

Uperolein, a physalaemin-like endecapeptide, has been shown to be selective for Neurokinin 1 receptor. As a first step towards understanding the structure-activity relationship, we report the membrane-induced structure of Uperolein with the aid of circular dichroism and 2D (1)H NMR spectroscopy. Sequence-specific resonance assignments of protons have been made using correlation spectroscopy (TOCSY, DQF-COSY) and NOESY spectroscopy. The interproton distance constraints and dihedral angle constraints have been utilized to generate a family of structures using torsion angle molecular dynamics within program DYANA. The conformational range of the peptide revealed by NMR and CD studies has been analysed in terms of characteristic secondary features. Analysis of NMR data indicates that the global fold of Uperolein can be explained in terms of equilibrium between 3(10)-helix and alpha-helix from residues 5 to 11. An extended highly flexible N-terminus displays some degree of order and a possible turn structure. A comparison between the structures of Uperolein and Substance P, a prototype and endogenous Neurokinin 1 receptor agonist, indicates several common features in the distribution of hydrophobic and hydrophilic residues. Both the peptides show an amphiphilic character towards the middle region. The similarities suggest that the molecules interact with the receptor in an analogous manner.     

High-resolution solution structure of reduced French bean plastocyanin and comparison with the crystal structure of poplar plastocyanin

The three-dimensional solution structure of reduced (CuI) plastocyanin from French bean leaves has been determined by distance geometry and restrained molecular dynamics methods using constraints obtained from 1H n.m.r. (nuclear magnetic resonance) spectroscopy. A total of 1244 experimental constraints were used, including 1120 distance constraints, 103 dihedral angle constraints and 21 hydrogen bond constraints. Stereospecific assignments were made for 26 methylene groups and the methyls of 11 valines. Additional constraints on copper co-ordination were included in the restrained dynamics calculations. The structures are well defined with average atomic root-mean-square deviations from the mean of 0.45 A for all backbone heavy atoms and 1.08 A for side-chain heavy atoms. French bean plastocyanin adopts a beta-sandwich structure in solution that is similar to the X-ray structure of reduced poplar plastocyanin; the average atomic root-mean-square difference between 16 n.m.r. structures and the X-ray structure is 0.76 A for all backbone heavy atoms. The conformations of the side-chains that constitute the hydrophobic core of French bean plastocyanin are very well defined. Of 47 conserved residues that populate a single chi 1 angle in solution, 43 have the same rotamer in the X-ray structure. Many surface side-chains adopt highly preferred conformations in solution, although the 3J alpha beta coupling constants often indicate some degree of conformational averaging. Some surface side-chains are disordered in both the solution and crystal structures of plastocyanin. There is a striking correlation between measures of side-chain disorder in solution and side-chain temperature factors in the X-ray structure. Side-chains that form a distinctive acidic surface region, believed to be important in binding other electron transfer proteins, appear to be disordered. Fifty backbone amide protons form hydrogen bonds to carbonyls in more than 60% of the n.m.r. structures; 45 of these amide protons exchange slowly with solvent deuterons. Ten hydrogen bonds are formed between side-chain and backbone atoms, eight of which are correlated with decreased proton exchange. Of the 60 hydrogen bonds formed in French bean plastocyanin, 56 occur in the X-ray structure of the poplar protein; two of the missing hydrogen bonds are absent as a result of mutations. It appears that molecular dynamics refinement of highly constrained n.m.r. structures allows accurate prediction of the pattern of hydrogen bonding.     

Solution structure of the calcium channel antagonist omega-conotoxin GVIA.

The three-dimensional solution structure is reported for omega-conotoxin GVIA, which is a potent inhibitor of presynaptic calcium channels in vertebrate neuromuscular junctions. Structures were generated by a hybrid distance geometry and restrained molecular dynamics approach using interproton distance, torsion angle, and hydrogen-bonding constraints derived from 1H NMR data. Conformations of GVIA with low constraint violations converged to a common peptide fold. The secondary structure in the peptide is an antiparallel triple-stranded beta-sheet containing a beta-hairpin and three tight turns. The NMR data are consistent with the region of the peptide from residues S9 to C16 being more dynamic than the rest of the peptide. The peptide has an amphiphilic structure with a positively charged hydrophilic side and an opposite side that contains a small hydrophobic region. Residues that are thought to be important in binding and function are located on the hydrophilic face of the peptide.     

Solution conformation of thymosin beta 4: a nuclear magnetic resonance and simulated annealing study

The conformation of the polypeptide thymosin beta 4 in solutions of 60% (v/v) trifluoroethanol-d3 and 50% (v/v) hexafluoroisopropyl-d2 alcohol in water is investigated by nuclear magnetic resonance (NMR) spectroscopy. Under these conditions thymosin beta 4 adopts an ordered structure. By use of a combination of two-dimensional NMR techniques, the 1H NMR spectrum of thymosin beta 4 is assigned. A set of 180 approximate interproton distance constraints is derived from nuclear Overhauser enhancement (NOE) measurements. These, together with 33 phi constraints obtained for JNH alpha coupling data and the 23 psi dihedral angles identified on the basis of the pattern of short-range NOEs, form the basis of a three-dimensional structure determination by dynamical simulated annealing. The calculations are carried out starting from three initial structures, an alpha-helix, an extended beta-strand, and a mixed alpha/beta structure. Ten independent structures are computed from each starting structure by using different random number seeds for the assignments of the initial velocities. All 30 calculated structures satisfy the experimental constraints, display very small deviations from idealized covalent geometry, and possess good nonbonded contacts. Analysis of the 30 converged structures indicates that there are two helical regions extending from residues 4-16 and from residues 30-40, which are well defined both in terms of atomic root mean square differences and backbone torsion angles. For the two helical regions individually the average backbone rms difference between all pairs of structures is approximately 2 A. The two helices exhibit typical amino acid preferences for specific locations at the ends of helices.(ABSTRACT TRUNCATED AT 250 WORDS)     

Three-dimensional solution structure of the B domain of staphylococcal protein A: comparisons of the solution and crystal structures.

The three-dimensional solution structure of the recombinant B domain (FB) of staphylococcal protein A, which specifically binds to the Fc portion of immunoglobulin G, was determined by NMR spectroscopy and hybrid distance geometry-dynamical simulated annealing calculations. On the basis of 692 experimental constraints including 587 distance constraints obtained from the nuclear Overhauser effect (NOE), 57 torsion angle (phi, chi 1) constraints, and 48 constraints associated with 24 hydrogen bonds, a total of 10 converged structures of FB were obtained. The atomic root mean square difference among the 10 converged structures is 0.52 +/- 0.10 A for the backbone atoms and 0.98 +/- 0.08 A for all heavy atoms (excluding the N-terminal segment from Thr1 to Glu9 and the C-terminal segment from Gln56 to Ala60, which are partially disordered). FB is composed of a bundle of three alpha-helices, i.e., helix I (Gln10-His19), helix II (Glu25-Asp37), and helix III (Ser42-Ala55). Helix II and helix III are antiparallel to each other, whereas the long axis of helix I is tilted at an angle of about 30 degrees with respect to those of helix II and helix III. Most of the hydrophobic residues of FB are buried in the interior of the bundle of the three helices. It is suggested that the buried hydrophobic residues form a hydrophobic core, contributing to the stability of FB.(ABSTRACT TRUNCATED AT 250 WORDS)