Biophysical characterization of the CXC chemokine receptor 2 ligands.

The chemokines of the immune system act as first responders by operating as chemoattractants, directing immune cells to specific locations of inflamed tissues. This promiscuous network is comprised of 50 ligands and 18 receptors where the ligands may interact with the receptors in various oligomeric states i.e., monomers, homodimers, and heterodimers. Chemokine receptors are G-protein coupled receptors (GPCRs) present in the membrane of immune cells. The migration of immune cells occurs in response to a concentration gradient of the ligands. Chemotaxis of neutrophils is directed by CXC-ligand (CXCL) activation of the membrane bound CXC chemokine receptor 2 (CXCR2). CXCR2 plays an important role in human health and is linked to disorders such as autoimmune disorders, inflammation, and cancer. Yet, despite their important role, little is known about the biophysical characteristics controlling ligand:ligand and ligand:receptor interaction essential for biological activity. In this work, we study the homodimers of three of the CXCR2 cognate ligands, CXCL1, CXCL5, and CXCL8. The ligands share high structural integrity but a low sequence identity. We show that the sequence diversity has evolved different binding affinities and stabilities for the CXC-ligands resulting in diverse agonist/antagonist behavior. Furthermore, CXC-ligands fold through a three-state mechanism, populating a folded monomeric state before associating into an active dimer.

Exploring the folding landscape of leptin: Insights into threading pathways.

The discovery of new protein topologies with entanglements and loop-crossings have shown the impact of local amino acid arrangement and global three-dimensional structures. This phenomenon plays a crucial role in understanding how protein structure relates to folding and function, affecting the global stability, and biological activity. Protein entanglements encompassing knots and non-trivial topologies add complexity to their folding free energy landscapes. However, the initial native contacts driving the threading event for entangled proteins remains elusive. The Pierced Lasso Topology (PLT) represents an entangled topology where a covalent linker creates a loop in which the polypeptide backbone is threaded through. Compared to true knotted topologies, PLTs are simpler topologies where the covalent-loop persists in all conformations. In this work, the PLT protein leptin, is used to visualize and differentiate the preference for slipknotting over plugging transition pathways along the folding route. We utilize the Energy Landscape Visualization Method (ELViM), a multidimensional projection technique, to visualize and distinguish early threaded conformations that cannot be observed in an in vitro experiment. Critical contacts for the leptin threading mechanisms were identified where the competing pathways are determined by the formation of a hairpin loop in the unfolded basin. Thus, prohibiting the dominant slipknotting pathway. Furthermore, ELViM offers insights into distinct folding pathways associated with slipknotting and plugging providing a novel tool for de novo design and in vitro experiments with residue specific information of threading events in silico.

Topological Reaction Coordinate Captures the Folding Transition State Ensemble in a Pierced Lasso Protein.

Proteins with a pierced lasso topology (PLT) have a covalent loop created by a disulfide bond, and the backbone circles back to thread the loop. This threaded topology has unique features compared to knotted topologies; notably, the topology is controlled by the chemical environment and the covalent loop remains intact even when denatured. In this work, we use the hormone leptin as our model PLT system and study its folding using molecular dynamics simulations that employ a structure-based (Go̅-like) model. We find that the reduced protein has a two-state folding mechanism with a transition state ensemble (TSE) that can be characterized by the reaction coordinate Q, the fraction of native contacts formed. In contrast, the oxidized protein, which must thread part of the polypeptide chain through a covalent loop, has a folding process that is poorly characterized by Q. Instead, we find that a topological coordinate that monitors the residue crossing the loop can identify the TSE of oxidized leptin. By precisely identifying the predicted TSE, one may now reliably calculate theoretical phi-values for the PLT protein, thereby enabling a comparison with experimental measurements. We find the loop-threading constraint leads to noncanonical phi-values that are uniformly small because this PLT protein has a flat energy landscape through the TSE.

A Small Contribution to a Large System: The Leptin Receptor Complex.

Obesity is a classified epidemic, increasing the risk of secondary diseases such as diabetes, inflammation, cardiovascular disease, and cancer. The pleiotropic hormone leptin is the proposed link for the gut-brain axis controlling nutritional status and energy expenditure. Research into leptin signaling provides great promise toward discovering therapeutics for obesity and its related diseases targeting leptin and its cognate leptin receptor (LEP-R). The molecular basis underlying the human leptin receptor complex assembly remains obscure, due to the lack of structural information regarding the biologically active complex. In this work, we investigate the proposed receptor binding sites in human leptin utilizing designed antagonist proteins combined with AlphaFold predictions. Our results show that binding site I has a more intricate role in the active signaling complex than previously described. We hypothesize that the hydrophobic patch in this region engages a third receptor forming a higher-order complex, or a new LEP-R binding site inducing allosteric rearrangement.

Menthol in electronic cigarettes causes biophysical inhibition of pulmonary surfactant.

With an increasing prevalence of electronic cigarette (e-cigarette) use, especially among youth, there is an urgent need to better understand the biological risks and pathophysiology of health conditions related to e-cigarettes. A majority of e-cigarette aerosols are in the submicron size and would deposit in the alveolar region of the lung, where they must first interact with the endogenous pulmonary surfactant. To date, little is known whether e-cigarette aerosols have an adverse impact on the pulmonary surfactant. We have systematically studied the effect of individual e-cigarette ingredients on an animal-derived clinical surfactant preparation, bovine lipid extract surfactant, using a combination of biophysical and analytical techniques, including in vitro biophysical simulations using constrained drop surfactometry, molecular imaging with atomic force microscopy, chemical assays using carbon nuclear magnetic resonance and circular dichroism, and in silico molecular dynamics simulations. All data collectively suggest that flavorings used in e-cigarettes, especially menthol, play a predominant role in inhibiting the biophysical function of the surfactant. The mechanism of biophysical inhibition appears to involve menthol interactions with both phospholipids and hydrophobic proteins of the natural surfactant. These results provide novel insights into the understanding of the health impact of e-cigarettes and may contribute to better regulation of e-cigarette products.

BAP1 forms a trimer with HMGB1 and HDAC1 that modulates gene × environment interaction with asbestos.

Carriers of heterozygous germline BAP1 mutations (BAP1+/-) are affected by the "BAP1 cancer syndrome." Although they can develop almost any cancer type, they are unusually susceptible to asbestos carcinogenesis and mesothelioma. Here we investigate why among all carcinogens, BAP1 mutations cooperate with asbestos. Asbestos carcinogenesis and mesothelioma have been linked to a chronic inflammatory process promoted by the extracellular release of the high-mobility group box 1 protein (HMGB1). We report that BAP1+/- cells secrete increased amounts of HMGB1, and that BAP1+/- carriers have detectable serum levels of acetylated HMGB1 that further increase when they develop mesothelioma. We linked these findings to our discovery that BAP1 forms a trimeric protein complex with HMGB1 and with histone deacetylase 1 (HDAC1) that modulates HMGB1 acetylation and its release. Reduced BAP1 levels caused increased ubiquitylation and degradation of HDAC1, leading to increased acetylation of HMGB1 and its active secretion that in turn promoted mesothelial cell transformation.

Topological Twists in Nature.

The first entangled protein was observed about 30 years ago, resulting in an increased interest for uncovering the biological functions and biophysical properties of these complex topologies. Recently, the Pierced Lasso Topology (PLT) was discovered in which a covalent bond forms an intramolecular loop, leaving one or both termini free to pierce the loop. This topology is related to knots and other entanglements. PLTs exist in many well-researched systems where the PLTs have previously been unnoticed. PLTs represents 18% of all disulfide containing proteins across all kingdoms of life. In this review, we investigate the biological implications of this specific topology in which the PLT-forming disulfide may act as a molecular switch for protein function and consequently human health.

The Pierced Lasso Topology Leptin has a Bolt on Dynamic Domain Composed by the Disordered Loops I and III.

Leptin is an important signaling hormone, mostly known for its role in energy expenditure and satiety. Furthermore, leptin plays a major role in other proteinopathies, such as cancer, marked hyperphagia, impaired immune function, and inflammation. In spite of its biological relevance in human health, there are no NMR resonance assignments of the human protein available, obscuring high-resolution characterization of the soluble protein and/or its conformational dynamics, suggested as being important for receptor interaction and biological activity. Here, we report the nearly complete backbone resonance assignments of human leptin. Chemical shift-based secondary structure prediction confirms that in solution leptin forms a four-helix bundle including a pierced lasso topology. The conformational dynamics, determined on several timescales, show that leptin is monomeric, has a rigid four-helix scaffold, and a dynamic domain, including a transiently formed helix. The dynamic domain is anchored to the helical scaffold by a secondary hydrophobic core, pinning down the long loops of leptin to the protein body, inducing motional restriction without a well-defined secondary or tertiary hydrogen bond stabilized structure. This dynamic region is well suited for and may be involved in functional allosteric dynamics upon receptor binding.

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.

Uncovering the molecular mechanisms behind disease-associated leptin variants.

The pleiotropic hormone leptin has a pivotal role in regulating energy balance by inhibiting hunger and increasing energy expenditure. Homozygous mutations found in the leptin gene are associated with extreme obesity, marked hyperphagia, and impaired immune function. Although these mutations have been characterized in vivo, a detailed understanding of how they affect leptin structure and function remains elusive. In the current work, we used NMR, differential scanning calorimetry, molecular dynamics simulations, and bioinformatics calculations to characterize the effects of these mutations on leptin structure and function and binding to its cognate receptor. We found that mutations identified in patients with congenital leptin deficiency not only cause leptin misfolding or aggregation, but also cause changes in the dynamics of leptin residues on the receptor-binding interface. Therefore, we infer that mutation-induced leptin deficiency may arise from several distinct mechanisms including (i) blockade of leptin receptor interface II, (ii) decreased affinity in the second step of leptin's interaction with its receptor, (iii) leptin destabilization, and (iv) unsuccessful threading through the covalent loop, leading to leptin misfolding/aggregation. We propose that this expanded framework for understanding the mechanisms underlying leptin deficiency arising from genetic mutations may be useful in designing therapeutics for leptin-associated disorders.

Pierced Lasso Topology Controls Function in Leptin.

Protein engineering is a powerful tool in drug design and therapeutics, where disulphide bridges are commonly introduced to stabilize proteins. However, these bonds also introduce covalent loops, which are often neglected. These loops may entrap the protein backbone on opposite sides, leading to a "knotted" topology, forming a so-called Pierced Lasso (PL). In this elegant system, the "knot" is held together with a single disulphide bridge where part of the polypeptide chain is threaded through. The size and position of these covalent loops can be manipulated through protein design in vitro, whereas nature uses polymorphism to switch the PL topology. The PL protein leptin shows genetic modification of an N-terminal residue, adding a third cysteine to the same sequence. In an effort to understand the mechanism of threading of these diverse topologies, we designed three loop variants to mimic the polymorphic sequence. This adds elegance to the system under study, as it allows the generation of three possible covalent loops; they are the original wild-type C-terminal loop protein, the fully circularized unthreaded protein, and the N-terminal loop protein, responsible for different lasso topologies. The size of the loop changes the threading mechanism from a slipknotting to a plugging mechanism, with increasing loop size. Interestingly, the ground state of the native protein structure is largely unaffected, but biological assays show that the activity is maximized by properly controlled dynamics in the threaded state. A threaded topology with proper conformational dynamics is important for receptor interaction and activation of the signaling pathways in vivo.

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.

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.

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 protein's 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.

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.

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.

Trimming down a protein structure to its bare foldons: spatial organization of the cooperative unit.

Folding of the ribosomal protein S6 is a malleable process controlled by two competing, and partly overlapping, folding nuclei. Together, these nuclei extend over most of the S6 structure, except the edge strand ß2, which is consistently missing in the folding transition states; despite being part of the S6 four-stranded sheet, ß2 seems not to be part of the cooperative unit of the protein. The question is then whether ß2 can be removed from the S6 structure without compromising folding cooperativity or native state integrity. To investigate this, we constructed a truncated variant of S6 lacking ß2, reducing the size of the protein from 96 to 76 residues (S6(Δß2)). The new S6 variant expresses well in Escherichia coli and has a well dispersed heteronuclear single quantum correlation spectrum and a perfectly wild-type-like crystal structure, but with a smaller three-stranded ß-sheet. Moreover, S6(Δß2) displays an archetypical v-shaped chevron plot with decreased slope of the unfolding limb, as expected from a protein with maintained folding cooperativity and reduced size. The results support the notion that foldons, as defined by the structural distribution of the folding nuclei, represent a property-based level of hierarchy in the build-up of larger protein structures and suggest that the role of ß2 in S6 is mainly in intermolecular binding, consistent with the position of this strand in the ribosomal assembly.

The HD-exchange motions of ribosomal protein S6 are insensitive to reversal of the protein-folding pathway.

An increasing number of protein structures are found to encompass multiple folding nuclei, allowing their structures to be formed by several competing pathways. A typical example is the ribosomal protein S6, which comprises two folding nuclei (sigma1 and sigma2) defining two competing pathways in the folding energy landscape: sigma1 --> sigma2 and sigma2 --> sigma1. The balance between the two pathways, and thus the order of folding events, is easily controlled by circular permutation. In this study, we make use of this ability to manipulate the folding pathway to demonstrate that the dynamic motions of the S6 structure are independent of how the protein folds. The HD-exchange protection factors remain the same upon complete reversal of the folding order. The phenomenon arises because the HD-exchange motions and the high-energy excitations controlling the folding pathway occur at separated free-energy levels: the Boltzmann distribution of unproductive unfolding attempts samples all unfolding channels in parallel, even those that end up in excessively high barriers. Accordingly, the findings provide a simple rationale for how to interpret native-state dynamics without the need to invoke fluctuations off the normal unfolding reaction coordinate.