The C-terminal Region of D-DT Regulates Molecular Recognition for Protein-Ligand Complexes.

Systematic analysis of molecular recognition is critical for understanding the biological function of macromolecules. For the immunomodulatory protein D-dopachrome tautomerase (D-DT), the mechanism of protein-ligand interactions is poorly understood. Here, 17 carefully designed protein variants and wild type (WT) D-DT were interrogated with an array of complementary techniques to elucidate the structural basis of ligand recognition. Utilization of a substrate and two selective inhibitors with distinct binding profiles offered previously unseen mechanistic insights into D-DT-ligand interactions. Our results demonstrate that the C-terminal region serves a key role in molecular recognition via regulation of the active site opening, protein-ligand interactions, and conformational flexibility of the pocket's environment. While our study is the first comprehensive analysis of molecular recognition for D-DT, the findings reported herein promote the understanding of protein functionality and enable the design of new structure-based drug discovery projects.

Design of Efficient Artificial Enzymes Using Crystallographically Enhanced Conformational Sampling.

The ability to create efficient artificial enzymes for any chemical reaction is of great interest. Here, we describe a computational design method for increasing the catalytic efficiency of de novo enzymes by several orders of magnitude without relying on directed evolution and high-throughput screening. Using structural ensembles generated from dynamics-based refinement against X-ray diffraction data collected from crystals of Kemp eliminases HG3 (kcat/KM 125 M-1 s-1) and KE70 (kcat/KM 57 M-1 s-1), we design from each enzyme ≤10 sequences predicted to catalyze this reaction more efficiently. The most active designs display kcat/KM values improved by 100-250-fold, comparable to mutants obtained after screening thousands of variants in multiple rounds of directed evolution. Crystal structures show excellent agreement with computational models, with catalytic contacts present as designed and transition-state root-mean-square deviations of ≤0.65 Å. Our work shows how ensemble-based design can generate efficient artificial enzymes by exploiting the true conformational ensemble to design improved active sites.

Changes in an enzyme ensemble during catalysis observed by high-resolution XFEL crystallography.

Enzymes populate ensembles of structures necessary for catalysis that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography at an x-ray free electron laser to observe catalysis in a designed mutant isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations, and formation of the thioimidate intermediate selects for catalytically competent substates. The influence of cysteine ionization on the ICH ensemble is validated by determining structures of the enzyme at multiple pH values. Large molecular dynamics simulations in crystallo and time-resolved electron density maps show that Asp17 ionizes during catalysis and causes conformational changes that propagate across the dimer, permitting water to enter the active site for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.

Structural characterization of ligand binding and pH-specific enzymatic activity of mouse Acidic Mammalian Chitinase.

Chitin is an abundant biopolymer and pathogen-associated molecular pattern that stimulates a host innate immune response. Mammals express chitin-binding and chitin-degrading proteins to remove chitin from the body. One of these proteins, Acidic Mammalian Chitinase (AMCase), is an enzyme known for its ability to function under acidic conditions in the stomach but is also active in tissues with more neutral pHs, such as the lung. Here, we used a combination of biochemical, structural, and computational modeling approaches to examine how the mouse homolog (mAMCase) can act in both acidic and neutral environments. We measured kinetic properties of mAMCase activity across a broad pH range, quantifying its unusual dual activity optima at pH 2 and 7. We also solved high resolution crystal structures of mAMCase in complex with oligomeric GlcNAcn, the building block of chitin, where we identified extensive conformational ligand heterogeneity. Leveraging these data, we conducted molecular dynamics simulations that suggest how a key catalytic residue could be protonated via distinct mechanisms in each of the two environmental pH ranges. These results integrate structural, biochemical, and computational approaches to deliver a more complete understanding of the catalytic mechanism governing mAMCase activity at different pH. Engineering proteins with tunable pH optima may provide new opportunities to develop improved enzyme variants, including AMCase, for therapeutic purposes in chitin degradation.

Design of efficient artificial enzymes using crystallographically-enhanced conformational sampling.

The ability to create efficient artificial enzymes for any chemical reaction is of great interest. Here, we describe a computational design method for increasing catalytic efficiency of de novo enzymes to a level comparable to their natural counterparts without relying on directed evolution. Using structural ensembles generated from dynamics-based refinement against X-ray diffraction data collected from crystals of Kemp eliminases HG3 (kcat/KM 125 M-1 s-1) and KE70 (kcat/KM 57 M-1 s-1), we design from each enzyme ≤10 sequences predicted to catalyze this reaction more efficiently. The most active designs display kcat/KM values improved by 100-250-fold, comparable to mutants obtained after screening thousands of variants in multiple rounds of directed evolution. Crystal structures show excellent agreement with computational models. Our work shows how computational design can generate efficient artificial enzymes by exploiting the true conformational ensemble to more effectively stabilize the transition state.

Computational remodeling of an enzyme conformational landscape for altered substrate selectivity.

Structural plasticity of enzymes dictates their function. Yet, our ability to rationally remodel enzyme conformational landscapes to tailor catalytic properties remains limited. Here, we report a computational procedure for tuning conformational landscapes that is based on multistate design of hinge-mediated domain motions. Using this method, we redesign the conformational landscape of a natural aminotransferase to preferentially stabilize a less populated but reactive conformation and thereby increase catalytic efficiency with a non-native substrate, resulting in altered substrate selectivity. Steady-state kinetics of designed variants reveals activity increases with the non-native substrate of approximately 100-fold and selectivity switches of up to 1900-fold. Structural analyses by room-temperature X-ray crystallography and multitemperature nuclear magnetic resonance spectroscopy confirm that conformational equilibria favor the target conformation. Our computational approach opens the door to targeted alterations of conformational states and equilibria, which should facilitate the design of biocatalysts with customized activity and selectivity.

Combining temperature perturbations with X-ray crystallography to study dynamic macromolecules: A thorough discussion of experimental methods.

Temperature is an important state variable that governs the behavior of microscopic systems, yet crystallographers rarely exploit temperature changes to study the structure and dynamics of biological macromolecules. In fact, approximately 90% of crystal structures in the Protein Data Bank were determined under cryogenic conditions, because sample cryocooling makes crystals robust to X-ray radiation damage and facilitates data collection. On the other hand, cryocooling can introduce artifacts into macromolecular structures, and can suppress conformational dynamics that are critical for function. Fortunately, recent advances in X-ray detector technology, X-ray sources, and computational data processing algorithms make non-cryogenic X-ray crystallography easier and more broadly applicable than ever before. Without the reliance on cryocooling, high-resolution crystallography can be combined with various temperature perturbations to gain deep insight into the conformational landscapes of macromolecules. This Chapter reviews the historical reasons for the prevalence of cryocooling in macromolecular crystallography, and discusses its potential drawbacks. Next, the Chapter summarizes technological developments and methodologies that facilitate non-cryogenic crystallography experiments. Finally, the chapter discusses the theoretical underpinnings and practical aspects of multi-temperature and temperature-jump crystallography experiments, which are powerful tools for understanding the relationship between the structure, dynamics, and function of proteins and other biological macromolecules.

Mapping protein dynamics at high spatial resolution with temperature-jump X-ray crystallography.

Understanding and controlling protein motion at atomic resolution is a hallmark challenge for structural biologists and protein engineers because conformational dynamics are essential for complex functions such as enzyme catalysis and allosteric regulation. Time-resolved crystallography offers a window into protein motions, yet without a universal perturbation to initiate conformational changes the method has been limited in scope. Here we couple a solvent-based temperature jump with time-resolved crystallography to visualize structural motions in lysozyme, a dynamic enzyme. We observed widespread atomic vibrations on the nanosecond timescale, which evolve on the submillisecond timescale into localized structural fluctuations that are coupled to the active site. An orthogonal perturbation to the enzyme, inhibitor binding, altered these dynamics by blocking key motions that allow energy to dissipate from vibrations into functional movements linked to the catalytic cycle. Because temperature jump is a universal method for perturbing molecular motion, the method demonstrated here is broadly applicable for studying protein dynamics.

Changes in an Enzyme Ensemble During Catalysis Observed by High Resolution XFEL Crystallography.

Enzymes populate ensembles of structures with intrinsically different catalytic proficiencies that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography (MISC) at an X-ray free electron laser (XFEL) to observe catalysis in a designed mutant (G150T) isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations and formation of the thioimidate catalytic intermediate selects for catalytically competent substates. A prior proposal for active site cysteine charge-coupled conformational changes in ICH is validated by determining structures of the enzyme over a range of pH values. A combination of large molecular dynamics simulations of the enzyme in crystallo and time-resolved electron density maps shows that ionization of the general acid Asp17 during catalysis causes additional conformational changes that propagate across the dimer interface, connecting the two active sites. These ionization-linked changes in the ICH conformational ensemble permit water to enter the active site in a location that is poised for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and x-ray diffraction at room temperature.

The nonstructural protein 3 (NSP3) macrodomain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Mac1) removes adenosine diphosphate (ADP) ribosylation posttranslational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the coronavirus disease 2019 pandemic. Here, we determined neutron and x-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase the potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a reevaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

Accurate positioning of functional residues with robotics-inspired computational protein design.

SignificanceComputational protein design promises to advance applications in medicine and biotechnology by creating proteins with many new and useful functions. However, new functions require the design of specific and often irregular atom-level geometries, which remains a major challenge. Here, we develop computational methods that design and predict local protein geometries with greater accuracy than existing methods. Then, as a proof of concept, we leverage these methods to design new protein conformations in the enzyme ketosteroid isomerase that change the protein's preference for a key functional residue. Our computational methods are openly accessible and can be applied to the design of other intricate geometries customized for new user-defined protein functions.

Generation of bright monomeric red fluorescent proteins via computational design of enhanced chromophore packing.

Red fluorescent proteins (RFPs) have found widespread application in chemical and biological research due to their longer emission wavelengths. Here, we use computational protein design to increase the quantum yield and thereby brightness of a dim monomeric RFP (mRojoA, quantum yield = 0.02) by optimizing chromophore packing with aliphatic residues, which we hypothesized would reduce torsional motions causing non-radiative decay. Experimental characterization of the top 10 designed sequences yielded mSandy1 (λ em = 609 nm, quantum yield = 0.26), a variant with equivalent brightness to mCherry, a widely used RFP. We next used directed evolution to further increase brightness, resulting in mSandy2 (λ em = 606 nm, quantum yield = 0.35), the brightest Discosoma sp. derived monomeric RFP with an emission maximum above 600 nm reported to date. Crystallographic analysis of mSandy2 showed that the chromophore p-hydroxybenzylidene moiety is sandwiched between the side chains of Leu63 and Ile197, a structural motif that has not previously been observed in RFPs, and confirms that aliphatic packing leads to chromophore rigidification. Our results demonstrate that computational protein design can be used to generate bright monomeric RFPs, which can serve as templates for the evolution of novel far-red fluorescent proteins.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and X-ray diffraction at room temperature.

The NSP3 macrodomain of SARS CoV 2 (Mac1) removes ADP-ribosylation post-translational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the COVID-19 pandemic. Here, we determined neutron and X-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site, and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a re-evaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

Identification of Isoprene Oxidation Reaction Products via Anion Photoelectron Spectroscopy.

We present a study on the oxidation of isoprene under several different conditions that may model both atmospheric and combustion chemistry. Anions, formed by passing isoprene/oxidant gas mixtures through a pulsed discharge generating a range of species, are separated via mass spectrometry and characterized by anion photoelectron (PE) spectroscopy supported by computations. Specifically, a UV-irradiated isoprene/O2 mixture, which additionally produces O3, and an isoprene/O2/H2 mixture, which generates •OH when passed through the discharge, were sampled. The mass spectra of ions generated under both conditions show the production of intact molecular ions, ion-molecule complexes (e.g., O2-, O4-, and O2-·isoprene), and singly deprotonated species (e.g., deprotonated isoprene, C5H7-). In addition, both smaller and oxidized fragments are observed using both gas mixtures, though relative abundances differ. From the UV-irradiated isoprene/O2 gas mixture, additional intact molecular products of reactions initiated by ozonolysis of isoprene, methylglyoxal, and dimethylglyoxal were observed. Fragmentation and oxidation of isoprene observed in both gas mixtures included species with m/z 39, 53, 67, 69, and 83 that we attribute to a series of alkyl- and alkenoxide-based anions. The coexistence of intact molecules and complexes with fragments and reaction products demonstrates the versatility of this ion source as a simple and efficient anion formation method for studying species that may be relevant in atmospheric and combustion chemistry.

Fragment binding to the Nsp3 macrodomain of SARS-CoV-2 identified through crystallographic screening and computational docking.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) macrodomain within the nonstructural protein 3 counteracts host-mediated antiviral adenosine diphosphate-ribosylation signaling. This enzyme is a promising antiviral target because catalytic mutations render viruses nonpathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of 2533 diverse fragments resulted in 214 unique macrodomain-binders. An additional 60 molecules were selected from docking more than 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several fragment hits were confirmed by solution binding using three biophysical techniques (differential scanning fluorimetry, homogeneous time-resolved fluorescence, and isothermal titration calorimetry). The 234 fragment structures explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

Fragment Binding to the Nsp3 Macrodomain of SARS-CoV-2 Identified Through Crystallographic Screening and Computational Docking.

The SARS-CoV-2 macrodomain (Mac1) within the non-structural protein 3 (Nsp3) counteracts host-mediated antiviral ADP-ribosylation signalling. This enzyme is a promising antiviral target because catalytic mutations render viruses non-pathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of diverse fragment libraries resulted in 214 unique macrodomain-binding fragments, out of 2,683 screened. An additional 60 molecules were selected from docking over 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several crystallographic and docking fragment hits were validated for solution binding using three biophysical techniques (DSF, HTRF, ITC). Overall, the 234 fragment structures presented explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

An ultrapotent synthetic nanobody neutralizes SARS-CoV-2 by stabilizing inactive Spike.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus enters host cells via an interaction between its Spike protein and the host cell receptor angiotensin-converting enzyme 2 (ACE2). By screening a yeast surface-displayed library of synthetic nanobody sequences, we developed nanobodies that disrupt the interaction between Spike and ACE2. Cryo-electron microscopy (cryo-EM) revealed that one nanobody, Nb6, binds Spike in a fully inactive conformation with its receptor binding domains locked into their inaccessible down state, incapable of binding ACE2. Affinity maturation and structure-guided design of multivalency yielded a trivalent nanobody, mNb6-tri, with femtomolar affinity for Spike and picomolar neutralization of SARS-CoV-2 infection. mNb6-tri retains function after aerosolization, lyophilization, and heat treatment, which enables aerosol-mediated delivery of this potent neutralizer directly to the airway epithelia.

Ensemble-based enzyme design can recapitulate the effects of laboratory directed evolution in silico.

The creation of artificial enzymes is a key objective of computational protein design. Although de novo enzymes have been successfully designed, these exhibit low catalytic efficiencies, requiring directed evolution to improve activity. Here, we use room-temperature X-ray crystallography to study changes in the conformational ensemble during evolution of the designed Kemp eliminase HG3 (kcat/KM 146 M-1s-1). We observe that catalytic residues are increasingly rigidified, the active site becomes better pre-organized, and its entrance is widened. Based on these observations, we engineer HG4, an efficient biocatalyst (kcat/KM 103,000 M-1s-1) containing key first and second-shell mutations found during evolution. HG4 structures reveal that its active site is pre-organized and rigidified for efficient catalysis. Our results show how directed evolution circumvents challenges inherent to enzyme design by shifting conformational ensembles to favor catalytically-productive sub-states, and suggest improvements to the design methodology that incorporate ensemble modeling of crystallographic data.

Expanding the space of protein geometries by computational design of de novo fold families.

Naturally occurring proteins vary the precise geometries of structural elements to create distinct shapes optimal for function. We present a computational design method, loop-helix-loop unit combinatorial sampling (LUCS), that mimics nature's ability to create families of proteins with the same overall fold but precisely tunable geometries. Through near-exhaustive sampling of loop-helix-loop elements, LUCS generates highly diverse geometries encompassing those found in nature but also surpassing known structure space. Biophysical characterization showed that 17 (38%) of 45 tested LUCS designs encompassing two different structural topologies were well folded, including 16 with designed non-native geometries. Four experimentally solved structures closely matched the designs. LUCS greatly expands the designable structure space and offers a new paradigm for designing proteins with tunable geometries that may be customizable for novel functions.

An ultra-potent synthetic nanobody neutralizes SARS-CoV-2 by locking Spike into an inactive conformation.

Without an effective prophylactic solution, infections from SARS-CoV-2 continue to rise worldwide with devastating health and economic costs. SARS-CoV-2 gains entry into host cells via an interaction between its Spike protein and the host cell receptor angiotensin converting enzyme 2 (ACE2). Disruption of this interaction confers potent neutralization of viral entry, providing an avenue for vaccine design and for therapeutic antibodies. Here, we develop single-domain antibodies (nanobodies) that potently disrupt the interaction between the SARS-CoV-2 Spike and ACE2. By screening a yeast surface-displayed library of synthetic nanobody sequences, we identified a panel of nanobodies that bind to multiple epitopes on Spike and block ACE2 interaction via two distinct mechanisms. Cryogenic electron microscopy (cryo-EM) revealed that one exceptionally stable nanobody, Nb6, binds Spike in a fully inactive conformation with its receptor binding domains (RBDs) locked into their inaccessible down-state, incapable of binding ACE2. Affinity maturation and structure-guided design of multivalency yielded a trivalent nanobody, mNb6-tri, with femtomolar affinity for SARS-CoV-2 Spike and picomolar neutralization of SARS-CoV-2 infection. mNb6-tri retains stability and function after aerosolization, lyophilization, and heat treatment. These properties may enable aerosol-mediated delivery of this potent neutralizer directly to the airway epithelia, promising to yield a widely deployable, patient-friendly prophylactic and/or early infection therapeutic agent to stem the worst pandemic in a century.