Molecular modeling of the multidomain structures of the proteoglycan binding region and the link protein of cartilage by neutron and synchrotron X-ray scattering

The interaction of proteoglycan monomers with hyaluronate in cartilage is mediated by a globular binding region at the N-terminus of the proteoglycan monomer; this interaction is stabilized by link protein. Sequences show that both the binding region (27% carbohydrate) and the link protein (6% carbohydrate) contain an immunoglobulin (Ig) fold domain and two proteoglycan tandem repeat (PTR) domains. Both proteins were investigated by neutron and synchrotron X-ray solution scattering, in which nonspecific aggregate formation was reduced by the use of citraconylation to modify surface lysine residues. The neutron and X-ray radius of gyration RG of native and citraconylated binding region is 5.1 nm, and the cross-sectional RG (RXS) is 1.9-2.0 nm. No neutron contrast dependence of the RG values was observed; however, a large contrast dependence was seen for the RXS values which is attributed to the high carbohydrate content of the binding region. The neutron RG for citraconylated link protein is 2.9 nm, its RXS is 0.8 nm, and these data are also independent of the neutron contrast. The scattering curves of binding region and link protein were modeled using small spheres. Both protein structures were defined initially by the representation of one domain by a crystal structure for a variable Ig fold and a fixed volume for the two PTR domains calculated from sequence data. The final models showed that the different dimensions and neutron contrast properties of binding region compared to link protein could be attributed to an extended glycosylated C-terminal peptide with extended carbohydrate structures in the binding region.(ABSTRACT TRUNCATED AT 250 WORDS)     

The immunoglobulin fold family: sequence analysis and 3D structure comparisons.

Fifty-two 3D structures of Ig-like domains covering the immunoglobulin fold family (IgFF) were compared and classified according to the conservation of their secondary structures. Members of the IgFF are distantly related proteins or evolutionarily unrelated proteins with a similar fold, the Ig fold. In this paper, a multiple structural alignment of the conserved common core is described and the correlation between corresponding sequences is discussed. While the members of the IgFF exhibit wide heterogeneity in terms of tissue and species distribution or functional implications, the 3D structures of these domains are far more conserved than their sequences. We define topologically equivalent residues in the Ig-like domains, describe the hydrophobic common cores and discuss the presence of additional strands. The disulfide bridges, not necessary for the stability of the Ig fold, may have an effect on the compactness of the domains. Based upon sequence and structure analysis, we propose the introduction of two new subtypes (C3 and C4) to the previous classifications, in addition to a new global structural classification. The very low mean sequence identity between subgroups of the IgFF suggests the occurrence of both divergent and convergent evolutionary processes, explaining the wide diversity of the superfamily. Finally, this review suggest that hydrophobic residues constituting the common hydrophobic cores are important clues to explain how highly divergent sequences can adopt a similar fold.     

Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.     

Mechanical stability and differentially conserved physical–chemical properties of titin Ig‐domains

The mechanisms that determine mechanical stabilities of protein folds remain elusive. Our understanding of these mechanisms is vital to both bioengineering efforts and to the better understanding and eventual treatment of pathogenic mutations affecting mechanically important proteins such as titin. We present a new approach to analyze data from single‐molecule force spectroscopy for different domains of the giant muscle protein titin. The region of titin found in the I‐band of a sarcomere is composed of about 40 Ig‐domains and is exposed to force under normal physiological conditions and connects the free‐hanging ends of the myosin filaments to the Z‐disc. Recent single‐molecule force spectroscopy data show a mechanical hierarchy in the I‐band domains. Domains near the C‐terminus in this region unfold at forces two to three times greater than domains near the beginning of the I‐band. Though all of these Ig‐domains are thought to share a fold and topology common to members of the Ig‐like fold family, the sequences of neighboring domains vary greatly with an average sequence identity of only 25%. We examine in this study the relation of these unique mechanical stabilities of each I‐band Ig domain to specific, conserved physical–chemical properties of amino acid sequences in related Ig domains. We find that the sequences of each individual titin Ig domain are very highly conserved, with an average sequence identity of 79% across species that are divergent as humans, chickens, and zebra fish. This indicates that the mechanical properties of each domain are well conserved and tailored to its unique position in the titin molecule. We used the PCPMer software to determine the conservation of amino acid properties in titin Ig domains grouped by unfolding forces into “strong” and “weak” families. We found two motifs unique to each family that may have some role in determining the mechanical properties of these Ig domains. A detailed statistical analysis of properties of individual residues revealed several positions that displayed differentially conserved properties in strong and weak families. In contrast to previous studies, we find evidence that suggests that the mechanical stability of Ig domains is determined by several residues scattered across the β‐sandwich fold, and force sensitive residues are not only confined to the A′‐G region. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Structure of the C-terminal domain from Trypanosoma brucei variant surface glycoprotein MITat1.2

The variant surface glycoprotein (VSG) of African trypanosomes has a structural role in protecting other cell surface proteins from effector molecules of the mammalian immune system and also undergoes antigenic variation necessary for a persistent infection in a host. Here we have reported the solution structure of a VSG type 2 C-terminal domain from MITat1.2, completing the first structure of both domains of a VSG. The isolated C-terminal domain is a monomer in solution and forms a novel fold, which commences with a short alpha-helix followed by a single turn of 3(10)-helix and connected by a short loop to a small anti-parallel beta-sheet and then a longer alpha-helix at the C terminus. This compact domain is flanked by two unstructured regions. The structured part of the domain contains 42 residues, and the core comprises 2 disulfide bonds and 2 hydrophobic residues. These cysteines and hydrophobic residues are conserved in other VSGs, and we have modeled the structures of two further VSG C-terminal domains using the structure of MITat1.2. The models suggest that the overall structure of the core is conserved in the different VSGs but that the C-terminal alpha-helix is of variable length and depends on the presence of charged residues. The results provided evidence for a conserved tertiary structure for all the type 2 VSG C-terminal domains, indicated that VSG dimers form through interactions between N-terminal domains, and showed that the selection pressure for sequence variation within a conserved tertiary structure acts on the whole of the VSG molecule.     

Identification of cDNA clones encoding different domains of the basement membrane heparan sulfate proteoglycan

We have used antibodies to the basement membrane proteoglycan to screen lambda gt11 expression vector libraries and have isolated two cDNA clones, termed BPG 5 and BPG 7, which encode different portions of the core protein of the heparan sulfate basement membrane proteoglycan. These clones hybridize to a single mRNA species of approximately 12 kilobases. Amino acid sequences obtained on peptides derived from protease digests of the core protein were found in the deduced sequence, confirming the identity of these clones. BPG 5 spanned 1986 base pairs and has an open reading frame of 662 amino acids. The amino acid sequence deduced from BPG 5 contains two cysteine-rich domains and two internally homologous domains lacking cysteine. The cysteine-rich domains show homology to the cysteine-rich domains of the laminin chains. A globule-rod structure, similar to that of the short arms of the laminin chains, is proposed for this region of the proteoglycan. The other clone, BPG 7, is 2193 base pairs long and has an open reading frame of 731 amino acids. The deduced sequence contains eight internal repeats with 2 cysteine residues in each repeat. These repeats show homology to the neural-cell adhesion molecule N-CAM and the plasma alpha 1B-glycoprotein. Looping structures similar to these proteins and to other proteins of the immunoglobulin gene superfamily are proposed for this region of the proteoglycan. The sequence DSGEY was found four times in this domain and could be heparan sulfate attachment sites.     

Conformational analysis of peptides corresponding to all the secondary structure elements of protein L B1 domain: secondary structure propensities are not conserved in proteins with the same fold.

The solution conformation of three peptides corresponding to the two beta-hairpins and the alpha-helix of the protein L B1 domain have been analyzed by circular dichroism (CD) and nuclear magnetic resonance spectroscopy (NMR). In aqueous solution, the three peptides show low populations of native and non-native locally folded structures, but no well-defined hairpin or helix structures are formed. In 30% aqueous trifluoroethanol (TFE), the peptide corresponding to the alpha-helix adopts a high populated helical conformation three residues longer than in the protein. The hairpin peptides aggregate in TFE, and no significant conformational change occurs in the NMR observable fraction of molecules. These results indicate that the helical peptide has a significant intrinsic tendency to adopt its native structure and that the hairpin sequences seem to be selected as non-helical. This suggests that these sequences favor the structure finally attained in the protein, but the contribution of the local interactions alone is not enough to drive the formation of a detectable population of native secondary structures. This pattern of secondary structure tendencies is different to those observed in two structurally related proteins: ubiquitin and the protein G B1 domain. The only common feature is a certain propensity of the helical segments to form the native structure. These results indicate that for a protein to fold, there is no need for large native-like secondary structure propensities, although a minimum tendency to avoid non-native structures and to favor native ones could be required.     

Structural analysis of the PsbQ protein of photosystem II by Fourier transform infrared and circular dichroic spectroscopy and by bioinformatic methods

The structure of PsbQ, one of the three main extrinsic proteins associated with the oxygen-evolving complex (OEC) of higher plants and green algae, is examined by Fourier transform infrared (FTIR) and circular dichroic (CD) spectroscopy and by computational structural prediction methods. This protein, together with two other lumenally bound extrinsic proteins, PsbO and PsbP, is essential for the stability and full activity of the OEC in plants. The FTIR spectra obtained in both H(2)O and D(2)O suggest a mainly alpha-helix structure on the basis of the relative areas of the constituents of the amide I and I' bands. The FTIR quantitative analyses indicate that PsbQ contains about 53% alpha-helix, 7% turns, 14% nonordered structure, and 24% beta-strand plus other beta-type extended structures. CD analyses indicate that PsbQ is a mainly alpha-helix protein (about 64%), presenting a small percentage assigned to beta-strand ( approximately 7%) and a larger amount assigned to turns and nonregular structures ( approximately 29%). Independent of the spectroscopic analyses, computational methods for protein structure prediction of PsbQ were utilized. First, a multiple alignment of 12 sequences of PsbQ was obtained after an extensive search in the public databases for protein and EST sequences. Based on this alignment, computational prediction of the secondary structure and the solvent accessibility suggest the presence of two different structural domains in PsbQ: a major C-terminal domain containing four alpha-helices and a minor N-terminal domain with a poorly defined secondary structure enriched in proline and glycine residues. The search for PsbQ analogues by fold recognition methods, not based on the secondary structure, also indicates that PsbQ is a four alpha-helix protein, most probably folding as an up-down bundle. The results obtained by both the spectroscopic and computational methods are in agreement, all indicating that PsbQ is mainly an alpha protein, and show the value of using both methodologies for protein structure investigation.     

Designed protein G core variants fold to native-like structures: sequence selection by ORBIT tolerates variation in backbone specification

The solution structures of two computationally designed core variants of the beta 1 domain of streptococcal protein G (G beta 1) were solved by (1)H NMR methods to assess the robustness of amino acid sequence selection by the ORBIT protein design package under changes in protein backbone specification. One variant has mutations at three of 10 core positions and corresponds to minimal perturbations of the native G beta 1 backbone. The other, with mutations at six of 10 positions, was calculated for a backbone in which the separation between G beta 1's alpha-helix and beta-sheet was increased by 15% relative to native G beta 1. Exchange broadening of some resonances and the complete absence of others in spectra of the sixfold mutant bespeak conformational heterogeneity in this protein. The NMR data were sufficiently abundant, however, to generate structures of similar, moderately high quality for both variants. Both proteins adopt backbone structures similar to their target folds. Moreover, the sequence selection algorithm successfully predicted all core chi(1) angles in both variants, five of six chi(2) angles in the threefold mutant and four of seven chi(2) angles in the sixfold mutant. We conclude that ORBIT calculates sequences that fold specifically to a geometry close to the template, even when the template is moderately perturbed relative to a naturally occurring structure. There are apparently limits to the size of acceptable perturbations: In this study, the larger perturbation led to undesired dynamic behavior.     

Relationship between sequence determinants of stability for two natural homologous proteins with different folds

In the Cro protein family, an evolutionary change in secondary structure has converted an alpha-helical fold to a mixture of alpha-helix and beta-sheet. P22 Cro and lambda Cro represent the ancestral all-alpha and descendant alpha+beta folds, respectively. The major structural differences between these proteins are at the C-terminal end of the domain (residues 34-56), where two alpha-helices in P22 Cro align with two beta-strands in lambda Cro. We sought to assess the possibility that smooth evolutionary transitions could have converted the all-alpha structure to the alpha+beta structure through sequences that could adopt both folds. First, we used scanning mutagenesis to identify and compare patterns of key stabilizing residues in the C-terminal regions of both P22 Cro and lambda Cro. These patterns exhibited little similarity to each other, with structurally important residues in the two proteins most often occurring at different sequence positions. Second, "hybrid scanning" studies, involving replacement of each wild-type residue in P22 Cro with the aligned wild-type residue in lambda Cro and vice versa, revealed five or six residues in each protein that strongly destabilized the other. These results suggest that key stability determinants for each Cro fold are quite different and that the P22 Cro sequence strongly favors the all-alpha structure while the lambda Cro sequence strongly favors the alpha+beta structure. Nonetheless, we were able to design a "structurally ambivalent" sequence fragment (SASF1), which corresponded to residues 39-56 and simultaneously incorporated most key stabilizing residues for both P22 Cro and lambda Cro. NMR experiments showed SASF1 to stably fold as a beta-hairpin when incorporated into the lambda Cro sequence but as a pair of alpha-helices when incorporated into P22 Cro.     

Identification of novel restriction endonuclease-like fold families among hypothetical proteins

Restriction endonucleases and other nucleic acid cleaving enzymes form a large and extremely diverse superfamily that display little sequence similarity despite retaining a common core fold responsible for cleavage. The lack of significant sequence similarity between protein families makes homology inference a challenging task and hinders new family identification with traditional sequence-based approaches. Using the consensus fold recognition method Meta-BASIC that combines sequence profiles with predicted protein secondary structure, we identify nine new restriction endonuclease-like fold families among previously uncharacterized proteins and predict these proteins to cleave nucleic acid substrates. Application of transitive searches combined with gene neighborhood analysis allow us to confidently link these unknown families to a number of known restriction endonuclease-like structures and thus assign folds to the uncharacterized proteins. Finally, our method identifies a novel restriction endonuclease-like domain in the C-terminus of RecC that is not detected with structure-based searches of the existing PDB database.     

Solution structure of the hypothetical protein TA0095 from Thermoplasma acidophilum: a novel superfamily with a two-layer sandwich architecture

TA0095 is a 96-residue hypothetical protein from Thermoplasma acidophilum that exhibits no sequence similarity to any protein of known structure. Also, TA0095 is a member of the COG4004 orthologous group of unknown function found in Archaea bacteria. We determined its three-dimensional structure by NMR methods. The structure displays an alpha/beta two-layer sandwich architecture formed by three alpha-helices and five beta-strands following the order beta1-alpha1-beta2-beta3-beta4-beta5-alpha2-alpha3. Searches for structural homologs indicate that the TA0095 structure belongs to the TBP-like fold, constituting a novel superfamily characterized by an additional C-terminal helix. The TA0095 structure provides a fold common to the COG4004 proteins that will obviously belong to this new superfamily. Most hydrophobic residues conserved in the COG4004 proteins are buried in the structure determined herein, thus underlying their importance for structure stability. Considering that the TA0095 surface shows a large positively charged patch with a high degree of residue conservation within the COG4004 domain, the biological function of TA0095 and the rest of COG4004 proteins might occur through binding a negatively charged molecule. Like other TBP-like fold proteins, the COG4004 proteins might be DNA-binding proteins. The fact that TA0095 is shown to interact with large DNA fragments is in favor of this hypothesis, although nonspecific DNA binding cannot be ruled out.     

Modularity and homology: modelling of the titin type I modules and their interfaces

Titin is a giant muscle protein with a highly modular architecture consisting of multiple repeats of two sequence motifs, named type I and type II. Type I motifs are homologous to members of the fibronectin type 3 (Fn3) superfamily, one of the motifs most widespread in modular proteins. Fn3 domains are thought to mediate protein-protein interactions and to act as spacers. In titin, Fn3 modules are present in two different super-repeated patterns, likely to be involved in sarcomere assembly through interactions with A-band proteins. Here, we discuss results from homology modelling the whole family of Fn3 domains in titin. Homology modelling is a powerful tool that will play an increasingly important role in the post-genomic era. It is particularly useful for extending experimental structure determinations of parts of multidomain proteins that contain multiple copies of the same motif. The 3D structures of a representative titin type I domain and of other extracellular Fn3 modules were used as a template to model the structures of the 132 copies in titin. The resulting models suggest residues that contribute to the fold stability and allow us to distinguish these from residues likely to have functional importance. In particular, analysis of the models and mapping of the consensus sequence onto the 3D structure suggest putative surfaces of interaction with other proteins. From the structures of isolated modules and the pattern of conservation in the multiple alignment of the whole titin Ig and Fn3 families, it is possible to address the question of how tandem modules are assembled. Our predictions can be validated experimentally.     

Similarities and differences among 105 members of the Int family of site-specific recombinases.

Alignments of 105 site-specific recombinases belonging to the Int family of proteins identified extended areas of similarity and three types of structural differences. In addition to the previously recognized conservation of the tetrad R-H-R-Y, located in boxes I and II, several newly identified sequence patches include charged amino acids that are highly conserved and a specific pattern of buried residues contributing to the overall protein fold. With some notable exceptions, unconserved regions correspond to loops in the crystal structures of the catalytic domains of lambda Int (Int c170) and HP1 Int (HPC) and of the recombinases XerD and Cre. Two structured regions also harbor some pronounced differences. The first comprises beta-sheets 4 and 5, alpha-helix D and the adjacent loop connecting it to alpha-helix E: two Ints of phages infecting thermophilic bacteria are missing this region altogether; the crystal structures of HPC, XerD and Cre reveal a lack of beta-sheets 4 and 5; Cre displays two additional beta-sheets following alpha-helix D; five recombinases carry large insertions. The second involves the catalytic tyrosine and is seen in a comparison of the four crystal structures. The yeast recombinases can theoretically be fitted to the Int fold, but the overall differences, involving changes in spacing as well as in motif structure, are more substantial than seen in most other proteins. The phenotypes of mutations compiled from several proteins are correlated with the available structural information and structure-function relationships are discussed. In addition, a few prokaryotic and eukaryotic enzymes with partial homology with the Int family of recombinases may be distantly related, either through divergent or convergent evolution. These include a restriction enzyme and a subgroup of eukaryotic RNA helicases (D-E-A-D proteins).     

Genes encoding three new members of the leukocyte antigen 6 superfamily and a novel member of Ig superfamily, together with genes encoding the regulatory nuclear chloride ion channel protein (hRNCC) and an N omega-N omega-dimethylarginine dimethylaminohydrolase homologue, are found in a 30-kb segment of the MHC class III region.

Many of the genes in the class III region of the human MHC encode proteins involved in the immune and inflammatory responses. We have sequenced a 30-kb segment of the MHC class III region lying between the heat shock protein 70 and TNF genes as part of a program aimed at identifying genes that could be involved in autoimmune disease susceptibility. The sequence analysis has revealed the localization of seven genes, whose precise position and order is cen-G7-G6-G6A-G6B-G6C-G6D-G6E-tel, five of which are fully encoded in the sequence, allowing their genomic structures to be defined. Three of them (G6C, G6D, and G6E) encode putative proteins that belong to the Ly-6 superfamily, known to be GPI-anchored proteins attached to the cell surface. Members of the family are specifically expressed and are important in leukocyte maturation. A fourth gene, G6B, encodes a novel member of the Ig superfamily containing a single Ig V-like domain and a cytoplasmic tail with several signal transduction features. The G6 gene encodes a regulatory nuclear chloride ion channel protein, while the G6A gene encodes a putative homologue of the enzyme N omega,N omega-dimethylarginine dimethylaminohydrolase, which is thought to be involved in regulating nitric oxide synthesis. In addition, three microsatellite markers, 9N-1, 82-2, and D6S273 are contained within the sequence, the last two of which have been reported to be strongly associated with the autoimmune disease ankylosing spondylitis.     

A simple motif for protein recognition in DNA secondary structures

DNA in a single-stranded form (ssDNA) exists transiently within the cell and comprises the telomeres of linear chromosomes and the genomes of some DNA viruses. As with RNA, in the single-stranded state, some DNA sequences are able to fold into complex secondary and tertiary structures that may be recognized by proteins and participate in gene regulation. To better understand how such DNA elements might fold and interact with proteins, and to compare recognition features to those of a structured RNA, we used in vitro selection to identify ssDNAs that bind an RNA-binding peptide from the HIV Rev protein with high affinity and specificity. The large majority of selected binders contain a non-Watson-Crick G.T base-pair and an adjacent C:G base-pair and both are essential for binding. This GT motif can be presented in different DNA contexts, including a nearly perfect duplex and a branched three-helix structure, and appears to be recognized in large part by arginine residues separated by one turn of an alpha-helix. Interestingly, a very similar GT motif is necessary also for protein binding and function of a well-characterized model ssDNA regulatory element from the proenkephalin promoter.     

Crystal structure of the major periplasmic domain of the bacterial membrane protein assembly facilitator YidC

The essential bacterial membrane protein YidC facilitates insertion and assembly of proteins destined for integration into the inner membrane. It has homologues in both mitochondria and chloroplasts. Here we report the crystal structure of the Escherichia coli YidC major periplasmic domain (YidCECP1) at 2.5A resolution. This domain is present in YidC from Gram-negative bacteria and is more than half the size of the full-length protein. The structure reveals that YidCECP1 is made up of a large twisted beta-sandwich protein fold with a C-terminal alpha-helix that packs against one face of the beta-sandwich. Our structure and sequence analysis reveals that the C-terminal alpha-helix and the beta-sheet that it lays against are the most conserved regions of the domain. The region corresponding to the C-terminal alpha-helix was previously shown to be important for the protein insertase function of YidC and is conserved in other YidC-like proteins. The structure reveals that a region of YidC that was previously shown to be involved in binding to SecF maps to one edge of the beta-sandwich. Electrostatic analysis of the molecular surface for this region of YidC reveals a predominantly charged surface and suggests that the SecF-YidC interaction may be electrostatic in nature. Interestingly, YidCECP1 has significant structural similarity to galactose mutarotase from Lactococcus lactis, suggesting that this domain may have another function besides its role in membrane protein assembly.