Full-matrix least-squares refinement of lysozymes and analysis of anisotropic thermal motion

Crystal structures of turkey egg lysozyme (TEL) and human lysozyme (HL) were refined by full-matrix least-squares method using anisotropic temperature factors. The refinement converged at the conventional R-values of 0.104 (TEL) and 0.115 (HL) for reflections with F_o > 0 to the resolution of 1.12 Å and 1.15 Å, respectively. The estimated r.m.s. coordinate errors for protein atoms were 0.031 Å (TEL) and 0.034 Å (HL). The introduction of anisotropic temperature factors markedly reduced the R-value but did not significantly affect the main chain coordinates. The degree of anisotropy of atomic thermal motion has strong positive correlation with the square of distance from the molecular centroid. The ratio of the radial component of thermal ellipsoid to the r.m.s. magnitude of three principal components has negative correlation with the distance from the molecular centroid, suggesting the domination of libration rather than breathing motion. The TLS model was applied to elucidate the characteristics of the rigid-body motion. The TLS tensors were determined by the least-squares fit to observed temperature factors. The profile of the magnitude of reproduced temperature factors by the TLS method well fitted to that of observed B_eqv. However, considerable disagreement was observed in the shape and orientation of thermal ellipsoid for atoms with large temperature factors, indicating the large contribution of local motion. The upper estimate of the external motion, 67% (TEL) and 61% (HL) of B_eqv, was deduced from the plot of the magnitude of TLS tensors determined for main chain atoms which were grouped into shells according to the distance from the center of libration. In the external motion, the translational portion is predominant and the contribution of libration and screw motion is relatively small. The internal motion, estimated by subtracting the upper estimate of the external motion from the observed temperature factor, is very similar between TEL and HL in spite of the difference in 54 of 130 amino acid residues and in crystal packing, being suggested to reflect the intrinsic internal motion of chicken-type lysozymes. Proteins 30:232–243, 1998. © 1998 Wiley-Liss, Inc.     

Functional implications of interleukin-1 beta based on the three-dimensional structure.

The molecular structure of interleukin-1 beta, a hormone-like cytokine with roles in several disease processes, has been determined at 2.0 A resolution and refined to a crystallographic R-factor of 0.19. The framework of this molecule consists of 12 antiparallel beta-strands exhibiting pseudo-3-fold symmetry. Six of the strands make up a beta-barrel with polar residues concentrated at either end. Analysis of the three-dimensional structure, together with results from site-directed mutagenesis and biochemical and immunological studies, suggest that the core of the beta-barrel plays an important functional role. A large patch of charged residues on one end of the barrel is proposed as the binding surface with which IL-1 interacts with its receptor.     

The crystal structure of staphylococcal nuclease refined at 1.7 A resolution

The crystal structure of staphylococcal nuclease has been determined to 1.7 A resolution with a final R-factor of 16.2% using stereochemically restrained Hendrickson-Konnert least-squares refinement. The structure reveals a number of conformational changes relative to the structure of the ternary complex of staphylococcal nuclease 1,2 bound with deoxythymidine-3',5'-diphosphate and Ca2+. Tyr-113 and Tyr-115, which pack against the nucleotide base in the nuclease complex, are rotated outward creating a more open binding pocket in the absence of nucleotide. The side chains of Ca2+ ligands Asp-21 and Asp-40 shift as does Glu-43, the proposed general base in the hydrolysis of the 5'-phosphodiester bond. The significance of some changes in the catalytic site is uncertain due to the intrusion of a symmetry related Lys-70 side chain which hydrogen bonds to both Asp-21 and Glu-43. The position of a flexible loop centered around residue 50 is altered, most likely due to conformational changes propagated from the Ca2+ site. The side chains of Arg-35, Lys-84, Tyr-85, and Arg-87, which hydrogen bond to the 3'- and 5'-phosphates of the nucleotide in the nuclease complex, are unchanged in conformation, with packing interactions with adjacent protein side chains sufficient to fix the geometry in the absence of ligand. The nuclease structure presented here, in combination with the stereochemically restrained refinement of the nuclease complex structure at 1.65 A, provides a wealth of structural information for the increasing number of studies using staphylococcal nuclease as a model system of protein structure and function.     

The crystal structure of the "open" and the "closed" conformation of the flexible loop of trypanosomal triosephosphate isomerase

Triosephosphate isomerase has an important loop near the active site which can exist in a "closed" and in an "open" conformation. Here we describe the structural properties of this "flexible" loop observed in two different structures of trypanosomal triosephosphate isomerase. Trypanosomal triosephosphate isomerase, crystallized in the presence of 2.4 M ammonium sulfate, packs as an asymmetric dimer of 54,000 Da in the crystallographic asymmetric unit. Due to different crystal contacts, peptide 167-180 (the flexible loop of subunit-1) is an open conformation, whereas in subunit-2, this peptide (residues 467-480) is in a closed conformation. In the closed conformation, a hydrogen bond exists between the tip of the loop and a well-defined sulfate ion which is bound to the active site of subunit-2. Such an active site sulfate is not present in subunit-1 due to crystal contacts. When the native (2.4 M ammonium sulfate) crystals are transferred to a sulfate-free mother liquor, the flexible loop of subunit-2 adopts the open conformation. From a closed starting model, this open conformation was discovered through molecular dynamics refinement without manual intervention, despite involving C alpha shifts of up to 7 A. The tip of the loop, residues 472, 473, 474, and 475, moves as a rigid body. Our analysis shows that in this crystal form the flexible loop of subunit-2 faces a solvent channel. Therefore the open and the closed conformations of this flexible loop are virtually unaffected by crystal contacts. The actual observed conformation depends only on the absence or presence of a suitable ligand in the active site.     

The dynamic NMR structure of the TΨC-loop: Implications for the specificity of tRNA methylation

tRNA (m5U54)-methyltransferase (RUMT) catalyzes the S-adenosylmethionine-dependentmethylation of uridine-54 in the TC-loop of all transfer RNAs in E. coli to form the 54-ribosylthymine residue. However, in all tRNA structures, residue 54 is completely buried andthe question arises as to how RUMT gains access to the methylation site. A 17-mer RNAhairpin consisting of nucleotides 49–65 of the T-loop is a substrate for RUMT.Homonuclear NMR methods in conjunction with restrained molecular dynamics (MD)methods were used to determine the solution structure of the 17-mer T-arm fragment. Theloop of the hairpin exhibits enhanced flexibility which renders the conventional NMR averagestructure less useful compared to the more commonly found situation where a molecule existsin predominantly one major conformation. However, when resorting to softer refinementmethods such as MD with time-averaged restraints, the conflicting restraints in the loop canbe satisfied much better. The dynamic structure of the T-arm is represented as an ensembleof 10 time-clusters. In all of these, U54 is completely exposed. The flexibility of the T-loop in solution in conjunction with extensive binding studies of RUMT with the TC-loop and tRNA suggest that the specificity of the RUMT/tRNA recognition is associated withtRNA tertiary structure elements. For the methylation, RUMT would simply have to breakthe tertiary interactions between the D- and T-loops, leading to a melting of the T-armstructure and making U54 available for methylation.     

Protein solution structure calculations in solution: Solvated molecular dynamics refinement of calbindin D9k

The three-dimensional solution structures of proteins determinedwith NMR-derived constraints are almost always calculated in vacuo. Thesolution structure of (Ca²⁺)_{_2}-calbindinD_9k has been redetermined by new restrained molecular dynamics(MD) calculations that include Ca²⁺ ions and explicit solventmolecules. Four parallel sets of MD refinements were run to provide accuratecomparisons of structures produced in vacuo, in vacuo withCa²⁺ ions, and with two different protocols in a solvent bathwith Ca²⁺ ions. The structural ensembles were analyzed interms of structural definition, molecular energies, packing density,solvent-accessible surface, hydrogen bonds, and the coordination of calciumions in the two binding loops. Refinement including Ca²⁺ ionsand explicit solvent results in significant improvements in the precisionand accuracy of the structure, particularly in the binding loops. Theseresults are consistent with results previously obtained in free MDsimulations of proteins in solution and show that the rMD refinedNMR-derived solution structures of proteins, especially metalloproteins, canbe significantly improved by these strategies.     

Structure refinement using time-averaged J-coupling constant restraints

Summary We describe a new penalty function for use in restrained molecular dynamics simulations which allows experimental J-coupling information to be enforced as a time-averaged, rather than instantaneous, quantity. The pseudo-energy term has been formulated in terms of a calculated J value (a measured quantity) rather than the relevant dihedral angle (a derived quantity). This accounts for the distinct non-linearity of the coupling constant with respect to either Cartesian coordinates or dihedral angles. Example simulations of the cyclic decapeptide antamanide show the procedure's ability to enforce experimental restraints while exploring a large region of conformational space and producing a relatively small disturbance of the physical force field.Abbreviations MD molecular dynamics - rms root-mean-square - NOE nuclear Overhauser effect     

Direct NOE refinement of biomolecular structures using 2D NMR data

Summary A set of computer programs called DINOSAUR has been developed, which allows the refinement of biomolecular structures directly from 2D NOE intensities. The NOE restraining potential implemented emphasises the weak intensities corresponding to larger distances which are more likely to determine the three-dimensional structure. An approximation based on a two-spin approximation is proposed for the gradient of the NOE intensities instead of the exact solution which is extremely time-consuming. The DINOSAUR routines have been implemented in various refinement programs (Distance bound Driven Dynamics, Molecular Dynamics and Energy Minimisation) and tested on an eight-residue model peptide.     

Refining homology models by combining replica-exchange molecular dynamics and statistical potentials

A protocol is presented for the global refinement of homology models of proteins. It combines the advantages of temperature-based replica-exchange molecular dynamics (REMD) for conformational sampling and the use of statistical potentials for model selection. The protocol was tested using 21 models. Of these 14 were models of 10 small proteins for which high-resolution crystal structures were available, the remainder were targets of the recent CASPR exercise. It was found that REMD in combination with currently available force fields could sample near-native conformational states starting from high-quality homology models. Conformations in which the backbone RMSD of secondary structure elements (SSE-RMSD) was lower than the starting value by 0.5-1.0 A were found for 15 out of the 21 cases (average 0.82 A). Furthermore, when a simple scoring function consisting of two statistical potentials was used to rank the structures, one or more structures with SSE-RMSD of at least 0.2 A lower than the starting value was found among the five best ranked structures in 11 out of the 21 cases. The average improvement in SSE-RMSD for the best models was 0.42 A. However, none of the scoring functions tested identified the structures with the lowest SSE-RMSD as the best models although all identified the native conformation as the one with lowest energy. This suggests that while the proposed protocol proved effective for the refinement of high-quality models of small proteins scoring functions remain one of the major limiting factors in structure refinement. This and other aspects by which the methodology could be further improved are discussed.