Prediction of protein interdomain linker regions by a hidden Markov model

MOTIVATION: Our aim was to predict protein interdomain linker regions using sequence alone, without requiring known homology. Identifying linker regions will delineate domain boundaries, and can be used to computationally dissect proteins into domains prior to clustering them into families. We developed a hidden Markov model of linker/non-linker sequence regions using a linker index derived from amino acid propensity. We employed an efficient Bayesian estimation of the model using Markov Chain Monte Carlo, Gibbs sampling in particular, to simulate parameters from the posteriors. Our model recognizes sequence data to be continuous rather than categorical, and generates a probabilistic output. RESULTS: We applied our method to a dataset of protein sequences in which domains and interdomain linkers had been delineated using the Pfam-A database. The prediction results are superior to a simpler method that also uses linker index.     

Armadillo: domain boundary prediction by amino acid composition

The identification and annotation of protein domains provides a critical step in the accurate determination of molecular function. Both computational and experimental methods of protein structure determination may be deterred by large multi-domain proteins or flexible linker regions. Knowledge of domains and their boundaries may reduce the experimental cost of protein structure determination by allowing researchers to work on a set of smaller and possibly more successful alternatives. Current domain prediction methods often rely on sequence similarity to conserved domains and as such are poorly suited to detect domain structure in poorly conserved or orphan proteins. We present here a simple computational method to identify protein domain linkers and their boundaries from sequence information alone. Our domain predictor, Armadillo (http://armadillo.blueprint.org), uses any amino acid index to convert a protein sequence to a smoothed numeric profile from which domains and domain boundaries may be predicted. We derived an amino acid index called the domain linker propensity index (DLI) from the amino acid composition of domain linkers using a non-redundant structure dataset. The index indicates that Pro and Gly show a propensity for linker residues while small hydrophobic residues do not. Armadillo predicts domain linker boundaries from Z-score distributions and obtains 35% sensitivity with DLI in a two-domain, single-linker dataset (within +/-20 residues from linker). The combination of DLI and an entropy-based amino acid index increases the overall Armadillo sensitivity to 56% for two domain proteins. Moreover, Armadillo achieves 37% sensitivity for multi-domain proteins, surpassing most other prediction methods. Armadillo provides a simple, but effective method by which prediction of domain boundaries can be obtained with reasonable sensitivity. Armadillo should prove to be a valuable tool for rapidly delineating protein domains in poorly conserved proteins or those with no sequence neighbors. As a first-line predictor, domain meta-predictors could yield improved results with Armadillo predictions.     

An improved profile-level domain linker propensity index for protein domain boundary prediction

Protein domain boundary prediction is critical for understanding protein structure and function. In this study, we present a novel method, an order profile domain linker propensity index (OPI), which uses the evolutionary information extracted from the protein sequence frequency profiles calculated from the multiple sequence alignments. A protein sequence is first converted into smooth and normalized numeric order profiles by OPI, from which the domain linkers can be predicted. By discriminating the different frequencies of the amino acids in the protein sequence frequency profiles, OPI clearly shows better performance than our previous method, a binary profile domain linker propensity index (PDLI). We tested our new method on two different datasets, SCOP-1 dataset and SCOP-2 dataset, and we were able to achieve a precision of 0.82 and 0.91 respectively. OPI also outperforms other residue-level, profile-level indexes as well as other state-of-the-art methods.     

Annealing control primer system for improving specificity of PCR amplification

A novel primer designed to improve the specificity of PCR amplification, called the annealing control primer (ACP), comprises a tripartite structure with a polydeoxyinosine [poly(dI)] linker between the 3' end target core sequence and the 5' end nontarget universal sequence. We show that this ACP linker prevents annealing of the 5' end nontarget sequence to the template and facilitates primer hybridization at the 3' end to the target sequence at specific temperatures, resulting in a dramatic improvement of annealing specificity. The effect of this linker is demonstrated by the incorporation of ACP sequences as primers during the amplification of target nucleotide sequence and as hybridization probes in the genotyping of single nucleotide polymorphisms. This is the first report to show that a poly(dI) linker between two different sequences of ACP forms a bubble-like structure and disrupts or destabilizes DNA duplex formation at certain annealing temperatures.     

Probing the primary structural determinants of streptokinase inter-domain linkers by site-specific substitution and deletion mutagenesis

The bacterial protein streptokinase (SK) contains three independently folded domains (α, β and γ), interconnected by two flexible linkers with noticeable sequence homology. To investigate their primary structure requirements, the linkers were swapped amongst themselves i.e. linker 1 (between α and β domains) was swapped with linker 2 (between β and γ domains) and vice versa. The resultant construct exhibited very low activity essentially due to an enhanced proteolytic susceptibility. However, a SK mutant with two linker 1 sequences, which was proteolytically as stable as WT-rSK retained about 10% of the plasminogen activator activity of rSK When the native sequence of each linker was substituted with 9 consecutive glycine sequences, in case of the linker 1 substitution mutant substantial activity was seen to survive, whereas the linker 2 mutant lost nearly all its activity. The optimal length of linkers was then studied through deletion mutagenesis experiments, which showed that deletion beyond three residues in either of the linkers resulted in virtually complete loss of activator activity. The effect of length of the linkers was then also examined by insertion of extraneous pentapeptide sequences having a propensity for adopting either an extended conformation or a relatively rigid conformation. The insertion of poly-Pro sequences into native linker 2 sequence caused up to 10-fold reduction in activity, whereas its effect in linker 1 was relatively minor. Interestingly, most of the linker mutants could form stable 1:1 complexes with human plasminogen. Taken together, these observations suggest that (i) the functioning of the inter-domain linkers of SK requires a critical minimal length, (ii) linker 1 is relatively more tolerant to insertions and sequence alterations, and appears to function primarily as a covalent connector between the α and β domains, and (iii) the native linker 2 sequence is virtually indispensable for the activity of SK probably because of structural and/or flexibility requirements in SK action during catalysis.     

LINKER: a program to generate linker sequences for fusion proteins

The construction of functional fusion proteins often requires a linker sequence that adopts an extended conformation to allow for maximal flexibility. Linker sequences are generally selected based on intuition. Without a reliable selection criterion, the design of such linkers is often difficult, particularly in situations where longer linker sequences are required. Here we describe a program called LINKER which can automatically generate a set of linker sequences that are known to adopt extended conformations as determined by X-ray crystallography and NMR. The only required input to the program is the desired linker sequence length. The program is specifically designed to assist in fusion protein construction. A number of optional input parameters have been incorporated so that users are able to enhance sequence selection based on specific applications. The program output simply contains a set of sequences with a specified length. This program should be a useful tool in both the biotechnology industry and biomedical research. It can be accessed through the Web page http://www.fccc. edu/research/labs/feng/linker.html.     

Activity of a two-domain antifreeze protein is not dependent on linker sequence

The reported NMR structure of RD3, a naturally occurring two-domain antifreeze protein, suggests that the two nearly identical domains are oriented to allow simultaneous binding of their active regions to the ice surface. It is implied that the nine residues linking the two domains play a role in this alignment, but this has not been established. We have designed and expressed a modified form of RD3 that replaces the nine-residue linker with a generic sequence of one serine and eight glycine residues to test the importance of the linker amino acid sequence. The modified linker is shown to have significantly different characteristics compared to the original linker. Heteronuclear nuclear Overhauser effect experiments show that the new linker residues have more mobility than the linker residues in the native protein. Further, NMR data show that the folding of the C-terminal domain is somewhat perturbed by the altered linker. Finally, distributions of residual dipolar couplings indicate that the two domains tumble and move independently of each other. Nevertheless, the thermal hysteresis activity of the modified protein is indistinguishable from that of native RD3, proving that increased activity of the two-domain antifreeze protein is not dependent on structure of the linker.     

Primary structure of a linker subunit of the tube worm 3000-kDa hemoglobin

The deep-sea tube worm Lamellibrachia contains two giant extracellular hemoglobins, a 3000-kDa hemoglobin and a 440-kDa hemoglobin. The former consists of four heme-containing chains (AI-AIV) and two linker chains (AV and AVI) for the assembly of the heme-containing chains. The 440-kDa hemoglobin consists of only four heme-containing chains (Suzuki, T., Takagi, T., and Ohta, S. (1988) Biochem. J. 255, 541-545). The complete amino acid sequence of a linker subunit (chain AV) has been determined by automated Edman sequencing of the peptides derived by digestions with lysyl endopeptidase and endoproteinase Asp-N. The chain is composed of 224 amino acid residues, and the molecular mass for the protein moiety was calculated to be 24,894 Da. An Asn-X-Thr sequence which is possible as a glycosylation site was suggested at positions 108-110. A computer-assisted homology search showed that the sequence shows no notable homology with any other globins and proteins. However a careful alignment of the linker sequence with a heme-containing chain sequence suggested that there is a slight, but significant homology between the two sequences. The alignment also suggested that the linker resulted from gene duplication of a heme-containing chain with a three exon-two intron structure, and that the first exon of domain 1 and the last exon of domain 2 had been lost during evolution. In our alignment, domain 1 has the heme-binding proximal histidine, but domain 2 does not. This is the first linker subunit to be sequenced completely.     

Endoglucanase A from Cellulomonas fimi in which the hinge sequence of human IgA1 is substituted for the linker connecting its two domains is hydrolyzed by IgA proteases from Neisseria gonorrhoeae

The hinge in IgA1 and the linker in endoglucanase A (CenA) are quite similar. The IgA1 hinge is 18 amino acids long and contains only proline, threonine and serine. The linker in CenA is 27 amino acids long and contains only proline, threonine and a single serine. IgA proteases from Neisseria gonorrhoeae cleave Pro-Ser and Pro-Thr bonds within the IgA1 hinge sequence, but they do not attack CenA. When the linker sequence of CenA is replaced with the hinge sequence of IgA1, the hybrid polypeptide is susceptible to the N. gonorrhoeae proteases. It is cleaved within the hinge sequence at the same sites as IgA1.     

Enzymatic activity of the Arabidopsis sulfurtransferase resides in the C-terminal domain but is boosted by the N-terminal domain and the linker peptide in the full-length enzyme

Sulfurtransferases/rhodaneses are a group of enzymes widely distributed in plants, animals, and bacteria that catalyze the transfer of sulfur from a donor molecule to a thiophilic acceptor substrate. Sulfurtransferases (STs) consist of two globular domains of nearly identical size and conformation connected by a short linker sequence. In plant STs this linker sequence is exceptionally longer than in sequences from other species. The Arabidopsis ST1 protein (AJ131404) contains five cysteine residues: one residue is universally conserved in all STs and considered to be catalytically essential; a second one, closely located in the primary sequence, is conserved only in sequences from eukaryotic species. Of the remaining three cysteine residues two are conserved in the so far known plant STs and one is unique to the Arabidopsis ST1. The aim of our study was to investigate the role of the two-domain structure, of the unique plant linker sequence and of each cysteine residue. The N- and C-terminal domains of the Arabidopsis ST1, the full-length protein with a shortened linker sequence and several point-mutated proteins were overexpressed in E. coli, purified and used for enzyme activity measurements. The C-terminal domain itself displayed ST activity which could be increased by adding the separately prepared N-terminal domain. The activity of an ST1 derivative with a shortened linker sequence was reduced by more than 60% of the wild-type activity, probably because of a drastically reduced protein stability. The replacement of each cysteine residue resulted in mutant forms which differed significantly in their stability, in the specific ST activities, and in their kinetic parameters which were determined for 3-mercaptopyruvate as well as thiosulfate as sulfur substrates: mutation of the putative active site cysteine (C332) essentially abolished activity; for C339 a crucial role at least for the turnover of thiosulfate could be identified.     

Modulation of beta-amyloid aggregation by engineering the sequence connecting beta-strand forming domains

Aggregation of beta-amyloid (Aβ) into oligomers and fibrils is associated with the pathology of Alzheimer's disease. The major structural characteristics of Aβ fibrils include the presence of β sheet-loop-β sheet conformations. Several lines of study suggested a potentially important role of the Aβ loop forming sequence (referred to as the Aβ linker region) in Aβ aggregation. Effects of mutations in several charged residues within the Aβ linker region on aggregation have been extensively studied. However, little is known about oligomerization effects of sequence variation in other residues within the Aβ linker region. Moreover, modulation effects of the Aβ linker mutants on Aβ aggregation have yet to be characterized. Here, we created and characterized Aβ linker variants containing sequences preferentially found in specific β turn conformations. Our results indicate that a propensity to form oligomers may be changed by local sequence variation in the Aβ linker region without mutating the charged residues. Strikingly, one Aβ linker variant rapidly formed protofibrillar oligomers, which did not convert to fibrillar aggregates in contrast to Aβ aggregating to fibrils under similar incubation conditions. Moreover, our results suggest that molecular forces critical in oligomerization and fibrillization may differ at least for those involved in the linker region. When co-incubated with Aβ, some Aβ linker variants were found to induce accumulation of Aβ oligomers. Our results suggest that engineering of the Aβ linker region as described in this paper may represent a novel approach to control Aβ oligomerization and create Aβ oligomerization modulators.     

A flexible C‐terminal linker is required for proper FtsZ assembly in vitro and cytokinetic ring formation in vivo

Assembly of the cytoskeletal protein FtsZ into a ring‐like structure is required for bacterial cell division. Structurally, FtsZ consists of four domains: the globular N‐terminal core, a flexible linker, 8–9 conserved residues implicated in interactions with modulatory proteins, and a highly variable set of 4–10 residues at its very C terminus. Largely ignored and distinguished by lack of primary sequence conservation, the linker is presumed to be intrinsically disordered. Here we employ genetics, biochemistry and cytology to dissect the role of the linker in FtsZ function. Data from chimeric FtsZs substituting the native linker with sequences from unrelated FtsZs as well as a helical sequence from human beta‐catenin indicate that while variations in the primary sequence are well tolerated, an intrinsically disordered linker is essential for Bacillus subtilis FtsZ assembly. Linker lengths ranging from 25 to 100 residues supported FtsZ assembly, but replacing the B. subtilis FtsZ linker with a 249‐residue linker from Agrobacterium tumefaciens FtsZ interfered with cell division. Overall, our results support a model in which the linker acts as a flexible tether allowing FtsZ to associate with the membrane through a conserved C‐terminal domain while simultaneously interacting with itself and modulatory proteins in the cytoplasm.     

Configuration of the catalytic GIY-YIG domain of intron endonuclease I-TevI: coincidence of computational and molecular findings

I-TevI is a member of the GIY-YIG family of homing endonucleases. It is folded into two structural and functional domains, an N-terminal catalytic domain and a C-terminal DNA-binding domain, separated by a flexible linker. In this study we have used genetic analyses, computational sequence analysis andNMR spectroscopy to define the configuration of theN-terminal domain and its relationship to the flexible linker. The catalytic domain is an alpha/beta structure contained within the first 92 amino acids of the 245-amino acid protein followed by an unstructured linker. Remarkably, this structured domain corresponds precisely to the GIY-YIG module defined by sequence comparisons of 57 proteins including more than 30 newly reported members of the family. Although much of the unstructured linker is not essential for activity, residues 93-116 are required, raising the possibility that this region may adopt an alternate conformation upon DNA binding. Two invariant residues of the GIY-YIG module, Arg27 and Glu75, located in alpha-helices, have properties of catalytic residues. Furthermore, the GIY-YIG sequence elements for which the module is named form part of a three-stranded antiparallel beta-sheet that is important for I-TevI structure and function.     

A new assay for functional screening of BRCA2 linker region mutations identifies variants that alter chemoresistance to cisplatin

Variants of unknown significance (VUS) complicate the assignment of risk to new DNA sequence variants found in at-risk populations. This study focused on the poorly studied linker region of the cancer-associated BRCA2 protein encoded by exons twelve through fourteen of BRCA2. To develop a new method to characterize VUS in this region of BRCA2, we first chose to study 4 reported VUS occurring on evolutionarily conserved residues within the linker region. To determine if these VUS represent neutral changes or if they impact the function of the BRCA2 protein, we stably transfected expression plasmids encoding wild-type or each mutant peptide into T47D breast cancer cells, which are wild-type for BRCA2. Four mutant peptide expressing cell lines and a wild-type linker region expressing cell line next were studied by challenging transfected cell lines with the DNA crosslinking compound cisplatin (10 μM) for 5 days. Expression of the wild-type linker region and certain mutant linker peptides (N2452D and I2285V) decreased apoptosis (as demonstrated by cell death detection assay) in transfected cell lines, indicating that the linker region peptide directly or indirectly affects the DNA damage repair pathway. By determining the cell survival and assaying the apoptotic index of treated cell lines, one could potentially use this screen to determine that a particular VUS has a functional impact on BRCA2 function, and hence is of functional significance. We conclude that this method is useful for screening the effect of linker region VUS on BRCA2 function, and to identify mutations for further testing. We also conclude that mutations in the linker region may have heretofore unappreciated roles in BRCA2 function.     

New type of linker useful for cloning experiments

A new class of linker oligodeoxynucleotide sequences is described which allows the original sequence of a foreign DNA to be restored after cloning. In the standard linker approach the terminal base pairs of the linker-sequences are irreversibly attached to the cloned DNA, thus altering the genetic information. The use of the new type of linker is demonstrated by the transformation of the unique Pvu II site in pBR322 into Bam HI site using the ISO ('in-site-out') linker d(GATCCGGATC). Possible further applications of these linker sequences are described.     

The carboxyl-terminal linker is important for chemoreceptor function

Sensory adaptation in bacterial chemotaxis is mediated by chemoreceptor methylation and demethylation. In Escherichia coli, methyltransferase CheR and methylesterase CheB bind both substrate sites and a carboxyl-terminal pentapeptide sequence carried by certain receptors. Pentapeptide binding enhances enzyme action, an enhancement required for effective adaptation and chemotaxis. Pentapeptides are linked to the conserved body of chemoreceptors through a notably variable sequence of 30-35 residues. We created nested deletions from the distal end of this linker in chemoreceptor Tar. Chemotaxis was eliminated by deletion of 20-40 residues and reduced by shorter deletions. This did not reflect generalized disruption, because all but the most extremely truncated receptors activated kinase, were substrates for adaptational modification and performed transmembrane signalling. In contrast, linker truncations reduced rates of adaptational modification in parallel with chemotaxis. We concluded the linker is important for chemotaxis because of its role in adaptational modification. Effects of linker truncations on CheR binding to receptor-borne pentapeptide implied linker (i) makes pentapeptide available to modification enzymes by separation from the helical receptor body, and (ii) is a flexible arm allowing dual binding of enzyme to pentapeptide and modification site. The data suggest linker and the helix from which it emerges are structurally dynamic.     

Single strand targeted triplex formation: targeting purine-pyrimidine mixed sequences using abasic linkers.

Foldback triplex-forming oligonucleotides (FTFOs) that contain an abasic linker, [2-(4-aminobutyr-1-yl)-1,3-propanediol] (APD linker), in the Hoogsteen domain against pyrimidine bases of a C:G and a T:A base pair were studied for their relative stability and sequence specificity of triplex formation. In general, the APD linker has less destabilizing effect against a C:G base pair than a T:A base pair. Incorporation of three APD linker moieties resulted in decreased binding to the target, which was comparable to results observed with three imperfectly matched natural base triplets. The APD linker incorporation did not result in the loss of sequence specificity of FTFOs, unlike in the case of normal triplex-forming oligonucleotides (TFOs). The introduction of a positively charged abasic linker, however, resulted in decreased stability of the triplex, because of loss of hydrogen bonding and stacking interactions in the major groove. The results of a molecular modeling study show that APD linker can be readily incorporated without any change in the conformation of the natural sugar-phosphate backbone conserving overall triple helix geometry. Further, the modeling study suggests a hydrogen bond formation between the amino group of linker and N4 of cytosine mediated by a solvent molecule (water) in the floor of the base triplet in addition to a contribution from the positive charge on the APD linker amino group. Either a direct or water-mediated hydrogen bond between the amino group of the APD linker and the O4 of thymine is unlikely when the linker is placed against a T:A base pair.     

Effects of length and position of an extended linker on sequence-selective DNA recognition of zinc finger peptides

Engineered zinc finger proteins revealed that a linker sequence connecting zinc finger units has a significant effect on the DNA binding property of the protein. The recognition for a noncontiguous DNA target beyond the current recognition code of zinc finger proteins has never been determined because of the limitation of a zinc finger framework. DNA recognition of zinc finger proteins is limited only to a contiguous subset of three base pairs. We propose the recognition for a noncontiguous DNA target by inserting amino acids into the canonical linker between zinc finger units. The sequence selectivity of the new zinc finger peptides was evaluated by gel mobility shift assays. DNase I footprinting analyses clearly showed different DNA binding of various linker-extended zinc finger peptides. The application of a SPR measurement also revealed a DNA sequence selectivity of peptides. Insertion of three amino acids is enough for recognition of a noncontiguous DNA target with sequence selectivity. An extended linker will be useful for expansion of the recognition code of zinc finger proteins and for development of a new role for linker sequences in DNA binding of zinc finger proteins.