Interleukin-2 druggability is modulated by global conformational transitions controlled by a helical capping switch.

Interleukin-2 (IL-2) is a small α-helical cytokine that regulates immune cell homeostasis through its recruitment to a high-affinity heterotrimeric receptor complex (IL-2Rα/IL-2Rβ/γc). IL-2 has been shown to have therapeutic efficacy for immune diseases by preferentially expanding distinct T cell compartments, and several regulatory T cell (Treg)-biasing anti-IL-2 antibodies have been developed for combination therapies. The conformational plasticity of IL-2 plays an important role in its biological actions by modulating the strength of receptor and drug interactions. Through an NMR analysis of milliseconds-timescale dynamics of free mouse IL-2 (mIL-2), we identify a global transition to a sparse conformation which is regulated by an α-helical capping "switch" at the loop between the A and B helices (AB loop). Binding to either an anti-mouse IL-2 monoclonal antibody (mAb) or a small molecule inhibitor near the loop induces a measurable response at the core of the structure, while locking the switch to a single conformation through a designed point mutation leads to a global quenching of core dynamics accompanied by a pronounced effect in mAb binding. By elucidating key details of the long-range allosteric communication between the receptor binding surfaces and the core of the IL-2 structure, our results offer a direct blueprint for designing precision therapeutics targeting a continuum of conformational states

Peptide exchange on MHC-I by TAPBPR is driven by a negative allostery release cycle.

Chaperones TAPBPR and tapasin associate with class I major histocompatibility complexes (MHC-I) to promote optimization (editing) of peptide cargo. Here, we use solution NMR to investigate the mechanism of peptide exchange. We identify TAPBPR-induced conformational changes on conserved MHC-I molecular surfaces, consistent with our independently determined X-ray structure of the complex. Dynamics present in the empty MHC-I are stabilized by TAPBPR and become progressively dampened with increasing peptide occupancy. Incoming peptides are recognized according to the global stability of the final pMHC-I product and anneal in a native-like conformation to be edited by TAPBPR. Our results demonstrate an inverse relationship between MHC-I peptide occupancy and TAPBPR binding affinity, wherein the lifetime and structural features of transiently bound peptides control the regulation of a conformational switch located near the TAPBPR binding site, which triggers TAPBPR release. These results suggest a similar mechanism for the function of tapasin in the peptide-loading complex

Role of water in Protein Aggregation and Amyloid Polymorphism

A variety of neurodegenerative diseases are associated with the formation of amyloid plaques. Our incomplete understanding of this process underscores the need to decipher the principles governing protein aggregation. Most experimental and simulation studies have been interpreted largely from the perspective of proteins: the role of solvent has been relatively overlooked. In this Account, we provide a perspective on how interactions with water affect folding landscapes of A$\beta$ monomers, A$\beta_{16-22}$ oligomer formation, and protofilament formation in a Sup35 peptide. Simulations show that the formation of aggregation-prone structures (N$^*$) similar to the structure in the fibril requires overcoming high desolvation barrier. The mechanism of protofilament formation in a polar Sup35 peptide fragment illustrates that water dramatically slows down self-assembly. Release of water trapped in the pores as water wires creates protofilament with a dry interface. Similarly, one of the main driving force for addition of a solvated monomer to a preformed fibril is the entropy gain of released water. We conclude by postulating that two-step model for protein crystallization must also hold for higher order amyloid structure formation starting from N$^*$. Multiple N$^*$ structures with varying water content results in a number of distinct water-laden polymorphic structures. In predominantly hydrophobic sequences, water accelerates fibril formation. In contrast, water-stabilized metastable intermediates dramatically slow down fibril growth rates in hydrophilic sequences.Comment: 27 pages, 4 figures; Accounts of Chemical Research, 201

In Silico Elucidation of the Recognition Dynamics of Ubiquitin

Elucidation of the mechanism of biomacromolecular recognition events has been a topic of intense interest over the past century. The inherent dynamic nature of both protein and ligand molecules along with the continuous reshaping of the energy landscape during the binding process renders it difficult to characterize this process at atomic detail. Here, we investigate the recognition dynamics of ubiquitin via microsecond all-atom molecular dynamics simulation providing both thermodynamic and kinetic information. The high-level of consistency found with respect to experimental NMR data lends support to the accuracy of the in silico representation of the conformational substates and their interconversions of free ubiquitin. Using an energy-based reweighting approach, the statistical distribution of conformational states of ubiquitin is monitored as a function of the distance between ubiquitin and its binding partner Hrs-UIM. It is found that extensive and dense sampling of conformational space afforded by the µs MD trajectory is essential for the elucidation of the binding mechanism as is Boltzmann sampling, overcoming inherent limitations of sparsely sampled empirical ensembles. The results reveal a population redistribution mechanism that takes effect when the ligand is at intermediate range of 1–2 nm from ubiquitin. This mechanism, which may be depicted as a superposition of the conformational selection and induced fit mechanisms, also applies to other binding partners of ubiquitin, such as the GGA3 GAT domain

An siRNA Screen in Pancreatic Beta Cells Reveals a Role for Gpr27 in Insulin Production

The prevalence of type 2 diabetes in the United States is projected to double or triple by 2050. We reasoned that the genes that modulate insulin production might be new targets for diabetes therapeutics. Therefore, we developed an siRNA screening system to identify genes important for the activity of the insulin promoter in beta cells. We created a subclone of the MIN6 mouse pancreatic beta cell line that expresses destabilized GFP under the control of a 362 base pair fragment of the human insulin promoter and the mCherry red fluorescent protein under the control of the constitutively active rous sarcoma virus promoter. The ratio of the GFP to mCherry fluorescence of a cell indicates its insulin promoter activity. As G protein coupled receptors (GPCRs) have emerged as novel targets for diabetes therapies, we used this cell line to screen an siRNA library targeting all known mouse GPCRs. We identified several known GPCR regulators of insulin secretion as regulators of the insulin promoter. One of the top positive regulators was Gpr27, an orphan GPCR with no known role in beta cell function. We show that knockdown of Gpr27 reduces endogenous mouse insulin promoter activity and glucose stimulated insulin secretion. Furthermore, we show that Pdx1 is important for Gpr27's effect on the insulin promoter and insulin secretion. Finally, the over-expression of Gpr27 in 293T cells increases inositol phosphate levels, while knockdown of Gpr27 in MIN6 cells reduces inositol phosphate levels, suggesting this orphan GPCR might couple to Gq/11. In summary, we demonstrate a MIN6-based siRNA screening system that allows rapid identification of novel positive and negative regulators of the insulin promoter. Using this system, we identify Gpr27 as a positive regulator of insulin production

Dimer Formation Enhances Structural Differences between Amyloid β-Protein (1–40) and (1–42): An Explicit-Solvent Molecular Dynamics Study

Amyloid -protein (A) is central to the pathology of Alzheimer's disease. A 5% difference in the primary structure of the two predominant alloforms, A and A, results in distinct assembly pathways and toxicity properties. Discrete molecular dynamics (DMD) studies of A and A assembly resulted in alloform-specific oligomer size distributions consistent with experimental findings. Here, a large ensemble of DMD–derived A and A monomers and dimers was subjected to fully atomistic molecular dynamics (MD) simulations using the OPLS-AA force field combined with two water models, SPCE and TIP3P. The resulting all-atom conformations were slightly larger, less compact, had similar turn and lower -strand propensities than those predicted by DMD. Fully atomistic A and A monomers populated qualitatively similar free energy landscapes. In contrast, the free energy landscape of A dimers indicated a larger conformational variability in comparison to that of A dimers. A dimers were characterized by an increased flexibility in the N-terminal region D1-R5 and a larger solvent exposure of charged amino acids relative to A dimers. Of the three positively charged amino acids, R5 was the most and K16 the least involved in salt bridge formation. This result was independent of the water model, alloform, and assembly state. Overall, salt bridge propensities increased upon dimer formation. An exception was the salt bridge propensity of K28, which decreased upon formation of A dimers and was significantly lower than in A dimers. The potential relevance of the three positively charged amino acids in mediating the A oligomer toxicity is discussed in the light of available experimental data

Insight into the Stability of Cross-β Amyloid Fibril from VEALYL Short Peptide with Molecular Dynamics Simulation

Amyloid fibrils are found in many fatal neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, type II diabetes, and prion disease. The VEALYL short peptide from insulin has been confirmed to aggregate amyloid-like fibrils. However, the aggregation mechanism of amyloid fibril is poorly understood. Here, we utilized molecular dynamics simulation to analyse the stability of VEALYL hexamer. The statistical results indicate that hydrophobic residues play key roles in stabilizing VEALYL hexamer. Single point and two linkage mutants confirmed that Val1, Leu4, and Tyr5 of VEALYL are key residues. The consistency of the results for the VEALYL oligomer suggests that the intermediate states might be trimer (3-0) and pentamer(3-2). These results can help us to obtain an insight into the aggregation mechanism of amyloid fibril. These methods can be used to study the stability of amyloid fibril from other short peptides

Conformational Preferences of a 14-Residue Fibrillogenic Peptide from Acetylcholinesterase†

A 14-residue fragment from near the C-terminus of the enzyme acetylcholinesterase (AChE) is believed to have a neurotoxic/neurotrophic effect acting via an unknown pathway. While the peptide is α-helical in the full-length enzyme, the structure and association mechanism of the fragment are unknown. Using multiple molecular dynamics simulations, starting from a tetrameric complex of the association domain of AChE and systematicall disassembled subsets that include the peptide fragment, we show that the fragment is incapable of retaining its helicity in solution. Extensive replica exchange Monte Carlo folding and unfolding simulations in implicit solvent with capped and uncappted termini failed to converge to any consistent cluster of structures, suggesting that the fragment remains largely unstructured in solution under the conditions considered. Furthermore, extended molecular dynamics simulations of two steric zipper models show that the peptide is likely to form a zipper with antiparallel sheets and that peptides with mutations known to prevent fibril formation likely do so by interfering with this packing. The results demonstrate how the local environment of a peptide can stabilize a particular conformation