Combining experimental information from crystal and solution studies: joint X-ray and NMR refinement

Therapeutic targeting of IL-6 trans-signaling

Interleukin-6 (IL-6) is a cytokine, which is involved in innate and acquired immunity, in neural cell maintenance and in metabolism. IL-6 can be synthesized by many different cells including myeloid cells, fibroblasts, endothelial cells and lymphocytes. The synthesis of IL-6 is strongly stimulated by Toll like receptors and by IL-1. Therefore, IL-6 levels in the body are high during infection and inflammatory processes. Moreover, IL-6 is a prominent growth factor of tumor cells and plays a major role in inflammation associated cancer. On target cells, IL-6 binds to an IL-6 receptor, which is not signaling competent. The complex of IL-6 and IL-6 receptor associate with a second receptor subunit, glycoprotein gp130, which dimerizes and initiates intracellular signaling. Cells, which do not express the IL-6 receptor are not responsive to IL-6. They can, however, be stimulated by the complex of IL-6 and a soluble form of the IL-6 receptor, which is generated by limited proteolysis and to a lesser extent by translation from an alternatively spliced mRNA. This process has been named IL-6 trans-signaling. This review article will explain the biology of IL-6 trans-signaling and the specific inhibition of this mode of signaling, which has been recognized to be fundamental in inflammation and cancer.

Common mechanism of thermostability in small α- and β-proteins studied by molecular dynamics

Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five β-proteins with known, distinct structures and thermostability (T_m ) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to T_m . Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and β-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH - TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much.

On the complementarity of X-ray and NMR data

X-ray crystallography and NMR contain complementary information for the structural characterization of biological macromolecules. X-ray diffraction is primarily sensitive to the overall shape of the molecule, whereas NMR is mostly sensitive to the atomic detail. Their combination can therefore provide a stronger justification for the resulting structure. For their combination we have recently proposed REFMAC-NMR, which relies on primary data from both techniques for joint refinement. This possibility raises the compelling question of how far the complementarity can be extended. In this paper, we describe an integrative approach to the refinement with NMR data of four X-ray structures of hen-egg-white lysozyme, solved at atomic resolution in four different crystal forms, and we demonstrate that the outcome critically depends on the crystal form itself, reflecting the sensitivity of NMR to fine details.

Integrated BioNMR - "getting by with a little help from my friends"

Single types of methodologies are insufficient to adequately describe complex biological structures. As a result, integrated approaches that combine complementary data are being developed. Here, I describe the benefits of integrating solution and magic angle spinning BioNMR approaches to characterize structure and dynamics of protein assemblies.

Joint X-ray/NMR structure refinement of multidomain/multisubunit systems

Data integration in structural biology has become a paradigm for the characterization of biomolecular systems, and it is now accepted that combining different techniques can fill the gaps in each other's blind spots. In this frame, one of the combinations, which we have implemented in REFMAC-NMR, is residual dipolar couplings from NMR together with experimental data from X-ray diffraction. The first are exquisitely sensitive to the local details but does not give any information about overall shape, whereas the latter encodes more the information about the overall shape but at the same time tends to miss the local details even at the highest resolutions. Once crystals are obtained, it is often rather easy to obtain a complete X-ray dataset, however it is time-consuming to obtain an exhaustive NMR dataset. Here, we discuss the effect of including a-priori knowledge on the properties of the system to reduce the number of experimental data needed to obtain a more complete picture. We thus introduce a set of new features of REFMAC-NMR that allow for improved handling of RDC data for multidomain proteins and multisubunit biomolecular complexes, and encompasses the use of pseudo-contact shifts as an additional source of NMR-based information. The new feature may either help in improving the refinement, or assist in spotting differences between the crystal and the solution data. We show three different examples where NMR and X-ray data can be reconciled to a unique structural model without invoking mobility.

BMI1 regulates PRC1 architecture and activity through homo- and hetero-oligomerization

BMI1 is a core component of the polycomb repressive complex 1 (PRC1) and emerging data support a role of BMI1 in cancer. The central domain of BMI1 is involved in protein–protein interactions and is essential for its oncogenic activity. Here, we present the structure of BMI1 bound to the polyhomeotic protein PHC2 illustrating that the central domain of BMI1 adopts an ubiquitin-like (UBL) fold and binds PHC2 in a β-hairpin conformation. Unexpectedly, we find that the UBL domain is involved in homo-oligomerization of BMI1. We demonstrate that both the interaction of BMI1 with polyhomeotic proteins and homo-oligomerization via UBL domain are necessary for H2A ubiquitination activity of PRC1 and for clonogenic potential of U2OS cells. Here, we also emphasize need for joint application of NMR spectroscopy and X-ray crystallography to determine the overall structure of the BMI1–PHC2 complex.

BMI1, a core element of the polycomb repressive complex 1, is suggested to have oncogenic activity in a variety of cancers. Here, the authors report the structure of BMI1 bound to the protein PHC2, identify BMI1 homo-oligomerization interfaces, and analyse the role of BMI1 protein-protein interactions in PRC1 function.

How to tackle protein structural data from solution and solid state: An integrated approach

Long-range NMR restraints, such as diamagnetic residual dipolar couplings and paramagnetic data, can be used to determine 3D structures of macromolecules. They are also used to monitor, and potentially to improve, the accuracy of a macromolecular structure in solution by validating or "correcting" a crystal model. Since crystal structures suffer from crystal packing forces they may not be accurate models for the macromolecular structures in solution. However, the presence of real differences should be tested for by simultaneous refinement of the structure using both crystal and solution NMR data. To achieve this, the program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic and paramagnetic NMR data and/or diamagnetic residual dipolar couplings. Inconsistencies between crystal structures and solution NMR data, if any, may be due either to structural rearrangements occurring on passing from the solution to solid state, or to a greater degree of conformational heterogeneity in solution with respect to the crystal. In the case of multidomain proteins, paramagnetic restraints can provide the correct mutual orientations and positions of domains in solution, as well as information on the conformational variability experienced by the macromolecule.

Simultaneous use of solution NMR and X-ray data in REFMAC5 for joint refinement/detection of structural differences

The program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic data and paramagnetic NMR data (pseudocontact shifts and self-orientation residual dipolar couplings) and/or diamagnetic residual dipolar couplings. Incorporation of these long-range NMR restraints in REFMAC5 can reveal differences between solid-state and solution conformations of molecules or, in their absence, can be used together with X-ray crystallographic data for structural refinement. Since NMR and X-ray data are complementary, when a single structure is consistent with both sets of data and still maintains reasonably `ideal' geometries, the reliability of the derived atomic model is expected to increase. The program was tested on five different proteins: the catalytic domain of matrix metalloproteinase 1, GB3, ubiquitin, free calmodulin and calmodulin complexed with a peptide. In some cases the joint refinement produced a single model consistent with both sets of observations, while in other cases it indicated, outside the experimental uncertainty, the presence of different protein conformations in solution and in the solid state.

Integrated description of protein dynamics from room-temperature X-ray crystallography and NMR

Detailed descriptions of atomic coordinates and motions are required for an understanding of protein dynamics and their relation to molecular recognition, catalytic function, and allostery. Historically, NMR relaxation measurements have played a dominant role in the determination of the amplitudes and timescales (picosecond-nanosecond) of bond vector fluctuations, whereas high-resolution X-ray diffraction experiments can reveal the presence of and provide atomic coordinates for multiple, weakly populated substates in the protein conformational ensemble. Here we report a hybrid NMR and X-ray crystallography analysis that provides a more complete dynamic picture and a more quantitative description of the timescale and amplitude of fluctuations in atomic coordinates than is obtainable from the individual methods alone. Order parameters (S(2)) were calculated from single-conformer and multiconformer models fitted to room temperature and cryogenic X-ray diffraction data for dihydrofolate reductase. Backbone and side-chain order parameters derived from NMR relaxation experiments are in excellent agreement with those calculated from the room-temperature single-conformer and multiconformer models, showing that the picosecond timescale motions observed in solution occur also in the crystalline state. These motions are quenched in the crystal at cryogenic temperatures. The combination of NMR and X-ray crystallography in iterative refinement promises to provide an atomic resolution description of the alternate conformational substates that are sampled through picosecond to nanosecond timescale fluctuations of the protein structure. The method also provides insights into the structural heterogeneity of nonmethyl side chains, aromatic residues, and ligands, which are less commonly analyzed by NMR relaxation measurements.

Characterization of three pro-inflammatory cytokines, TNFα1, TNFα2 and IL-1β, in cage-reared Atlantic bluefin tuna Thunnus thynnus

Atlantic bluefin tuna (BFT) (Thunnus thynnus) is of great economic significance for world aquaculture and therefore it is necessary to ensure optimal and sustainable conditions for the farming of this species. Intensive culture of fish may be limited by infectious diseases that can impact on growth performance and cause heavy losses. However, to date there are no reports of cloning and expression analysis of any major immune genes of Atlantic BFT although some immune genes are known in other BFT species. Therefore the aim of this study was to characterize the first cytokine molecules in Atlantic BFT, through: 1) Isolation of full-length cDNA and gene sequences of TNFα1, TNFα2 and IL-1β, 2) comparison of these molecules to known sequences in other vertebrates, especially teleost fish, by multiple sequence alignment, phylogenetic tree analysis and homology modeling; 3) Quantification of in vivo expression of these cytokines in selected tissues in reared BFT over the duration of the farming process. The results indicated that these three cytokines could have value for monitoring Atlantic BFT health status. Curiously, the liver seemed to be an important site of cytokine production during poor health conditions in this species, perhaps reflecting its role as an important organ involved in fish defenses.

The structure of TAX1BP1 UBZ1+2 provides insight into target specificity and adaptability

TAX1BP1 is a novel ubiquitin-binding adaptor protein involved in the negative regulation of the NF-kappaB transcription factor, which is a key player in inflammatory responses, immunity and tumorigenesis. TAX1BP1 recruits A20 to the ubiquitinated signaling proteins TRAF6 and RIP1, leading to their A20-mediated deubiquitination and the disruption of IL-1-induced and TNF-induced NF-kappaB signaling, respectively. The two zinc fingers localized at its C-terminus function as novel ubiquitin-binding domains (UBZ, ubiquitin-binding zinc finger). Here we present for the first time both the solution and crystal structures of two classical UBZ domains in tandem within the human TAX1BP1. The relative orientation of the two domains is slightly different in the X-ray structure with respect to the NMR structure, indicating some degree of conformational flexibility, which is rationalized by NMR relaxation data. The observed degree of flexibility and stability between the two UBZ domains might have consequences on the recognition mechanism of interacting partners.

STAP Refinement of the NMR database: a database of 2405 refined solution NMR structures

According to several studies, some nuclear magnetic resonance (NMR) structures are of lower quality, less reliable and less suitable for structural analysis than high-resolution X-ray crystallographic structures. We present a public database of 2405 refined NMR solution structures [statistical torsion angle potentials (STAP) refinement of the NMR database, http://psb.kobic.re.kr/STAP/refinement] from the Protein Data Bank (PDB). A simulated annealing protocol was employed to obtain refined structures with target potentials, including the newly developed STAP. The refined database was extensively analysed using various quality indicators from several assessment programs to determine the nuclear Overhauser effect (NOE) completeness, Ramachandran appearance, χ₁-χ₂ rotamer normality, various parameters for protein stability and other indicators. Most quality indicators are improved in our protocol mainly due to the inclusion of the newly developed knowledge-based potentials. This database can be used by the NMR structure community for further development of research and validation tools, structure-related studies and modelling in many fields of research.

High-resolution membrane protein structure by joint calculations with solid-state NMR and X-ray experimental data

X-ray diffraction and nuclear magnetic resonance spectroscopy (NMR) are the staple methods for revealing atomic structures of proteins. Since crystals of biomolecular assemblies and membrane proteins often diffract weakly and such large systems encroach upon the molecular tumbling limit of solution NMR, new methods are essential to extend structures of such systems to high resolution. Here we present a method that incorporates solid-state NMR restraints alongside of X-ray reflections to the conventional model building and refinement steps of structure calculations. Using the 3.7 Å crystal structure of the integral membrane protein complex DsbB-DsbA as a test case yielded a significantly improved backbone precision of 0.92 Å in the transmembrane region, a 58% enhancement from using X-ray reflections alone. Furthermore, addition of solid-state NMR restraints greatly improved the overall quality of the structure by promoting 22% of DsbB transmembrane residues into the most favored regions of Ramachandran space in comparison to the crystal structure. This method is widely applicable to any protein system where X-ray data are available, and is particularly useful for the study of weakly diffracting crystals.

Consensus structure of Pf1 filamentous bacteriophage from X-ray fibre diffraction and solid-state NMR

Filamentous bacteriophages (filamentous bacterial viruses or Inovirus) are simple and well-characterised macromolecular assemblies that are widely used in molecular biology and biophysics, both as paradigms for studying basic biological questions and as practical tools in areas as diverse as immunology and solid-state physics. The strains fd, M13 and f1 are virtually identical filamentous phages that infect bacteria expressing F-pili, and are sometimes grouped as the Ff phages. For historical reasons fd has often been used for structural studies, but M13 and f1 are more often used for biological experiments. Many other strains have been identified that are genetically quite distinct from Ff and yet have a similar molecular structure and life cycle. One of these, Pf1, gives the highest resolution X-ray fibre diffraction patterns known for filamentous bacteriophage. These diffraction patterns have been used in the past to derive a molecular model for the structure of the phage. Solid-state NMR experiments have been used in separate studies to derive a significantly different model of Pf1. Here we combine previously published X-ray fibre diffraction data and solid-state NMR data to give a consensus structure model for Pf1 filamentous bacteriophage, and we discuss the implications of this model for assembly of the phage at the bacterial membrane.

Water in the polar and nonpolar cavities of the protein interleukin-1β

Water in the protein interior serves important structural and functional roles and is also increasingly recognized as a relevant factor in drug binding. The nonpolar cavity in the protein interleukin-1β has been reported to be filled by water on the basis of some experiments and simulations and to be empty on the basis of others. Here we study the thermodynamics of filling the central nonpolar cavity and the four polar cavities of interleukin-1β by molecular dynamics simulation. We use different water models (TIP3P and SPC/E) and protein force fields (amber94 and amber03) to calculate the semigrand partition functions term by term that quantify the hydration equilibria. We consistently find that water in the central nonpolar cavity is thermodynamically unstable, independent of force field and water model. The apparent reason is the relatively small size of the cavity, with a volume less than ∼80 Å(3). Our results are consistent with the most recent X-ray crystallographic and simulation studies but disagree with an earlier interpretation of nuclear magnetic resonance (NMR) experiments probing protein-water interactions. We show that, at least semiquantitatively, the measured nuclear Overhauser effects indicating the proximity of water to the methyl groups lining the nonpolar cavity can, in all likelihood, be attributed to interactions with buried and surface water molecules near the cavity. The same methods applied to determine the occupancy of the polar cavities show that they are filled by the same number of water molecules observed in crystallography, thereby validating the theoretical and simulation methods used to study the water occupancy in the nonpolar protein cavity.

Molecular dynamics free energy calculations to assess the possibility of water existence in protein nonpolar cavities

Are protein nonpolar cavities filled with water molecules? Although many experimental and theoretical investigations have been done, particularly for the nonpolar cavity of IL-1 beta, the results are still conflicting. To study this problem from the thermodynamic point of view, we calculated hydration free energies of four protein nonpolar cavities by means of the molecular dynamics thermodynamic integration method. In addition to the IL-1 beta cavity (69 A(3)), we selected the three largest nonpolar cavities of AvrPphB (81 A(3)), Trp repressor (87 A(3)), and hemoglobin (108 A(3)) from the structural database, in view of the simulation result from another study that showed larger nonpolar cavities are more likely to be hydrated. The calculations were performed with flexible and rigid protein models. The calculated free energy changes were all positive; hydration of the nonpolar cavities was energetically unfavorable for all four cases. Because hydration of smaller cavities should happen more rarely, we conclude that existing protein nonpolar cavities are not likely to be hydrated. Although a possibility remains for much larger nonpolar cavities, such cases are not found experimentally. We present a hypothesis to explain this: hydrated nonpolar cavities are quite unstable and the conformation could not be maintained.

Systematic comparison of crystal and NMR protein structures deposited in the protein data bank

Nearly all the macromolecular three-dimensional structures deposited in Protein Data Bank were determined by either crystallographic (X-ray) or Nuclear Magnetic Resonance (NMR) spectroscopic methods. This paper reports a systematic comparison of the crystallographic and NMR results deposited in the files of the Protein Data Bank, in order to find out to which extent these information can be aggregated in bioinformatics. A non-redundant data set containing 109 NMR – X-ray structure pairs of nearly identical proteins was derived from the Protein Data Bank. A series of comparisons were performed by focusing the attention towards both global features and local details. It was observed that: (1) the RMDS values between NMR and crystal structures range from about 1.5 Å to about 2.5 Å; (2) the correlation between conformational deviations and residue type reveals that hydrophobic amino acids are more similar in crystal and NMR structures than hydrophilic amino acids; (3) the correlation between solvent accessibility of the residues and their conformational variability in solid state and in solution is relatively modest (correlation coefficient = 0.462); (4) beta strands on average match better between NMR and crystal structures than helices and loops; (5) conformational differences between loops are independent of crystal packing interactions in the solid state; (6) very seldom, side chains buried in the protein interior are observed to adopt different orientations in the solid state and in solution.

A refinement protocol to determine structure, topology, and depth of insertion of membrane proteins using hybrid solution and solid-state NMR restraints

To fully describe the fold space and ultimately the biological function of membrane proteins, it is necessary to determine the specific interactions of the protein with the membrane. This property of membrane proteins that we refer to as structural topology cannot be resolved using X-ray crystallography or solution NMR alone. In this article, we incorporate into XPLOR-NIH a hybrid objective function for membrane protein structure determination that utilizes solution and solid-state NMR restraints, simultaneously defining structure, topology, and depth of insertion. Distance and angular restraints obtained from solution NMR of membrane proteins solubilized in detergent micelles are combined with backbone orientational restraints (chemical shift anisotropy and dipolar couplings) derived from solid-state NMR in aligned lipid bilayers. In addition, a supplementary knowledge-based potential, E (z) (insertion depth potential), is used to ensure the correct positioning of secondary structural elements with respect to a virtual membrane. The hybrid objective function is minimized using a simulated annealing protocol implemented into XPLOR-NIH software for general use.

Synergy of NMR, computation, and X-ray crystallography for structural biology

NMR spectroscopy and X-ray crystallography are currently the two most widely applied methods for the determination of macromolecular structures at high resolution. More recently, significant advances have been made in algorithms for the de novo prediction of protein structure, and, in favorable cases, the predicted models agree extremely well with experimentally determined structures. Here, we demonstrate a synergistic combination of NMR spectroscopy, de novo structure prediction, and X-ray crystallography in an effective overall strategy for rapidly determining the structure of the coat protein C-terminal domain from the Sulfolobus islandicus rod-shaped virus (SIRV). This approach takes advantage of the most accessible aspects of each structural technique and may be widely applicable for structure determination.

Improving NMR protein structure quality by Rosetta refinement: a molecular replacement study

The structure of human protein HSPC034 has been determined by both solution nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography. Refinement of the NMR structure ensemble, using a Rosetta protocol in the absence of NMR restraints, resulted in significant improvements not only in structure quality, but also in molecular replacement (MR) performance with the raw X-ray diffraction data using MOLREP and Phaser. This method has recently been shown to be generally applicable with improved MR performance demonstrated for eight NMR structures refined using Rosetta (Qian et al., Nature 2007;450:259-264). Additionally, NMR structures of HSPC034 calculated by standard methods that include NMR restraints have improvements in the RMSD to the crystal structure and MR performance in the order DYANA, CYANA, XPLOR-NIH, and CNS with explicit water refinement (CNSw). Further Rosetta refinement of the CNSw structures, perhaps due to more thorough conformational sampling and/or a superior force field, was capable of finding alternative low energy protein conformations that were equally consistent with the NMR data according to the Recall, Precision, and F-measure (RPF) scores. On further examination, the additional MR-performance shortfall for NMR refined structures as compared with the X-ray structure were attributed, in part, to crystal-packing effects, real structural differences, and inferior hydrogen bonding in the NMR structures. A good correlation between a decrease in the number of buried unsatisfied hydrogen-bond donors and improved MR performance demonstrates the importance of hydrogen-bond terms in the force field for improving NMR structures. The superior hydrogen-bond network in Rosetta-refined structures demonstrates that correct identification of hydrogen bonds should be a critical goal of NMR structure refinement. Inclusion of nonbivalent hydrogen bonds identified from Rosetta structures as additional restraints in the structure calculation results in NMR structures with improved MR performance.

Production and comprehensive quality control of recombinant human Interleukin-1beta: a case study for a process development strategy

We describe an efficient strategy to produce high-quality proteins by using a single large IMAC chromatography column and enzymatic His-tag removal via the TAGZyme system in pilot scale. Numerous quality assays demonstrated a high purity of the final product, the human cytokine Interleukin-1beta (IL-1beta). The protein preparation was apparently free of host cell proteins, endotoxins, protease, and aggregates. The N-terminal amino acid sequence of IL-1beta was in full agreement with the natural mature form of IL-1beta. The homogeneity of the product was further shown by X-ray structure determination which confirmed the previously solved structure of the protein. We propose the applied workflow as a strategy for industrial production of protein-based biopharmaceuticals.

Structural differences between the putative carbohydrate-recognition domains of human IL-1 alpha, IL-1 beta and IL-1 receptor antagonist obtained by in silico modeling

In a previous report (Cebo et al. J Biol Chem 276 (2001) 5685-5691), it was established that biologically active recombinant human IL-1alpha and IL-1beta had different carbohydrate-binding properties. IL-1alpha recognized a di-antennary N-glycan with two alpha2-3-linked sialic acid residues, whereas IL-1beta recognized the GM(4), a alpha2-3-linked sialylated glycosphingolipid. These different carbohydrate-binding properties of two interleukins binding to the same receptor (IL-1R) could explain why these molecules had different biological effects and cell specificities. Molecular modeling of the ligands and in silico docking experiments defined putative carbohydrate-recognition domains localized in the same area of the two molecules, a domain different from that defined as the type I IL-1R binding domain. The calculated pattern of hydrogen bonding and of van der Waals interactions fulfilled the essential features observed for calcium-independent lectins (mammalian, viral or bacterial). The analysis of the same domain of the third members of this family of molecules, the IL-1R-antagonist, indicated it did not fulfill the criteria for carbohydrate-recognition domains. It is proposed that its role as a pure antagonist is due to the absence of lectin activity and consequently explained its inability to associate IL-1R with other surface molecular complexes necessary for signaling.

A general method for the unbiased improvement of solution NMR structures by the use of related X-ray data, the AUREMOL-ISIC algorithm

Background

Rapid and accurate three-dimensional structure determination of biological macromolecules is mandatory to keep up with the vast progress made in the identification of primary sequence information. During the last few years the amount of data deposited in the protein data bank has substantially increased providing additional information for novel structure determination projects. The key question is how to combine the available database information with the experimental data of the current project ensuring that only relevant information is used and a correct structural bias is produced. For this purpose a novel fully automated algorithm based on Bayesian reasoning has been developed. It allows the combination of structural information from different sources in a consistent way to obtain high quality structures with a limited set of experimental data. The new ISIC (Intelligent Structural Information Combination) algorithm is part of the larger AUREMOL software package.

Results

Our new approach was successfully tested on the improvement of the solution NMR structures of the Ras-binding domain of Byr2 from Schizosaccharomyces pombe, the Ras-binding domain of RalGDS from human calculated from a limited set of NMR data, and the immunoglobulin binding domain from protein G from Streptococcus by their corresponding X-ray structures. In all test cases clearly improved structures were obtained. The largest danger in using data from other sources is a possible bias towards the added structure. In the worst case instead of a refined target structure the structure from the additional source is essentially reproduced. We could clearly show that the ISIC algorithm treats these difficulties properly.

Conclusion

In summary, we present a novel fully automated method to combine strongly coupled knowledge from different sources. The combination with validation tools such as the calculation of NMR R-factors strengthens the impact of the method considerably since the improvement of the structures can be assessed quantitatively. The ISIC method can be applied to a large number of similar problems where the quality of the obtained three-dimensional structures is limited by the available experimental data like the improvement of large NMR structures calculated from sparse experimental data or the refinement of low resolution X-ray structures. Also structures may be refined using other available structural information such as homology models.

The Xplor-NIH NMR molecular structure determination package

We announce the availability of the Xplor-NIH software package for NMR biomolecular structure determination. This package consists of the pre-existing XPLOR program, along with many NMR-specific extensions developed at the NIH. In addition to many features which have been developed over the last 20 years, the Xplor-NIH package contains an interface with a new programmatic framework written in C++. This interface currently supports the general purpose scripting languages Python and TCL, enabling rapid development of new tools, such as new potential energy terms and new optimization methods. Support for these scripting languages also facilitates interaction with existing external programs for structure analysis, structure manipulation, visualization, and spectral analysis.

A single amino acid difference between human and monkey interleukin (IL)-1beta dictates effective binding to soluble type II IL-1 receptor

Soluble type II interleukin (IL)-1 receptor (sIL1R-II) binds human IL-1beta with high affinity and neutralizes its activity. Recombinant sIL1R-II is considered a potentially useful anti-IL-1 therapeutic, and preclinical studies have been undertaken with this molecule in primates. To better understand the cytokine-receptor interactions occurring in this nonhuman context, monkey IL-1 and IL1R-II were cloned, and their binding abilities were examined in vitro. IL-1beta from cynomolgus monkey was capable of binding and activating the human type I IL-1 receptor. However, unlike human IL-1beta, it was unable to effectively bind and become neutralized by sIL1R-II. Human and cynomolgus IL-1beta proteins are 96% identical, differing by only six amino acids. Structural and mutational analysis revealed that the unique sIL1R-II binding ability of human IL-1beta is due to a single amino acid difference compared with monkey IL-1beta.

Ribosomal protein L9: a structure determination by the combined use of X-ray crystallography and NMR spectroscopy

The structure of protein L9 from the Bacillus stearothernophilus ribosome has been determined at 2.5 A resolution by refinement against single crystal X-ray diffraction data with additional constraints provided by NMR data. This highly elongated protein consists of two domains separated by a nine-turn connecting helix. Conserved aromatic and positively charged amino acid residues on the surface of each domain are likely to be directly involved in binding 23 S ribosomal RNA. The shape of the protein, with its two widely spaced RNA-binding sites, suggests that it may serve as a "molecular strut", most likely playing a role in ribosome assembly and/or maintaining the catalytically active conformation of the ribosomal RNA. The combined use of X-ray and NMR data in the refinement procedure was essential in defining the N-terminal domain of the protein, which was relatively poorly determined by the X-ray data alone. In addition to resolving the ambiguities in defining the hydrophobic core and side-chain conformations with the N-terminal domain, this combined NMR-X-ray analysis provides the first detailed and accurate view of the N-terminal RNA-binding site. NMR data also showed that the N-terminal domain is stable in solution, indicated by amide protons that are protected from solvent exchange. The lack of definition of the N-terminal domain in the X-ray structure is therefore likely due to packing disorder within the crystal rather than structural instability. This combined NMR-X-ray analysis provides a useful model as to how X-ray and NMR data can be practically and logically combined in the determination of the structure of a single protein molecule.

The oligomerization domain of p53: crystal structure of the trigonal form

The structure of the oligomerization domain of the p53 tumor suppressor protein was determined in the trigonal crystal form, using a refined NMR structure as a model. A synthetic peptide comprising residues 319-360 of human p53 crystallized in the space group P3(1)21. There is one biologically relevant tetrameric domain in the crystallographic asymmetric unit. The structure was refined jointly with NMR data, only the third such case (the previous examples being IL-1beta (Shaanan, B., Gronenborn, A.M., Cohen, G.H., Gilliland, G.L., Veerapandian, B., Davies, D.R. and Clore, G.M. (1992) Science 257, 961-964 [1]) and BPTI (Schiffer, C., Huber, R., Wuthrich, K. and Van Gunsteren, W.F. (1994) J. Mol. Biol. 241, 588-599 [21)), to 2.5 A resolution with an R factor of 0.207. The distribution of tumor-derived mutations in the oligomerization region together with structural and biological data suggest a strategy for the design of antitumor therapeutics.

Formation of an active dimer during storage of interleukin-1 receptor antagonist in aqueous solution

The degradation products of recombinant human interleukin-1 receptor antagonist (rhIL-1ra) formed during storage at 30 degrees C in aqueous solution were characterized. Cationic exchange chromatography of the stored sample showed two major, new peaks eluting before (P1) and after (L2) the native protein, which were interconvertible. Size-exclusion chromatography and electrophoresis documented that both the P1 and L2 fractions were irreversible dimers, formed by noncovalent interactions. A competition assay with interleukin-1 indicated that on a per monomer basis the P1 and L2 dimers retained about two-thirds of the activity of the native monomer. Infrared and far-UV circular dichroism spectroscopies showed that only minor alterations in secondary structure arose upon the formation of the P1 dimer. However, alteration in the near-UV circular dichroism spectrum suggested the presence of disulfide bonds in the P1 dimer, which are absent in the native protein. Mass spectroscopy and tryptic mapping, before and after carboxymethylation, demonstrated that the P1 dimer contained an intramolecular disulfide bond between Cys-66 and Cys-69. Although conversion of native protein to the P1 dimer was irreversible in buffer alone, the native monomer could be regained by denaturing the P1 dimer with guanidine hydrochloride and renaturing it by dialysis, suggesting that the intramolecular disulfide bond does not interfere with refolding. Analysis of the time course of P1 formation during storage at 30 degrees C indicated that the process followed first-order, and not second-order, kinetics, suggesting that the rate-limiting step was not dimerization. It is proposed that a conformational change in the monomer is the rate-limiting step in the formation of the P1 dimer degradation product. Sucrose stabilized the native monomer against this process. This result can be explained by the general stabilization mechanism for this additive, which is due to its preferential exclusion from the protein surface.

A comparison of X-ray and NMR structures for human endothelin-1

Direct comparisons between the recently solved X-ray and NMR structures of human endothelin-1 with respect to secondary structure, RMS deviations, surface accessibilities, and side-chain conformers indicate important differences in conformation, especially in the C-terminus, but also in the central loop region, that are important for defining the specificity of binding. These differences are larger than seen for other X-ray and NMR structures that have been compared. Comparisons between the X-ray structure and the NMR NOE constraints highlight the regions of flexibility and environment-induced diversity in the endothelin structures.

Comparison of four independently determined structures of human recombinant interleukin-4

Four independent structures of human interleukin-4, two determined by nuclear magnetic resonance techniques and two by X-ray diffraction, have been compared in detail. The core of this four helix bundle protein is very similar in all the structures but there are some differences in loop regions that are known to be mobile in solution. Careful comparison of the experimental data sets and the methods of analysis of the different laboratories has provided clues to the sources of most of the differences, and also answered some general questions about the accuracy of protein structure determination by these two techniques.

Young Investigator Award Lecture. Structures of larger proteins, protein-ligand and protein-DNA complexes by multidimensional heteronuclear NMR

The recent development of a whole panoply of multidimensional heteronuclear-edited and -filtered NMR experiments has revolutionized the field of protein structure determination by NMR, making it possible to extend the methodology from the 10-kDa limit of conventional 2-dimensional NMR to systems up to potentially 35-40 kDa. The basic strategy for solving 3-dimensional structures of larger proteins and protein-ligand complexes in solution using 3- and 4-dimensional NMR spectroscopy is summarized, and the power of these methods is illustrated using 3 examples: interleukin-1 beta, the complex of calmodulin with a target peptide, and the specific complex of the transcription factor GATA-1 with its cognate DNA target site.

Hematopoietic cytokines: similarities and differences in the structures, with implications for receptor binding

Crystal and NMR structures of helical cytokines--interleukin-4 (IL-4), granulocyte-macrophage colony-stimulating factor (GM-CSF), and interleukin-2 (IL-2)--have been compared. Root mean square deviations in the C alpha coordinates for the conserved regions of the helices were 1-2 A between different cytokines, about twice the differences observed for independently determined crystal and solution structures of IL-4. Considerable similarity in amino acid sequence in the areas expected to interact with the receptors was detected, and the available mutagenesis data for these cytokines were correlated with structure conservation. Models of cytokine-receptor interactions were postulated for IL-4 based on its structure as well as on the published structure of human growth hormone interacting with its receptors (de Vos, A.M., Ultsch, M., & Kossiakoff, A.A., 1992, Science 255, 306-312). Patches of positively charged residues on the surfaces of helices C and D of IL-4 may be responsible for the interactions with the negatively charged residues found in the complementary parts of the IL-4 receptors.

Kinetics of folding of the all-beta sheet protein interleukin-1 beta

The folding of the all-beta sheet protein, interleukin-1 beta, was studied with nuclear magnetic resonance (NMR) spectroscopy, circular dichroism, and fluorescence. Ninety percent of the beta structure present in the native protein, as monitored by far-ultraviolet circular dichroism, was attained within 25 milliseconds, correlating with the first kinetic phase determined by tryptophan and 1-anilinonaphthalene-8-sulfonate fluorescence. In contrast, formation of stable native secondary structure, as measured by quenched-flow deuterium-hydrogen exchange experiments, began after only 1 second. Results from the NMR experiments indicated the formation of at least two intermediates with half-lives of 0.7 to 1.5 and 15 to 25 seconds. The final stabilization of the secondary structure, however, occurs on a time scale much greater than 25 seconds. These results differ from previous results on mixed alpha helix-beta sheet proteins in which both the alpha helices and beta sheets were stabilized very rapidly (less than 10 to 20 milliseconds).

Computational challenges for macromolecular structure determination by X-ray crystallography and solution NMR-spectroscopy

Macromolecular structure determination by X-ray crystallography and solution NMR spectroscopy has experienced unprecedented growth during the past decade.