Low-resolution structures of proteins in solution retrieved from X-ray scattering with a genetic algorithm.

Small-angle x-ray solution scattering (SAXS) is analyzed with a new method to retrieve convergent model structures that fit the scattering profiles. An arbitrary hexagonal packing of several hundred beads containing the problem object is defined. Instead of attempting to compute the Debye formula for all of the possible mass distributions, a genetic algorithm is employed that efficiently searches the configurational space and evolves best-fit bead models. Models from different runs of the algorithm have similar or identical structures. The modeling resolution is increased by reducing the bead radius together with the search space in successive cycles of refinement. The method has been tested with protein SAXS (0.001 < S < 0.06 A(-1)) calculated from x-ray crystal structures, adding noise to the profiles. The models obtained closely approach the volumes and radii of gyration of the known structures, and faithfully reproduce the dimensions and shape of each of them. This includes finding the active site cavity of lysozyme, the bilobed structure of gamma-crystallin, two domains connected by a stalk in betab2-crystallin, and the horseshoe shape of pancreatic ribonuclease inhibitor. The low-resolution solution structure of lysozyme has been directly modeled from its experimental SAXS profile (0.003 < S < 0.03 A(-1)). The model describes lysozyme size and shape to the resolution of the measurement. The method may be applied to other proteins, to the analysis of domain movements, to the comparison of solution and crystal structures, as well as to large macromolecular assemblies.     

Structural Insights into the Globular Tails of the Human Type V Myosins Myo5a,Myo5b,and Myo5c

Vertebrate type V myosins (MyoV) Myo5a, Myo5b, and Myo5c mediate transport of several different cargoes. All MyoV paralogs bind to cargo complexes mainly by their C-terminal globular domains. In absence of cargo, the globular domain of Myo5a inhibits its motor domain. Here, we report low-resolution SAXS models for the globular domains from human Myo5a, Myo5b, and Myo5c, which suggest very similar overall shapes of all three paralogs. We determined the crystal structures of globular domains from Myo5a and Myo5b, and provide a homology model for human Myo5c. When we docked the Myo5a crystal structure into a previously published electron microscopy density of the autoinhibited full-length Myo5a, only one domain orientation resulted in a good fit. This structural arrangement suggests the participation of additional region of the globular domain in autoinhibition. Quantification of the interaction of the Myo5a globular domain with its motor complex revealed a tight binding with dissociation half-life in the order of minutes, suggesting a rather slow transition between the active and inactive states.     

Lumbricus terrestris hemoglobin: A comparison of small-angle X-ray scattering and cryoelectron microscopy data

The quaternary structure of Lumbricus terrestris hemoglobin was investigated by small-angle x-ray scattering (SAXS). Based on the SAXS data from several independent experiments, a three-dimensional (3D) consensus model was established to simulate the solution structure of this complex protein at low resolution (about 3 nm) and to yield the particle dimensions. The model is built up from a large number of small spheres of different weights, a result of the two-step procedure used to calculate the SAXS model. It accounts for the arrangement of 12 subunits in a hexagonal bilayer structure and for an additional central unit of cylinder-like shape. This model provides an excellent fit of the experimental scattering curve of the protein up to h = 1 nm⁻¹ and a nearly perfect fit of the experimental distance distribution function p(r) in the whole range. Scattering curves and p(r) functions were also calculated for low-resolution models based on 3D reconstructions obtained by cryoelectron microscopy (EM). The calculated functions of these models also provide a very good fit of the experimental scattering curve (even at h > 1 nm⁻¹) and p(r) function, if hydration is taken into account and the original model coordinates are slightly rescaled. The comparison of models reveals that both the SAXS-based and the EM-based model lead to a similar simulation of the protein structure and to similar particle dimensions. The essential differences between the models concern the hexagonal bilayer arrangement (eclipsed in the SAXS model, one layer slightly rotated in the EM model), and the mass distribution, mainly on the surface and in the central part of the protein complex. © John Wiley & Sons, Inc. Biopoly 45: 289–298, 1998     

Conformation of the Drosophila motor protein non-claret disjunctional in solution from X-ray and neutron scattering

The quaternary structures of monomeric and dimeric Drosophila non-claret disjunctional (ncd) constructs were investigated using synchrotron x-ray and neutron solution scattering, and their low resolution shapes were restored ab initio from the scattering data. The experimental curves were further compared with those computed from crystallographic models of one monomeric and three available dimeric ncd structures in the microtubule-independent ADP-bound state. These comparisons indicate that accounting for the missing parts in the crystal structures for all these constructs is indispensable to obtain reasonable fits to the scattering patterns. A ncd construct (MC6) lacking the coiled-coil region is monomeric in solution, but the calculated scattering from the crystallographic monomer yields a poor fit to the data. A tentative configuration of the missing C-terminal residues in the form of an antiparallel beta-sheet was found that significantly improves the fit. The atomic model of a short dimeric ncd construct (MC5) without 2-fold symmetry is found to fit the data better than the symmetric models. Addition of the C-terminal residues to both head domains gives an excellent fit to the x-ray and neutron experimental data, although the orientation of the beta-sheet differs from that of the monomer. The solution structure of the long ncd construct (MC1) including complete N-terminal coiled-coil and motor domains is modeled by adding a straight coiled-coil section to the model of MC5.     

Solution structure of C-terminal Escherichia coli translation initiation factor IF2 by small-angle X-ray scattering

Initiation of protein synthesis in bacteria involves the combined action of three translation initiation factors, including translation initiation factor IF2. Structural knowledge of this bacterial protein is scarce. A fragment consisting of the four C-terminal domains of IF2 from Escherichia coli was expressed, purified, and characterized by small-angle X-ray scattering (SAXS), and from the SAXS data, a radius of gyration of 43 +/- 1 A and a maximum dimension of approximately 145 A were obtained for the molecule. Furthermore, the SAXS data revealed that E. coli IF2 in solution adopts a structure that is significantly different from the crystal structure of orthologous aIF5B from Methanobacterium thermoautotrophicum. This crystal structure constitutes the only atomic resolution structural knowledge of the full-length factor. Computer programs were applied to the SAXS data to provide an initial structural model for IF2 in solution. The low-resolution nature of SAXS prevents the elucidation of a complete and detailed structure, but the resulting model for C-terminal E. coli IF2 indicates important structural differences between the aIF5B crystal structure and IF2 in solution. The chalice-like structure with a highly exposed alpha-helical stretch observed for the aIF5B crystal structure was not found in the structural model of IF2 in solution, in which domain VI-2 is moved closer to the rest of the protein.     

A model for the three-dimensional structure of human plasma vitronectin from small-angle scattering measurements

Small-angle X-ray scattering (SAXS) measurements were used to characterize vitronectin, a circulatory protein found in human plasma that functions in regulating cell adhesion and migration, as well as proteolytic cascades that affect blood coagulation, fibrinolysis, and pericellular proteolysis. SAXS measurements were taken over a 3-fold range of protein concentrations, yielding data that characterize a monodisperse system of particles with an average radius of gyration of 30.3 +/- 0.6 A and a maximum linear dimension of 110 A. Shape restoration was applied to the data to produce two models of the solution structure of the ligand-free protein. A low-resolution model of the protein was generated that indicates the protein to be roughly peanut-shaped. A better understanding of the domain structure of vitronectin resulted from low-resolution models developed from available high-resolution structures of the domains. These domains include the N-terminal domain that was determined experimentally by NMR [Mayasundari, A., Whittemore, N. A., Serpersu, E. H., and Peterson, C. B. (2004) J. Biol. Chem. 279, 29359-29366] and the docked structure of the central and C-terminal domains that were determined by computational threading [Xu, D., Baburaj, K., Peterson, C. B., and Xu, Y. (2001) Proteins: Struct., Funct., Genet. 44, 312-320]. This model provides an indication of the disposition of the central domain and C-terminal heparin-binding domains of vitronectin with respect to the N-terminal somatomedin B (SMB) domain. This model constructed from the available domain structures, which agrees with the low-resolution model produced from the SAXS data, shows the SMB domain well separated from the central and heparin-binding domains by a disordered linker (residues 54-130). Also, binding sites within the SMB domain are predicted to be well exposed to the surrounding solvent for ease of access to its various ligands.     

Integration of small-angle X-ray scattering data into structural modeling of proteins and their assemblies

A major challenge in structural biology is to determine the configuration of domains and proteins in multidomain proteins and assemblies, respectively. All available data should be considered to maximize the accuracy and precision of these models. Small-angle X-ray scattering (SAXS) efficiently provides low-resolution experimental data about the shapes of proteins and their assemblies. Thus, we integrated SAXS profiles into our software for modeling proteins and their assemblies by satisfaction of spatial restraints. Specifically, we modeled the quaternary structures of multidomain proteins with structurally defined rigid domains as well as quaternary structures of binary complexes of structurally defined rigid proteins. In addition to SAXS profiles and the component structures, we used stereochemical restraints and an atomic distance-dependent statistical potential. The scoring function is optimized by a biased Monte Carlo protocol, including quasi-Newton and simulated annealing schemes. The final prediction corresponds to the best scoring solution in the largest cluster of many independently calculated solutions. To quantify how well the quaternary structures are determined based on their SAXS profiles, we used a benchmark of 12 simulated examples as well as an experimental SAXS profile of the homotetramer d-xylose isomerase. Optimization of the SAXS-dependent scoring function generally results in accurate models if sufficiently precise approximations for the constituent rigid bodies are available; otherwise, the best scoring models can have significant errors. Thus, SAXS profiles can play a useful role in the structural characterization of proteins and assemblies if they are combined with additional data and used judiciously. Our integration of a SAXS profile into modeling by satisfaction of spatial restraints will facilitate further integration of different kinds of data for structure determination of proteins and their assemblies.     

Accurate flexible fitting of high-resolution protein structures to small-angle x-ray scattering data using a coarse-grained model with implicit hydration shell

Small-angle x-ray scattering (SAXS) is a powerful technique widely used to explore conformational states and transitions of biomolecular assemblies in solution. For accurate model reconstruction from SAXS data, one promising approach is to flexibly fit a known high-resolution protein structure to low-resolution SAXS data by computer simulations. This is a highly challenging task due to low information content in SAXS data. To meet this challenge, we have developed what we believe to be a novel method based on a coarse-grained (one-bead-per-residue) protein representation and a modified form of the elastic network model that allows large-scale conformational changes while maintaining pseudobonds and secondary structures. Our method optimizes a pseudoenergy that combines the modified elastic-network model energy with a SAXS-fitting score and a collision energy that penalizes steric collisions. Our method uses what we consider a new implicit hydration shell model that accounts for the contribution of hydration shell to SAXS data accurately without explicitly adding waters to the system. We have rigorously validated our method using five test cases with simulated SAXS data and three test cases with experimental SAXS data. Our method has successfully generated high-quality structural models with root mean-squared deviation of 1 ∼ 3 Å from the target structures.     

Addition of missing loops and domains to protein models by x-ray solution scattering

Inherent flexibility and conformational heterogeneity in proteins can often result in the absence of loops and even entire domains in structures determined by x-ray crystallographic or NMR methods. X-ray solution scattering offers the possibility of obtaining complementary information regarding the structures of these disordered protein regions. Methods are presented for adding missing loops or domains by fixing a known structure and building the unknown regions to fit the experimental scattering data obtained from the entire particle. Simulated annealing was used to minimize a scoring function containing the discrepancy between the experimental and calculated patterns and the relevant penalty terms. In low-resolution models where interface location between known and unknown parts is not available, a gas of dummy residues represents the missing domain. In high-resolution models where the interface is known, loops or domains are represented as interconnected chains (or ensembles of residues with spring forces between the C(alpha) atoms), attached to known position(s) in the available structure. Native-like folds of missing fragments can be obtained by imposing residue-specific constraints. After validation in simulated examples, the methods have been applied to add missing loops or domains to several proteins where partial structures were available.     

The SAXS solution structure of RF1 differs from its crystal structure and is similar to its ribosome bound cryo-EM structure

Bacterial class I release factors (RFs) are seen by cryo-electron microscopy (cryo-EM) to span the distance between the ribosomal decoding and peptidyl transferase centers during translation termination. The compact conformation of bacterial RF1 and RF2 observed in crystal structures will not span this distance, and large structural rearrangements of RFs have been suggested to play an important role in termination. We have collected small-angle X-ray scattering (SAXS) data from E. coli RF1 and from a functionally active truncated RF1 derivative. Theoretical scattering curves, calculated from crystal and cryo-EM structures, were compared with the experimental data, and extensive analyses of alternative conformations were made. Low-resolution models were constructed ab initio, and by rigid-body refinement using RF1 domains. The SAXS data were compatible with the open cryo-EM conformation of ribosome bound RFs and incompatible with the crystal conformation. These conclusions obviate the need for assuming large conformational changes in RFs during termination.     

Low-resolution structure and fluorescence anisotropy analysis of protein tyrosine phosphatase eta catalytic domain

The rat protein tyrosine phosphatase eta, rPTPeta, is a class I "classical" transmembrane RPTP, with an intracellular portion composed of a unique catalytic region. The rPTPeta and the human homolog DEP-1 are downregulated in rat and human neoplastic cells, respectively. However, the malignant phenotype is reverted after exogenous reconstitution of rPTPeta, suggesting that its function restoration could be an important tool for gene therapy of human cancers. Using small-angle x-ray scattering (SAXS) and biophysical techniques, we characterized the intracellular catalytic domain of rat protein tyrosine phosphatase eta (rPTPetaCD) in solution. The protein forms dimers in solution as confirmed by SAXS data analysis. The SAXS data also indicated that rPTPetaCD dimers are elongated and have an average radius of gyration of 2.65 nm and a D(max) of 8.5 nm. To further study the rPTPetaCD conformation in solution, we built rPTPetaCD homology models using as scaffolds the crystallographic structures of RPTPalpha-D1 and RPTPmicro-D1 dimers. These models were, then, superimposed onto ab initio low-resolution SAXS structures. The structural comparisons and sequence alignment analysis of the putative dimerization interfaces provide support to the notion that the rPTPetaCD dimer architecture is more closely related to the crystal structure of autoinhibitory RPTPalpha-D1 dimer than to the dimeric arrangement exemplified by RPTPmicro-D1. Finally, the characterization of rPTPetaCD by fluorescence anisotropy measurements demonstrates that the dimer dissociation is concentration dependent with a dissociation constant of 21.6 +/- 2.0 microM.     

Structural features of the full-length adaptor protein GADS in solution determined using small-angle X-ray scattering

The Grb2-related adaptor protein GADS plays a central role during the initial phases of signal transduction in T lymphocytes. GADS possesses N- and C-terminal Src homology 3 (SH3) domains flanking a central Src homology 2 (SH2) domain and a 126-residue region rich in glutamine and proline residues, presumed to be largely unstructured. The SH2 domain of GADS binds the adaptor protein LAT; the C-terminal SH3 domain pairs GADS to the adaptor protein SLP-76, whereas the function of the central region is unknown. High-resolution three-dimensional models are available for the isolated SH2 and C-terminal SH3 domains in complex with their respective binding partners, LAT and SLP-76. However, in part because of its intrinsic instability, there is no structural information for the entire GADS molecule. Here, we report the low-resolution structure of full-length GADS in solution using small-angle x-ray scattering (SAXS). Based on the SAXS data, complemented by gel filtration experiments, we show that full-length GADS is monomeric in solution and that its overall structural parameters are smaller than those expected for a protein with a long unstructured region. Ab initio and rigid body modeling of the SAXS data reveal that full-length GADS is a relatively compact molecule and that the potentially unstructured region retains a significant degree of structural order. The biological function of GADS is discussed based on its overall structure.     

Biophysical characterisation of human LincRNA-p21 sense and antisense Alu inverted repeats

Human Long Intergenic Noncoding RNA-p21 (LincRNA-p21) is a regulatory noncoding RNA that plays an important role in promoting apoptosis. LincRNA-p21 is also critical in down-regulating many p53 target genes through its interaction with a p53 repressive complex. The interaction between LincRNA-p21 and the repressive complex is likely dependent on the RNA tertiary structure. Previous studies have determined the two-dimensional secondary structures of the sense and antisense human LincRNA-p21 AluSx1 IRs using SHAPE. However, there were no insights into its three-dimensional structure. Therefore, we in vitro transcribed the sense and antisense regions of LincRNA-p21 AluSx1 Inverted Repeats (IRs) and performed analytical ultracentrifugation, size exclusion chromatography, light scattering, and small angle X-ray scattering (SAXS) studies. Based on these studies, we determined low-resolution, three-dimensional structures of sense and antisense LincRNA-p21. By adapting previously known two-dimensional information, we calculated their sense and antisense high-resolution models and determined that they agree with the low-resolution structures determined using SAXS. Thus, our integrated approach provides insights into the structure of LincRNA-p21 Alu IRs. Our study also offers a viable pipeline for combining the secondary structure information with biophysical and computational studies to obtain high-resolution atomistic models for long noncoding RNAs.     

Determination of protein oligomeric structure from small‐angle X‐ray scattering

Small‐angle X‐ray scattering (SAXS) is useful for determining the oligomeric states and quaternary structures of proteins in solution. The average molecular mass in solution can be calculated directly from a single SAXS curve collected on an arbitrary scale from a sample of unknown protein concentration without the need for beamline calibration or protein standards. The quaternary structure in solution can be deduced by comparing the experimental SAXS curve to theoretical curves calculated from proposed models of the oligomer. This approach is especially robust when the crystal structure of the target protein is known, and the candidate oligomer models are derived from the crystal lattice. When SAXS data are obtained at multiple protein concentrations, this analysis can provide insight into dynamic self‐association equilibria. Herein, we summarize the computational methods that are used to determine protein molecular mass and quaternary structure from SAXS data. These methods are organized into a workflow and demonstrated with four case studies using experimental SAXS data from the published literature.     

High-resolution crystal structure reveals a HEPN domain at the C-terminal region of S. cerevisiae RNA endonuclease Swt1

Swt1 is an RNA endonuclease that plays an important role in quality control of nuclear messenger ribonucleoprotein particles (mRNPs) in eukaryotes; however, its structural details remain to be elucidated. Here, we report the crystal structure of the C-terminal (CT) domain of Swt1 from Saccharomyces cerevisiae, which shares common characteristics of higher eukaryotes and prokaryotes nucleotide binding (HEPN) domain superfamily. To study in detail the full-length protein structure, we analyzed the low-resolution architecture of Swt1 in solution using small angle X-ray scattering (SAXS) method. Both the CT domain and middle domain exhibited a good fit upon superimposing onto the molecular envelope of Swt1. Our study provides the necessary structural information for detailed analysis of the functional role of Swt1, and its importance in the process of nuclear mRNP surveillance.     

Mechanism of the interaction of beta(2)-glycoprotein I with negatively charged phospholipid membranes.

In an attempt to understand the multifunctional involvement of beta(2)-glycoprotein I (beta(2)GPI) in autoimmune diseases, thrombosis, atherosclerosis, and inflammatory processes, substantial interest is focused on the interaction of beta(2)GPI with negatively charged ligands, in particular, with acidic phospholipids. In this study, unilamellar vesicles composed of cardiolipin were used as in vitro membrane system to test and further refine a model of interaction based on the crystal structure of beta(2)GPI. The data suggest that beta(2)GPI anchors to the membrane surface with its hydrophobic loop adjacent to the positively charged lysine rich region in domain V. Subsequently, beta(2)GPI penetrates the membrane interfacial headgroup region as indicated by a restriction of the lipid side chain mobility, but without formation of a nonbilayer lipid phase. A structural rearrangement of beta(2)GPI upon lipid binding was detected by microcalorimetry and may result in the exposure of cryptic epitopes located in the complement control protein domains. This lipid-dependent conformational change may induce oligomerization of beta(2)GPI and promote intermolecular associations. Thus, the aggregation tendency of beta(2)GPI may serve as the basis for the formation of a molecular link between cells but may also be an essential feature for binding of autoantibodies and hence determine the role of beta(2)GPI in autoimmune diseases.     

A covalent linker allows for membrane targeting of an oxylipin biosynthetic complex

A naturally occurring bifunctional protein from Plexaura homomalla links sequential catalytic activities in an oxylipin biosynthetic pathway. The C-terminal lipoxygenase (LOX) portion of the molecule catalyzes the transformation of arachidonic acid (AA) to the corresponding 8 R-hydroperoxide, and the N-terminal allene oxide synthase (AOS) domain promotes the conversion of the hydroperoxide intermediate to the product allene oxide (AO). Small-angle X-ray scattering data indicate that in the absence of a covalent linkage the two catalytic domains that transform AA to AO associate to form a complex that recapitulates the structure of the bifunctional protein. The SAXS data also support a model for LOX and AOS domain orientation in the fusion protein inferred from a low-resolution crystal structure. However, results of membrane binding experiments indicate that covalent linkage of the domains is required for Ca (2+)-dependent membrane targeting of the sequential activities, despite the noncovalent domain association. Furthermore, membrane targeting is accompanied by a conformational change as monitored by specific proteolysis of the linker that joins the AOS and LOX domains. Our data are consistent with a model in which Ca (2+)-dependent membrane binding relieves the noncovalent interactions between the AOS and LOX domains and suggests that the C2-like domain of LOX mediates both protein-protein and protein-membrane interactions.     

Ceruloplasmin: Macromolecular Assemblies with Iron-Containing Acute Phase Proteins

Copper-containing ferroxidase ceruloplasmin (Cp) forms binary and ternary complexes with cationic proteins lactoferrin (Lf) and myeloperoxidase (Mpo) during inflammation. We present an X-ray crystal structure of a 2Cp-Mpo complex at 4.7 Å resolution. This structure allows one to identify major protein–protein interaction areas and provides an explanation for a competitive inhibition of Mpo by Cp and for the activation of p-phenylenediamine oxidation by Mpo. Small angle X-ray scattering was employed to construct low-resolution models of the Cp-Lf complex and, for the first time, of the ternary 2Cp-2Lf-Mpo complex in solution. The SAXS-based model of Cp-Lf supports the predicted 1∶1 stoichiometry of the complex and demonstrates that both lobes of Lf contact domains 1 and 6 of Cp. The 2Cp-2Lf-Mpo SAXS model reveals the absence of interaction between Mpo and Lf in the ternary complex, so Cp can serve as a mediator of protein interactions in complex architecture. Mpo protects antioxidant properties of Cp by isolating its sensitive loop from proteases. The latter is important for incorporation of Fe³⁺ into Lf, which activates ferroxidase activity of Cp and precludes oxidation of Cp substrates. Our models provide the structural basis for possible regulatory role of these complexes in preventing iron-induced oxidative damage.     

Structure, dynamics and domain organization of the repeat protein Cin1 from the apple scab fungus

Venturia inaequalis is a hemi-biotrophic fungus that causes scab disease of apple. A recently-identified gene from this fungus, cin1 (cellophane-induced 1), is up-regulated over 1000-fold in planta and considerably on cellophane membranes, and encodes a cysteine-rich secreted protein of 523 residues with eight imperfect tandem repeats of ~60 amino acids. The Cin1 sequence has no homology to known proteins and appears to be genus-specific; however, Cin1 repeats and other repeat domains may be structurally similar. An NMR-derived structure of the first two repeat domains of Cin1 (Cin1-D1D2) and a low-resolution model of the full-length protein (Cin1-FL) using SAXS data were determined. The structure of Cin1-D1D2 reveals that each domain comprises a core helix-loop-helix (HLH) motif as part of a three-helix bundle, and is stabilized by two intra-domain disulfide bonds. Cin1-D1D2 adopts a unique protein fold as DALI and PDBeFOLD analysis identified no structural homology. A (15)N backbone NMR dynamic analysis of Cin1-D1D2 showed that a short stretch of the inter-domain linker has large amplitude motions that give rise to reciprocal domain-domain mobility. This observation was supported by SAXS data modeling, where the scattering length density envelope remains thick at the domain-domain boundary, indicative of inter-domain dynamics. Cin1-FL SAXS data models a loosely-packed arrangement of domains, rather than the canonical parallel packing of adjacent HLH repeats observed in α-solenoid repeat proteins. Together, these data suggest that the repeat domains of Cin1 display a "beads-on-a-string" organization with inherent inter-domain flexibility that is likely to facilitate interactions with target ligands.