The impact of molecular variants, crystallization conditions and the space group on ligand-protein complexes: a case study on bacterial phosphotriesterase.

A bacterial phosphotriesterase was employed as an experimental paradigm to examine the effects of multiple factors, such as the molecular constructs, the ligands used during protein expression and purification, the crystallization conditions and the space group, on the visualization of molecular complexes of ligands with a target enzyme. In this case, the ligands used were organophosphates that are fragments of the nerve agents and insecticides on which the enzyme acts as a bioscavenger. 12 crystal structures of various phosphotriesterase constructs obtained by directed evolution were analyzed, with resolutions of up to 1.38 Å. Both apo forms and holo forms, complexed with the organophosphate ligands, were studied. Crystals obtained from three different crystallization conditions, crystallized in four space groups, with and without N-terminal tags, were utilized to investigate the impact of these factors on visualizing the organophosphate complexes of the enzyme. The study revealed that the tags used for protein expression can lodge in the active site and hinder ligand binding. Furthermore, the space group in which the protein crystallizes can significantly impact the visualization of bound ligands. It was also observed that the crystallization precipitants can compete with, and even preclude, ligand binding, leading to false positives or to the incorrect identification of lead drug candidates. One of the co-crystallization conditions enabled the definition of the spaces that accommodate the substituents attached to the P atom of several products of organophosphate substrates after detachment of the leaving group. The crystal structures of the complexes of phosphotriesterase with the organophosphate products reveal similar short interaction distances of the two partially charged O atoms of the P-O bonds with the exposed ß-Zn2+ ion and the buried α-Zn2+ ion. This suggests that both Zn2+ ions have a role in stabilizing the transition state for substrate hydrolysis. Overall, this study provides valuable insights into the challenges and considerations involved in studying the crystal structures of ligand-protein complexes, highlighting the importance of careful experimental design and rigorous data analysis in ensuring the accuracy and reliability of the resulting phosphotriesterase-organophosphate structures.

Specific quinone reductase 2 inhibitors reduce metabolic burden and reverse Alzheimer's disease phenotype in mice.

Biological aging can be described as accumulative, prolonged metabolic stress and is the major risk factor for cognitive decline and Alzheimer's disease (AD). Recently, we identified and described a quinone reductase 2 (QR2) pathway in the brain, in which QR2 acts as a removable memory constraint and metabolic buffer within neurons. QR2 becomes overexpressed with age, and it is possibly a novel contributing factor to age-related metabolic stress and cognitive deficit. We found that, in human cells, genetic removal of QR2 produced a shift in the proteome opposing that found in AD brains while simultaneously reducing oxidative stress. We therefore created highly specific QR2 inhibitors (QR2is) to enable evaluation of chronic QR2 inhibition as a means to reduce biological age-related metabolic stress and cognitive decline. QR2is replicated results obtained by genetic removal of QR2, while local QR2i microinjection improved hippocampal and cortical-dependent learning in rats and mice. Continuous consumption of QR2is in drinking water improved cognition and reduced pathology in the brains of AD-model mice (5xFAD), with a noticeable between-sex effect on treatment duration. These results demonstrate the importance of QR2 activity and pathway function in the healthy and neurodegenerative brain and what we believe to be the great therapeutic potential of QR2is as first-in-class drugs.

Allosteric activation of vinculin by talin.

The talin-vinculin axis is a key mechanosensing component of cellular focal adhesions. How talin and vinculin respond to forces and regulate one another remains unclear. By combining single-molecule magnetic tweezers experiments, Molecular Dynamics simulations, actin-bundling assays, and adhesion assembly experiments in live cells, we here describe a two-ways allosteric network within vinculin as a regulator of the talin-vinculin interaction. We directly observe a maturation process of vinculin upon talin binding, which reinforces the binding to talin at a rate of 0.03 s-1. This allosteric transition can compete with force-induced dissociation of vinculin from talin only at forces up to 10 pN. Mimicking the allosteric activation by mutation yields a vinculin molecule that bundles actin and localizes to focal adhesions in a force-independent manner. Hence, the allosteric switch confines talin-vinculin interactions and focal adhesion build-up to intermediate force levels. The 'allosteric vinculin mutant' is a valuable molecular tool to further dissect the mechanical and biochemical signalling circuits at focal adhesions and elsewhere.

Design of a stable human acid-ß-glucosidase: towards improved Gaucher disease therapy and mutation classification.

Acid-ß-glucosidase (GCase, EC3.2.1.45), the lysosomal enzyme which hydrolyzes the simple glycosphingolipid, glucosylceramide (GlcCer), is encoded by the GBA1 gene. Biallelic mutations in GBA1 cause the human inherited metabolic disorder, Gaucher disease (GD), in which GlcCer accumulates, while heterozygous GBA1 mutations are the highest genetic risk factor for Parkinson's disease (PD). Recombinant GCase (e.g., Cerezyme® ) is produced for use in enzyme replacement therapy for GD and is largely successful in relieving disease symptoms, except for the neurological symptoms observed in a subset of patients. As a first step toward developing an alternative to the recombinant human enzymes used to treat GD, we applied the PROSS stability-design algorithm to generate GCase variants with enhanced stability. One of the designs, containing 55 mutations compared to wild-type human GCase, exhibits improved secretion and thermal stability. Furthermore, the design has higher enzymatic activity than the clinically used human enzyme when incorporated into an AAV vector, resulting in a larger decrease in the accumulation of lipid substrates in cultured cells. Based on stability-design calculations, we also developed a machine learning-based approach to distinguish benign from deleterious (i.e., disease-causing) GBA1 mutations. This approach gave remarkably accurate predictions of the enzymatic activity of single-nucleotide polymorphisms in the GBA1 gene that are not currently associated with GD or PD. This latter approach could be applied to other diseases to determine risk factors in patients carrying rare mutations.

The identification and characterization of an oxalyl-CoA synthetase from grass pea (Lathyrus sativus L.).

Oxalic acid is a small metabolite found in many plants. It serves as protection from herbivores, a chelator of metal ions, a regulator of calcium levels, and additional tasks. However, it is also a strong di-carboxylic acid that can compromise plant viability by reducing cellular pH. Several metabolic pathways have evolved to control oxalate levels in plants by enzymatic degradation. Among them is the pathway that utilizes oxalyl-CoA synthetase (OCS, EC 6.2.1.8) and ATP to convert oxalate to oxalyl-CoA. Oxalyl-CoA can then be degraded to CO2 or utilized as a precursor for the synthesis of other compounds. In grass pea (Lathyrus sativus L.), a grain legume grown in Asia and Africa for human and animal consumption, the neurotoxic compound ß-N-oxalyl-l-α,ß-diaminopropionic acid (ß-ODAP) is synthesized from oxalyl-CoA and l-α,ß-diaminopropionic acid (l-DAPA). Here, we report on the identification and characterization of oxalyl CoA-synthetase from grass pea (LsOCS). The gene encoding LsOCS was amplified from grass pea, and then expressed and purified from E. coli cells as an untagged, monomeric protein of 56 kDa. Its catalytic efficiency with oxalate, K oxalate M = 71.5 ± 13.3 µM, V max = 8.2 ± 0.8 µmole min-1 mg-1, was similar to that of OCS homologs from Arabidopsis thaliana (AtAAE3) and Medicago truncatula (MtAAE3). The enzyme was crystalized in complex with AMP and is the first OCS whose structure was determined in the thioester-forming conformation. Finally, we propose that substituting LsOCS with an oxalate oxidase or decarboxylase may reduce the levels of ß-ODAP in grass pea.

Identification and characterization of the key enzyme in the biosynthesis of the neurotoxin ß-ODAP in grass pea.

Grass pea (Lathyrus sativus L.) is a grain legume commonly grown in Asia and Africa for food and forage. It is a highly nutritious and robust crop, capable of surviving both droughts and floods. However, it produces a neurotoxic compound, ß-N-oxalyl-L-α,ß-diaminopropionic acid (ß-ODAP), which can cause a severe neurological disorder when consumed as a primary diet component. While the catalytic activity associated with ß-ODAP formation was demonstrated more than 50 years ago, the enzyme responsible for this activity has not been identified. Here, we report on the identity, activity, 3D structure, and phylogenesis of this enzyme-ß-ODAP synthase (BOS). We show that BOS belongs to the benzylalcohol O-acetyltransferase, anthocyanin O-hydroxycinnamoyltransferase, anthranilate N-hydroxycinnamoyl/benzoyltransferase, deacetylvindoline 4-O-acetyltransferase superfamily of acyltransferases and is structurally similar to hydroxycinnamoyl transferase. Using molecular docking, we propose a mechanism for its catalytic activity, and using heterologous expression in tobacco leaves (Nicotiana benthamiana), we demonstrate that expression of BOS in the presence of its substrates is sufficient for ß-ODAP production in vivo. The identification of BOS may pave the way toward engineering ß-ODAP-free grass pea cultivars, which are safe for human and animal consumption.

Computationally designed dual-color MRI reporters for noninvasive imaging of transgene expression.

Imaging of gene-expression patterns in live animals is difficult to achieve with fluorescent proteins because tissues are opaque to visible light. Imaging of transgene expression with magnetic resonance imaging (MRI), which penetrates to deep tissues, has been limited by single reporter visualization capabilities. Moreover, the low-throughput capacity of MRI limits large-scale mutagenesis strategies to improve existing reporters. Here we develop an MRI system, called GeneREFORM, comprising orthogonal reporters for two-color imaging of transgene expression in deep tissues. Starting from two promiscuous deoxyribonucleoside kinases, we computationally designed highly active, orthogonal enzymes ('reporter genes') that specifically phosphorylate two MRI-detectable synthetic deoxyribonucleosides ('reporter probes'). Systemically administered reporter probes exclusively accumulate in cells expressing the designed reporter genes, and their distribution is displayed as pseudo-colored MRI maps based on dynamic proton exchange for noninvasive visualization of transgene expression. We envision that future extensions of GeneREFORM will pave the way to multiplexed deep-tissue mapping of gene expression in live animals.

Evolution of CPEB4 Dynamics Across its Liquid-Liquid Phase Separation Transition.

Knowledge about the structural and dynamic properties of proteins that form membrane-less organelles in cells via liquid-liquid phase separation (LLPS) is required for understanding the process at a molecular level. We used spin labeling and electron paramagnetic resonance (EPR) spectroscopy to investigate the dynamic properties (rotational diffusion) of the low complexity N-terminal domain of cytoplasmic polyadenylation element binding-4 protein (CPEB4NTD) across its LLPS transition, which takes place with increasing temperature. We report the coexistence of three spin labeled CPEB4NTD (CPEB4*) populations with distinct dynamic properties representing different conformational spaces, both before and within the LLPS state. Monomeric CPEB4* exhibiting fast motion defines population I and shows low abundance prior to and following LLPS. Populations II and III are part of CPEB4* assemblies where II corresponds to loose conformations with intermediate range motions and population III represents compact conformations with strongly attenuated motions. As the temperature increased the population of component II increased reversibly at the expense of component III, indicating the existence of an III ⇌ II equilibrium. We correlated the macroscopic LLPS properties with the III ⇌ II exchange process upon varying temperature and CPEB4* and salt concentrations. We hypothesized that weak transient intermolecular interactions facilitated by component II lead to LLPS, with the small assemblies integrated within the droplets. The LLPS transition, however, was not associated with a clear discontinuity in the correlation times and populations of the three components. Importantly, CPEB4NTD exhibits LLPS properties where droplet formation occurs from a preformed microscopic assembly rather than the monomeric protein molecules.

Community-Wide Experimental Evaluation of the PROSS Stability-Design Method.

Recent years have seen a dramatic improvement in protein-design methodology. Nevertheless, most methods demand expert intervention, limiting their widespread adoption. By contrast, the PROSS algorithm for improving protein stability and heterologous expression levels has been successfully applied to a range of challenging enzymes and binding proteins. Here, we benchmark the application of PROSS as a stand-alone tool for protein scientists with no or limited experience in modeling. Twelve laboratories from the Protein Production and Purification Partnership in Europe (P4EU) challenged the PROSS algorithm with 14 unrelated protein targets without support from the PROSS developers. For each target, up to six designs were evaluated for expression levels and in some cases, for thermal stability and activity. In nine targets, designs exhibited increased heterologous expression levels either in prokaryotic and/or eukaryotic expression systems under experimental conditions that were tailored for each target protein. Furthermore, we observed increased thermal stability in nine of ten tested targets. In two prime examples, the human Stem Cell Factor (hSCF) and human Cadherin-Like Domain (CLD12) from the RET receptor, the wild type proteins were not expressible as soluble proteins in E. coli, yet the PROSS designs exhibited high expression levels in E. coli and HEK293 cells, respectively, and improved thermal stability. We conclude that PROSS may improve stability and expressibility in diverse cases, and that improvement typically requires target-specific expression conditions. This study demonstrates the strengths of community-wide efforts to probe the generality of new methods and recommends areas for future research to advance practically useful algorithms for protein science.

Correction: Optimizing antibody affinity and stability by the automated design of the variable light-heavy chain interfaces.

[This corrects the article DOI: 10.1371/journal.pcbi.1007207.].

Highly active rubiscos discovered by systematic interrogation of natural sequence diversity.

CO2 is converted into biomass almost solely by the enzyme rubisco. The poor carboxylation properties of plant rubiscos have led to efforts that made it the most kinetically characterized enzyme, yet these studies focused on < 5% of its natural diversity. Here, we searched for fast-carboxylating variants by systematically mining genomic and metagenomic data. Approximately 33,000 unique rubisco sequences were identified and clustered into ≈ 1,000 similarity groups. We then synthesized, purified, and biochemically tested the carboxylation rates of 143 representatives, spanning all clusters of form-II and form-II/III rubiscos. Most variants (> 100) were active in vitro, with the fastest having a turnover number of 22 ± 1 s-1 -sixfold faster than the median plant rubisco and nearly twofold faster than the fastest measured rubisco to date. Unlike rubiscos from plants and cyanobacteria, the fastest variants discovered here are homodimers and exhibit a much simpler folding and activation kinetics. Our pipeline can be utilized to explore the kinetic space of other enzymes of interest, allowing us to get a better view of the biosynthetic potential of the biosphere.

Efficient Targeted Degradation via Reversible and Irreversible Covalent PROTACs.

Proteolysis targeting chimeras (PROTACs) represent an exciting inhibitory modality with many advantages, including substoichiometric degradation of targets. Their scope, though, is still limited to date by the requirement for a sufficiently potent target binder. A solution that proved useful in tackling challenging targets is the use of electrophiles to allow irreversible binding to the target. However, such binding will negate the catalytic nature of PROTACs. Reversible covalent PROTACs potentially offer the best of both worlds. They possess the potency and selectivity associated with the formation of the covalent bond, while being able to dissociate and regenerate once the protein target is degraded. Using Bruton's tyrosine kinase (BTK) as a clinically relevant model system, we show efficient degradation by noncovalent, irreversible covalent, and reversible covalent PROTACs, with <10 nM DC50's and >85% degradation. Our data suggest that part of the degradation by our irreversible covalent PROTACs is driven by reversible binding prior to covalent bond formation, while the reversible covalent PROTACs drive degradation primarily by covalent engagement. The PROTACs showed enhanced inhibition of B cell activation compared to ibrutinib and exhibit potent degradation of BTK in patient-derived primary chronic lymphocytic leukemia cells. The most potent reversible covalent PROTAC, RC-3, exhibited enhanced selectivity toward BTK compared to noncovalent and irreversible covalent PROTACs. These compounds may pave the way for the design of covalent PROTACs for a wide variety of challenging targets.

Two closed ATP- and ADP-dependent conformations in yeast Hsp90 chaperone detected by Mn(II) EPR spectroscopic techniques.

Hsp90 plays a central role in cell homeostasis by assisting folding and maturation of a large variety of clients. It is a homo-dimer, which functions via hydrolysis of ATP-coupled to conformational changes. Hsp90's conformational cycle in the absence of cochaperones is currently postulated as apo-Hsp90 being an ensemble of "open"/"closed" conformations. Upon ATP binding, Hsp90 adopts an active ATP-bound closed conformation where the N-terminal domains, which comprise the ATP binding site, are in close contact. However, there is no consensus regarding the conformation of the ADP-bound Hsp90, which is considered important for client release. In this work, we tracked the conformational states of yeast Hsp90 at various stages of ATP hydrolysis in frozen solutions employing electron paramagnetic resonance (EPR) techniques, particularly double electron-electron resonance (DEER) distance measurements. Using rigid Gd(III) spin labels, we found the C domains to be dimerized with same distance distribution at all hydrolysis states. Then, we substituted the ATPase Mg(II) cofactor with paramagnetic Mn(II) and followed the hydrolysis state using hyperfine spectroscopy and measured the inter-N-domain distance distributions via Mn(II)-Mn(II) DEER. The point character of the Mn(II) spin label allowed us resolve 2 different closed states: The ATP-bound (prehydrolysis) characterized by a distance distribution having a maximum of 4.3 nm, which broadened and shortened, shifting the mean to 3.8 nm at the ADP-bound state (posthydrolysis). This provides experimental evidence to a second closed conformational state of Hsp90 in solution, referred to as "compact." Finally, the so-called high-energy state, trapped by addition of vanadate, was found structurally similar to the posthydrolysis state.

Optimizing antibody affinity and stability by the automated design of the variable light-heavy chain interfaces.

Antibodies developed for research and clinical applications may exhibit suboptimal stability, expressibility, or affinity. Existing optimization strategies focus on surface mutations, whereas natural affinity maturation also introduces mutations in the antibody core, simultaneously improving stability and affinity. To systematically map the mutational tolerance of an antibody variable fragment (Fv), we performed yeast display and applied deep mutational scanning to an anti-lysozyme antibody and found that many of the affinity-enhancing mutations clustered at the variable light-heavy chain interface, within the antibody core. Rosetta design combined enhancing mutations, yielding a variant with tenfold higher affinity and substantially improved stability. To make this approach broadly accessible, we developed AbLIFT, an automated web server that designs multipoint core mutations to improve contacts between specific Fv light and heavy chains (http://AbLIFT.weizmann.ac.il). We applied AbLIFT to two unrelated antibodies targeting the human antigens VEGF and QSOX1. Strikingly, the designs improved stability, affinity, and expression yields. The results provide proof-of-principle for bypassing laborious cycles of antibody engineering through automated computational affinity and stability design.

Tracking Conformational Changes in Calmodulin in vitro, in Cell Extract, and in Cells by Electron Paramagnetic Resonance Distance Measurements.

It is an open question whether the conformations of proteins sampled in dilute solutions are the same as in the cellular environment. Here we address this question by double electron-electron resonance (DEER) distance measurements with Gd(III) spin labels to probe the conformations of calmodulin (CaM) in vitro, in cell extract, and in human HeLa cells. Using the CaM mutants N53C/T110C and T34C/T117C labeled with maleimide-DOTA-Gd(III) in the N- and C-terminal domains, we observed broad and varied interdomain distance distributions. The in vitro distance distributions of apo-CaM and holo-CaM in the presence and absence of the IQ target peptide can be described by combinations of closed, open, and collapsed conformations. In cell extract, apo- and holo-CaM bind to target proteins in a similar way as apo- and holo-CaM bind to IQ peptide in vitro. In HeLa cells, however, in the presence or absence of elevated in-cell Ca2+ levels CaM unexpectedly produced more open conformations and very broad distance distributions indicative of many different interactions with in-cell components. These results show-case the importance of in-cell analyses of protein structures.

Covalent Docking Identifies a Potent and Selective MKK7 Inhibitor.

The c-Jun NH2-terminal kinase (JNK) signaling pathway is central to the cell response to stress, inflammatory signals, and toxins. While selective inhibitors are known for JNKs and for various upstream MAP3Ks, no selective inhibitor is reported for MKK7--one of two direct MAP2Ks that activate JNK. Here, using covalent virtual screening, we identify selective MKK7 covalent inhibitors. We optimized these compounds to low-micromolar inhibitors of JNK phosphorylation in cells. The crystal structure of a lead compound bound to MKK7 demonstrated that the binding mode was correctly predicted by docking. We asserted the selectivity of our inhibitors on a proteomic level and against a panel of 76 kinases, and validated an on-target effect using knockout cell lines. Lastly, we show that the inhibitors block activation of primary mouse B cells by lipopolysaccharide. These MKK7 tool compounds will enable better investigation of JNK signaling and may serve as starting points for therapeutics.

Ultrahigh specificity in a network of computationally designed protein-interaction pairs.

Protein networks in all organisms comprise homologous interacting pairs. In these networks, some proteins are specific, interacting with one or a few binding partners, whereas others are multispecific and bind a range of targets. We describe an algorithm that starts from an interacting pair and designs dozens of new pairs with diverse backbone conformations at the binding site as well as new binding orientations and sequences. Applied to a high-affinity bacterial pair, the algorithm results in 18 new ones, with cognate affinities from pico- to micromolar. Three pairs exhibit 3-5 orders of magnitude switch in specificity relative to the wild type, whereas others are multispecific, collectively forming a protein-interaction network. Crystallographic analysis confirms design accuracy, including in new backbones and polar interactions. Preorganized polar interaction networks are responsible for high specificity, thus defining design principles that can be applied to program synthetic cellular interaction networks of desired affinity and specificity.

Structure of Type-I Mycobacterium tuberculosis fatty acid synthase at 3.3 Å resolution.

Tuberculosis (TB) is a devastating and rapidly spreading disease caused by Mycobacterium tuberculosis (Mtb). Therapy requires prolonged treatment with a combination of multiple agents and interruptions in the treatment regimen result in emergence and spread of multi-drug resistant (MDR) Mtb strains. MDR Mtb poses a significant global health problem, calling for urgent development of novel drugs to combat TB. Here, we report the 3.3 Å resolution structure of the ~2 MDa type-I fatty acid synthase (FAS-I) from Mtb, determined by single particle cryo-EM. Mtb FAS-I is an essential enzymatic complex that contributes to the virulence of Mtb, and thus a prime target for anti-TB drugs. The structural information for Mtb FAS-I we have obtained enables computer-based drug discovery approaches, and the resolution achieved by cryo-EM is sufficient for elucidating inhibition mechanisms by putative small molecular weight inhibitors.