Local sequence‐structure relationships in proteins

We seek to understand the interplay between amino acid sequence and local structure in proteins. Are some amino acids unique in their ability to fit harmoniously into certain local structures? What is the role of sequence in sculpting the putative native state folds from myriad possible conformations? In order to address these questions, we represent the local structure of each C_α atom of a protein by just two angles, θ and μ, and we analyze a set of more than 4,000 protein structures from the PDB. We use a hierarchical clustering scheme to divide the 20 amino acids into six distinct groups based on their similarity to each other in fitting local structural space. We present the results of a detailed analysis of patterns of amino acid specificity in adopting local structural conformations and show that the sequence‐structure correlation is not very strong compared with a random assignment of sequence to structure. Yet, our analysis may be useful to determine an effective scoring rubric for quantifying the match of an amino acid to its putative local structure.     

From Discrete to Continuum Models of Three-Dimensional Deformations in Epithelial Sheets

Metamorphic proteins, including proteins with high levels of sequence identity but different folds, are exceptions to the long-standing rule-of-thumb that proteins with as little as 30% sequence identity adopt the same fold. Which topologies can be bridged by these highly identical sequences remains an open question. Here we bridge two 3-α-helix bundle proteins with two radically different folds. Using a straightforward approach, we engineered the sequences of one subdomain within maltose binding protein (MBP, α/β/α-sandwich) and another within outer surface protein A (OspA, β-sheet) to have high sequence identity (80 and 77%, respectively) with engineered variants of protein G (G_A, 3-α-helix bundle). Circular dichroism and nuclear magnetic resonance spectra of all engineered variants demonstrate that they maintain their native conformations despite substantial sequence modification. Furthermore, the MBP variant (80% identical to G_A) remained active. Thermodynamic analysis of numerous G_A and MBP variants suggests that the key to our approach involved stabilizing the modified MBP and OspA subdomains via external interactions with neighboring substructures, indicating that subdomain interactions can stabilize alternative folds over a broad range of sequence variation. These findings suggest that it is possible to bridge one fold with many other topologies, which has implications for protein folding, evolution, and misfolding diseases.     

Harnessing the Unique Structural Properties of Isolated α-Helices

The α-helix is a ubiquitous secondary structural element that is almost exclusively observed in proteins when stabilized by tertiary or quaternary interactions. However, beginning with the unexpected observations of α-helix formation in the isolated C-peptide in ribonuclease A, there is growing evidence that a significant percentage (0.2%) of all proteins contain isolated stable single α-helical domains (SAH). These SAH domains provide unique structural features essential for normal protein function. A subset of SAH domains contain a characteristic ER/K motif, composed of a repeating sequence of ∼4 consecutive glutamic acids followed by ∼4 consecutive basic arginine or lysine (R/K) residues. The ER/K α-helix, also termed the ER/K linker, has been extensively characterized in the context of the myosin family of molecular motors and is emerging as a versatile structural element for protein and cellular engineering applications. Here, we review the structure and function of SAH domains, as well as the tools to identify them in natural proteins. We conclude with a discussion of recent studies that have successfully used the modular ER/K linker for engineering chimeric myosin proteins with altered mechanical properties, as well as synthetic polypeptides that can be used to monitor and systematically modulate protein interactions within cells.     

Effect of Short- and Long-Range Interactions on Trp Rotamer Populations Determined by Site-Directed Tryptophan Fluorescence of Tear Lipocalin

In the lipocalin family, the conserved interaction between the main α-helix and the β-strand H is an ideal model to study protein side chain dynamics. Site-directed tryptophan fluorescence (SDTF) has successfully elucidated tryptophan rotamers at positions along the main alpha helical segment of tear lipocalin (TL). The rotamers assigned by fluorescent lifetimes of Trp residues corroborate the restriction expected based on secondary structure. Steric conflict constrains Trp residues to two (t, g ⁻) of three possible χ₁ (t, g ⁻, g ⁺) canonical rotamers. In this study, investigation focused on the interplay between rotamers for a single amino acid position, Trp 130 on the α-helix and amino acids Val 113 and Leu 115 on the H strand, i.e. long range interactions. Trp130 was substituted for Phe by point mutation (F130W). Mutations at positions 113 and 115 with combinations of Gly, Ala, Phe residues alter the rotamer distribution of Trp130. Mutations, which do not distort local structure, retain two rotamers (two lifetimes) populated in varying proportions. Replacement of either long range partner with a small amino acid, V113A or L115A, eliminates the dominance of the t rotamer. However, a mutation that distorts local structure around Trp130 adds a third fluorescence lifetime component. The results indicate that the energetics of long-range interactions with Trp 130 further tune rotamer populations. Diminished interactions, evident in W130G113A115, result in about a 22% increase of α-helix content. The data support a hierarchic model of protein folding. Initially the secondary structure is formed by short-range interactions. TL has non-native α-helix intermediates at this stage. Then, the long-range interactions produce the native fold, in which TL shows α-helix to β-sheet transitions. The SDTF method is a valuable tool to assess long-range interaction energies through rotamer distribution as well as the characterization of low-populated rotameric states of functionally important excited protein states.     

The Effect of a ΔK280 Mutation on the Unfolded State of a Microtubule-Binding Repeat in Tau

Tau is a natively unfolded protein that forms intracellular aggregates in the brains of patients with Alzheimer''s disease. To decipher the mechanism underlying the formation of tau aggregates, we developed a novel approach for constructing models of natively unfolded proteins. The method, energy-minima mapping and weighting (EMW), samples local energy minima of subsequences within a natively unfolded protein and then constructs ensembles from these energetically favorable conformations that are consistent with a given set of experimental data. A unique feature of the method is that it does not strive to generate a single ensemble that represents the unfolded state. Instead we construct a number of candidate ensembles, each of which agrees with a given set of experimental constraints, and focus our analysis on local structural features that are present in all of the independently generated ensembles. Using EMW we generated ensembles that are consistent with chemical shift measurements obtained on tau constructs. Thirty models were constructed for the second microtubule binding repeat (MTBR2) in wild-type (WT) tau and a ΔK280 mutant, which is found in some forms of frontotemporal dementia. By focusing on structural features that are preserved across all ensembles, we find that the aggregation-initiating sequence, PHF6*, prefers an extended conformation in both the WT and ΔK280 sequences. In addition, we find that residue K280 can adopt a loop/turn conformation in WT MTBR2 and that deletion of this residue, which can adopt nonextended states, leads to an increase in locally extended conformations near the C-terminus of PHF6*. As an increased preference for extended states near the C-terminus of PHF6* may facilitate the propagation of β-structure downstream from PHF6*, these results explain how a deletion at position 280 can promote the formation of tau aggregates.     

pH responsiveness of fibrous assemblies of repeat-sequence amphipathic α-helix polypeptides

We reported previously that our designed polypeptide α3 (21 residues), which has three repeats of a seven-amino-acid sequence (LETLAKA)₃, forms not only an amphipathic α-helix structure but also long fibrous assemblies in aqueous solution. To address the relationship between the electrical states of the polypeptide and its α-helix and fibrous assembly formation, we characterized mutated polypeptides in which charged amino acid residues of α3 were replaced with Ser. We prepared the following polypeptides: 2Sα3 (LSTLAKA)₃, in which all Glu residues were replaced with Ser residues; 6Sα3 (LETLASA)₃, in which all Lys residues were replaced with Ser; and 2S6Sα3 (LSTLASA)₃; in which all Glu and Lys residues were replaced with Ser. In 0.1M KCl, 2Sα3 formed an α-helix under basic conditions and 6Sα3 formed an α-helix under acid conditions. In 1M KCl, they both formed α-helices under a wide pH range. In addition, 2Sα3 and 6Sα3 formed fibrous assemblies under the same buffer conditions in which they formed α-helices. α-Helix and fibrous assembly formation by these polypeptides was reversible in a pH-dependent manner. In contrast, 2S6Sα3 formed an α-helix under basic conditions in 1M KCl. Taken together, these findings reveal that the charge states of the charged amino acid residues and the charge state of the Leu residue located at the terminus play an important role in α-helix formation.     

Dynamic Motions of the HIV-1 Frameshift Site RNA

The HIV-1 frameshift site (FS) plays a critical role in viral replication. During translation, the HIV-1 FS transitions from a 3-helix to a 2-helix junction RNA secondary structure. The 2-helix junction structure contains a GGA bulge, and purine-rich bulges are common motifs in RNA secondary structure. Here, we investigate the dynamics of the HIV-1 FS 2-helix junction RNA. Interhelical motions were studied under different ionic conditions using NMR order tensor analysis of residual dipolar couplings. In 150 mM potassium, the RNA adopts a 43°(±4°) interhelical bend angle (β) and displays large amplitude, anisotropic interhelical motions characterized by a 0.52(±0.04) internal generalized degree of order (GDO_int) and distinct order tensor asymmetries for its two helices (η = 0.26(±0.04) and 0.5(±0.1)). These motions are effectively quenched by addition of 2 mM magnesium (GDO_int = 0.87(±0.06)), which promotes a near-coaxial conformation (β = 15°(±6°)) of the two helices. Base stacking in the bulge was investigated using the fluorescent purine analog 2-aminopurine. These results indicate that magnesium stabilizes extrahelical conformations of the bulge nucleotides, thereby promoting coaxial stacking of helices. These results are highly similar to previous studies of the HIV transactivation response RNA, despite a complete lack of sequence similarity between the two RNAs. Thus, the conformational space of these RNAs is largely determined by the topology of their interhelical junctions.     

Conformational contagion in a protein: structural properties of a chameleon sequence

Certain sequences, known as chameleon sequences, take both alpha- and beta-conformations in natural proteins. We demonstrate that a wild chameleon sequence fused to the C-terminal alpha-helix or beta-sheet in foreign stable proteins from hyperthermophiles forms the same conformation as the host secondary structure. However, no secondary structural formation is observed when the sequence is attached to the outside of the secondary structure. These results indicate that this sequence inherently possesses an ability to make either alpha- or beta-conformation, depending on the sequentially neighboring secondary structure if little other nonlocal interaction occurs. Thus, chameleon sequences take on a satellite state through contagion by the power of a secondary structure. We propose this "conformational contagion" as a new nonlocal determinant factor in protein structure and misfolding related to protein conformational diseases.     

Analysis of Chameleon Sequences by Energy Decomposition on a Pairwise Per-residue Basis

The conversion from α-helix to β-strand that has been widely observed in so-called chameleon sequences has received considerable attention since such a structural change may induce many amyloidogenic proteins to self-assemble into fibrils thus causing fatal diseases. Here we report a large scale-analysis of the energetics of secondary structural conversions in a collection of chameleon sequences retrieved from the Protein Data Bank. Major energetic contributions to the secondary structural conversion were analyzed by carrying out energy decomposition on a pairwise per-residue basis, i.e., (i,i), (i,i ± 1), (i,i ± 2), (i,i ± 3), (i,i ± 4) and > (i,i ± 4) intra-/inter-residual interactions. While the overall potential energy differences were subtle, individual residue-based interacting energy differences were observed to vary significantly depending on the specific type of secondary structural conversion. The average energy difference between α-helix and β-strand, <ΔE _α→β>, in the chameleon sequences varied significantly in (i,i), (i,i ± 1) and > (i,i ± 4) interactions. The major energetic factors in secondary structure conversions were electrostatic interactions and the polar term for solvation energy. In addition, residue-based average energy differences in α-helix → β-strand conversions were well-correlated to those in α-helix → random coil → β-strand conversions (R ² = 0.92). Assuming that three secondary structural elements can transform in either direction, this strong correlation indicates that the present energy decomposition method using database structures of chameleon sequences provides a reliable tool for the characterization of secondary structure fluctuations in amino acid sequences.     

Loop L5 Assumes Three Distinct Orientations during the ATPase Cycle of the Mitotic Kinesin Eg5: A TRANSIENT AND TIME-RESOLVED FLUORESCENCE STUDY*

Members of the kinesin superfamily of molecular motors differ in several key structural domains, which probably allows these molecular motors to serve the different physiologies required of them. One of the most variable of these is a stem-loop motif referred to as L5. This loop is longest in the mitotic kinesin Eg5, and previous structural studies have shown that it can assume different conformations in different nucleotide states. However, enzymatic domains often consist of a mixture of conformations whose distribution shifts in response to substrate binding or product release, and this information is not available from the “static” images that structural studies provide. We have addressed this issue in the case of Eg5 by attaching a fluorescent probe to L5 and examining its fluorescence, using both steady state and time-resolved methods. This reveals that L5 assumes an equilibrium mixture of three orientations that differ in their local environment and segmental mobility. Combining these studies with transient state kinetics demonstrates that there is a major shift in this distribution during transitions that interconvert weak and strong microtubule binding states. Finally, in conjunction with previous cryo-EM reconstructions of Eg5·microtubule complexes, these fluorescence studies suggest a model in which L5 regulates both nucleotide and microtubule binding through a set of reversible interactions with helix α3. We propose that these features facilitate the production of sustained opposing force by Eg5, which underlies its role in supporting formation of a bipolar spindle in mitosis.     

Aromaticity at position 39 in α‐synuclein: A modulator of amyloid fibril assembly and membrane‐bound conformations

Recent studies revealed that molecular events related with the physiology and pathology of αS might be regulated by specific sequence motifs in the primary sequence of αS. The importance of individual residues in these motifs remains an important open avenue of investigation. In this work, we have addressed the structural details related to the amyloid fibril assembly and lipid‐binding features of αS through the design of site‐directed mutants at position 39 of the protein and their study by in vitro and in vivo assays. We demonstrated that aromaticity at position 39 of αS primary sequence influences strongly the aggregation properties and the membrane‐bound conformations of the protein, molecular features that might have important repercussions for the function and dysfunction of αS. Considering that aggregation and membrane damage is an important driver of cellular toxicity in amyloid diseases, future work is needed to link our findings with studies based on toxicity and neuronal cell death.Brief statement outlining significanceModulation by distinct sequential motifs and specific residues of αS on its physiological and pathological states is an active area of research. Here, we demonstrated that aromaticity at position 39 of αS modulates the membrane‐bound conformations of the protein, whereas removal of aromatic functionality at position 39 reduces strongly the amyloid assembly in vitro and in vivo. Our study provides new evidence for the modulation of molecular events related with the physiology and pathology of αS.     

Structure of a Conserved Golgi Complex-targeting Signal in Coronavirus Envelope Proteins

Coronavirus envelope (CoV E) proteins are ∼100-residue polypeptides with at least one channel-forming α-helical transmembrane (TM) domain. The extramembrane C-terminal tail contains a completely conserved proline, at the center of a predicted β-coil-β motif. This hydrophobic motif has been reported to constitute a Golgi-targeting signal or a second TM domain. However, no structural data for this or other extramembrane domains in CoV E proteins is available. Herein, we show that the E protein in the severe acute respiratory syndrome virus has only one TM domain in micelles, whereas the predicted β-coil-β motif forms a short membrane-bound α-helix connected by a disordered loop to the TM domain. However, complementary results suggest that this motif is potentially poised for conformational change or in dynamic exchange with other conformations.     

The Common Structural Architecture of Shigella flexneri and Salmonella typhimurium Type Three Secretion Needles

The Type Three Secretion System (T3SS), or injectisome, is a macromolecular infection machinery present in many pathogenic Gram-negative bacteria. It consists of a basal body, anchored in both bacterial membranes, and a hollow needle through which effector proteins are delivered into the target host cell. Two different architectures of the T3SS needle have been previously proposed. First, an atomic model of the Salmonella typhimurium needle was generated from solid-state NMR data. The needle subunit protein, PrgI, comprises a rigid-extended N-terminal segment and a helix-loop-helix motif with the N-terminus located on the outside face of the needle. Second, a model of the Shigella flexneri needle was generated from a high-resolution 7.7-Å cryo-electron microscopy density map. The subunit protein, MxiH, contains an N-terminal α-helix, a loop, another α-helix, a 14-residue-long β-hairpin (Q51–Q64) and a C-terminal α-helix, with the N-terminus facing inward to the lumen of the needle. In the current study, we carried out solid-state NMR measurements of wild-type Shigella flexneri needles polymerized in vitro and identified the following secondary structure elements for MxiH: a rigid-extended N-terminal segment (S2-T11), an α-helix (L12-A38), a loop (E39-P44) and a C-terminal α-helix (Q45-R83). Using immunogold labeling in vitro and in vivo on functional needles, we located the N-terminus of MxiH subunits on the exterior of the assembly, consistent with evolutionary sequence conservation patterns and mutagenesis data. We generated a homology model of Shigella flexneri needles compatible with both experimental data: the MxiH solid-state NMR chemical shifts and the state-of-the-art cryoEM density map. These results corroborate the solid-state NMR structure previously solved for Salmonella typhimurium PrgI needles and establish that Shigella flexneri and Salmonella typhimurium subunit proteins adopt a conserved structure and orientation in their assembled state. Our study reveals a common structural architecture of T3SS needles, essential to understand T3SS-mediated infection and develop treatments.     

Edge Strand Dissociation and Conformational Changes in Transthyretin under Amyloidogenic Conditions

During amyloidogenesis, proteins undergo conformational changes that allow them to aggregate and assemble into insoluble, fibrillar structures. Soluble oligomers that form during this process typically contain 2–24 monomeric subunits and are cytotoxic. Before the formation of these soluble oligomers, monomeric species first adopt aggregation-competent conformations. Knowledge of the structures of these intermediate states is invaluable to the development of molecular strategies to arrest pathological amyloid aggregation. However, the highly dynamic and interconverting nature of amyloidogenic species limits biophysical characterization of their structures during amyloidogenesis. Here, we use molecular dynamics simulations to probe conformations sampled by monomeric transthyretin under amyloidogenic conditions. We show that certain β-strands in transthyretin tend to unfold and sample nonnative conformations and that the edge strands in one β-sheet (the DAGH sheet) are particularly susceptible to conformational changes in the monomeric state. We also find that changes in the tertiary structure of transthyretin can be associated with disruptions to the secondary structure. We evaluated the conformations produced by molecular dynamics by calculating how well molecular-dynamics-derived structures reproduced NMR-derived interatomic distances. Finally, we leverage our computational results to produce experimentally testable hypotheses that may aid experimental explorations of pathological conformations of transthyretin.     

Forcing Switch from Short- to Intermediate- and Long-lived States of the αA Domain Generates LFA-1/ICAM-1 Catch Bonds

Binding of lymphocyte function-associated antigen-1 (LFA-1) to intercellular adhesion molecule-1 (ICAM-1) mediates leukocyte adhesion under force. Using a biomembrane force probe capable of measuring single bond interactions, we showed ICAM-1 binding to LFA-1 at different conformations, including the bent conformation with the lowest affinity. We quantify how force and conformations of LFA-1 regulate its kinetics with ICAM-1. At zero-force, on-rates were substantially changed by conditions that differentially favor a bent or extended LFA-1 with a closed or open headpiece; but off-rates were identical. With increasing force, LFA-1/ICAM-1 bond lifetimes (reciprocal off-rates) first increased (catch bonds) and then decreased (slip bonds). Three states with distinct off-rates were identified from lifetime distributions. Force shifted the associated fractions from the short- to intermediate- and long-lived states, producing catch bonds at low forces, but increased their off-rates exponentially, converting catch to slip bonds at high forces. An internal ligand antagonist that blocks pulling of the α₇-helix suppressed the intermediate-/long-lived states and eliminated catch bonds, revealing an internal catch bond between the αA and βA domains. These results elucidate an allosteric mechanism for the mechanochemistry of LFA-1/ICAM-1 binding.     

Choreography for nucleosomes: the conformational freedom of the nucleosomal filament and its limitations

Eukaryotic DNA is organized into nucleosomes by coiling around core particles of histones, forming a nucleosomal filament. The significance for the conformation of the filament of the DNA entry/exit angle (α) at the nucleosome, the angle of rotation (β) of nucleosomes around their interconnecting DNA (linker DNA) and the length of the linker DNA, has been studied by means of wire models with straight linkers. It is shown that variations in α and β endow the filament with an outstanding conformational freedom when α is increased beyond 60–90°, owing to the ability of the filament to change between forward right-handed and backward left-handed coiling. A wealth of different helical and looped conformations are formed in response to repeated β sequences, and helical conformations are shown to be able to contract to a high density and to associate pairwise into different types of double fibers. Filaments with random β sequences are characterized by relatively stable loop clusters connected by segments of higher flexibility. Displacement of core particles along the DNA in such fibers, combined with limited twisting of the linkers, can generate the β sequence necessary for compaction into a regular helix, thus providing a model for heterochromatinization.     

Solution Structure of a Phytocystatin from Ananas comosus and Its Molecular Interaction with Papain

The structure of a recombinant pineapple cystatin (AcCYS) was determined by NMR with the RMSD of backbone and heavy atoms of twenty lowest energy structures of 0.56 and 1.11 Å, respectively. It reveals an unstructured N-terminal extension and a compact inhibitory domain comprising a four-stranded antiparallel β-sheet wrapped around a central α-helix. The three structural motifs (G⁴⁵, Q⁸⁹XVXG, and W¹²⁰) putatively responsible for the interaction with papain-like proteases are located in one side of AcCYS. Significant chemical shift perturbations in two loop regions, residues 45 to 48 (GIYD) and residues 89 to 91 (QVV), of AcCYS strongly suggest their involvement in the binding to papain, consistent with studies on other members of the cystatin family. However, the highly conserved W120 appears not to be involved in the binding with papain as no chemical shift perturbation was observed. Chemical shift index analysis further indicates that the length of the α-helix is shortened upon association with papain. Collectively, our data suggest that AcCYS undergoes local secondary structural rearrangements when papain is brought into close contact. A molecular model of AcCYS/papain complex is proposed to illustrate the interaction between AcCYS and papain, indicating a complete blockade of the catalytic triad by AcCYS.