5-HT6 receptor antagonists. Design, synthesis, and structure-activity relationship of substituted 2-(1-methyl-4-piperazinyl)pyridines.

In this study, we present the discovery and pharmacological characterization of a new series of 6-piperazinyl-7-azaindoles. These compounds demonstrate potent antagonism and selectivity against the 5-HT6 receptor. Our research primarily focuses on optimizing the lead structure and investigating the structure-activity relationship (SAR) of these compounds. Our main objective is to improve their activity and selectivity against off-target receptors. Overall, our findings contribute to the advancement of novel compounds targeting the 5-HT6 receptor. Compound 29 exhibits significant promise in terms of pharmacological, physicochemical, and ADME (Absorption, Distribution, Metabolism, and Excretion) properties. Consequently, it merits thorough exploration as a potential drug candidate due to its favorable activity profile and successful outcomes in a range of in vivo experiments.

1-Sulfonyl-6-Piperazinyl-7-Azaindoles as potent and pseudo-selective 5-HT6 receptor antagonists.

A series of 1-Sulfonyl-6-Piperazinyl-7-Azaindoles, showing strong antagonistic activity to 5-HT6 receptor (5-HT6R) was synthesized and characterized. The series was optimized to reduce activity on D2 receptor. Based on the selectivity against this off-target and the analysis of the ADME-tox profile, compound 1c was selected for in vivo efficacy assessment, which demonstrated procognitive effects as shown in reversal of scopolamine induced amnesia in an elevated plus maze test in mice. Compound 3, the demethylated version of compound 1c, was profiled against a panel of 106 receptors, channels and transporters, indicating only D3 receptor as a major off-target. Compound 3 has been selected for this study over compound 1c because of the higher 5-HT6R/D2R binding ratio. These results have defined a new direction for the design of our pseudo-selective 5-HT6R antagonists.

Rational engineering of L-asparaginase reveals importance of dual activity for cancer cell toxicity.

Using proteins in a therapeutic context often requires engineering to modify functionality and enhance efficacy. We have previously reported that the therapeutic antileukemic protein macromolecule Escherichia coli L-asparaginase is degraded by leukemic lysosomal cysteine proteases. In the present study, we successfully engineered L-asparaginase to resist proteolytic cleavage and at the same time improve activity. We employed a novel combination of mutant sampling using a genetic algorithm in tandem with flexibility studies using molecular dynamics to investigate the impact of lid-loop and mutations on drug activity. Applying these methods, we successfully predicted the more active L-asparaginase mutants N24T and N24A. For the latter, a unique hydrogen bond network contributes to higher activity. Furthermore, interface mutations controlling secondary glutaminase activity demonstrated the importance of this enzymatic activity for drug cytotoxicity. All selected mutants were expressed, purified, and tested for activity and for their ability to form the active tetrameric form. By introducing the N24A and N24A R195S mutations to the drug L-asparaginase, we are a step closer to individualized drug design.

Molecular basis of reduced glucosylceramidase activity in the most common Gaucher disease mutant, N370S.

Gaucher disease is caused by the defective activity of the lysosomal hydrolase, glucosylceramidase. Although the x-ray structure of wild type glucosylceramidase has been resolved, little is known about the structural features of any of the >200 mutations. Various treatments for Gaucher disease are available, including enzyme replacement and chaperone therapies. The latter involves binding of competitive inhibitors at the active site to enable correct folding and transport of the mutant enzyme to the lysosome. We now use molecular dynamics, a set of structural analysis tools, and several statistical methods to determine the flexible behavior of the N370S Gaucher mutant at various pH values, with and without binding the chaperone, N-butyl-deoxynojirimycin. We focus on the effect of the chaperone on the whole protein, on the active site, and on three important structural loops, and we demonstrate how the chaperone modifies the behavior of N370S in such a way that it becomes more active at lysosomal pH. Our results suggest a mechanism whereby the binding of N-butyl-deoxynojirimycin helps target correctly folded glucosylceramidase to the lysosome, contributes to binding with saposin C, and explains the initiation of the substrate-enzyme complex. Such analysis provides a new framework for determination of the structure of other Gaucher disease mutants and suggests new approaches for rational drug design.

Influence of the electric field on supramolecular structure and properties of amyloid-specific reagent Congo red.

Among specific amyloid ligands, Congo red and its analogues are often considered potential therapeutic compounds. However, the results of the studies so far have not been univocal because the properties of this dye, derived mostly from its supramolecular nature, are still poorly understood. The supramolecular structure of Congo red, formed by π-π stacking of dye molecules, is susceptible to the influence of the electric field, which may significantly facilitate electron delocalization. Consequently, the electric field may generate altered physico-chemical properties of the dye. Enhanced electron delocalization, induced by the electric field, alters the total charge of Congo red, making the dye more acidic (negatively charged). This is a consequence of withdrawing electrons from polar substituents of aromatic rings-sulfonic and amino groups-thus increasing their tendency to dissociate protons. The electric field-induced charge alteration observed in electrophoresis depends on dye concentration. This concentration-dependent charge alteration effect disappears when the supramolecular structure disintegrates in DMSO. Dipoles formed from supramolecular fibrillar species in the electric field become ordered in the solution, introducing the modified arrangement to liquid crystalline phase. Experimental results and theoretical studies provide evidence confirming predictions that the supramolecular character of Congo red is the main reason for its specific properties and reactivity.

Discovery of 5-{2-[5-Chloro-2-(5-ethoxyquinoline-8-sulfonamido)phenyl]ethynyl}-4-methoxypyridine-2-carboxylic Acid, a Highly Selective in Vivo Useable Chemical Probe to Dissect MCT4 Biology.

Due to increased lactate production during glucose metabolism, tumor cells heavily rely on efficient lactate transport to avoid intracellular lactate accumulation and acidification. Monocarboxylate transporter 4 (MCT4/SLC16A3) is a lactate transporter that plays a central role in tumor pH modulation. The discovery and optimization of a novel class of MCT4 inhibitors (hit 9a), identified by a cellular screening in MDA-MB-231, is described. Direct target interaction of the optimized compound 18n with the cytosolic domain of MCT4 was shown after solubilization of the GFP-tagged transporter by fluorescence cross-correlation spectroscopy and microscopic studies. In vitro treatment with 18n resulted in lactate efflux inhibition and reduction of cellular viability in MCT4 high expressing cells. Moreover, pharmacokinetic properties of 18n allowed assessment of lactate modulation and antitumor activity in a mouse tumor model. Thus, 18n represents a valuable tool for investigating selective MCT4 inhibition and its effect on tumor biology.

How to design potent and selective DYRK1B inhibitors? Molecular modeling study.

DYRK1B protein kinase is an emerging anticancer target due to its overexpression in a variety of cancers and its role in cancer chemoresistance through maintaining cancer cells in the G0 (quiescent) state. Consequently, there is a growing interest in the development of potent and selective DYRK1B inhibitors for anticancer therapy. One of the major off-targets is another protein kinase, GSK3ß, which phosphorylates an important regulator of cell cycle progression on the same residue as DYRK1B and is involved in multiple signaling pathways. In the current work, we performed a detailed comparative structural analysis of DYRK1B and GSK3ß ATP-binding sites and identified key regions responsible for selectivity. As the crystal structure of DYRK1B has never been reported, we built and optimized a homology model by comparative modeling and metadynamics simulations. Calculation of interaction energies between docked ligands in the ATP-binding sites of both kinases allowed us to pinpoint key residues responsible for potency and selectivity. Specifically, the role of the gatekeeper residues in DYRK1B and GSK3ß is discussed in detail, and two other residues are identified as key to selectivity of DYRK1B inhibition versus GSK3ß. The analysis presented in this work was used to support the design of potent and selective azaindole-quinoline-based DYRK1B inhibitors and can facilitate development of more selective inhibitors for DYRK kinases.

Flexible relaxation of rigid-body docking solutions.

Molecular Dynamics (MD) simulations have been performed on a set of rigid-body docking poses, carried out over 25 protein-protein complexes. The results show that fully flexible relaxation increases the fraction of native contacts (NC) by up to 70% for certain docking poses. The largest increase in the fraction of NC is observed for docking poses where anchor residues are able to sample their bound conformation. For each MD simulation, structural snap-shots were clustered and the centre of each cluster used as the MD-relaxed docking pose. A comparison between two energy-based scoring schemes, the first calculated for the MD-relaxed poses, the second for energy minimized poses, shows that the former are better in ranking complexes with large hydrophobic interfaces. Furthermore, complexes with large interfaces are generally ranked well, regardless of the type of relaxation method chosen, whereas complexes with small hydrophobic interfaces remain difficult to rank. In general, the results indicate that current force-fields are able to correctly describe direct intermolecular interactions between receptor and ligand molecules. However, these force-fields still fail in cases where protein-protein complexes are stabilized by subtle energy contributions.

Implicit flexibility in protein docking: cross-docking and local refinement.

In previous CAPRI rounds (3-5) we showed that using MD-generated ensembles, as inputs for a rigid-body docking algorithm, increased our success rate, especially for targets exhibiting substantial amounts of induced fit. In recent rounds (6-11), our cross-docking was followed by a short MD-based local refinement for the subset of solutions with the lowest interaction energies after minimization. The above approach showed promising results for target 20, where we were able to recover 30% of native contacts for one of our submitted models. Further tests, performed a posteriori, revealed that cross-docking approach produces more near-native (NN) solutions but only for targets with large conformational changes upon binding. However, at the time of the blind docking experiment, these improved solutions were not chosen for the subsequent refinement, as their interaction energies after minimization ranked poorly compared with other solutions. This indicates deficiencies in the present scoring schemes that are based on interaction energies of minimized structures. Refinement MD simulations substantially increase the fraction of native contacts for NN docked solutions, but generally worsen interface and ligand RMSD. Further analysis shows that although MD simulations are able to improve sidechain packing across the interface, which results in an increased fraction of native contacts, they are not capable of improving interface and ligand backbone RMSD for NN structures beyond 1.5 and 3.5 A, respectively, even if explicit solvent is used.

The increased flexibility of CDR loops generated in antibodies by Congo red complexation favors antigen binding.

The dye Congo red and related self-assembling compounds were found to stabilize immune complexes by binding to antibodies currently engaged in complexation to antigen. In our simulations, it was shown that the site that becomes accessible for binding the supramolecular dye ligand is located in the V domain, and is normally occupied by the N-terminal polypeptide chain fragment. The binding of the ligand disrupts the beta-structure in the domain, increasing the plasticity of the antigen-binding site. The higher fluctuation of CDR-bearing loops enhances antigen binding, and allows even low-affinity antibodies to be engaged in immune complexes. Experimental observations of the enhancement effect were supported by theoretical studies using L lambda chain (4BJL-PDB identification) and the L chain from the complex of IgM-rheumatoid factor bound to the CH3 domain of the Fc fragment (1ADQ-PDB identification) as the initial structures for theoretical studies of dye-induced changes. Commercial IgM-type rheumatoid factor (human) and sheep red blood cells with coupled IgG (human) were used for experimental tests aimed to reveal the dye-enhancement effect in this system. The specificity of antigen-antibody interaction enhanced by dye binding was studied using rabbit anti-sheep red cell antibodies to agglutinate red cells of different species. Red blood cells of hoofed mammals (horse, goat) showed weak enhancement of agglutination in the presence of Congo red. Neither agglutination nor enhancement were observed in the case of human red cells. The dye-enhancement capability in the SRBC-antiSRBC system was lost after pepsin-digestion of antibodies producing (Fab)2 fragments still agglutinating red cells. Monoclonal (myeloma) IgG, L lambda chain and ovoalbumin failed to agglutinate red cells, as expected, and showed no enhancement effect. This indicates that the enhancement effect is specific.

The use of the Congo red-related dye DBACR to recognize the heavy chain-derived abnormality of myeloma immunoglobulins.

INTRODUCTION: The aim of this study was to differentiate heavy and light chain-derived instability of monoclonal myeloma immunoglobulins by complexation of matched supramolecular dyes. These are composed of several micellar pieces of self-assembled dye molecules which may penetrate the protein interior of the binding locus with polypeptide chains. These dyes were used to elicit, by precipitation, the postulated higher aggregation tendency of the heavy chain derived from its higher hydrophobicity. MATERIALS AND METHODS: Agarose gel electrophoresis was used to create conditions for dye complexation and to reveal the precipitation. RESULTS: Congo red derivatives with aromatic ring substitutes, BACR and DBACR, of increased penetrating capability were chosen to provoke the precipitation of abnormal immunoglobulins by displacing association-prone polypeptide chains from the protein interior. CONCLUSIONS: The results of this study confirm the heavy chain-related propensity of some monoclonal immunoglobulins to aggregate and precipitate. The simplicity of the technique may improve clinical diagnosis and facilitate predictions of disease complications.

Albumin binds self-assembling dyes as specific polymolecular ligands.

Self-assembling dyes with a structure related to Congo red (e.g. Evans blue) form polymolecular complexes with albumin. The dyes, which are lacking a self-assembling property (Trypan blue, ANS) bind as single molecules. The supramolecular character of dye ligands bound to albumin was demonstrated by indicating the complexation of dye molecules outnumbering the binding sites in albumin and by measuring the hydrodynamic radius of albumin which is growing upon complexation of self-assembling dye in contrast to dyes lacking this property. The self-assembled character of Congo red was also proved using it as a carrier introducing to albumin the intercalated nonbonding foreign compounds. Supramolecular, ordered character of the dye in the complex with albumin was also revealed by finding that self-assembling dyes become chiral upon complexation. Congo red complexation makes albumin less resistant to low pH as concluded from the facilitated N-F transition, observed in studies based on the measurement of hydrodynamic radius. This particular interference with protein stability and the specific changes in digestion resulted from binding of Congo red suggest that the self-assembled dye penetrates the central crevice of albumin.

Analysis of correlated domain motions in IgG light chain reveals possible mechanisms of immunological signal transduction.

It was shown experimentally that binding of a micelle composed of Congo red molecules to immunological complexes leads to the enhanced stability of the latter, and simultaneously prevents binding of a complement molecule (C1q). The dye binds in a cavity created by the removal of N-terminal polypeptide chain, as observed experimentally in a model system-immunoglobulin G (IgG) light chain dimer. Molecular Dynamics (MD) simulations of three forms of IgG light chain dimer, with and without the dye, were performed to investigate the role of N-terminal fragment and self-assembled ligand in coupling between V and C domains. Root-mean-square distance (RMSD) time profiles show that removal of N-terminal fragment leads to destabilization of V domain. A micelle composed of four self-assembled dye molecules stabilizes and fixes the domain. Analysis of root-mean-square fluctuation (RMSF) values and dynamic cross-correlation matrices (DCCM) reveals that removal of N-terminal fragment results in complete decoupling between V and C domains. Binding of self-assembled Congo red molecules improves the coupling, albeit slightly. The disruption of a small beta-sheet composed of N- and C-terminal fragments of the domain (NC sheet) is the most likely reason for the decoupling. Self-assembled ligand, bound in the place originally occupied by N-terminal fragment, is not able to take over the function of the beta-sheet. Lack of correlation of motions between residues in V and C domains denotes that light chain-Congo red complexes have hampered ability to transmit conformational changes between domains. This is a likely explanation of the lack of complement binding by immunological complexes, which bind Congo red, and supports the idea that the NC sheet is the key structural fragment taking part in immunological signal transduction.

Instability of monoclonal myeloma protein may be identified as susceptibility to penetration and binding by newly synthesized Congo red derivatives.

Monoclonal myeloma proteins often have an abnormal, unstable structure, and tend to aggregate with fatal clinical consequences. A method for early clinical identification of this aggregation tendency is impatiently awaited. This work proposes the use of supramolecular dyes as specific ligands to reveal protein instability. Disclosure of excessive polypeptide chain flexibility in unstable monoclonal proteins, leading to increased susceptibility to penetration by foreign compounds, appeared possible when new supramolecular Congo red-derived dyes with different protein-binding capabilities were used for complexation. Two basic protein instability levels, local and global, were differentiated by comparing the extent of protein loading with dye and the subsequent electrophoretic migration rate of the complexes. A simple electrophoretic test is proposed for assessment of the instability of monoclonal proteins in clinical conditions.

The structure and protein binding of amyloid-specific dye reagents.

The self-assembling tendency and protein complexation capability of dyes related to Congo red and also some dyes of different structure were compared to explain the mechanism of Congo red binding and the reason for its specific affinity for beta-structure. Complexation with proteins was measured directly and expressed as the number of dye molecules bound to heat-aggregated IgG and to two light chains with different structural stability. Binding of dyes to rabbit antibodies was measured indirectly as the enhancement effect of the dye on immune complex formation. Self-assembling was tested using dynamic light scattering to measure the size of the supramolecular assemblies. In general the results show that the supramolecular form of a dye is the main factor determining its complexation capability. Dyes that in their compact supramolecular organization are ribbon-shaped may adhere to polypeptides of beta-conformation due to the architectural compatibility in this unique structural form. The optimal fit in complexation seems to depend on two contradictory factors involving, on the one hand, the compactness of the non-covalently stabilized supramolecular ligand, and the dynamic character producing its plasticity on the other. As a result, the highest protein binding capability is shown by dyes with a moderate self-assembling tendency, while those arranging into either very rigid or very unstable supramolecular entities are less able to bind.

The use of supramolecular structures as protein ligands.

Congo red dye as well as other eagerly self-assembling organic molecules which form rod-like or ribbon-like supramolecular structures in water solutions, appears to represent a new class of protein ligands with possible wide-ranging medical applications. Such molecules associate with proteins as integral clusters and preferentially penetrate into areas of low molecular stability. Abnormal, partly unfolded proteins are the main binding target for such ligands, while well packed molecules are generally inaccessible. Of particular interest is the observation that local susceptibility for binding supramolecular ligands may be promoted in some proteins as a consequence of function-derived structural changes, and that such complexation may alter the activity profile of target proteins. Examples are presented in this paper.

Comparison of a molecular dynamics model with the X-ray structure of the N370S acid-beta-glucosidase mutant that causes Gaucher disease.

Recently, two studies were published that examined the structure of the acid-ß-glucosidase N370S mutant, the most common mutant that causes Gaucher disease. One study used the experimental tool of X-ray crystallography, and the other utilized molecular dynamics (MD). The two studies reinforced each other through the similarities in their findings, but each approach also added some unique information. Both studies report that the conformation of active site loop 3 changes, due to an altered hydrogen bonding network; however, the MD study produced additional data concerning the flexibility of loop 1 and the catalytic residues that are not observed in the other study.

A dyad of lymphoblastic lysosomal cysteine proteases degrades the antileukemic drug L-asparaginase.

l-Asparaginase is a key therapeutic agent for treatment of childhood acute lymphoblastic leukemia (ALL). There is wide individual variation in pharmacokinetics, and little is known about its metabolism. The mechanisms of therapeutic failure with l-asparaginase remain speculative. Here, we now report that 2 lysosomal cysteine proteases present in lymphoblasts are able to degrade l-asparaginase. Cathepsin B (CTSB), which is produced constitutively by normal and leukemic cells, degraded asparaginase produced by Escherichia coli (ASNase) and Erwinia chrysanthemi. Asparaginyl endopeptidase (AEP), which is overexpressed predominantly in high-risk subsets of ALL, specifically degraded ASNase. AEP thereby destroys ASNase activity and may also potentiate antigen processing, leading to allergic reactions. Using AEP-mediated cleavage sequences, we modeled the effects of the protease on ASNase and created a number of recombinant ASNase products. The N24 residue on the flexible active loop was identified as the primary AEP cleavage site. Sole modification at this site rendered ASNase resistant to AEP cleavage and suggested a key role for the flexible active loop in determining ASNase activity. We therefore propose what we believe to be a novel mechanism of drug resistance to ASNase. Our results may help to identify alternative therapeutic strategies with the potential of further improving outcome in childhood ALL.

The use of rigid, fibrillar Congo red nanostructures for scaffolding protein assemblies and inducing the formation of amyloid-like arrangement of molecules.

The ordered amyloid-like organization of protein aggregates was obtained using for their formation the rigid fibrillar nanostructures of Congo red as the scaffolding. The higher rigidity of used dye nanoparticles resulted from the stronger stacking of molecules at low pH (near the pK of the dye amino group) because of the decreased charge repulsion. The polylysine, human globin, and immunoglobulin L chain were arranged in this way to form deposits of amyloid properties. The scaffolding was introduced simply by mixing the dye and proteins at a low pH or the dye was used in the preorganized form by maintaining it in the electric field before and during protein addition. The polarization and electron microscopy studies confirmed the unidirectional organization of the complex. The precipitate of the complex was used for studies directly or after the partial or complete removal of the dye. The results suggest that the process of formation of amyloid-like deposits may bypass the nucleation step. It is possible if the protein aggregation occurs in unidirectionally organized (because of scaffolding) assembly of molecules, arranged prior to self-association. The recognition of the structure of amphoteric Congo red nanoparticles used for the scaffolding was based on the molecular dynamics simulation.

Macroscopic pKa Calculations for Fluorescein and Its Derivatives.

This study describes the calculation of the microscopic dissociation and tautomerization constants of fluorescein and its derivatives, 2',7'-dichlorofluorescein (DCF) and 2',7'-difluorofluorescein (DFF), in an aqueous environment. In vacuo free energies were obtained using complete basis set (CBS) and DFT-based methods, while free energies of solvation were calculated with the CPCM implicit solvation protocol using the UAHF, UAKS, and Pauling radii sets. Our results indicate that the different vacuum protocols give free energy changes upon dissociation within 1 kcal/mol of each other for a given molecule. Therefore, we suggest that the computationally less intensive PBE1PBE/6-311+G(2d,2p)//PBE1PBE/6-31+G(d) model chemistry may reasonably be used in pKa calculations of larger molecules. The calculations also provided a rigorous test of the implicit solvation models. Relative calculations of dissociation constants gave results in good agreement with experiment; absolute values deviated from experimental data by 1-3 pKa units. Consistently better results were obtained with the Pauling radii set. The influence of geometry relaxation on going from vacuum to solvent is negligible for pKa2 and larger for pKa1 but still smaller than the variation due to the radii set. Calculation of tautomerization constants gave more variable results, with none of the solvation methods able to reproduce experimental values consistently, although certain individual constants were correctly calculated.