Backbone dynamics measurements on leukemia inhibitory factor, a rigid four-helical bundle cytokine

The structure of the extracellular domains of human interleukin 11α receptor reveals mechanisms of cytokine engagement

Interleukin (IL) 11 activates multiple intracellular signaling pathways by forming a complex with its cell surface α-receptor, IL-11Rα, and the β-subunit receptor, gp130. Dysregulated IL-11 signaling has been implicated in several diseases, including some cancers and fibrosis. Mutations in IL-11Rα that reduce signaling are also associated with hereditary cranial malformations. Here we present the first crystal structure of the extracellular domains of human IL-11Rα and a structure of human IL-11 that reveals previously unresolved detail. Disease-associated mutations in IL-11Rα are generally distal to putative ligand-binding sites. Molecular dynamics simulations showed that specific mutations destabilize IL-11Rα and may have indirect effects on the cytokine-binding region. We show that IL-11 and IL-11Rα form a 1:1 complex with nanomolar affinity and present a model of the complex. Our results suggest that the thermodynamic and structural mechanisms of complex formation between IL-11 and IL-11Rα differ substantially from those previously reported for similar cytokines. This work reveals key determinants of the engagement of IL-11 by IL-11Rα that may be exploited in the development of strategies to modulate formation of the IL-11-IL-11Rα complex.

GRPY: An Accurate Bead Method for Calculation of Hydrodynamic Properties of Rigid Biomacromolecules

Two main problems that arise in the context of hydrodynamic bead modeling are an inaccurate treatment of bead overlaps and the necessity of using volume corrections when calculating intrinsic viscosity. We present a formalism based on the generalized Rotne-Prager-Yamakawa approximation that successfully addresses both of these issues. The generalized Rotne-Prager-Yamakawa method is shown to be highly effective for the calculation of transport properties of rigid biomolecules represented as assemblies of spherical beads of different sizes, both overlapping and nonoverlapping. We test the method on simple molecular shapes as well as real protein structures and compare its performance with other computational approaches.

Molecular Dynamics Simulations Reveal Key Roles of the Interleukin-6 Alpha Receptor in the Assembly of the Human Interleukin-6 Receptor Complex

Human interleukin-6 (hIL-6) is a pleiotropic cytokine with three distinct receptor epitopes, termed sites I, II, and III, which function to assemble a signaling complex. hIL-6 signals via a glycoprotein 130 (gp130) homodimer after initially forming a heterodimer with the nonsignaling α-receptor (IL-6Rα). The molecular description of the assembly of the hIL-6 signaling complex remains elusive because available structures provide descriptions of hIL-6 in its free and fully bound receptor forms, but not for intermediate steps that are crucial in the stepwise assembly of the signaling complex. In this report, molecular dynamics simulations provide atomic details describing the functional role of the initial hIL-6/IL-6Rα complex in facilitating subsequent interactions with gp130, which have not been previously shown. IL-6Rα binding to hIL-6 rigidifies the flexible N-terminus of the hIL-6 AB-loop through interactions with the D2 domain of IL-6Rα. This rigidification combined with repositioning of residues involved in gp130 receptor recognition promotes gp130 binding at site III. Binding of gp130 receptors at sites II and III is coupled with the release of the hIL-6 N-terminal AB-loop interaction and a pivoting of IL-6Rα around the hIL-6 helix bundle to the state of the hIL-6/IL-6Rα/gp130 complex.

A comparative NMR study of the polypeptide backbone dynamics of hemoglobin in the deoxy and carbonmonoxy forms

Model-free-based NMR dynamics studies have been undertaken for polypeptide backbone amide N-H bond vectors for both the deoxy and carbonmonoxy forms of chain-specific, isotopically (15N and 2H) labeled tetrameric hemoglobin (Hb) using 15N-relaxation parameters [longitudinal relaxation rate (R1), transverse relaxation rate (R2), and heteronuclear nuclear Overhauser effect (NOE)] measured at two temperatures (29 and 34 degrees C) and two magnetic field strengths (11.7 and 14.1 T). In both deoxy and carbonmonoxy forms of human normal adult hemoglobin (Hb A), the amide N-H bonds of most amino acid residues are rigid on the fast time scale (nanosecond to picosecond), except for the loop regions and certain helix-helix connections. Although rigid in deoxy-Hb A, beta146His has been found to be free from restriction of its backbone motions in the CO form, presumably due to the rupture of its hydrogen bond/salt bridge network. We now have direct dynamics evidence for this structural transition of Hb in solution. While remarkably flexible in the deoxy state, alpha31Arg and beta123Thr, neighbors in the intradimer (alpha1beta1) interface, exhibit stiffening upon CO binding. These findings imply a role for alpha31Arg and beta123Thr in the intradimer communication but contradict the results from X-ray crystallography. We have also found that there is considerable flexibility in the intradimer (alpha1beta1) interface (i.e., B, G, and H helices and the GH corner) and possible involvement of several amino acid residues (e.g., alpha31Arg, beta3Leu, beta41Phe, beta123Thr, and beta146His) in the allosteric pathway. Several amino acid residues at the intradimer interfaces, such as beta109Val, appear to be involved in possible conformational exchange processes. The dynamic picture derived from the present study provides new insights into the traditional description of the stereochemical mechanism for the cooperative oxygenation of Hb A based on X-ray crystallographic results.

Dynamics of the SPRY domain-containing SOCS box protein 2: flexibility of key functional loops

The SPRY domain was identified originally as a sequence repeat in the dual-specificity kinase splA and ryanodine receptors and subsequently found in many other distinct proteins, including more than 70 encoded in the human genome. It is a subdomain of the B30.2/SPRY domain and is believed to function as a protein-protein interaction module. Three-dimensional structures of several B30.2/SPRY domain-containing proteins have been reported recently: murine SSB-2 in solution by NMR spectroscopy, a Drosophila SSB (GUSTAVUS), and human PRYSPRY protein by X-ray crystallography. The three structures share a core of two antiparallel beta-sheets for the B30.2/SPRY domain but show differences located mainly at one end of the beta-sandwich. Analysis of SSB-2 residues required for interactions with its intracellular ligands has provided insights into B30.2/SPRY binding specificity and identified loop residues critical for the function of this domain. We have investigated the backbone dynamics of SSB-2 by means of Modelfree analysis of its backbone (15)N relaxation parameters and carried out coarse-grained dynamics simulation of B30.2/SPRY domain-containing proteins using normal mode analysis. Translational self-diffusion coefficients of SSB-2 measured using pulsed field gradient NMR were used to confirm the monomeric state of SSB-2 in solution. These results, together with previously reported amide exchange data, highlight the underlying flexibility of the loop regions of B30.2/SPRY domain-containing proteins that have been shown to be important for protein-protein interactions. The underlying flexibility of certain regions of the B30.2/SPRY domain-containing proteins may also contribute to some apparent structural differences observed between GUSTAVUS or PRYSPRY and SSB-2.

The SPRY domain of SSB-2 adopts a novel fold that presents conserved Par-4-binding residues

The four mammalian SPRY domain-containing SOCS box proteins (SSB-1 to SSB-4) are characterized by a C-terminal SOCS box and a central SPRY domain. We have determined the first SPRY-domain structure, as part of SSB-2, by NMR. This domain adopts a novel fold consisting of a beta-sandwich structure formed by two four-stranded antiparallel beta-sheets with a unique topology. We demonstrate that SSB-1, SSB-2 and SSB-4, but not SSB-3, bind prostate apoptosis response protein-4 (Par-4). Mutational analysis of SSB-2 loop regions identified conserved structural determinants for its interaction with Par-4 and the hepatocyte growth factor receptor, c-Met. Mutations in analogous loop regions of pyrin and midline-1 SPRY domains have been shown to cause Mediterranean fever and Opitz syndrome, respectively. Our findings provide a template for SPRY-domain structure and an insight into the mechanism of SPRY-protein interaction.

Affinity maturation of leukemia inhibitory factor and conversion to potent antagonists of signaling

Leukemia inhibitory factor (LIF)-induced cell signaling occurs following sequential binding to the LIF receptor alpha-chain (LIFR), then to the gp130 co-receptor used by all members of the interleukin-6 family of cytokines. By monovalently displaying human LIF on the surface of M13 phage and randomizing clusters of residues in regions predicted to be important for human LIFR binding, we have identified mutations, which lead to significant increases in affinity for binding to LIFR. Six libraries were constructed in which regions of 4-6 amino acids were randomized then panned against LIFR. Mutations identified in three distinct clusters, residues 53-57, 102-103, and 150-155, gave rise to proteins with significantly increased affinity for binding to both human and mouse LIFR. Combining the mutations for each of these regions further increased the affinity, such that the best mutants bound to human LIFR with >1000-fold higher affinity than wild-type human LIF. NMR analysis indicated that the mutations did not alter the overall structure of the molecule relative to the native protein, although some local changes occurred in the vicinity of the substituted residues. Despite increases in LIFR binding affinity, these mutants did not show any increase in activity as agonists of LIF-induced proliferation of Ba/F3 cells expressing human LIFR and gp130 compared with wild-type LIF. Incorporation of two additional mutations (Q29A and G124R), which were found to abrogate cell signaling, led to the generation of highly potent antagonists of both human and murine LIF-induced bioactivity.

NMR studies of the backbone flexibility and structure of human growth hormone: a comparison of high and low pH conformations

(15)N NMR relaxation parameters and amide (1)H/(2)H-exchange rates have been used to characterize the structural flexibility of human growth hormone (rhGH) at neutral and acidic pH. Our results show that the rigidity of the molecule is strongly affected by the solution conditions. At pH 7.0 the backbone dynamics parameters of rhGH are uniform along the polypeptide chain and their values are similar to those of other folded proteins. In contrast, at pH 2.7 the overall backbone flexibility increases substantially compared to neutral pH and the average order parameter approaches the lower limit expected for a folded protein. However, a significant variation of the backbone dynamics through the molecule indicates that under acidic conditions the mobility of the residues becomes more dependent on their location within the secondary structure units. In particular, the order parameters of certain loop regions decrease dramatically and become comparable to those found in unfolded proteins. Furthermore, the HN-exchange rates at low pH reveal that the residues most protected from exchange are clustered at one end of the helical bundle, forming a stable nucleus. We suggest that this nucleus maintains the overall fold of the protein under destabilizing conditions. We therefore conclude that the acid state of rhGH consists of a structurally conserved, but dynamically more flexible helical core surrounded by an aura of highly mobile, unstructured loops. However, in spite of its prominent flexibility the acid state of rhGH cannot be considered a "molten globule" state because of its high stability. It appears from our work that under certain conditions, a protein can tolerate a considerable increase in flexibility of its backbone, along with an increased penetration of water into its core, while still maintaining a stable folded conformation.

Solution structure of the eukaryotic pore-forming cytolysin equinatoxin II: implications for pore formation

Sea anemones produce a family of 18-20 kDa proteins, the actinoporins, that lyse cells by forming pores in cell membranes. Sphingomyelin plays an important role in their lytic activity, with membranes lacking this lipid being largely refractory to these toxins. The structure of the actinoporin equinatoxin II in aqueous solution, determined from NMR data, consists of two short helices packed against opposite faces of a beta-sandwich structure formed by two five-stranded beta-sheets. The protein core has extensive hydrophobic interfaces formed by residues projecting from the internal faces of the two beta-sheets. 15N relaxation data show uniform backbone dynamics, implying that equinatoxin II in solution is relatively rigid, except at the N terminus; its inferred rotational correlation time is consistent with values for monomeric proteins of similar mass. Backbone amide exchange rate data also support the view of a stable structure, even though equinatoxin II lacks disulfide bonds. As monitored by NMR, it unfolds at around 70 degrees C at pH 5.5. At 25 degrees C the structure is stable over the pH range 2.5-7.3 but below pH 2.5 it undergoes a slow transition to an incompletely unfolded structure resembling a molten globule. Equinatoxin II has two significant patches of positive electrostatic potential formed by surface-exposed Lys and Arg residues, which may assist its interaction with charged regions of the lipid head groups. Tyr and Trp residues on the surface may also contribute by interacting with the carbonyl groups of the acyl chains of target membranes. Data from mutational studies and truncated analogues identify two regions of the protein involved in membrane interactions, the N-terminal helix and the Trp-rich region. Once the protein is anchored, the N-terminal helix may penetrate the membrane, with up to four helices lining the pore, although other mechanisms of pore formation cannot be ruled out.

A unified picture of protein hydration: prediction of hydrodynamic properties from known structures

Hydration is essential for the structural and functional integrity of globular proteins. How much hydration water is required for that integrity? A number of techniques such as X-ray diffraction, nuclear magnetic resonance (NMR) spectroscopy, calorimetry, infrared spectroscopy, and molecular dynamics (MD) simulations indicate that the hydration level is 0.3-0.5 g of water per gram of protein for medium sized proteins. Hydrodynamic properties, when accounted for by modeling proteins as ellipsoids, appear to give a wide range of hydration levels. In this paper we describe an alternative numerical technique for hydrodynamic calculations that takes account of the detailed protein structures. This is made possible by relating hydrodynamic properties (translational and rotational diffusion constants and intrinsic viscosity) to electrostatic properties (capacitance and polarizability). We show that the use of detailed protein structures in predicting hydrodynamic properties leads to hydration levels in agreement with other techniques. A unified picture of protein hydration emerges. There are preferred hydration sites around a protein surface. These sites are occupied nearly all the time, but by different water molecules at different times. Thus, though a given water molecule may have a very short residence time (approximately 100-500 ps from NMR spectroscopy and MD simulations) in a particular site, the site appears fully occupied in experiments in which time-averaged properties are measured.