Characterization of proteins with wide-angle X-ray solution scattering (WAXS)

X-ray solution scattering in both the small-angle (SAXS) and wide-angle (WAXS) regimes is making an increasing impact on our understanding of biomolecular complexes. The accurate calculation of WAXS patterns from atomic coordinates has positioned the approach for rapid growth and integration with existing Structural Genomics efforts. WAXS data are sensitive to small structural changes in proteins; useful for calculation of the pair-distribution function at relatively high resolution; provides a means to characterize the breadth of the structural ensemble in solution; and can be used to identify proteins with similar folds. WAXS data are often used to test structural models, identify structural similarities and characterize structural changes. WAXS is highly complementary to crystallography and NMR. It holds great potential for the testing of structural models of proteins; identification of proteins that may exhibit novel folds; characterization of unfolded or natively disordered proteins; and detection of structural changes associated with protein function.     

X-ray solution scattering studies of the structural diversity intrinsic to protein ensembles

It is becoming increasingly clear that characterization of the protein ensemble-the collection of all conformations of which the protein is capable-will be a critical step in developing a full understanding of the linkage between structure, dynamics, and function. X-ray solution scattering in the small angle (SAXS) and wide-angle (WAXS) regimes represents an important new window to exploring the behavior of ensembles. The characteristics of the ensemble express themselves in X-ray solution scattering data in predictable ways. Here we present an overview of the effect that structural diversity intrinsic to protein ensembles has on scattering data. We then demonstrate the observation of these effects in scattering from four molecular systems; myoglobin; ubiquitin; alcohol dehydrogenase; and HIV protease; and demonstrate the modulation of these ensembles by ligand binding, mutation, and environmental factors. The observations are analyzed quantitatively in terms of the average spatial extent of structural fluctuations occurring within these proteins under different experimental conditions. The insights which these analyses support are discussed in terms of the function of the various proteins.     

Conformational Changes of Three Periplasmic Receptors for Bacterial Chemotaxis and Transport: The Maltose-, Glucose/Galactose- and Ribose-binding Proteins.

Small-angle X-ray scattering experiments were carried out for the maltose-, glucose/galactose- and ribose-binding proteins of Gram negative bacteria. All were shown to be monomers that decrease in radius of gyration on ligand binding.The results obtained for the maltose-binding protein agree well with crystal structures of the closed, ligand-bound, and open, ligand-free protein, suggesting that these are indeed the primary forms in solution. The closed form is stabilized by protein – sugar interactions, while the open conformation is stabilized by close contacts between the two domains. Since it is the proper spacial relationship of the domains in the closed form that is most important for interaction with chemotaxis and transport partners, the stabilization of the open form would help keep ligand-free molecules from interfering in function.The scattering results also provide evidence that a large conformational change takes place in association with ligand binding to the glucose/galactose- and ribose-binding proteins, and that the two changes are similar. Modeling suggests that the open forms resemble those found in the related leucine and leucine/isoleucine/valine-binding proteins, but are different from those observed for the maltose-binding protein and the related purine repressor.     

Subdividing repressor function: DNA binding affinity, selectivity, and allostery can be altered by amino acid substitution of nonconserved residues in a LacI/GalR homologue

Many mutations that impact protein function occur at residues that do not directly contact ligand. To understand the functional contributions from the sequence that links the DNA-binding and regulatory domains of the LacI/GalR homologues, we have created a chimeric protein (LLhP), which comprises the LacI DNA-binding domain, the LacI linker, and the PurR regulatory domain. Although DNA binding site residues are identical in LLhP and LacI, thermodynamic measurements of DNA binding affinity show that LLhP does not discriminate between alternative DNA ligands as well as LacI. In addition, small-angle scattering experiments show that LLhP is more compact than LacI. When DNA is released, LacI shows a 20 A increase in length that was previously attributed to unfolding of the linker. This change is not seen in apo-LLhP, even though the linker sequences of the two proteins are identical. Together, results indicate that long-range functional and structural changes are propagated across the interface that forms between the linker and regulatory domain. These changes could be mediated via the side chains of several linker residues that contact the regulatory domains of the naturally occurring proteins, LacI and PurR. Substitution of these residues in LLhP leads to a range of functional effects. Four variants exhibit altered affinity for DNA, with no changes in selectivity or allosteric response. Another two result in proteins that bind operator DNA with very low affinity and no allosteric response, similar to LacI binding nonspecific DNA sequences. Two more substitutions simultaneously diminish affinity, enhance allostery, and profoundly alter DNA ligand selectivity. Thus, positions within the linker can be varied to modulate different aspects of repressor function.     

Characterization of protein fold by wide-angle X-ray solution scattering

Wide-angle X-ray solution scattering (WAXS) patterns contain substantial information about the three-dimensional structure of a protein. Although WAXS data have far less information than is required for determination of a full three-dimensional structure, the actual amount of information contained in a WAXS pattern has not been carefully quantified. Here we carry out an analysis of the amount of information that can be extracted from a WAXS pattern and demonstrate that it is adequate to estimate the secondary-structure content of a protein and to strongly limit its possible tertiary structures. WAXS patterns computed from the atomic coordinates of a set of 498 protein domains representing all of known fold space were used as the basis for constructing a multidimensional space of all corresponding WAXS patterns (‘WAXS space’). Within WAXS space, each scattering pattern is represented by a single vector. A principal components analysis was carried out to identify those directions in WAXS space that provide the greatest discrimination among patterns. The number of dimensions that provide significant discrimination among protein folds agrees well with the number of independent parameters estimated from a naïve Shannon sampling theorem approach. Estimates of the relative abundances of secondary structures were made using training/test sets derived from this data set. The average error in the estimate of α-helical content was 11%, and of β-sheet content was 9%. The distribution of proteins that are members of the four structure classes, α, β, α/β and α+β, are well separated in WAXS space when data extending to a spacing of 2.2 Å are used. Quantification of the information embedded within a WAXS pattern indicates that these data can be used as a powerful constraint in homology modeling of protein structures.     

WAXS studies of the structural diversity of hemoglobin in solution

Specific ligation states of hemoglobin are, when crystallized, capable of taking on multiple quaternary structures. The relationship between these structures, captured in crystal lattices, and hemoglobin structure in solution remains uncertain. Wide-angle X-ray solution scattering (WAXS) is a sensitive probe of protein structure in solution that can distinguish among similar structures and has the potential to contribute to these issues. We used WAXS to assess the relationships among the structures of human and bovine hemoglobins in different liganded forms in solution. WAXS data readily distinguished among the various forms of hemoglobins. WAXS patterns confirm some of the relationships among hemoglobin structures that have been defined through crystallography and NMR and extend others. For instance, methemoglobin A in solution is, as expected, nearly indistinguishable from HbCO A. Interestingly, for bovine hemoglobin, the differences between deoxy-Hb, methemoglobin and HbCO are smaller than the corresponding differences in human hemoglobin. WAXS data were also used to assess the spatial extent of structural fluctuations of various hemoglobins in solution. Dynamics has been implicated in allosteric control of hemoglobin, and increased dynamics has been associated with lowered oxygen affinity. Consistent with that notion, WAXS patterns indicate that deoxy-Hb A exhibits substantially larger structural fluctuations than HbCO A. Comparisons between the observed WAXS patterns and those predicted on the basis of atomic coordinate sets suggest that the structures of Hb in different liganded forms exhibit clear differences from known crystal structures.     

EF loop conformational change triggers ligand binding in beta-lactoglobulins

Beta-lactoglobulins, belonging to the lipocalin family, are a widely studied group of proteins, characterized by the ability to solubilize and transport hydrophobic ligands, especially fatty acids. Despite many reports, the mechanism of ligand binding and the functional role of these proteins is still unclear, and many contradicting concepts are often encountered in the literature. In the present paper the comparative analysis of the binding properties of beta-lactoglobulins has been performed using sequence-derived information, structure-based electrostatic calculations, docking simulations, and NMR experiments. Our results reveal for the first time the mechanism of beta-lactoglobulin ligand binding, which is completely determined by the opening-closing of EF loop, triggered by Glu89 protonation. The alkaline shift observed for Glu89 pKa in porcine beta-lactoglobulin (pKa 9.7) with respect to the bovine species (pKa 5.5) depends upon the interplay of electrostatic effects of few nearby key residues. Porcine protein is therefore able to bind fatty acids provided that the appropriate pH solution conditions are met (pH > 8.6), where the EF loop conformational change can take place. The unusually high pH of binding detected for porcine beta-lactoglobulin seems to be functional to lipases activity. Theoretical pKa calculations extended to representative beta-lactoglobulins allowed the identification of key residues involved in structurally and functionally important electrostatic interactions. The results presented here provide a strong indication that the described conformational change is a common feature of all beta-lactoglobulins.     

The solution structure of the N-terminal domain of riboflavin synthase

The structure of the amino-terminal domain of Escherichia coli riboflavin synthase (RiSy) has been determined by NMR spectroscopy with riboflavin as a bound ligand. RiSy is functional as a 75 kDa homotrimer, each subunit of which consists of two domains which share very similar sequences and structures. The N-terminal domain (RiSy-N; 97 residues) forms a 20 kDa homodimer in solution which binds riboflavin with high affinity. The structure features a six-stranded antiparallel beta-barrel with a Greek-key fold, both ends of which are closed by an alpha-helix. One riboflavin molecule is bound per monomer in a site at one end of the barrel which is comprised of elements of both monomers. The structure and ligand binding are similar to that of the FAD binding domains of ferrodoxin reductase family proteins. The structure provides insights into the structure of the whole enzyme, the organisation of the functional trimer and the mechanism of riboflavin synthesis. C48 from the N-terminal domain is identified as the free cysteine implicated in a nucleophilic role in the synthesis mechanism, while H102 from the C-terminal domains is also likely to play a key role. Both are invariant in all known riboflavin synthase sequences.     

Validating Solution Ensembles from Molecular Dynamics Simulation by Wide-Angle X-ray Scattering Data

Wide-angle x-ray scattering (WAXS) experiments of biomolecules in solution have become increasingly popular because of technical advances in light sources and detectors. However, the structural interpretation of WAXS profiles is problematic, partly because accurate calculations of WAXS profiles from structural models have remained challenging. In this work, we present the calculation of WAXS profiles from explicit-solvent molecular dynamics (MD) simulations of five different proteins. Using only a single fitting parameter that accounts for experimental uncertainties because of the buffer subtraction and dark currents, we find excellent agreement to experimental profiles both at small and wide angles. Because explicit solvation eliminates free parameters associated with the solvation layer or the excluded solvent, which would require fitting to experimental data, we minimize the risk of overfitting. We further find that the influence from water models and protein force fields on calculated profiles are insignificant up to q ≈ 15 nm^?1. Using a series of simulations that allow increasing flexibility of the proteins, we show that incorporating thermal fluctuations into the calculations significantly improves agreement with experimental data, demonstrating the importance of protein dynamics in the interpretation of WAXS profiles. In addition, free MD simulations up to one microsecond suggest that the calculated profiles are highly sensitive with respect to minor conformational rearrangements of proteins, such as an increased flexibility of a loop or an increase of the radius of gyration by  <  1%. The present study suggests that quantitative comparison between MD simulations and experimental WAXS profiles emerges as an accurate tool to validate solution ensembles of biomolecules.     

The polyprotein and FAR lipid binding proteins of nematodes: shape and monomer/dimer states in ligand-free and bound forms

Nematodes produce two classes of small, helix-rich fatty acid- and retinol-binding proteins whose structures and in vivo functions remain to be elucidated. These are the polyprotein allergens (NPA) and the FAR proteins. The solution properties of recombinant forms of these proteins from parasitic [Ascaris suum (As) and Onchocerca volvulus (Ov)] and free-living [Caenorhabditis elegans (Ce)] nematodes have been examined. Analytical ultracentrifugation (AUC) showed that, contrary to previous findings, the rAs-NPA-1A polyprotein unit (approximately 15 kDa) is a monomer, and this stoichiometry is unaltered by ligand (oleic acid) binding. The rOv-FAR-1 and rCe-FAR-5 proteins differ in that the former forms a tight dimer and the latter a monomer, and these oligomeric states are also unaffected by ligand binding or protein concentration. Sedimentation equilibrium experiments showed that the partial specific volume v of the unliganded proteins agree well with the value calculated from amino acid composition extrapolated to experimental temperature, and was unaffected upon ligand binding. Data from small-angle X-ray scattering (SAXS) indicated that both of the monomeric proteins rAs-NPA-1A and rCe-FAR-5 are globular, although slightly elongated and flattened. These data are in good agreement with shapes predicted from sedimentation velocity experiments and hydrodynamic bead modelling. On the basis of functional and secondary structural homology with the ligand-binding domain of the retinoic acid receptor RXRalpha, de novo atomic resolution structures for rAs-NPA-1A and rCe-FAR-5 have been constructed which are consistent with the SAXS and hydrodynamic data.     

Solution structure and functional ligand screening of HI0719, a highly conserved protein from bacteria to humans in the YjgF/YER057c/UK114 family

HI0719 belongs to a large family of highly conserved proteins with no definitive molecular function and is found in organisms ranging from bacteria to humans. We describe the NMR structure of HI0719, the first solution structure for a member of this family. The overall fold is similar to the crystal structures of two homologues, YabJ from Bacillus subtilis and YjgF from Escherichia coli, and all three structures are similar to that of chorismate mutase, although there is little sequence homology and no apparent functional connection. HI0719 is a homotrimer with a distinct cavity located at the subunit interface. Six of the seven invariant residues in the high identity group of proteins are located in this cavity, suggesting that this may be a binding site for small molecules. Using previously published observations about the biological role of HI0719 family members as a guide, over 100 naturally occurring small molecules or structural analogues were screened for ligand binding using NMR spectroscopy. The targeted screening approach identified six compounds that bind to HI0719 at the putative active site. Five of these compounds are either alpha-keto acids or alpha,beta-unsaturated acids, while the sixth compound is structurally similar. Previous studies have proposed that some HI0719 homologues may act on small molecules in the isoleucine biosynthetic path and, if this is correct, the ligand screening results presented here suggest that the interaction most likely occurs with 2-ketobutyrate and/or its unstable enamine precursor.     

Small-angle X-ray scattering studies of enzymes

Enzyme function requires conformational changes to achieve substrate binding, domain rearrangements, and interactions with partner proteins, but these movements are difficult to observe. Small-angle X-ray scattering (SAXS) is a versatile structural technique that can probe such conformational changes under solution conditions that are physiologically relevant. Although it is generally considered a low-resolution structural technique, when used to study conformational changes as a function of time, ligand binding, or protein interactions, SAXS can provide rich insight into enzyme behavior, including subtle domain movements. In this perspective, we highlight recent uses of SAXS to probe structural enzyme changes upon ligand and partner-protein binding and discuss tools for signal deconvolution of complex protein solutions.     

Ligand binding to cytochrome c and other related haem proteins and peptides. Part III. Temperature dependence studies

The temperature dependency of ligand binding processes lend support to the proposed mechanisms and the factors affecting ligand binding reported earlier in this series. The free energy contribution from each factor affecting ligand binding was estimated for a number of haem proteins. The structures of the haem proteins used, as conveyed from ligand binding data, are in agreement with the structures of these haem proteins as determined by other methods (e.g. X-ray crystallography, NMR, etc.). Therefore, ligand binding could be used as a facile probe to investigate some of the structural and functional properties of haem proteins. In this respect, it was concluded that the structure of native cytochrome c at pH 10 is similar to the structure of carboxymethyl-Met 80 cytochrome c between pH 7 and 10.     

The role of FhuD2 in iron(III)-hydroxamate transport in Staphylococcus aureus. Demonstration that FhuD2 binds iron(III)-hydroxamates but with minimal conformational change and implication of mutations on transport

The fhuD2 gene encodes a lipoprotein that has previously been shown to be important for the utilization of iron(III)-hydroxamates by Staphylococcus aureus. We have studied the function of the FhuD2 protein in greater detail, and demonstrate here that the protein binds several iron(III)-hydroxamates. Mutagenesis of FhuD2 identified several residues that were important for the ability of the protein to function in iron(III)-hydroxamate transport. Several residues, notably Tyr-191, Trp-197, and Glu-202, were found to be critical for ligand binding. Moreover, mutation of two highly conserved glutamate residues, Glu-97 and Glu-231, had no affect on ligand binding, but did impair iron(III)-hydroxamate transport. Interestingly, the transport defect was not equivalent for all iron(III)-hydroxamates. We modeled FhuD2 against the high resolution structures of Escherichia coli FhuD and BtuF, two structurally related proteins, and showed that the three proteins share a similar overall structure. FhuD2 Glu-97 and Glu-231 were positioned on the surface of the N and C domains, respectively. Characterization of E97A, E231A, or E97A/E231A mutants suggests that these residues, along with the ligand itself, play a cumulative role in recognition by the ABC transporter FhuBGC2. In addition, small angle x-ray scattering was used to demonstrate that, in solution, FhuD2 does not undergo a detectable change in conformation upon binding iron(III)-hydroxamates. Therefore, the mechanism of binding and transport of ligands for binding proteins within this family is significantly different from that of other well studied binding protein families, such as that represented by maltose-binding protein.     

Non-catalytic Binding Sites Induce Weaker Long-Range Evolutionary Rate Gradients than Catalytic Sites in Enzymes

Enzymes exhibit a strong long-range evolutionary constraint that extends from their catalytic site and affects even distant sites, where site-specific evolutionary rate increases monotonically with distance. While protein–protein sites in enzymes were previously shown to induce only a weak conservation gradient, a comprehensive relationship between different types of functional sites in proteins and the magnitude of evolutionary rate gradients they induce has yet to be established. Here, we systematically calculate the evolutionary rate (dN/dS) of sites as a function of distance from different types of binding sites in enzymes and other proteins: catalytic sites, non-catalytic ligand binding sites, allosteric binding sites, and protein–protein interaction sites. We show that catalytic sites indeed induce significantly stronger evolutionary rate gradient than all other types of non-catalytic binding sites. In addition, catalytic sites in enzymes with no known allosteric function still induce strong long-range conservation gradients. Notably, the weak long-range conservation gradients induced by non-catalytic binding sites in enzymes is nearly identical in magnitude to those induced by ligand binding sites in non-enzymes. Finally, we show that structural determinants such as local solvent exposure of sites cannot explain the observed difference between catalytic and non-catalytic functional sites. Our results suggest that enzymes and non-enzymes share similar evolutionary constraints only when examined from the perspective of non-catalytic functional sites. Hence, the unique evolutionary rate gradient from catalytic sites in enzymes is likely driven by the optimization of catalysis rather than ligand binding and allosteric functions.     

Protein function annotation by local binding site surface similarity

Hundreds of protein crystal structures exist for proteins whose function cannot be confidently determined from sequence similarity. Surflex‐PSIM, a previously reported surface‐based protein similarity algorithm, provides an alternative method for hypothesizing function for such proteins. The method now supports fully automatic binding site detection and is fast enough to screen comprehensive databases of protein binding sites. The binding site detection methodology was validated on apo/holo cognate protein pairs, correctly identifying 91% of ligand binding sites in holo structures and 88% in apo structures where corresponding sites existed. For correctly detected apo binding sites, the cognate holo site was the most similar binding site 87% of the time. PSIM was used to screen a set of proteins that had poorly characterized functions at the time of crystallization, but were later biochemically annotated. Using a fully automated protocol, this set of 8 proteins was screened against ～60,000 ligand binding sites from the PDB. PSIM correctly identified functional matches that predated query protein biochemical annotation for five out of the eight query proteins. A panel of 12 currently unannotated proteins was also screened, resulting in a large number of statistically significant binding site matches, some of which suggest likely functions for the podorly characterized proteins. Proteins 2014; 82:679–694. © 2013 Wiley Periodicals, Inc.     

Solution Structure of Ectodomains of the Insulin Receptor Family: The Ectodomain of the Type 1 Insulin-Like Growth Factor Receptor Displays Asymmetry of Ligand Binding Accompanied by Limited Conformational Change

The insulin receptor (IR) and the homologous Type 1 insulin-like growth factor receptor (IGF-1R) are cell-surface tyrosine kinase receptors that effect signaling within the respective pathways of glucose metabolism and normal human growth. While ligand binding to these receptors is assumed to result in a structural transition within the receptor ectodomain that then effects signal transduction across the cell membrane, little is known about the molecular detail of these events. Presented here are small-angle X-ray scattering data obtained from the IR and IGF-1R ectodomains in solution. We show that, in solution, the ectodomains of IR and IGF-1R have a domain disposition that is very similar to that seen in the crystal structure of the ectodomain of IR, despite the constituent domains being in relatively sparse contact and potentially mobile. We also show that the IGF-1R ectodomain is capable of binding up to three molecules of IGF-1 in solution, with surprisingly little apparent change in relative domain disposition compared to the apo form. While the observed 3:1 ligand-binding stoichiometry appears to contradict earlier explanations of the absence of a bell-shaped dose-response curve for IGF-1R in ligand displacement assays, it is readily understood in the context of the harmonic oscillator model of the negative cooperativity of ligand binding to IGF-1R. Taken together, our findings suggest that the structural movements within these receptors upon ligand binding are small and are possibly limited to local rotation of domains.     

Coarse‐grained and all‐atom modeling of structural states and transitions in hemoglobin

Hemoglobin (Hb), an oxygen‐binding protein composed of four subunits (α1, α2, β1, and β2), is a well‐known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse‐grained (CG) modeling, all‐atom simulation, and structural data from X‐ray crystallography and wide‐angle X‐ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot‐spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all‐atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all‐atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Promotion of Enzyme Flexibility by Dephosphorylation and Coupling to the Catalytic Mechanism of a Phosphohexomutase

The enzyme phosphomannomutase/phosphoglucomutase (PMM/PGM) from Pseudomonas aeruginosa catalyzes an intramolecular phosphoryl transfer across its phosphosugar substrates, which are precursors in the synthesis of exoproducts involved in bacterial virulence. Previous structural studies of PMM/PGM have established a key role for conformational change in its multistep reaction, which requires a dramatic 180° reorientation of the intermediate within the active site. Here hydrogen-deuterium exchange by mass spectrometry and small angle x-ray scattering were used to probe the conformational flexibility of different forms of PMM/PGM in solution, including its active, phosphorylated state and the unphosphorylated state that occurs transiently during the catalytic cycle. In addition, the effects of ligand binding were assessed through use of a substrate analog. We found that both phosphorylation and binding of ligand produce significant effects on deuterium incorporation. Phosphorylation of the conserved catalytic serine has broad effects on residues in multiple domains and is supported by small angle x-ray scattering data showing that the unphosphorylated enzyme is less compact in solution. The effects of ligand binding are generally manifested near the active site cleft and at a domain interface that is a site of conformational change. These results suggest that dephosphorylation of the enzyme may play two critical functional roles: a direct role in the chemical step of phosphoryl transfer and secondly through propagation of structural flexibility. We propose a model whereby increased enzyme flexibility facilitates the reorientation of the reaction intermediate, coupling changes in structural dynamics with the unique catalytic mechanism of this enzyme.