SAXSDom: Modeling multidomain protein structures using small-angle X-ray scattering data

Many proteins are composed of several domains that pack together into a complex tertiary structure. Multidomain proteins can be challenging for protein structure modeling, particularly those for which templates can be found for individual domains but not for the entire sequence. In such cases, homology modeling can generate high quality models of the domains but not for the orientations between domains. Small-angle X-ray scattering (SAXS) reports the structural properties of entire proteins and has the potential for guiding homology modeling of multidomain proteins. In this article, we describe a novel multidomain protein assembly modeling method, SAXSDom that integrates experimental knowledge from SAXS with probabilistic Input-Output Hidden Markov model to assemble the structures of individual domains together. Four SAXS-based scoring functions were developed and tested, and the method was evaluated on multidomain proteins from two public datasets. Incorporation of SAXS information improved the accuracy of domain assembly for 40 out of 46 critical assessment of protein structure prediction multidomain protein targets and 45 out of 73 multidomain protein targets from the ab initio domain assembly dataset. The results demonstrate that SAXS data can provide useful information to improve the accuracy of domain-domain assembly. The source code and tool packages are available at https://github.com/jianlin-cheng/SAXSDom .     

Heterologous expression, isotopic-labeling and immuno-characterization of Cin1, a novel protein secreted by the phytopathogenic fungus Venturia inaequalis

The phytopathogenic fungus Venturia inaequalis causes scab of apple. Once this fungus penetrates the plant surface, it forms a specialized body called a stroma between the inner cuticle surface and the epidermal cell wall. A novel V. inaequalis gene, cin1, is strongly up-regulated in the early stages of infection. This gene codes for a 523 residue secreted protein, containing eight imperfect repeats of 60 amino acids. Cin1 was expressed in the methanolytic yeast Pichia pastoris using the pPICZ vector system. A protein of 57 kDa was secreted by these transformants and peptide fingerprinting indicated that it was the Cin1 protein product. Multiple angle laser light scattering confirmed the predicted mass of Cin1, showing it was not glycosylated by Pichia and was monomeric in solution. Through measurements of the hydrodynamic properties of Cin1, the experimental Stokes radius of Cin1 was calculated and corresponded to the theoretical value for a natively folded globular protein of size 57 kDa. The mobility of recombinant Cin1 on native PAGE was also consistent with that of a folded protein. To simplify future structural analyses, a two-domain truncated version, Cin1-2D, consisting of domains one and two, was also expressed using the same vector system. Both proteins were purified to homogeneity. Conditions for maximal (>98%) incorporation of ¹³C and ¹⁵N were determined. A mouse polyclonal antibody and three monoclonal antibodies (MAbs) were raised against the full-length version of Cin1. Analysis of the three MAbs using surface plasmon resonance indicated binding to distinct epitopes on the Cin1 protein. Western blots confirmed the different specificities of each MAb.     

Inter-domain interactions influence the stability and catalytic activity of the bi-domain protein tyrosine phosphatase PTP99A

The two protein tyrosine phosphatase (PTP) domains in bi-domain PTPs share high sequence and structural similarity. However, only one of the two PTP domains is catalytically active. Here we describe biochemical studies on the two tandem PTP domains of the bi-domain PTP, PTP99A. Phosphatase activity, monitored using small molecule as well as peptide substrates, revealed that the inactive (D2) domain activates the catalytic (D1) domain. Thermodynamic measurements suggest that the inactive D2 domain stabilizes the bi-domain (D1-D2) protein. The mechanism by which the D2 domain activates and stabilizes the bi-domain protein is governed by few interactions at the inter-domain interface. In particular, mutating Lys990 at the interface attenuates inter-domain communication. This residue is located at a structurally equivalent location to the so-called allosteric site of the canonical single domain PTP, PTP1B. These observations suggest functional optimization in bi-domain PTPs whereby the inactive PTP domain modulates the catalytic activity of the bi-domain enzyme.     

Homology of the precursor of pulmonary surfactant-associated protein SP-B with prosaposin and sulfated glycoprotein 1.

The precursor of pulmonary surfactant-associated protein, SP-B, is composed of an NH2-terminal domain of 30 residues (a-type domain) and three tandem repeats of about 90 residues (b-type domain); biophysically active mature SP-B corresponds to the second b-type repeat. Consensus sequences constructed for the b-type repeats were used to search the data base for homologous sequences, and the search has revealed that prosaposin and sulfated glycoprotein 1 show a remarkable homology with these repeats. The domain organizations of the latter proteins, however, differ from that of SP-B precursor inasmuch as they contain four tandem copies of the b-type domain and a-type domains are present both in the NH2-terminal and COOH-terminal parts of the proteins. The implications of the homology of saposins and SP-B for their structure and function are discussed.     

Pleiotropic function of intersectin homologue Cin1 in Cryptococcus neoformans

The manifestation of virulence traits in Cryptococcus neoformans is thought to rely on intracellular transport, a process not fully explored in this pathogenic fungus. Through interaction cloning, we identified a multi‐modular protein, Cin1 (cryptococcal intersectin 1), whose domain structure is similar to that of the human endocytic protein ITSN1. Cin1 contains an N‐terminal EH domain, a central coiled‐coil region, a WH2 domain, two SH3 domains and a C‐terminal RhoGEF (DH)‐PH domain. Interestingly, alternative mRNA splicing resulted in two Cin1 isoforms, and Cin1 homologues are also restricted to basidiomycetous fungi. Disruption of the CIN1 gene had a pleiotropic effect on growth, normal cytokinesis, intracellular transports and the production of several virulence factors. Additionally, Cin1 interacts with cryptococcal Cdc42 and Wsp1 (a WASP homologue) proteins in vitro, suggesting a conserved role in the regulation of the actin cytoskeleton. However, deletion of RhoGEF or SH3 and RhoGEF domains did not result in any phenotypic changes, suggesting that functional redundancy exists in proteins containing similar domains or that the activities by other domains are necessary for Cin1 function. Our study presents the first evidence of a multi‐modular protein whose function in intracellular transport underlies the growth, differentiation and virulence of a pathogenic microorganism.     

The structure of the unliganded extracellular domain of the interleukin-6 signal transducer gp130 in solution

Interleukin-6 (IL-6) plays an important role in immune responses and signals via two different pathways. When IL-6 binds to its non-signalling membrane-bound receptor (IL-6R), a non-covalent dimer of the ubiquitous interleukin-6 signal transducer gp130 is recruited to initiate intracellular signalling cascades. This so-called classical signalling pathway is restricted to cells expressing the membrane-bound IL-6R, such as hepatocytes and certain leukocytes. In addition, an alternative trans-signalling pathway uses soluble forms of IL-6R (sIL-6R) in complex with IL-6 to activate cells expressing gp130, but not membrane-bound IL-6R. In both cases, a tetrameric or hexameric signalling complex consisting of two gp130 molecules and one or two molecules each of IL-6 and (s)IL-6R is formed. The structure of the hexameric complex of the ligand-binding domains of gp130 (D1-D3) with IL-6 and sIL-6R has been solved by X-ray crystallography as well as the full-length extracellular part of gp130 (D1-D6) as a monomer. Since gp130 exists as a preformed dimer on the cell surface, we used a sgp130Fc fusion protein - consisting of two extracellular gp130 regions (D1-D6) dimerised by an IgG1-Fc part - to study the structure of unliganded gp130 extracellular domains in solution by small-angle X-ray scattering (SAXS). The SAXS data indicated that sgp130Fc forms a rigid molecule in solution. The low resolution structural model reveals an elongated assembly with an Fc base and two gp130 arms, whereby the orientation of the ligand-binding domains D1-D3 with respect to the membrane-proximal domains D4-D6 differs from that in the crystallographic monomer. Functional implications of these findings are discussed.     

Identification of residues involved in the interaction of Staphylococcus aureus fibronectin-binding protein with the (4)F1(5)F1 module pair of human fibronectin using heteronuclear NMR spectroscopy

Many pathogenic Gram-positive bacteria express cell surface proteins that bind to components of the extracellular matrix. This paper describes studies of the interaction between ligand binding repeats (D3 and D1-D4) of a fibronectin-binding protein from Staphylococcus aureus with a module pair ((4)F1(5)F1) from the N-terminal region of fibronectin. When D3 was added to isotope-labeled (4)F1(5)F1, (1)H, (15)N, and (13)C NMR chemical shift changes indicate that binding is primarily via residues in (4)F1, although a few residues in (5)F1 are also affected. Both hydrophobic and electrostatic interactions appear to be involved. The NMR data indicate that part of the D3 repeat converts from a disordered to a more ordered, extended conformation on binding to (4)F1(5)F1. In further NMR experiments, selective reduction of the intensity of D1-D4 resonances was observed on binding to (4)F1(5)F1, consistent with previous suggestions that in each of D1, D2, and D3 repeats, the main fibronectin binding site is in the C-terminal region of the repeat. In D1-D4, these regions also appear to go from a disordered to a more ordered conformation of fibronectin binding. Although the regions of the two proteins which interact had been previously identified, the findings presented here identify, for the first time, the specific residues in both proteins that are likely to be involved in the interaction.     

Structure and evolution of the Cinful retrotransposon family of maize.

A maize cDNA clone was isolated by virtue of its intense hybridization to total maize genomic DNA, indicating homology to highly repetitive sequences. Genomic homologues were identified and subcloned from an adh1-bearing maize yeast artificial chromosome (YAC). Sequencing revealed that the expressed sequence was part of a Ty3-gypsy-type retrotransposon. We discovered and sequenced two complete retrotransposons of this family, and named them Cinful elements because they are members of a family of maize retrotransposons including Zeon-1 and the first plant transposable element sequenced, the solo long terminal repeat (LTR) called Cin1. All are defective, as Cinful-1 and Cinful-2 elements lack gag and Zeon-1 lacks pol homology. Despite the apparent lack of an intact "autonomous" element, the Cinful family has expanded to a copy number of about 18 000, representing just under 9% of the maize genome. Both point mutations and major rearrangements, including possible gene acquisition, differentiate members of the Cinful family. Cinful family members were found to have an unusual feature that we also observed in two other Ty3-class retrotransposons of teosinte and tobacco: related tandem repeats that separate their internal domains with a gag- or pol-containing homology from a 3' segment of unknown function. The conserved and variable features identified provide insights into the origin, mutational history, and functional components of this major constituent of the maize genome.     

The trypanosome leucine repeat gene in the variant surface glycoprotein expression site encodes a putative metal-binding domain and a region resembling protein-binding domains of yeast, Drosophila, and mammalian proteins.

We have identified a new variant surface glycoprotein expression site-associated gene (ESAG) in Trypanosoma brucei, the trypanosome leucine repeat (T-LR) gene. Like most other ESAGs, it is expressed in a life cycle stage-specific manner. The N-terminal 20% of the predicted T-LR protein resembles the metal-binding domains of nucleic acid-binding proteins. The remainder is composed of leucine-rich repeats that are characteristic of protein-binding domains found in a variety of other eucaryote proteins. This is the first report of leucine-rich repeats and potential nucleic acid-binding domains on the same protein. The T-LR gene is adjacent to ESAG 4, which has homology to the catalytic domain of adenylate cyclase. This is intriguing, since yeast adenylate cyclase has a leucine-rich repeat regulatory domain. The leucine-rich repeat and putative metal-binding domains suggest a possible regulatory role that may involve adenylate cyclase activity or nucleic acid binding.     

Conformational rearrangement within the soluble domains of the CD4 receptor is ligand-specific

Ligand binding induces shape changes within the four modular ectodomains (D1-D4) of the CD4 receptor, an important receptor in immune signaling. Small angle x-ray scattering (SAXS) on both a two-domain and a four-domain construct of the soluble CD4 (sCD4) is consistent with known crystal structures demonstrating a bilobal and a semi-extended tetralobal Z conformation in solution, respectively. Detection of conformational changes within sCD4 as a result of ligand binding was followed by SAXS on sCD4 bound to two different glycoprotein ligands: the tick saliva immunosuppressor Salp15 and the HIV-1 envelope protein gp120. Ab initio modeling of these data showed that both Salp15 and gp120 bind to the D1 domain of sCD4 and yet induce drastically different structural rearrangements. Upon binding, Salp15 primarily distorts the characteristic lobal architecture of the sCD4 without significantly altering the semi-extended shape of the sCD4 receptor. In sharp contrast, the interaction of gp120 with sCD4 induces a shape change within sCD4 that can be described as a Z-to-U bi-fold closure of the four domains across its flexible D2-D3 linker. Placement of known crystal structures within the boundaries of the SAXS-derived models suggests that the ligand-induced shape changes could be a result of conformational changes within this D2-D3 linker. Functionally, the observed shape changes in CD4 receptor causes dissociation of lymphocyte kinase from the cytoplasmic domain of Salp15-bound CD4 and facilitates an interaction between the exposed V3 loops of CD4-bound gp120 molecule to the extracellular loops of its co-receptor, a step essential for HIV-1 viral entry.     

Detection of breaking points in helices linking separate domains

A novel method is proposed to predict whether two domains connected by a helical link can mutually reorient themselves, as well as where the helix can be distorted to allow the domain-domain movements. The method, based on analysis of the variation of the a.d.p. values along the helix link, is applied to three proteins--calmodulin, lysozyme, and hemagglutinin--for which both the domain-domain flexibility and the helix fragment responsible for it are well documented. The helix regions that are variously distorted to permit domain-domain reorientation are well predicted. The method is also applied to colicin Ia and shows that an inter-domain rearrangement can take place as previously postulated. The prediction of the helix breaking point should prove useful in interpreting structural data and in defining the domain borders automatically for proteins built by domains connected by helical links.     

Two novel Venturia inaequalis genes induced upon morphogenetic differentiation during infection and in vitro growth on cellophane

Venturia inaequalis is a hemibiotrophic ascomycete that causes apple scab. Germ tubes, from conidia or ascospores, penetrate the leaf or fruit surface directly via appressoria-like swellings; subsequently the hyphae divide laterally to form a stroma between the cuticle and the outer wall of the epidermal cells. This morphological switch can be mimicked by growing the fungus in vitro on cellophane discs. The aim of this work was to identify genes upregulated in planta using growth on cellophane as a model. Four cDNA clones were found to be induced by growth on cellophane, and qRT-PCR showed two of these genes were up-regulated over a thousand fold in infected apple leaves compared to liquid culture. The predicted proteins for both genes possess putative signal peptides for secretion but have no similarity to sequences in publicly available databases. Both genes encode proteins with novel, imperfect repeat domain structures, the number of which vary in an isolate-specific fashion. Cin1 has seven or eight repeats of about 60 amino acids with four conserved cysteine residues per repeat, while Cin3 has four or five repeats of 32 amino acids with no cysteines. Both proteins appear to have evolved through internal duplication. Cin3, in particular, shows considerable between-strain variation in domain structure, indicating a high degree of recombination at this locus and revealing that the repeat structure has most likely arisen by unequal crossing-over. Results of this study support the hypothesis that cellophane-grown V. inaequalis mimics aspects of biotrophic infection and provide the first insights into novel fungal genes expressed during apple scab infection and their mechanisms of evolution.     

A simple genetic algorithm for the optimization of multidomain protein homology models driven by NMR residual dipolar coupling and small angle <Emphasis Type="Italic">X</Emphasis>-ray scattering data

Most proteins comprise several domains and/or participate in functional complexes. Owing to ongoing structural genomic projects, it is likely that it will soon be possible to predict, with reasonable accuracy, the conserved regions of most structural domains. Under these circumstances, it will be important to have methods, based on simple-to-acquire experimental data, that allow to build and refine structures of multi-domain proteins or of protein complexes from homology models of the individual domains/proteins. It has been recently shown that small angle X-ray scattering (SAXS) and NMR residual dipolar coupling (RDC) data can be combined to determine the architecture of such objects when the X-ray structures of the domains are known and can be considered as rigid objects. We developed a simple genetic algorithm to achieve the same goal, but by using homology models of the domains considered as deformable objects. We applied it to two model systems, an S1KH bi-domain of the NusA protein and the γS-crystallin protein. Despite its simplicity our algorithm is able to generate good solutions when driven by SAXS and RDC data.     

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.     

Assembly and routing of von Willebrand factor variants: the requirements for disulfide-linked dimerization reside within the carboxy-terminal 151 amino acids

The precursor protein of von Willebrand factor (pro-vWF) consists of four different repeated domains, denoted D1-D2-D'-D3-A1-A2-A3-D4-B1-B2-B3-C1-C2, followed by a carboxy-terminal region of 151 amino acids without obvious internal homology. Previously, we have shown the requirement of the domains D1, D2, D', and D3 of pro-vWF in the assembly of pro-vWF dimers into multimers. Here, we define the domains of vWF involved in dimerization, using deletion mutants of full-length vWF cDNA transiently expressed in monkey kidney COS-1 cells. It is shown that only the carboxy-terminal 151 amino acid residues of vWF are required for dimerization. In addition, by analyzing a construct, encoding only the carboxy-terminal 151 amino acids of vWF, we find that the formation of dimers is an event independent of other domains present on pro-vWF, such as the domains C1 and C2 previously suggested to be involved in dimerization. Furthermore, it is shown that a deletion mutant of vWF, lacking the carboxy-terminal 151 amino acid residues and thus unable to dimerize, is proteolytically degraded in the ER. In contrast, a mutant protein, composed only of the carboxy-terminal 151 amino acids of vWF, and able to dimerize, is transported from the ER in a similar fashion as wild-type vWF. The role of the ER in the assembly of vWF is discussed with regard to the data presented in this paper on the intracellular fate of several vWF mutant proteins.     

Structure of a putative ancestral protein encoded by a single sequence repeat from a multidomain proteinase inhibitor gene from Nicotiana alata.

BACKGROUND: The ornamental tobacco Nicotiana alata produces a series of proteinase inhibitors (PIs) that are derived from a 43 kDa precursor protein, NaProPI. NaProPI contains six highly homologous repeats that fold to generate six separate structural domains, each corresponding to one of the native PIs. An unusual feature of NaProPI is that the structural domains lie across adjacent repeats and that the sixth PI domain is generated from fragments of the first and sixth repeats. Although the homology of the repeats suggests that they may have arisen from gene duplication, the observed folding does not appear to support this. This study of the solution structure of a single NaProPI repeat (aPI1) forms a basis for unravelling the mechanism by which this protein may have evolved. RESULTS: The three-dimensional structure of aPI1 closely resembles the triple-stranded antiparallel beta sheet observed in each of the native PIs. The five-residue sequence Glu-Glu-Lys-Lys-Asn, which forms the linker between the six structural domains in NaProPI, exists as a disordered loop in aPI1. The presence of this loop in aPI1 results in a loss of the characteristically flat and disc-like topography of the native inhibitors. CONCLUSIONS: A single repeat from NaProPI is capable of folding into a compact globular domain that displays native-like PI activity. Consequently, it is possible that a similar single-domain inhibitor represents the ancestral protein from which NaProPI evolved.     

The A5 antigen, a candidate for the neuronal recognition molecule, has homologies to complement components and coagulation factors

The A5 antigen is a neuronal cell surface protein of Xenopus presumed to be involved in the neuronal recognition between the optic nerve fibers and the visual centers. Analyses of cDNA clones revealed that the A5 antigen is a class I membrane protein containing two different internal repeats in the extracellular segment. The first repeat bears homology to domain III of complement components C1r and C1s, and the second repeat is homologous to the C1 and C2 domains of coagulation factors V and VIII. The mRNA for the A5 antigen was present in retinal ganglion cells and visual center neurons. Nonneuronal cells in the peripheral and central nervous systems did not express the mRNA for the A5 antigen.     

Domains involved in multimer assembly of von willebrand factor (vWF): multimerization is independent of dimerization.

The precursor protein of von Willebrand factor (pro-vWF) consist of four repeated domains, denoted D1-D2-D'-D3-A1-A2-A3-D4-B1-B2-B3-C1-C2. The domains D1 and D2 constitute the amino-terminal pro-polypeptide and the remaining domains mature vWF, generated upon proteolytic processing. We have shown previously that the pro-polypeptide of pro-vWF is obligatory for assembly of pro-vWF dimers into multimers, a process vital for efficient adhesion of platelets to an injured vessel wall. Here, we have employed full length vWF cDNA to construct a series of deletion mutants, based on the homology between the various domains. Specifically, the domains D', D3 or both were deleted and the multimeric pattern of the mutant vWF proteins was analysed after transient expression in COS-1 cells. It is demonstrated that in addition to the pro-polypeptide, both the D' and the D3 domain are required for multimer assembly. Furthermore, by analysing a construct containing only the domains D' and D3 next to the pro-polypeptide it is shown that this is the only part of the vWF protein involved in multimer assembly. Since, the formation of pro-vWF dimers relies on the carboxy-terminal area of mature vWF, it is concluded that multimer assembly is a process independent of dimerization.     

A helical string of alternately connected three-helix bundles for the cell wall-associated adhesion protein Ebh from Staphylococcus aureus

The 1.1 MDa cell-wall-associated adhesion protein of staphylococci, Ebh, consists of several distinct regions, including a large central region with 52 imperfect repeats of 126 amino acid residues. We investigated the structure of this giant molecule by X-ray crystallography, circular dichroism (CD) spectrometry, and small-angle X-ray scattering (SAXS). The crystal structure of two repeats showed that each repeat consists of two distinct three-helix bundles, and two such repeats are connected along the long axis, resulting in a rod-like structure that is 120 A in length. CD and SAXS analyses of the samples with longer repeats suggested that each repeat has an identical structure, and that such repeats are connected tandemly to form a rod-like structure in solution, the length of which increased proportionately with the number of repeating units. On the basis of these results, it was proposed that Ebh is a 320 nm rod-like molecule with some plasticity at module junctions.     

Prediction of the structure of human Janus kinase 2 (JAK2) comprising JAK homology domains 1 through 7

A theoretical model of human Janus kinase 2 (JAK2) comprising all seven Janus homology domains is presented. The model was generated by application of homology modelling approaches. The three-dimensional structure contains, starting from the N-terminus, FERM (4.1, ezrin, radixin, moesin), SH2 (Src homology region 2), tyrosine kinase-like, and tyrosine kinase domains. The predicted inter-domain orientation in JAK2 is discussed and the currently existing mutational data for Janus kinases are evaluated. Structural details of the SH2 and the FERM domains are presented. The predictions indicate that the SH2 domain is not fully functional. A number of hydrophobic amino acids of the FERM domain that are predicted to be involved in the constitutive association with the cytokine receptors are highlighted. The model gives new insights into the structure-function relationship of this important protein, and areas that could be investigated by mutation studies are highlighted.