Search

The de novo design of small molecule-binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase-1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from

The SARS-CoV-2 envelope (E) protein forms a five-helix bundle in lipid bilayers whose cation-conducting activity is associated with the inflammatory response and respiratory distress symptoms of COVID-19. E channel activity is inhibited by the drug 5-(N,N-hexamethylene) amiloride (HMA). However, the binding site of HMA in E has not been determined. Here we use solid-state NMR to measure distances between HMA and the E transmembrane domain (ETM) in lipid bilayers. 13 C, 15 N-labeled HMA is combined with fluorinated or 13 C-labeled ETM. Conversely, fluorinated HMA is combined with 13 C, 15 N-labeled ETM. These orthogonal isotopic labeling patterns allow us to conduct dipolar recoupling NMR experiments to determine the HMA binding stoichiometry to ETM as well as HMA-protein distances. We find that HMA binds ETM with a stoichiometry of one drug per pentamer. Unexpectedly, the bound HMA is not centrally located within the channel pore, but lies on the lipid-facing surface in the middle of the TM domain. This result suggests that HMA may inhibit the E channel activity by interfering with the gating function of an aromatic network. These distance data are obtained under much lower drug concentrations than in previous chemical shift perturbation data, which showed the largest perturbation for N-terminal residues. This difference suggests that HMA has higher affinity for the protein-lipid interface than the channel pore. These results give insight into the inhibition mechanism of HMA for SARS-CoV-2 E.

Fungal infections cause more than 1.5 million deaths a year. Due to emerging antifungal drug resistance, novel strategies are urgently needed to combat life-threatening fungal diseases. Here, we identify the host defense peptide mimetic, brilacidin (BRI) as a synergizer with caspofungin (CAS) against CAS-sensitive and CAS-resistant isolates of Aspergillus fumigatus, Candida albicans, C. auris, and CAS-intrinsically resistant Cryptococcus neoformans. BRI also potentiates azoles against A. fumigatus and several Mucorales fungi. BRI acts in A. fumigatus by affecting cell wall integrity pathway and cell membrane potential. BRI combined with CAS significantly clears A. fumigatus lung infection in an immunosuppressed murine model of invasive pulmonary aspergillosis. BRI alone also decreases A. fumigatus fungal burden and ablates disease development in a murine model of fungal keratitis. Our results indicate that combinations of BRI and antifungal drugs in clinical use are likely to improve the treatment outcome of aspergillosis and other fungal infections.

In neurodegenerative diseases, proteins fold into amyloid structures with distinct conformations (strains) that are characteristic of different diseases. However, there is a need to rapidly identify amyloid conformations in situ. Here, we use machine learning on the full information available in fluorescent excitation/emission spectra of amyloid-binding dyes to identify six distinct different conformational strains in vitro, as well as amyloid-β (Aβ) deposits in different transgenic mouse models. Our EMBER (excitation multiplexed bright emission recording) imaging method rapidly identifies conformational differences in Aβ and tau deposits from Down syndrome, sporadic and familial Alzheimer's disease human brain slices. EMBER has in situ identified distinct conformational strains of tau inclusions in astrocytes, oligodendrocytes, and neurons from Pick's disease. In future studies, EMBER should enable high-throughput measurements of the fidelity of strain transmission in cellular and animal neurodegenerative diseases models, time course of amyloid strain propagation, and identification of pathogenic versus benign strains.

Accumulation of filamentous aggregates of tau protein in the brain is a pathological hallmark of Alzheimer's disease (AD) and many other neurodegenerative tauopathies. The filaments adopt disease-specific cross-β amyloid conformations that self-propagate and are implicated in neuronal loss. Development of molecular diagnostics and therapeutics is of critical importance. However, mechanisms of small molecule binding to the amyloid core is poorly understood. We used cryo-electron microscopy to determine a 2.7 Å structure of AD patient-derived tau paired-helical filaments bound to the PET ligand GTP-1. The compound is bound stoichiometrically at a single site along an exposed cleft of each protofilament in a stacked arrangement matching the fibril symmetry. Multiscale modeling reveals pi-pi aromatic interactions that pair favorably with the small molecule-protein contacts, supporting high specificity and affinity for the AD tau conformation. This binding mode offers critical insight into designing compounds to target different amyloid folds found across neurodegenerative diseases.

The SARS-CoV-2 envelope (E) protein is a viroporin associated with the acute respiratory symptoms of COVID-19. E forms cation-selective ion channels that assemble in the lipid membrane of the endoplasmic reticulum Golgi intermediate compartment. The channel activity of E is linked to the inflammatory response of the host cell to the virus. Like many viroporins, E is thought to oligomerize with a well-defined stoichiometry. However, attempts to determine the E stoichiometry have led to inconclusive results and suggested mixtures of oligomers whose exact nature might vary with the detergent used. Here, we employ 19F solid-state nuclear magnetic resonance and the centerband-only detection of exchange (CODEX) technique to determine the oligomeric number of E's transmembrane domain (ETM) in lipid bilayers. The CODEX equilibrium value, which corresponds to the inverse of the oligomeric number, indicates that ETM assembles into pentamers in lipid bilayers, without any detectable fraction of low-molecular-weight oligomers. Unexpectedly, at high peptide concentrations and in the presence of the lipid phosphatidylinositol, the CODEX data indicate that more than five 19F spins are within a detectable distance of about 2 nm, suggesting that the ETM pentamers cluster in the lipid bilayer. Monte Carlo simulations that take into account peptide-peptide and peptide-lipid interactions yielded pentamer clusters that reproduced the CODEX data. This supramolecular organization is likely important for E-mediated virus assembly and budding and for the channel function of the protein.

Chemical probes for chromatin reader proteins are valuable tools for investigating epigenetic regulatory mechanisms and evaluating whether the target of interest holds therapeutic potential. Developing potent inhibitors for the plant homeodomain (PHD) family of methylation readers remains a difficult task due to the charged, shallow and extended nature of the histone binding site that precludes effective engagement of conventional small molecules. Herein, we describe the development of novel proximity-reactive cyclopeptide inhibitors for PHD3-a trimethyllysine reader domain of histone demethylase KDM5A. Guided by the PHD3-histone co-crystal structure, we designed a sidechain-to-sidechain linking strategy to improve peptide proteolytic stability whilst maintaining binding affinity. We have developed an operationally simple solid-phase macrocyclization pathway, capitalizing on the inherent reactivity of the dimethyllysine ε-amino group to generate scaffolds bearing charged tetraalkylammonium functionalities that effectively engage the shallow aromatic 'groove' of PHD3. Leveraging a surface-exposed lysine residue on PHD3 adjacent to the ligand binding site, cyclic peptides were rendered covalent through installation of an arylsulfonyl fluoride warhead. The resulting lysine-reactive cyclic peptides demonstrated rapid and efficient labeling of the PHD3 domain in HEK293T lysates, showcasing the feasibility of employing proximity-induced reactivity for covalent labeling of this challenging family of reader domains.

Every voltage-gated ion channel (VGIC) has a pore domain (PD) made from four subunits, each comprising an antiparallel transmembrane helix pair bridged by a loop. The extent to which PD subunit structure requires quaternary interactions is unclear. Here, we present crystal structures of a set of bacterial voltage-gated sodium channel (BacNaV) 'pore only' proteins that reveal a surprising collection of non-canonical quaternary arrangements in which the PD tertiary structure is maintained. This context-independent structural robustness, supported by molecular dynamics simulations, indicates that VGIC-PD tertiary structure is independent of quaternary interactions. This fold occurs throughout the VGIC superfamily and in diverse transmembrane and soluble proteins. Strikingly, characterization of PD subunit-binding Fabs indicates that non-canonical quaternary PD conformations can occur in full-length VGICs. Together, our data demonstrate that the VGIC-PD is an autonomously folded unit. This property has implications for VGIC biogenesis, understanding functional states, de novo channel design, and VGIC structural origins.

Natural killer (NK) cells can detect antibody-coated cells through recognition by the CD16 Fc receptor. The importance of CD16 in human NK cell biology has long been appreciated, but how CD16 functions in mouse NK cells remains poorly understood. Here, we report drastic differences between human and mouse CD16 functions in NK cells. We demonstrate that one of the adaptor molecules that CD16 associates with and signals through, CD3ζ, plays a critical role in these functional differences. Using a systematic approach, we demonstrate that residues in the transmembrane domain of the mouse CD3ζ molecule prevent efficient complex formation with mouse CD16, thereby dampening receptor function. Mutating these residues in mouse CD3ζ to those encoded by human CD3ζ resulted in rescue of CD16 receptor function. We reveal that the mouse CD3ζ transmembrane domain adopts a tightly packed confirmation, preventing association with CD16, whereas human CD3ζ adopts a versatile configuration that accommodates receptor assembly.

Abstract: Natural killer (NK) cells are innate lymphocytes capable of distinguishing between healthy and pathogenic cells using a myriad of receptors expressed on their cell surface. The CD16 receptor detects the Fc portion of IgG on antibody-coated cells resulting in NK cell activation. The importance of this receptor on human NK cell biology has long been appreciated; however, how CD16 functions in mouse NK cells remains poorly understood. Here, we report drastic differences between human and mouse CD16 functions in NK cells. We demonstrate that one of the adaptor molecules that CD16 associates with and signals through, CD3ζ, plays a critical role in these functional differences. Using a systematic approach, we demonstrate that key residues in the transmembrane domain of the mouse CD3ζ molecule prevent efficient complex formation with mouse CD16, thereby dampening receptor function. Mutating these residues in mouse CD3ζ to those encoded by human CD3ζ resulted in rescue of CD16 receptor function. Supported by grants from the Cancer Research Institute (CRI), Burroughs Wellcome Fund (BWF), the Parker Institute for Cancer Immunotherapy (PICI) and NIH grants AI068129 and AI146581.

Triggering receptor expressed on myeloid cells 2 (TREM2) is a single-pass transmembrane receptor of the immunoglobulin superfamily that is secreted in a soluble (sTREM2) form. Mutations in TREM2 have been linked to increased risk of Alzheimer's disease (AD). A prominent neuropathological component of AD is deposition of the amyloid-β (Aβ) into plaques, particularly Aβ40 and Aβ42. While the membrane-bound form of TREM2 is known to facilitate uptake of Aβ fibrils and the polarization of microglial processes toward amyloid plaques, the role of its soluble ectodomain, particularly in interactions with monomeric or fibrillar Aβ, has been less clear. Our results demonstrate that sTREM2 does not bind to monomeric Aβ40 and Aβ42, even at a high micromolar concentration, while it does bind to fibrillar Aβ42 and Aβ40 with equal affinities (2.6 ± 0.3 µM and 2.3 ± 0.4 µM). Kinetic analysis shows that sTREM2 inhibits the secondary nucleation step in the fibrillization of Aβ, while having little effect on the primary nucleation pathway. Furthermore, binding of sTREM2 to fibrils markedly enhanced uptake of fibrils into human microglial and neuroglioma derived cell lines. The disease-associated sTREM2 mutant, R47H, displayed little to no effect on fibril nucleation and binding, but it decreased uptake and functional responses markedly. We also probed the structure of the WT sTREM2-Aβ fibril complex using integrative molecular modeling based primarily on the cross-linking mass spectrometry data. The model shows that sTREM2 binds fibrils along one face of the structure, leaving a second, mutation-sensitive site free to mediate cellular binding and uptake.

Amyloid beta (Aβ) is thought to play a critical role in the pathogenesis of Alzheimer's disease (AD). Prion-like Aβ polymorphs, or "strains", can have varying pathogenicity and may underlie the phenotypic heterogeneity of the disease. In order to develop effective AD therapies, it is critical to identify the strains of Aβ that might arise prior to the onset of clinical symptoms and understand how they may change with progressing disease. Down syndrome (DS), as the most common genetic cause of AD, presents promising opportunities to compare such features between early and advanced AD. In this work, we evaluate the neuropathology and Aβ strain profile in the post-mortem brain tissues of 210 DS, AD, and control individuals. We assayed the levels of various Aβ and tau species and used conformation-sensitive fluorescent probes to detect differences in Aβ strains among individuals and populations. We found that these cohorts have some common but also some distinct strains from one another, with the most heterogeneous populations of Aβ emerging in subjects with high levels of AD pathology. The emergence of distinct strains in DS at these later stages of disease suggests that the confluence of aging, pathology, and other DS-linked factors may favor conditions that generate strains that are unique from sporadic AD.

Multicomponent immune receptors are essential complexes in which distinct ligand-recognition and signaling subunits are held together by interactions between acidic and basic residues of their transmembrane helices. A 2:1 acidic-to-basic motif in the transmembrane domains of the subunits is necessary and sufficient to assemble these receptor complexes. Here, we study a prototype for these receptors, a DAP12-NKG2C 2:1 heterotrimeric complex, in which the two DAP12 subunits each contribute a single transmembrane Asp residue, and the NKG2C subunit contributes a Lys to form the complex. DAP12 can also associate with 20 other subunits using a similar motif. Here, we use molecular-dynamics simulations to understand the basis for the high affinity and diversity of interactions in this group of receptors. Simulations of the transmembrane helices with differing protonation states of the Asp-Asp-Lys triad identified a structurally stable interaction in which a singly-protonated Asp-Asp pair forms a hydrogen-bonded carboxyl-carboxylate clamp that clasps onto a charged Lys side chain. This polar motif was also supported by density functional theory and a Protein Data Bank-wide search. In contrast, the helices are dynamic at sites distal to the stable carboxyl-carboxylate clamp motif. Such a locally stable but globally dynamic structure is well suited to accommodate the sequence and structural variations in the transmembrane helices of multicomponent receptors, which mix and match subunits to create combinatorial functional diversity from a limited number of subunits. It also supports a signaling mechanism based on multisubunit clustering rather than propagation of rigid conformational changes through the membrane.

The randomization and screening of combinatorial DNA libraries is a powerful technique for understanding sequence-function relationships and optimizing biosynthetic pathways. Although it can be difficult to predict a priori which sequence combinations encode functional units, it is often possible to omit undesired combinations that inflate library size and screening effort. However, defined library generation is difficult when a complex scan through sequence space is needed. To overcome this challenge, we designed a hybrid valve- and droplet-based microfluidic system that deterministically assembles DNA parts in picoliter droplets, reducing reagent consumption and bias. Using this system, we built a combinatorial library encoding an engineered histidine kinase (HK) based on bacterial CpxA. Our library encodes designed transmembrane (TM) domains that modulate the activity of the cytoplasmic domain of CpxA and variants of the structurally distant "S helix" located near the catalytic domain. We find that the S helix sets a basal activity further modulated by the TM domain. Surprisingly, we also find that a given TM motif can elicit opposing effects on the catalytic activity of different S-helix variants. We conclude that the intervening HAMP domain passively transmits signals and shapes the signaling response depending on subtle changes in neighboring domains. This flexibility engenders a richness in functional outputs as HKs vary in response to changing evolutionary pressures.

We describe the de novo design of an allosterically regulated protein, which comprises two tightly coupled domains. One domain is based on the DF (Due Ferri in Italian or two-iron in English) family of de novo proteins, which have a diiron cofactor that catalyzes a phenol oxidase reaction, while the second domain is based on PS1 (Porphyrin-binding Sequence), which binds a synthetic Zn-porphyrin (ZnP). The binding of ZnP to the original PS1 protein induces changes in structure and dynamics, which we expected to influence the catalytic rate of a fused DF domain when appropriately coupled. Both DF and PS1 are four-helix bundles, but they have distinct bundle architectures. To achieve tight coupling between the domains, they were connected by four helical linkers using a computational method to discover the most designable connections capable of spanning the two architectures. The resulting protein, DFP1 (Due Ferri Porphyrin), bound the two cofactors in the expected manner. The crystal structure of fully reconstituted DFP1 was also in excellent agreement with the design, and it showed the ZnP cofactor bound over 12 Å from the dimetal center. Next, a substrate-binding cleft leading to the diiron center was introduced into DFP1. The resulting protein acts as an allosterically modulated phenol oxidase. Its Michaelis-Menten parameters were strongly affected by the binding of ZnP, resulting in a fourfold tighter K m and a 7-fold decrease in k cat These studies establish the feasibility of designing allosterically regulated catalytic proteins, entirely from scratch.

Here we present a platform for discovery of protease-activated prodrugs and apply it to antibiotics that target Gram-negative bacteria. Because cleavable linkers for prodrugs had not been developed for bacterial proteases, we used substrate phage to discover substrates for proteases found in the bacterial periplasm. Rather than focusing on a single protease, we used a periplasmic extract of E. coli to find sequences with the greatest susceptibility to the endogenous mixture of periplasmic proteases. Using a fluorescence assay, candidate sequences were evaluated to identify substrates that release native amine-containing payloads. We next designed conjugates consisting of (1) an N-terminal siderophore to facilitate uptake, (2) a protease-cleavable linker, and (3) an amine-containing antibiotic. Using this strategy, we converted daptomycin-which by itself is active only against Gram-positive bacteria-into an antibiotic capable of targeting Gram-negative Acinetobacter species. We similarly demonstrated siderophore-facilitated delivery of oxazolidinone and macrolide antibiotics into a number of Gram-negative species. These results illustrate this platform's utility for development of protease-activated prodrugs, including Trojan horse antibiotics.

Prevalent resistance to inhibitors that target the influenza A M2 proton channel has necessitated a continued drug design effort, supported by a sustained study of the mechanism of channel function and inhibition. Recent high-resolution X-ray crystal structures present the first opportunity to see how the adamantyl amine class of inhibitors bind to M2 and disrupt and interact with the channel's water network, providing insight into the critical properties that enable their effective inhibition in wild-type M2. In this work, we examine the hypothesis that these drugs act primarily as mechanism-based inhibitors by comparing hydrated excess proton stabilization during proton transport in M2 with the interactions revealed in the crystal structures, using the Multiscale Reactive Molecular Dynamics (MS-RMD) methodology. MS-RMD, unlike classical molecular dynamics, models the hydrated proton (hydronium-like cation) as a dynamic excess charge defect and allows bonds to break and form, capturing the intricate interactions between the hydrated excess proton, protein atoms, and water. Through this, we show that the ammonium group of the inhibitors is effectively positioned to take advantage of the channel's natural ability to stabilize an excess protonic charge and act as a hydronium mimic. Additionally, we show that the channel is especially stable in the drug binding region, highlighting the importance of this property for binding the adamantane group. Finally, we characterize an additional hinge point near Val27, which dynamically responds to charge and inhibitor binding. Altogether, this work further illuminates a dynamic understanding of the mechanism of drug inhibition in M2, grounded in the fundamental properties that enable the channel to transport and stabilize excess protons, with critical implications for future drug design efforts.

The de novo design of proteins that bind highly functionalized small molecules represents a great challenge. To enable computational design of binders, we developed a unit of protein structure-a van der Mer (vdM)-that maps the backbone of each amino acid to statistically preferred positions of interacting chemical groups. Using vdMs, we designed six de novo proteins to bind the drug apixaban; two bound with low and submicromolar affinity. X-ray crystallography and mutagenesis confirmed a structure with a precisely designed cavity that forms favorable interactions in the drug-protein complex. vdMs may enable design of functional proteins for applications in sensing, medicine, and catalysis.

Protein conformations are shaped by cellular environments, but how environmental changes alter the conformational landscapes of specific proteins in vivo remains largely uncharacterized, in part due to the challenge of probing protein structures in living cells. Here, we use deep mutational scanning to investigate how a toxic conformation of α-synuclein, a dynamic protein linked to Parkinson's disease, responds to perturbations of cellular proteostasis. In the context of a course for graduate students in the UCSF Integrative Program in Quantitative Biology, we screened a comprehensive library of α-synuclein missense mutants in yeast cells treated with a variety of small molecules that perturb cellular processes linked to α-synuclein biology and pathobiology. We found that the conformation of α-synuclein previously shown to drive yeast toxicity-an extended, membrane-bound helix-is largely unaffected by these chemical perturbations, underscoring the importance of this conformational state as a driver of cellular toxicity. On the other hand, the chemical perturbations have a significant effect on the ability of mutations to suppress α-synuclein toxicity. Moreover, we find that sequence determinants of α-synuclein toxicity are well described by a simple structural model of the membrane-bound helix. This model predicts that α-synuclein penetrates the membrane to constant depth across its length but that membrane affinity decreases toward the C terminus, which is consistent with orthogonal biophysical measurements. Finally, we discuss how parallelized chemical genetics experiments can provide a robust framework for inquiry-based graduate coursework.

Defining the biologically active structures of proteins in their cellular environments remains challenging for proteins with multiple conformations and functions, where only a minor conformer might be associated with a given function. Here, we use deep mutational scanning to probe the structure and dynamics of α-synuclein, a protein known to adopt disordered, helical and amyloid conformations. We examined the effects of 2,600 single-residue substitutions on the ability of intracellularly expressed α-synuclein to slow the growth of yeast. Computational analysis of the data showed that the conformation responsible for this phenotype is a long, uninterrupted, amphiphilic helix with increasing dynamics toward the C terminus. Deep mutational scanning can therefore determine biologically active conformations in cellular environments, even for a highly dynamic multi-conformational protein.

The V27A mutation confers adamantane resistance on the influenza A matrix 2 (M2) proton channel and is becoming more prevalent in circulating populations of influenza A virus. We have used X-ray crystallography to determine structures of a spiro-adamantyl amine inhibitor bound to M2(22-46) V27A and also to M2(21-61) V27A in the Inwardclosed conformation. The spiro-adamantyl amine binding site is nearly identical for the two crystal structures. Compared to the M2 "wild type" (WT) with valine at position 27, we observe that the channel pore is wider at its N-terminus as a result of the V27A mutation and that this removes V27 side chain hydrophobic interactions that are important for binding of amantadine and rimantadine. The spiro-adamantyl amine inhibitor blocks proton conductance in the WT and V27A mutant channels by shifting its binding site in the pore depending on which residue is present at position 27. Additionally, in the structure of the M2(21-61) V27A construct, the C-terminus of the channel is tightly packed relative to that of the M2(22-46) construct. We observe that residues Asp44, Arg45, and Phe48 face the center of the channel pore and would be well-positioned to interact with protons exiting the M2 channel after passing through the His37 gate. A 300 ns molecular dynamics simulation of the M2(22-46) V27A-spiro-adamantyl amine complex predicts with accuracy the position of the ligands and waters inside the pore in the X-ray crystal structure of the M2(22-46) V27A complex.

Proteins are molecular machines whose function depends on their ability to achieve complex folds with precisely defined structural and dynamic properties. The rational design of proteins from first-principles, or de novo, was once considered to be impossible, but today proteins with a variety of folds and functions have been realized. We review the evolution of the field from its earliest days, placing particular emphasis on how this endeavor has illuminated our understanding of the principles underlying the folding and function of natural proteins, and is informing the design of macromolecules with unprecedented structures and properties. An initial set of milestones in de novo protein design focused on the construction of sequences that folded in water and membranes to adopt folded conformations. The first proteins were designed from first-principles using very simple physical models. As computers became more powerful, the use of the rotamer approximation allowed one to discover amino acid sequences that stabilize the desired fold. As the crystallographic database of protein structures expanded in subsequent years, it became possible to construct proteins by assembling short backbone fragments that frequently recur in Nature. The second set of milestones in de novo design involves the discovery of complex functions. Proteins have been designed to bind a variety of metals, porphyrins, and other cofactors. The design of proteins that catalyze hydrolysis and oxygen-dependent reactions has progressed significantly. However, de novo design of catalysts for energetically demanding reactions, or even proteins that bind with high affinity and specificity to highly functionalized complex polar molecules remains an importnant challenge that is now being achieved. Finally, the protein design contributed significantly to our understanding of membrane protein folding and transport of ions across membranes. The area of membrane protein design, or more generally of biomimetic polymers that function in mixed or non-aqueous environments, is now becoming increasingly possible.

Misfolding of the microtubule-binding protein tau into filamentous aggregates is characteristic of many neurodegenerative diseases such as Alzheimer's disease and progressive supranuclear palsy. Determining the structures and dynamics of these tau fibrils is important for designing inhibitors against tau aggregation. Tau fibrils obtained from patient brains have been found by cryo-electron microscopy to adopt disease-specific molecular conformations. However, in vitro heparin-fibrillized 2N4R tau, which contains all four microtubule-binding repeats (4R), was recently found to adopt polymorphic structures. Here we use solid-state NMR spectroscopy to investigate the global fold and dynamics of heparin-fibrillized 0N4R tau. A single set of 13C and 15N chemical shifts was observed for residues in the four repeats, indicating a single β-sheet conformation for the fibril core. This rigid core spans the R2 and R3 repeats and adopts a hairpin-like fold that has similarities to but also clear differences from any of the polymorphic 2N4R folds. Obtaining a homogeneous fibril sample required careful purification of the protein and removal of any proteolytic fragments. A variety of experiments and polarization transfer from water and mobile side chains indicate that 0N4R tau fibrils exhibit heterogeneous dynamics: Outside the rigid R2-R3 core, the R1 and R4 repeats are semirigid even though they exhibit β-strand character and the proline-rich domains undergo large-amplitude anisotropic motions, whereas the two termini are nearly isotropically flexible. These results have significant implications for the structure and dynamics of 4R tau fibrils in vivo.

Integrin αIIbβ3, a transmembrane heterodimer, mediates platelet aggregation when it switches from an inactive to an active ligand-binding conformation following platelet stimulation. Central to regulating αIIbβ3 activity is the interaction between the αIIb and β3 extracellular stalks, which form a tight heterodimer in the inactive state and dissociate in the active state. Here, we demonstrate that alanine replacements of sensitive positions in the heterodimer stalk interface destabilize the inactive conformation sufficiently to cause constitutive αIIbβ3 activation. To determine the structural basis for this effect, we performed a structural bioinformatics analysis and found that perturbing intersubunit contacts with favorable interaction geometry through substitutions to alanine quantitatively accounted for the degree of constitutive αIIbβ3 activation. This mutational study directly assesses the relationship between favorable interaction geometry at mutation-sensitive positions and the functional activity of those mutants, giving rise to a simple model that highlights the importance of interaction geometry in contributing to the stability between protein-protein interactions.

The amantadine-resistant S31N mutant of the influenza A M2 proton channel has become prevalent in currently circulating viruses. Here, we have solved an X-ray crystal structure of M2(22-46) S31N that contains two distinct conformational states within its asymmetric unit. This structure reveals the mechanism of adamantane resistance in both conformational states of the M2 channel. In the Inwardopen conformation, the mutant Asn31 side chain faces the channel pore and sterically blocks the adamantane binding site. In the Inwardclosed conformation, Asn31 forms hydrogen bonds with carbonyls at the monomer-monomer interface, which twists the monomer helices and constricts the channel pore at the drug binding site. We also examine M2(19-49) WT and S31N using solution NMR spectroscopy and show that distribution of the two conformational states is dependent on both detergent choice and experimental pH.

The influenza A M2 protein is an acid-activated proton channel responsible for acidification of the inside of the virus, a critical step in the viral life cycle. This channel has four central histidine residues that form an acid-activated gate, binding protons from the outside until an activated state allows proton transport to the inside. While previous work has focused on proton transport through the channel, the structural and dynamic changes that accompany proton flux and enable activation have yet to be resolved. In this study, extensive Multiscale Reactive Molecular Dynamics simulations with explicit Grotthuss-shuttling hydrated excess protons are used to explore detailed molecular-level interactions that accompany proton transport in the +0, + 1, and +2 histidine charge states. The results demonstrate how the hydrated excess proton strongly influences both the protein and water hydrogen-bonding network throughout the channel, providing further insight into the channel's acid-activation mechanism and rectification behavior. We find that the excess proton dynamically, as a function of location, shifts the protein structure away from its equilibrium distributions uniquely for different pH conditions consistent with acid-activation. The proton distribution in the xy-plane is also shown to be asymmetric about the channel's main axis, which has potentially important implications for the mechanism of proton conduction and future drug design efforts.

Lateral transmembrane (TM) helix-helix interactions between single-span membrane proteins play an important role in the assembly and signaling of many cell-surface receptors. Often, these helices contain two highly conserved yet distinct interaction motifs, arranged such that the motifs cannot be engaged simultaneously. However, there is sparse experimental evidence that dual-engagement mechanisms play a role in biological signaling. Here, we investigate the function of the two conserved interaction motifs in the TM domain of the integrin β3-subunit. The first motif uses reciprocating "large-large-small" amino acid packing to mediate the interaction of the β3 and αIIb TM domains and maintain the inactive resting conformation of the platelet integrin αIIbβ3. The second motif, S-x3-A-x3-I, is a variant of the classical "G-x3-G" motif. Using site-directed mutagenesis, optical trap-based force spectroscopy, and molecular modeling, we show that S-x3-A-x3-I does not engage αIIb but rather mediates the interaction of the β3 TM domain with the TM domain of the αv-subunit of the integrin αvβ3. Like αIIbβ3, αvβ3 on circulating platelets is inactive, and in the absence of platelet stimulation is unable to interact with components of the subendothelial matrix. However, disrupting any residue in the β3 S-x3-A-x3-I motif by site-directed mutations is sufficient to induce αvβ3 binding to the αvβ3 ligand osteopontin and to the monoclonal antibody WOW-1. Thus, the β3-integrin TM domain is able to engage in two mutually exclusive interactions that produce alternate α-subunit pairing, creating two integrins with distinct biological functions.

De novo protein design represents an attractive approach for testing and extending our understanding of metalloprotein structure and function. Here, we describe our work on the design of DF (Due Ferri or two-iron in Italian), a minimalist model for the active sites of much larger and more complex natural diiron and dimanganese proteins. In nature, diiron and dimanganese proteins protypically bind their ions in 4-Glu, 2-His environments, and they catalyze diverse reactions, ranging from hydrolysis, to O2-dependent chemistry, to decarbonylation of aldehydes. In the design of DF, the position of each atom-including the backbone, the first-shell ligands, the second-shell hydrogen-bonded groups, and the well-packed hydrophobic core-was bespoke using precise mathematical equations and chemical principles. The first member of the DF family was designed to be of minimal size and complexity and yet to display the quintessential elements required for binding the dimetal cofactor. After thoroughly characterizing its structural, dynamic, spectroscopic, and functional properties, we added additional complexity in a rational stepwise manner to achieve increasingly sophisticated catalytic functions, ultimately demonstrating substrate-gated four-electron reduction of O2 to water. We also briefly describe the extension of these studies to the design of proteins that bind nonbiological metal cofactors (a synthetic porphyrin and a tetranuclear cluster), and a Zn2+/proton antiporting membrane protein. Together these studies demonstrate a successful and generally applicable strategy for de novo metalloprotein design, which might indeed mimic the process by which primordial metalloproteins evolved. We began the design process with a highly symmetrical backbone and binding site, by using point-group symmetry to assemble the secondary structures that position the amino acid side chains required for binding. The resulting models provided a rough starting point and initial parameters for the subsequent precise design of the final protein using modern methods of computational protein design. Unless the desired site is itself symmetrical, this process requires reduction of the symmetry or lifting it altogether. Nevertheless, the initial symmetrical structure can be helpful to restrain the search space during assembly of the backbone. Finally, the methods described here should be generally applicable to the design of highly stable and robust catalysts and sensors. There is considerable potential in combining the efficiency and knowledge base associated with homogeneous metal catalysis with the programmability, biocompatibility, and versatility of proteins. While the work reported here focuses on testing and learning the principles of natural metalloproteins by designing and studying proteins one at a time, there is also considerable potential for using designed proteins that incorporate both biological and nonbiological metal ion cofactors for the evolution of novel catalysts.

Expansion of the genetic code with unnatural amino acids (Uaas) has significantly increased the chemical space available to proteins for exploitation. Due to the inherent limitation of translational machinery and the required compatibility with biological settings, function groups introduced via Uaas to date are restricted to chemically inert, bioorthogonal, or latent bioreactive groups. To break this barrier, here we report a new strategy enabling the specific incorporation of biochemically reactive amino acids into proteins. A latent bioreactive amino acid is genetically encoded at a position proximal to the target natural amino acid; they react via proximity-enabled reactivity, selectively converting the latter into a reactive residue in situ. Using this Genetically Encoded Chemical COnversion (GECCO) strategy and harnessing the sulfur-fluoride exchange (SuFEx) reaction between fluorosulfate-l-tyrosine and serine or threonine, we site-specifically generated the reactive dehydroalanine and dehydrobutyrine into proteins. GECCO works both inter- and intramolecularly, and is compatible with various proteins. We further labeled the resultant dehydroalanine-containing protein with thiol-saccharide to generate glycoprotein mimetics. GECCO represents a new solution for selectively introducing biochemically reactive amino acids into proteins and is expected to open new avenues for exploiting chemistry in live systems for biological research and engineering.

Infrared (IR) spectroscopy has provided considerable insight into the structures, dynamics, and formation mechanisms of amyloid fibrils. IR probes, such as main chain 13C═18O, have been widely employed to obtain site-specific structural information, yet only secondary structures and strand-to-strand arrangements can be probed. Very few nonperturbative IR probes are available to report on the side-chain conformation and environments, which are critical to determining sheet-to-sheet arrangements in steric zippers within amyloids. Polar residues, such as glutamine, contribute significantly to the stability of amyloids and thus are frequently found in core regions of amyloid peptides/proteins. Furthermore, polyglutamine (polyQ) repeats form toxic aggregates in several neurodegenerative diseases. Here we report the synthesis and application of a new nonperturbative IR probe-glutamine side chain 13C═18O. We use side chain 13C═18O labeling and isotope dilution to detect the presence of intermolecularly hydrogen-bonded arrays of glutamine side chains (Gln ladders) in amyloid-forming peptides. Moreover, the line width of the 13C═18O peak is highly sensitive to its local hydration environment. The IR data from side chain labeling allows us to unambiguously determine the sheet-to-sheet arrangement in a short amyloid-forming peptide, GNNQQNY, providing insight that was otherwise inaccessible through main chain labeling. With several different fibril samples, we also show the versatility of this IR probe in studying the structures and aggregation kinetics of amyloids. Finally, we demonstrate the capability of modeling amyloid structures with IR data using the integrative modeling platform (IMP) and the potential of integrating IR with other biophysical methods for more accurate structural modeling. Together, we believe that side chain 13C═18O will complement main chain isotope labeling in future IR studies of amyloids and integrative modeling using IR data will significantly expand the power of IR spectroscopy to elucidate amyloid assemblies.

The hallmarks of Alzheimer's disease (AD) are the accumulation of Aβ plaques and neurofibrillary tangles composed of hyperphosphorylated tau. We developed sensitive cellular assays using human embryonic kidney-293T cells to quantify intracellular self-propagating conformers of Aβ in brain samples from patients with AD or other neurodegenerative diseases. Postmortem brain tissue from patients with AD had measurable amounts of pathological Aβ conformers. Individuals over 80 years of age had the lowest amounts of prion-like Aβ and phosphorylated tau. Unexpectedly, the longevity-dependent decrease in self-propagating tau conformers occurred in spite of increasing amounts of total insoluble tau. When corrected for the abundance of insoluble tau, the ability of postmortem AD brain homogenates to induce misfolded tau in the cellular assays showed an exponential decrease with longevity, with a half-life of about one decade over the age range of 37 to 99 years. Thus, our findings demonstrate an inverse correlation between longevity in patients with AD and the abundance of pathological tau conformers. Our cellular assays can be applied to patient selection for clinical studies and the development of new drugs and diagnostics for AD.

Recently, l-4-cyanotryptophan has been shown to be an efficient blue fluorescence emitter, with the potential to enable novel applications in biological spectroscopy and microscopy. However, lack of facile synthetic routes to this unnatural amino acid limits its wide use. Herein, we describe an expedient approach to synthesize Fmoc protected l-4-cyanotryptophan with high optical purity (>99%). Additionally, we test the utility of this blue fluorophore in imaging cell-membrane-bound peptides and in determining peptide-membrane binding constants.

Site-specific protein cleavage is essential for many protein-production protocols and typically requires proteases. We report the development of a chemical protein-cleavage method that is achieved through the use of a sequence-specific nickel-assisted cleavage (SNAC)-tag. We demonstrate that the SNAC-tag can be inserted before both water-soluble and membrane proteins to achieve fusion protein cleavage under biocompatible conditions with efficiency comparable to that of enzymes, and that the method works even when enzymatic cleavages fail.

The features that stabilize the structures of membrane proteins remain poorly understood. Polar interactions contribute modestly, and the hydrophobic effect contributes little to the energetics of apolar side-chain packing in membranes. Disruption of steric packing can destabilize the native folds of membrane proteins, but is packing alone sufficient to drive folding in lipids? If so, then membrane proteins stabilized by this feature should be readily designed and structurally characterized-yet this has not been achieved. Through simulation of the natural protein phospholamban and redesign of variants, we define a steric packing code underlying its assembly. Synthetic membrane proteins designed using this code and stabilized entirely by apolar side chains conform to the intended fold. Although highly stable, the steric complementarity required for their folding is surprisingly stringent. Structural informatics shows that the designed packing motif recurs across the proteome, emphasizing a prominent role for precise apolar packing in membrane protein folding, stabilization, and evolution.

One of the fundamental events in protein folding is α-helix formation, which involves sequential development of a series of helical hydrogen bonds between the backbone C═O group of residues i and the -NH group of residues i + 4. While we now know a great deal about α-helix folding dynamics, a key question that remains to be answered is where the productive helical nucleation event occurs. Statistically, a helical nucleus (or the first helical hydrogen-bond) can form anywhere within the peptide sequence in question; however, the one that leads to productive folding may only form at a preferred location. This consideration is based on the fact that the α-helical structure is inherently asymmetric, due to the specific alignment of the helical hydrogen bonds. While this hypothesis is plausible, validating it is challenging because there is not an experimental observable that can be used to directly pinpoint the location of the productive nucleation process. Therefore, in this study we combine several techniques, including peptide cross-linking, laser-induced temperature-jump infrared spectroscopy, and molecular dynamics simulations, to tackle this challenge. Taken together, our experimental and simulation results support an α-helix folding mechanism wherein the productive nucleus is formed at the N-terminus, which propagates toward the C-terminal end of the peptide to yield the folded structure. In addition, our results show that incorporation of a cross-linker can lead to formation of differently folded conformations, underscoring the need for all-atom simulations to quantitatively assess the proposed cross-linking design.

Unnatural nucleosides possessing unique spectroscopic properties that mimic natural nucleobases in both size and chemical structure are ideally suited for spectroscopic measurements of DNA/RNA structure and dynamics in a site-specific manner. However, such unnatural nucleosides are scarce, which prompts us to explore the utility of a recently found unnatural nucleoside, 4-cyanoindole-2'-deoxyribonucleoside (4CNI-NS), as a site-specific spectroscopic probe of DNA. A recent study revealed that 4CNI-NS is a universal nucleobase that maintains the high fluorescence quantum yield of 4-cyanoindole and that among the four natural nucleobases, only guanine can significantly quench its fluorescence. Herein, we further show that the C≡N stretching frequency of 4CNI-NS is sensitive to the local environment, making it a useful site-specific infrared probe of oligonucleotides. In addition, we demonstrate that the fluorescence-quencher pair formed by 4CNI-NS and guanine can be used to quantitatively assess the binding affinity of a single-stranded DNA to the protein system of interest via fluorescence spectroscopy, among other applications. We believe that this fluorescence binding assay is especially useful as its potentiality allows high-throughput screening of DNA⁻protein interactions.

Water-mediated interactions play key roles in drug binding. In protein sites with sparse polar functionality, a small-molecule approach is often viewed as insufficient to achieve high affinity and specificity. Here we show that small molecules can enable potent inhibition by targeting key waters. The M2 proton channel of influenza A is the target of the antiviral drugs amantadine and rimantadine. Structural studies of drug binding to the channel using X-ray crystallography have been limited because of the challenging nature of the target, with the one previously solved crystal structure limited to 3.5 Å resolution. Here we describe crystal structures of amantadine bound to M2 in the Inwardclosed conformation (2.00 Å), rimantadine bound to M2 in both the Inwardclosed (2.00 Å) and Inwardopen (2.25 Å) conformations, and a spiro-adamantyl amine inhibitor bound to M2 in the Inwardclosed conformation (2.63 Å). These X-ray crystal structures of the M2 proton channel with bound inhibitors reveal that ammonium groups bind to water-lined sites that are hypothesized to stabilize transient hydronium ions formed in the proton-conduction mechanism. Furthermore, the ammonium and adamantyl groups of the adamantyl-amine class of drugs are free to rotate in the channel, minimizing the entropic cost of binding. These drug-bound complexes provide the first high-resolution structures of drugs that interact with and disrupt networks of hydrogen-bonded waters that are widely utilized throughout nature to facilitate proton diffusion within proteins.

Deregulated HER2 is a target of many approved cancer drugs. We analyzed 111,176 patient tumors and identified recurrent mutations in HER2 transmembrane domain (TMD) and juxtamembrane domain (JMD) that include G660D, R678Q, E693K, and Q709L. Using a saturation mutagenesis screen and testing of patient-derived mutations we found several activating TMD and JMD mutations. Structural modeling and analysis showed that the TMD/JMD mutations function by improving the active dimer interface or stabilizing an activating conformation. Further, we found that HER2 G660D employed asymmetric kinase dimerization for activation and signaling. Importantly, anti-HER2 antibodies and small-molecule kinase inhibitors blocked the activity of TMD/JMD mutants. Consistent with this, a G660D germline mutant lung cancer patient showed remarkable clinical response to HER2 blockade.

De novo design provides an attractive approach, which allows one to test and refine the principles guiding metalloproteins in defining the geometry and reactivity of their metal ion cofactors. Although impressive progress has been made in designing proteins that bind transition metal ions including iron-sulfur clusters, the design of tetranuclear clusters with oxygen-rich environments remains in its infancy. In previous work, we described the design of homotetrameric four-helix bundles that bind tetra-Zn2+ clusters. The crystal structures of the helical proteins were in good agreement with the overall design, and the metal-binding and conformational properties of the helical bundles in solution were consistent with the crystal structures. However, the corresponding apo-proteins were not fully folded in solution. In this work, we design three peptides, based on the crystal structure of the original bundles. One of the peptides forms tetramers in aqueous solution in the absence of metal ions as assessed by CD and NMR. It also binds Zn2+ in the intended stoichiometry. These studies strongly suggest that the desired structure has been achieved in the apo state, providing evidence that the peptide is able to actively impart the designed geometry to the metal cluster.

Chemical cross-linking mass spectrometry (CXMS) is being increasingly used to study protein assemblies and complex protein interaction networks. Existing CXMS chemical cross-linkers target only Lys, Cys, Glu, and Asp residues, limiting the information measurable. Here we report a "plant-and-cast" cross-linking strategy that employs a heterobifunctional cross-linker that contains a highly reactive succinimide ester as well as a less reactive sulfonyl fluoride. The succinimide ester reacts rapidly with surface Lys residues "planting" the reagent at fixed locations on protein. The pendant aryl sulfonyl fluoride is then "cast" across a limited range of the protein surface, where it can react with multiple weakly nucleophilic amino acid sidechains in a proximity-enhanced sulfur-fluoride exchange (SuFEx) reaction. Using proteins of known structures, we demonstrated that the heterobifunctional agent formed cross-links between Lys residues and His, Ser, Thr, Tyr, and Lys sidechains. This geometric specificity contrasts with current bis-succinimide esters, which often generate nonspecific cross-links between lysines brought into proximity by rare thermal fluctuations. Thus, the current method can provide diverse and robust distance restraints to guide integrative modeling. This work provides a chemical cross-linker targeting unactivated Ser, Thr, His, and Tyr residues using sulfonyl fluorides. In addition, this methodology yielded a variety of cross-links when applied to the complex Escherichia coli cell lysate. Finally, in combination with genetically encoded chemical cross-linking, cross-linking using this reagent markedly increased the identification of weak and transient enzyme-substrate interactions in live cells. Proximity-dependent cross-linking will dramatically expand the scope and power of CXMS for defining the identities and structures of protein complexes.

The pathways that G protein-coupled receptor (GPCR) ligands follow as they bind to or dissociate from their receptors are largely unknown. Protease-activated receptor-1 (PAR1) is a GPCR activated by intramolecular binding of a tethered agonist peptide that is exposed by thrombin cleavage. By contrast, the PAR1 antagonist vorapaxar is a lipophilic drug that binds in a pocket almost entirely occluded from the extracellular solvent. The binding and dissociation pathway of vorapaxar is unknown. Starting with the crystal structure of vorapaxar bound to PAR1, we performed temperature-accelerated molecular dynamics simulations of ligand dissociation. In the majority of simulations, vorapaxar exited the receptor laterally into the lipid bilayer through openings in the transmembrane helix (TM) bundle. Prior to full dissociation, vorapaxar paused in metastable intermediates stabilized by interactions with the receptor and lipid headgroups. Derivatives of vorapaxar with alkyl chains predicted to extend between TM6 and TM7 into the lipid bilayer inhibited PAR1 with apparent on rates similar to that of the parent compound in cell signaling assays. These data are consistent with vorapaxar binding to PAR1 via a pathway that passes between TM6 and TM7 from the lipid bilayer, in agreement with the most consistent pathway observed by molecular dynamics. While there is some evidence of entry of the ligand into rhodopsin and lipid-activated GPCRs from the cell membrane, our study provides the first such evidence for a peptide-activated GPCR and suggests that metastable intermediates along drug binding and dissociation pathways can be stabilized by specific interactions between lipids and the ligand.

Analogous to reversible post-translational protein modifications, the ability to attach and subsequently remove modifications on proteins would be valuable for protein and biological research. Although bioorthogonal functionalities have been developed to conjugate or cleave protein modifications, they are introduced into proteins on separate residues and often with bulky side chains, limiting their use to one type of control and primarily protein surface. Here we achieved dual control on one residue by genetically encoding S-propargyl-cysteine (SprC), which has bioorthogonal alkyne and propargyl groups in a compact structure, permitting usage in protein interior in addition to surface. We demonstrated its incorporation at the dimer interface of glutathione transferase for in vivo crosslinking via thiol-yne click chemistry, and at the active site of human rhinovirus 3C protease for masking and then turning on enzyme activity via Pd-cleavage of SprC into Cys. In addition, we installed biotin onto EGFP via Sonogashira coupling of SprC and then tracelessly removed it via Pd cleavage. SprC is small in size, commercially available, nontoxic, and allows for bond building and breaking on a single residue. Genetically encoded SprC will be valuable for chemically controlling proteins with an essential Cys and for reversible protein modifications.