Determination of the three-dimensional solution structure of the C-terminal domain of cellobiohydrolase I from Trichoderma reesei. A study using nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing

The solution structure of a synthetic 36-residue polypeptide comprising the C-terminal cellulose binding domain of cellobiohydrolase I (CT-CBH I) from Trichoderma reesei was investigated by nuclear magnetic resonance (NMR) spectroscopy. The 1H NMR spectrum was completely assigned in a sequential manner by two-dimensional NMR techniques. A large number of stereospecific assignments for beta-methylene protons, as well as ranges for the phi, psi, and chi 1 torsion angles, were obtained on the basis of sequential and intraresidue nuclear Overhauser enhancement (NOE) and coupling constant data in combination with a conformational data base search. The structure calculations were carried out in an iterative manner by using the hybrid distance geometry-dynamical simulated annealing method. This involved computing a series of initial structures from a subset of the experimental data in order to resolve ambiguities in the assignments of some NOE cross-peaks arising from chemical shift degeneracy. Additionally, this permitted us to extend the stereospecific assignments to the alpha-methylene protons of glycine using information on phi torsion angles derived from the initial structure calculations. The final experimental data set consisted of 554 interproton distance restraints, 24 restraints for 12 hydrogen bonds, and 33 phi, 24 psi, and 25 chi 1 torsion angle restraints. CT-CBH I has two disulfide bridges whose pairing was previously unknown. Analysis of structures calculated with all three possible combinations of disulfide bonds, as well as without disulfide bonds, indicated that the correct disulfide bridge pairing was 8-25 and 19-35. Forty-one structures were computed with the 8-25 and 19-35 disulfide bridges, and the average atomic rms difference between the individual structures and the mean structure obtained by averaging their coordinates was 0.33 +/- 0.04 A for the backbone atoms and 0.52 +/- 0.06 A for all atoms. The protein has a wedgelike shape with an amphiphilic character, one face being predominantly hydrophilic and the other mainly hydrophobic. The principal element of secondary structure is made up of an irregular triple-stranded antiparallel beta-sheet composed of residues 5-9 (beta 1), 24-28 (beta 2), and 33-36 (beta 3) in which strand beta 3 is hydrogen bonded to the other two strands.(ABSTRACT TRUNCATED AT 400 WORDS)     

Solution structure of a lipid transfer protein extracted from rice seeds. Comparison with homologous proteins.

Nuclear magnetic resonance (NMR) spectroscopy was used to determine the three dimensional structure of rice nonspecific lipid transfer protein (ns-LTP), a 91 amino acid residue protein belonging to the broad family of plant ns-LTP. Sequence specific assignment was obtained for all but three HN backbone 1H resonances and for more than 95% of the 1H side-chain resonances using a combination of 1H 2D NOESY; TOCSY and COSY experiments at 293 K. The structure was calculated on the basis of four disulfide bridge restraints, 1259 distance constraints derived from 1H-1H Overhauser effects, 72 phi angle restraints and 32 hydrogen-bond restraints. The final solution structure involves four helices (H1: Cys3-Arg18, H2: Ala25-Ala37, H3: Thr41-Ala54 and H4: Ala66-Cys73) followed by a long C-terminal tail (T) with no observable regular structure. N-capping residues (Thr2, Ser24, Thr40), whose side-chain oxygen atoms are involved in hydrogen bonds with i + 3 amide proton additionally stabilize the N termini of the first three helices. The fourth helix involving Pro residues display a mixture of alpha and 3(10) conformation. The rms deviation of 14 final structures with respect to the average structure is 1.14 +/- 0.16 A for all heavy atoms (C, N, O and S) and 0.72 +/- 0.01 A for the backbone atoms. The global fold of rice ns-LTP is close to the previously published structures of wheat, barley and maize ns-LTPs exhibiting nearly identical pattern of the numerous sequence specific interactions. As reported previously for different four-helix topology proteins, hydrophobic, hydrogen bonding and electrostatic mechanisms of fold stabilization were found for the rice ns-LTP. The sequential alignment of 36 ns-LTP primary structures strongly suggests that there is a uniform pattern of specific long-range interactions (in terms of sequence), which stabilize the fold of all plant ns-LTPs.     

Reduced representation model of protein structure prediction: statistical potential and genetic algorithms.

A reduced representation model, which has been described in previous reports, was used to predict the folded structures of proteins from their primary sequences and random starting conformations. The molecular structure of each protein has been reduced to its backbone atoms (with ideal fixed bond lengths and valence angles) and each side chain approximated by a single virtual united-atom. The coordinate variables were the backbone dihedral angles phi and psi. A statistical potential function, which included local and nonlocal interactions and was computed from known protein structures, was used in the structure minimization. A novel approach, employing the concepts of genetic algorithms, has been developed to simultaneously optimize a population of conformations. With the information of primary sequence and the radius of gyration of the crystal structure only, and starting from randomly generated initial conformations, I have been able to fold melittin, a protein of 26 residues, with high computational convergence. The computed structures have a root mean square error of 1.66 A (distance matrix error = 0.99 A) on average to the crystal structure. Similar results for avian pancreatic polypeptide inhibitor, a protein of 36 residues, are obtained. Application of the method to apamin, an 18-residue polypeptide with two disulfide bonds, shows that it folds apamin to native-like conformations with the correct disulfide bonds formed.     

Method of determining protein conformations by the two-dimensional nuclear Overhauser enhancement spectroscopy data

A method is suggested to determine the most probable values of the angles phi, psi of the protein backbone by the data on the availability and absence of d connectivities in the two-dimensional nuclear Overhauser enhancement spectra. In view of this, the dependences of the proton-proton distances in dipeptide units of L-amino acid residues on the dihedral angles phi, psi, chi1 are considered and the conformational states of amino acid residues of the proteins with the known spatial structure are analysed statistically. The potentialities of the method are assessed with the aid of model spectral nuclear magnetic resonance (NMR) parameters obtained from the X-ray data for the bovine pancreatic trypsin inhibitor and avian pancreatic polypeptide. It is shown that the developed procedure of structural interpretation of the NMR data allows one to correctly reproduce the local conformation of the protein backbone. The obtained backbone conformation may serve as a starting point to build and refine molecular three-dimensional structure.     

Combined use of homo- and heteronuclear coupling constants as restraints in molecular dynamics simulations.

A penalty function for scalar coupling constants has been applied in molecular dynamics simulations as an experimental constraint. The function is based on the difference between the coupling constant calculated from the dihedral angle and the experimentally measured coupling constant. The method is illustrated on a model cyclic pentapeptide for which 3JHN-H alpha and 3JHN-C beta, both about the phi backbone dihedral angle, have been measured. The function is efficient in producing structures consistent with the scalar couplings, but removed from the conformation observed in solution. This arises from the lack of J restraints for the psi dihedral angle. Simulations with both nuclear Overhauser effect (NOE) and J-coupling restraints illustrates small but significant differences from simulations using only NOEs.     

High-resolution three-dimensional structure of reduced recombinant human thioredoxin in solution

The solution structure of recombinant human thioredoxin (105 residues) has been determined by nuclear magnetic resonance (NMR) spectroscopy combined with hybrid distance geometry-dynamical simulated annealing calculations. Approximate interproton distance restraints were derived from nuclear Overhauser effect (NOE) measurements. In addition, a large number of stereospecific assignments for beta-methylene protons and torsion angle restraints for phi, psi, and chi 1 were obtained by using a conformational grid search on the basis of the intraresidue and sequential NOE data in conjunction with 3JHN alpha and 3J alpha beta coupling constants. The structure calculations were based on 1983 approximate interproton distance restraints, 52 hydrogen-bonding restraints for 26 hydrogen bonds, and 98 phi, 71 psi, and 72 chi 1 torsion angle restraints. The 33 final simulated annealing structures obtained had an average atomic rms distribution of the individual structures about the mean coordinate positions of 0.40 +/- 0.06 A for the backbone atoms and 0.78 +/- 0.05 A for all atoms. The solution structure of human thioredoxin consists of a five-stranded beta-sheet surrounded by four alpha-helices, with an active site protrusion containing the two redox-active cysteines. The overall structure is similar to the crystal and NMR structures of oxidized [Katti, S. K., LeMaster, D. M., & Eklund, H. (1990) J. Mol. Biol. 212, 167-184] and reduced [Dyson, J. H., Gippert, G. P., Case, D. A., Holmgren, A., & Wright, P. (1990) Biochemistry 29, 4129-4136] Escherichia coli thioredoxin, respectively, despite the moderate 25% amino acid sequence homology. Several differences, however, can be noted. The human alpha 1 helix is a full turn longer than the corresponding helix in E. coli thioredoxin and is characterized by a more regular helical geometry. The helix labeled alpha 3 in human thioredoxin has its counterpart in the 3(10) helix of the E. coli protein and is also longer in the human protein. In contrast to these structural differences, the conformation of the active site loop in both proteins is very similar, reflecting the perfect sequence identity for a stretch of eight amino acid residues around the redox-active cysteines.     

Flexibility of bovine pancreatic trypsin inhibitor

The native conformation of a protein may be expressed in terms of the dihedral angles, phi's and psi's for the backbone, and kappa's for the side chains, for a given geometry (bond lengths and bond angles). We have developed a method to obtain the dihedral angles for a low-energy structure of a protein, starting with the X-ray structure; it is applied here to examine the degree of flexibility of bovine pancreatic trypsin inhibitor. Minimization of the total energy of the inhibitor (including nonbonded, electrostatic, torsional, hydrogen bonding, and disulfide loop energies) yields a conformation having a total energy of -221 kcal/mol and a root mean square deviation between all atoms of the computed and experimental structures of 0.63 A. The optimal conformation is not unique, however, there being at least two other conformations of low-energy (-222 and -220 kcal/mol), which resemble the experimental one (root mean square deviations of 0.66 and 0.64 A, respectively). These three conformations are located in different positions in phi, psi space, i.e., with a total deviation of 81 degrees, 100 degrees and 55 degrees from each other (with a root mean square deviation of several degrees per dihedral angle from each other). The nonbonded energies of the backbones, calculated along lines in phi, psi space connecting these three conformations, are all negative, without any intervening energy barriers (on an energy contour map in the phi, psi plane). Side chains were attached at several representative positions in this plane, and the total energy was minimized by varying the kappa's. The energies were of approximately the same magnitude as the previous ones, indicating that the conformation of low energy is flexible to some extent in a restricted region of phi, psi space. Interestingly, the difference delta phi i+1 in phi i+1 for the (i + 1)th residue from one conformation to another is approximately the same as -delta psi i for the ith residue; i.e., the plane of the peptide group between the ith and (i + 1)th residues re-orient without significant changes in the positions of the other atoms. The flexibility of the orientations of the planes of the peptide groups is probably coupled in a cooperative manner to the flexibility of the positions of the backbone and side-chain atoms.     

Secondary structure and dynamics of an intrinsically unstructured linker domain

The transient secondary structure and dynamics of an intrinsically unstructured linker domain from the 70 kDa subunit of human replication protein A was investigated using solution state NMR. Stable secondary structure, inferred from large secondary chemical shifts, was observed for a segment of the intrinsically unstructured linker domain when it is attached to an N-terminal protein interaction domain. Results from NMR relaxation experiments showed the rotational diffusion for this segment of the intrinsically unstructured linker domain to be correlated with the N-terminal protein interaction domain. When the N-terminal domain is removed, the stable secondary structure is lost and faster rotational diffusion is observed. The large secondary chemical shifts were used to calculate phi and psi dihedral angles and these dihedral angles were used to build a backbone structural model. Restrained molecular dynamics were performed on this new structure using the chemical shift based dihedral angles and a single NOE distance as restraints. In the resulting family of structures a large, solvent exposed loop was observed for the segment of the intrinsically unstructured linker domain that had large secondary chemical shifts.     

Theoretical conformational analysis of methylamide of N-acetyl-L-lysine]

The spatial structure of the methylamide of N-acetyl-L-lysine has been analysed taking into account non-bonded and electrostatic interactions, torsional energy, bond angles distortion and hydrogen bonding. Conformational capacities of the backbone and mutual dependence of spatial structures of the backbone and the side chain was described by conformational maps obtained by energy minimisation, the dihedral angles and the bond angles of the side chain being varied for every phi, psi point. Every possible combination for phi, psi, x1-x5-angles was used corresponding to the stable form of the backbone and to torsion potential minima of the initial approximations in the calculation of preferred conformations of the molecule. Comparisons are made between stable forms of the methylamide of N-acetyl-L-lysine and Lys residues in proteins with known structure.     

Modeling protein loops using a phi i + 1, psi i dimer database.

We present an automated method for modeling backbones of protein loops. The method samples a database of phi i + 1 and psi i angles constructed from a nonredundant version of the Protein Data Bank (PDB). The dihedral angles phi i + 1 and psi i completely define the backbone conformation of a dimer when standard bond lengths, bond angles, and a trans planar peptide configuration are used. For the 400 possible dimers resulting from 20 natural amino acids, a list of allowed phi i + 1, psi i pairs for each dimer is created by pooling all such pairs from the loop segments of each protein in the nonredundant version of the PDB. Starting from the N-terminus of the loop sequence, conformations are generated by assigning randomly selected pairs of phi i + 1, psi i for each dimer from the respective pool using standard bond lengths, bond angles, and a trans peptide configuration. We use this database to simulate protein loops of lengths varying from 5 to 11 amino acids in five proteins of known three-dimensional structures. Typically, 10,000-50,000 models are simulated for each protein loop and are evaluated for stereochemical consistency. Depending on the length and sequence of a given loop, 50-80% of the models generated have no stereochemical strain in the backbone atoms. We demonstrate that, when simulated loops are extended to include flanking residues from homologous segments, only very few loops from an ensemble of sterically allowed conformations orient the flanking segments consistent with the protein topology. The presence of near-native backbone conformations for loops from five different proteins suggests the completeness of the dimeric database for use in modeling loops of homologous proteins. Here, we take advantage of this observation to design a method that filters near-native loop conformations from an ensemble of sterically allowed conformations. We demonstrate that our method eliminates the need for a loop-closure algorithm and hence allows for the use of topological constraints of the homologous proteins or disulfide constraints to filter near-native loop conformations.     

Exact Solutions for Internuclear Vectors and Backbone Dihedral Angles from NH Residual Dipolar Couplings in Two Media, and their Application in a Systematic Search Algorithm for Determining Protein Backbone Structure

We have derived a quartic equation for computing the direction of an internuclear vector from residual dipolar couplings (RDCs) measured in two aligning media, and two simple trigonometric equations for computing the backbone (phi,psi) angles from two backbone vectors in consecutive peptide planes. These equations make it possible to compute, exactly and in constant time, the backbone (phi,psi) angles for a residue from RDCs in two media on any single backbone vector type. Building upon these exact solutions we have designed a novel algorithm for determining a protein backbone substructure consisting of alpha-helices and beta-sheets. Our algorithm employs a systematic search technique to refine the conformation of both alpha-helices and beta-sheets and to determine their orientations using exclusively the angular restraints from RDCs. The algorithm computes the backbone substructure employing very sparse distance restraints between pairs of alpha-helices and beta-sheets refined by the systematic search. The algorithm has been demonstrated on the protein human ubiquitin using only backbone NH RDCs, plus twelve hydrogen bonds and four NOE distance restraints. Further, our results show that both the global orientations and the conformations of alpha-helices and beta-strands can be determined with high accuracy using only two RDCs per residue. The algorithm requires, as its input, backbone resonance assignments, the identification of alpha-helices and beta-sheets as well as sparse NOE distance and hydrogen bond restraints.     

Low resolution structure of interleukin-1 beta in solution derived from 1H-15N heteronuclear three-dimensional nuclear magnetic resonance spectroscopy

A low resolution solution structure of the cytokine interleukin-1 beta, a 153 residue protein of molecular weight 17,400, has been determined on the basis of 446 nuclear Overhauser effect (NOE) derived approximate interproton distance restraints involving solely NH, C alpha H and C beta H protons, supplemented by 90 distance restraints for 45 hydrogen bonds, and 79 phi torsion angle restraints. With the exception of 27 C alpha H-C alpha H NOEs, all the NOEs were assigned from a three-dimensional 1H-1H NOE 15N-1H heteronuclear multiple quantum coherence (HMQC) spectrum. The torsion angle restraints were obtained from accurate 3JHN alpha coupling constants measured from a HMQC-J spectrum, while the hydrogen bonds were derived from a qualitative analysis of the NOE, coupling constant and amide exchange data. A total of 20 simulated annealing (SA) structures was computed using the hybrid distance geometry-dynamical simulated annealing method. The solution structure of IL-1 beta comprises 12 beta-strands arranged in three pseudo-symmetrical topological units (each consisting of 5 anti-parallel beta-strands), joined by turns, short loops and long loops. The core of the structure, which is made up of the 12 beta-strands, together with the turns joining strands I and II, strands VIII and IX and strands X and XI, is well determined with a backbone atomic root-mean-square (r.m.s.) distribution about the mean co-ordinate positions of 1.2(+/- 0.1) A. The loop conformations, on the other hand, are poorly determined by the current data. A comparison of the core of the low resolution solution structure of IL-1 beta with that of the X-ray structure indicates that they are similar, with a backbone atomic r.m.s. difference of only 1.5 A between the co-ordinates of the restrained minimized mean of the SA structures and the X-ray structure.     

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.     

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

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

Application of cross-correlated NMR spin relaxation to the zinc-finger protein CRP2(LIM2): evidence for collective motions in LIM domains

The solution structure of quail CRP2(LIM2) was significantly improved by using an increased number of NOE constraints obtained from a 13C,15N-labeled protein sample and by applying a recently developed triple-resonance cross-correlated relaxation experiment for the determination of the backbone dihedral angle psi. Additionally, the relative orientation of the 15N(i)-1HN(i) dipole and the 13CO(i) CSA tensor, which is related to both backbone angles phi and psi, was probed by nitrogen-carbonyl multiple-quantum relaxation and used as an additional constraint for the refinement of the local geometry of the metal-coordination sites in CRP2(LIM2). The backbone dynamics of residues located in the folded part of CRP2(LIM2) have been characterized by proton-detected 13C'(i-1)-15N(i) and 15N(i)-1HN(i) multiple-quantum relaxation, respectively. We show that regions having cross-correlated time modulation of backbone isotropic chemical shifts on the millisecond to microsecond time scale correlate with residues that are structurally altered in the mutant protein CRP2(LIM2)R122A (disruption of the CCHC zinc-finger stabilizing side-chain hydrogen bond) and that these residues are part of an extended hydrogen-bonding network connecting the two zinc-binding sites. This indicates the presence of long-range collective motions in the two zinc-binding subdomains. The conformational plasticity of the LIM domain may be of functional relevance for this important protein recognition motif.     

The influence of stereospecific assignments on the determination of three-dimensional structures of proteins by nuclear magnetic resonance spectroscopy. Application to the sea anemone protein BDS-I

The influence of the stereospecific assignments of beta-methylene protons and the classification of chi 1 torsion angles on the definition of the three-dimensional structures of proteins determined from NMR data is investigated using the sea anemone protein BDS-I (43 residues) as a model system. Two sets of structures are computed. The first set comprises 42 converged structures (denoted STEREO structures) calculated on the basis of the complete list of restraints derived from the NMR data, consisting of 489 interproton and 24 hydrogen bonding distance restraints, supplemented by 23 phi backbone and 21 chi 1 side chain torsion angle restraints. The second set comprises 31 converged structures (denoted NOSTEREO structures) calculated from a reduced data set in which those restraints arising from stereospecific assignments, and the corresponding chi 1 torsion angle restraints, are explicitly omitted. The results show that the inclusion of the stereospecific restraints leads to a significant improvement in the definition of the structure of BDS-I, both with respect to the backbone and the detailed arrangement of the side chains. Average atomic rms differences between the individual structures and the mean structures for the backbone atoms are 0.67 +/- 0.12 A and 0.93 +/- 0.16 A for the STEREO and NOSTEREO structures, respectively; the corresponding values for all atoms are 0.90 +/- 0.17 A and 1.17 +/- 0.17 A, respectively. In addition, while the overall fold remains unchanged, there is a small but significant atomic displacement between the two sets of structures.     

The geometry of alpha-sheet: Implications for its possible function as amyloid precursor in proteins

alpha-sheet has been proposed as the main constituent of the prefibrillar intermediate during amyloid formation. Here the helical parameters of the alpha-sheet strand are calculated from average main-chain dihedral angles reported from molecular dynamics simulations. It is an almost linear polypeptide that forms a right-handed helix of about 100 A diameter, with 100 residues and a rise of 30 A per turn. The strands are curved but untwisted, which implies that neighboring strands need not coil to make interstrand hydrogen bonds. This suggests that compared to beta-sheets in native folded proteins, alpha-sheets can be larger and stack more easily to create extensive 3D blocks. It is shown that alpha-sheet is related to a category of structures termed "mirror" structures. Mirror structures have repetitive pairs of main-chain dihedral angles at residues i and i+1 that satisfy the condition phi(i) (+1) = -psi(i), psi(i) (+1) = -phi(i). They are uniquely identified by the two orientations of their peptide planes, specified by phi(i) and psi(i). Their side chains point alternately in opposite directions. Interestingly, their conformations are insensitive to phi(i) and psi(i) in that the pseudo dihedral angle formed by four consecutive C(alpha) atoms is always close to 180 degrees . There are two types: "beta-mirror" and "alpha-mirror" structure; beta-mirror structures relate to beta-sheet by small peptide plane rotations, of less than 90 degrees , while alpha-mirror structures are close to alpha-sheet and relate to beta-sheet by approximately 180 degrees peptide plane flips. Most mirror structures, and in particular the alpha-mirror, form wide helices with diameters 50-70 A. Their gentle curvature, and therefore that of the alpha-sheet, arises from the orientation of successive peptide units causing the difference in the bond angles at the C and N atoms of the peptide unit to gradually change the direction of the chain.     

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.     

Real-value prediction of backbone torsion angles

The backbone structure of a protein is largely determined by the phi and psi torsion angles. Thus, knowing these angles, even if approximately, will be very useful for protein-structure prediction. However, in a previous work, a sequence-based, real-value prediction of psi angle could only achieve a mean absolute error of 54 degrees (83 degrees, 35 degrees, 33 degrees for coil, strand, and helix residues, respectively) between predicted and actual angles. Moreover, a real-value prediction of phi angle is not yet available. This article employs a neural-network based approach to improve psi prediction by taking advantage of angle periodicity and apply the new method to the prediction to phi angles. The 10-fold-cross-validated mean absolute error for the new method is 38 degrees (58 degrees, 33 degrees, 22 degrees for coil, strand, and helix, respectively) for psi and 25 degrees (35 degrees, 22 degrees, 16 degrees for coil, strand, and helix, respectively) for phi. The accuracy of real-value prediction is comparable to or more accurate than the predictions based on multistate classification of the phi-psi map. More accurate prediction of real-value angles will likely be useful for improving the accuracy of fold recognition and ab initio protein-structure prediction. The Real-SPINE 2.0 server is available on the website http://sparks.informatics.iupui.edu.     

Solution structure of a polypeptide dimer comprising the fourth Ca(2+)-binding site of troponin C by nuclear magnetic resonance spectroscopy

The structure of a 39 amino acid proteolytic fragment of rabbit skeletal troponin C containing the fourth Ca(2+)-binding site has been determined by an approach involving nuclear magnetic resonance (NMR) spectroscopy combined with hybrid distance geometry-dynamical simulated annealing calculations. Hydrodynamic and NMR evidence establishes unambiguously that the fragment forms a stable dimer in solution in the presence of excess Ca2+. The calculation of the dimeric structure is based on a total of 1056 experimental restraints comprising 422 interproton distances, 35 phi, 28 psi, and 28 chi 1 torsion angle restraints within each subunit, 30 intermonomer distance restraints, and 6 Ca2+ restraints per subunit. A total of 48 final structures were calculated having an rms deviation about the mean atomic backbone coordinate positions of 1.0 A for residues Asp128-Glu156. The solution structure consists of a dimer of helix-loop-helix motifs related by a 2-fold axis of symmetry. The overall architecture of the dimer is very similar to the C-terminal domain in the crystal structure of chicken skeletal troponin C.