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.

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.

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.

Elucidation of the molecular interactions that enable stable assembly and structural diversity in multicomponent immune receptors.

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.

Acyl Transfer Catalytic Activity in De Novo Designed Protein with N-Terminus of α-Helix As Oxyanion-Binding Site.

The design of catalytic proteins with functional sites capable of specific chemistry is gaining momentum and a number of artificial enzymes have recently been reported, including hydrolases, oxidoreductases, retro-aldolases, and others. Our goal is to develop a peptide ligase for robust catalysis of amide bond formation that possesses no stringent restrictions to the amino acid composition at the ligation junction. We report here the successful completion of the first step in this long-term project by building a completely de novo protein with predefined acyl transfer catalytic activity. We applied a minimalist approach to rationally design an oxyanion hole within a small cavity that contains an adjacent thiol nucleophile. The N-terminus of the α-helix with unpaired hydrogen-bond donors was exploited as a structural motif to stabilize negatively charged tetrahedral intermediates in nucleophilic addition-elimination reactions at the acyl group. Cysteine acting as a principal catalytic residue was introduced at the second residue position of the α-helix N-terminus in a designed three-α-helix protein based on structural informatics prediction. We showed that this minimal set of functional elements is sufficient for the emergence of catalytic activity in a de novo protein. Using peptide-αthioesters as acyl-donors, we demonstrated their catalyzed amidation concomitant with hydrolysis and proved that the environment at the catalytic site critically influences the reaction outcome. These results represent a promising starting point for the development of efficient catalysts for protein labeling, conjugation, and peptide ligation.

Modulating Integrin αIIbß3 Activity through Mutagenesis of Allosterically Regulated Intersubunit Contacts.

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.

Extending chemical perturbations of the ubiquitin fitness landscape in a classroom setting reveals new constraints on sequence tolerance.

Although the primary protein sequence of ubiquitin (Ub) is extremely stable over evolutionary time, it is highly tolerant to mutation during selection experiments performed in the laboratory. We have proposed that this discrepancy results from the difference between fitness under laboratory culture conditions and the selective pressures in changing environments over evolutionary timescales. Building on our previous work (Mavor et al., 2016), we used deep mutational scanning to determine how twelve new chemicals (3-Amino-1,2,4-triazole, 5-fluorocytosine, Amphotericin B, CaCl2, Cerulenin, Cobalt Acetate, Menadione, Nickel Chloride, p-Fluorophenylalanine, Rapamycin, Tamoxifen, and Tunicamycin) reveal novel mutational sensitivities of ubiquitin residues. Collectively, our experiments have identified eight new sensitizing conditions for Lys63 and uncovered a sensitizing condition for every position in Ub except Ser57 and Gln62. By determining the ubiquitin fitness landscape under different chemical constraints, our work helps to resolve the inconsistencies between deep mutational scanning experiments and sequence conservation over evolutionary timescales.

Computation-Guided Design of a Stimulus-Responsive Multienzyme Supramolecular Assembly.

The construction of stimulus-responsive supramolecular complexes of metabolic pathway enzymes, inspired by natural multienzyme assemblies (metabolons), provides an attractive avenue for efficient and spatiotemporally controllable one-pot biotransformations. We have constructed a phosphorylation- and optically responsive metabolon for the biodegradation of the environmental pollutant 1,2,3-trichloropropane.