DMSO Solubility Assessment for Fragment-Based Screening.

In this paper, we report comprehensive experimental and chemoinformatics analyses of the solubility of small organic molecules ("fragments") in dimethyl sulfoxide (DMSO) in the context of their ability to be tested in screening experiments. Here, DMSO solubility of 939 fragments has been measured experimentally using an NMR technique. A Support Vector Classification model was built on the obtained data using the ISIDA fragment descriptors. The analysis revealed 34 outliers: experimental issues were retrospectively identified for 28 of them. The updated model performs well in 5-fold cross-validation (balanced accuracy = 0.78). The datasets are available on the Zenodo platform (DOI:10.5281/zenodo.4767511) and the model is available on the website of the Laboratory of Chemoinformatics.

The One-carbon Carrier Methylofuran from Methylobacterium extorquens AM1 Contains a Large Number of α- and γ-Linked Glutamic Acid Residues.

Methylobacterium extorquens AM1 uses dedicated cofactors for one-carbon unit conversion. Based on the sequence identities of enzymes and activity determinations, a methanofuran analog was proposed to be involved in formaldehyde oxidation in Alphaproteobacteria. Here, we report the structure of the cofactor, which we termed methylofuran. Using an in vitro enzyme assay and LC-MS, methylofuran was identified in cell extracts and further purified. From the exact mass and MS-MS fragmentation pattern, the structure of the cofactor was determined to consist of a polyglutamic acid side chain linked to a core structure similar to the one present in archaeal methanofuran variants. NMR analyses showed that the core structure contains a furan ring. However, instead of the tyramine moiety that is present in methanofuran cofactors, a tyrosine residue is present in methylofuran, which was further confirmed by MS through the incorporation of a (13)C-labeled precursor. Methylofuran was present as a mixture of different species with varying numbers of glutamic acid residues in the side chain ranging from 12 to 24. Notably, the glutamic acid residues were not solely γ-linked, as is the case for all known methanofurans, but were identified by NMR as a mixture of α- and γ-linked amino acids. Considering the unusual peptide chain, the elucidation of the structure presented here sets the basis for further research on this cofactor, which is probably the largest cofactor known so far.

Local and Global Dynamics in Klebsiella pneumoniae Outer Membrane Protein a in Lipid Bilayers Probed at Atomic Resolution.

The role of membrane proteins in cellular mechanism strongly depends on their dynamics, and solid-state magic-angle spinning (MAS) nuclear magnetic resonance (NMR) is a unique method to exhaustively characterize motions of proteins in a lipid environment. Herein, we make use of advances in 1H-detected MAS NMR to describe the dynamics of the membrane domain of the Outer membrane protein A of Klebsiella pneumoniae (KpOmpA). By measuring 1H-15N dipolar-coupling as well as 15N R1 and R1ρ relaxation rates at fast (60 kHz) MAS and high magnetic field (1 GHz), we were able to describe the motions of the residues of the ß-barrel as a collective rocking of low amplitude and of hundreds of nanoseconds time scale. Residual local motions at the edges of the strands, underscored by enhanced 15N R1ρ relaxation rates, report on the mobility of the connected loops. In agreement with MAS NMR data, proteolysis experiments performed on the full length KpOmpA as well as on its membrane domain, reconstituted in liposomes or in detergent micelles, revealed in all cases the existence of a unique trypsin cleavage site within the membrane domain (out of 16 potential Lys and Arg sites). This site is located in the extracellular loop L3, showing that it is highly accessible to protein-protein interactions. KpOmpA is involved in cell-cell recognition, for adhesion and immune response mechanisms. The L3 region may therefore play a key role in pathogenicity.

The C-terminal region of the transcriptional regulator THAP11 forms a parallel coiled-coil domain involved in protein dimerization.

Thanatos associated protein 11 (THAP11) is a cell cycle and cell growth regulator differentially expressed in cancer cells. THAP11 belongs to a distinct family of transcription factors recognizing specific DNA sequences via an atypical zinc finger motif and regulating diverse cellular processes. Outside the extensively characterized DNA-binding domain, THAP proteins vary in size and predicted domains, for which structural data are still lacking. We report here the crystal structure of the C-terminal region of human THAP11 protein, providing the first 3D structure of a coiled-coil motif from a THAP family member. We further investigate the stability, dynamics and oligomeric properties of the determined structure combining molecular dynamics simulations and biophysical experiments. Our results show that the C-ter region of THAP11 forms a left-handed parallel homo-dimeric coiled-coil structure possessing several unusual features.

Mutation of a pH-modulating residue in a GH51 α-l-arabinofuranosidase leads to a severe reduction of the secondary hydrolysis of transfuranosylation products.

BACKGROUND: The development of enzyme-mediated glycosynthesis using glycoside hydrolases is still an inexact science, because the underlying molecular determinants of transglycosylation are not well understood. In the framework of this challenge, this study focused on the family GH51 α-l-arabinofuranosidase from Thermobacillus xylanilyticus, with the aim to understand why the mutation of position 344 provokes a significant modification of the transglycosylation/hydrolysis partition. METHODS: Detailed kinetic analysis (kcat, KM, pKa determination and time-course NMR kinetics) and saturation transfer difference nuclear magnetic resonance spectroscopy was employed to determine the synthetic and hydrolytic ability modification induced by the redundant N344 mutation disclosed in libraries from directed evolution. RESULTS: The mutants N344P and N344Y displayed crippled hydrolytic abilities, and thus procured improved transglycosylation yields. This behavior was correlated with an increased pKa of the catalytic nucleophile (E298), the pKa of the acid/base catalyst remaining unaffected. Finally, mutations at position 344 provoked a pH-dependent product inhibition phenomenon, which is likely to be the result of a significant modification of the proton sharing network in the mutants. CONCLUSIONS AND GENERAL SIGNIFICANCE: Using a combination of biochemical and biophysical methods, we have studied TxAbf-N344 mutants, thus revealing some fundamental details concerning pH modulation. Although these results concern a GH51 α-l-arabinofuranosidase, it is likely that the general lessons that can be drawn from them will be applicable to other glycoside hydrolases. Moreover, the effects of mutations at position 344 on the transglycosylation/hydrolysis partition provide clues as to how TxAbf can be further engineered to obtain an efficient transfuranosidase.

Combined experimental and simulation studies suggest a revised mode of action of the anti-Alzheimer disease drug NQ-Trp.

Inhibition of the aggregation of the monomeric peptide ß-amyloid (Aß) into oligomers is a widely studied therapeutic approach in Alzheimer's disease (AD). Many small molecules have been reported to work in this way, including 1,4-naphthoquinon-2-yl-L-tryptophan (NQ-Trp). NQ-Trp has been reported to inhibit aggregation, to rescue cells from Aß toxicity, and showed complete phenotypic recovery in an in vivo AD model. In this work we investigated its molecular mechanism by using a combined approach of experimental and theoretical studies, and obtained converging results. NQ-Trp is a relatively weak inhibitor and the fluorescence data obtained by employing the fluorophore widely used to monitor aggregation into fibrils can be misinterpreted due to the inner filter effect. Simulations and NMR experiments showed that NQ-Trp has no specific "binding site"-type interaction with mono- and dimeric Aß, which could explain its low inhibitory efficiency. This suggests that the reported anti-AD activity of NQ-Trp-type molecules in in vivo models has to involve another mechanism. This study has revealed the potential pitfalls in the development of aggregation inhibitors for amyloidogenic peptides, which are of general interest for all the molecules studied in the context of inhibiting the formation of toxic aggregates.

Two classes of cholesterol binding sites for the ß2AR revealed by thermostability and NMR.

Cholesterol binding to G protein-coupled receptors (GPCRs) and modulation of their activities in membranes is a fundamental issue for understanding their function. Despite the identification of cholesterol binding sites in high-resolution x-ray structures of the ?2 adrenergic receptor (ß2AR) and other GPCRs, the binding affinity of cholesterol for this receptor and exchange rates between the free and bound cholesterol remain unknown. In this study we report the existence of two classes of cholesterol binding sites in ß2AR. By analyzing the ß2AR unfolding temperature in lipidic cubic phase (LCP) as a function of cholesterol concentration we observed high-affinity cooperative binding of cholesterol with sub-nM affinity constant. In contrast, saturation transfer difference (STD) NMR experiments revealed the existence of a second class of cholesterol binding sites, in fast exchange on the STD NMR timescale. Titration of the STD signal as a function of cholesterol concentration provided a lower limit of 100 mM for their dissociation constant. However, these binding sites are specific for both cholesterol and ß2AR, as shown with control experiments using ergosterol and a control membrane protein (KpOmpA). We postulate that this specificity is mediated by the high-affinity bound cholesterol molecules and propose the formation of transient cholesterol clusters around the high-affinity binding sites.

Cord factor (trehalose 6,6'-dimycolate) forms fully stable and non-permeable lipid bilayers required for a functional outer membrane.

Cord factor (trehalose 6,6'-dimycolate, TDM) is the major lipid in the outer membrane of Corynebacteria and Mycobacteria. Although its role is well recognized in the immune response phenomena, its membrane biophysical properties remained largely unexplored and TDM has often been described as a detergent. We purified the main components of the outer membrane from Corynebacterium glutamicum and analyzed their membrane forming properties. In mixture with endogenous cardiolipin, but not alone, the spontaneous hydration of TDM produces liposomes. As a pure component, TDM formed vesicles only by the detergent dialysis method. Perdeuterated cardiolipin-TDM mixtures were shown by deuterium nuclear magnetic resonance (NMR) to exhibit a gel to liquid crystalline phase transition over a 273-295K temperature range, for cells grown at 303K, and thus to be in a liquid crystalline state at physiological temperature. Molecular dynamics simulations of hydrated TDM bilayers provided the trehalose average orientation and conformation, the chain order parameters, the area per lipid and the bilayer thickness which was confirmed by electron microscopy. Finally the Porin A-Porin H ion channel from the Corynebacterial outer membrane was reconstituted in TDM liposomes. With properly mycoloylated proteins, it manifested the typical voltage dependent ion channel properties of an outer membrane porin.

Kinetic solubility: Experimental and machine-learning modeling perspectives.

Kinetic aqueous or buffer solubility is important parameter measuring suitability of compounds for high throughput assays in early drug discovery while thermodynamic solubility is reserved for later stages of drug discovery and development. Kinetic solubility is also considered to have low inter-laboratory reproducibility because of its sensitivity to protocol parameters [1]. Presumably, this is why little efforts have been put to build QSPR models for kinetic in comparison to thermodynamic aqueous solubility. Here, we investigate the reproducibility and modelability of kinetic solubility assays. We first analyzed the relationship between kinetic and thermodynamic solubility data, and then examined the consistency of data from different kinetic assays. In this contribution, we report differences between kinetic and thermodynamic solubility data that are consistent with those reported by others [1, 2] and good agreement between data from different kinetic solubility campaigns in contrast to general expectations. The latter is confirmed by achieving high performing QSPR models trained on merged kinetic solubility datasets. The poor performance of QSPR model trained on thermodynamic solubility when applied to kinetic solubility dataset reinforces the conclusion that kinetic and thermodynamic solubilities do not correlate: one cannot be used as an ersatz for the other. This encourages for building predictive models for kinetic solubility. The kinetic solubility QSPR model developed in this study is freely accessible through the Predictor web service of the Laboratory of Chemoinformatics (https://chematlas.chimie.unistra.fr/cgi-bin/predictor2.cgi).

Dynamics of Klebsiella pneumoniae OmpA transmembrane domain: the four extracellular loops display restricted motion behavior in micelles and in lipid bilayers.

The transmembrane domain of Klebsiella pneumoniae OmpA (KpOmpA) possesses four long extracellular loops that exhibit substantial sequence variability throughout OmpA homologs in Enterobacteria, in comparison with the highly conserved membrane-embedded ß-barrel core. These loops are responsible for the immunological properties of the protein, including cellular and humoral recognition. In addition to key features revealed by structural elucidation of the KpOmpA transmembrane domain in detergent micelles, studies of protein dynamics provide insight into its function and/or mechanism of action. We have investigated the dynamics of KpOmpA in a lipid bilayer, using magic angle spinning solid-state NMR. The dynamics of the ß-barrel and loop regions were probed by the spin-lattice relaxation times of the C(α) and C(ß) atoms of the serine and threonine residues, and by cross-polarization dynamics. The ß-barrel core of the protein is rigid; the C-terminal halves of two of the four extracellular loops (L1 and L3), which are particularly long in KpOmpA, are highly mobile. The other two loops (L2 and L4), which are very similar to their homologs in Escherichia coli OmpA, and the N-terminal halves of L1 and L3 exhibit more restricted motions. We suggest a correlation between the sequence variability and the dynamics of certain loop regions, which accounts for their respective contributions to the structural and immunological properties of the protein.

Hydrogen bonding of cholesterol in the lipidic cubic phase.

The addition of cholesterol to the monoolein-based lipidic cubic phase (LCP) has been instrumental in obtaining high-resolution crystal structures of several G protein-coupled receptors. Here, we report the use of high-resolution magic angle spinning NMR spectroscopy to record and assign the isotropic (13)C chemical shifts of cholesterol in lipidic lamellar and cubic phases at different hydration levels with monoolein and chain-deuterated DMPC as host lipids. The hydrogen-bonding patterns of cholesterol in these phases were determined from the NMR data by quantum chemical calculations. The results are consistent with the normal orientation of cholesterol in lipid bilayers and with the cholesterol hydroxyl group located at the hydrophobic/hydrophilic interface. The (13)C chemical shifts of cholesterol are mostly affected by the host lipid identity with little or no dependency on the hydration (20% vs 40%) or the phase identity (lamellar vs LCP). In chain-deuterated DMPC bilayers, the hydroxyl group of cholesterol forms most of its hydrogen bonds with water, while in monoolein bilayers it predominately interacts with monoolein. Such differences in the hydrogen-bonding network of cholesterol may have implications for the design of experiments in monoolein-based LCP.

Sodium is a negative allosteric regulator of the ghrelin receptor.

The functional properties of G protein-coupled receptors (GPCRs) are intimately associated with the different components in their cellular environment. Among them, sodium ions have been proposed to play a substantial role as endogenous allosteric modulators of GPCR-mediated signaling. However, this sodium effect and the underlying mechanisms are still unclear for most GPCRs. Here, we identified sodium as a negative allosteric modulator of the ghrelin receptor GHSR (growth hormone secretagogue receptor). Combining 23Na-nuclear magnetic resonance (NMR), molecular dynamics, and mutagenesis, we provide evidence that, in GHSR, sodium binds to the allosteric site conserved in class A GPCRs. We further leveraged spectroscopic and functional assays to show that sodium binding shifts the conformational equilibrium toward the GHSR-inactive ensemble, thereby decreasing basal and agonist-induced receptor-catalyzed G protein activation. All together, these data point to sodium as an allosteric modulator of GHSR, making this ion an integral component of the ghrelin signaling machinery.

Functional expression of the PorAH channel from Corynebacterium glutamicum in cell-free expression systems: implications for the role of the naturally occurring mycolic acid modification.

PorA and PorH are two small membrane proteins from the outer membrane of Corynebacterium glutamicum, which have been shown to form heteromeric ion channels and to be post-translationally modified by mycolic acids. Any structural details of the channel could not be analyzed so far due to tremendous difficulties in the production of sufficient amounts of protein samples. Cell-free (CF) expression is a new and remarkably successful strategy for the production of membrane proteins for which toxicity, membrane targeting, and degradation are key issues. In addition, reaction conditions can easily be modified to modulate the quality of synthesized protein samples. We developed an efficient CF expression strategy to produce the channel subunits devoid of post-translational modifications. (15)N-labeled PorA and PorH samples were furthermore characterized by NMR and gave well resolved spectra, opening the way for structural studies. The comparison of ion channel activities of CF-expressed proteins with channels isolated from C. glutamicum gave clear insights on the influence of the mycolic acid modification of the two subunits on their functional properties.

Structural determinants of specific DNA-recognition by the THAP zinc finger.

Human THAP1 is the prototype of a large family of cellular factors sharing an original THAP zinc-finger motif responsible for DNA binding. Human THAP1 regulates endothelial cell proliferation and G1/S cell-cycle progression, through modulation of pRb/E2F cell-cycle target genes including rrm1. Recently, mutations in THAP1 have been found to cause DYT6 primary torsion dystonia, a human neurological disease. We report here the first 3D structure of the complex formed by the DNA-binding domain of THAP1 and its specific DNA target (THABS) found within the rrm1 target gene. The THAP zinc finger uses its double-stranded beta-sheet to fill the DNA major groove and provides a unique combination of contacts from the beta-sheet, the N-terminal tail and surrounding loops toward the five invariant base pairs of the THABS sequence. Our studies reveal unprecedented insights into the specific DNA recognition mechanisms within this large family of proteins controlling cell proliferation, cell cycle and pluripotency.

Mycobacterial resistance to zinc poisoning requires assembly of P-ATPase-containing membrane metal efflux platforms.

The human pathogen Mycobacterium tuberculosis requires a P1B-ATPase metal exporter, CtpC (Rv3270), for resistance to zinc poisoning. Here, we show that zinc resistance also depends on a chaperone-like protein, PacL1 (Rv3269). PacL1 contains a transmembrane domain, a cytoplasmic region with glutamine/alanine repeats and a C-terminal metal-binding motif (MBM). PacL1 binds Zn2+, but the MBM is required only at high zinc concentrations. PacL1 co-localizes with CtpC in dynamic foci in the mycobacterial plasma membrane, and the two proteins form high molecular weight complexes. Foci formation does not require flotillin nor the PacL1 MBM. However, deletion of the PacL1 Glu/Ala repeats leads to loss of CtpC and sensitivity to zinc. Genes pacL1 and ctpC appear to be in the same operon, and homologous gene pairs are found in the genomes of other bacteria. Furthermore, PacL1 colocalizes and functions redundantly with other PacL orthologs in M. tuberculosis. Overall, our results indicate that PacL proteins may act as scaffolds that assemble P-ATPase-containing metal efflux platforms mediating bacterial resistance to metal poisoning.

NMR-based structural glycomics for high-throughput screening of carbohydrate-active enzyme specificity.

We report here the development of a straightforward, sensitive, and quantitative NMR-based method for high-throughput characterization of carbohydrate structure and screening of carbohydrate active enzyme (CAZyme) specificity. Automated assays starting from gene library expression to carbohydrate structure determination directly from crude reaction media have been established and successfully used to screen a library of 4032 CAZymes obtained by combinatorial engineering, at a rate of 480 enzyme variants per day. This allowed one to accurately discriminate 303 enzyme variants with altered specificity. The results demonstrate the potential of high-throughput NMR technology in glycomics, to mine artificial and natural enzyme diversity for novel biocatalysts.

Order parameters of a transmembrane helix in a fluid bilayer: case study of a WALP peptide.

A new solid-state NMR-based strategy is established for the precise and efficient analysis of orientation and dynamics of transmembrane peptides in fluid bilayers. For this purpose, several dynamically averaged anisotropic constraints, including (13)C and (15)N chemical shift anisotropies and (13)C-(15)N dipolar couplings, were determined from two different triple-isotope-labeled WALP23 peptides ((2)H, (13)C, and (15)N) and combined with previously published quadrupolar splittings of the same peptide. Chemical shift anisotropy tensor orientations were determined with quantum chemistry. The complete set of experimental constraints was analyzed using a generalized, four-parameter dynamic model of the peptide motion, including tilt and rotation angle and two associated order parameters. A tilt angle of 21 degrees was determined for WALP23 in dimyristoylphosphatidylcholine, which is much larger than the tilt angle of 5.5 degrees previously determined from (2)H NMR experiments. This approach provided a realistic value for the tilt angle of WALP23 peptide in the presence of hydrophobic mismatch, and can be applied to any transmembrane helical peptide. The influence of the experimental data set on the solution space is discussed, as are potential sources of error.

GATEWAY technology and E. coli recombinant system produce a properly folded and functional recombinant allergen of the lipid transfer protein of apple (Mal d 3).

The lipid transfer protein of apple fruit, Mal d 3, has been produced as a soluble recombinant protein in transformed Escherichia coli cells using the GATEWAY technology. Circular dichroism spectra showing the protein essentially consists of alpha-helices indicate that the rMal d 3 is properly folded. The (1)H NMR spectra also indicates a correct fold for the recombinant allergen. The reactivity of rMal d 3 towards IgE from apple allergic patients and in vitro degranulation activity measured on transformed rat basophil leukemia cells expressing the human Fc epsilon RI alpha-subunit of rMal d 3 is similar to those of the native allergen purified from apple fruits. The expression of active rMal d 3 in E. coli is readily feasible and offers an interesting alternative to the production of recombinant allergen in the yeast Pichia pastoris. This expression in E. coli open the way to the modification of Mal d 3 by site-directed mutagenesis for immunotherapy purposes.

Conformational flexibility of coenzyme A and its impact on the post-translational modification of acyl carrier proteins by 4'-phosphopantetheinyl transferases.

One central question surrounding the biosynthesis of fatty acids and polyketide-derived natural products is how the 4'-phosphopantetheinyl transferase (PPTase) interrogates the essential acyl carrier protein (ACP) domain to fulfill the initial activation step. The triggering factor of this study was the lack of structural information on PPTases at physiological pH, which could bias our comprehension of the mechanism of action of these important enzymes. Structural and functional studies on the family II PPTase PptAb of Mycobacterium abscessus show that pH has a profound effect on the coordination of metal ions and on the conformation of endogenously bound coenzyme A (CoA). The observed conformational flexibility of CoA at physiological pH is accompanied by a disordered 4'-phosphopantetheine (Ppant) moiety. Finally, structural and dynamical information on an isolated mycobacterial ACP domain, in its apo form and in complex with the activator PptAb, suggests an alternate mechanism for the post-translational modification of modular megasynthases.

Structural insights on the pamoic acid and the 8 kDa domain of DNA polymerase beta complex: towards the design of higher-affinity inhibitors.

BACKGROUND: DNA polymerase beta (pol beta), the error-prone DNA polymerase of single-stranded DNA break repair as well as base excision repair pathways, is overexpressed in several tumors and takes part in chemotherapeutic agent resistance, like that of cisplatin, through translesion synthesis. For this reason pol beta has become a therapeutic target. Several inhibitors have been identified, but none of them presents a sufficient affinity and specificity to become a drug. The fragment-based inhibitor design allows an important improvement in affinity of small molecules. The initial and critical step for setting up the fragment-based strategy consists in the identification and structural characterization of the first fragment bound to the target. RESULTS: We have performed docking studies of pamoic acid, a 9 micromolar pol beta inhibitor, and found that it binds in a single pocket at the surface of the 8 kDa domain of pol beta. However, docking studies provided five possible conformations for pamoic acid in this site. NMR experiments were performed on the complex to select a single conformation among the five retained. Chemical Shift Mapping data confirmed pamoic acid binding site found by docking while NOESY and saturation transfer experiments provided distances between pairs of protons from the pamoic acid and those of the 8 kDa domain that allowed the identification of the correct conformation. CONCLUSION: Combining NMR experiments on the complex with docking results allowed us to build a three-dimensional structural model. This model serves as the starting point for further structural studies aimed at improving the affinity of pamoic acid for binding to DNA polymerase beta.