Expanding chalcogen bonds in thiophenes to interactions with halogens.

Compounds containing the thiophene moiety find several applications in physics and chemistry, such as electrical conduction, which depends on specific conformations to properly exhibiting the desired properties. In turn, chalcogen bonding has found to modulate the conformation of some N-thiophen-2-ylfomamides. Since halogens participate in a kin interaction (halogen bonding) and are abundant in agrochemicals, pharmaceuticals, and materials, we have quantum-chemically explored the interaction between organic halogen and thiophene as a conformational modulator in some model compounds. Although such interaction indeed appears, as demonstrated by atoms in molecules and natural bond orbital analysis, it is inefficient to control the conformational equilibrium. An energy decomposition analysis scheme demonstrated that halomethane and thiophene tend to move away from one another due to a core component (Pauli repulsion and exchange), which is mainly due to a deformation term. Therefore, chalcogen bonds with halogens appear weaker than with other chalcogens.

Fluorinated benzoxazinones designed via MIA-QSAR, docking and molecular dynamics as protoporphyrinogen IX oxidase inhibitors.

BACKGROUND: Fluorine plays a significant role in agrochemical science because approximately 25% of herbicides licensed worldwide contain this element. In a pool of previously synthesized benzoxazinones, some compounds contained fluorine and demonstrated inhibitory activities against protoporphyrinogen IX oxidase (PPO). Therefore, three data sets of benzoxazinone derivatives with known inhibitory activity against PPO were employed to build a multivariate image analysis applied to a quantitative structure-activity relationships (MIA-QSAR) model to identify improved analogs with at least one fluorine substituent. RESULTS: The QSAR model was vigorously validated and demonstrated to be highly predictive (r2 = 0.85, q2 = 0.71, and r2 pred = 0.88); thus, the model can provide reliable estimations for the PPO inhibitory activity of unknown derivatives. From these compounds, a couple of N-substituted benzoxazinones that contained the -CH2 CHF2 group were found with predicted pKi values larger than 8 (Ki in mol L-1 ) and higher lipophilicity than the most active data set compounds. In addition, we carried out a systematic investigation of the binding mode of PPO by performing computational docking followed by molecular dynamics simulations. The proposed binding mode was consistent with experimental studies, and several potential key residues were identified. CONCLUSION: Two new proposed benzoxazinones exhibited better performance than compounds of the data set, and fluorine substituents played pivotal roles in describing the biological activities. © 2024 Society of Chemical Industry.

Multivariate image analysis applied to quantitative structure-activity relationships and docking studies of recent hydroxyphenylpyruvate deoxygenase inhibitors.

BACKGROUND: Mesotrione is a triketone widely used as an inhibitor of the hydroxyphenylpyruvate deoxygenase (HPPD) enzyme. However, new agrochemicals should be developed continuously to tackle the problem of herbicide resistance. Two sets of mesotrione analogs have been synthesized recently and they have demonstrated successful phytotoxicity against weeds. In this study, these compounds were joined to form a single data set and the HPPD inhibition of this enlarged library of triketones was modeled using multivariate image analysis applied to quantitative structure-activity relationships (MIA-QSAR). Docking studies were also carried out to validate the MIA-QSAR findings and to aid the interpretation of ligand-enzyme interactions responsible for the bioactivity (pIC50 ). RESULTS: The MIA-QSAR models based on van der Waals radii (rvdW ), electronegativity (ε), and the rvdW /ε ratio as molecular descriptors were both predictive to an acceptable degree (r2 ≥ 0.80, q2 ≥ 0.68 and r2 pred ≥ 0.68). Subsequently, partial least squares (PLS) regression parameters were applied to predict the pIC50 values of newly proposed derivatives, yielding a few promising agrochemical candidates. The calculated log P for most of these derivatives was found to be higher than that of mesotrione and the library compounds, indicating that they should be less prone to leach out and contaminate groundwater. CONCLUSION: Multivariate image analysis descriptors corroborated by docking studies were capable of modeling the herbicidal activities of 68 triketones reliably. Due to the substituent effects at the triketone framework, particularly of a nitro group in R3 , promising analogs could be designed. The P9 proposal demonstrated higher calculated activity and log P than commercial mesotrione. © 2023 Society of Chemical Industry.

Synergistic effect of intra- and intermolecular hydrogen bond in 2-haloethanols probed by infrared.

Fluorine is often considered the only halogen to effectively engage in hydrogen bonds, while the other halogens, particularly iodine, are not electronegative enough to participate as hydrogen bond acceptors in electrostatic interactions. 2-Fluoroethanol and 2-iodoethanol have been studied herein to test this assumption, since a highly stable gauche conformation can experience the intramolecular hydrogen bond. However, the infrared O  H stretching frequency indicates that the hydroxyl group in 2-fluoroethanol is not engaged in intramolecular hydrogen bond, while the corresponding vibration mode for 2-iodoethanol suggests that not only the O  H is engaged in such interaction, but also that intramolecular hydrogen bond may drive the conformational equilibrium in this molecule. Theoretical calculations support the covalent nature of this interaction, and provide evidence that intermolecular hydrogen bond with a water molecule, and probably with the polar solvents tested experimentally, occurs with the hydroxyl rather than with the iodine substituent, as conventionally, in order to keep the intramolecular hydrogen bond effective.

MIA-QSAR study of the structural merging of (thio)benzamide herbicides with photosynthetic system II inhibitory activities.

Benzamide herbicides consist of a class of photosynthetic system II (PSII) inhibitors widely used for weed control. However, the development of resistance by these weeds to the known herbicides requires an ongoing search for new agrochemicals. We report the combination of two congeneric series of (thio)benzamide herbicides into a single data set and subsequent modeling of their herbicidal activities against PSII using MIA-QSAR. The robust and predictive models were used to estimate the pIC50 of new agrochemical candidates, which were proposed based on a chemical mixing of the substructures of the most active compounds present in the data set. The chemical features affecting the herbicidal activities were analyzed using MIA contour maps, whereas the ligand-enzyme interactions responsible for the binding affinities were rationalized through docking studies. The proposed compound possessing a thiobenzamide moiety and C-11 chain, H, NO2, OH, and OH as variable substituents was the most promising alternative.Communicated by Ramaswamy H. Sarma.

Theoretical Exploitation of 1,2,3,4,5,6-Hexachloro- and 1,2,3,4,5,6-Hexafluorocyclohexane Isomers as Biologically Active Compounds.

Hexachlorocyclohexanes (HCHs) have been widely explored as biological compounds during the last century. However, most of them were banned due to their potential toxicity in humans, animals, and the environment. Revisiting HCHs to explore their biological activity while improving key features is valuable and may lead to a new class of pesticides that utilizes the biological response of HCHs without their toxic characteristics. In this sense, the fluorine atom can be a possible alternative since a large number of therapeutics and agrochemicals have been developed with this halogen in their structure. We have evaluated herein the conformational behavior of HCHs and their bioisosteric fluorinated compounds, namely, hexafluorocyclohexanes (HFHs), through quantum-chemical calculations. We also explored the potential of the HCH and HFH isomers as biological compounds by docking them inside three possible targets. It was demonstrated that HCH and HFH have similar ligand-protein interactions with three pockets: the picrotoxin and barbiturate sites of the GABAA receptor and the ryanodine receptor. The results support HFHs as possible alternatives for HCHs since the replacement of Cl with F does not forfeit the main ligand-protein interactions. Finally, we demonstrated that HFHs have a lower log P than HCHs by almost two logarithmic units. This result highlights the role of fluorine in distribution and bioaccumulation.

Iodine Gauche Effect Induced by an Intramolecular Hydrogen Bond.

The gauche conformer in 1-X,2-Y-disubstituted ethanes, that is, the staggered orientation in which X and Y are in closer contact, is only favored for relatively small substituents that do not give rise to large X···Y steric repulsion. For more diffuse substituents, weakly attractive orbital interactions between antiperiplanar bonds (i.e., hyperconjugation) cannot overrule the repulsive forces between X and Y. Our quantum chemical analyses of the rotational isomerism of XCH2CH2Y (X = F, OH; Y = I) at ZORA-BP86-D3(BJ)/QZ4P reveal that indeed the anti conformer is generally favored due to a less destabilizing I···F and I···O-H steric repulsion. The only case when the gauche conformer is preferred is when the hydroxyl hydrogen is oriented toward the iodine atom in the 2-iodoethanol. This is because of the significantly stabilizing covalent component of the I···H-O intramolecular hydrogen bond. Therefore, we show that strong intramolecular interactions can overcome the steric repulsion between bulky substituents in 1,2-disubstituted ethanes and cause the gauche effect. Our quantum chemical computations have guided nuclear magnetic resonance experiments that confirm the increase in the gauche population as X goes from F to OH.

In Silico Interactions of the Components from the Schinus terebinthifolius Extract with Human Tyrosinase.

The anti-tyrosinase activity of the leaf extract of Schinus terebinthifolius, also known as Brazilian peppertree, was evaluated using multiple in silico approaches, such as molecular homology, molecular docking, MM-GBSA, molecular dynamics, MM-PBSA, QSAR, and skin permeability predictions. With these computational tools, the compounds that downregulate tyrosinase enzyme activity could be evaluated, and more potent molecules could be identified. The results indicated that various compounds, especially luteolin, are accountable for the anti-tyrosinase activity of S. terebinthifolius. For cosmetic application, further studies with luteolin are especially recommended, for having presented a good performance both in theoretical inhibition (30.92 kJ mol-1 ) and skin permeability (LogKp=-6.62 cm-1 ).

Computer-Assisted Proposition of Promising Aryloxyacetic Acid Derivatives as HPPD Inhibitors.

A series of aryloxyacetic acid derivatives have demonstrated promising herbicidal performance by inhibition of the hydroxyphenylpyruvate deoxygenase (HPPD) enzyme. We hereby applied quantitative structure-activity relationship (QSAR) and docking strategies to model and chemically understand the bioactivities of these compounds and subsequently propose unprecedented analogues aiming at improving the herbicidal and environmental properties. Bulky halogens at the 2-, 3-, 4-, and 6-positions of an aromatic ring, CF3 in 4-position, and the 2-NO2 group in a phenyl ring appear to favor the HPPD inhibition. At the same time, Me and OMe substituents contribute to decreasing the pKi values. Accordingly, a few compounds were proposed and the candidate with 2,4,6-triBr substituents demonstrated an estimated pKi similar to those of the best library compounds. This finding was corroborated by the docking scores of the ligand-enzyme interactions. In addition, the high calculated lipophilicity of some proposed agrochemicals suggests that they should have low soil mobility and, therefore, are not prone to easily leach out and reach groundwater, despite causing other ecological issues.

In silico modeling of the AHAS inhibition of an augmented series of pyrimidine herbicides and design of novel derivatives.

Pyrimidine compounds comprise a class of acetohydroxyacid synthase (AHAS) inhibitors, thus possessing herbicidal activity. Due to the ongoing development of resistance by weeds to current herbicides, the design of new agrochemical candidates is often required. This work reports the proposition of unprecedented pyrimidines as herbicides guided by quantitative structure-activity relationship (QSAR) modeling. Multivariate image analysis (MIA) descriptors for 66 pyrimidine derivatives obtained from different sources were regressed against inhibitory activity data, and the resulting QSAR models were found to be reliable and predictive (r2 = 0.88 ± 0.07, q2 = 0.53 ± 0.06, and r2pred = 0.51 ± 0.10 in a bootstrapping experiment using electronegativity-based descriptors). The chemical features responsible for the herbicidal activities were analyzed through MIA contour maps that describe the substituent effects on the response variables, whereas the interaction between the proposed compounds and AHAS was analyzed through docking studies. From the proposed compounds, at least five pyrimidine derivatives exhibited promising performance as AHAS inhibitors compared to the known analogs.

Is conformation relevant for QSAR purposes? 2D Chemical representation in a 3D-QSAR perspective.

Conformation has a key role in the mechanism of interaction between small molecules and biological receptors. However, encoding this type of information in molecular descriptors for the construction of robust quantitative structure-activity relationships (QSAR) models is not an easy task and, so far, the dependence of these models on such feature has not been thoroughly investigated. In the present study, the authors explore the effects of conformational information on a 3D-QSAR technique by comparing models built with descriptors that encode fully described tridimensional aspects (structures docked inside a biological target), with descriptors in which this information is suppressed (flat structures) or not fully described (structures with quantum-chemically optimized geometries). As a result, the validation parameters indicate that the robustness of the models seems to be more related to the alignment aspect of the structures than to how well their tridimensional features are described.

Computer-Assisted Improvement of Sulfonylureas with Antifungal Properties and Limited Herbicidal Activity: Potential Application in Forage Conservation.

This work reports studies at the molecular level of a series of modified sulfonylureas to determine the chemophoric sites responsible for their antifungal and herbicidal activities. For forage conservation, high antifungal potency and low phytotoxicity are required. A molecular modeling study based on multivariate image analysis applied to quantitative structure-activity relationship (MIA-QSAR) was performed to model these properties, as well as to guide the design of new agrochemical candidates. As a result, the MIA-QSAR models were reliable, robust, and predictive; for antifungal activity, the averages of the main validation parameters were r2 = 0.936, q2 = 0.741, and r2pred = 0.720, and for herbicidal activity, the model was very predictive (r2pred = 0.981 and r2m = 0.944). From the interpretation of the MIA-plots, 46 novel sulfonylureas with likely improved performance were proposed, from which 9 presented promising calculated selectivity indexes. Docking studies were performed to validate the QSAR predictions and to understand the interaction mode of the proposed ligands with the acetohydroxyacid synthase enzyme.

QSAR-Guided Proposition of N-(4-methanesulfonyl)Benzoyl-N'-(Pyrimidin-2-yl)Thioureas as Effective and Safer Herbicides.

Chlorinated agrochemicals play a major role in toxicity due especially to the labile C - Cl bond and high lipophilicity of organochlorines. In turn, urea and thiourea herbicides are widely used for weed control. A series of substituted N-benzoyl-N'-pyrimidin-2-yl thioureas has been recently synthesized and tested against Brassica napus L., demonstrating promising herbicidal activities, particularly for chlorinated derivatives. We have therefore modeled these activities using multivariate image analysis applied to quantitative structure-activity relationships (MIA-QSAR) to find out a significant and reliable correlation between measured and predicted inhibition of B. napus L. root growth (%) and, ultimately, to propose effective, non-chlorinated and/or less lipophilic N-(4-methanesulfonyl)benzoyl-N'-(pyrimidin-2-yl)thiourea candidates. The model was found to be predictive, giving an average r2pred in the external validation of 0.833. The predicted data for the proposed herbicides, interpreted in terms of MIA-plots of the chemical moieties responsible for bioactivity and supported by docking studies towards the photosystem II enzyme, suggest that substituents at both R1 and R2 positions modulate the agrochemical (R1 = Cl increases and R2 = OR decreases bioactivity) and environmental friendship (particularly with R2 = OH) performances of this class of compounds.

Theoretical and X-Ray Evidence of Electrostatic Phosphonium anti and gauche Effects.

Sulphur, not phosphorus, is the only known third-row element capable of experiencing an electrostatic gauche effect with fluorine. Some six-membered rings containing an endocyclic phosphorus atom and a ß-fluorine substituent that can interconvert to axial (gauche relative to phosphorus) and equatorial positions were then analysed. While phosphines do not establish an electrostatic attraction between fluorine and phosphorus, some oxidised forms exhibit surprising stability for the sterically disfavoured axial orientation. Because the nature of this behaviour was not obvious, since an intramolecular hydrogen bond can appear, a phosphonium derivative was further studied and its axial conformation was found to be highly stable. A preference for the gauche arrangement appears even for the acyclic and sterically hindered (2-fluoroethyl)triphenylphosphonium cation. On the other hand, (ethane-1,2-diyl)bis(phosphonium) cations are exclusively in anti conformation due to an (+/+)-electrostatic repulsion between the positively charged phosphonium groups.

Fluorinated dihydropyridines as candidates to block L-type voltage-dependent calcium channels.

Voltage-gated calcium (Cav) channels malfunction may lead to Alzheimer's and cardiovascular disorders, thus a critical protein target for drug development and treatment against several diseases. Indeed, dihydropyridines (DHPs) as nifedipine and amlodipine are top-selling pharmaceuticals and, respectively, the 121st and 5th most prescribed drugs in the United States that have been used as successful selective blockers for L-type Ca2+ channels (LCC) and may be helpful model structures to compare with new DHP analogs. In this context, we have performed a structure-based drug design (SBDD) study of several fluorinated DHPs by using homology modeling, molecular docking, quantitative structure activity relationship (QSAR) and molecular dynamics calculations. Such approaches combined with molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) interaction energy results and screening of ADMET (absorption, distribution, metabolism, excretion and toxicity) properties indicate that all ligands in this study are potential new candidates to be tested experimentally for inhibition of LCC and may have higher affinities than the commonly used drugs, being convenient synthetic routes proposed for 11-16, which are among the ligands that showed the best theoretical results concerning LCC inhibition. Furthermore, the ligand interactions with the binding site were carefully examined using the topological noncovalent interactions (NCI) method, which highlighted specifically responsible amino acid residues that increase the spontaneity of the new proposed DHP ligands.Communicated by Ramaswamy H. Sarma.

An examination of the relationship between molecular dipole moment and blood-gas partition for common anaesthetic gases.

The solubility of inhalational anaesthetics in the bloodstream is related to the minimum alveolar concentration (MAC), which is an indicator of anaesthetic gas potency. The blood-gas partition coefficient (Kbg) is a measure of how much anaesthetics bind to plasma proteins in the blood compared to air. Just like the octanol-water partition coefficient, the Kbg may be related to the molecular dipole moment (µ), which can be modulated by the molecular conformation. Our quantum-chemical calculations demonstrated that subtle stereoelectronic interactions, namely those responsible for the anomeric and gauche effects, control the conformational equilibrium of some widely used volatile fluorinated anaesthetics and, consequently, of their molecular dipole moments. Since a remarkable correlation between empirical Kbg and calculated µ was found for these anaesthetics, the average molecular dipole moments may be used to predict the anaesthetic gas potency and other properties, such as lipid solubility, of inhalational anaesthetic analogs.

The Gauche Effect in XCH2 CH2 X Revisited.

We have quantum chemically investigated the rotational isomerism of 1,2-dihaloethanes XCH2 CH2 X (X = F, Cl, Br, I) at ZORA-BP86-D3(BJ)/QZ4P. Our Kohn-Sham molecular orbital (KS-MO) analyses reveal that hyperconjugative orbital interactions favor the gauche conformation in all cases (X = F-I), not only for X = F as in the current model of this so-called gauche effect. We show that, instead, it is the interplay of hyperconjugation with Pauli repulsion between lone-pair-type orbitals on the halogen substituents that constitutes the causal mechanism for the gauche effect. Thus, only in the case of the relatively small fluorine atoms, steric Pauli repulsion is too weak to overrule the gauche preference of the hyperconjugative orbital interactions. For the larger halogens, X⋅⋅⋅X steric Pauli repulsion becomes sufficiently destabilizing to shift the energetic preference from gauche to anti, despite the opposite preference of hyperconjugation.

Conformational preferences of fluorine-containing agrochemicals and their implications for lipophilicity prediction.

Molecular polarity governs lipophilicity, which in turn determines important agrochemical and environmental properties, such as soil sorption and bioconcentration of organic compounds. Since the C-F bond is the most polar in organic chemistry, the orientation of fluorine substituents originating from the rotation around C-C(F) bonds should affect the polarity and, consequently, the physicochemical and biological properties of fluorine-containing agrochemicals. Accordingly, this study aims to determine the most likely conformers of some fluorine-containing agrochemicals and to correlate their molecular dipole moments with the respective n-octanol/water partition coefficients (log P), in order to investigate the dependence of the lipophilicity with the molecular conformation.

Theoretical study of fluorinated bioisosteres of organochlorine compounds as effective and eco-friendly pesticides.

Chlordane is a worldwide banned organochlorine insecticide because of its hazard to animal and human health. It is also a persistent organic pollutant, which can affect either the soil or the aquatic life. The same applies to other chlorinated cyclodiene insecticides, such as dieldrin and aldrin. In turn, organofluorine compounds have a widespread use in agriculture. Therefore, density functional calculations and docking studies showed that the bioisosteric replacement of chlorines in the above-mentioned compounds by fluorines improves some physicochemical parameters used to estimate the toxicity and environmental risk of these compounds, as well as the ligand-enzyme (GABAA receptor-chloride channel complex) interactions related to their insecticidal activity. This work is an effort to provide an improved new class of organofluorine pesticides.

Bioconformational modulation of a thymidine kinase enzyme ligand through F⋯HO intramolecular hydrogen bond.

While the induced-fit of a ligand towards an enzyme is pivotally dictated by intermolecular hydrogen bonds between the small molecule and amino acid residues in the binding site, the role of intramolecular hydrogen bond as contributing interaction for a bioactive conformation is not well understood. This work reports a theoretical conformational analysis of a thymidine kinase enzyme ligand (NMF) that is prone to experience an F⋯HO intramolecular hydrogen bond, inside and outside the biological binding site. This interaction stabilizes the most favorable conformations of NMF in the gas phase and, although it can be disrupted in a biological environment due to intermolecular hydrogen bonds in some cases, these interactions are competitive in other systems. Therefore, an intramolecular hydrogen bond can affect the conformational likeliness most related to the bioactivity. Moreover, isolated conformations governed by this interaction cannot be unequivocally used to generate molecular descriptors in 3D-QSAR (Quantitative Structure-Activity Relationships), since the bioactive conformation may not be determined only by intramolecular interactions.