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.     

A study on the effects of linker flexibility on acid phosphatase PhoC-GFP fusion protein using a novel linker library

Fusion strategy has been widely used to construct artificial multifunction proteins. The flexibility or rigidity of linkers between two fused partners is an important parameter that affects the function of fusion proteins. By combining the flexible unit GGGGS (F) and rigid unit EAAAK (R), ten linkers consisting of five elementary units that cover the fully rigid RRRRR linker to the fully flexible FFFFF linker were used to construct acid phosphatase-green fluorescence protein fusion protein (PhoC-GFP). By varying the linker flexibility in PhoC-GFPs, the relative specific activity of phosphotransferase and phosphatase varied from ∼19.0% to 100% and ∼9.35% to 100%, respectively. There exists an optimal linker capable of achieving the highest phosphotransferase/phosphatase activity and GFP fluorescence intensity. We found that the highest activities were achieved neither with the rigid RRRRR linker nor with the flexible FFFFF linker, but with the FFFRR linker. Linker flexibility could adjust the activity ratio between phosphotransferase and phosphatase and varied between ∼30% to 100%. PhoC-GFP with FRRRR linker achieved the highest relative specific phosphotransferase activity/relative specific phosphatase activity (T/P) value. Our results show that applying a linker library with controllable flexibility to the fusion proteins will be an efficient way to adjust the function of fusion enzymes.     

Nanostructure design using protein building blocks enhanced by conformationally constrained synthetic residues

Increasing efforts are being invested in the construction of nanostructures with desired shapes and physical and chemical properties. Our strategy involves nanostructure design using naturally occurring protein building blocks. Inspection of the protein structural database (PDB) reveals the richness of the conformations, shapes, and chemistries of proteins and their building blocks. To increase the population of the native fold in the selected building block, we mutate natural residues by engineered, constrained residues that restrict the conformational freedom at the targeted site and have favorable interactions, geometry, and size. Here, as a model system, we construct nanotubes using building blocks from left-handed beta-helices which are commonly occurring repeat protein architectures. We pick two-turn beta-helical segments, duplicate and stack them, and using all-atom molecular dynamics simulations (MD) with explicit solvent probe the structural stability of these nanotubular structures as indicated by their capacity to retain the initial organization and their conformational dynamics. Comparison of the results for the wild-type and mutated sequences shows that the introduction of the conformationally restricted 1-aminocyclopropanecarboxylic acid (Ac3c) residue in loop regions greatly enhances the stability of beta-helix nanotubes. The Ac3c geometrical confinement effect is sequence-specific and position-specific. The achievement of high stability of nanotubular structures originates not only from the reduction of mobility at the mutation site induced by Ac3c but also from stabilizing association forces between building blocks such as hydrogen bonds and hydrophobic contacts. For the selected synthetic residue, similar size, hydrophobicity, and backbone conformational tendencies are desirable as in the Ac3c.     

Dynamic allostery: linkers are not merely flexible

Most proteins consist of multiple domains. How do linkers efficiently transfer information between sites that are on different domains to activate the protein? Mere flexibility only implies that the conformations would be sampled. For fast timescales between triggering events and cellular response, which often involves large conformational change, flexibility on its own may not constitute a good solution. We posit that successive conformational states along major allosteric propagation pathways are pre-encoded in linker sequences where each state is encoded by the previous one. The barriers between these states that are hierarchically populated are lower, achieving faster timescales even for large conformational changes. We further propose that evolution has optimized the linker sequences and lengths for efficiency, which explains why mutations in linkers may affect protein function and review the literature in this light.     

Conformations of variably linked chimeric proteins evaluated by synchrotron X-ray small-angle scattering

We constructed chimeric proteins that consist of two green fluorescent protein variants, EBFP and EGFP, connected by flexible linkers, (GGGGS)n (n = 3 approximately 4), and helical linkers, (EAAAK)n (n = 2 approximately 5). The conformations of the chimeric proteins with the various linkers were evaluated using small-angle X-ray scattering (SAXS). The SAXS experiments showed that introducing the short helical linkers (n = 2 approximately 3) causes multimerization, while the longer linkers (n = 4 approximately 5) solvate monomeric chimeric proteins. With the moderate-length linkers (n = 4), the observed radius of gyration (R(g)) and maximum dimension (D(max)) were 38.8 A and 120 A with the flexible linker, and 40.2 A and 130 A with the helical linker, respectively. The chimeric protein with the helical linker assumed a more elongated conformation as compared to that with the flexible linker. When the length of the helical linker increased (n = 5), R(g) and D(max) increased to 43.2 A and 140 A, respectively. These results suggest that the longer helix effectively separates the two domains of the chimeric protein. Considering the connectivity of the backbone peptide of the protein, the helical linker seems to connect the two domains diagonally. Surprisingly, the chimeric proteins with the flexible linker exhibited an elongated conformation, rather than the most compact side-by-side conformation expected from the fluorescence resonance energy transfer (FRET) analysis. Furthermore, the SAXS analyses suggest that destabilization of the short helical linker causes multimerization of the chimeric proteins. Information about the global conformation of the chimeric protein is thus be necessary for optimization of the linker design.     

Designing a nanotube using naturally occurring protein building blocks

Here our goal is to carry out nanotube design using naturally occurring protein building blocks. Inspection of the protein structural database reveals the richness of the conformations of proteins, their parts, and their chemistry. Given target functional protein nanotube geometry, our strategy involves scanning a library of candidate building blocks, combinatorially assembling them into the shape and testing its stability. Since self-assembly takes place on time scales not affordable for computations, here we propose a strategy for the very first step in protein nanotube design: we map the candidate building blocks onto a planar sheet and wrap the sheet around a cylinder with the target dimensions. We provide examples of three nanotubes, two peptide and one protein, in atomistic model detail for which there are experimental data. The nanotube models can be used to verify a nanostructure observed by low-resolution experiments, and to study the mechanism of tube formation.     

NMR studies on domain diffusion and alignment in modular GB1 repeats

Modular proteins contain individual domains that are often connected by flexible, unstructured linkers. Using a model system based on the GB1 domain, we constructed tandem repeat proteins and investigated the rotational diffusion and long-range angular ordering behavior of individual domains by measuring NMR relaxation parameters and residual dipolar couplings. Although they display almost identical protein-solvent interfaces, each domain exhibits distinct rotational diffusion and alignment properties. The diffusion tensor anisotropy of the N-terminal domain (NTD) is D_‖/D_⊥ = 1.5-1.6, similar to that of single-GB1 domains (D_‖/D_⊥ = 1.6-1.7), whereas the value for the C-terminal domain (CTD) is D_‖/D_⊥ = 2.0-2.2. In addition, the two domains have different rotational correlation times. These effects are observed for linkers of three to 24 residues, irrespective of linker length. The NTD and CTD also differ in their degree of magnetic alignment, even with a flexible linker of 18 residues, exhibiting D_a values of 7.7 Hz and 9.7 Hz, respectively. Our results suggest that diffusion differences and long-range influences may persist in modular protein systems, even for systems that have highly flexible linkers and exhibit no domain-domain or domain-linker interactions.     

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.     

Stability of tubular structures based on beta-helical proteins: self-assembled versus polymerized nanoconstructs and wild-type versus mutated sequences

In this work we used atomistic molecular dynamics simulations to examine different aspects of tubular nanostructures constructed using protein building blocks with a beta-helical conformation. Initially, we considered two different natural protein building blocks, which were extracted from the protein data base, to compare the relative stabilities of the nanotubes obtained made of self-assembled and covalently linked repeats. Results show nanotubes constructed by linking building blocks through covalent bonds are very stable suggesting that the basic principles of polymer physics are valid when the repeating units are made of large fragments of proteins. In contrast, the stability of self-assembled nanostructures strongly depends on the attractive nonbonding interactions associated to building blocks aligned in a complementary manner. On the other hand, we investigated the ability of a conformationally constrained synthetic amino acid to enhance the stability of both self-assembled and polymerized nanotubes when it is used to substitute natural residues. Specifically, we considered 1-aminocyclopentane-1-caboxylic acid, which involves strong stereochemical constraints produced by the cyclopentane side chain. We found that the incorporation of this amino acid within the more flexible regions of the beta-helical building blocks is an excellent strategy to enhance the stability of the nanotubes. Thus, when a single mutation is performed in the loop region of the beta-helix, the bend architecture of the whole loop is stabilized since the conformational mobility is reduced not only at the mutated position but also at the adjacent positions.     

Glycosylation effects on the structure and dynamics of a full-length Cel7A cellulase

Fungi cellulases are used to degrade cellulose-containing biomass for bioethanol production. Industrial cellulases such as Cel7A from Trichoderma reesei (TrCel7A) are critical in this process. Thus, the understanding of structure and dynamics is crucial for engineering variants with improved cellulolytic activity. This cellulase consists of two domains connected by a flexible and highly glycosylated linker. However, the linker flexibility has hindered the determination of Cel7A complete structure. Herein, based on atomic and sparse data, we applied integrative modelling to build a model of the complete enzyme structure. Next, through simulations, we studied the glycosylation effects on the structure and dynamics of a solubilized TrCel7A. Essential dynamics analysis showed that O-glycosylation in the linker led to the stabilization of protein overall dynamics. O-linked glycans seem to restrict protein dihedral angles distribution in this region, selecting more elongated conformations. Besides the reduced flexibility, functional interdomain motions occurred in a more concerted way in the glycosylated system. In contrast, in the absence of glycosylation, we observed vast conformational plasticity with the functional domains frequently collapsing. We report here evidence that targeting Cel7A linker flexibility by point mutations including modification of glycosylation sites could be a promising design strategy to improve cellulase activity.     

Domain cooperativity in multidomain proteins: what can we learn from molecular alignment in anisotropic media?

Many proteins have modular design with multiple globular domains connected via flexible linkers. As a simple model of such system, we study a tandem construct consisting of two identical SH3 domains and a variable-length Gly/Ser linker. When the linker is short, this construct represents a dumbbell-shaped molecule with limited amount of domain-domain mobility. Due to its elongated shape, this molecule efficiently aligns in steric alignment media. As the length of the linker increases, the two domains become effectively uncoupled and begin to behave as independent entities. Consequently, their degree of alignment drops, approaching that found in the (near-spherical) isolated SH3 domains. To model the dependence of alignment parameters on the length of the interdomain linker, we have generated in silico a series of conformational ensembles representing SH3 tandems with different linker length. These ensembles were subsequently used as input for alignment prediction software PALES. The predicted alignment tensors were compared with the results of experimental measurements using a series of tandem-SH3 samples in PEG/hexanol alignment media. This comparison broadly confirmed the expected trends. At the same time, it has been found that the isolated SH3 domain aligns much stronger than expected. This finding can be attributed to complex morphology of the PEG/hexanol media and/or to weak site-specific interactions between the protein and the media. In the latter case, there are strong indications that electrostatic interactions may play a role. The fact that PEG/hexanol does not behave as a simple steric media should serve as a caution for studies that use PALES as a quantitative prediction tool (especially for disordered proteins). Further progress in this area depends on our ability to accurately model the anisotropic media and its site-specific interactions with protein molecules. Once this ability is improved, it should be possible to use the alignment parameters as a measure of domain-domain cooperativity, thus identifying the situations where two domains transiently interact with each other or become coupled through a partially structured linker.     

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.     

Artificial proteins as allosteric modulators of PDZ3 and SH3 in two‐domain constructs: A computational characterization of novel chimeric proteins

Artificial multidomain proteins with enhanced structural and functional properties can be utilized in a broad spectrum of applications. The design of chimeric fusion proteins utilizing protein domains or one‐domain miniproteins as building blocks is an important advancement for the creation of new biomolecules for biotechnology and medical applications. However, computational studies to describe in detail the dynamics and geometry properties of two‐domain constructs made from structurally and functionally different proteins are lacking. Here, we tested an in silico design strategy using all‐atom explicit solvent molecular dynamics simulations. The well‐characterized PDZ3 and SH3 domains of human zonula occludens (ZO‐1) (3TSZ), along with 5 artificial domains and 2 types of molecular linkers, were selected to construct chimeric two‐domain molecules. The influence of the artificial domains on the structure and dynamics of the PDZ3 and SH3 domains was determined using a range of analyses. We conclude that the artificial domains can function as allosteric modulators of the PDZ3 and SH3 domains. Proteins 2016; 84:1358–1374. © 2016 Wiley Periodicals, Inc.     

A systematic study of the vibrational free energies of polypeptides in folded and random states

Molecular vibrations, especially low frequency motions, may be used as an indication of the rigidity or the flatness of the protein folding energy landscape. We have studied the vibrational properties of native folded as well as random coil structures of more than 60 polypeptides. The picture we obtain allows us to perceive how and why the energy landscape progressively rigidifies while still allowing potential flexibility. Compared with random coil structures, both alpha-helices and beta-hairpins are vibrationally more flexible. The vibrational properties of loop structures are similar to those of the corresponding random coil structures. Inclusion of an alpha-helix tends to rigidify peptides and so-called building blocks of the structure, whereas the addition of a beta-structure has less effect. When small building blocks coalesce to form larger domains, the protein rigidifies. However, some folded native conformations are still found to be vibrationally more flexible than random coil structures, for example, beta(2)-microglobulin and the SH3 domain. Vibrational free energy contributes significantly to the thermodynamics of protein folding and affects the distribution of the conformational substates. We found a weak correlation between the vibrational folding energy and the protein size, consistent with both previous experimental estimates and theoretical partition of the heat capacity change in protein folding.     

Thermal and chemical denaturation of the BRCT functional module of human 53BP1

BRCTs are protein-docking modules involved in eukaryotic DNA repair. They are characterized by low sequence homology with generally well-conserved structure organization. In a considerable number of proteins, a pair of BRCT structural repeats occurs, connected with inter-BRCT linkers, variable in length, sequence and structure. Linkers may separate and control the relative position of BRCT domains as well as protect and stabilize the hydrophobic inter-BRCT interface region. Their vital role in protein function has been demonstrated by recent findings associating missense mutations in the inter-repeat linker region of the BRCT domain of BRCA1 (BRCA1-BRCT) to hereditary breast/ovarian cancer. The interaction of 53BP1 with the core domain of the p53 tumor suppressor involves the C-terminal BRCT repeat as well as the inert-BRCT linker of the tandem BRCT domain of 53BP1 (53BP1-BRCT). High-accuracy differential scanning calorimetry (DSC) and circular dichroism (CD) have been employed to characterize the heat-induced unfolding of 53BP1-BRCT domain. The calorimetric results provide evidence for unfolding to an intermediate, only partly unfolded state, which, based on the CD results, retains the secondary structural characteristics of the native protein. A direct comparison with the corresponding thermal processes for BRAC1-BRCT and BARD1-BRCT provides evidence that the observed behavior is analogous to BRCA1-BRCT even though the two domains differ substantially in the linker structure. Moreover, chemical denaturation experiments of the untagged 53BP1-BRCT and comparison with BRCA1 and BARD1 BRCTs show that no clear association can be drawn between the structural organization of the inter-BRCT linkers and the overall stability of the BRCT domains.     

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.     

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.     

Quantitative relation between intermolecular and intramolecular binding of pro-rich peptides to SH3 domains

Flexible linkers are often found to tether binding sequence motifs or connect protein domains. Here we analyze three usages of flexible linkers: 1), intramolecular binding of proline-rich peptides (PRPs) to SH3 domains for kinase regulation; 2), intramolecular binding of PRP for increasing the folding stability of SH3 domains; and 3), covalent linking of PRPs and other ligands for high-affinity bivalent binding. The basis of these analyses is a quantitative relation between intermolecular and intramolecular binding constants. This relation has the form K(i) = K(e0)p for intramolecular binding and K(e) = K(e01)K(e02)p for bivalent binding. The effective concentration p depends on the length of the linker and the distance between the linker attachment points in the bound state. Several applications illustrate the usefulness of the quantitative relation. These include intramolecular binding to the Itk SH3 domain by an internal PRP and to a circular permutant of the alpha-spectrin SH3 domain by a designed PRP, and bivalent binding to the two SH3 domains of Grb2 by two linked PRPs. These and other examples suggest that flexible linkers and sequence motifs tethered to them, like folded protein domains, are also subject to tight control during evolution.     

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.