Stable single α-helices are constant force springs in proteins

In silico and in vivo analyses of a novel variant in MYO6 identified in a family with postlingual non-syndromic hearing loss from Argentina

Hereditary hearing loss stands as the most prevalent sensory disorder, with over 124 non-syndromic genes and approximately 400 syndromic forms of deafness identified in humans. The clinical presentation of these conditions spans a spectrum, ranging from mild to profound hearing loss. The aim of this study was to identify the genetic cause of hearing loss in a family and functionally validate a novel variant identified in the MYO6 gene. After Whole Exome Sequencing analysis, the variant c.2775G>C p.Arg925Ser in MYO6 was detected in a family with postlingual non-syndromic hearing loss. By protein modeling a change in the electrostatic charge of the single alpha helix domain surface was revealed. Through a knockdown phenotype rescue assay in zebrafish, the detrimental effects of the identified variant on the auditory system was determined. These findings underscore the significance of a comprehensive approach, integrating both in silico and in vivo strategies, to ascertain the pathogenicity of this candidate variant. Such an approach has demonstrated its effectiveness in achieving an accurate genetic diagnosis and in promoting a more profound comprehension of the mechanisms that underlie the pathophysiology of hearing.

Backbone Hydration of α-Helical Peptides: Hydrogen-Bonding and Surface Hydrophobicity/Hydrophilicity

The stability of proteins and small peptides depends on the way they interact with the surrounding water molecules. For small peptides, such as α-helical polyalanine (polyALA), water molecules can weaken the intramolecular hydrogen-bonds (HB) formed between the peptide backbone O and NH groups which are responsible for the α-helix structure. Here, we perform molecular dynamics simulations to study the hydration of polyALA, polyserine (polySER), and other homopolymer peptide α-helices at different temperatures and pressures. We find that water molecules form HB with most polyALA carbonyl O atoms, despite ALA hydrophobic CH3 side chain. Similar water-peptide backbone HB are found in other (hydrophobic and hydrophilic) homopolymer α-helices with large side chains, including polyvaline, polyleucine, and polyphenyalanine. A novel hydration mechanism is observed in polyserine (polySER): the backbone peptide rarely forms HB with water and, instead, the carbonyl O atoms tend to form HB with polySER side chain OH groups. We also quantify the hydrophobicity/hydrophilicity of polyALA and polySER by calculating the contact angle θ c of a water droplet pierced by a long polyALA/polySER α-helix. Unexpectedly, even when polyALA α-helix is supposed to be hydrophobic (θ c > 90°), we find that θ c ≈ 79°. For polySER, θ c ≈ 70°, consistent with α-helical polySER being hydrophilic.

Multi-integrated approach for unraveling small open reading frames potentially associated with secondary metabolism in Streptomyces

IMPORTANCE

Due to their small size and special chemical features, small open reading frame (smORF)-encoding peptides (SEPs) are often neglected. However, they may play critical roles in regulating gene expression, enzyme activity, and metabolite production. Studies on bacterial microproteins have mainly focused on pathogenic bacteria, which are importance to systematically investigate SEPs in streptomycetes and are rich sources of bioactive secondary metabolites. Our study is the first to perform a global identification of smORFs in streptomycetes. We established a peptidogenomic workflow for non-model microbial strains and identified multiple novel smORFs that are potentially linked to secondary metabolism in streptomycetes. Our multi-integrated approach in this study is meaningful to improve the quality and quantity of the detected smORFs. Ultimately, the workflow we established could be extended to other organisms and would benefit the genome mining of microproteins with critical functions for regulation and engineering useful microorganisms.

A novel single alpha-helix DNA-binding domain in CAF-1 promotes gene silencing and DNA damage survival through tetrasome-length DNA selectivity and spacer function

The histone chaperone chromatin assembly factor 1 (CAF-1) deposits two nascent histone H3/H4 dimers onto newly replicated DNA forming the central core of the nucleosome known as the tetrasome. How CAF-1 ensures there is sufficient space for the assembly of tetrasomes remains unknown. Structural and biophysical characterization of the lysine/glutamic acid/arginine-rich (KER) region of CAF-1 revealed a 128-Å single alpha-helix (SAH) motif with unprecedented DNA-binding properties. Distinct KER sequence features and length of the SAH drive the selectivity of CAF-1 for tetrasome-length DNA and facilitate function in budding yeast. In vivo, the KER cooperates with the DNA-binding winged helix domain in CAF-1 to overcome DNA damage sensitivity and maintain silencing of gene expression. We propose that the KER SAH links functional domains within CAF-1 with structural precision, acting as a DNA-binding spacer element during chromatin assembly.

Diversity of viscoelastic properties of an engineered muscle-inspired protein hydrogel

Folded protein hydrogels are prime candidates as tuneable biomaterials but it is unclear to what extent their mechanical properties have mesoscopic, as opposed to molecular origins. To address this, we probe hydrogels inspired by the muscle protein titin and engineered to the polyprotein I27₅, using a multimodal rheology approach. Across multiple protocols, the hydrogels consistently exhibit power-law viscoelasticity in the linear viscoelastic regime with an exponent β = 0.03, suggesting a dense fractal meso-structure, with predicted fractal dimension d_f = 2.48. In the nonlinear viscoelastic regime, the hydrogel undergoes stiffening and energy dissipation, indicating simultaneous alignment and unfolding of the folded proteins on the nanoscale. Remarkably, this behaviour is highly reversible, as the value of β, d_f and the viscoelastic moduli return to their equilibrium value, even after multiple cycles of deformation. This highlights a previously unrevealed diversity of viscoelastic properties that originate on both at the nanoscale and the mesoscopic scale, providing powerful opportunities for engineering novel biomaterials.

PK1 from Drosophila obscurin is an inactive pseudokinase with scaffolding properties

Obscurins are large filamentous proteins with crucial roles in the assembly, stability and regulation of muscle. Characteristic of these proteins is a tandem of two C-terminal kinase domains, PK1 and PK2, that are separated by a long intrinsically disordered sequence. The significance of this conserved domain arrangement is unknown. Our study of PK1 from Drosophila obscurin shows that this is a pseudokinase with features typical of the CAM-kinase family, but which carries a minimalistic regulatory tail that no longer binds calmodulin or has mechanosensory properties typical of other sarcomeric kinases. PK1 binds ATP with high affinity, but in the absence of magnesium and lacks detectable phosphotransfer activity. It also has a highly diverged active site, strictly conserved across arthropods, that might have evolved to accommodate an unconventional binder. We find that PK1 interacts with PK2, suggesting a functional relation to the latter. These findings lead us to speculate that PK1/PK2 form a pseudokinase/kinase dual system, where PK1 might act as an allosteric regulator of PK2 and where mechanosensing properties, akin to those described for regulatory tails in titin-like kinases, might now reside on the unstructured interkinase segment. We propose that the PK1-interkinase-PK2 region constitutes an integrated functional unit in obscurin proteins.

Quantification of carboxylate-bridged di-zinc site stability in protein due ferri by single-molecule force spectroscopy

Carboxylate-bridged diiron proteins belong to a protein family involved in different physiological processes. These proteins share the conservative EXXH motif, which provides the carboxylate bridge and is critical for metal binding. Here, we choose de novo-designed single-chain due ferri protein (DFsc), a four-helical protein with two EXXH motifs as a model protein, to study the stability of the carboxylate-bridged di-metal binding site. The mechanical and kinetic properties of the di-Zn site in DFsc were obtained by atomic force microscopy-based single-molecule force spectroscopy. Zn-DFsc showed a considerable rupture force of ~200 pN, while the apo-protein is mechanically labile. In addition, multiple rupture pathways were observed with different probabilities, indicating the importance of the EXXH-based carboxylate-bridged metal site. These results demonstrate carboxylate-bridged di-metal site is mechanically stable and improve our understanding of this important type of metalloprotein.

Contribution of Hydrophobic Interactions to Protein Mechanical Stability

The role of hydrophobic and polar interactions in providing thermodynamic stability to folded proteins has been intensively studied, but the relative contribution of these interactions to the mechanical stability is less explored. We used steered molecular dynamics simulations with constant-velocity pulling to generate force-extension curves of selected protein domains and monitor hydrophobic surface unravelling upon extension. Hydrophobic contribution was found to vary between one fifth and one third of the total force while the rest of the contribution is attributed primarily to hydrogen bonds. Moreover, hydrophobic force peaks were shifted towards larger protein extensions with respect to the force peaks attributed to hydrogen bonds. The higher importance of hydrogen bonds compared to hydrophobic interactions in providing mechanical resistance is in contrast with the relative importance of the hydrophobic interactions in providing thermodynamic stability of proteins. The different contributions of these interactions to the mechanical stability are explained by the steeper free energy dependence of hydrogen bonds compared to hydrophobic interactions on the relative positions of interacting atoms. Comparative analyses for several protein domains revealed that the variation of hydrophobic forces is modest, while the contribution of hydrogen bonds to the force peaks becomes increasingly important for mechanically resistant protein domains.

Switching an active site helix in dihydrofolate reductase reveals limits to subdomain modularity

To what degree are individual structural elements within proteins modular such that similar structures from unrelated proteins can be interchanged? We study subdomain modularity by creating 20 chimeras of an enzyme, Escherichia coli dihydrofolate reductase (DHFR), in which a catalytically important, 10-residue α-helical sequence is replaced by α-helical sequences from a diverse set of proteins. The chimeras stably fold but have a range of diminished thermal stabilities and catalytic activities. Evolutionary coupling analysis indicates that the residues of this α-helix are under selection pressure to maintain catalytic activity in DHFR. Reversion to phenylalanine at key position 31 was found to partially restore catalytic activity, which could be explained by evolutionary coupling values. We performed molecular dynamics simulations using replica exchange with solute tempering. Chimeras with low catalytic activity exhibit nonhelical conformations that block the binding site and disrupt the positioning of the catalytically essential residue D27. Simulation observables and in vitro measurements of thermal stability and substrate-binding affinity are strongly correlated. Several E. coli strains with chromosomally integrated chimeric DHFRs can grow, with growth rates that follow predictions from a kinetic flux model that depends on the intracellular abundance and catalytic activity of DHFR. Our findings show that although α-helices are not universally substitutable, the molecular and fitness effects of modular segments can be predicted by the biophysical compatibility of the replacement segment.

Titin kinase ubiquitination aligns autophagy receptors with mechanical signals in the sarcomere

Striated muscle undergoes remodelling in response to mechanical and physiological stress, but little is known about the integration of such varied signals in the myofibril. The interaction of the elastic kinase region from sarcomeric titin (A168-M1) with the autophagy receptors Nbr1/p62 and MuRF E3 ubiquitin ligases is well suited to link mechanosensing with the trophic response of the myofibril. To investigate the mechanisms of signal cross-talk at this titin node, we elucidated its 3D structure, analysed its response to stretch using steered molecular dynamics simulations and explored its functional relation to MuRF1 and Nbr1/p62 using cellular assays. We found that MuRF1-mediated ubiquitination of titin kinase promotes its scaffolding of Nbr1/p62 and that the process can be dynamically down-regulated by the mechanical unfolding of a linker sequence joining titin kinase with the MuRF1 receptor site in titin. We propose that titin ubiquitination is sensitive to the mechanical state of the sarcomere, the regulation of sarcomere targeting by Nbr1/p62 being a functional outcome. We conclude that MuRF1/Titin Kinase/Nbr1/p62 constitutes a distinct assembly that predictably promotes sarcomere breakdown in inactive muscle.

When Familial Hearing Loss Means Genetic Heterogeneity: A Model Case Report

We describe a family with both hearing loss (HL) and thrombocytopenia, caused by pathogenic variants in three genes. The proband was a child with neonatal thrombocytopenia, childhood-onset HL, hyper-laxity and severe myopia. The child’s mother (and some of her relatives) presented with moderate thrombocytopenia and adulthood-onset HL. The child’s father (and some of his relatives) presented with adult-onset HL. An HL panel analysis, completed by whole exome sequencing, was performed in this complex family. We identified three pathogenic variants in three different genes: MYH9, MYO7A and ACTG1. The thrombocytopenia in the child and her mother is explained by the MYH9 variant. The post-lingual HL in the paternal branch is explained by the MYO7A variant, absent in the proband, while the congenital HL of the child is explained by a de novo ACTG1 variant. This family, in which HL segregates, illustrates that multiple genetic conditions coexist in individuals and make patient care more complex than expected.

Disentangling the complexity of low complexity proteins

There are multiple definitions for low complexity regions (LCRs) in protein sequences, with all of them broadly considering LCRs as regions with fewer amino acid types compared to an average composition. Following this view, LCRs can also be defined as regions showing composition bias. In this critical review, we focus on the definition of sequence complexity of LCRs and their connection with structure. We present statistics and methodological approaches that measure low complexity (LC) and related sequence properties. Composition bias is often associated with LC and disorder, but repeats, while compositionally biased, might also induce ordered structures. We illustrate this dichotomy, and more generally the overlaps between different properties related to LCRs, using examples. We argue that statistical measures alone cannot capture all structural aspects of LCRs and recommend the combined usage of a variety of predictive tools and measurements. While the methodologies available to study LCRs are already very advanced, we foresee that a more comprehensive annotation of sequences in the databases will enable the improvement of predictions and a better understanding of the evolution and the connection between structure and function of LCRs. This will require the use of standards for the generation and exchange of data describing all aspects of LCRs.

Short abstract

There are multiple definitions for low complexity regions (LCRs) in protein sequences. In this critical review, we focus on the definition of sequence complexity of LCRs and their connection with structure. We present statistics and methodological approaches that measure low complexity (LC) and related sequence properties. Composition bias is often associated with LC and disorder, but repeats, while compositionally biased, might also induce ordered structures. We illustrate this dichotomy, plus overlaps between different properties related to LCRs, using examples.

Structure of the shutdown state of myosin-2

Myosin-2 is essential for processes as diverse as cell division and muscle contraction. Dephosphorylation of its regulatory light chain (RLC) promotes an inactive, ‘shutdown’ state with the filament-forming tail folded onto the two heads¹, preventing filament formation and inactivating the motors². The mechanism by which this happens is obscure. Here we report a cryo-electron microscopy structure of shutdown smooth muscle myosin, with a resolution of 6 Å in the head region. A pseudo-atomic model, obtained by flexible fitting of crystal structures into the density and molecular dynamics simulations, describes interaction interfaces at the atomic level. The N-terminal extension of one RLC interacts with the tail and the other with the partner head, revealing how the RLCs stabilise the shutdown state in different ways and how their phosphorylation would allow myosin activation. Additional interactions between the three segments of the coiled coil, the motor domains and LCs stabilise the shutdown molecule. The structure of the lever in each head is competent to generate force upon activation. This shutdown structure is relevant to all myosin-2 isoforms and provides a framework for understanding their disease-causing mutations.

ER/K-link-Leveraging a native protein linker to probe dynamic cellular interactions

ER/K α-helices are a subset of single alpha helical domains, which exhibit unusual stability as isolated protein secondary structures. They adopt an elongated structural conformation, while regulating the frequency of interactions between proteins or polypeptides fused to their ends. Here we review recent advances on the structure, stability and function of ER/K α-helices as linkers (ER/K linkers) in native proteins. We describe methodological considerations in the molecular cloning, protein expression and measurement of interaction strengths, using sensors incorporating ER/K linkers. We highlight biological insights obtained over the last decade by leveraging distinct biophysical features of ER/K-linked sensors. We conclude with the outlook for the use of ER/K linkers in the selective modulation of dynamic cellular interactions.

Unfolding mechanism and free energy landscape of single, stable, alpha helices at low pull speeds

Single alpha helices (SAHs) stable in isolated form are often found in motor proteins where they bridge functional domains. Understanding the mechanical response of SAHs is thus critical to understand their function. The quasi-static force-extension relation of a small number of SAHs is known from single-molecule experiments. Unknown, or still controversial, are the molecular scale details behind those observations. We show that the deformation mechanism of SAHs pulled from the termini at pull speeds approaching the quasi-static limit differs from that of typical helices found in proteins, which are stable only when interacting with other protein domains. Using molecular dynamics simulations with atomistic resolution at low pull speeds previously inaccessible to simulation, we show that SAHs start unfolding from the termini at all pull speeds we investigated. Unfolding proceeds residue-by-residue and hydrogen bond breaking is not the main event determining the barrier to unfolding. We use the molecular simulation data to test the cooperative sticky chain model. This model yields excellent fits of the force-extension curves and quantifies the distance, xE = 0.13 nm, to the transition state, the natural frequency of bond vibration, ν0 = 0.82 ns-1, and the height, V0 = 2.9 kcal mol-1, of the free energy barrier associated with the deformation of single residues. Our results demonstrate that the sticky chain model could advantageously be used to analyze experimental force-extension curves of SAHs and other biopolymers.

Single-Molecule Force Spectroscopy on the N2A Element of Titin: Effects of Phosphorylation and CARP

Titin is a large filamentous protein that forms a sarcomeric myofilament with a molecular spring region that develops force in stretched sarcomeres. The molecular spring has a complex make-up that includes the N2A element. This element largely consists of a 104-residue unique sequence (N2A-Us) flanked by immunoglobulin domains (I80 and I81). The N2A element is of interest because it assembles a signalosome with CARP (Cardiac Ankyrin Repeat Protein) as an important component; CARP both interacts with the N2A-Us and I81 and is highly upregulated in response to mechanical stress. The mechanical properties of the N2A element were studied using single-molecule force spectroscopy, including how these properties are affected by CARP and phosphorylation. Three protein constructs were made that consisted of 0, 1, or 2 N2A-Us elements with flanking I80 and I81 domains and with specific handles at their ends for study by atomic force microscopy (AFM). The N2A-Us behaved as an entropic spring with a persistence length (Lp) of ∼0.35 nm and contour length (Lc) of ∼39 nm. CARP increased the Lp of the N2A-Us and the unfolding force of the Ig domains; force clamp experiments showed that CARP reduced the Ig domain unfolding kinetics. These findings suggest that CARP might function as a molecular chaperone that protects I81 from unfolding when mechanical stress is high. The N2A-Us was found to be a PKA substrate, and phosphorylation was blocked by CARP. Mass spectrometry revealed a PKA phosphosite (Ser-9895 in NP_001254479.2) located at the border between the N2A-Us and I81. AFM studies showed that phosphorylation affected neither the Lp of the N2A-Us nor the Ig domain unfolding force (Funfold). Simulating the force-sarcomere length relation of a single titin molecule containing all spring elements showed that the compliance of the N2A-Us only slightly reduces passive force (1.4%) with an additional small reduction by CARP (0.3%). Thus, it is improbable that the compliance of the N2A element has a mechanical function per se. Instead, it is likely that this compliance has local effects on binding of signaling molecules and that it contributes thereby to strain- and phosphorylation- dependent mechano-signaling.

Modulating the Stiffness of the Myosin VI Single α-Helical Domain

Highly charged, single α-helical (SAH) domains contain a high percentage of Arg, Lys, and Glu residues. Their dynamic salt bridge pairing creates the exceptional stiffness of these helical rods, with a persistence length of more than 200 Å for the myosin VI SAH domain. With the aim of modulating the stiffness of the helical structure, we investigated the effect, using NMR spectroscopy, of substituting key charged Arg, Lys, Glu, and Asp residues by Gly or His. Results indicate that such mutations result in the transient breaking of the helix at the site of mutation but with noticeable impact on amide hydrogen exchange rates extending as far as ±2 helical turns, pointing to a substantial degree of cooperativity in SAH stability. Whereas a single Gly substitution caused transient breaks ∼20% of the time, two consecutive Gly substitutions break the helix ∼65% of the time. NMR relaxation measurements indicate that the exchange rate between an intact and a broken helix is fast (>300,000 s^-1) and that for the wild-type sequence, the finite persistence length is dominated by thermal fluctuations of backbone torsion angles and H-bond lengths, not by transient helix breaking. The double mutation D27H/E28H causes a pH-dependent fraction of helix disruption, in which the helix breakage increases from 26% at pH 7.5 to 53% at pH 5.5. The ability to modulate helical integrity by pH may enable incorporation of externally tunable dynamic components in the design of molecular machines.

Helix 8 is the essential structural motif of mechanosensitive GPCRs

G-protein coupled receptors (GPCRs) are versatile cellular sensors for chemical stimuli, but also serve as mechanosensors involved in various (patho)physiological settings like vascular regulation, cardiac hypertrophy and preeclampsia. However, the molecular mechanisms underlying mechanically induced GPCR activation have remained elusive. Here we show that mechanosensitive histamine H₁ receptors (H₁Rs) are endothelial sensors of fluid shear stress and contribute to flow-induced vasodilation. At the molecular level, we observe that H₁Rs undergo stimulus-specific patterns of conformational changes suggesting that mechanical forces and agonists induce distinct active receptor conformations. GPCRs lacking C-terminal helix 8 (H8) are not mechanosensitive, and transfer of H8 to non-responsive GPCRs confers, while removal of H8 precludes, mechanosensitivity. Moreover, disrupting H8 structural integrity by amino acid exchanges impairs mechanosensitivity. Altogether, H8 is the essential structural motif endowing GPCRs with mechanosensitivity. These findings provide a mechanistic basis for a better understanding of the roles of mechanosensitive GPCRs in (patho)physiology.

GPCRs are versatile cellular sensors for chemical stimuli but the molecular mechanisms underlying mechanically induced GPCR activation have remained elusive. Here authors identify the C-terminal helix 8 (H8) as the essential structural motif endowing H₁R and other GPCRs with mechanosensitivity.

The hierarchical emergence of worm-like chain behaviour from globular domain polymer chains

Biological organisms make use of hierarchically organised structures to modulate mechanical behaviour across multiple lengthscales, allowing microscopic objects to generate macroscopic effects. Within these structural hierarchies, the resultant physical behaviour of the entire system is determined not only by the intrinsic mechanical properties of constituent subunits, but also by their organisation in three-dimensional space. When these subunits are polyproteins, colloidal chains or other globular domain polymers, the Kratky-Porod model is often assumed for the individual subunits. Hence, it is implicitly asserted that the polymeric object has an intrinsic parameter, the persistence length, that defines its flexibility. However, the persistence lengths extracted from experiment vary, and are often relatively small. Through a series of simulations on polymer chains formed of globular subunits, we show that the persistence length itself is a hierarchical structural property, related not only to the intrinsic mechanical properties of the underlying monomeric subunits, but emerging due to the organisation of inhomogenous geometry along the polymer contour.

Probing the Basis of α-Synuclein Aggregation by Comparing Simulations to Single-Molecule Experiments

Intrinsically disordered proteins often play an important role in protein aggregation. However, it is challenging to determine the structures and interactions that drive the early stages of aggregation because they are transient and obscured in a heterogeneous mixture of disordered states. Even computational methods are limited because the lack of ordered structure makes it difficult to ensure that the relevant conformations are sampled. We address these challenges by integrating atomistic simulations with high-resolution single-molecule measurements reported previously, using the measurements to help discern which parts of the disordered ensemble of structures in the simulations are most probable while using the simulations to identify residues and interactions that are important for oligomer stability. This approach was applied to α-synuclein, an intrinsically disordered protein that aggregates in the context of Parkinson's disease. We simulated single-molecule pulling experiments on dimers, the minimal oligomer, and compared them to force spectroscopy measurements. Force-extension curves were simulated starting from a set of 66 structures with substantial structured content selected from the ensemble of dimer structures generated at zero force via Monte Carlo simulations. The pattern of contour length changes as the structures unfolded through intermediate states was compared to the results from optical trapping measurements on the same dimer to discern likely structures occurring in the measurements. Simulated pulling curves were generally consistent with experimental data but with a larger number of transient intermediates. We identified an ensemble of β-rich dimer structures consistent with the experimental data from which dimer interfaces could be deduced. These results suggest specific druggable targets in the structural motifs of α-synuclein that may help prevent the earliest steps of oligomerization.

Remarkable Rigidity of the Single α-Helical Domain of Myosin-VI As Revealed by NMR Spectroscopy

Although the α-helix has long been recognized as an all-important element of secondary structure, it generally requires stabilization by tertiary interactions with other parts of a protein's structure. Highly charged single α-helical (SAH) domains, consisting of a high percentage (>75%) of Arg, Lys, and Glu residues, are exceptions to this rule but have been difficult to characterize structurally. Our study focuses on the 68-residue medial tail domain of myosin-VI, which is found to contain a highly ordered α-helical structure extending from Glu-6 to Lys-63. High hydrogen exchange protection factors (15-150), small (ca. 4 Hz) 3 JHNHα couplings, and a near-perfect fit to an ideal model α-helix for its residual dipolar couplings (RDCs), measured in a filamentous phage medium, support the high regularity of this helix. Remarkably, the hydrogen exchange rates are far more homogeneous than the protection factors derived from them, suggesting that for these transiently broken helices the intrinsic exchange rates derived from the amino acid sequence are not appropriate reference values. 15N relaxation data indicate a very high degree of rotational diffusion anisotropy ( D∥/ D⊥ ≈ 7.6), consistent with the hydrodynamic behavior predicted for such a long, nearly straight α-helix. Alignment of the helix by a paramagnetic lanthanide ion attached to its N-terminal region shows a decrease in alignment as the distance from the tagging site increases. This decrease yields a precise measure for the persistence length of 224 ± 10 Å at 20 °C, supporting the idea that the role of the SAH helix is to act as an extension of the myosin-VI lever arm.

Dynamic ion pair behavior stabilizes single α-helices in proteins

Ion pairs are key stabilizing interactions between oppositely charged amino acid side chains in proteins. They are often depicted as single conformer salt bridges (hydrogen-bonded ion pairs) in crystal structures, but it is unclear how dynamic they are in solution. Ion pairs are thought to be particularly important in stabilizing single α-helix (SAH) domains in solution. These highly stable domains are rich in charged residues (such as Arg, Lys, and Glu) with potential ion pairs across adjacent turns of the helix. They provide a good model system to investigate how ion pairs can contribute to protein stability. Using NMR spectroscopy, small-angle X-ray light scattering (SAXS), and molecular dynamics simulations, we provide here experimental evidence that ion pairs exist in a SAH in murine myosin 7a (residues 858-935), but that they are not fixed or long lasting. In silico modeling revealed that the ion pairs within this α-helix exhibit dynamic behavior, rapidly forming and breaking and alternating between different partner residues. The low-energy helical state was compatible with a great variety of ion pair combinations. Flexible ion pair formation utilizing a subset of those available at any one time avoided the entropic penalty of fixing side chain conformations, which likely contributed to helix stability overall. These results indicate the dynamic nature of ion pairs in SAHs. More broadly, thermodynamic stability in other proteins is likely to benefit from the dynamic behavior of multi-option solvent-exposed ion pairs.

Helical Polyampholyte Sequences Have Unique Thermodynamic Properties

Helices are the most common structural pattern observed in structured proteins. Polypeptide sequences that form helices in isolation have been identified and extensively studied. These are generally rich in alanine, the amino acid with strongest helical propensity. Insertion of charged or polar amino acids has been shown to be necessary to make alanine-rich peptides soluble and sometimes even increase the helicity of the peptides. More recently sequences that contain mostly charged residues (E-R/K rich) have been found in naturally occurring proteins that are highly helical, soluble, and extended regardless their length. Artificial sequences composed mostly or exclusively of charged amino acids have been designed that are also highly helical, depending on the specific pattern of oppositely charged residues. Here we explore the thermodynamic properties of a number of 16-residue long peptides with varying helical propensity by performing equilibrium simulations over a broad range of temperatures. We observe quantitative differences in the peptides' helical propensities that can be related to qualitative differences in the free energy landscape, depending on the ampholytic patterns in the sequence. The results provide hints on how the specific physical properties of naturally occurring long sequences with similar patterns of charged residues may relate to their biological function.

Detection of single alpha-helices in large protein sequence sets using hardware acceleration

Single alpha-helices (SAHs) are increasingly recognized as important structural and functional elements of proteins. Comprehensive identification of SAH segments in large protein datasets was largely hindered by the slow speed of the most restrictive prediction tool for their identification, FT_CHARGE on common hardware. We have previously implemented an FPGA-based version of this tool allowing fast analysis of a large number of sequences. Using this implementation, we have set up of a semi-automated pipeline capable of analyzing full UniProt releases in reasonable time and compiling monthly updates of a comprehensive database of SAH segments. Releases of this database, denoted CSAHDB, is available on the CSAHserver 2 website at csahserver.itk.ppke.hu. An overview of human SAH-containing sequences combined with a literature survey suggests specific roles of SAH segments in proteins involved in RNA-based regulation processes as well as cytoskeletal proteins, a number of which is also linked to the development and function of synapses.

Deciphering the Mechanical Properties of Type III Secretion System EspA Protein by Single Molecule Force Spectroscopy

Bacterial pathogens inject virulence factors into host cells during bacterial infections using type III secretion systems. In enteropathogenic Escherichia coli, this system contains an external filament, formed by a self-oligomerizing protein called E. coli secreted protein A (EspA). The EspA filament penetrates the thick viscous mucus layer to facilitate the attachment of the bacteria to the gut-epithelium. To do that, the EspA filament requires noteworthy mechanical endurance considering the mechanical shear stresses found within the intestinal tract. To date, the mechanical properties of the EspA filament and the structural and biophysical knowledge of monomeric EspA are very limited, mostly due to the strong tendency of the protein to self-oligomerize. To overcome this limitation, we employed a single molecule force spectroscopy (SMFS) technique and studied the mechanical properties of EspA. Force extension dynamic of (I91)₄-EspA-(I91)₄ chimera revealed two structural unfolding events occurring at low forces during EspA unfolding, thus indicating no unique mechanical stability of the monomeric protein. SMFS examination of purified monomeric EspA protein, treated by a gradually refolding protocol, exhibited similar mechanical properties as the EspA protein within the (I91)₄-EspA-(I91)₄ chimera. Overall, our results suggest that the mechanical integrity of the EspA filament likely originates from the interactions between EspA monomers and not from the strength of an individual monomer.

Determining Stable Single Alpha Helical (SAH) Domain Properties by Circular Dichroism and Atomic Force Microscopy

Stable, single α-helical (SAH) domains exist in a number of unconventional myosin isoforms, as well as other proteins. These domains are formed from sequences rich in charged residues (Arg, Lys, and Glu), they can be hundreds of residues long, and in isolation they can tolerate significant changes in pH and salt concentration without loss in helicity. Here we describe methods for the preparation and purification of SAH domains and SAH domain-containing constructs, using the myosin 10 SAH domain as an example. We go on to describe the use of circular dichroism spectroscopy and force spectroscopy with the atomic force microscope for the elucidation of structural and mechanical properties of these unusual helical species.

Distribution and evolution of stable single α-helices (SAH domains) in myosin motor proteins

Stable single-alpha helices (SAHs) are versatile structural elements in many prokaryotic and eukaryotic proteins acting as semi-flexible linkers and constant force springs. This way SAH-domains function as part of the lever of many different myosins. Canonical myosin levers consist of one or several IQ-motifs to which light chains such as calmodulin bind. SAH-domains provide flexibility in length and stiffness to the myosin levers, and may be particularly suited for myosins working in crowded cellular environments. Although the function of the SAH-domains in human class-6 and class-10 myosins has well been characterised, the distribution of the SAH-domain in all myosin subfamilies and across the eukaryotic tree of life remained elusive. Here, we analysed the largest available myosin sequence dataset consisting of 7919 manually annotated myosin sequences from 938 species representing all major eukaryotic branches using the SAH-prediction algorithm of Waggawagga, a recently developed tool for the identification of SAH-domains. With this approach we identified SAH-domains in more than one third of the supposed 79 myosin subfamilies. Depending on the myosin class, the presence of SAH-domains can range from a few to almost all class members indicating complex patterns of independent and taxon-specific SAH-domain gain and loss.

Characterization of long and stable de novo single alpha-helix domains provides novel insight into their stability

Naturally-occurring single α-helices (SAHs), are rich in Arg (R), Glu (E) and Lys (K) residues, and stabilized by multiple salt bridges. Understanding how salt bridges promote their stability is challenging as SAHs are long and their sequences highly variable. Thus, we designed and tested simple de novo 98-residue polypeptides containing 7-residue repeats (AEEEXXX, where X is K or R) expected to promote salt-bridge formation between Glu and Lys/Arg. Lys-rich sequences (EK3 (AEEEKKK) and EK2R1 (AEEEKRK)) both form SAHs, of which EK2R1 is more helical and thermo-stable suggesting Arg increases stability. Substituting Lys with Arg (or vice versa) in the naturally-occurring myosin-6 SAH similarly increased (or decreased) its stability. However, Arg-rich de novo sequences (ER3 (AEEERRR) and EK1R2 (AEEEKRR)) aggregated. Combining a PDB analysis with molecular modelling provides a rational explanation, demonstrating that Glu and Arg form salt bridges more commonly, utilize a wider range of rotamer conformations, and are more dynamic than Glu-Lys. This promiscuous nature of Arg helps explain the increased propensity of de novo Arg-rich SAHs to aggregate. Importantly, the specific K:R ratio is likely to be important in determining helical stability in de novo and naturally-occurring polypeptides, giving new insight into how single α-helices are stabilized.

Scaffolding in the Spliceosome via Single α Helices

The spliceosomal B complex-specific protein Prp38 forms a complex with the intrinsically unstructured proteins MFAP1 and Snu23. Our binding and crystal structure analyses show that MFAP1 and Snu23 contact Prp38 via ER/K motif-stabilized single α helices, which have previously been recognized only as rigid connectors or force springs between protein domains. A variant of the Prp38-binding single α helix of MFAP1, in which ER/K motifs not involved in Prp38 binding were mutated, was less α-helical in isolation and showed a reduced Prp38 affinity, with opposing tendencies in interaction enthalpy and entropy. Our results indicate that the strengths of single α helix-based interactions can be tuned by the degree of helix stabilization in the unbound state. MFAP1, Snu23, and several other spliceosomal proteins contain multiple regions that likely form single α helices via which they might tether several binding partners and act as intermittent scaffolds that facilitate remodeling steps during assembly of an active spliceosome.

Nanomechanical properties of MscL α helices: A steered molecular dynamics study

Gating of mechanosensitive (MS) channels is driven by a hierarchical cascade of movements and deformations of transmembrane helices in response to bilayer tension. Determining the intrinsic mechanical properties of the individual transmembrane helices is therefore central to understanding the intricacies of the gating mechanism of MS channels. We used a constant-force steered molecular dynamics (SMD) approach to perform unidirectional pulling tests on all the helices of MscL in M. tuberculosis and E. coli homologs. Using this method, we could overcome the issues encountered with the commonly used constant-velocity SMD simulations, such as low mechanical stability of the helix during stretching and high dependency of the elastic properties on the pulling rate. We estimated Young's moduli of the α-helices of MscL to vary between 0.2 and 12.5 GPa with TM2 helix being the stiffest. We also studied the effect of water on the properties of the pore-lining TM1 helix. In the absence of water, this helix exhibited a much stiffer response. By monitoring the number of hydrogen bonds, it appears that water acts like a 'lubricant' (softener) during TM1 helix elongation. These data shed light on another physical aspect underlying hydrophobic gating of MS channels, in particular MscL.

CARP interacts with titin at a unique helical N2A sequence and at the domain Ig81 to form a structured complex

The cardiac ankyrin repeat protein (CARP) is up-regulated in the myocardium during cardiovascular disease and in response to mechanical or toxic stress. Stress-induced CARP interacts with the N2A spring region of the titin filament to modulate muscle compliance. We characterize the interaction between CARP and titin-N2A and show that the binding site in titin spans the dual domain UN2A-Ig81. We find that the unique sequence UN2A is not structurally disordered, but that it has a stable, elongated α-helical fold that possibly acts as a constant force spring. Our findings portray CARP/titin-N2A as a structured node and help to rationalize the molecular basis of CARP mechanosensing in the sarcomeric I-band.

The myosin X motor is optimized for movement on actin bundles

Myosin X has features not found in other myosins. Its structure must underlie its unique ability to generate filopodia, which are essential for neuritogenesis, wound healing, cancer metastasis and some pathogenic infections. By determining high-resolution structures of key components of this motor, and characterizing the in vitro behaviour of the native dimer, we identify the features that explain the myosin X dimer behaviour. Single-molecule studies demonstrate that a native myosin X dimer moves on actin bundles with higher velocities and takes larger steps than on single actin filaments. The largest steps on actin bundles are larger than previously reported for artificially dimerized myosin X constructs or any other myosin. Our model and kinetic data explain why these large steps and high velocities can only occur on bundled filaments. Thus, myosin X functions as an antiparallel dimer in cells with a unique geometry optimized for movement on actin bundles.

A Combination of Diffusion and Active Translocation Localizes Myosin 10 to the Filopodial Tip

Myosin 10 is an actin-based molecular motor that localizes to the tips of filopodia in mammalian cells. To understand how it is targeted to this distinct region of the cell, we have used total internal reflection fluorescence microscopy to study the movement of individual full-length and truncated GFP-tagged molecules. Truncation mutants lacking the motor region failed to localize to filopodial tips but still bound transiently at the plasma membrane. Deletion of the single α-helical and anti-parallel coiled-coil forming regions, which lie between the motor and pleckstrin homology domains, reduced the instantaneous velocity of intrafilopodial movement but did not affect the number of substrate adherent filopodia. Deletion of the anti-parallel coiled-coil forming region, but not the EKR-rich region of the single α-helical domain, restored intrafilopodial trafficking, suggesting this region is important in determining myosin 10 motility. We propose a model by which myosin 10 rapidly targets to the filopodial tip via a sequential reduction in dimensionality. Molecules first undergo rapid diffusion within the three-dimensional volume of the cell body. They then exhibit periods of slower two-dimensional diffusion in the plane of the plasma membrane. Finally, they move in a unidimensional, highly directed manner along the polarized actin filament bundle within the filopodium becoming confined to a single point at the tip. Here we have observed directly each phase of the trafficking process using single molecule fluorescence imaging of live cells and have quantified our observations using single particle tracking, autocorrelation analysis, and kymographs.

Functional and Structural Characterization of Novel Type of Linker Connecting Capsid and Nucleocapsid Protein Domains in Murine Leukemia Virus

The assembly of immature retroviral particles is initiated in the cytoplasm by the binding of the structural polyprotein precursor Gag with viral genomic RNA. The protein interactions necessary for assembly are mediated predominantly by the capsid (CA) and nucleocapsid (NC) domains, which have conserved structures. In contrast, the structural arrangement of the CA-NC connecting region differs between retroviral species. In HIV-1 and Rous sarcoma virus, this region forms a rod-like structure that separates the CA and NC domains, whereas in Mason-Pfizer monkey virus, this region is densely packed, thus holding the CA and NC domains in close proximity. Interestingly, the sequence connecting the CA and NC domains in gammaretroviruses, such as murine leukemia virus (MLV), is unique. The sequence is called a charged assembly helix (CAH) due to a high number of positively and negatively charged residues. Although both computational and deletion analyses suggested that the MLV CAH forms a helical conformation, no structural or biochemical data supporting this hypothesis have been published. Using an in vitro assembly assay, alanine scanning mutagenesis, and biophysical techniques (circular dichroism, NMR, microcalorimetry, and electrophoretic mobility shift assay), we have characterized the structure and function of the MLV CAH. We provide experimental evidence that the MLV CAH belongs to a group of charged, E(R/K)-rich, single α-helices. This is the first single α-helix motif identified in viral proteins.

Differential Effects of Hydrophobic Core Packing Residues for Thermodynamic and Mechanical Stability of a Hyperthermophilic Protein

Proteins from organisms that have adapted to environmental extremes provide attractive systems to explore and determine the origins of protein stability. Improved hydrophobic core packing and decreased loop-length flexibility can increase the thermodynamic stability of proteins from hyperthermophilic organisms. However, their impact on protein mechanical stability is not known. Here, we use protein engineering, biophysical characterization, single-molecule force spectroscopy (SMFS), and molecular dynamics (MD) simulations to measure the effect of altering hydrophobic core packing on the stability of the cold shock protein TmCSP from the hyperthermophilic bacterium Thermotoga maritima. We make two variants of TmCSP in which a mutation is made to reduce the size of aliphatic groups from buried hydrophobic side chains. In the first, a mutation is introduced in a long loop (TmCSP L40A); in the other, the mutation is introduced on the C-terminal β-strand (TmCSP V62A). We use MD simulations to confirm that the mutant TmCSP L40A shows the most significant increase in loop flexibility, and mutant TmCSP V62A shows greater disruption to the core packing. We measure the thermodynamic stability (ΔGD-N) of the mutated proteins and show that there is a more significant reduction for TmCSP L40A (ΔΔG = 63%) than TmCSP V62A (ΔΔG = 47%), as might be expected on the basis of the relative reduction in the size of the side chain. By contrast, SMFS measures the mechanical stability (ΔG*) and shows a greater reduction for TmCSP V62A (ΔΔG* = 8.4%) than TmCSP L40A (ΔΔG* = 2.5%). While the impact on the mechanical stability is subtle, the results demonstrate the power of tuning noncovalent interactions to modulate both the thermodynamic and mechanical stability of a protein. Such understanding and control provide the opportunity to design proteins with optimized thermodynamic and mechanical properties.

The physics of pulling polyproteins: a review of single molecule force spectroscopy using the AFM to study protein unfolding

One of the most exciting developments in the field of biological physics in recent years is the ability to manipulate single molecules and probe their properties and function. Since its emergence over two decades ago, single molecule force spectroscopy has become a powerful tool to explore the response of biological molecules, including proteins, DNA, RNA and their complexes, to the application of an applied force. The force versus extension response of molecules can provide valuable insight into its mechanical stability, as well as details of the underlying energy landscape. In this review we will introduce the technique of single molecule force spectroscopy using the atomic force microscope (AFM), with particular focus on its application to study proteins. We will review the models which have been developed and employed to extract information from single molecule force spectroscopy experiments. Finally, we will end with a discussion of future directions in this field.

Tuning protein mechanics through an ionic cluster graft from an extremophilic protein

Proteins from extremophilic organisms provide excellent model systems to determine the role of non-covalent interactions in defining protein stability and dynamics as well as being attractive targets for the development of robust biomaterials. Hyperthermophilic proteins have a prevalence of salt bridges, relative to their mesophilic homologues, which are thought to be important for enhanced thermal stability. However, the impact of salt bridges on the mechanical properties of proteins is far from understood. Here, a combination of protein engineering, biophysical characterisation, single molecule force spectroscopy (SMFS) and molecular dynamics (MD) simulations directly investigates the role of salt bridges in the mechanical stability of two cold shock proteins; BsCSP from the mesophilic organism Bacillus subtilis and TmCSP from the hyperthermophilic organism Thermotoga maritima. Single molecule force spectroscopy shows that at ambient temperatures TmCSP is mechanically stronger yet, counter-intuitively, its native state can withstand greater deformation before unfolding (i.e. it is mechanically soft) compared with BsCSP. MD simulations were used to identify the location and quantify the population of salt bridges, and reveal that TmCSP contains a larger number of highly occupied salt bridges than BsCSP. To test the hypothesis that salt-bridges endow these mechanical properties on the hyperthermophilic CSP, a charged triple mutant (CTM) variant of BsCSP was generated by grafting an ionic cluster from TmCSP into the BsCSP scaffold. As expected CTM is thermodynamically more stable and mechanically softer than BsCSP. We show that a grafted ionic cluster can increase the mechanical softness of a protein and speculate that it could provide a mechanical recovery mechanism and that it may be a design feature applicable to other proteins.

A conserved charged single α-helix with a putative steric role in paraspeckle formation

Paraspeckles are subnuclear particles involved in the regulation of mRNA expression. They are formed by the association of DBHS family proteins and the NEAT1 long noncoding RNA. Here, we show that a recently identified structural motif, the charged single α-helix, is largely conserved in the DBHS family. Based on the available structural data and a previously suggested multimerization scheme of DBHS proteins, we built a structural model of a (PSPC1/NONO)(n) multimer that might have relevance in paraspeckle formation. Our model contains an extended coiled-coil region that is followed by and partially overlaps with the predicted charged single α-helix. We suggest that the charged single α-helix can act as an elastic ruler governing the exact positioning of the dimeric core structures relative to each other during paraspeckle assembly along the NEAT1 noncoding RNA.

The Inner Centromere Protein (INCENP) Coil Is a Single α-Helix (SAH) Domain That Binds Directly to Microtubules and Is Important for Chromosome Passenger Complex (CPC) Localization and Function in Mitosis

The chromosome passenger complex (CPC) is a master regulator of mitosis. Inner centromere protein (INCENP) acts as a scaffold regulating CPC localization and activity. During early mitosis, the N-terminal region of INCENP forms a three-helix bundle with Survivin and Borealin, directing the CPC to the inner centromere where it plays essential roles in chromosome alignment and the spindle assembly checkpoint. The C-terminal IN box region of INCENP is responsible for binding and activating Aurora B kinase. The central region of INCENP has been proposed to comprise a coiled coil domain acting as a spacer between the N- and C-terminal domains that is involved in microtubule binding and regulation of the spindle checkpoint. Here we show that the central region (213 residues) of chicken INCENP is not a coiled coil but a ∼ 32-nm-long single α-helix (SAH) domain. The N-terminal half of this domain directly binds to microtubules in vitro. By analogy with previous studies of myosin 10, our data suggest that the INCENP SAH might stretch up to ∼ 80 nm under physiological forces. Thus, the INCENP SAH could act as a flexible "dog leash," allowing Aurora B to phosphorylate dynamic substrates localized in the outer kinetochore while at the same time being stably anchored to the heterochromatin of the inner centromere. Furthermore, by achieving this flexibility via an SAH domain, the CPC avoids a need for dimerization (required for coiled coil formation), which would greatly complicate regulation of the proximity-induced trans-phosphorylation that is critical for Aurora B activation.

Myosin tails and single α-helical domains

The human genome contains 39 myosin genes, divided up into 12 different classes. The structure, cellular function and biochemical properties of many of these isoforms remain poorly characterized and there is still some controversy as to whether some myosin isoforms are monomers or dimers. Myosin isoforms 6 and 10 contain a stable single α-helical (SAH) domain, situated just after the canonical lever. The SAH domain is stiff enough to be able to lengthen the lever allowing the myosin to take a larger step. In addition, atomic force microscopy and atomistic simulations show that SAH domains unfold at relatively low forces and have a high propensity to refold. These properties are likely to be important for protein function, enabling motors to carry cargo in dense actin networks, and other proteins to remain attached to binding partners in the crowded cell.

Effect of external pulling forces on the length distribution of peptides

BACKGROUND

The distribution of the length of a polypeptide, or that of the distance between any two of its atoms, is an important property as it can be analytically or numerically estimated for a number of polymer models. Importantly, it is directly measurable through a number of different experimental techniques. Length distributions can be straightforwardly assessed from molecular dynamics simulation; however, true convergence through full accurate coverage of the length range is difficult to achieve.

METHODS

The application of external constant force combined with the weighted-histogram analysis method (WHAM) is used to enhance sampling of unlikely 'long' or 'short' conformations and obtain the potential of mean force, while also collecting dynamic properties of the chain under variable tension.

RESULTS

We demonstrate the utility of constant force to enhance the sampling efficiency and obtain experimentally measurable quantities on a series of short peptides, including charge-rich sequences that are known to be highly helical but whose properties are distinct from those of helical peptides undergoing helix-coil transitions.

CONCLUSIONS

Force-enhanced sampling enhances the range and accuracy of the length-based potential of mean force of the peptide, in particular those sequences that contain increased numbers of charged residues.

GENERAL SIGNIFICANCE

This approach allows users to simultaneously probe the force-dependent behaviour of peptides directly, enhance the range and accuracy of the length-based PMF of the peptide and also test the convergence of simulations by comparing the overlap of PMF profiles from different constant forces. This article is part of a special issue entitled Recent developments of molecular dynamics.