Ellinor Haglund

The unique cysteine knot regulates the pleotropic hormone leptin.

Leptin plays a key role in regulating energy intake/expenditure, metabolism and hypertension. It folds into a four-helix bundle that binds to the extracellular receptor to initiate signaling. Our work on leptin revealed a hidden complexity in the formation of a previously un-described, cysteine-knotted topology in leptin. We hypothesized that this unique topology could offer new mechanisms in regulating the protein activity. A combination of in silico simulation and in vitro experiments was used to probe the role of the knotted topology introduced by the disulphide-bridge on leptin folding and function. Our results surprisingly show that the free energy landscape is conserved between knotted and unknotted protein, however the additional complexity added by the knot formation is structurally important. Native state analyses led to the discovery that the disulphide-bond plays an important role in receptor binding and thus mediate biological activity by local motions on distal receptor-binding sites, far removed from the disulphide-bridge. Thus, the disulphide-bridge appears to function as a point of tension that allows dissipation of stress at a distance in leptin.

Pierced Lasso Bundles Are a New Class of Knot-like Motifs

A four-helix bundle is a well-characterized motif often used as a target for designed pharmaceutical therapeutics and nutritional supplements. Recently, we discovered a new structural complexity within this motif created by a disulphide bridge in the long-chain helical bundle cytokine leptin. When oxidized, leptin contains a disulphide bridge creating a covalent-loop through which part of the polypeptide chain is threaded (as seen in knotted proteins). We explored whether other proteins contain a similar intriguing knot-like structure as in leptin and discovered 11 structurally homologous proteins in the PDB. We call this new helical family class the Pierced Lasso Bundle (PLB) and the knot-like threaded structural motif a Pierced Lasso (PL). In the current study, we use structure-based simulation to investigate the threading/folding mechanisms for all the PLBs along with three unthreaded homologs as the covalent loop (or lasso) in leptin is important in folding dynamics and activity. We find that the presence of a small covalent loop leads to a mechanism where structural elements slipknot to thread through the covalent loop. Larger loops use a piercing mechanism where the free terminal plugs through the covalent loop. Remarkably, the position of the loop as well as its size influences the native state dynamics, which can impact receptor binding and biological activity. This previously unrecognized complexity of knot-like proteins within the helical bundle family comprises a completely new class within the knot family, and the hidden complexity we unraveled in the PLBs is expected to be found in other protein structures outside the four-helix bundles. The insights gained here provide critical new elements for future investigation of this emerging class of proteins, where function and the energetic landscape can be controlled by hidden topology, and should be take into account in ab initio predictions of newly identified protein targets.

Complex lasso: new entangled motifs in proteins

We identify new entangled motifs in proteins that we call complex lassos. Lassos arise in proteins with disulfide bridges (or in proteins with amide linkages), when termini of a protein backbone pierce through an auxiliary surface of minimal area, spanned on a covalent loop. We find that as much as 18% of all proteins with disulfide bridges in a non-redundant subset of PDB form complex lassos, and classify them into six distinct geometric classes, one of which resembles supercoiling known from DNA. Based on biological classification of proteins we find that lassos are much more common in viruses, plants and fungi than in other kingdoms of life. We also discuss how changes in the oxidation/reduction potential may affect the function of proteins with lassos. Lassos and associated surfaces of minimal area provide new, interesting and possessing many potential applications geometric characteristics not only of proteins, but also of other biomolecules.

Engineering covalent loops in proteins can serve as an on/off switch to regulate threaded topologies.

Knots in proteins are under active investigation motivating refinements of current techniques and the development of tools to better understand the knotted topology. A strong focus is to identify new knots and expand upon the current understanding of their complex topology. Previous work has shown that the knotted topology, even in the simplest case of knots, encompasses a variety of unique challenges in folding and tying a chain. To bypass many of the in vitro experimental complications involved in working with knots, it is useful to apply methodologies to a more simplified system. The pierced lasso bundles (PLB), we discovered where a single disulphide bridge holds the threaded topology together, presents a simpler system to study knots in vitro. Having a disulphide bridge as an on/off switch between the threaded/unthreaded topology is advantageous because a covalent loop allows manipulation of the knot without directly altering affecting secondary and tertiary structure. Because disulphide bridges are commonly used in protein engineering, a pierced lasso (PL) topology can be easily introduced into a protein of interest to form a knotted topology within a given secondary structure. It is also important to take into account that if formed, disulphides can inadvertently introduce an unwanted PL. This was found upon determination of the crystal structure (PDB code 2YHG) of the recently de novo designed nucleoside hydrolase. Our detailed investigations of the PL presented here will allow researchers to look at the introduction of disulphide bridges in a larger context with respect to potential geometrical consequences on the structure and functional properties of proteins.

Heterogeneous side chain conformation highlights a network of interactions implicated in hysteresis of the knotted protein, minimal tied trefoil.

Hysteresis is a signature for a bistability in the native landscape of a protein with significant transition state barriers for the interconversion of stable species. Large global stability, as in GFP, contributes to the observation of this rare hysteretic phenomenon in folding. The signature for such behavior is non-coincidence in the unfolding and refolding transitions, despite waiting significantly longer than the time necessary for complete denaturation. Our work indicates that hysteresis in the knotted protein, the minimal tied trefoil from Thermotoga maritma (MTTTm), is mediated by a network of side chain interactions within a tightly packed core. These initially identified interactions include proline 62 from a tight β-like turn, phenylalanine 65 at the beginning of the knotting loop, and histidine 114 that initiates the threading element. It is this tightly packed region and the knotting element that we propose is disrupted with prolonged incubation in the denatured state, and is involved in the observed hysteresis. Interestingly, the disruption is not linked to backbone interactions, but rather to the packing of side chains in this critical region.

Constructing a folding model for protein S6 guided by native fluctuations deduced from NMR structures

The diversity in a set of protein nuclear magnetic resonance (NMR) structures provides an estimate of native state fluctuations that can be used to refine and enrich structure-based protein models (SBMs). Dynamics are an essential part of a proteins functional native state. The dynamics in the native state are controlled by the same funneled energy landscape that guides the entire folding process. SBMs apply the principle of minimal frustration, drawn from energy landscape theory, to construct a funneled folding landscape for a given protein using only information from the native structure. On an energy landscape smoothed by evolution towards minimal frustration, geometrical constraints, imposed by the native structure, control the folding mechanism and shape the native dynamics revealed by the model. Native-state fluctuations can alternatively be estimated directly from the diversity in the set of NMR structures for a protein. Based on this information, we identify a highly flexible loop in the ribosomal protein S6 and modify the contact map in a SBM to accommodate the inferred dynamics. By taking into account the probable native state dynamics, the experimental transition state is recovered in the model, and the correct order of folding events is restored. Our study highlights how the shared energy landscape connects folding and function by showing that a better description of the native basin improves the prediction of the folding mechanism.

Geometrical Frustration in Interleukin-33 Decouples the Dynamics of the Functional Element from the Folding Transition State Ensemble

Interleukin-33 (IL-33) is currently the focus of multiple investigations into targeting pernicious inflammatory disorders. This mediator of inflammation plays a prevalent role in chronic disorders such as asthma, rheumatoid arthritis, and progressive heart disease. In order to better understand the possible link between the folding free energy landscape and functional regions in IL-33, a combined experimental and theoretical approach was applied. IL-33 is a pseudo- symmetrical protein composed of three distinct structural elements that complicate the folding mechanism due to competition for nucleation on the dominant folding route. Trefoil 1 constitutes the majority of the binding interface with the receptor whereas Trefoils 2 and 3 provide the stable scaffold to anchor Trefoil 1. We identified that IL-33 folds with a three-state mechanism, leading to a rollover in the refolding arm of its chevron plots in strongly native conditions. In addition, there is a second slower refolding phase that exhibits the same rollover suggesting similar limitations in folding along parallel routes. Characterization of the intermediate state and the rate limiting steps required for folding suggests that the rollover is attributable to a moving transition state, shifting from a post- to pre-intermediate transition state as you move from strongly native conditions to the midpoint of the transition. On a structural level, we found that initially, all independent Trefoil units fold equally well until a QCA of 0.35 when Trefoil 1 will backtrack in order to allow Trefoils 2 and 3 to fold in the intermediate state, creating a stable scaffold for Trefoil 1 to fold onto during the final folding transition. The formation of this intermediate state and subsequent moving transition state is a result of balancing the difficulty in folding the functionally important Trefoil 1 onto the remainder of the protein. Taken together our results indicate that the functional element of the protein is geometrically frustrated, requiring the more stable elements to fold first, acting as a scaffold for docking of the functional element to allow productive folding to the native state.

Guiding the design of bacterial signaling interactions using a coevolutionary landscape

The selection of mutations that encode new interactions between bacterial two-component signaling (TCS) proteins remains a significant challenge. Recent work constructed a coevolutionary landscape where mutations can readily be selected to maintain signal transfer interactions between partner TCS proteins without introducing unwanted crosstalk. A bigger challenge is to select mutations for a TCS protein from the landscape to enhance, suppress, or have a neutral effect on its basal signal transfer with a non-partner. This study focuses on the computational selection of 12 single-point mutations to a response regulator from Bacillus subtilis and its effect on phosphotransfer with a histidine kinase from Escherichia Coli. These mutations are experimentally expressed to directly test the theoretical predictions, of which seven mutants successfully perturb phosphotransfer in the predicted manner. Furthermore, Differential Scanning Calorimetry is used to monitor any protein stability effects caused by the mutations, which could be detrimental to proper protein function.

Designing bacterial signaling interactions with coevolutionary landscapes

Selecting amino acids to design novel protein-protein interactions that facilitate catalysis is a daunting challenge. We propose that a computational coevolutionary landscape based on sequence analysis alone offers a major advantage over expensive, time-consuming brute-force approaches currently employed. Our coevolutionary landscape allows prediction of single amino acid substitutions that produce functional interactions between non-cognate, interspecies signaling partners. In addition, it can also predict mutations that maintain segregation of signaling pathways across species. Specifically, predictions of phosphotransfer activity between the Escherichia coli histidine kinase EnvZ to the non-cognate receiver Spo0F from Bacillus subtilis were compiled. Twelve mutations designed to enhance, suppress, or have a neutral effect on kinase phosphotransfer activity to a non-cognate partner were selected. We experimentally tested the ability of the kinase to relay phosphate to the respective designed Spo0F receiver proteins against the theoretical predictions. Our key finding is that the coevolutionary landscape theory, with limited structural data, can significantly reduce the search-space for successful prediction of single amino acid substitutions that modulate phosphotransfer between the two-component His-Asp relay partners in a predicted fashion. This combined approach offers significant improvements over large-scale mutations studies currently used for protein engineering and design.

Structure Based Protein Models and In Vitro Experiments to Investigate the Folding Landscape of the Human Obesity Protein Leptin

Jeffery Friedman discovered the human obesity gen in 1994. The expressed hormone leptin plays a key role in regulating energy intake/expenditure through the JAK/STAT signaling pathway. Even though, biologically, it is a well-studied protein not much is known about the biophysical properties.The mature protein is 146 residues (16kDa) buildup of 5 α-helixes and one stabilizing disulfide bond. This disulfide bridge encloses the C-terminal cysteine with residue 96 into a 50 residues long loop, creating a ‘so-called’ cysteine knot.In this research we explore the free energy landscape through structure-based protein models (Cα- and all-atom) and experiments (CD and fluorescence). The result shows that leptin can fold both in its reduced- and oxidized state. Even though the oxidized protein has a non-trivial pathway, where the N-terminal has to ‘thread’ through the loop to reach the native state, the folding landscape looks almost the same as for the reduced protein. The first event of folding is to form the loop or horseshoe (in the reduced state) constructed of α4 and α5 to stabilize the structure. Secondly, α3 is pulled into position: Either by a slipknot, backing into the loop in the oxidized state or by being pulled back to the correct conformation behind the horseshoe in the reduced state. The dragging force for this slipknot event is pulled from the energetic gain of zipping α2 together with α5. The last step is to fold α1 into its native conformation in the back of the structure between α5 and α2.