De novo design of drug-binding proteins with predictable binding energy and specificity.

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 <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule-binding proteins with tuned interaction energies is feasible entirely from computation.

De novo-designed transmembrane proteins bind and regulate a cytokine receptor.

Transmembrane (TM) domains as simple as a single span can perform complex biological functions using entirely lipid-embedded chemical features. Computational design has the potential to generate custom tool molecules directly targeting membrane proteins at their functional TM regions. Thus far, designed TM domain-targeting agents have been limited to mimicking the binding modes and motifs of natural TM interaction partners. Here, we demonstrate the design of de novo TM proteins targeting the erythropoietin receptor (EpoR) TM domain in a custom binding topology competitive with receptor homodimerization. The TM proteins expressed in mammalian cells complex with EpoR and inhibit erythropoietin-induced cell proliferation. In vitro, the synthetic TM domain complex outcompetes EpoR homodimerization. Structural characterization reveals that the complex involves the intended amino acids and agrees with our designed molecular model of antiparallel TM helices at 1:1 stoichiometry. Thus, membrane protein TM regions can now be targeted in custom-designed topologies.

An enhanced broad-spectrum peptide inhibits Omicron variants in vivo.

The continual emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) poses a major challenge to vaccines and antiviral therapeutics due to their extensive evasion of immunity. Aiming to develop potent and broad-spectrum anticoronavirus inhibitors, we generated A1-(GGGGS)7-HR2m (A1L35HR2m) by introducing an angiotensin-converting enzyme 2 (ACE2)-derived peptide A1 to the N terminus of the viral HR2-derived peptide HR2m through a long flexible linker, which showed significantly improved antiviral activity. Further cholesterol (Chol) modification at the C terminus of A1L35HR2m greatly enhanced the inhibitory activities against SARS-CoV-2, SARS-CoV-2 VOCs, SARS-CoV, and Middle East respiratory syndrome coronavirus (MERS-CoV) pseudoviruses, with IC50 values ranging from 0.16 to 5.53 nM. A1L35HR2m-Chol also potently inhibits spike-protein-mediated cell-cell fusion and the replication of authentic Omicron BA.2.12.1, BA.5, and EG.5.1. Importantly, A1L35HR2m-Chol distributed widely in respiratory tract tissue and had a long half-life (>10 h) in vivo. Intranasal administration of A1L35HR2m-Chol to K18-hACE2 transgenic mice potently inhibited Omicron BA.5 and EG.5.1 infection both prophylactically and therapeutically.

Docking for Molecules That Bind in a Symmetric Stack with SymDOCK.

Discovering ligands for amyloid fibrils, such as those formed by the tau protein, is an area of great current interest. In recent structures, ligands bind in stacks in the tau fibrils to reflect the rotational and translational symmetry of the fibril itself; in these structures, the ligands make few interactions with the protein but interact extensively with each other. To exploit this symmetry and stacking, we developed SymDOCK, a method to dock molecules that follow the protein's symmetry. For each prospective ligand pose, we apply the symmetry operation of the fibril to generate a self-interacting and fibril-interacting stack, checking that doing so will not cause a clash between the original molecule and its image. Absent a clash, we retain that pose and add the ligand-ligand van der Waals energy to the ligand's docking score (here using DOCK3.8). We can check these geometries and energies using an implementation of ANI, a neural-network-based quantum-mechanical evaluation of the ligand stacking energies. In retrospective calculations, symmetry docking can reproduce the poses of three tau PET tracers whose structures have been determined. More convincingly, in a prospective study, SymDOCK predicted the structure of the PET tracer MK-6240 bound in a symmetrical stack to AD PHF tau before that structure was determined; the docked pose was used to determine how MK-6240 fit the cryo-EM density. In proof-of-concept studies, SymDOCK enriched known ligands over property-matched decoys in retrospective screens without sacrificing docking speed and can address large library screens that seek new symmetrical stackers. Future applications of this approach will be considered.

Docking for molecules that bind in a symmetric stack with SymDOCK.

Discovering ligands for amyloid fibrils, such as those formed by the tau protein, is an area of much current interest. In recent structures, ligands bind in stacks in the tau fibrils to reflect the rotational and translational symmetry of the fibril itself; in these structures the ligands make few interactions with the protein but interact extensively with each other. To exploit this symmetry and stacking, we developed SymDOCK, a method to dock molecules that follow the protein's symmetry. For each prospective ligand pose, we apply the symmetry operation of the fibril to generate a self-interacting and fibril-interacting stack, checking that doing so will not cause a clash between the original molecule and its image. Absent a clash, we retain that pose and add the ligand-ligand van der Waals energy to the ligand's docking score (here using DOCK3.8). We can check these geometries and energies using an implementation of ANI, a neural network-based quantum-mechanical evaluation of the ligand stacking energies. In retrospective calculations, symmetry docking can reproduce the poses of three tau PET tracers whose structures have been determined. More convincingly, in a prospective study SymDOCK predicted the structure of the PET tracer MK-6240 bound in a symmetrical stack to AD PHF tau before that structure was determined; the docked pose was used to determine how MK-6240 fit the cryo-EM density. In proof-of-concept studies, SymDOCK enriched known ligands over property-matched decoys in retrospective screens without sacrificing docking speed, and can address large library screens that seek new symmetrical stackers. Future applications of this approach will be considered.

De novo Design of Peptides that Bind Specific Conformers of α-Synuclein.

Insoluble amyloids rich in cross-ß fibrils are observed in a number of neurodegenerative diseases. Depending on the clinicopathology, the amyloids can adopt distinct supramolecular assemblies, termed conformational strains. However, rapid methods to study amyloid in a conformationally specific manner are lacking. We introduce a novel computational method for de novo design of peptides that tile the surface of α-synuclein fibrils in a conformationally specific manner. Our method begins by identifying surfaces that are unique to the conformational strain of interest, which becomes a "target backbone" for the design of a peptide binder. Next, we interrogate structures in the PDB database with high geometric complementarity to the target. Then, we identify secondary structural motifs that interact with this target backbone in a favorable, highly occurring geometry. This method produces monomeric helical motifs with a favorable geometry for interaction with the strands of the underlying amyloid. Each motif is then symmetrically replicated to form a monolayer that tiles the amyloid surface. Finally, amino acid sequences of the peptide binders are computed to provide a sequence with high geometric and physicochemical complementarity to the target amyloid. This method was applied to a conformational strain of α-synuclein fibrils, resulting in a peptide with high specificity for the target relative to other amyloids formed by α-synuclein, tau, or Aß40. This designed peptide also markedly slowed the formation of α-synuclein amyloids. Overall, this method offers a new tool for examining conformational strains of amyloid proteins.

Hexamethylene amiloride binds the SARS-CoV-2 envelope protein at the protein-lipid interface.

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.

Patterning and folding of intestinal villi by active mesenchymal dewetting.

Tissue folding generates structural motifs critical to organ function. In the intestine, bending of a flat epithelium into a periodic pattern of folds gives rise to villi, the numerous finger-like protrusions that are essential for nutrient absorption. However, the molecular and mechanical mechanisms driving the initiation and morphogenesis of villi remain a matter of debate. Here, we identify an active mechanical mechanism that simultaneously patterns and folds intestinal villi. We find that PDGFRA+ subepithelial mesenchymal cells generate myosin II-dependent forces sufficient to produce patterned curvature in neighboring tissue interfaces. At the cell-level, this occurs through a process dependent upon matrix metalloproteinase-mediated tissue fluidization and altered cell-ECM adhesion. By combining computational models with in vivo experiments, we reveal these cellular features manifest at the tissue-level as differences in interfacial tensions that promote mesenchymal aggregation and interface bending through a process analogous to the active de-wetting of a thin liquid film.

Transient water wires mediate selective proton transport in designed channel proteins.

Selective proton transport through proteins is essential for forming and using proton gradients in cells. Protons are conducted along hydrogen-bonded 'wires' of water molecules and polar side chains, which, somewhat surprisingly, are often interrupted by dry apolar stretches in the conduction pathways, inferred from static protein structures. Here we hypothesize that protons are conducted through such dry spots by forming transient water wires, often highly correlated with the presence of the excess protons in the water wire. To test this hypothesis, we performed molecular dynamics simulations to design transmembrane channels with stable water pockets interspersed by apolar segments capable of forming flickering water wires. The minimalist designed channels conduct protons at rates similar to viral proton channels, and they are at least 106-fold more selective for H+ over Na+. These studies inform the mechanisms of biological proton conduction and the principles for engineering proton-conductive materials.

Stacked binding of a PET ligand to Alzheimer's tau paired helical filaments.

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.

A host defense peptide mimetic, brilacidin, potentiates caspofungin antifungal activity against human pathogenic fungi.

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.

EMBER multidimensional spectral microscopy enables quantitative determination of disease- and cell-specific amyloid strains.

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.

EMBER multi-dimensional spectral microscopy enables quantitative determination of disease- and cell-specific amyloid strains.

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 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. Significance: In neurodegenerative diseases proteins fold into amyloid structures with distinct conformations (strains) that are characteristic of different diseases. There is a need to rapidly identify these 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 Aß deposits in different transgenic mouse models. Our imaging method rapidly identifies conformational differences in Aß and tau deposits from Down syndrome, sporadic and familial Alzheimer's disease human brain slices. We also identified distinct conformational strains of tau inclusions in astrocytes, oligodendrocytes, and neurons from Pick's disease. These findings will facilitate the identification of pathogenic protein aggregates to guide research and treatment of protein misfolding diseases.

Designed Transmembrane Proteins Inhibit the Erythropoietin Receptor in a Custom Binding Topology.

Transmembrane (TM) domains as simple as a single span can perform complex biological functions using entirely lipid-embedded chemical features. Computational design has potential to generate custom tool molecules directly targeting membrane proteins at their functional TM regions. Thus far, designed TM domain-targeting agents have been limited to mimicking binding modes and motifs of natural TM interaction partners. Here, we demonstrate the design of de novo TM proteins targeting the erythropoietin receptor (EpoR) TM domain in a custom binding topology competitive with receptor homodimerization. The TM proteins expressed in mammalian cells complex with EpoR and inhibit erythropoietin-induced cell proliferation. In vitro, the synthetic TM domain complex outcompetes EpoR homodimerization. Structural characterization reveals that the complex involves the intended amino acids and agrees with our designed molecular model of antiparallel TM helices at 1:1 stoichiometry. Thus, membrane protein TM regions can now be targeted in custom designed topologies.

Aß and Tau Prions Causing Alzheimer's Disease.

Studies show that patients with Alzheimer's disease (AD) have both Aß and tau prions, and thus, AD is a double-prion disease. AD patients with the greatest longevity exhibited low levels of both Aß and tau prions; tau prions were nearly absent in the brains of almost half of the patients who lived beyond 80 years of age. Using cellular bioassays for prions in postmortem samples, we found that both Aß and tau proteins misfold into prions leading to AD, which is either a sporadic or familial dementing disorder. Although AD is transmissible experimentally, there is no evidence that AD is either communicable or contagious. Since the progression of AD correlates poorly with insoluble Aß in the central nervous system (CNS), it was difficult to distinguish between inert amyloids and Aß prions. To measure the progression of AD, we devised rapid bioassays to measure the abundance of isoform-specific Aß prions in the brains of transgenic (Tg) mice and in postmortem human CNS samples from AD victims and people who died of other neurodegenerative diseases (NDs). We found significant correlations between the longevity of individuals with AD, sex, and genetic background, despite the fact that all postmortem brain tissue had essentially the same confirmed neuropathology.Although brains from all AD patients had measurable levels of Aß prions at death, the oldest individuals had lower Aß prion levels than the younger ones. Additionally, the long-lived individuals had low tau prion levels that correlated with the extent of phosphorylated tau (p-tau). Unexpectedly, a longevity-dependent decrease in tau prions was found in spite of increasing amounts of total insoluble tau. When corrected for the abundance of insoluble tau, the tau prion levels decreased exponentially with respect to the age at death with a half-time of approximately one decade, and this correlated with the abundance of phosphorylated tau.Even though our findings with tau prions were not unexpected, they were counterintuitive; thus, tau phosphorylation and tau prion activity decreased exponentially with longevity in patients with AD ranging from ages 37 to 99 years. Our findings demonstrated an inverse correlation between longevity in AD patients and the abundance of neurotoxic tau prions. Moreover, our discovery may have profound implications for the selection of phenotypically distinct patient populations and the development of diagnostics and effective therapeutics for AD.

De novo design of drug-binding proteins with predictable binding energy and specificity.

The de novo design of small-molecule-binding proteins has seen exciting recent progress; however, the ability to achieve exquisite affinity for binding small molecules while tuning specificity has not yet been demonstrated directly from computation. Here, we develop a computational procedure that results in the highest affinity binders to date with predetermined relative affinities, targeting a series of PARP1 inhibitors. Two of four designed proteins bound with affinities ranging from < 5 nM to low µM, in a predictable manner. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free-energy calculations performed directly on the designed models are in excellent agreement with the experimentally measured affinities, suggesting that the de novo design of small-molecule-binding proteins with tuned interaction energies is now feasible entirely from computation. We expect these methods to open many opportunities in biomedicine, including rapid sensor development, antidote design, and drug delivery vehicles.

SARS-CoV-2 Envelope Protein Forms Clustered Pentamers in Lipid Bilayers.

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.

De novo metalloprotein design.

Natural metalloproteins perform many functions - ranging from sensing to electron transfer and catalysis - in which the position and property of each ligand and metal, is dictated by protein structure. De novo protein design aims to define an amino acid sequence that encodes a specific structure and function, providing a critical test of the hypothetical inner workings of (metallo)proteins. To date, de novo metalloproteins have used simple, symmetric tertiary structures - uncomplicated by the large size and evolutionary marks of natural proteins - to interrogate structure-function hypotheses. In this Review, we discuss de novo design applications, such as proteins that induce complex, increasingly asymmetric ligand geometries to achieve function, as well as the use of more canonical ligand geometries to achieve stability. De novo design has been used to explore how proteins fine-tune redox potentials and catalyse both oxidative and hydrolytic reactions. With an increased understanding of structure-function relationships, functional proteins including O2-dependent oxidases, fast hydrolases, and multi-proton/multi-electron reductases, have been created. In addition, proteins can now be designed using xeno-biological metals or cofactors and principles from inorganic chemistry to derive new-to-nature functions. These results and the advances in computational protein design suggest a bright future for the de novo design of diverse, functional metalloproteins.

Covalent labeling of a chromatin reader domain using proximity-reactive cyclic peptides.

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.

Quaternary structure independent folding of voltage-gated ion channel pore domain subunits.

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.