Symmetry-based self-assembled nanotubes constructed using native protein structures: the key role of flexible linkers

We construct nanotubes using native protein structures and their native associations from structural databases. The construction is based on a shape-guided symmetric self-assembly concept. Our strategy involves fusing judiciously-selected oligomerization domains via peptide linkers. Linkers are inherently flexible, hence their choice is critical: they should position the domains in three-dimensional space in the desired orientation while retaining their own natural conformational tendencies; however, at the same time, retain the construct stability. Here we outline a design scheme which accounts for linker flexibility considerations, and present two examples. The first is HIV-1 capsid protein, which in vitro self-assembles into nanotubes and conical capsids, and its linker exists as a short flexible loop. The second involves novel nanotubes construction based on antimicrobial homodimer Magainin 2, employing linkers of distinct lengths and flexibility levels. Our strategy utilizes the abundance of unique shapes and sizes of proteins and their building blocks which can assemble into a vast number of combinations, and consequently, nanotubes of distinct morphologies and diameters. Computational design and assessment methodologies can help reduce the number of candidates for experimental validation. This is an invited paper for a special issue on protein dynamics, here focusing on flexibility in nanotube design based on protein building blocks.     

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.     

Extended conformational states dominate the Hsp90 chaperone dynamics

The heat shock protein 90 (Hsp90) is a molecular chaperone central to client protein folding and maturation in eukaryotic cells. During its chaperone cycle, Hsp90 undergoes ATPase-coupled large-scale conformational changes between open and closed states, where the N-terminal and middle domains of the protein form a compact dimerized conformation. However, the molecular principles of the switching motion between the open and closed states remain poorly understood. Here we show by integrating atomistic and coarse-grained molecular simulations with small-angle X-ray scattering experiments and NMR spectroscopy data that Hsp90 exhibits rich conformational dynamics modulated by the charged linker, which connects the N-terminal with the middle domain of the protein. We show that the dissociation of these domains is crucial for the conformational flexibility of the open state, with the separation distance controlled by a β-sheet motif next to the linker region. Taken together, our results suggest that the conformational ensemble of Hsp90 comprises highly extended states, which could be functionally crucial for client processing.     

Probing Interdomain Linkers and Protein Supertertiary Structure In Vitro and in Live Cells with Fluorescent Protein Resonance Energy Transfer

Many proteins are composed of independently-folded domains connected by flexible linkers. The primary sequence and length of such linkers can set the effective concentration for the tethered domains, which impacts rates of association and enzyme activity. The length of such linkers can be sensitive to environmental conditions, which raises questions as to how studies in dilute buffer relate to the highly-crowded cellular environment. To examine the role of linkers in domain separation, we measured Fluorescent Protein-Fluorescence Resonance Energy Transfer (FP-FRET) for a series of tandem FPs that varied in the length of their interdomain linkers. We used discrete molecular dynamics to map the underlying conformational distribution, which revealed intramolecular contact states that we confirmed with single molecule FRET. Simulations found that attached FPs increased linker length and slowed conformational dynamics relative to the bare linkers. This makes the CLYs poor sensors of inherent linker properties. However, we also showed that FP-FRET in CLYs was sensitive to solvent quality and macromolecular crowding making them potent environmental sensors. Finally, we targeted the same proteins to the plasma membrane of living mammalian cells to measure FP-FRET in cellulo. The measured FP-FRET when tethered to the plasma membrane was the same as that in dilute buffer. While caveats remain regarding photophysics, this suggests that the supertertiary conformational ensemble of these CLY proteins may not be affected by this specific cellular environment.     

Characterization and prediction of linker sequences of multi-domain proteins by a neural network

In this paper, we describe a neural network analysis of sequences connecting two protein domains (domain linkers). The neural network was trained to distinguish between domain linker sequences and non-linker sequences, using a SCOP-defined domain library. The analysis indicated that a significant difference existed between domain linkers and non-linker regions, including intra-domain loop regions. Moreover, the resulting Hinton diagram showed a position-dependent amino acid preference of the domain linker sequences, and implied their non-random nature. We then applied the neural network to predict domain linkers in multi-domain protein sequences. As the result of a Jack-knife test, 58% of the predicted regions matched actual linker regions (specificity), and 36% of the SCOP-derived domain linkers were predicted (sensitivity). This prediction efficiency is superior to simpler methods derived from secondary structure prediction that assume that long loop regions are putative domain linkers. Altogether, these results suggest that domain linkers possess local characteristics different from those of loop regions.     

Single-molecule observations of neck linker conformational changes in the kinesin motor protein

Kinesin-1 is a dimeric motor protein that moves cargo processively along microtubules. Kinesin motility has been proposed to be driven by the coordinated forward extension of the neck linker (a approximately 12-residue peptide) in one motor domain and the rearward positioning of the neck linker in the partner motor domain. To test this model, we have introduced fluorescent dyes selectively into one subunit of the kinesin dimer and performed 'half-molecule' fluorescence resonance energy transfer to measure conformational changes of the neck linker. We show that when kinesin binds with both heads to the microtubule, the neck linkers in the rear and forward heads extend forward and backward, respectively. During ATP-driven motility, the neck linkers switch between these conformational states. These results support the notion that neck linker movements accompany the 'hand-over-hand' motion of the two motor domains.     

Cellulase Linkers Are Optimized Based on Domain Type and Function: Insights from Sequence Analysis,Biophysical Measurements,and Molecular Simulation

Cellulase enzymes deconstruct cellulose to glucose, and are often comprised of glycosylated linkers connecting glycoside hydrolases (GHs) to carbohydrate-binding modules (CBMs). Although linker modifications can alter cellulase activity, the functional role of linkers beyond domain connectivity remains unknown. Here we investigate cellulase linkers connecting GH Family 6 or 7 catalytic domains to Family 1 or 2 CBMs, from both bacterial and eukaryotic cellulases to identify conserved characteristics potentially related to function. Sequence analysis suggests that the linker lengths between structured domains are optimized based on the GH domain and CBM type, such that linker length may be important for activity. Longer linkers are observed in eukaryotic GH Family 6 cellulases compared to GH Family 7 cellulases. Bacterial GH Family 6 cellulases are found with structured domains in either N to C terminal order, and similar linker lengths suggest there is no effect of domain order on length. O-glycosylation is uniformly distributed across linkers, suggesting that glycans are required along entire linker lengths for proteolysis protection and, as suggested by simulation, for extension. Sequence comparisons show that proline content for bacterial linkers is more than double that observed in eukaryotic linkers, but with fewer putative O-glycan sites, suggesting alternative methods for extension. Conversely, near linker termini where linkers connect to structured domains, O-glycosylation sites are observed less frequently, whereas glycines are more prevalent, suggesting the need for flexibility to achieve proper domain orientations. Putative N-glycosylation sites are quite rare in cellulase linkers, while an N-P motif, which strongly disfavors the attachment of N-glycans, is commonly observed. These results suggest that linkers exhibit features that are likely tailored for optimal function, despite possessing low sequence identity. This study suggests that cellulase linkers may exhibit function in enzyme action, and highlights the need for additional studies to elucidate cellulase linker functions.     

An analysis of protein domain linkers: their classification and role in protein folding

Recent advances in protein engineering have come from creating multi-functional chimeric proteins containing modules from various proteins. These modules are typically joined via an oligopeptide linker, the correct design of which is crucial for the desired function of the chimeric protein. Here we analyse the properties of naturally occurring inter-domain linkers with the aim to design linkers for domain fusion. Two main types of linker were identified; helical and non-helical. Helical linkers are thought to act as rigid spacers separating two domains. Non-helical linkers are rich in prolines, which also leads to structural rigidity and isolation of the linker from the attached domains. This means that both linker types are likely to act as a scaffold to prevent unfavourable interactions between folding domains. Based on these results we have constructed a linker database intended for the rational design of linkers for domain fusion, which can be accessed via the Internet at http://mathbio.nimr.mrc.ac.uk.     

Exploring the Structure-Function Loop Adaptability of a (β/α)(8)-Barrel Enzyme through Loop Swapping and Hinge Variability

Cellulosomes are large extracellular multi-enzyme complexes that exhibit elevated activity on plant cell-wall polysaccharides. In the present study, the relationships between the conformational flexibility and efficacy of cellulosomes, and the inter-modules linkers of their scaffold protein were investigated. For this purpose, the length of the intrinsically disordered Ser/Thr-rich 50-residue linker connecting a Clostridium thermocellum and a Clostridium cellulolyticum cohesin in a hybrid scaffoldin (Scaf4) was changed by sequences ranging from 4 to 128 residues. The composition was also modified and new linkers composed of series of N, S or repeats of the EPPV motif were generated. Two model cellulases (Cel48F and Cel9G) appended with appropriate dockerins were subsequently bound to the engineered scaffoldins. All the resulting minicomplexes displayed the same activity on crystalline cellulose as the complex based on the initial Scaf4, and were found to be 2-fold more active than Cel48F and Cel9G bound to separate cohesins. Small-angle X-ray scattering assays of the engineered scaffoldins confirmed, however, that the size and the conformational flexibility of some of the new inter-cohesins linkers differed significantly from that of the initial 50 residue linker displayed by the parental Scaf4. Our data suggest that the synergy induced by proximity does not require a specific inter-cohesins sequence or distance. The present study reveals that complexation onto the hybrid scaffoldins modifies the type of soluble sugars released from crystalline cellulose by the selected cellulases, compared to the free enzyme system.     

Role of the S3-S4 Linker in Shaker Potassium Channel Activation

Structural models of voltage-gated channels suggest that flexibility of the S3-S4 linker region may be important in allowing the S4 region to undergo large conformational changes in its putative voltage-sensing function. We report here the initial characterization of 18 mutations in the S3-S4 linker of the Shaker channel, including deletions, insertions, charge changes, substitution of prolines, and chimeras replacing the 25-residue Shaker linker with 7- or 9-residue sequences from Shab, Shaw, or Shal. As measured in Xenopus oocytes with a two-microelectrode voltage clamp, each mutant construct yielded robust currents. Changes in the voltage dependence of activation were small, with activation voltage shifts of 13 mV or less. Substitution of linkers from the slowly activating Shab and Shaw channels resulted in a three- to fourfold slowing of activation and deactivation. It is concluded that the S3-S4 linker is unlikely to participate in a large conformational change during channel activation. The linker, which in some channel subfamilies has highly conserved sequences, may however be a determinant of activation kinetics in potassium channels, as previously has been suggested in the case of calcium channels.     

Quantitative understanding of the energy transfer between fluorescent proteins connected via flexible peptide linkers

The fusion of different protein domains via peptide linkers is a powerful, modular approach to obtain proteins with new functions. A detailed understanding of the conformational behavior of peptide linkers is important for applications such as fluorescence resonance energy transfer (FRET)-based sensor proteins and multidomain proteins involved in multivalent interactions. To investigate the conformational behavior of flexible glycine- and serine-containing peptide linkers, we constructed a series of fusion proteins of enhanced cyan and yellow fluorescent proteins (ECFP-linker-EYFP) in which the linker length was systematically varied by incorporating between 1 and 9 GGSGGS repeats. As expected, both steady-state and time-resolved fluorescence measurements showed a decrease in energy transfer with increasing linker length. The amount of energy transfer observed in these fusion proteins can be quantitatively understood by simple models that describe the flexible linker as a worm-like chain with a persistence length of 4.5 A or a Gaussian chain with a characteristic ratio of 2.3. The implications of our results for understanding the properties of FRET-based sensors and other fusion proteins with Gly/Ser linkers are discussed.     

Characteristics and prediction of domain linker sequences in multi-domain proteins

To facilitate swift structural characterizations, structural genomic/proteomic projects need to divide large multi-domain proteins into structural domains and to determine their structures separately. Thus, the assignment of structural domains based solely on sequence information, especially on the physico-chemical properties of the amino acid sequences, could be very helpful for such projects. In this study, we examined the characteristics of domain linker sequences, which are loop sequences connecting two structural domains. To this end, we prepared a set of 101 non-redundant multi-domain protein sequences with known structures, and performed an analysis of the linker sequences. The analysis revealed that the frequencies of five (Pro, Gly, Asp, Asn, Lys) amino acid residues differed significantly between the linker and non-linker loop sequences. Moreover, we observed a similar deviation for the residue pair frequencies between the two types of loop sequences. Finally, we describe an automated method, based on the above analysis, to detect loops that have high probabilities of being domain linkers in a protein sequence.     

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.     

Optimized linker sequences for the expression of monomeric and dimeric bispecific single-chain diabodies.

Bispecific single-chain diabodies (scDb) consist of the variable heavy and light chain domains of two antibodies connected by three linkers. The structure of an scDb in the V(H)-V(L) orientation is V(H)A-linkerA-V(L)B-linkerM-V(H)B-linkerB-V(L)A, with linkers A and B routinely chosen to be 5-6 residues and linker M 15-20 residues. Here, we applied display of scDb on filamentous phage to analyse the composition of optimal linker sequences. The three linkers were randomized in length and sequence using degenerated triplets coding for only six hydrophilic or aliphatic amino acids (Thr, Ser, Asp, Asn, Gly, Ala). Antigen-binding clones were then isolated by one to two rounds of selection on the two different antigens recognized by the bispecific scDb. Using an scDb directed against carcinoembryonic antigen (CEA) and beta-galactosidase (Gal), we found that monomeric scDb had a preferred length of 15 or more amino acid residues for the middle linker M and of 3-6 residues for the linkers A and B. No obvious bias towards a preferred linker sequence was observed. Reduction of the middle linker below 13 residues led to the formation of dimeric scDb, which most likely results from interchain pairing between all the V(H) and V(L) domains. Dimeric scDb were also formed by fragments possessing a long linker M and linkers A and B of 0 or 1 residue. We assume that these dimeric scDb are formed by intrachain pairing of the central variable domains and interchain pairing of the flanking variable domains. Thus, the latter molecules represent a novel format of bispecific and tetravalent molecules. The described strategy allows for the isolation of both optimized and minimal linker sequences for the assembly of monomeric or dimeric single-chain diabodies.     

Structure and Dynamics of Ribosomal Protein L12: An Ensemble Model Based on SAXS and NMR Relaxation

Ribosomal protein L12 is a two-domain protein that forms dimers mediated by its N-terminal domains. A 20-residue linker separates the N- and C-terminal domains. This linker results in a three-lobe topology with significant flexibility, known to be critical for efficient translation. Here we present an ensemble model of spatial distributions and correlation times for the domain reorientations of L12 that reconciles experimental data from small-angle x-ray scattering and nuclear magnetic resonance. We generated an ensemble of L12 conformations in which the structure of each domain is fixed but the domain orientations are variable. The ensemble reproduces the small-angle x-ray scattering data and the optimized correlation times of its reorientational eigenmodes fit the ¹⁵N relaxation data. The ensemble model reveals intrinsic conformational properties of L12 that help explain its function on the ribosome. The two C-terminal domains sample a large volume and extend further away from the ribosome anchor than expected for a random-chain linker, indicating that the flexible linker has residual order. Furthermore, the distances between each C-terminal domain and the anchor are anticorrelated, indicating that one of them is more retracted on average. We speculate that these properties promote the function of L12 to recruit translation factors and control their activity on the ribosome.     

Linkers in the structural biology of protein–protein interactions

Linkers or spacers are short amino acid sequences created in nature to separate multiple domains in a single protein. Most of them are rigid and function to prohibit unwanted interactions between the discrete domains. However, Gly‐rich linkers are flexible, connecting various domains in a single protein without interfering with the function of each domain. The advent of recombinant DNA technology made it possible to fuse two interacting partners with the introduction of artificial linkers. Often, independent proteins may not exist as stable or structured proteins until they interact with their binding partner, following which they gain stability and the essential structural elements. Gly‐rich linkers have been proven useful for these types of unstable interactions, particularly where the interaction is weak and transient, by creating a covalent link between the proteins to form a stable protein–protein complex. Gly‐rich linkers are also employed to form stable covalently linked dimers, and to connect two independent domains that create a ligand‐binding site or recognition sequence. The lengths of linkers vary from 2 to 31 amino acids, optimized for each condition so that the linker does not impose any constraints on the conformation or interactions of the linked partners. Various structures of covalently linked protein complexes have been described using X‐ray crystallography, nuclear magnetic resonance and cryo‐electron microscopy techniques. In this review, we evaluate several structural studies where linkers have been used to improve protein quality, to produce stable protein–protein complexes, and to obtain protein dimers.     

Ligand-binding properties and conformational dynamics of autolysin repeat domains in staphylococcal cell wall recognition

The bifunctional major autolysin Atl plays a key role in staphylococcal cell separation. Processing of Atl yields catalytically active amidase (AM) and glucosaminidase (GL) domains that are each fused to repeating units. The two repeats of AM (R1 and R2) target the enzyme to the septum, where it cleaves murein between dividing cells. We have determined the crystal structure of R2, which reveals that each repeat folds into two half-open β-barrel subunits. We further demonstrate that lipoteichoic acid serves as a receptor for the repeats and that this interaction depends on conserved surfaces in each subunit. Small-angle X-ray scattering of the mature amidase reveals the presence of flexible linkers separating the AM, R1, and R2 units. Different levels of flexibility for each linker provide mechanistic insights into the conformational dynamics of the full-length protein and the roles of its components in cell wall association and catalysis. Our analysis supports a model in which the repeats direct the catalytic AM domain to the septum, where it can optimally perform the final step of cell division.     

Crystal structure and NMR studies of the apo SH2 domains of ZAP-70: two bikes rather than a tandem

The protein kinase ZAP-70 is involved in T-cell activation, and interacts with tyrosine-phosphorylated peptide sequences known as immunoreceptor tyrosine activation motifs (ITAMs), which are present in three of the subunits of the T-cell receptor. We have studied the tandem SH2 (tSH2) domains of ZAP-70, by both X-ray and NMR. Here, we present the crystal structure of the apoprotein, i.e., the tSH2 domain in the absence of ITAM. Comparison with the previously reported complex structure reveals that binding to the ITAM peptide induces surprisingly large movements between the two SH2 domains and within the actual binding sites. The conformation of the ITAM-free protein is partly governed by a hydrophobic cluster between the linker region and the C-terminal SH2 domain. Our data suggest that the two SH2 domains are able to undergo large interdomain movements. The proposed relative flexibility of the SH2 domains is further supported by the finding that no NMR signals could be detected for the two helices connecting the SH2 domains; these are likely to be broadened beyond detection due to conformational exchange. It is likely that this conformational reorientation induced by ITAM binding is the main signaling event activating the kinase domain in ZAP-70. Another NMR observation was that the N-terminal SH2 domain could bind tetrapeptides derived from the ITAM sequence, apparently without the need to interact with the C-terminal domain. In contrast, the C-terminal domain has little affinity for tetrapeptides. The opposite situation is true for binding to plain phosphotyrosine, where the C-terminal domain has a higher affinity. Distinct features in the crystal structure, showing the interdependence of both domains, explain these binding data.     

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.