Multistructure 3D-QSAR studies on a series of conformationally constrained butyrophenones docked into a new homology model of the 5-HT2A receptor

Molecular Docking Assessment of Cathinones as 5-HT2AR Ligands: Developing of Predictive Structure-Based Bioactive Conformations and Three-Dimensional Structure-Activity Relationships Models for Future Recognition of Abuse Drugs

Commercially available cathinones are drugs of long-term abuse drugs whose pharmacology is fairly well understood. While their psychedelic effects are associated with 5-HT2AR, the enclosed study summarizes efforts to shed light on the pharmacodynamic profiles, not yet known at the receptor level, using molecular docking and three-dimensional quantitative structure-activity relationship (3-D QSAR) studies. The bioactive conformations of cathinones were modeled by AutoDock Vina and were used to build structure-based (SB) 3-D QSAR models using the Open3DQSAR engine. Graphical inspection of the results led to the depiction of a 3-D structure analysis-activity relationship (SAR) scheme that could be used as a guideline for molecular determinants by which any untested cathinone molecule can be predicted as a potential 5-HT2AR binder prior to experimental evaluation. The obtained models, which showed a good agreement with the chemical properties of co-crystallized 5-HT2AR ligands, proved to be valuable for future virtual screening campaigns to recognize unused cathinones and similar compounds, such as 5-HT2AR ligands, minimizing both time and financial resources for the characterization of their psychedelic effects.

Difluorodiazoethane as a masked acetylene equivalent in formal [3 + 2] cycloadditions with ketones to access 2,3-functionalized furans

A transition-metal-free [3+2] cycloaddition of CF2HCHN2 with β-ketones is reported, which enables the synthesis of 2,3-functionalized furans. Sequential defluorination, nucleophilic addition, and cyclization are key elemental steps of the reaction.

Aminothiaindanone as an Accessible Scaffold for a Three-Point Chemical Diversity

Aminothiaindanone heterocycle appears to be a scaffold of interest in medicinal chemistry. To increase the chemical diversity in this series, the introduction of three-point chemical diversity on the cyclopenta[b]thiophen-4-one scaffold was explored. About thirty newly functionalized thiophene-containing bicycles were obtained using various chemical reactions, paving the way for novel possibilities in medicinal chemistry projects.

Two Binding Geometries for Risperidone in Dopamine D3 Receptors: Insights on the Fast-Off Mechanism through Docking, Quantum Biochemistry, and Molecular Dynamics Simulations

Risperidone is an atypical antipsychotic used in the treatment of schizophrenia and of symptoms of irritability associated with autism spectrum disorder (ASD). Its main action mechanism is the blockade of D2-like receptors acting over positive and negative symptoms of schizophrenia with small risk of extrapyramidal symptoms (EPS) at doses corresponding to low/moderate D2 occupancy. Such a decrease in the side effect incidence can be associated with its fast unbinding from D2 receptors in the nigrostriatal region allowing the recovery of dopamine signaling pathways. We performed docking essays using risperidone and the D3 receptor crystallographic data and results suggested two possible distinct orientations for risperidone at the binding pocket. Orientation 1 is more close to the opening of the binding site and has the 6-fluoro-1,2 benzoxazole fragment toward the bottom of the D3 receptor cleft, while orientation 2 is deeper inside the binding pocket with the same fragment toward to the receptor surface. In order to unveil the implications of these two binding orientations, classical molecular dynamics and quantum biochemistry computations within the density functional theory formalism and the molecular fractionation with conjugate caps framework were performed. Quantum mechanics/molecular mechanics suggests that orientation 2 (considering the contribution of Glu90) is slightly more energetically stable than orientation 1 with the main contribution coming from residue Asp110. The residue Glu90, positioned at the opening of the binding site, is closer to orientation 1 than 2, suggesting that it may have a key role in stability through attractive interaction with risperidone. Therefore, although orientations 1 and 2 are both likely to occur, we suggest that the occurrence of the first may contribute to the reduction of side effects in patients taking risperidone due to the reduction of dopamine receptor occupancy in the nigrostriatal region through a mechanism of fast dissociation. The atypical effect may be obtained simply by either delaying D3R full blockage by spatial hindrance of orientation 1 at the binding site or through an effective blockade followed by orientation 1 fast dissociation. While the molecular interpretation suggested in this work shed some light on the potential molecular mechanisms accounting for the reduced extrapyramidal symptoms observed during risperidone treatment, further studies are necessary in order to evaluate the implications of both orientations during the receptor activation/inhibition. Altogether these data highlight important hot spots in the dopamine receptor binding site bringing relevant information for the development of novel/derivative agents with atypical profile.

An application of two MIFs-based tools (Volsurf+ and Pentacle) to binary QSAR: the case of a palinurin-related data set of non-ATP competitive glycogen synthase kinase 3β (GSK-3β) inhibitors

VolSurf+ and GRIND descriptors extract the information present in MIFs calculated by GRID: the first are simpler to interpret and generally applied to ADME-Tox topics, whereas the latter are more sophisticated and thus more suited for pharmacodynamics events. Here we present a study which compares binary QSAR models obtained with VolSurf+ descriptors and GRIND for a data set of non-ATP competitive GSK-3β inhibitors chemically related to palinurin for which the biological activity is expressed in binary format. Results suggest not only that the simpler Volsurf+ descriptors are good enough to predict and chemically interpret the investigated phenomenon but also a bioactive conformation of palinurin which may guide future design of ATP non-competitive GSK-3 inhibitors.

Affinity of aporphines for the human 5-HT2A receptor: insights from homology modeling and molecular docking studies

Analogs of nantenine were docked into a modeled structure of the human 5-HT(2A) receptor using ICM Pro, GLIDE, and GOLD docking methods. The resultant docking scores were used to correlate with observed in vitro apparent affinity (K(e)) data. The GOLD docking algorithm when used with a homology model of 5-HT(2A), based on a bovine rhodopsin template and built by the program MODELLER, gives results which are most in agreement with the in vitro results. Further analysis of the docking poses among members of a C1 alkyl series of nantenine analogs, indicate that they bind to the receptor in a similar orientation, but differently than nantenine. Besides an important interaction between the protonated nitrogen of the C1 alkyl analogs and residue Asp155, we identified Ser242, Phe234, and Gly238 as key residues responsible for the affinity of these compounds for the 5-HT(2A) receptor. Specifically, the ability of some of these analogs to establish a H-bond with Ser242 and hydrophobic interactions with Phe234 and Gly238 appears to explain their enhanced affinity as compared to nantenine.

The development and validation of a novel virtual screening cascade protocol to identify potential serotonin 5-HT(7)R antagonists

In an attempt to identify new ligands for the 5-HT(7) receptor (5-HT(7)R), we developed and tested a hierarchical multi-step strategy of virtual screening (VS) based on two-dimensional (2D) pharmacophore similarity, physicochemical scalar descriptors, an ADME/Tox filter, three-dimensional (3D) pharmacophore searches and a docking protocol. Six chemical classes of 5-HT(7)R antagonists were used as query structures in a double-path virtual screening scheme. The Enamine screening database, consisting of approximately 730,000 commercially available drug-like compounds, was adopted and used as a source of structures. A biological evaluation of 26 finally selected virtual hits resulted in finding two benzodioxane derivatives with significant affinity (K(i)=197 and 265 nM). The approach described in this case study can be easily used as a general rational drug design tool for other biological targets.

Phe369(7.38) at human 5-HT(7) receptors confers interspecies selectivity to antagonists and partial agonists

BACKGROUND AND PURPOSE

Human and rat 5-HT(7) receptors were studied with a particular emphasis on the molecular interactions involved in ligand binding, searching for an explanation to the interspecies selectivity observed for a set of compounds. We performed affinity studies, molecular modelling and site-directed mutagenesis, with special focus on residue Phe(7.38) of the human 5-HT(7) receptor [Cys(7.38) in rat].

EXPERIMENTAL APPROACH

Competition binding studies were performed for seven 5-HT(7) receptor ligands at three different 5-HT(7) receptors. The functional behaviour was evaluated by measuring 5-carboxytryptamine-stimulated cAMP production. Computational simulations were carried out to explore the structural bases in ligand binding observed for these compounds.

KEY RESULTS

Competition experiments showed a remarkable selectivity for the human receptor when compared with the rat receptor. These results indicate that mutating Cys to Phe at position 7.38 profoundly affects the binding affinities at the 5-HT(7) receptor. Computational simulations provide a structural interpretation for this key finding. Pharmacological characterization of compounds mr25020, mr25040 and mr25053 revealed a competitive antagonistic behaviour. Compounds mr22423, mr22433, mr23284 and mr25052 behaved as partial agonists.

CONCLUSIONS AND IMPLICATIONS

We propose that the interspecies difference in binding affinities observed for the compounds at human and rat 5-HT(7) receptors is due to the nature of the residue at position 7.38. Our molecular modelling simulations suggest that Phe(7.38) in the human receptor is integrated in the hydrophobic pocket in the central part of the binding site [Phe(6.51)-Phe(6.52)] and allows a tighter binding of the ligands when compared with the rat receptor.

Cyndi: a multi-objective evolution algorithm based method for bioactive molecular conformational generation

BACKGROUND

Conformation generation is a ubiquitous problem in molecule modelling. Many applications require sampling the broad molecular conformational space or perceiving the bioactive conformers to ensure success. Numerous in silico methods have been proposed in an attempt to resolve the problem, ranging from deterministic to non-deterministic and systemic to stochastic ones. In this work, we described an efficient conformation sampling method named Cyndi, which is based on multi-objective evolution algorithm.

RESULTS

The conformational perturbation is subjected to evolutionary operation on the genome encoded with dihedral torsions. Various objectives are designated to render the generated Pareto optimal conformers to be energy-favoured as well as evenly scattered across the conformational space. An optional objective concerning the degree of molecular extension is added to achieve geometrically extended or compact conformations which have been observed to impact the molecular bioactivity (J Comput -Aided Mol Des 2002, 16: 105-112). Testing the performance of Cyndi against a test set consisting of 329 small molecules reveals an average minimum RMSD of 0.864 A to corresponding bioactive conformations, indicating Cyndi is highly competitive against other conformation generation methods. Meanwhile, the high-speed performance (0.49 +/- 0.18 seconds per molecule) renders Cyndi to be a practical toolkit for conformational database preparation and facilitates subsequent pharmacophore mapping or rigid docking. The copy of precompiled executable of Cyndi and the test set molecules in mol2 format are accessible in Additional file 1.

CONCLUSION

On the basis of MOEA algorithm, we present a new, highly efficient conformation generation method, Cyndi, and report the results of validation and performance studies comparing with other four methods. The results reveal that Cyndi is capable of generating geometrically diverse conformers and outperforms other four multiple conformer generators in the case of reproducing the bioactive conformations against 329 structures. The speed advantage indicates Cyndi is a powerful alternative method for extensive conformational sampling and large-scale conformer database preparation.

Potential modes of interaction of 9-aminomethyl-9,10-dihydroanthracene (AMDA) derivatives with the 5-HT2A receptor: a ligand structure-affinity relationship, receptor mutagenesis and receptor modeling investigation

The effects of 3-position substitution of 9-aminomethyl-9,10-dihydroanthracene (AMDA) on 5-HT 2A receptor affinity were determined and compared to a parallel series of DOB-like 1-(2,5-dimethoxyphenyl)-2-aminopropanes substituted at the 4-position. The results were interpreted within the context of 5-HT 2A receptor models that suggest that members of the DOB-like series can bind to the receptor in two distinct modes that correlate with the compounds' functional activity. Automated ligand docking and molecular dynamics suggest that all of the AMDA derivatives, the parent of which is a 5-HT 2A antagonist, bind in a fashion analogous to that for the sterically demanding antagonist DOB-like compounds. The failure of the F340 (6.52)L mutation to adversely affect the affinity of AMDA and the 3-bromo derivative is consistent with the proposed modes of orientation. Evaluation of ligand-receptor complex models suggest that a valine/threonine exchange between the 5-HT 2A and D 2 receptors may be the origin of selectivity for AMDA and two substituted derivatives.

Methoxy-substituted 9-aminomethyl-9,10-dihydroanthracene (AMDA) derivatives exhibit differential binding affinities at the 5-HT(2A) receptor

The effects of methoxy-substitution at the 1-, 2-, 3-, and 4-positions of 9-aminomethyl-9,10-dihydroanthracene (AMDA) on h5-HT(2A) receptor affinity were determined. Racemic mixtures of these compounds were found to show the following affinity trend: 3-MeO > 4-MeO > 1-MeO approximately 2-MeO. Comparison of the effects of these substitutions, with the aid of computational molecular modeling techniques, suggest that the various positional and stereochemical isomers of the methoxy-substituted AMDA compounds interact differently with the h5-HT(2A) receptor. It is predicted that for the compounds with higher affinities, the methoxy oxygen atom is able to interact with hydrogen bond-donating sidechains within alternative h5-HT(2A) receptor binding sites, whereas the lower-affinity isomers lack this ability.