The dopamine transporter gene SLC6A3: multidisease risks

The human dopamine transporter gene SLC6A3 has been consistently implicated in several neuropsychiatric diseases but the disease mechanism remains elusive. In this risk synthesis, we have concluded that SLC6A3 represents an increasingly recognized risk with a growing number of familial mutants associated with neuropsychiatric and neurological disorders. At least five loci were related to common and severe diseases including alcohol use disorder (high activity variant), attention-deficit/hyperactivity disorder (low activity variant), autism (familial proteins with mutated networking) and movement disorders (both regulatory variants and familial mutations). Association signals depended on genetic markers used as well as ethnicity examined. Strong haplotype selection and gene-wide epistases support multimarker assessment of functional variations and phenotype associations. Inclusion of its promoter region’s functional markers such as DNPi (rs67175440) and 5’VNTR (rs70957367) may help delineate condensate-based risk action, testing a locus-pathway-phenotype hypothesis for one gene-multidisease etiology

Challenges Predicting Ligand-Receptor Interactions of Promiscuous Proteins: The Nuclear Receptor PXR

Transcriptional regulation of some genes involved in xenobiotic detoxification and apoptosis is performed via the human pregnane X receptor (PXR) which in turn is activated by structurally diverse agonists including steroid hormones. Activation of PXR has the potential to initiate adverse effects, altering drug pharmacokinetics or perturbing physiological processes. Reliable computational prediction of PXR agonists would be valuable for pharmaceutical and toxicological research. There has been limited success with structure-based modeling approaches to predict human PXR activators. Slightly better success has been achieved with ligand-based modeling methods including quantitative structure-activity relationship (QSAR) analysis, pharmacophore modeling and machine learning. In this study, we present a comprehensive analysis focused on prediction of 115 steroids for ligand binding activity towards human PXR. Six crystal structures were used as templates for docking and ligand-based modeling approaches (two-, three-, four- and five-dimensional analyses). The best success at external prediction was achieved with 5D-QSAR. Bayesian models with FCFP_6 descriptors were validated after leaving a large percentage of the dataset out and using an external test set. Docking of ligands to the PXR structure co-crystallized with hyperforin had the best statistics for this method. Sulfated steroids (which are activators) were consistently predicted as non-activators while, poorly predicted steroids were docked in a reverse mode compared to 5α-androstan-3β-ol. Modeling of human PXR represents a complex challenge by virtue of the large, flexible ligand-binding cavity. This study emphasizes this aspect, illustrating modest success using the largest quantitative data set to date and multiple modeling approaches

An NMR-based scoring function improves the accuracy of binding pose predictions by docking by two orders of magnitude

Low-affinity ligands can be efficiently optimized into high-affinity drug leads by structure based drug design when atomic-resolution structural information on the protein/ligand complexes is available. In this work we show that the use of a few, easily obtainable, experimental restraints improves the accuracy of the docking experiments by two orders of magnitude. The experimental data are measured in nuclear magnetic resonance spectra and consist of protein-mediated NOEs between two competitively binding ligands. The methodology can be widely applied as the data are readily obtained for low-affinity ligands in the presence of non-labelled receptor at low concentration. The experimental inter-ligand NOEs are efficiently used to filter and rank complex model structures that have been pre-selected by docking protocols. This approach dramatically reduces the degeneracy and inaccuracy of the chosen model in docking experiments, is robust with respect to inaccuracy of the structural model used to represent the free receptor and is suitable for high-throughput docking campaigns

Pyrazoleamide compounds are potent antimalarials that target Na+ homeostasis in intraerythrocytic Plasmodium falciparum

The quest for new antimalarial drugs, especially those with novel modes of action, is essential in the face of emerging drug-resistant parasites. Here we describe a new chemical class of molecules, pyrazoleamides, with potent activity against human malaria parasites and showing remarkably rapid parasite clearance in an in vivo model. Investigations involving pyrazoleamide-resistant parasites, whole-genome sequencing and gene transfers reveal that mutations in two proteins, a calcium-dependent protein kinase (PfCDPK5) and a P-type cation-ATPase (PfATP4), are necessary to impart full resistance to these compounds. A pyrazoleamide compound causes a rapid disruption of Na+ regulation in blood-stage Plasmodium falciparum parasites. Similar effect on Na+ homeostasis was recently reported for spiroindolones, which are antimalarials of a chemical class quite distinct from pyrazoleamides. Our results reveal that disruption of Na+ homeostasis in malaria parasites is a promising mode of antimalarial action mediated by at least two distinct chemical classes

Development of a Novel Virtual Screening Cascade Protocol to Identify Potential Trypanothione Reductase Inhibitors

The implementation of a novel sequential computational approach that can be used effectively for virtual screening and identification of prospective ligands that bind to trypanothione reductase (TryR) is reported. The multistep strategy combines a ligand-based virtual screening for building an enriched library of small molecules with a docking protocol (AutoDock, X-Score) for screening against the TryR target. Compounds were ranked by an exhaustive conformational consensus scoring approach that employs a rank-by-rank strategy by combining both scoring functions. Analysis of the predicted ligand-protein interactions highlights the role of bulky quaternary amine moieties for binding affinity. The scaffold hopping (SHOP) process derived from this computational approach allowed the identification of several chemotypes, not previously reported as antiprotozoal agents, which includes dibenzothiepine, dibenzooxathiepine, dibenzodithiepine, and polycyclic cationic structures like thiaazatetracyclo-nonadeca-hexaen-3-ium. Assays measuring the inhibiting effect of these compounds on T. cruzi and T. brucei TryR confirm their potential for further rational optimization

Inhibitors of Helicobacter pylori Protease HtrA Found by ‘Virtual Ligand’ Screening Combat Bacterial Invasion of Epithelia

Background: The human pathogen Helicobacter pylori (H. pylori) is a main cause for gastric inflammation and cancer. Increasing bacterial resistance against antibiotics demands for innovative strategies for therapeutic intervention. Methodology/Principal Findings: We present a method for structure-based virtual screening that is based on the comprehensive prediction of ligand binding sites on a protein model and automated construction of a ligand-receptor interaction map. Pharmacophoric features of the map are clustered and transformed in a correlation vector (‘virtual ligand’) for rapid virtual screening of compound databases. This computer-based technique was validated for 18 different targets of pharmaceutical interest in a retrospective screening experiment. Prospective screening for inhibitory agents was performed for the protease HtrA from the human pathogen H. pylori using a homology model of the target protein. Among 22 tested compounds six block E-cadherin cleavage by HtrA in vitro and result in reduced scattering and wound healing of gastric epithelial cells, thereby preventing bacterial infiltration of the epithelium. Conclusions/Significance: This study demonstrates that receptor-based virtual screening with a permissive (‘fuzzy’) pharmacophore model can help identify small bioactive agents for combating bacterial infection

Qualitative prediction of blood–brain barrier permeability on a large and refined dataset

The prediction of blood–brain barrier permeation is vitally important for the optimization of drugs targeting the central nervous system as well as for avoiding side effects of peripheral drugs. Following a previously proposed model on blood–brain barrier penetration, we calculated the cross-sectional area perpendicular to the amphiphilic axis. We obtained a high correlation between calculated and experimental cross-sectional area (r = 0.898, n = 32). Based on these results, we examined a correlation of the calculated cross-sectional area with blood–brain barrier penetration given by logBB values. We combined various literature data sets to form a large-scale logBB dataset with 362 experimental logBB values. Quantitative models were calculated using bootstrap validated multiple linear regression. Qualitative models were built by a bootstrapped random forest algorithm. Both methods found similar descriptors such as polar surface area, pKa, logP, charges and number of positive ionisable groups to be predictive for logBB. In contrast to our initial assumption, we were not able to obtain models with the cross-sectional area chosen as relevant parameter for both approaches. Comparing those two different techniques, qualitative random forest models are better suited for blood-brain barrier permeability prediction, especially when reducing the number of descriptors and using a large dataset. A random forest prediction system (ntrees = 5) based on only four descriptors yields a validated accuracy of 88%

Ligand-Dependent Conformations and Dynamics of the Serotonin 5-HT2A Receptor Determine Its Activation and Membrane-Driven Oligomerization Properties

From computational simulations of a serotonin 2A receptor (5-HT2AR) model complexed with pharmacologically and structurally diverse ligands we identify different conformational states and dynamics adopted by the receptor bound to the full agonist 5-HT, the partial agonist LSD, and the inverse agonist Ketanserin. The results from the unbiased all-atom molecular dynamics (MD) simulations show that the three ligands affect differently the known GPCR activation elements including the toggle switch at W6.48, the changes in the ionic lock between E6.30 and R3.50 of the DRY motif in TM3, and the dynamics of the NPxxY motif in TM7. The computational results uncover a sequence of steps connecting these experimentally-identified elements of GPCR activation. The differences among the properties of the receptor molecule interacting with the ligands correlate with their distinct pharmacological properties. Combining these results with quantitative analysis of membrane deformation obtained with our new method (Mondal et al, Biophysical Journal 2011), we show that distinct conformational rearrangements produced by the three ligands also elicit different responses in the surrounding membrane. The differential reorganization of the receptor environment is reflected in (i)-the involvement of cholesterol in the activation of the 5-HT2AR, and (ii)-different extents and patterns of membrane deformations. These findings are discussed in the context of their likely functional consequences and a predicted mechanism of ligand-specific GPCR oligomerization