Conformational analysis of methylphenidate: comparison of molecular orbital and molecular mechanics methods

Synthesis, herbicidal activity and soil degradation of novel 5-substituted sulfonylureas as AHAS inhibitors

BACKGROUND

Chlorsulfuron, metsulfuron-methyl and ethametsulfuron can damage sensitive crops in rotation pattern as a result of their long persistence in soil. To explore novel sulfonylurea (SU) herbicides with favorable soil degradation rates, four series of SUs were synthesized through a structure-based drug design (SBDD) strategy.

RESULTS

The target compounds, especially Ia, Id and Ie, exhibited prospective herbicidal activity against dicotyledon oil seed rape (Brassica campestris), amaranth (Amaranthus retroflexus), monocotyledon barnyard grass (Echinochloa crusgalli) and crab grass (Digitaria sanguinalis) at a concentration of 15 a.i. g ha^-1 . Additionally, Ia, Id and Ig displayed excellent inhibitory effects against AtAHAS, with Kapp i values of 59.1, 34.5 and 71.8 μm, respectively, which were much lower than that of chlorsulfuron at 149.4 μm. The π-π stack and H-bonds between the Ia conformation and AtAHAS in the molecular docking results confirmed the series of compounds to be conventional AHAS inhibitors. In alkaline soil (pH = 8.46), compounds Ia-Ig revealed various degrees of acceleration in the degradation rate compared with chlorsulfuron. Besides, compound Ia showed considerable wheat and corn safety under postemergence at the concentration of 30, 60 and even 120 a.i. g ha^-1 .

CONCLUSION

Overall, based on the synthetic procedure, herbicidal activity, soil degradation and crop safety, the Ia sulfonylureas series were chosen to be investigated as prospective AHAS inhibitors. The 5-dimethylamino group on SUs accelerated the degradation rate at different levels in alkaline soils which seems to be controllable in conventional cropping systems in their further application. © 2022 Society of Chemical Industry.

Controllable Soil Degradation Rate of 5-Substituted Sulfonylurea Herbicides as Novel AHAS Inhibitors

Chlorsulfuron has been applied in wheat fields as a recognized herbicide worldwide, yet it was officially banned in China since 2014 for its soil persistence problem. On the basis of our previous research that 5-dimethylamino distinctively accelerated degradation rate in soils, a modified amino moiety (Ia-c) and monosubstituted amino group (Id-e) were introduced onto the fifth position of the benzene ring in sulfonylurea structures, as well as heterocyclic amino substituents (If-g) to seek a suitable soil degradation rate during such an in situ crop rotation system. Referring to the biological data and ScAHAS inhibition and ScAHAS docking results, they turned out to be AHAS inhibitors with high potent herbicidal activities. The various influence on soil degradation rate along with crop safety indicated that different substituents on the fifth position have exerted an apparent impact. Their united study of structure-activity-safety-degradation relationship has great potential to provide valuable information for further development of eco-friendly agrochemicals.

Cholesteryl ester transfer protein inhibitors in coronary heart disease: Validated comparative QSAR modeling of N, N-disubstituted trifluoro-3-amino-2-propanols

Cholesteryl ester transfer protein (CETP) converts high density lipoprotein cholesterol to low density lipoproteins. It is a promising target for treatment of coronary heart disease. Two dimensional quantitative structure activity relationship (2D-QSAR), hologram QSAR (HQSAR) studies and comparative molecular field analysis (CoMFA) as well as comparative molecular similarity analysis (CoMSIA) were performed on 104 CETP inhibitors. The statistical qualities of generated models were justified by internal and external validation, i.e., q(2) and R(2)pred respectively. The best 2D-QSAR model was obtained with q(2) and R(2)pred values of 0.794 and 0.796 respectively. The 2D-QSAR study suggests that unsaturation, branching and van der Waals volumes may play important roles. The HQSAR model showed q(2) and R(2)pred values of 0.628 and 0.550 respectively. Similarly, CoMFA model showed q(2) and R(2)pred values of 0.707 and 0.755 respectively whereas CoMSIA model was obtained with q(2) and R(2)pred values of 0.696 and 0.703 respectively. CoMFA and CoMSIA studies indicate that steric factors are important at substituted phenoxy and tetrafluoroethoxy groups whereas electropositive factors play important role at difluoromethyl group. The results of 3D-QSAR studies validate those of 2D-QSAR and HQSAR studies as well as the earlier observed SAR data. Current work may help to develop better CETP inhibitors.

Conformational analysis of piperazine and piperidine analogs of GBR 12909: stochastic approach to evaluating the effects of force fields and solvent

Analogs of the flexible dopamine reuptake inhibitor, GBR 12909 (1), may have potential utility in the treatment of cocaine abuse. As a first step in the 3D-QSAR modeling of the dopamine transporter (DAT)/serotonin transporter (SERT) selectivity of these compounds, we carried out conformational analyses of two analogs of 1: a piperazine (2) and a related piperidine (3). Ensembles of conformers consisting of local minima on the potential energy surface of the molecule were generated in the vacuum phase and in implicit solvent by random search conformational analysis using the Tripos and MMFF94 force fields. Some differences were noted in the conformer populations due to differences in the treatment of the tertiary amine nitrogen and ether oxygen atom types by the force fields. The force fields also differed in their descriptions of internal rotation around the C(sp³)-O(sp³) bond proximal to the bisphenyl moiety. Molecular orbital calculations at the HF/6-31G(d) and B3LYP/6-31G(d) levels of C-O internal rotation in model compound (5), designed to model the effect of the proximity of the bisphenyl group on C-O internal rotation, showed a broad region of low energy between -60° to 60° with minima at both -60° and 30° and a low rotational barrier at 0°, in closer agreement with the MMFF94 results than the Tripos results. Molecular mechanics calculations on model compound (6) showed that the MMFF94 force field was much more sensitive than the Tripos force field to the effects of the bisphenyl moiety on C-O internal rotation.

Quantitative conformationally sampled pharmacophore for delta opioid ligands: reevaluation of hydrophobic moieties essential for biological activity

Recent studies have indicated several therapeutic applications for delta opioid agonists and antagonists. To exploit the therapeutic potential of delta opioids developing a structural basis for the activity of ligands at the delta opioid receptor is essential. The conformationally sampled pharmacophore (CSP) method (Bernard et al. J. Am. Chem. Soc. 2003, 125, 3103-3107) is extended here to obtain quantitative models of delta opioid ligand efficacy and affinity. Quantification is performed via overlap integrals of the conformational space sampled by ligands with respect to a reference compound. Iterative refinement of the CSP model identified hydrophobic groups other than the traditional phenylalanine residues as important for efficacy and affinity in DSLET and ICI 174 864. The obtained models for a structurally diverse set of peptidic and nonpeptidic delta opioid ligands offer good predictions with R2 values>0.9, and the predicted efficacy for a set of test compounds was consistent with the experimental values.

DAT/SERT selectivity of flexible GBR 12909 analogs modeled using 3D-QSAR methods

The dopamine reuptake inhibitor GBR 12909 (1-{2-[bis(4-fluorophenyl)methoxy]ethyl}-4-(3-phenylpropyl)piperazine, 1) and its analogs have been developed as tools to test the hypothesis that selective dopamine transporter (DAT) inhibitors will be useful therapeutics for cocaine addiction. This 3D-QSAR study focuses on the effect of substitutions in the phenylpropyl region of 1. CoMFA and CoMSIA techniques were used to determine a predictive and stable model for the DAT/serotonin transporter (SERT) selectivity (represented by pK(i) (DAT/SERT)) of a set of flexible analogs of 1, most of which have eight rotatable bonds. In the absence of a rigid analog to use as a 3D-QSAR template, six conformational families of analogs were constructed from six pairs of piperazine and piperidine template conformers identified by hierarchical clustering as representative molecular conformations. Three models stable to y-value scrambling were identified after a comprehensive CoMFA and CoMSIA survey with Region Focusing. Test set correlation validation led to an acceptable model, with q(2)=0.508, standard error of prediction=0.601, two components, r(2)=0.685, standard error of estimate=0.481, F value=39, percent steric contribution=65, and percent electrostatic contribution=35. A CoMFA contour map identified areas of the molecule that affect pK(i) (DAT/SERT). This work outlines a protocol for deriving a stable and predictive model of the biological activity of a set of very flexible molecules.

Hierarchical clustering analysis of flexible GBR 12909 dialkyl piperazine and piperidine analogs

Pharmacophore modeling of large, drug-like molecules, such as the dopamine reuptake inhibitor GBR 12909, is complicated by their flexibility. A comprehensive hierarchical clustering study of two GBR 12909 analogs was performed to identify representative conformers for input to three-dimensional quantitative structure-activity relationship studies of closely-related analogs. Two data sets of more than 700 conformers each produced by random search conformational analysis of a piperazine and a piperidine GBR 12909 analog were studied. Several clustering studies were carried out based on different feature sets that include the important pharmacophore elements. The distance maps, the plot of the effective number of clusters versus actual number of clusters, and the novel derived clustering statistic, percentage change in the effective number of clusters, were shown to be useful in determining the appropriate clustering level. Six clusters were chosen for each analog, each representing a different region of the torsional angle space that determines the relative orientation of the pharmacophore elements. Conformers of each cluster that are representative of these regions were identified and compared for each analog. This study illustrates the utility of using hierarchical clustering for the classification of conformers of highly flexible molecules in terms of the three-dimensional spatial orientation of key pharmacophore elements.

Singular value decomposition of torsional angles of analogs of the dopamine reuptake inhibitor GBR 12909

Analysis of large, flexible molecules, such as the dopamine reuptake inhibitor GBR 12909 (1), is complicated by the fact that they can take on a wide range of closely related conformations. The first step in the analysis is to classify the conformers into groups. Here, Singular Value Decomposition (SVD) was used to group conformations of GBR 12909 analogs by the similarity of their nonring torsional angles. The significance of the present work, the first application of SVD to the analysis of very flexible molecules, lies in the development of a novel scaling technique for circular data and in the grouping of molecular conformations using a technique that is independent of molecular alignment. Over 700 conformers each of a piperazine (2) and piperidine (3) analog of 1 were studied. Analysis of the score and loading plots showed that the conformers of 2 separate into three large groups due to torsional angles on the naphthalene side of the molecule, whereas those of 3 separate into nine groups due to torsional angles on the bisphenyl side of the molecule. These differences are due to nitrogen inversion at the unprotonated piperazinyl nitrogen of 2, which results in a different ensemble of conformers than those of 3, where no inversion is possible at the corresponding piperidinyl carbon.