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.     

Small-angle X-ray scattering constraints and local geometry like secondary structures can construct a coarse-grained protein model at amino acid residue resolution

We have developed a new methodology that determines protein structures using small-angle X-ray scattering (SAXS) data. The current bottlenecks in determining the protein structures require a new strategy using the simple design of an experiment, and SAXS is suitable for this purpose in spite of its low information content. First we demonstrated that SAXS constraints work additively to NMR-derived information in calculating structures. Next, structure calculations for nine proteins taking different folds were performed using the SAXS constraints combined with the NMR-derived distance restraints for local geometry such as secondary structures or those for tertiary structure. The results show that the SAXS constraints complemented the tertiary-structural information for all the proteins, and that accuracy of the structures thus obtained with SAXS constraints and local geometrical restraints ranged from 1.85 to 4.33 Å. Based on these results, we were able to construct a coarse-grained protein model at amino acid residue resolution.     

Solution conformation of adenovirus virus associated RNA-I and its interaction with PKR

Adenovirus virus-associated RNA (VA_I) provides protection against the host antiviral response in part by inhibiting the interferon-induced double stranded RNA-activated protein kinase (PKR). VA_I consists of three base-paired regions; the apical stem responsible for the interaction with double-stranded RNA binding motifs (dsRBMs) of PKR, the central stem required for inhibition, and the terminal stem. The solution conformation of VA_I and VA_I lacking the terminal stem were determined using SAXS that suggested extended conformations that are in agreement with their secondary structures. Solution conformations of VA_I lacking the terminal stem in complex with the dsRBMs of PKR indicated that the apical stem interacts with both dsRNA-binding motifs whereas the central stem does not. Hydrodynamic properties calculated from ab initio models were compared to experimentally determined parameters for model validation. Furthermore, SAXS envelopes were used as a constraint for the in silico modeling of tertiary structure for RNA and RNA–protein complex. Finally, full-length PKR was also studied, but concentration-dependent changes in hydrodynamic parameters prevented ab initio shape determination. Taken together, results provide an improved structural framework that further our understanding of the role VA_I plays in evading host innate immune responses.     

Insulin-like growth factor binding protein-2: NMR analysis and structural characterization of the N-terminal domain

The insulin-like growth factor binding proteins are a family of six proteins (IGFBP-1 to -6) that bind insulin-like growth factors-I and -II (IGF-I/II) with high affinity. In addition to regulating IGF actions, IGFBPs have IGF-independent functions. IGFBP-2, the largest member of this family, is over-expressed in many cancers and has been proposed as a possible target for the development of novel anti-cancer therapeutics. The IGFBPs have a common architecture consisting of conserved N- and C-terminal domains joined by a variable linker domain. The solution structure and dynamics of the C-terminal domain of human IGFBP-2 have been reported (Kuang Z. et al. J. Mol. Biol. 364, 690-704, 2006) but neither the N-domain (N-BP-2) nor the linker domain have been characterised. Here we present NMR resonance assignments for human N-BP-2, achieved by recording spectra at low protein concentration using non-uniform sampling and maximum entropy reconstruction. Analysis of secondary chemical shifts shows that N-BP-2 possesses a secondary structure similar to that of other IGFBPs. Although aggregation hampered determination of the solution structure for N-BP-2, a homology model was generated based on the high degree of sequence and structure homology exhibited by the IGFBPs. This model was consistent with experimental NMR and SAXS data and displayed some unique features such as a Pro/Ala-rich non-polar insert, which formed a flexible solvent-exposed loop on the surface of the protein opposite to the IGF-binding interface. NMR data indicated that this loop could adopt either of two alternate conformations in solution - an entirely flexible conformation and one containing nascent helical structure. This loop and an adjacent poly-proline sequence may comprise a potential SH3 domain interaction site for binding to other proteins.     

Structural insights into the beta-mannosidase from T. reesei obtained by synchrotron small-angle X-ray solution scattering enhanced by X-ray crystallography

A molecular envelope of the beta-mannosidase from Trichoderma reesei has been obtained by combined use of solution small-angle X-ray scattering (SAXS) and protein crystallography. Crystallographic data at 4 A resolution have been used to enhance informational content of the SAXS data and to obtain an independent, more detailed protein shape. The phased molecular replacement technique using a low resolution SAXS model, building, and refinement of a free atom model has been employed successfully. The SAXS and crystallographic free atom models exhibit a similar globular form and were used to assess available crystallographic models of glycosyl hydrolases. The structure of the beta-galactosidase, a member of a family 2, clan GHA glycosyl hydrolases, shows an excellent fit to the experimental molecular envelope and distance distribution function of the beta-mannosidase, indicating gross similarities in their three-dimensional structures. The secondary structure of beta-mannosidase quantified by circular dichroism measurements is in a good agreement with that of beta-galactosidase. We show that a comparison of distance distribution functions in combination with 1D and 2D sequence alignment techniques was able to restrict the number of possible structurally homologous proteins. The method could be applied as a general method in structural genomics and related fields once protein solution scattering data are available.     

A solution SAXS study of Borrelia burgdorferi OspA, a protein containing a single-layer beta-sheet.

The crystal structure of a soluble form of Borrelia burgdorferi outer surface protein A (OspA) complexed with the Fab fragment of a monoclonal antibody has revealed an unusual structure that has a repetitive antiparallel beta topology with a nonglobular, single layer beta-sheet connecting the globular N- and C-terminal domains. Earlier NMR studies have shown that the local structure of OspA including the single layer beta-sheet is similar to the crystal structure. Here we report a small angle X-ray scattering (SAXS) study of the global conformation of OspA in solution. The radius of gyration (Rg) and the length distribution function (P(r)) of OspA measured by SAXS in solution are nearly identical to the calculated ones from the crystal structure, respectively. The NMR and SAXS experiments complement each other to show that OspA including the central single-layer beta-sheet is a stable structure in solution, and that the OspA crystal structure represents the predominant solution conformation of the protein.     

Evolutionary model of protein secondary structure capable of revealing new biological relationships

Ancestral sequence reconstruction has had recent success in decoding the origins and the determinants of complex protein functions. However, phylogenetic analyses of remote homologues must handle extreme amino acid sequence diversity resulting from extended periods of evolutionary change. We exploited the wealth of protein structures to develop an evolutionary model based on protein secondary structure. The approach follows the differences between discrete secondary structure states observed in modern proteins and those hypothesized in their immediate ancestors. We implemented maximum likelihood-based phylogenetic inference to reconstruct ancestral secondary structure. The predictive accuracy from the use of the evolutionary model surpasses that of comparative modeling and sequence-based prediction; the reconstruction extracts information not available from modern structures or the ancestral sequences alone. Based on a phylogenetic analysis of a sequence-diverse protein family, we showed that the model can highlight relationships that are evolutionarily rooted in structure and not evident in amino acid-based analysis.     

Determination of a molecular shape for netrin-4 from hydrodynamic and small angle X-ray scattering measurements

As part of a continuing investigation of netrins, an emerging class of extracellular matrix proteins that are involved in axon guidance activity, we have used dynamic light scattering (DLS) and small angle X-ray scattering to investigate the solution conformation of a truncated version of netrin-4 (Δnetrin-4) that lacks the C-terminal portion. The protein is characterized by a hydrodynamic (Stokes) radius (r(H)) of 4.60 (±0.20) nm, a radius of gyration (r(G)) of 4.42 (±0.20) nm and a maximum particle dimension (D(max)) of 16nm. More detailed ab initio modeling of the SAXS data indicates an extended rod like conformation for Δnetrin-4 in solution-a concept supported by the excellent agreement observed between experimental parameter estimates and those calculated for the ab initio models for Δnetrin-4 by the HYDROPRO program.     

Structural studies of the manganese stabilizing subunit in photosystem II

Photosystem II (PSII) is the plant photosynthetic reaction center that carries out the light driven oxidation of water. The water splitting reactions are catalyzed at a tetranuclear manganese cluster. The manganese stabilizing protein (MSP) of PSII stabilizes the manganese cluster and accelerates the rate of oxygen evolution. MSP can be removed from PSII, with an accompanying decrease in activity. Either an Escherichia coli expressed version of MSP or native, plant MSP can be rebound to the PSII reaction center; MSP reconstitution reverses the deleterious effects associated with MSP removal. We have employed Fourier transform infrared (FTIR) spectroscopy and solution small angle x-ray scattering (SAXS) techniques to investigate the structure of MSP in solution and to define the structural changes that occur before and after reconstitution to PSII. FTIR and SAXS are complementary, because FTIR spectroscopy detects changes in MSP secondary structure and SAXS detects changes in MSP size/shape. From the SAXS data, we conclude that the size/shape and domain structure of MSP do not change when MSP binds to PSII. From FTIR data acquired before and after reconstitution, we conclude that the reconstitution-induced increase in beta-sheet content, which was previously reported, persists after MSP is removed from the PSII reaction center. However, the secondary structural change in MSP is metastable after removal from PSII, which indicates that this form of MSP is not the lowest energy conformation in solution.     

Production and characterization of recombinant P1 adhesin essential for adhesion,gliding, and antigenic variation in the human pathogenic bacterium,Mycoplasma pneumoniae

Mycoplasma pneumoniae forms an attachment organelle at one cell pole, binds to the host cell surface, and glides via a unique mechanism. A 170-kDa protein, P1 adhesin, present on the organelle surface plays a critical role in the binding and gliding process. In this study, we obtained a recombinant P1 adhesin comprising 1476 amino acid residues, excluding the C-terminal domain of 109 amino acids that carried the transmembrane segment, that were fused to additional 17 amino acid residues carrying a hexa-histidine (6?×?His) tag using an Escherichia coli expression system. The recombinant protein showed solubility, and chirality in circular dichroism (CD). The results of analytical gel filtration, ultracentrifugation, negative-staining electron microscopy, and small-angle X-ray scattering (SAXS) showed that the recombinant protein exists in a monomeric form with a uniformly folded structure. SAXS analysis suggested the presence of a compact and ellipsoidal structure rather than random or molten globule-like conformation. Structure model based on SAXS results fitted well with the corresponding structure obtained with cryo-electron tomography from a closely related species, M. genitalium. This recombinant protein may be useful for structural and functional studies as well as for the preparation of antibodies for medical applications.     

Solution structure of subunit a, a 104-363, of the Saccharomyces cerevisiae V-ATPase and the importance of its C-terminus in structure formation

The 95 kDa subunit a of eukaryotic V-ATPases consists of a C-terminal, ion-translocating part and an N-terminal cytosolic domain. The latter's N-terminal domain (~40 kDa) is described to bind in an acidification-dependent manner with cytohesin-2 (ARNO), giving the V-ATPase the putative function as pH-sensing receptor. Recently, the solution structure of the very N-terminal segment of the cytosolic N-terminal domain has been solved. Here we produced the N-terminal truncated form SCa??????? of the N-terminal domain (SCa?????) of the Saccharomyces cerevisiae V-ATPase and determined its low resolution solution structure, derived from SAXS data. SCa??????? shows an extended S-like conformation with a width of about 3.88 nm and a length of 11.4 nm. The structure has been superimposed into the 3D reconstruction of the related A?A? ATP synthase from Pyrococcus furiosus, revealing that the SCa??????? fits well into the density of the collar structure of the enzyme complex. To understand the importance of the C-terminus of the protein SCa?????, and to determine the localization of the N- and C-termini in SCa???????, the C-terminal truncated form SCa??????? was produced and analyzed by SAXS. Comparison of the SCa??????? and SCa??????? shapes showed that the additional loop region in SCa??????? consists of the C-terminal residues. Whereas SCa??????? is monomeric in solution, SCa??????? forms a dimer, indicating the importance of the very C-terminus in structure formation. Finally, the solution structure of SCa??????? and SCa??????? will be discussed in terms of the topological arrangement of subunit a and cytoheisn-2 in V-ATPases.     

Limited tolerance towards folded elements during secretion of the autotransporter Hbp

Many virulence factors secreted by pathogenic Gram-negative bacteria belong to the autotransporter (AT) family. ATs consist of a passenger domain, which is the actual secreted moiety, and a beta-domain that facilitates the transfer of the passenger domain across the outer membrane. Here, we analysed folding and translocation of the AT passenger, using Escherichia coli haemoglobin protease (Hbp) as a model protein. Dual cysteine mutagenesis, instigated by the unique crystal structure of the Hbp passenger, resulted in intramolecular disulphide bond formation dependent on the periplasmic enzyme DsbA. A small loop tied off by a disulphide bond did not interfere with secretion of Hbp. In contrast, a bond between different domains of the Hbp passenger completely blocked secretion resulting in degradation by the periplasmic protease DegP. In the absence of DegP, a translocation intermediate accumulated in the outer membrane. A similar jammed intermediate was formed upon insertion of a calmodulin folding moiety into Hbp. The data suggest that Hbp can fold in the periplasm but must retain a certain degree of flexibility and/or modest width to allow translocation across the outer membrane.     

Structure and flexibility of the complete periplasmic domain of BamA: the protein insertion machine of the outer membrane

Folding and insertion of β-barrel outer membrane proteins (OMPs) is essential for Gram-negative bacteria. This process is mediated by the multiprotein complex BAM, composed of the essential β-barrel OMP BamA and four lipoproteins (BamBCDE). The periplasmic domain of BamA is key for its function and contains five "polypeptide transport-associated" (POTRA) repeats. Here, we report the crystal structure of the POTRA4-5 tandem, containing the essential for BAM complex formation and cell viability POTRA5. The domain orientation observed in the crystal is validated by solution NMR and SAXS. Using previously determined structures of BamA POTRA1-4, we present a spliced model of the entire BamA periplasmic domain validated by SAXS. Solution scattering shows that conformational flexibility between POTRA2 and 3 gives rise to compact and extended conformations. The length of BamA in its extended conformation suggests that the protein may bridge the inner and outer membranes across the periplasmic space.     

Structural stability and solution structure of chaperonin GroES heptamer studied by synchrotron small-angle X-ray scattering

The GroES protein from Escherichia coli is a well-known member of the molecular chaperones. GroES consists of seven identical 10 kDa subunits, and forms a dome-like oligomeric structure. In order to obtain information on the structural stability and unfolding-refolding mechanism of GroES protein, especially at protein concentrations (0.4-1.2 mM GroES monomer) that would mimic heat stress conditions in vivo, we have performed synchrotron small-angle X-ray scattering (SAXS) experiments. Surprisingly, in spite of the high protein concentration, reversibility in the unfolding-refolding reaction was confirmed by SAXS experiments structurally. Although the unfolding-refolding reaction showed an apparent single transition with a Cm of 1.1 M guanidium hydrochloride, a more detailed analysis of this transition demonstrated that the unfolding mechanism could be best explained by a sequential three-state model, which consists of native heptamer, dissociated monomer, and unfolded monomer. Together with our previous result that GroES unfolded completely via a partially folded monomer according to a three-state model at low protein concentration (5 microM monomer), the unfolding-refolding mechanism of GroES protein could be explained uniformly by the three-state model from low to high protein concentrations. Furthermore, to clarify an ambiguity of the native GroES structure in solution, especially mobile loop structures, we have estimated a solution structure of GroES using SAXS profiles obtained from experiments and simulation analysis. The result suggested that the native structure of GroES in solution was very similar to that seen in GroES-GroEL complex determined by crystallography.     

Alterations in the protein-synthesis, -degradation and/or -secretion rates in hepatic subcellular fractions of selenium-deficient mice.

Conformational studies of myelin basic protein (MBP) in solution generally have used protein purified in organic solvents and acid. The use of such conditions raises the possibility that the secondary structure reported for the basic protein represents a denatured state. Therefore we have purified this protein by using a procedure that avoids denaturants. Bovine myelin was extracted with 0.2 M-CaCl2 and the protein was purified from the supernatant by chromatography on Sephadex G-75. The conformation of the basic protein was characterized by using c.d. and 1H-n.m.r. spectroscopy. In solution, it appeared to be predominantly randomly coiled, with only small segments of persistent structure. However, in the presence of myristoyl lysophosphatidylcholine the secondary structure of MBP became more ordered, and sedimentation-velocity experiments showed that MBP aggregated. Comparison of our results with published data indicates that Ca2+-extracted basic protein behaves similarly to the protein purified by traditional methods with respect to its ordered conformation in solution in the absence and in the presence of lipid and with respect to its self-association. Thus its thermodynamically stable structure in aqueous solution appears to be a highly flexible coil.     

Ligand-induced structural transitions in ErbB receptor extracellular domains

Crystallographic studies showed that epidermal growth factor (EGF) receptor activation involves major domain rearrangements. Without bound ligand, the extracellular region of the receptor (sEGFR) adopts a "tethered" configuration with its dimerization site occluded by apparently autoinhibitory intramolecular interactions. Ligand binding causes the receptor to become "extended," breaking the tether and exposing the dimerization site. Using small-angle X-ray scattering (SAXS), we confirm that the tethered and extended conformations are also adopted in solution, and we describe low-resolution molecular envelopes for an intact sEGFR dimer. We also use SAXS to monitor directly the transition from a tethered to extended configuration in the monomeric extracellular regions of ErbB3 and a dimerization-defective EGFR mutant. Finally, we show that mutating every intramolecular tether interaction in sEGFR does not greatly alter its conformation. These findings explain why tether mutants fail to activate EGF receptor and provide new insight into regulation of ErbB receptor conformation.     

Collagen VI, conformation of A-domain arrays and microfibril architecture

Collagen VI is a ubiquitous extracellular matrix protein that assembles into beaded microfibrils that form networks linking cells to the matrix. Collagen VI microfibrils are typically formed from a heterotrimer of the α1, α2, and α3 chains. The α3 chain is distinct as it contains an extended N terminus with up to 10 consecutive von Willebrand factor type A-domains (VWA). Here, we use solution small angle x-ray scattering (SAXS) and single particle analysis EM to determine the nanostructure of nine of these contiguous A-domains. Both techniques reveal a tight C-shape conformation for the A-domains. Furthermore, using biophysical approaches, we demonstrate that the N-terminal region undergoes a conformational change and a proportion forms dimers in the presence of Zn(2+). This is the first indication that divalent cations interact with collagen VI A-domains. A three-dimensional reconstruction of tissue-purified collagen VI microfibrils was generated using EM and single particle image analysis. The reconstruction showed the intricate architecture of the collagen VI globular regions, in particular the highly structurally conserved C-terminal region and variations in the appearance of the N-terminal region. The N-terminal domains project out from the globular beaded region like angled radial spokes. These could potentially provide interactive surfaces for other cell matrix molecules.     

Structural stability of soybean lipoxygenase-1 in solution as probed by small angle X-ray scattering

Soybean lipoxygenase-1 (LOX-1) is used widely as a model for studying the structural and functional properties of the homologous family of lipoxygenases. The crystallographic structure revealed that LOX-1 is organized in a beta-sheet N-terminal domain and a larger, mostly helical, C-terminal domain. Here, we describe the overall structural characterization of native unliganded LOX-1 in solution, using small angle X-ray scattering (SAXS). We show that the scattering pattern of the unliganded enzyme in solution does not display any significant difference compared with that calculated from the crystal structure, and that models of the overall shape of the protein calculated ab initio from the SAXS pattern provide a close envelope to the crystal structure. These data, demonstrating that LOX-1 has a compact structure also in solution, rule out any major motional flexibility of the LOX-1 molecule in aqueous solutions. In addition we show that eicosatetraynoic acid, an irreversible inhibitor of lipoxygenase used to mimic the effect of substrate binding, does not alter the overall conformation of LOX-1 nor its ability to bind to membranes. In contrast, the addition of glycerol (to 5%, v/v) causes an increase in the binding of the enzyme to membranes without altering its catalytic efficiency towards linoleic acid nor its SAXS pattern, suggesting that the global conformation of the enzyme is unaffected. Therefore, the compact structure determined in the crystal appears to be essentially preserved in these various solution conditions. During the preparation of this article, a paper by M. Hammel and co-workers showed instead a sharp difference between crystal and solution conformations of rabbit 15-LOX-1. The possible cause of this difference might be the presence of oligomers in the rabbit lipoxygenase preparations.     

Contribution of the copper ions in the dinuclear active site to the stability of Carcinus aestuarii hemocyanin

We have investigated the effect of copper binding on the structural properties of hemocyanin (Hc). To this aim, we have studied the holo- and apo-form of the protein, both in the hexameric and in the monomeric state (CaeSS2 subunit), with experimental approaches that report on the protein aggregation and conformational stability. The results of gel-filtration chromatography and small angle X-ray scattering (SAXS) provide evidence that the hydrodynamic and gyration radius (R(g)) of Hc in the hexameric form only slightly increase upon copper removal, whereas a remarkable enhancement in the R(g) value is observed for the CaeSS2 monomer. CD measurements in the far- and near-UV region indicate that removal of copper only marginally affects the conformation of the hexameric Hc. Instead, copper depletion in the CaeSS2 strongly alters the tertiary structure of the monomer (near-UV CD), even though it is almost inconsequential on the secondary structure content (far-UV CD). These findings are fully consistent with the results of limited proteolysis experiments showing that the hexameric Hc is similarly resistant to proteolysis by trypsin both in the holo- and apo-form. Conversely, the apo-form of CaeSS2 monomer is much more susceptible to proteolytic attack by trypsin than the holo-form. Based on SAXS measurements, the concentration-dependent oligomerization process for apo-CaeSS2 has been analyzed on the basis of a thermodynamic model involving a concentration-dependent equilibrium between a monomer in a native-like and an hexameric aggregate of monomers.     

The structure of the ATP-bound state of S. cerevisiae phosphofructokinase determined by cryo-electron microscopy

Phosphofructokinase (Pfk1, EC 2.7.1.11) plays a key regulatory role in the glycolytic pathway. The combination of X-ray crystallographic and biochemical data has provided an understanding of the different conformational changes that occur between the active and inhibited states of the bacterial enzyme, and of the role of the two bacterial effectors. Eukaryotic phosphofructokinases exhibit a far more sophisticated regulatory mechanism, they are more complex structures regulated by a large number of effectors (around 20). Saccharomyces cerevisiae Pfk1 is an 835 kDa hetero-octamer which shows cooperative binding for fructose-6-phosphate (F6P) and non-cooperative binding for ATP. The 3D structure of the F6P-bound state was obtained by cryo-electron microscopy to 1.1 nm resolution. This electron microscopy structure, in combination with molecular replacement using the bacterial enzyme has helped provide initial phases to solve the X-ray structure of the F6P-bound state 12S yeast truncated-tetramer. Biochemical and small-angle X-ray scattering (SAXS) studies had indicated that Pfk1 underwent a large conformational change upon Mg-ATP binding. We have calculated a reconstruction using reference-based 3D projection alignment methods from 0 degrees images acquired from frozen-hydrated preparations of the enzyme in the presence of Mg-ATP. The ATP-bound structure is more extended or open, and the calculated radius of gyration of 7.33 nm (7.0 nm for F6P) is in good agreement with the SAXS data. There is a substantial decrease in the rotational angle between the top and bottom tetramers. Interestingly, all these changes have arisen from a reorientation of the alpha- and beta-subunits in the dimers. The interface region between the alpha- and beta-subunits is now approximately half the size of the one in the F6P-bound structure. This is the first time that the 3D structure of a eukaryotic Pfk1 has been visualized in its T-state (inhibited-state).