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.

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.

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.

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.

Allosteric mechanism of signal transduction in the two-component system histidine kinase PhoQ.

Transmembrane signaling proteins couple extracytosolic sensors to cytosolic effectors. Here, we examine how binding of Mg2+ to the sensor domain of an E. coli two component histidine kinase (HK), PhoQ, modulates its cytoplasmic kinase domain. We use cysteine-crosslinking and reporter-gene assays to simultaneously and independently probe the signaling state of PhoQ's sensor and autokinase domains in a set of over 30 mutants. Strikingly, conservative single-site mutations distant from the sensor or catalytic site strongly influence PhoQ's ligand-sensitivity as well as the magnitude and direction of the signal. Data from 35 mutants are explained by a semi-empirical three-domain model in which the sensor, intervening HAMP, and catalytic domains can adopt kinase-promoting or inhibiting conformations that are in allosteric communication. The catalytic and sensor domains intrinsically favor a constitutively 'kinase-on' conformation, while the HAMP domain favors the 'off' state; when coupled, they create a bistable system responsive to physiological concentrations of Mg2+. Mutations alter signaling by locally modulating domain intrinsic equilibrium constants and interdomain couplings. Our model suggests signals transmit via interdomain allostery rather than propagation of a single concerted conformational change, explaining the diversity of signaling structural transitions observed in individual HK domains.

Allosteric cooperation in a de novo-designed two-domain protein.

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 Km and a 7-fold decrease in kcat These studies establish the feasibility of designing allosterically regulated catalytic proteins, entirely from scratch.

A defined structural unit enables de novo design of small-molecule-binding proteins.

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.

Excitation energy-dependent photocurrent switching in a single-molecule photodiode.

The direction of electron flow in molecular optoelectronic devices is dictated by charge transfer between a molecular excited state and an underlying conductor or semiconductor. For those devices, controlling the direction and reversibility of electron flow is a major challenge. We describe here a single-molecule photodiode. It is based on an internally conjugated, bichromophoric dyad with chemically linked (porphyrinato)zinc(II) and bis(terpyridyl)ruthenium(II) groups. On nanocrystalline, degenerately doped indium tin oxide electrodes, the dyad exhibits distinct frequency-dependent, charge-transfer characters. Variations in the light source between red-light (∼1.9 eV) and blue-light (∼2.7 eV) excitation for the integrated photodiode result in switching of photocurrents between cathodic and anodic. The origin of the excitation frequency-dependent photocurrents lies in the electronic structure of the chromophore excited states, as shown by the results of theoretical calculations, laser flash photolysis, and steady-state spectrophotometric measurements.

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.

Engineering opposite electronic polarization of singlet and triplet states increases the yield of high-energy photoproducts.

Efficient photosynthetic energy conversion requires quantitative, light-driven formation of high-energy, charge-separated states. However, energies of high-lying excited states are rarely extracted, in part because the congested density of states in the excited-state manifold leads to rapid deactivation. Conventional photosystem designs promote electron transfer (ET) by polarizing excited donor electron density toward the acceptor ("one-way" ET), a form of positive design. Curiously, negative design strategies that explicitly avoid unwanted side reactions have been underexplored. We report here that electronic polarization of a molecular chromophore can be used as both a positive and negative design element in a light-driven reaction. Intriguingly, prudent engineering of polarized excited states can steer a "U-turn" ET-where the excited electron density of the donor is initially pushed away from the acceptor-to outcompete a conventional one-way ET scheme. We directly compare one-way vs. U-turn ET strategies via a linked donor-acceptor (DA) assembly in which selective optical excitation produces donor excited states polarized either toward or away from the acceptor. Ultrafast spectroscopy of DA pinpoints the importance of realizing donor singlet and triplet excited states that have opposite electronic polarizations to shut down intersystem crossing. These results demonstrate that oppositely polarized electronically excited states can be employed to steer photoexcited states toward useful, high-energy products by routing these excited states away from states that are photosynthetic dead ends.

Genetically Introducing Biochemically Reactive Amino Acids Dehydroalanine and Dehydrobutyrine in Proteins.

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.

Inhibitor binding mode and allosteric regulation of Na+-glucose symporters.

Sodium-dependent glucose transporters (SGLTs) exploit sodium gradients to transport sugars across the plasma membrane. Due to their role in renal sugar reabsorption, SGLTs are targets for the treatment of type 2 diabetes. Current therapeutics are phlorizin derivatives that contain a sugar moiety bound to an aromatic aglycon tail. Here, we develop structural models of human SGLT1/2 in complex with inhibitors by combining computational and functional studies. Inhibitors bind with the sugar moiety in the sugar pocket and the aglycon tail in the extracellular vestibule. The binding poses corroborate mutagenesis studies and suggest a partial closure of the outer gate upon binding. The models also reveal a putative Na+ binding site in hSGLT1 whose disruption reduces the transport stoichiometry to the value observed in hSGLT2 and increases inhibition by aglycon tails. Our work demonstrates that subtype selectivity arises from Na+-regulated outer gate closure and a variable region in extracellular loop EL5.

Inhibitors of the M2 Proton Channel Engage and Disrupt Transmembrane Networks of Hydrogen-Bonded Waters.

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.

De novo design of a hyperstable non-natural protein-ligand complex with sub-Å accuracy.

Protein catalysis requires the atomic-level orchestration of side chains, substrates and cofactors, and yet the ability to design a small-molecule-binding protein entirely from first principles with a precisely predetermined structure has not been demonstrated. Here we report the design of a novel protein, PS1, that binds a highly electron-deficient non-natural porphyrin at temperatures up to 100 °C. The high-resolution structure of holo-PS1 is in sub-Å agreement with the design. The structure of apo-PS1 retains the remote core packing of the holoprotein, with a flexible binding region that is predisposed to ligand binding with the desired geometry. Our results illustrate the unification of core packing and binding-site definition as a central principle of ligand-binding protein design.

Molecular Road Map to Tuning Ground State Absorption and Excited State Dynamics of Long-Wavelength Absorbers.

Realizing chromophores that simultaneously possess substantial near-infrared (NIR) absorptivity and long-lived, high-yield triplet excited states is vital for many optoelectronic applications, such as optical power limiting and triplet-triplet annihilation photon upconversion (TTA-UC). However, the energy gap law ensures such chromophores are rare, and molecular engineering of absorbers having such properties has proven challenging. Here, we present a versatile methodology to tackle this design issue by exploiting the ethyne-bridged (polypyridyl)metal(II) (M; M = Ru, Os)-(porphinato)metal(II) (PM'; M' = Zn, Pt, Pd) molecular architecture (M-(PM')n-M), wherein high-oscillator-strength NIR absorptivity up to 850 nm, near-unity intersystem crossing (ISC) quantum yields (ΦISC), and triplet excited-state (T1) lifetimes on the microseconds time scale are simultaneously realized. By varying the extent to which the atomic coefficients of heavy metal d orbitals contribute to the one-electron excitation configurations describing the initially prepared singlet and triplet excited-state wave functions, we (i) show that the relative magnitudes of fluorescence (k0F), S1 → S0 nonradiative decay (knr), S1 → T1 ISC (kISC), and T1 → S0 relaxation (kT1→S0) rate constants can be finely tuned in M-(PM')n-M compounds and (ii) demonstrate designs in which the kISC magnitude dominates singlet manifold relaxation dynamics but does not give rise to T1 → S0 conversion dynamics that short-circuit a microseconds time scale triplet lifetime. Notably, the NIR spectral domain absorptivities of M-(PM')n-M chromophores far exceed those of classic coordination complexes and organic materials possessing similarly high yields of triplet-state formation: in contrast to these benchmark materials, this work demonstrates that these M-(PM')n-M systems realize near unit ΦISC at extraordinarily modest S1-T1 energy gaps (∼0.25 eV). This study underscores the photophysical diversity of the M-(PM')n-M platform and presents a new library of long-wavelength absorbers that efficiently populate long-lived T1 states.

Zinc-binding structure of a catalytic amyloid from solid-state NMR.

Throughout biology, amyloids are key structures in both functional proteins and the end product of pathologic protein misfolding. Amyloids might also represent an early precursor in the evolution of life because of their small molecular size and their ability to self-purify and catalyze chemical reactions. They also provide attractive backbones for advanced materials. When ß-strands of an amyloid are arranged parallel and in register, side chains from the same position of each chain align, facilitating metal chelation when the residues are good ligands such as histidine. High-resolution structures of metalloamyloids are needed to understand the molecular bases of metal-amyloid interactions. Here we combine solid-state NMR and structural bioinformatics to determine the structure of a zinc-bound metalloamyloid that catalyzes ester hydrolysis. The peptide forms amphiphilic parallel ß-sheets that assemble into stacked bilayers with alternating hydrophobic and polar interfaces. The hydrophobic interface is stabilized by apolar side chains from adjacent sheets, whereas the hydrated polar interface houses the Zn2+-binding histidines with binding geometries unusual in proteins. Each Zn2+ has two bis-coordinated histidine ligands, which bridge adjacent strands to form an infinite metal-ligand chain along the fibril axis. A third histidine completes the protein ligand environment, leaving a free site on the Zn2+ for water activation. This structure defines a class of materials, which we call metal-peptide frameworks. The structure reveals a delicate interplay through which metal ions stabilize the amyloid structure, which in turn shapes the ligand geometry and catalytic reactivity of Zn2.

Engineering High-Potential Photo-oxidants with Panchromatic Absorption.

Challenging photochemistry demands high-potential visible-light-absorbing photo-oxidants. We report (i) a highly electron-deficient Ru(II) complex (eDef-Rutpy) bearing an E1/20/+ potential more than 300 mV more positive than that of any established Ru(II) bis(terpyridyl) derivative, and (ii) an ethyne-bridged eDef-Rutpy-(porphinato)Zn(II) (eDef-RuPZn) supermolecule that affords both panchromatic UV-vis spectral domain absorptivity and a high E1/20/+ potential, comparable to that of Ce(NH4)2(NO3)6 [E1/2(Ce3+/4+) = 1.61 V vs NHE], a strong and versatile ground-state oxidant commonly used in organic functional group transformations. eDef-RuPZn exhibits ∼8-fold greater absorptive oscillator strength over the 380-700 nm range relative to conventional Ru(II) polypyridyl complexes, and impressive excited-state reduction potentials (1E-/* = 1.59 V; 3E-/* = 1.26 V). eDef-RuPZn manifests electronically excited singlet and triplet charge-transfer state lifetimes more than 2 orders of magnitude longer than those typical of conventional Ru(II) bis(terpyridyl) chromophores, suggesting new opportunities in light-driven oxidation reactions for energy conversion and photocatalysis.

Where Is the Electronic Oscillator Strength? Mapping Oscillator Strength across Molecular Absorption Spectra.

The effectiveness of solar energy capture and conversion materials derives from their ability to absorb light and to transform the excitation energy into energy stored in free carriers or chemical bonds. The Thomas-Reiche-Kuhn (TRK) sum rule mandates that the integrated (electronic) oscillator strength of an absorber equals the total number of electrons in the structure. Typical molecular chromophores place only about 1% of their oscillator strength in the UV-vis window, so individual chromophores operate at about 1% of their theoretical limit. We explore the distribution of oscillator strength as a function of excitation energy to understand this circumstance. To this aim, we use familiar independent-electron model Hamiltonians as well as first-principles electronic structure methods. While model Hamiltonians capture the qualitative electronic spectra associated with π electron chromophores, these Hamiltonians mistakenly focus the oscillator strength in the fewest low-energy transitions. Advanced electronic structure methods, in contrast, spread the oscillator strength over a very wide excitation energy range, including transitions to Rydberg and continuum states, consistent with experiment. Our analysis rationalizes the low oscillator strength in the UV-vis spectral region in molecules, a step toward the goal of oscillator strength manipulation and focusing.

Photoinduced Electron Transfer Elicits a Change in the Static Dielectric Constant of a de Novo Designed Protein.

We provide a direct measure of the change in effective dielectric constant (ε(S)) within a protein matrix after a photoinduced electron transfer (ET) reaction. A linked donor-bridge-acceptor molecule, PZn-Ph-NDI, consisting of a (porphinato)Zn donor (PZn), a phenyl bridge (Ph), and a naphthalene diimide acceptor (NDI), is shown to be a "meter" to indicate protein dielectric environment. We calibrated PZn-Ph-NDI ET dynamics as a function of solvent dielectric, and computationally de novo designed a protein SCPZnI3 to bind PZn-Ph-NDI in its interior. Mapping the protein ET dynamics onto the calibrated ET catalogue shows that SCPZnI3 undergoes a switch in the effective dielectric constant following photoinduced ET, from ε(S) ≈ 8 to ε(S) ≈ 3.

Mean First-Passage Times in Biology.

Many biochemical processes, such as charge hopping or protein folding, can be described by an average timescale to reach a final state, starting from an initial state. Here, we provide a pedagogical treatment of the mean first-passage time (MFPT) of a physical process, which depends on the number of intervening states between the initial state and the target state. Our aim in this tutorial review is to provide a clear development of the mean first-passage time formalism and to show some of its practical utility. The MFPT treatment can provide a useful link between microscopic rates and the average timescales often probed by experiment.