Mapping the human brain proteome: opportunities, challenges, and clinical potential.

INTRODUCTION: Due to the segmented functions and complexity of the human brain, the characterization of molecular profiles within specific areas such as brain structures and biofluids is essential to unveil the molecular basis for structure specialization as well as the molecular imbalance associated with neurodegenerative and psychiatric diseases. AREAS COVERED: Much of our knowledge about brain functionality derives from neurophysiological, anatomical, and transcriptomic approaches. More recently, laser capture and imaging proteomics, technological and computational developments in LC-MS/MS, as well as antibody/aptamer-based platforms have allowed the generation of novel cellular, spatial, and posttranslational dimensions as well as innovative facets in biomarker validation and druggable target identification. EXPERT OPINION: Proteomics is a powerful toolbox to functionally characterize, quantify, and localize the extensive protein catalog of the human brain across physiological and pathological states. Brain function depends on multi-dimensional protein homeostasis, and its elucidation will help us to characterize biological pathways that are essential to properly maintain cognitive functions. In addition, comprehensive human brain pathological proteomes may be the basis in computational drug-repositioning methods as a strategy for unveiling potential new therapies in neurodegenerative and psychiatric disorders.

Use of the Novel Site-Directed Enzyme Enhancement Therapy (SEE-Tx) Drug Discovery Platform to Identify Pharmacological Chaperones for Glutaric Acidemia Type 1.

Allosteric regulators acting as pharmacological chaperones hold promise for innovative therapeutics since they target noncatalytic sites and stabilize the folded protein without competing with the natural substrate, resulting in a net gain of function. Exogenous allosteric regulators are typically more selective than active site inhibitors and can be more potent than competitive inhibitors when the natural substrate levels are high. To identify novel structure-targeted allosteric regulators (STARs) that bind to and stabilize the mitochondrial enzyme glutaryl-CoA dehydrogenase (GCDH), the computational site-directed enzyme enhancement therapy (SEE-Tx) technology was applied. SEE-Tx is an innovative drug discovery platform with the potential to identify drugs for treating protein misfolding disorders, such as glutaric acidemia type 1 (GA1) disease. Putative allosteric regulators were discovered using structure- and ligand-based virtual screening methods and validated using orthogonal biophysical and biochemical assays. The computational approach presented here could be used to discover allosteric regulators of other protein misfolding disorders.

Discovery of allosteric regulators with clinical potential to stabilize alpha-L-iduronidase in mucopolysaccharidosis type I.

Mucopolysaccharidosis type I (MPS I) is an inherited lysosomal disease caused by lowered activity of the enzyme alpha-L-iduronidase (IDUA). Current therapeutic options show limited efficacy and do not treat some important aspects of the disease. Therefore, it may be advantageous to identify strategies that could improve the efficacy of existing treatments. Pharmacological chaperones are small molecules that protect proteins from degradation, and their use in combination with enzyme replacement therapy (ERT) has been proposed as an alternative therapeutic strategy. Using the SEE-Tx® proprietary computational drug discovery platform, a new allosteric ligand binding cavity in IDUA was identified distal from the active site. Virtual high-throughput screening of approximately 5 million compounds using the SEE-Tx® docking platform identified a subset of small molecules that bound to the druggable cavity and functioned as novel allosteric chaperones of IDUA. Experimental validation by differential scanning fluorimetry showed an overall hit rate of 11.4%. Biophysical studies showed that one exemplary hit molecule GT-01803 bound to (Kd = 22 µM) and stabilized recombinant human IDUA (rhIDUA) in a dose-dependent manner. Co-administration of rhIDUA and GT-01803 increased IDUA activity in patient-derived fibroblasts. Preliminary in vivo studies have shown that GT-01803 improved the pharmacokinetic (PK) profile of rhIDUA, increasing plasma levels in a dose-dependent manner. Furthermore, GT-01803 also increased IDUA enzymatic activity in bone marrow tissue, which benefits least from standard ERT. Oral bioavailability of GT-01803 was found to be good (50%). Overall, the discovery and validation of a novel allosteric chaperone for rhIDUA presents a promising strategy to enhance the efficacy of existing treatments for MPS I. The compound's ability to increase rhIDUA activity in patient-derived fibroblasts and its good oral bioavailability underscore its potential as a potent adjunct to ERT, particularly for addressing aspects of the disease less responsive to standard treatment.

Neuropathological stage-dependent proteome mapping of the olfactory tract in Alzheimer's disease: From early olfactory-related omics signatures to computational repurposing of drug candidates.

Alzheimer's disease (AD) is the most common form of dementia, characterized by an early olfactory dysfunction, progressive memory loss, and behavioral deterioration. Albeit substantial progress has been made in characterizing AD-associated molecular and cellular events, there is an unmet clinical need for new therapies. In this study, olfactory tract proteotyping performed in controls and AD subjects (n = 17/group) showed a Braak stage-dependent proteostatic impairment accompanied by the progressive modulation of amyloid precursor protein and tau functional interactomes. To implement a computational repurposing of drug candidates with the capacity to reverse early AD-related olfactory omics signatures (OMSs), we generated a consensual OMSs database compiling differential omics datasets obtained by mass-spectrometry or RNA-sequencing derived from initial AD across the olfactory axis. Using the Connectivity Map-based drug repurposing approach, PKC, EGFR, Aurora kinase, Glycogen synthase kinase, and CDK inhibitors were the top pharmacologic classes capable to restore multiple OMSs, whereas compounds with targeted activity to inhibit PI3K, Insulin-like growth factor 1 (IGF-1), microtubules, and Polo-like kinase (PLK) represented a family of drugs with detrimental potential to induce olfactory AD-associated gene expression changes. To validate the potential therapeutic effects of the proposed drugs, in vitro assays were performed. These validation experiments revealed that pretreatment of human neuron-like SH-SY5Y cells with the EGFR inhibitor AG-1478 showed a neuroprotective effect against hydrogen peroxide-induced damage while the pretreatment with the Aurora kinase inhibitor Reversine reduced amyloid-beta (Aß)-induced neurotoxicity. Taken together, our data pointed out that OMSs may be useful as substrates for drug repurposing to propose novel neuroprotective treatments against AD.

Involvement of Glucosamine 6 Phosphate Isomerase 2 (GNPDA2) Overproduction in ß-Amyloid- and Tau P301L-Driven Pathomechanisms.

Alzheimer's disease (AD) is a neurodegenerative olfactory disorder affecting millions of people worldwide. Alterations in the hexosamine- or glucose-related pathways have been described through AD progression. Specifically, an alteration in glucosamine 6 phosphate isomerase 2 (GNPDA2) protein levels has been observed in olfactory areas of AD subjects. However, the biological role of GNPDA2 in neurodegeneration remains unknown. Using mass spectrometry, multiple GNPDA2 interactors were identified in human nasal epithelial cells (NECs) mainly involved in intraciliary transport. Moreover, GNPDA2 overexpression induced an increment in NEC proliferation rates, accompanied by transcriptomic alterations in Type II interferon signaling or cellular stress responses. In contrast, the presence of beta-amyloid or mutated Tau-P301L in GNPDA2-overexpressing NECs induced a slowdown in the proliferative capacity in parallel with a disruption in protein processing. The proteomic characterization of Tau-P301L transgenic zebrafish embryos demonstrated that GNPDA2 overexpression interfered with collagen biosynthesis and RNA/protein processing, without inducing additional changes in axonal outgrowth defects or neuronal cell death. In humans, a significant increase in serum GNPDA2 levels was observed across multiple neurological proteinopathies (AD, Lewy body dementia, progressive supranuclear palsy, mixed dementia and amyotrophic lateral sclerosis) (n = 215). These data shed new light on GNPDA2-dependent mechanisms associated with the neurodegenerative process beyond the hexosamine route.

Validation of a highly sensitive HaloTag-based assay to evaluate the potency of a novel class of allosteric ß-Galactosidase correctors.

Site-directed Enzyme Enhancement Therapy (SEE-Tx®) technology is a disease-agnostic drug discovery tool that can be applied to any protein target of interest with a known three-dimensional structure. We used this proprietary technology to identify and characterize the therapeutic potential of structurally targeted allosteric regulators (STARs) of the lysosomal hydrolase ß-galactosidase (ß-Gal), which is deficient due to gene mutations in galactosidase beta 1 (GLB1)-related lysosomal storage disorders (LSDs). The biochemical HaloTag cleavage assay was used to monitor the delivery of wildtype (WT) ß-Gal and four disease-related ß-Gal variants (p.Ile51Thr, p.Arg59His, p.Arg201Cys and p.Trp273Leu) in the presence and absence of two identified STAR compounds. In addition, the ability of STARs to reduce toxic substrate was assessed in a canine fibroblast cell model. In contrast to the competitive pharmacological chaperone N-nonyl-deoxygalactonojirimycin (NN-DGJ), the two identified STAR compounds stabilized and substantially enhanced the lysosomal transport of wildtype enzyme and disease-causing ß-Gal variants. In addition, the two STAR compounds reduced the intracellular accumulation of exogenous GM1 ganglioside, an effect not observed with the competitive chaperone NN-DGJ. This proof-of-concept study demonstrates that the SEE-Tx® platform is a rapid and cost-effective drug discovery tool for identifying STARs for the treatment of LSDs. In addition, the HaloTag assay developed in our lab has proved valuable in investigating the effect of STARs in promoting enzyme transport and lysosomal delivery. Automatization and upscaling of this assay would be beneficial for screening STARs as part of the drug discovery process.

Extracting Atomic Contributions to Binding Free Energy Using Molecular Dynamics Simulations with Mixed Solvents (MDmix).

BACKGROUND: Mixed solvents MD (MDmix) simulations have proved to be a useful and increasingly accepted technique with several applications in structure-based drug discovery. One of the assumptions behind the methodology is the transferability of free energy values from the simulated cosolvent molecules to larger drug-like molecules. However, the binding free energy maps (ΔGbind) calculated for the different moieties of the cosolvent molecules (e.g. a hydroxyl map for the ethanol) are largely influenced by the rest of the solvent molecule and do not reflect the intrinsic affinity of the moiety in question. As such, they are hardly transferable to different molecules. METHOD: To achieve transferable energies, we present here a method for decomposing the molecular binding free energy into accurate atomic contributions. RESULT: We demonstrate with two qualitative visual examples how the corrected energy maps better match known binding hotspots and how they can reveal hidden hotspots with actual drug design potential. CONCLUSION: Atomic decomposition of binding free energies derived from MDmix simulations provides transferable and quantitative binding free energy maps.

Structural properties of g,t-parallel duplexes.

The structure of G,T-parallel-stranded duplexes of DNA carrying similar amounts of adenine and guanine residues is studied by means of molecular dynamics (MD) simulations and UV- and CD spectroscopies. In addition the impact of the substitution of adenine by 8-aminoadenine and guanine by 8-aminoguanine is analyzed. The presence of 8-aminoadenine and 8-aminoguanine stabilizes the parallel duplex structure. Binding of these oligonucleotides to their target polypyrimidine sequences to form the corresponding G,T-parallel triplex was not observed. Instead, when unmodified parallel-stranded duplexes were mixed with their polypyrimidine target, an interstrand Watson-Crick duplex was formed. As predicted by theoretical calculations parallel-stranded duplexes carrying 8-aminopurines did not bind to their target. The preference for the parallel-duplex over the Watson-Crick antiparallel duplex is attributed to the strong stabilization of the parallel duplex produced by the 8-aminopurines. Theoretical studies show that the isomorphism of the triads is crucial for the stability of the parallel triplex.

Emerging opportunities and concerns for drug discovery at serotonin 5-HT2B receptors.

5-HT(2B) receptors have been reported to play an important role at cardiac, intestinal and central levels, as well as in bone marrow formation and growth. In the last decade, 5-HT(2B) receptors have also gained much attention as new targets in therapeutics, but also as off-targets because their activation along with the inhibition of serotonin transporters plays a significant role in the pathogenesis of 5-HT induced valvulopathy. Taking this into account, the present review focuses on the new therapeutic applications of 5-HT(2B) receptor ligands as well as on the potential concerns.

Discovery of 5-HT6 receptor ligands based on virtual HTS.

Based on a pharmacophore alignment on a 5-HT(6) ligand applying 4SCan technology, a new lead series was identified and further structurally investigated. K(i)s down to 8 nM were achieved.

Theoretical study of the Hoogsteen-Watson-Crick junctions in DNA.

A series of d (AT)(n) oligonucleotides containing mixtures of normal B-type Watson-Crick and antiparallel Hoogsteen helices have been studied using molecular dynamics simulation techniques to analyze the structural and thermodynamic impact of the junction between Watson-Crick and antiparallel Hoogsteen structures. Analysis of molecular dynamics simulations strongly suggests that for all oligonucleotides studied the antiparallel Hoogsteen appears as a reasonable conformation, only slightly less stable than the canonical B-type Watson-Crick one. The junctions between the Watson-Crick and Hoogsteen structures introduces a priori a sharp discontinuity in the helix, because the properties of each type of conformation are very well preserved in the corresponding fragments. However, and quite counterintuitively, junctions do not largely distort the duplex in structural, dynamics or energetic terms. Our results strongly support the possibility that small fragments of antiparallel Hoogsteen duplex might be embedded into large fragments of B-type Watson-Crick helices, making possible protein-DNA interactions that are specific of the antiparallel Hoogsteen conformation.

Theoretical study of the guanine --> 6-thioguanine substitution in duplexes, triplexes, and tetraplexes.

Molecular dynamics and thermodynamic integration calculations have been carried out on a set of G-rich single-strand, duplex, triplex, and quadruplex DNAs to study the structural and stability changes connected with the guanine --> 6-thioguanine (G --> S) mutation. The presence of 6-thioguanine leads to a shift of the geometry from the B/A intermediate to the pure B-form in duplex DNA. The G --> S mutation does not largely affect the structure of the antiparallel triplex when it is located at the reverse-Hoogsteen position, but leads to a non-negligible local distortion in the structure when it is located at the Watson-Crick position. The G --> S mutation leads to destabilization of all studied structures: the lowest effect has been observed for the G --> S mutation in the reverse-Hoogsteen strand of the triplex, a medium effect has been observed in the Watson-Crick strand of the triplex and duplex, and the highest influence of the G -->S mutation has been found for the quadruplex structures.

Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations?

The counterion distribution around a DNA dodecamer (5'-CGCGAATTCGCG-3') is analyzed using both standard and novel techniques based on state of the art molecular dynamics simulations. Specifically, we have explored the population of Na(+) in the minor groove of DNA duplex, and whether or not a string of Na(+) can replace the spine of hydration in the narrow AATT minor groove. The results suggest that the insertion of Na(+) in the minor groove is a very rare event, but that when once the ion finds specific sites deep inside the groove it can reside there for very long periods of time. According to our simulation the presence of Na(+) inside the groove does not have a dramatic influence in the structure or dynamics of the duplex DNA. The ability of current MD simulations to obtain equilibrated pictures of the counterion atmosphere around DNA is critically discussed.

Theoretical study of alkyl-pi and aryl-pi interactions. Reconciling theory and experiment.

Quantum mechanical and quantum mechanical/molecular mechanical calculations in conjunction with continuum solvation models have been used to analyze CH-pi interactions in model systems of aryl- and alkyl-aromatic interactions, as well as in a model folding system designed to study those interactions. High level calculations reproduced accurately the interaction of CH-pi interactions in both alkyl- and aryl-based model systems. Dispersion effects dominate the interaction, but the electrostatics term is also relevant for aryl CH-pi interactions. Theoretical calculations were also used to examine the influence of CH-pi interactions in determining the conformational flexibility of folding models. Finally, a critical comparison of the results obtained from high level calculations on model systems and the experimental data derived for folding models in apolar solvents was carried out, which allowed us to reconcile the apparent discrepancy found between both data.

Theoretical study of a new DNA structure: the antiparallel Hoogsteen duplex.

The structure of a new form of duplex DNA, the antiparallel Hoogsteen duplex, is studied in polyd(AT) sequences by means of state-of-the-art molecular dynamics simulations in aqueous solution. The structure, which was found to be stable in all of the simulations, has many similarities with the standard Watson-Crick duplex in terms of general structure, flexibility, and molecular recognition patterns. Accurate MM-PB/SA (and MM-GB/SA) analysis shows that the new structure has an effective energy similar to that of the B-type duplex, while it is slightly disfavored by intramolecular entropic considerations. Overall, MD simulations strongly suggest that the antiparallel Hoogsteen duplex is an accessible structure for a polyd(AT) sequence, which might compete under proper experimental conditions with normal B-DNA. MD simulations also suggest that chimeras containing Watson-Crick duplex and Hoogsteen antiparallel helices might coexist in a common structure, but with the differential characteristics of both type of structures preserved.

Antiparallel triple helices. Structural characteristics and stabilization by 8-amino derivatives.

The structural, dynamical, and recognition properties of antiparallel DNA triplexes formed by the antiparallel d(G#G.C), d(A#A.T), and d(T#A.T) motifs (the pound sign and dot mean reverse-Hoogsteen and Watson-Crick hydrogen bonds, respectively) are studied by means of "state of the art" molecular dynamics simulations. Once the characteristics of the helix are defined, molecular dynamics and thermodynamic integration calculations are used to determine the expected stabilization of the antiparallel triplex caused by the introduction of 8-aminopurines. Finally, oligonucleotides containing 8-aminopurine derivatives are synthesized and tested experimentally using several approaches in a variety of systems. A very large stabilization of the triplex is found experimentally, as predicted by simulations. These results open the possibility for the use of oligonucleotides carrying 8-aminopurines to bind single-stranded nucleic acids by formation of antiparallel triplexes.

Effect of bulky lesions on DNA: solution structure of a DNA duplex containing a cholesterol adduct.

The three-dimensional solution structure of two DNA decamers of sequence d(CCACXGGAAC)-(GTTCCGGTGG) with a modified nucleotide containing a cholesterol derivative (X) in its C1 '(chol)alpha or C1 '(chol)beta diastereoisomer form has been determined by using NMR and restrained molecular dynamics. This DNA derivative is recognized with high efficiency by the UvrB protein, which is part of the bacterial nucleotide excision repair, and the alpha anomer is repaired more efficiently than the beta one. The structures of the two decamers have been determined from accurate distance constraints obtained from a complete relaxation matrix analysis of the NOE intensities and torsion angle constraints derived from J-coupling constants. The structures have been refined with molecular dynamics methods, including explicit solvent and applying the particle mesh Ewald method to properly evaluate the long range electrostatic interactions. These calculations converge to well defined structures whose conformation is intermediate between the A- and B-DNA families as judged by the root mean square deviation but with sugar puckerings and groove shapes corresponding to a distorted B-conformation. Both duplex adducts exhibit intercalation of the cholesterol group from the major groove of the helix and displacement of the guanine base opposite the modified nucleotide. Based on these structures and molecular dynamics calculations, we propose a tentative model for the recognition of damaged DNA substrates by the UvrB protein.

Hoogsteen-based parallel-stranded duplexes of DNA. Effect of 8-amino-purine derivatives.

The structure of parallel-stranded duplexes of DNA-containing a mixture of guanines (G) and adenines (A) is studied by means of molecular dynamics (MD) simulation, as well as NMR and circular dichroism (CD) spectroscopy. Results demonstrate that the structure is based on the Hoogsteen motif rather than on the reverse Watson-Crick one. Molecular dynamics coupled to thermodynamic integration (MD/TI) calculations and melting experiments allowed us to determine the effect of 8-amino derivatives of A and G and of 8-amino-2'-deoxyinosine on the stability of parallel-stranded duplexes. The large stabilization of the parallel-stranded helix upon 8-amino substitution agrees with a Hoogsteen pairing, confirming MD, NMR, and CD data, and suggests new methods to obtain DNA triplexes for antigene and antisense purposes.