Binding of retinol induces changes in rat cellular retinol-binding protein II conformation and backbone dynamics.

The structure and backbone dynamics of rat holo cellular retinol-binding protein II (holo-CRBP II) in solution has been determined by multidimensional NMR. The final structure ensemble was based on 3980 distance and 30 dihedral angle restraints, and was calculated using metric matrix distance geometry with pairwise Gaussian metrization followed by simulated annealing. The average RMS deviation of the backbone atoms for the final 25 structures relative to their mean coordinates is 0.85(+/-0.09) A. Comparison of the solution structure of holo-CRBP II with apo-CRBP II indicates that the protein undergoes conformational changes not previously observed in crystalline CRBP II, affecting residues 28-35 of the helix-turn-helix, residues 37-38 of the subsequent linker, as well as the beta-hairpin C-D, E-F and G-H loops. The bound retinol is completely buried inside the binding cavity and oriented as in the crystal structure. The order parameters derived from the (15)N T(1), T(2) and steady-state NOE parameters show that the backbone dynamics of holo-CRBP II is restricted throughout the polypeptide. The T(2) derived apparent backbone exchange rate and amide (1)H exchange rate both indicate that the microsecond to second timescale conformational exchange occurring in the portal region of the apo form has been suppressed in the holo form.

Ab initio fold prediction of small helical proteins using distance geometry and knowledge-based scoring functions.

The problem of protein tertiary structure prediction from primary sequence can be separated into two subproblems: generation of a library of possible folds and specification of a best fold given the library. A distance geometry procedure based on random pairwise metrization with good sampling properties was used to generate a library of 500 possible structures for each of 11 small helical proteins. The input to distance geometry consisted of sets of restraints to enforce predicted helical secondary structure and a generic range of 5 to 11 A between predicted contact residues on all pairs of helices. For each of the 11 targets, the resulting library contained structures with low RMSD versus the native structure. Near-native sampling was enhanced by at least three orders of magnitude compared to a random sampling of compact folds. All library members were scored with a combination of an all-atom distance-dependent function, a residue pair-potential, and a hydrophobicity function. In six of the 11 cases, the best-ranking fold was considered to be near native. Each library was also reduced to a final ab initio prediction via consensus distance geometry performed over the 50 best-ranking structures from the full set of 500. The consensus results were of generally higher quality, yielding six predictions within 6.5 A of the native fold. These favorable predictions corresponded to those for which the correlation between the RMSD and the scoring function were highest. The advantage of the reported methodology is its extreme simplicity and potential for including other types of structural restraints.

The structure and dynamics of rat apo-cellular retinol-binding protein II in solution: comparison with the X-ray structure.

The structure and dynamics of rat apo-cellular retinol binding protein II (apo-CRBP II) in solution has been determined by multidimensional NMR analysis of uniformly enriched recombinant rat 13C, 15N-apo-CRBP II and 15N-apo-CRBP II. The final ensemble of 24 NMR structures has been calculated from 3274 conformational restraints or 24.4 restraints/residue. The average root-mean-square deviation of the backbone atoms for the final 24 structures relative to their mean structure is 1.06 A. Although the average solution structure is very similar to the crystal structure, it differs at the putative entrance to the binding cavity, which is formed by the helix-turn-helix motif, the betaC-betaD turn and the betaE-betaF turn. The mean coordinates of the main-chain atoms of amino acid residues 28-38 are displaced in the solution structure relative to the crystal structure. The side-chain of F58, located on the betaC-betaD turn, is reoriented such that it interacts with L37 and no longer blocks entry into the ligand-binding pocket. Residues 28-35, which form the second helix of the helix-turn-helix motif in the crystal structure, do not exhibit a helical conformation in the solution structure. The solution structure of apo-CRBP II exhibits discrete regions of backbone disorder which are most pronounced at residues 28-32, 37-38 and 73-76 in the betaE-betaF turn as evaluated by the consensus chemical shift index, the root-mean-square deviation, amide 1H exchange rates and 15N relaxation studies. These studies indicate that fluctuations in protein conformation occur on the microseconds to ms time-scale in these regions of the protein. Some of these exchange processes can be directly observed in the three-dimensional 15N-resolved NOESY spectrum. These results suggest that in solution, apo-CRBP II undergoes conformational changes on the microseconds to ms time-scale which result in increased access to the binding cavity.

A potential smoothing algorithm accurately predicts transmembrane helix packing.

Potential smoothing, a deterministic analog of stochastic simulated annealing, is a powerful paradigm for the solution of conformational search problems that require extensive sampling, and should be a useful tool in computational approaches to structure prediction and refinement. A novel potential smoothing and search (PSS) algorithm has been developed and applied to predict the packing of transmembrane helices. The highlight of this method is the efficient manner in which it circumvents the combinatorial explosion associated with the large number of minima on multidimensional potential energy surfaces in order to converge to the global energy minimum. Here we show how our potential smoothing and search method succeeds in finding the global minimum energy structure for the glycophorin A (GpA) transmembrane helix dimer by optimizing interhelical van der Waals interactions over rigid and semi-rigid helices. Structures obtained from our ab initio predictions are in close agreement with recent experimental data.

Accuracy of side-chain prediction upon near-native protein backbones generated by Ab initio folding methods.

The ab initio folding problem can be divided into two sequential tasks of approximately equal computational complexity: the generation of native-like backbone folds and the positioning of side chains upon these backbones. The prediction of side-chain conformation in this context is challenging, because at best only the near-native global fold of the protein is known. To test the effect of displacements in the protein backbones on side-chain prediction for folds generated ab initio, sets of near-native backbones (< or = 4 A C alpha RMS error) for four small proteins were generated by two methods. The steric environment surrounding each residue was probed by placing the side chains in the native conformation on each of these decoys, followed by torsion-space optimization to remove steric clashes on a rigid backbone. We observe that on average 40% of the chi1 angles were displaced by 40 degrees or more, effectively setting the limits in accuracy for side-chain modeling under these conditions. Three different algorithms were subsequently used for prediction of side-chain conformation. The average prediction accuracy for the three methods was remarkably similar: 49% to 51% of the chi1 angles were predicted correctly overall (33% to 36% of the chi1+2 angles). Interestingly, when the inter-side-chain interactions were disregarded, the mean accuracy increased. A consensus approach is described, in which side-chain conformations are defined based on the most frequently predicted chi angles for a given method upon each set of near-native backbones. We find that consensus modeling, which de facto includes backbone flexibility, improves side-chain prediction: chi1 accuracy improved to 51-54% (36-42% of chi1+2). Implications of a consensus method for ab initio protein structure prediction are discussed.

Distance geometry generates native-like folds for small helical proteins using the consensus distances of predicted protein structures.

For successful ab initio protein structure prediction, a method is needed to identify native-like structures from a set containing both native and non-native protein-like conformations. In this regard, the use of distance geometry has shown promise when accurate inter-residue distances are available. We describe a method by which distance geometry restraints are culled from sets of 500 protein-like conformations for four small helical proteins generated by the method of Simons et al. (1997). A consensus-based approach was applied in which every inter-Calpha distance was measured, and the most frequently occurring distances were used as input restraints for distance geometry. For each protein, a structure with lower coordinate root-mean-square (RMS) error than the mean of the original set was constructed; in three cases the topology of the fold resembled that of the native protein. When the fold sets were filtered for the best scoring conformations with respect to an all-atom knowledge-based scoring function, the remaining subset of 50 structures yielded restraints of higher accuracy. A second round of distance geometry using these restraints resulted in an average coordinate RMS error of 4.38 A.

Light-activated rhodopsin induces structural binding motif in G protein alpha subunit.

A large superfamily of transmembrane receptors control cellular responses to diverse extracellular signals by catalyzing activation of specific types of heterotrimeric GTP-binding proteins. How these receptors recognize and promote nucleotide exchange on G protein alpha subunits to initiate signal amplification is unknown. The three-dimensional structure of the transducin (Gt) alpha subunit C-terminal undecapeptide Gtalpha(340-350) IKENLKDCGLF was determined by transferred nuclear Overhauser effect spectroscopy while it was bound to photoexcited rhodopsin. Light activation of rhodopsin causes a dramatic shift from a disordered conformation of Gtalpha(340-350) to a binding motif with a helical turn followed by an open reverse turn centered at Gly-348, a helix-terminating C capping motif of an alphaL type. Docking of the NMR structure to the GDP-bound x-ray structure of Gt reveals that photoexcited rhodopsin promotes the formation of a continuous helix over residues 325-346 terminated by the C-terminal helical cap with a unique cluster of crucial hydrophobic side chains. A molecular mechanism by which activated receptors can control G proteins through reversible conformational changes at the receptor-G protein interface is demonstrated.

Type 2M von Willebrand disease: F606I and I662F mutations in the glycoprotein Ib binding domain selectively impair ristocetin- but not botrocetin-mediated binding of von Willebrand factor to platelets.

von Willebrand disease (vWD) is a common, autosomally inherited, bleeding disorder caused by quantitative and/or qualitative deficiency of von Willebrand factor (vWF). We describe two families with a variant form of vWD where affected members of both families have borderline or low vWF antigen levels, normal vWF multimer patterns, disproportionately low ristocetin cofactor activity, and significant bleeding symptoms. Whereas ristocetin-induced binding of plasma vWF from affected members of both families to fixed platelets was reduced, botrocetin-induced platelet binding was normal. The sequencing of genomic DNA identified unique missense mutations in each family in the vWF exon 28. In Family A, a missense mutation at nucleotide 4105T --> A resulted in a Phe606Ile amino acid substitution (F606I) and in Family B, a missense mutation at nucleotide 4273A --> T resulted in an Ile662Phe amino acid substitution (I662F). Both mutations are within the large disulfide loop between Cys509 and Cys695 in the A1 domain that mediates vWF interaction with platelet glycoprotein Ib. Expression of recombinant vWF containing either F606I or I662F mutations resulted in mutant recombinant vWF with decreased ristocetin-induced platelet binding, but normal multimer structure, botrocetin-induced platelet binding, collagen binding, and binding to the conformation-sensitive monoclonal antibody, AvW-3. Both mutations are phenotypically distinct from the previously reported variant type 2MMilwaukee-1 because of the presence of normal botrocetin-induced platelet binding, collagen binding, and AvW-3 binding, as well as the greater frequency and intensity of clinical bleeding. When the reported type 2M mutations are mapped on the predicted three-dimensional structure of the A1 loop of vWF, the mutations cluster in one region that is distinct from the region in which the type 2B mutations cluster.

The NMR solution structure of intestinal fatty acid-binding protein complexed with palmitate: application of a novel distance geometry algorithm.

The three-dimensional solution structure of rat intestinal fatty acid-binding protein (I-FABP) complexed with palmitate has been determined using multidimensional triple-resonance NMR methods. The structure is based on 3889 conformational restraints derived mostly from 3-D 13C- and 15N-resolved nuclear Overhauser (NOESY) experiments. The 3-D NOESY data for this 15.4 kDa complex contained an average of nine possible interpretations per cross-peak. To circumvent this ambiguity, an eight-stage iterative procedure was employed to gradually interpret and introduce unambiguous distance restraints during subsequent rounds of structure calculations. The first stage of this procedure relied critically upon an initial structural model based on the consensus 1H/13C chemical shift-derived secondary structure and a set of symmetry-checked restraints derived from the 3-D 13C-resolved NOESY spectrum. The structures were calculated using DISTGEOM, a program that implements a novel distance geometry algorithm with pairwise Gaussian metrization. A central feature of this algorithm is the use of an iteratively optimized Gaussian distribution for the selection of trial distances, which overcomes the tendency of metrization to produce crushed structures. In addition, this algorithm randomly selects pairwise elements of the distance matrix, which results in an improved sampling of conformational space for a given computational effort. The final family of 20 distance geometry/simulated annealing structures exhibited an average pairwise C(alpha) root-mean-square deviation of 0.98 A, and their stereochemical quality, as assessed by PROCHECK, was comparable to that of 2.5 A X-ray crystal structures. The NMR structure was compared with the X-ray crystal structure of the same ligand/protein complex and was found to be essentially identical within the precision of the results. The NMR structure was also compared with that of the palmitate complex with bovine heart FABP, which shares 30% sequence identity with rat I-FABP. The overall folds were the same, but differences were noted with respect to the presence or absence of apparent conformational heterogeneity and the location and conformation of the bound fatty acid.

Definition of a metal-dependent/Li(+)-inhibited phosphomonoesterase protein family based upon a conserved three-dimensional core structure.

Inositol polyphosphate 1-phosphatase, inositol monophosphate phosphatase, and fructose 1,6-bisphosphatase share a sequence motif, Asp-Pro-(Ile or Leu)-Asp-(Gly or Ser)-(Thr or Ser), that has been shown by crystallographic and mutagenesis studies to bind metal ions and participate in catalysis. We compared the six alpha-carbon coordinates of this motif from the crystal structures of these three phosphatases and found that they are superimposable with rms deviations ranging from 0.27 to 0.60 A. Remarkably, when these proteins were aligned by this motif a common core structure emerged, defined by five alpha-helices and 11 beta-strands comprising 155 residues having rms deviations ranging from 1.48 to 2.66 A. We used the superimposed structures to align the sequences within the common core, and a distant relationship was observed suggesting a common ancestor. The common core was used to align the sequences of several other proteins that share significant similarity to inositol monophosphate phosphatase, including proteins encoded by fungal qa-X and qutG, bacterial suhB and cysQ (identical to amtA), and yeast met22 (identical to hal2). Evolutionary comparison of the core sequences indicate that five distinct branches exist within this family. These proteins share metal-dependent/Li(+)-sensitive phosphomonoesterase activity, and each predicted tree branch exhibits unique substrate specificity. Thus, these proteins define an ancient structurally conserved family involved in diverse metabolic pathways including inositol signaling, gluconeogenesis, sulfate assimilation, and possibly quinone metabolism. Furthermore, we suggest that this protein family identifies candidate enzymes to account for both the therapeutic and toxic actions of Li+ as it is used in patients treated for manic depressive disease.

Crystal structure of inositol polyphosphate 1-phosphatase at 2.3-A resolution.

Bovine inositol polyphosphate 1-phosphatase (1-ptase), M(r) = 44,000, is a Mg(2+)-dependent/Li(+)-sensitive enzyme that catalyzes the hydrolysis of the 1-position phosphate from inositol 1,4-bisphosphate and inositol 1,3,4-trisphosphate. We have determined the crystal structure of recombinant bovine 1-ptase in the presence of Mg2+ by multiple isomorphous replacement. The structure is currently refined to an R value of 0.198 for 15,563 reflections within a resolution range of 8.0-2.3 A. 1-Ptase is monomeric in the crystal, consistent with biochemical data, and folds into an alternatively layered alpha/beta/alpha/beta sandwich. The central core of 1-ptase consists of a six-stranded antiparallel beta sheet perpendicular to two parallel three-turn alpha-helices. The beta sheet is flanked by two antiparallel six-turn alpha-helices aligned parallel to the beta sheet, and the central helices are flanked by a five-stranded largely parallel beta sheet. Two neighboring metal binding sites are located in adjacent acidic pockets formed by the intersection of several secondary structure elements including an unusual kink structure formed by the "DPIDST" sequence motif. The fold of 1-ptase is similar to that of two other metal-dependent/Li(+)-sensitive phosphatases, inositol monophosphate phosphatase and fructose 1,6-bisphosphatase despite minimal amino acid identity. Comparison of the active-site pockets of these proteins will likely provide insight into substrate binding and the mechanisms of metal-dependent catalysis and Li+ inhibition.

Crystallization and initial X-ray crystallographic characterization of recombinant bovine inositol polyphosphate 1-phosphatase produced in Spodoptera frugiperda cells.

Bovine inositol polyphosphate 1-phosphatase, a monomeric protein with a molecular mass of 44,000 Da, hydrolyzes the 1-position phosphate from inositol 1,3,4-trisphosphate and inositol 1,4-bisphosphate. The low abundance of inositol polyphosphate 1-phosphatase in tissues has precluded structural studies requiring large quantities of enzyme. We used recombinant Baculovirus harboring the cDNA of bovine inositol polyphosphate 1-phosphatase to infect Spodoptera frugiperda (Sf9) insect cells. Recombinant protein (25 mg per 1 x 10(9) cells) was purified to homogeneity. The enzyme produced in Sf9 cells was similar to the native purified protein as determined by immunoblotting catalytic properties, and inhibition by lithium ions. Crystals of the purified recombinant enzyme were grown by vapor diffusion. Precession photography was used to determine the parameters of inositol polyphosphate 1-phosphatase crystals. The tetragonal crystals belong to the space group P4(1) or P4(3), have unit cell dimensions of a = b = 51.6 A, c = 143.3 A, alpha = beta = gamma = 90 degrees, and contain one molecule per asymmetric unit. We have collected a complete diffraction data set extending to 2.3 A and are currently attempting to solve the three-dimensional structure of bovine inositol polyphosphate 1-phosphatase using a multiple isomorphous replacement strategy.

Tertiary templates for proteins. Use of packing criteria in the enumeration of allowed sequences for different structural classes.

We assume that each class of protein has a core structure that is defined by internal residues, and that the external, solvent-contacting residues contribute to the stability of the structure, are of primary importance to function, but do not determine the architecture of the core portions of the polypeptide chain. An algorithm has been developed to supply a list of permitted sequences of internal residues compatible with a known core structure. This list is referred to as the tertiary template for that structure. In general the positions in the template are not sequentially adjacent and are distributed throughout the polypeptide chain. The template is derived using the fixed positions for the main-chain and beta-carbon atoms in the test structure and selected stereochemical rules. The focus of this paper is on the use of two packing criteria: avoidance of steric overlap and complete filling of available space. The program also notes potential polar group interactions and disulfide bonds as well as possible burial of formal charges. Central to the algorithm is the side-chain rotamer library. In an update of earlier studies by others, we show that 17 of the 20 amino acids (omitting Met, Lys and Arg) can be represented adequately by 67 side-chain rotamers. A list of chi angles and their standard deviations is given. The newer, high-resolution, refined structures in the Brookhaven Protein Data Bank show similar mean chi values, but have much smaller deviations than those of earlier studies. This suggests that a rotamer library may be a better structural approximation than was previously thought. In using packing constraints, it has been found essential to include all hydrogen atoms specifically. The "unified atom" representation is not adequate. The permitted rotamer sequences are severely restricted by the main-chain plus beta-carbon atoms of the test structure. Further restriction is introduced if the full set of atoms of the external residues are held fixed, the full-chain model. The space-filling requirement has a major role in restricting the template lists. The preliminary tests reported here make it appear likely that templates prepared from the currently known core structures will be able to discriminate between these structures. The templates should thus be useful in deciding whether a sequence of unknown tertiary structure fits any of the known core classes and, if a fit is found, how the sequence should be aligned in three dimensions to fit the core of that class.(ABSTRACT TRUNCATED AT 400 WORDS)