Determination of a high-quality nuclear magnetic resonance solution structure of the bovine pancreatic trypsin inhibitor and comparison with three crystal structures.

A high-quality three-dimensional structure of the bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution was determined by 1H nuclear magnetic resonance (n.m.r.) spectroscopy and compared to the three available high-resolution X-ray crystal structures. A newly collected input of 642 distance constraints derived from nuclear Overhauser effects and 115 dihedral angle constraints was used for the structure calculations with the program DIANA, followed by restrained energy minimization with the program AMBER. The BPTI solution structure is represented by a group of 20 conformers with an average root-mean-square deviation (RMSD) relative to the mean solution structure of 0.43 A for backbone atoms and 0.92 A for all heavy atoms of residues 2 to 56. The pairwise RMSD values of the three crystal structures relative to the mean solution structure are 0.76 to 0.85 A for the backbone atoms and 1.24 to 1.33 A for all heavy atoms of residues 2 to 56. Small local differences in backbone atom positions between the solution structure and the X-ray structures near residues 9, 25 to 27, 46 to 48 and 52 to 58, and conformational differences for individual amino acid side-chains were analyzed for possible correlations with intermolecular protein-protein contacts in the crystal lattices, using the pairwise RMSD values among the three crystal structures as a reference.     

Solution structure of the catalytic domain of human stromelysin complexed with a hydrophobic inhibitor.

Stromelysin, a representative matrix metalloproteinase and target of drug development efforts, plays a prominent role in the pathological proteolysis associated with arthritis and secondarily in that of cancer metastasis and invasion. To provide a structural template to aid the development of therapeutic inhibitors, we have determined a medium-resolution structure of a 20-kDa complex of human stromelysin's catalytic domain with a hydrophobic peptidic inhibitor using multinuclear, multidimensional NMR spectroscopy. This domain of this zinc hydrolase contains a mixed beta-sheet comprising one antiparallel strand and four parallel strands, three helices, and a methionine-containing turn near the catalytic center. The ensemble of 20 structures was calculated using, on average, 8 interresidue NOE restraints per residue for the 166-residue protein fragment complexed with a 4-residue substrate analogue. The mean RMS deviation (RMSD) to the average structure for backbone heavy atoms is 0.91 A and for all heavy atoms is 1.42 A. The structure has good stereochemical properties, including its backbone torsion angles. The beta-sheet and alpha-helices of the catalytic domains of human stromelysin (NMR model) and human fibroblast collagenase (X-ray crystallographic model of Lovejoy B et al., 1994b, Biochemistry 33:8207-8217) superimpose well, having a pairwise RMSD for backbone heavy atoms of 2.28 A when three loop segments are disregarded. The hydroxamate-substituted inhibitor binds across the hydrophobic active site of stromelysin in an extended conformation. The first hydrophobic side chain is deeply buried in the principal S'1 subsite, the second hydrophobic side chain is located on the opposite side of the inhibitor backbone in the hydrophobic S'2 surface subsite, and a third hydrophobic side chain (P'3) lies at the surface.     

Solution structure and proposed domain domain recognition interface of an acyl carrier protein domain from a modular polyketide synthase

Polyketides are a medicinally important class of natural products. The architecture of modular polyketide synthases (PKSs), composed of multiple covalently linked domains grouped into modules, provides an attractive framework for engineering novel polyketide-producing assemblies. However, impaired domain-domain interactions can compromise the efficiency of engineered polyketide biosynthesis. To facilitate the study of these domain-domain interactions, we have used nuclear magnetic resonance (NMR) spectroscopy to determine the first solution structure of an acyl carrier protein (ACP) domain from a modular PKS, 6-deoxyerythronolide B synthase (DEBS). The tertiary fold of this 10-kD domain is a three-helical bundle; an additional short helix in the second loop also contributes to the core helical packing. Superposition of residues 14-94 of the ensemble on the mean structure yields an average atomic RMSD of 0.64 +/- 0.09 Angstrom for the backbone atoms (1.21 +/- 0.13 Angstrom for all non-hydrogen atoms). The three major helices superimpose with a backbone RMSD of 0.48 +/- 0.10 Angstrom (0.99 +/- 0.11 Angstrom for non-hydrogen atoms). Based on this solution structure, homology models were constructed for five other DEBS ACP domains. Comparison of their steric and electrostatic surfaces at the putative interaction interface (centered on helix II) suggests a model for protein-protein recognition of ACP domains, consistent with the previously observed specificity. Site-directed mutagenesis experiments indicate that two of the identified residues influence the specificity of ACP recognition.     

A large-scale comparison of computational models on the residue flexibility for NMR-derived proteins

As an alternative to X-ray crystallography, nuclear magnetic resonance (NMR) has also emerged as the method of choice for studying both protein structure and dynamics in solution. However, little work using computational models such as Gaussian network model (GNM) and machine learning approaches has focused on NMR-derived proteins to predict the residue flexibility, which is represented by the root mean square deviation (RMSD) with respect to the average structure. We provide a large-scale comparison of computational models, including GNM, parameter-free GNM and several linear regression models using local solvent exposures as inputs, based on a dataset of 1609 protein chains whose structures were resolved by NMR. The result again confirmed that the correlation of GNM outputs with raw RMSD values was better than that using B-factors of X-ray data. Nevertheless, it was also concluded that the parameter-free GNM and the solvent exposure based linear regression models performed worse than GNM when predicting RMSD, contrary to results using X-ray data. The discrepancy of residue flexibility prediction between NMR and X-ray data is likely attributable to a combination of their physical and methodological differences.     

NMR structure of tissue inhibitor of metalloproteinases-1 implicates localized induced fit in recognition of matrix metalloproteinases

A high quality solution structure of the matrix metalloproteinase inhibitory N-terminal domain of recombinant human tissue inhibitor of metalloproteinases-1 (N-TIMP-1) has been determined. For the rigidly packed residues, the average RMSD to the mean structure is 0. 57 A for the backbone atoms and 1.00 A for all heavy atoms. Comparison of the solution structure of free N-TIMP-1 with the crystal structure of TIMP-1 bound to the catalytic domain of MMP-3 ( Gomis-R]uth et al., 1997 ) shows that the structural core of the beta barrel flanked by helices is nearly unchanged by the association with MMP-3, evident from a backbone RMSD of 1.15 A. However, clear differences in the conformation of the MMP-binding ridge of free and MMP-bound TIMP-1 suggest induced fit throughout the ridge. The MMP-dependent conformational changes in the ridge include a dramatic bending of AB loop residues Glu28 through Leu34, moderate hinge bending of the CD-loop about residues Ala65 and Cys70, and modest bending of the Cys1 through Pro6 segment. A large number of interresidue Nuclear Overhauser enhancements (NOEs) augmented by stereospecific assignments, torsion restraints, and dipolar couplings (an average of 18 non-trivial restraints per residue) engender confidence in these structural inferences. A tight cluster of three lysine residues and one arginine residue atop beta-strands A and B, and identical among TIMP sequences, form the heart of a highly conserved electropositive patch that may interact with anionic components of the extracellular matrix.     

Refined solution structure of human profilin I.

Profilin is a ubiquitous eukaryotic protein that binds to both cytosolic actin and the phospholipid phosphatidylinositol-4,5-bisphosphate. These dual competitive binding capabilities of profilin suggest that profilin serves as a link between the phosphatidyl inositol cycle and actin polymerization, and thus profilin may be an essential component in the signaling pathway leading to cytoskeletal rearrangement. The refined three-dimensional solution structure of human profilin I has been determined using multidimensional heteronuclear NMR spectroscopy. Twenty structures were selected to represent the solution conformational ensemble. This ensemble of structures has root-mean-square distance deviations from the mean structure of 0.58 A for the backbone atoms and 0.98 A for all non-hydrogen atoms. Comparison of the solution structure of human profilin to the crystal structure of bovine profilin reveals that, although profilin adopts essentially identical conformations in both states, the solution structure is more compact than the crystal structure. Interestingly, the regions that show the most structural diversity are located at or near the actin-binding site of profilin. We suggest that structural differences are reflective of dynamical properties of profilin that facilitate favorable interactions with actin. The global folding pattern of human profilin also closely resembles that of Acanthamoeba profilin I, reflective of the 22% sequence identity and approximately 45% sequence similarity between these two proteins.     

Defining and computing optimum RMSD for gapped and weighted multiple-structure alignment

Pairwise structure alignment commonly uses root mean square deviation (RMSD) to measure the structural similarity, and methods for optimizing RMSD are well established. We extend RMSD to weighted RMSD for multiple structures. By using multiplicative weights, we show that weighted RMSD for all pairs is the same as weighted RMSD to an average of the structures. Thus, using RMSD or weighted RMSD implies that the average is a consensus structure. Although we show that in general, the two tasks of finding the optimal translations and rotations for minimizing weighted RMSD cannot be separated for multiple structures like they can for pairs, an inherent difficulty and a fact ignored by previous work, we develop a near-linear iterative algorithm to converge weighted RMSD to a local minimum. 10,000 experiments of gapped alignment done on each of 23 protein families from HOMSTRAD (where each structure starts with a random translation and rotation) converge rapidly to the same minimum. Finally we propose a heuristic method to iteratively remove the effect of outliers and find well-aligned positions that determine the structural conserved region by modeling B-factors and deviations from the average positions as weights and iteratively assigning higher weights to better aligned atoms.     

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.     

Packing of transmembrane helices in bacteriorhodopsin folding: structure and thermodynamics

We propose a coarse-grained (CG) model to study the native structure and physical properties of helical membrane proteins (HMPs) using off-lattice computer simulations. Instead of considering sequence heterogeneity explicitly, we model its effect on the packing of helices by employing a mean packing parameter r(0), which is calculated from an all-atom (AA) model. Specifically, this CG model is applied to investigate the packing of helices in bacteriorhodopsin (BR), and predicts the seven helix bundle structure of BR with a root mean square deviation (RMSD) in coordinates of helix backbone atoms (N, C, C(alpha)) of 3.99 A from its crystal structure. This predicted structure is further refined in an AA model by Amber and the refined structure has a RMSD (in coordinates of helix backbone atoms) of 2.64 A. The predicted packing position, tilting angle, and orientation angle of each helix in the refined structure are consistent with experimental data and their physical origins can be well understood in our model. Our results show that a reasonably good structure of BR can be predicted by using such a dual-scale approach, provided that its secondary structure is known. Starting from a random initial configuration, the folded structure can be obtained in days using a regular desktop computer. Various thermodynamic properties of helix packing of BR are also investigated in this CG model.     

Residue-level global and local ensemble-ensemble comparisons of protein domains

Many methods of protein structure generation such as NMR-based solution structure determination and template-based modeling do not produce a single model, but an ensemble of models consistent with the available information. Current strategies for comparing ensembles lose information because they use only a single representative structure. Here, we describe the ENSEMBLATOR and its novel strategy to directly compare two ensembles containing the same atoms to identify significant global and local backbone differences between them on per-atom and per-residue levels, respectively. The ENSEMBLATOR has four components: eePREP (ee for ensemble-ensemble), which selects atoms common to all models; eeCORE, which identifies atoms belonging to a cutoff-distance dependent common core; eeGLOBAL, which globally superimposes all models using the defined core atoms and calculates for each atom the two intraensemble variations, the interensemble variation, and the closest approach of members of the two ensembles; and eeLOCAL, which performs a local overlay of each dipeptide and, using a novel measure of local backbone similarity, reports the same four variations as eeGLOBAL. The combination of eeGLOBAL and eeLOCAL analyses identifies the most significant differences between ensembles. We illustrate the ENSEMBLATOR''s capabilities by showing how using it to analyze NMR ensembles and to compare NMR ensembles with crystal structures provides novel insights compared to published studies. One of these studies leads us to suggest that a “consistency check” of NMR-derived ensembles may be a useful analysis step for NMR-based structure determinations in general. The ENSEMBLATOR 1.0 is available as a first generation tool to carry out ensemble-ensemble comparisons.     

TASSER-based refinement of NMR structures

The TASSER structure prediction algorithm is employed to investigate whether NMR structures can be moved closer to their corresponding X-ray counterparts by automatic refinement procedures. The benchmark protein dataset includes 61 nonhomologous proteins whose structures have been determined by both NMR and X-ray experiments. Interestingly, by starting from NMR structures, the majority (79%) of TASSER refined models show a structural shift toward their X-ray structures. On average, the TASSER refined models have a root-mean-square-deviation (RMSD) from the X-ray structure of 1.785 A (1.556 A) over the entire chain (aligned region), while the average RMSD between NMR and X-ray structures (RMSD(NMR_X-ray)) is 2.080 A (1.731 A). For all proteins having a RMSD(NMR_X-ray) >2 A, the TASSER refined structures show consistent improvement. However, for the 34 proteins with a RMSD(NMR_X-ray) <2 A, there are only 21 cases (60%) where the TASSER model is closer to the X-ray structure than NMR, which may be due to the inherent resolution of TASSER. We also compare the TASSER models with 12 NMR models in the RECOORD database that have been recalculated recently by Nederveen et al. from original NMR restraints using the newest molecular dynamics tools. In 8 of 12 cases, TASSER models show a smaller RMSD to X-ray structures; in 3 of 12 cases, where RMSD(NMR_X-ray) <1 A, RECOORD does better than TASSER. These results suggest that TASSER can be a useful tool to improve the quality of NMR structures.     

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.     

Efficient molecular docking of NMR structures: application to HIV-1 protease

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.     

Solution structure of inhibitor-free human metalloelastase (MMP-12) indicates an internal conformational adjustment

Macrophage metalloelastase or matrix metalloproteinase-12 (MMP-12) appears to exacerbate atherosclerosis, emphysema, aortic aneurysm, rheumatoid arthritis, and inflammatory bowel disease. An inactivating E219A mutation, validated by crystallography and NMR spectra, prevents autolysis of MMP-12 and allows us to determine its NMR structure without an inhibitor. The structural ensemble of the catalytic domain without an inhibitor is based on 2813 nuclear Overhauser effects (NOEs) and has an average RMSD to the mean structure of 0.25 Å for the backbone and 0.61 Å for all heavy atoms for residues Trp109-Gly263. Compared to crystal structures of MMP-12, helix B (hB) at the active site is unexpectedly more deeply recessed under the β-sheet. This opens a pocket between hB and β-strand IV in the active-site cleft. Both hB and an internal cavity are shifted toward β-strand I, β-strand III, and helix A on the back side of the protease. About 25 internal NOE contacts distinguish the inhibitor-free solution structure and indicate hB's greater depth and proximity to the sheet and helix A. Line broadening and multiplicity of amide proton NMR peaks from hB are consistent with hB undergoing a slow conformational exchange among subtly different environments. Inhibitor-binding-induced perturbations of the NMR spectra of MMP-1 and MMP-3 map to similar locations across MMP-12 and encompass the internal conformational adjustments. Evolutionary trace analysis suggests a functionally important network of residues that encompasses most of the locations adjusting in conformation, including 18 residues with NOE contacts unique to inhibitor-free MMP-12. The conformational change, sequence analysis, and inhibitor perturbations of NMR spectra agree on the network they identify between structural scaffold and the active site of MMPs.     

Nuclear magnetic resonance solution structure of hirudin(1-51) and comparison with corresponding three-dimensional structures determined using the complete 65-residue hirudin polypeptide chain.

The three-dimensional structure of the N-terminal 51-residue domain of recombinant hirudin in aqueous solution was determined by 1H nuclear magnetic resonance (NMR) spectroscopy, and the resulting high-quality solution structure was compared with corresponding structures obtained from studies with the intact, 65-residue polypeptide chain of hirudin. On the basis of 580 distance constraints derived from nuclear Overhauser effects and 109 dihedral angle constraints, a group of 20 conformers representing the solution structure of hirudin(1-51) was computed with the program DIANA and energy-minimized with a modified version of the program AMBER. Residues 3 to 30 and 37 to 48 form a well-defined molecular core with two antiparallel beta-sheets composed of residues 14 to 16 and 20 to 22, and 27 to 31 and 36 to 40, and three reverse turns at residues 8 to 11 (type II), 17 to 20 (type II') and 23 to 26 (type II). The average root-mean-square deviation of the individual NMR conformers relative to their mean co-ordinates is 0.38 A for the backbone atoms and 0.77 A for all heavy atoms of these residues. Increased structural disorder was found for the N-terminal dipeptide segment, the loop at residues 31 to 36, and the C-terminal tripeptide segment. The solution structure of hirudin(1-51) has the same molecular architecture as the corresponding polypeptide segment in natural hirudin and recombinant desulfatohirudin. It is also closely similar to the crystal structure of the N-terminal 51-residue segment of hirudin in a hirudin-thrombin complex, with root-mean-square deviations of the crystal structure relative to the mean solution structure of 0.61 A for the backbone atoms and 0.91 A for all heavy atoms of residues 3 to 30 and 37 to 48. Further coincidence is found for the loop formed by residues 31 to 36, which shows increased structural disorder in all available solution structures of hirudin, and of which residues 32 to 35 are not observable in the electron density map of the thrombin complex. Significant local structural differences between hirudin(1-51) in solution and hirudin in the crystalline thrombin complex were identified mainly for the N-terminal tripeptide segment and residues 17 to 21. These are further analyzed in an accompanying paper.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

Heterogeneity and inaccuracy in protein structures solved by X-ray crystallography

Proteins are dynamic molecules, exhibiting structural heterogeneity in the form of anisotropic motion and discrete conformational substates, often of functional importance. In protein structure determination by X-ray crystallography, the observed diffraction pattern results from the scattering of X-rays by an ensemble of heterogeneous molecules, ordered and oriented by packing in a crystal lattice. The majority of proteins diffract to resolutions where heterogeneity is difficult to identify and model, and are therefore approximated by a single, average conformation with isotropic variance. Here we show that disregarding structural heterogeneity introduces degeneracy into the structure determination process, as many single, isotropic models exist that explain the diffraction data equally well. The large differences among these models imply that the accuracy of crystallographic structures has been widely overestimated. Further, it suggests that analyses that depend on small differences in the relative positions of atoms may be flawed.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

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.