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).

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.

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.

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.

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.

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.

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.

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.

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.

NMR localization of the O-mycoloylation on PorH, a channel forming peptide from Corynebacterium glutamicum.

PorH and PorA are two small peptides that, in complex, form a voltage-dependent ion channel in the outer membrane of Corynebacterium glutamicum. Specific post-translational modifications on PorA and PorH are required for the formation of a functional ion channel. The assignment of PorH proton NMR chemical shifts in DMSO, allowed identifying unambiguously the exact position of the PorH O-mycoloylation on Ser 56 side chain. This was further confirmed by site directed mutagenesis and mass spectrometry. Together with the previously published localization of PorA mycoloylation, this provides the complete primary structure characterization of this outer membrane porin.

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.

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.

Engineering transglycosidase activity into a GH51 α-l-arabinofuranosidase.

Directed evolution was applied to the α-l-arabinofuranosidase from Thermobacillus xylanilyticus to confer better transglycosylation ability, particularly for the synthesis of benzyl α-l-arabinofuranosyl-(1,2)-α-d-xylopyranoside, starting from p-nitrophenyl α-l-arabinofuranoside (donor) and benzyl α-d-xylopyranoside (acceptor). The aim was to obtain mutants displaying both lower hydrolytic and greater transglycosylation activities to favour the stable production of the target disaccharide. The implementation of a simple chromogenic screen ultimately provided three mutant enzymes whose properties correspond to those sought after. These all displayed lowered hydrolytic activity and conserved or slightly improved transfer activity, while one of them also displayed lowered secondary hydrolysis of the transglycosylation product. DNA sequence analysis of the mutants revealed between three and seven point mutations and biochemical analysis combined with STD-NMR experiments indicated that distinct molecular mechanisms were active among the three mutants.

Functional roles of H98 and W99 and ß2α2 loop dynamics in the α-l-arabinofuranosidase from Thermobacillus xylanilyticus.

This study is focused on the elucidation of the functional role of the mobile ß2α2 loop in the α-L-arabinofuranosidase from Thermobacillus xylanilyticus, and particularly on the roles of loop residues H98 and W99. Using site-directed mutagenesis, coupled to characterization methods including isothermal titration calorimetry (ITC) and saturation transfer difference nuclear magnetic resonance (STD-NMR) spectroscopy, and molecular dynamics simulations, it has been possible to provide a molecular level view of interactions and the consequences of mutations. Binding of para-nitrophenyl α-L-arabinofuranoside (pNP-α-l-Araf) to the wild-type arabinofuranosidase was characterized by K(d) values (0.32 and 0.16 mm, from ITC and STD-NMR respectively) that highly resembled that of the arabinoxylo-oligosaccharide XA(3)XX (0.21 mm), and determination of the thermodynamic parameters of enzyme : pNP-α-L-Araf binding revealed that this process is driven by favourable entropy, which is linked to the movement of the ß2α2 loop. Loop closure relocates the solvent-exposed W99 into a buried location, allowing its involvement in substrate binding and in the formation of a functional active site. Similarly, the data underline the role of H98 in the 'dynamic' formation and definition of a catalytically operational active site, which may be a specific feature of a subset of GH51 arabinofuranosidases. Substitution of H98 and W99 by alanine or phenylalanine revealed that mutations affected K(M) and/or k(cat). Molecular dynamics performed on W99A implied that this mutation causes the loss of a hydrogen bond and leads to an alternative binding mode that is detrimental for catalysis. STD-NMR experiments revealed altered binding of the aglycon motif in the active site, combined with reduced STD intensities of the α-L-arabinofuranosyl moiety for W99 substitutions.

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.

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.

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.

Solution-state NMR spectroscopy of membrane proteins in detergent micelles: structure of the Klebsiella pneumoniae outer membrane protein A, KpOmpA.

Structure determination of integral membrane proteins is one of the most important challenges of structural biology. Over the last 7 years, solution-state NMR spectroscopy has become an increasingly useful approach for 3D structure determination and dynamical analysis of membrane proteins solubilized in detergent micelles. We describe herein an ensemble of methods applied in this context, including isotopic labelling, in vitro refolding procedure, and state-of-the-art NMR experiments required for the structure determination of high molecular weight molecular complexes. Furthermore, the basic principles of spectrum interpretation and 3D structure calculation are reported. This approach is illustrated with the structural study of the transmembrane domain of the outer membrane protein A from Klebsiella pneumoniae (KpOmpA).