A universal deep-learning model for zinc finger design enables transcription factor reprogramming.

Cys2His2 zinc finger (ZF) domains engineered to bind specific target sequences in the genome provide an effective strategy for programmable regulation of gene expression, with many potential therapeutic applications. However, the structurally intricate engagement of ZF domains with DNA has made their design challenging. Here we describe the screening of 49 billion protein-DNA interactions and the development of a deep-learning model, ZFDesign, that solves ZF design for any genomic target. ZFDesign is a modern machine learning method that models global and target-specific differences induced by a range of library environments and specifically takes into account compatibility of neighboring fingers using a novel hierarchical transformer architecture. We demonstrate the versatility of designed ZFs as nucleases as well as activators and repressors by seamless reprogramming of human transcription factors. These factors could be used to upregulate an allele of haploinsufficiency, downregulate a gain-of-function mutation or test the consequence of regulation of a single gene as opposed to the many genes that a transcription factor would normally influence.

Basal oxidation of conserved cysteines modulates cardiac titin stiffness and dynamics.

Titin, as the main protein responsible for the passive stiffness of the sarcomere, plays a key role in diastolic function and is a determinant factor in the etiology of heart disease. Titin stiffness depends on unfolding and folding transitions of immunoglobulin-like (Ig) domains of the I-band, and recent studies have shown that oxidative modifications of cryptic cysteines belonging to these Ig domains modulate their mechanical properties in vitro. However, the relevance of this mode of titin mechanical modulation in vivo remains largely unknown. Here, we describe the high evolutionary conservation of titin mechanical cysteines and show that they are remarkably oxidized in murine cardiac tissue. Mass spectrometry analyses indicate a similar landscape of basal oxidation in murine and human myocardium. Monte Carlo simulations illustrate how disulfides and S-thiolations on these cysteines increase the dynamics of the protein at physiological forces, while enabling load- and isoform-dependent regulation of titin stiffness. Our results demonstrate the role of conserved cysteines in the modulation of titin mechanical properties in vivo and point to potential redox-based pathomechanisms in heart disease.

Engineered dual selection for directed evolution of SpCas9 PAM specificity.

The widely used Streptococcus pyogenes Cas9 (SpCas9) nuclease derives its DNA targeting specificity from protein-DNA contacts with protospacer adjacent motif (PAM) sequences, in addition to base-pairing interactions between its guide RNA and target DNA. Previous reports have established that the PAM specificity of SpCas9 can be altered via positive selection procedures for directed evolution or other protein engineering strategies. Here we exploit in vivo directed evolution systems that incorporate simultaneous positive and negative selection to evolve SpCas9 variants with commensurate or improved activity on NAG PAMs relative to wild type and reduced activity on NGG PAMs, particularly YGG PAMs. We also show that the PAM preferences of available evolutionary intermediates effectively determine whether similar counterselection PAMs elicit different selection stringencies, and demonstrate that negative selection can be specifically increased in a yeast selection system through the fusion of compensatory zinc fingers to SpCas9.

Dissecting the Structural and Chemical Determinants of the "Open-to-Closed" Motion in the Mannosyltransferase PimA from Mycobacteria.

The phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential peripheral membrane glycosyltransferase that initiates the biosynthetic pathway of phosphatidyl-myo-inositol mannosides (PIMs), key structural elements and virulence factors of Mycobacterium tuberculosis. PimA undergoes functionally important conformational changes, including (i) α-helix-to-ß-strand and ß-strand-to-α-helix transitions and (ii) an "open-to-closed" motion between the two Rossmann-fold domains, a conformational change that is necessary to generate a catalytically competent active site. In previous work, we established that GDP-Man and GDP stabilize the enzyme and facilitate the switch to a more compact active state. To determine the structural contribution of the mannose ring in such an activation mechanism, we analyzed a series of chemical derivatives, including mannose phosphate (Man-P) and mannose pyrophosphate-ribose (Man-PP-RIB), and additional GDP derivatives, such as pyrophosphate ribose (PP-RIB) and GMP, by the combined use of X-ray crystallography, limited proteolysis, circular dichroism, isothermal titration calorimetry, and small angle X-ray scattering methods. Although the ß-phosphate is present, we found that the mannose ring, covalently attached to neither phosphate (Man-P) nor PP-RIB (Man-PP-RIB), does promote the switch to the active compact form of the enzyme. Therefore, the nucleotide moiety of GDP-Man, and not the sugar ring, facilitates the "open-to-closed" motion, with the ß-phosphate group providing the high-affinity binding to PimA. Altogether, the experimental data contribute to a better understanding of the structural determinants involved in the "open-to-closed" motion not only observed in PimA but also visualized and/or predicted in other glycosyltransfeases. In addition, the experimental data might prove to be useful for the discovery and/or development of PimA and/or glycosyltransferase inhibitors.

The geometric influence on the Cys2His2 zinc finger domain and functional plasticity.

The Cys2His2 zinc finger is the most common DNA-binding domain expanding in metazoans since the fungi human split. A proposed catalyst for this expansion is an arms race to silence transposable elements yet it remains poorly understood how this domain is able to evolve the required specificities. Likewise, models of its DNA binding specificity remain error prone due to a lack of understanding of how adjacent fingers influence each other's binding specificity. Here, we use a synthetic approach to exhaustively investigate binding geometry, one of the dominant influences on adjacent finger function. By screening over 28 billion protein-DNA interactions in various geometric contexts we find the plasticity of the most common natural geometry enables more functional amino acid combinations across all targets. Further, residues that define this geometry are enriched in genomes where zinc fingers are prevalent and specificity transitions would be limited in alternative geometries. Finally, these results demonstrate an exhaustive synthetic screen can produce an accurate model of domain function while providing mechanistic insight that may have assisted in the domains expansion.

Conformational entropy of a single peptide controlled under force governs protease recognition and catalysis.

An immense repertoire of protein chemical modifications catalyzed by enzymes is available as proteomics data. Quantifying the impact of the conformational dynamics of the modified peptide remains challenging to understand the decisive kinetics and amino acid sequence specificity of these enzymatic reactions in vivo, because the target peptide must be disordered to accommodate the specific enzyme-binding site. Here, we were able to control the conformation of a single-molecule peptide chain by applying mechanical force to activate and monitor its specific cleavage by a model protease. We found that the conformational entropy impacts the reaction in two distinct ways. First, the flexibility and accessibility of the substrate peptide greatly increase upon mechanical unfolding. Second, the conformational sampling of the disordered peptide drives the specific recognition, revealing force-dependent reaction kinetics. These results support a mechanism of peptide recognition based on conformational selection from an ensemble that we were able to quantify with a torsional free-energy model. Our approach can be used to predict how entropy affects site-specific modifications of proteins and prompts conformational and mechanical selectivity.

The force-dependent mechanism of DnaK-mediated mechanical folding.

It is well established that chaperones modulate the protein folding free-energy landscape. However, the molecular determinants underlying chaperone-mediated mechanical folding remain largely elusive, primarily because the force-extended unfolded conformation fundamentally differs from that characterized in biochemistry experiments. We use single-molecule force-clamp spectroscopy, combined with molecular dynamics simulations, to study the effect that the Hsp70 system has on the mechanical folding of three mechanically stiff model proteins. Our results demonstrate that, when working independently, DnaJ (Hsp40) and DnaK (Hsp70) work as holdases, blocking refolding by binding to distinct substrate conformations. Whereas DnaK binds to molten globule-like forms, DnaJ recognizes a cryptic sequence in the extended state in an unanticipated force-dependent manner. By contrast, the synergetic coupling of the Hsp70 system exhibits a marked foldase behavior. Our results offer unprecedented molecular and kinetic insights into the mechanisms by which mechanical force finely regulates chaperone binding, directly affecting protein elasticity.

Disulfide isomerization reactions in titin immunoglobulin domains enable a mode of protein elasticity.

The response of titin to mechanical forces is a major determinant of the function of the heart. When placed under a pulling force, the unstructured regions of titin uncoil while its immunoglobulin (Ig) domains unfold and extend. Using single-molecule atomic force microscopy, we show that disulfide isomerization reactions within Ig domains enable a third mechanism of titin elasticity. Oxidation of Ig domains leads to non-canonical disulfide bonds that stiffen titin while enabling force-triggered isomerization reactions to more extended states of the domains. Using sequence and structural analyses, we show that 21% of titin's I-band Ig domains contain a conserved cysteine triad that can engage in disulfide isomerization reactions. We propose that imbalance of the redox status of myocytes can have immediate consequences for the mechanical properties of the sarcomere via alterations of the oxidation state of titin domains.

Single-molecule Force Spectroscopy Predicts a Misfolded, Domain-swapped Conformation in human γD-Crystallin Protein.

Cataract is a protein misfolding disease where the size of the aggregate is directly related to the severity of the disorder. However, the molecular mechanisms that trigger the onset of aggregation remain unknown. Here we use a combination of protein engineering techniques and single-molecule force spectroscopy using atomic force microscopy to study the individual unfolding pathways of the human γD-crystallin, a multidomain protein that must remain correctly folded during the entire lifetime to guarantee lens transparency. When stretching individual polyproteins containing two neighboring HγD-crystallin monomers, we captured an anomalous misfolded conformation in which the ß1 and ß2 strands of the N terminus domain of two adjacent monomers swap. This experimentally elusive domain-swapped conformation is likely to be responsible for the increase in molecular aggregation that we measure in vitro. Our results demonstrate the power of force spectroscopy at capturing rare misfolded conformations with potential implications for the understanding of the molecular onset of protein aggregation.

Secondary structure reshuffling modulates glycosyltransferase function at the membrane.

Secondary structure refolding is a key event in biology as it modulates the conformation of many proteins in the cell, generating functional or aberrant states. The crystal structures of mannosyltransferase PimA reveal an exceptional flexibility of the protein along the catalytic cycle, including ß-strand-to-α-helix and α-helix-to-ß-strand transitions. These structural changes modulate catalysis and are promoted by interactions of the protein with anionic phospholipids in the membrane.

A novel Plasmodium-specific prodomain fold regulates the malaria drug target SUB1 subtilase.

The Plasmodium subtilase SUB1 plays a pivotal role during the egress of malaria parasites from host hepatocytes and erythrocytes. Here we report the crystal structure of full-length SUB1 from the human-infecting parasite Plasmodium vivax, revealing a bacterial-like catalytic domain in complex with a Plasmodium-specific prodomain. The latter displays a novel architecture with an amino-terminal insertion that functions as a 'belt', embracing the catalytic domain to further stabilize the quaternary structure of the pre-protease, and undergoes calcium-dependent autoprocessing during subsequent activation. Although dispensable for recombinant enzymatic activity, the SUB1 'belt' could not be deleted in Plasmodium berghei, suggesting an essential role of this domain for parasite development in vivo. The SUB1 structure not only provides a valuable platform to develop new anti-malarial candidates against this promising drug target, but also defines the Plasmodium-specific 'belt' domain as a key calcium-dependent regulator of SUB1 during parasite egress from host cells.

S-glutathionylation of cryptic cysteines enhances titin elasticity by blocking protein folding.

The giant elastic protein titin is a determinant factor in how much blood fills the left ventricle during diastole and thus in the etiology of heart disease. Titin has been identified as a target of S-glutathionylation, an end product of the nitric-oxide-signaling cascade that increases cardiac muscle elasticity. However, it is unknown how S-glutathionylation may regulate the elasticity of titin and cardiac tissue. Here, we show that mechanical unfolding of titin immunoglobulin (Ig) domains exposes buried cysteine residues, which then can be S-glutathionylated. S-glutathionylation of cryptic cysteines greatly decreases the mechanical stability of the parent Ig domain as well as its ability to fold. Both effects favor a more extensible state of titin. Furthermore, we demonstrate that S-glutathionylation of cryptic cysteines in titin mediates mechanochemical modulation of the elasticity of human cardiomyocytes. We propose that posttranslational modification of cryptic residues is a general mechanism to regulate tissue elasticity.

Structure-function relationships of membrane-associated GT-B glycosyltransferases.

Membrane-associated GT-B glycosyltransferases (GTs) comprise a large family of enzymes that catalyze the transfer of a sugar moiety from nucleotide-sugar donors to a wide range of membrane-associated acceptor substrates, mostly in the form of lipids and proteins. As a consequence, they generate a significant and diverse amount of glycoconjugates in biological membranes, which are particularly important in cell-cell, cell-matrix and host-pathogen recognition events. Membrane-associated GT-B enzymes display two "Rossmann-fold" domains separated by a deep cleft that includes the catalytic center. They associate permanently or temporarily to the phospholipid bilayer by a combination of hydrophobic and electrostatic interactions. They have the remarkable property to access both hydrophobic and hydrophilic substrates that reside within chemically distinct environments catalyzing their enzymatic transformations in an efficient manner. Here, we discuss the considerable progress that has been made in recent years in understanding the molecular mechanism that governs substrate and membrane recognition, and the impact of the conformational transitions undergone by these GTs during the catalytic cycle.

Conformational plasticity of the essential membrane-associated mannosyltransferase PimA from mycobacteria.

Phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential glycosyltransferase (GT) that initiates the biosynthetic pathway of phosphatidyl-myo-inositol mannosides, lipomannan, and lipoarabinomannan, which are key glycolipids/lipoglycans of the mycobacterial cell envelope. PimA belongs to a large family of peripheral membrane-associated GTs for which the understanding of the molecular mechanism and conformational changes that govern substrate/membrane recognition and catalysis remains a major challenge. Here we used single molecule force spectroscopy techniques to study the mechanical and conformational properties of PimA. In our studies, we engineered a polyprotein containing PimA flanked by four copies of the well characterized I27 protein, which provides an unambiguous mechanical fingerprint. We found that PimA exhibits weak mechanical stability albeit displaying ß-sheet topology expected to unfold at much higher forces. Notably, PimA unfolds following heterogeneous multiple step mechanical unfolding pathways at low force akin to molten globule states. Interestingly, the ab initio low resolution envelopes obtained from small angle x-ray scattering of the unliganded PimA and the PimA·GDP complexed forms clearly demonstrate that not only the "open" and "closed" conformations of the GT-B enzyme are largely present in solution, but in addition, PimA experiences remarkable flexibility that undoubtedly corresponds to the N-terminal "Rossmann fold" domain, which has been proved to participate in protein-membrane interactions. Based on these results and on our previous experimental data, we propose a model wherein the conformational transitions are important for the mannosyltransferase to interact with the donor and acceptor substrates/membrane.

In Silico screening on the three-dimensional model of the Plasmodium vivax SUB1 protease leads to the validation of a novel anti-parasite compound.

Widespread drug resistance calls for the urgent development of new antimalarials that target novel steps in the life cycle of Plasmodium falciparum and Plasmodium vivax. The essential subtilisin-like serine protease SUB1 of Plasmodium merozoites plays a dual role in egress from and invasion into host erythrocytes. It belongs to a new generation of attractive drug targets against which specific potent inhibitors are actively searched. We characterize here the P. vivax SUB1 enzyme and show that it displays a typical auto-processing pattern and apical localization in P. vivax merozoites. To search for small PvSUB1 inhibitors, we took advantage of the similarity of SUB1 with bacterial subtilisins and generated P. vivax SUB1 three-dimensional models. The structure-based virtual screening of a large commercial chemical compounds library identified 306 virtual best hits, of which 37 were experimentally confirmed inhibitors and 5 had Ki values of <50 µM for PvSUB1. Interestingly, they belong to different chemical families. The most promising competitive inhibitor of PvSUB1 (compound 2) was equally active on PfSUB1 and displayed anti-P. falciparum and Plasmodium berghei activity in vitro and in vivo, respectively. Compound 2 inhibited the endogenous PfSUB1 as illustrated by the inhibited maturation of its natural substrate PfSERA5 and inhibited parasite egress and subsequent erythrocyte invasion. These data indicate that the strategy of in silico screening of three-dimensional models to select for virtual inhibitors combined with stringent biological validation successfully identified several inhibitors of the PvSUB1 enzyme. The most promising hit proved to be a potent cross-inhibitor of PlasmodiumSUB1, laying the groundwork for the development of a globally active small compound antimalarial.

Elasticity, structure, and relaxation of extended proteins under force.

Force spectroscopies have emerged as a powerful and unprecedented tool to study and manipulate biomolecules directly at a molecular level. Usually, protein and DNA behavior under force is described within the framework of the worm-like chain (WLC) model for polymer elasticity. Although it has been surprisingly successful for the interpretation of experimental data, especially at high forces, the WLC model lacks structural and dynamical molecular details associated with protein relaxation under force that are key to the understanding of how force affects protein flexibility and reactivity. We use molecular dynamics simulations of ubiquitin to provide a deeper understanding of protein relaxation under force. We find that the WLC model successfully describes the simulations of ubiquitin, especially at higher forces, and we show how protein flexibility and persistence length, probed in the force regime of the experiments, are related to how specific classes of backbone dihedral angles respond to applied force. Although the WLC model is an average, backbone model, we show how the protein side chains affect the persistence length. Finally, we find that the diffusion coefficient of the protein's end-to-end distance is on the order of 10(8) nm(2)/s, is position and side-chain dependent, but is independent of the length and independent of the applied force, in contrast with other descriptions.

Biological and structural characterization of the Mycobacterium smegmatis nitroreductase NfnB, and its role in benzothiazinone resistance.

Tuberculosis is still a leading cause of death in developing countries, for which there is an urgent need for new pharmacological agents. The synthesis of the novel antimycobacterial drug class of benzothiazinones (BTZs) and the identification of their cellular target as DprE1 (Rv3790), a component of the decaprenylphosphoryl-ß-d-ribose 2'-epimerase complex, have been reported recently. Here, we describe the identification and characterization of a novel resistance mechanism to BTZ in Mycobacterium smegmatis. The overexpression of the nitroreductase NfnB leads to the inactivation of the drug by reduction of a critical nitro-group to an amino-group. The direct involvement of NfnB in the inactivation of the lead compound BTZ043 was demonstrated by enzymology, microbiological assays and gene knockout experiments. We also report the crystal structure of NfnB in complex with the essential cofactor flavin mononucleotide, and show that a common amino acid stretch between NfnB and DprE1 is likely to be essential for the interaction with BTZ. We performed docking analysis of NfnB-BTZ in order to understand their interaction and the mechanism of nitroreduction. Although Mycobacterium tuberculosis seems to lack nitroreductases able to inactivate these drugs, our findings are valuable for the design of new BTZ molecules, which may be more effective in vivo.

Comparative evaluation of 3D virtual ligand screening methods: impact of the molecular alignment on enrichment.

In the early stage of drug discovery programs, when the structure of a complex involving a target and a small molecule is available, structure-based virtual ligand screening methods are generally preferred. However, ligand-based strategies like shape-similarity search methods can also be applied. Shape-similarity search methods consist in exploring a pseudo-binding-site derived from the known small molecule used as a reference. Several of these methods use conformational sampling algorithms which are also shared by corresponding docking methods: for example Surflex-dock/Surflex-sim, FlexX/FlexS, ICM, and OMEGA-FRED/OMEGA-ROCS. Using 11 systems issued from the challenging "own" subsets of the Directory of Useful Decoys (DUD-own), we evaluated and compared the performance of the above-cited programs in terms of molecular alignment accuracy, enrichment in active compounds, and enrichment in different chemotypes (scaffold-hopping). Since molecular alignment is a crucial aspect of performance for the different methods, we have assessed its impact on enrichment. We have also illustrated the paradox of retrieving active compounds with good scores even if they are inaccurately positioned. Finally, we have highlighted possible positive aspects of using shape-based approaches in drug-discovery protocols when the structure of the target in complex with a small molecule is known.

Substrate-induced conformational changes in the essential peripheral membrane-associated mannosyltransferase PimA from mycobacteria: implications for catalysis.

Phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential glycosyltransferase (GT) involved in the biosynthesis of phosphatidyl-myo-inositol mannosides (PIMs), which are key components of the mycobacterial cell envelope. PimA is the paradigm of a large family of peripheral membrane-binding GTs for which the molecular mechanism of substrate/membrane recognition and catalysis is still unknown. Strong evidence is provided showing that PimA undergoes significant conformational changes upon substrate binding. Specifically, the binding of the donor GDP-Man triggered an important interdomain rearrangement that stabilized the enzyme and generated the binding site for the acceptor substrate, phosphatidyl-myo-inositol (PI). The interaction of PimA with the beta-phosphate of GDP-Man was essential for this conformational change to occur. In contrast, binding of PI had the opposite effect, inducing the formation of a more relaxed complex with PimA. Interestingly, GDP-Man stabilized and PI destabilized PimA by a similar enthalpic amount, suggesting that they formed or disrupted an equivalent number of interactions within the PimA complexes. Furthermore, molecular docking and site-directed mutagenesis experiments provided novel insights into the architecture of the myo-inositol 1-phosphate binding site and the involvement of an essential amphiphatic alpha-helix in membrane binding. Altogether, our experimental data support a model wherein the flexibility and conformational transitions confer the adaptability of PimA to the donor and acceptor substrates, which seems to be of importance during catalysis. The proposed mechanism has implications for the comprehension of the peripheral membrane-binding GTs at the molecular level.

Molecular recognition and interfacial catalysis by the essential phosphatidylinositol mannosyltransferase PimA from mycobacteria.

Mycobacterial phosphatidylinositol mannosides (PIMs) and metabolically derived cell wall lipoglycans play important roles in host-pathogen interactions, but their biosynthetic pathways are poorly understood. Here we focus on Mycobacterium smegmatis PimA, an essential enzyme responsible for the initial mannosylation of phosphatidylinositol. The structure of PimA in complex with GDP-mannose shows the two-domain organization and the catalytic machinery typical of GT-B glycosyltransferases. PimA is an amphitrophic enzyme that binds mono-disperse phosphatidylinositol, but its transferase activity is stimulated by high concentrations of non-substrate anionic surfactants, indicating that the early stages of PIM biosynthesis involve lipid-water interfacial catalysis. Based on structural, calorimetric, and mutagenesis studies, we propose a model wherein PimA attaches to the membrane through its N-terminal domain, and this association leads to enzyme activation. Our results reveal a novel mode of phosphatidylinositol recognition and provide a template for the development of potential antimycobacterial compounds.