Human cardiac myosin binding protein C: structural flexibility within an extended modular architecture

New insights into the modular organization and flexibility of the N-terminal half of human cardiac myosin binding protein C (cMyBP-C) and information on the association state of the full-length protein have been deduced from a combined small-angle X-ray scattering (SAXS) and NMR study. SAXS data show that the first five immunoglobulin domains of cMyBP-C, which include those implicated in interactions with both myosin and actin, remain monodisperse and monomeric in solution and have a highly extended yet distinctively ‘bent’ modular arrangement that is similar to the giant elastic muscle protein titin. Analyses of the NMR and SAXS data indicate that a proline/alanine-rich linker connecting the cardiac-specific N-terminal C0 domain to the C1 domain provides significant structural flexibility at the N-terminus of the human isoform, while the modular arrangement of domains C1–C2–C3–C4 is relatively fixed. Domain fragments from the C-terminal half of the protein have a propensity to self-associate in vitro, while full-length bacterially expressed cMyBP-C forms flexible extended dimers at micromolar protein concentrations. In summary, our studies reveal that human cMyBP-C combines a distinctive modular architecture with regions of flexibility and that the N-terminal half of the protein is sufficiently extended to span the range of interfilament distances sampled within the dynamic environment of heart muscle. These structural features of cMyBP-C could facilitate its putative role as a molecular switch between actin and myosin and may contribute to modulating the transverse pliancy of the C-zone of the A-band across muscle sarcomeres.     

Stability and folding rates of domains spanning the large A-band super-repeat of titin

Titin is a very large (>3 MDa) protein found in striated muscle where it is believed to participate in myogenesis and passive tension. A prominent feature in the A-band portion of titin is the presence of an 11-domain super-repeat of immunoglobulin superfamily and fibronectin-type-III-like domains. Seven overlapping constructs from human cardiac titin, each consisting of two or three domains and together spanning the entire 11-domain super-repeat, have been expressed in Escherichia coli. Fluorescence unfolding experiments and circular dichroism spectroscopy have been used to measure folding stabilities for each of the constructs and to assign unfolding rates for each super-repeat domain. Immunoglobulin superfamily domains were found to fold correctly only in the presence of their C-terminal fibronectin type II domain, suggesting close and possibly rigid association between these units. The domain stabilities, which range from 8.6 to 42 kJ mol(-1) under physiological conditions, correlate with previously reported mechanical forces required to unfold titin domains. Individual domains vary greatly in their rates of unfolding, with a range of unfolding rate constants between 2.6 x 10(-6) and 1.2 s(-1). This variation in folding behavior is likely to be an important determinant in ensuring independent folding of domains in multi-domain proteins such as titin.     

1H, 15N and 13C backbone chemical shift assignment of the titin A67-A68 domain tandem

Single molecules of the giant protein titin extend across half of the muscle sarcomere, from the Z-line to the M-line, and have roles in muscle assembly and elasticity. In the A-band titin is attached to thick filaments and here the domain arrangement occurs in regular patterns of eleven called the large super-repeat. The large super-repeat itself occurs eleven times and forms nearly half the titin molecule. Interactions of the large super-repeats with myosin are consistent with a role in thick filament assembly. Here we report backbone assignments of the titin A67-A68 domain tandem (Fn-Ig) from the third super-repeat (A65-A75) completed using triple resonance NMR experiments.     

Structural and functional studies of titin's fn3 modules reveal conserved surface patterns and binding to myosin S1--a possible role in the Frank-Starling mechanism of the heart.

The A-band part of titin, a striated-muscle specific protein spanning from the Z-line to the M-line, mainly consists of a well-ordered super-repeat array of immunoglobulin-like and fibronectin-type III (fn3)-like domains. Since it has been suspected that the fn3 domains might represent titin's binding sites to myosin, we have developed structural models for all of titin's 132 fn3-like domains. A subset of eight experimentally determined fn3 structures from a range of proteins, including titin itself, was used as homology templates. After grouping the models according to their position within the super-repeat segment of the central A-band titin region, we analyzed the models with respect to side-chain conservation. This showed that conserved residues form an extensive surface pattern predominantly at one side of the domains, whereas domains outside the central C-zone super-repeat region show generally less conserved surfaces. Since the conserved surface residues may function as protein-binding sites, we experimentally studied the binding properties of expressed multi-domain fn3 fragments. This revealed that fn3 fragments specifically bind to the sub-fragment 1 of myosin. We also measured the effect of fn3 fragments on the contractile properties of single cardiac myocytes. At sub-maximal Ca(2+) concentrations, fn3 fragments significantly enhance active tension. This effect is most pronounced at short sarcomere length, and as a result the length-dependence of Ca(2+) activation is reduced. A model of how titin's fn3-like domains may influence actomyosin interaction is proposed.     

Two immunoglobulin tandem proteins with a linking β-strand reveal unexpected differences in cooperativity and folding pathways

The study of the folding of single domains, in the context of their multidomain environment, is important because more than 70% of eukaryotic proteins are composed of multiple domains. The structures of the tandem immunoglobulin (Ig) domain pairs A164-A165 and A168-A169, from the A-band of the giant muscle protein titin, reveal that they form tightly associated domain arrangements, connected by a continuous β-strand. We investigate the thermodynamic and kinetic properties of these tandem domain pairs. While A164-A165 apparently behaves as a single cooperative unit at equilibrium, unfolding without the accumulation of a large population of intermediates, domains in A168-A169 behave independently. Although A169 appears to be stabilized in the tandem protein, we show that this is due to nonspecific stabilization by extension. We elucidate the folding and unfolding pathways of both tandem pairs and show that cooperativity in A164-A165 is a manifestation of the relative refolding and unfolding rate constants of each individual domain. We infer that the differences between the two tandem pairs result from a different pattern of interactions at the domain/domain interface.     

Unraveling evolutionary constraints: A heterogeneous conservation in dynamics of the titin Ig domains

The giant protein titin, which comprises immunoglobulin (Ig) domains, acts as a bidirectional spring in muscle. The unfolding of Ig domains has been extensively studied, but their dynamics under native states have not been well-characterized. We performed molecular dynamics simulation on a single titin Ig domain and multi-domains. Mobile regions displaying concerted motions were identified. The dynamics of Ig domains are constrained by evolutionary pressures, in such a way that global dominant motion is conserved, yet different flexibilities within Ig domains and in linkers connecting neighbouring domains were observed. We explain these heterogeneous conserved dynamics in relation to sequence conservation across species and the sequence diversity among neighbouring Ig domains.     

Extended and flexible domain solution structure of the extracellular matrix protein anosmin-1 by X-ray scattering, analytical ultracentrifugation and constrained modelling

Kallmann's syndrome corresponds to a loss of sense of smell and hypogonadotrophic hypogonadism. Defects in anosmin-1 result in the X-linked inherited form of Kallmann's syndrome. Anosmin-1 is an extracellular matrix protein comprised of an N-terminal, cysteine-rich (Cys-box) domain and a whey acidic protein-like (WAP) domain, followed by four fibronectin type III (FnIII) domains. The solution structures of recombinant proteins containing the first three domains (PIWF1) and all six domains (PIWF4) were determined by X-ray scattering and analytical ultracentrifugation. Guinier analyses showed that PIWF1 and PIWF4 have different radii of gyration (R(G)) values of 3.1 nm and 6.7 nm, respectively, but similar cross-sectional radii of gyration (R(XS)) values of 1.5 nm and 1.9 nm, respectively. Distance distribution functions showed that the maximum lengths of PIWF1 and PIWF4 were 11 nm and 23 nm, respectively. Analytical ultracentrifugation gave sedimentation coefficients of 2.52 S and 3.55 S for PIWF1 and PIWF4, respectively. The interpretation of the scattering data by constrained modelling requires homology models for all six domains in anosmin-1. While models were already available for the WAP and FnIII domains, searches suggested the Cys-box domain may resemble the cysteine-rich region of the insulin-like growth factor receptor. Automated constrained molecular modelling based on joining the anosmin-1 domains with structurally randomised linkers resulted in 10,000 models for anosmin-1. A trial-and-error search showed that about 0.1-1.4% of these models fitted the X-ray data. The best models showed that the three domains and six domains in PIWF1 and PIWF4, respectively, were extended. The inter-domain linkers in anosmin-1 could not all be extended at the same time, and there was evidence for inter-domain flexibility. Models with folded-back domain arrangements do not fit the data. These solution structures account for the known biological function of anosmin-1, in particular its ability to interact with its three macromolecular ligands.     

Towards a molecular understanding of titin.

Titin is at present the largest known protein (M(r) 3000 kDa) and its expression is restricted to vertebrate striated muscle. Single molecules span from M- to Z-lines and therefore over 1 micron. We have isolated cDNAs encoding five distant titin A-band epitopes, extended their sequences and determined 30 kb (1000 kDa) of the primary structure of titin. Sequences near the M-line encode a kinase domain and are closely related to the C-terminus of twitchin from Caenorhabditis elegans. This suggests that the function of this region in the titin/twitchin family is conserved throughout the animal kingdom. All other A-band sequences consist of 100 amino acid (aa) repeats predicting immunoglobulin-C2 and fibronectin type III globular domains. These domains are arranged into highly ordered 11 domain super-repeat patterns likely to match the myosin helix repeat in the thick filament. Expressed titin fragments bind to the LMM part of myosin and C-protein. Binding strength increases with the number of domains involved, indicating a cumulative effect of multiple binding sites for myosin along the titin molecule. We conclude that A-band titin is likely to be involved in the ordered assembly of the vertebrate thick filament.     

1H, 15N and 13C backbone chemical shift assignment of titin domains A59–A60 and A60 alone

The giant protein titin is the third most abundant protein of vertebrate striated muscle. The titin molecule is >1 μm long and spans half the sarcomere, from the Z-disk to the M-line, and has important roles in sarcomere assembly, elasticity and intracellular signaling. In the A-band of the sarcomere titin is attached to the thick filaments and mainly consists immunoglobulin-like and fibronectin type III-like domains. These are mostly arranged in long-range patterns or ‘super-repeats’. The large super-repeats each contain 11 domains and are repeated 11 times, thus forming nearly half the titin molecule. Through interactions with myosin and C-protein, they are involved in thick filament assembly. The importance of titin in muscle assembly is highlighted by the effect of mutations in the A-band portion, which are the commonest cause of dilated cardiomyopathy, affecting ~1 in 250 (Herman et al. in N Engl J Med 366:619–628, 2012). Here we report backbone ¹⁵N, ¹³C and ¹H chemical shift and ¹³Cβ assignments for the A59–A60 domain tandem from the titin A59–A69 large super-repeat, completed using triple resonance NMR. Since, some regions of the backbone remained unassigned in A60 domain of the complete A59–A60 tandem, a construct containing a single A60 domain, A60sd, was also studied using the same methods. Considerably improved assignment coverage was achieved using A60sd due to its lower mass and improved molecular tumbling rate; these assignments also allowed the analysis of inter-domain interactions using chemical shift mapping against A59–A60.     

Molecular structure of the sarcomeric M band: mapping of titin and myosin binding domains in myomesin and the identification of a potential regulatory phosphorylation site in myomesin.

The M band of sarcomeric muscle is a highly complex structure which contributes to the maintenance of the regular lattice of thick filaments. We propose that the spatial coordination of this assembly is regulated by specific interactions of myosin filaments, the M band protein myomesin and the large carboxy-terminal region of titin. Corresponding binding sites between these proteins were identified. Myomesin binds myosin in the central region of light meromyosin (LMM, myosin residues 1506-1674) by its unique amino-terminal domain My1. A single titin immunoglobulin domain, m4, interacts with a myomesin fragment spanning domains My4-My6. This interaction is regulated by phosphorylation of Ser482 in the linker between myomesin domains My4 and My5. Myomesin phosphorylation at this site by cAMP-dependent kinase and similar or identical activities in muscle extracts block the association with titin. We propose that this demonstration of a phosphorylation-controlled interaction in the sarcomeric cytoskeleton is of potential relevance for sarcomere formation and/or turnover. It also reveals how binding affinities of modular proteins can be regulated by modifications of inter-domain linkers.     

Determination of the Individual Roles of the Linker Residues in the Interdomain Motions of Calmodulin Using NMR Chemical Shifts

Many protein molecules are formed by two or more domains whose structures and dynamics are closely related to their biological functions. It is thus important to develop methods to determine the structural properties of these multidomain proteins. Here, we characterize the interdomain motions in the calcium-bound state of calmodulin (Ca^2 +-CaM) using NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations. We find that the conformational fluctuations of the interdomain linker, which are largely responsible for the overall interdomain motions of CaM, can be well described by exploiting the information provided by chemical shifts. We thus identify 10 residues in the interdomain linker region that change their conformations upon substrate binding. Five of these residues (Met76, Lys77, Thr79, Asp80 and Ser81) are highly flexible and cover the range of conformations observed in the substrate-bound state, while the remaining five (Arg74, Lys75, Asp78, Glu82 and Glu83) are much more rigid and do not populate conformations typical of the substrate-bound form. The ensemble of conformations representing the Ca^2 +-CaM state obtained in this study is in good agreement with residual dipolar coupling, paramagnetic resonance enhancement, small-angle X-ray scattering and fluorescence resonance energy transfer measurements, which were not used as restraints in the calculations. These results provide initial evidence that chemical shifts can be used to characterize the conformational fluctuations of multidomain proteins.     

The Ig doublet Z1Z2: a model system for the hybrid analysis of conformational dynamics in Ig tandems from titin

Titin is a gigantic elastic filament that determines sarcomere ultrastructure and stretch response in vertebrate muscle. It folds into numerous Ig and FnIII domains connected in tandem. Data on interdomain arrangements and dynamics are essential for understanding the function of this filament. Here, we report a mechanistic analysis of the conformational dynamics of two Ig domains from the N terminus of titin, Z1Z2, by using X-ray crystallography, SAXS, NMR relaxation data, and residual dipolar couplings in combination. Z1Z2 preferentially adopts semiextended conformations in solution, with close-hinge arrangements representing low-probability states. Although interdomain contacts are not observed, the linker appears to acquire moderate rigidity via small, local hydrophobic interactions. Thus, Z1Z2 constitutes an adaptable modular system with restricted dynamics. We speculate that its preexistent conformation contributes to the selective recruitment of the binding partner telethonin onto the repetitive surface of the filament. The structural interconversion of four Z1Z2 conformers is analyzed.     

Structural investigation of zymogenic and activated forms of human blood coagulation factor VIII: a computational molecular dynamics study

Background

Human blood coagulation factor VIII (fVIII) is a large plasma glycoprotein with sequential domain arrangement in the order A1-a1-A2-a2-B-a3-A3-C1-C2. The A1, A2 and A3 domains are interconnected by long linker peptides (a1, a2 and a3) that possess the activation sites. Proteolysis of fVIII zymogen by thrombin or factor Xa results in the generation of the activated form (fVIIIa) which serves as a critical co-factor for factor IXa (fIXa) enzyme in the intrinsic coagulation pathway.

Results

In our efforts to elucidate the structural differences between fVIII and fVIIIa, we developed the solution structural models of both forms, starting from an incomplete 3.7 Å X-ray crystal structure of fVIII zymogen, using explicit solvent MD simulations. The full assembly of B-domainless single-chain fVIII was built between the A1-A2 (Ala1-Arg740) and A3-C1-C2 (Ser1669-Tyr2332) domains. The structural dynamics of fVIII and fVIIIa, simulated for over 70 ns of time scale, enabled us to evaluate the integral motions of the multi-domain assembly of the co-factor and the possible coordination pattern of the functionally important calcium and copper ion binding in the protein.

Conclusions

MD simulations predicted that the acidic linker peptide (a1) between the A1 and A2 domains is largely flexible and appears to mask the exposure of putative fIXa enzyme binding loop (Tyr555-Asp569) region in the A2 domain. The simulation of fVIIIa, generated from the zymogen structure, predicted that the linker peptide (a1) undergoes significant conformational reorganization upon activation by relocating completely to the A1-domain. The conformational transition led to the exposure of the Tyr555-Asp569 loop and the surrounding region in the A2 domain. While the proposed linker peptide conformation is predictive in nature and warrants further experimental validation, the observed conformational differences between the zymogen and activated forms may explain and support the large body of experimental data that implicated the critical importance of the cleavage of the peptide bond between the Arg372 and Ser373 residues for the full co-factor activity of fVIII.     

Mechanical Design of the Third FnIII Domain of Tenascin-C

By combining single-molecule atomic force microscopy (AFM), proline mutagenesis and steered molecular dynamics (SMD) simulations, we investigated the mechanical unfolding dynamics and mechanical design of the third fibronectin type III domain of tenascin-C (TNfn3) in detail. We found that the mechanical stability of TNfn3 is similar to that of other constituting FnIII domains of tenascin-C, and the unfolding process of TNfn3 is an apparent two-state process. By employing proline mutagenesis to block the formation of backbone hydrogen bonds and introduce structural disruption in β sheet, we revealed that in addition to the important roles played by hydrophobic core packing, backbone hydrogen bonds in β hairpins are also responsible for the overall mechanical stability of TNfn3. Furthermore, proline mutagenesis revealed that the mechanical design of TNfn3 is robust and the mechanical stability of TNfn3 is very resistant to structural disruptions caused by proline substitutions in β sheets. Proline mutant F88P is one exception, as the proline mutation at position 88 reduced the mechanical stability of TNfn3 significantly and led to unfolding forces of < 20 pN. This result suggests that Phe88 is a weak point of the mechanical resistance for TNfn3. We used SMD simulations to understand the molecular details underlying the mechanical unfolding of TNfn3. The comparison between the AFM results and SMD simulations revealed similarities and discrepancies between the two. We compared the mechanical unfolding and design of TNfn3 and its structural homologue, the tenth FnIII domain from fibronectin. These results revealed the complexity underlying the mechanical design of FnIII domains and will serve as a starting point for systematically analyzing the mechanical architecture of other FnIII domains in tenascins-C, and will help to gain a better understanding of some of the complex features observed for the stretching of native tenascin-C.     

Crystal structure of the C1 domain of cardiac myosin binding protein-C: implications for hypertrophic cardiomyopathy

C-protein is a major component of skeletal and cardiac muscle thick filaments. Mutations in the gene encoding cardiac C-protein [cardiac myosin binding protein-C (cMyBP-C)] are one of the principal causes of hypertrophic cardiomyopathy. cMyBP-C is a string of globular domains including eight immunoglobulin-like and three fibronectin-like domains termed C0-C10. It binds to myosin and titin, and probably to actin, and may have both a structural and a regulatory role in muscle function. To help to understand the pathology of the known mutations, we have solved the structure of the immunoglobulin-like C1 domain of MyBP-C by X-ray crystallography to a resolution of 1.55 Å. Mutations associated with hypertrophic cardiomyopathy are clustered at one end towards the C-terminus, close to the important C1C2 linker, where they alter the structural integrity of this region and its interactions.     

Functional chimeras of the phosphodiesterase 5 and 10 tandem GAF domains

The tandem GAF domain of hPDE10A uses cAMP as an allosteric ligand (Gross-Langenhoff, M., Hofbauer, K., Weber, J., Schultz, A., and Schultz, J. E. (2006) J. Biol. Chem. 281, 2841-2846). We used a two-pronged approach to study how discrimination of ligand is achieved in human (h)PDE10A and how domain selection in the phosphodiesterase GAF tandems is determined. First, we examined which functional groups of cAMP are responsible for purine ring discrimination. Changes at the C-6 ring position (removal of the amino group; chloride substitution) and at the N-1 ring position reduced stimulation efficacy by 80%, i.e. marking those positions as decisive for nucleotide discrimination. Second, we generated a GAF tandem chimera that consisted of the cGMP-binding GAF-A unit from hPDE5A1, which signals through cGMP in PDE5, and the GAF-B from hPDE10A1, which signals through cAMP in PDE10. Stimulation of the reporter enzyme exclusively was through the GAF-B domain of hPDE10A1 (EC(50) = 7 mum cAMP) as shown by respective point mutations. The PDE5 GAF-A domain in the chimera did not signal, and its function was reduced to a strictly structural role. Signaling was independent of the origin of the N terminus. Generating 10 additional PDE5/10 tandem GAF chimeras surprisingly demonstrated that the length-conserved linker in GAF tandems between GAF-A and GAF-B played an unforeseen decisive role in intramolecular signaling. Swapping the linker sections between PDE5 and PDE10 GAF tandem domains abrogated signaling completely pointing to specific domain interactions within GAF tandems, which are not visible in the available crystal structures with bound ligands.     

Molecular Characterisation of Titin N2A and Its Binding of CARP Reveals a Titin/Actin Cross-linking Mechanism

Striated muscle responds to mechanical overload by rapidly up-regulating the expression of the cardiac ankyrin repeat protein, CARP, which then targets the sarcomere by binding to titin N2A in the I-band region. To date, the role of this interaction in the stress response of muscle remains poorly understood. Here, we characterise the molecular structure of the CARP-receptor site in titin (UN2A) and its binding of CARP. We find that titin UN2A contains a central three-helix bundle fold (ca 45 residues in length) that is joined to N- and C-terminal flanking immunoglobulin domains by long, flexible linkers with partial helical content. CARP binds titin by engaging an α-hairpin in the three-helix fold of UN2A, the C-terminal linker sequence, and the BC loop in Ig81, which jointly form a broad binding interface. Mutagenesis showed that the CARP/N2A association withstands sequence variations in titin N2A and we use this information to evaluate 85 human single nucleotide variants. In addition, actin co-sedimentation, co-transfection in C2C12 cells, proteomics on heart lysates, and the mechanical response of CARP-soaked myofibrils imply that CARP induces the cross-linking of titin and actin myofilaments, thereby increasing myofibril stiffness. We conclude that CARP acts as a regulator of force output in the sarcomere that preserves muscle mechanical performance upon overload stress.     

Studies on the interaction between titin and myosin.

This study examines the interaction of titin and myosin. In order to analyze the domains of myosin contributing to the binding for titin, we conducted a solid phase binding assay. Different portions of myosin (heavy chains, light chains and myosin fragments) were coated on the microtiter wells and reacted with biotinylated titin. Then the binding of biotinylated titin to these polypeptides was detected by using the avidinbiotin-peroxidase method. The results demonstrated that light meromyosin and subfragment 1 were the major domains of myosin interacting with titin. Titin fragments obtained by trypsin digestion were allowed to react with myosin in an affinity column, and the bound fragments were isolated by an acidic elution. Immunoblot analysis of myosin-bound titin fragments revealed that an A-band domain of titin was responsible for the binding of myosin. In addition, biotinylated titin labelled the outer A-bands and Z-bands in intact myofibrils, thus confirming the in situ binding of titin to myosin.     

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.     

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.