Dissecting C-H∙∙∙π and N-H∙∙∙π Interactions in Two Proteins Using a Combined Experimental and Computational Approach

Discovery of a New Class of Lipophilic Pyrimidine-Biphenyl Herbicides Using an Integrated Experimental-Computational Approach

Herbicides are useful tools for managing weeds and promoting food production and sustainable agriculture. In this study, we report on the development of a novel class of lipophilic pyrimidine-biphenyl (PMB) herbicides. Firstly, three PMBs, Ia, IIa, and IIIa, were rationally designed via a scaffold hopping strategy and were determined to inhibit acetohydroxyacid synthase (AHAS). Computational simulation was carried out to investigate the molecular basis for the efficiency of PMBs against AHAS. With a rational binding mode, and the highest in vitro as well as in vivo potency, Ia was identified as a preferable hit. Furthermore, these integrated analyses guided the design of eighteen new PMBs, which were synthesized via a one-step Suzuki-Miyaura cross-coupling reaction. These new PMBs, Iba-ic, were more effective in post-emergence control of grass weeds compared with Ia. Interestingly, six of the PMBs displayed 98-100% inhibition in the control of grass weeds at 750 g ai/ha. Remarkably, Ica exhibited ≥ 80% control against grass weeds at 187.5 g ai/ha. Overall, our comprehensive and systematic investigation revealed that a structurally distinct class of lipophilic PMB herbicides, which pair excellent herbicidal activities with new interactions with AHAS, represent a noteworthy development in the pursuit of sustainable weed control solutions.

Rational Design of TDP-43 Derived α-Helical Peptide Inhibitors: An In Silico Strategy to Prevent TDP-43 Aggregation in Neurodegenerative Disorders

TDP-43, an essential RNA/DNA-binding protein, is central to the pathology of neurodegenerative diseases, such as amyotrophic lateral sclerosis and frontotemporal dementia. Pathological mislocalization and aggregation of TDP-43 disrupt RNA splicing, mRNA stability, and mRNA transport, thereby impairing neuronal function and survival. The formation of amyloid-like TDP-43 filaments is largely facilitated by the destabilization of an α-helical segment within the disordered C-terminal region. In this study, we hypothesized that preventing the destabilization of the α-helical domain could potentially halt the growth of these pathological filaments. To explore this, we utilized a range of in silico techniques to design and evaluate peptide-based therapeutics that bind to pathological TDP-43 amyloid-like filament crystal structures and resist β sheet conversion. Our computational approaches, including biophysical and secondary structure property prediction, molecular docking, 3D structure prediction, and molecular dynamics simulations, were used to assess the structure, stability, and binding affinity of these peptides in relation to pathological TDP-43 filaments. The results of our in silico analyses identified a selection of promising peptides which displayed a stable α-helical structure, exhibited an increased number of intramolecular hydrogen bonds within the helical domain, and demonstrated high binding affinities for pathological TDP-43 amyloid-like filaments. Molecular dynamics simulations provided further support for the structural and thermodynamic stability of these peptides, as they exhibited lower root-mean-square deviation and more favorable free energy landscapes over 300 ns. These findings establish α-helical propensity peptides as potential lead molecules for the development of novel therapeutics against TDP-43 aggregation. This structure-based computational approach for the rational design of peptide inhibitors opens a new direction in the search for effective interventions for ALS, FTD, and other related neurodegenerative diseases. The peptides identified as the most promising candidates in this study are currently subject to further testing and validation through both in vitro and in vivo experiments.

Structure-activity relationship of 2,4-D correlates auxinic activity with the induction of somatic embryogenesis in Arabidopsis thaliana

2,4-dichlorophenoxyacetic acid (2,4-D) is a synthetic analogue of the plant hormone auxin that is commonly used in many in vitro plant regeneration systems, such as somatic embryogenesis (SE). Its effectiveness in inducing SE, compared to the natural auxin indole-3-acetic acid (IAA), has been attributed to the stress triggered by this compound rather than its auxinic activity. However, this hypothesis has never been thoroughly tested. Here we used a library of forty 2,4-D analogues to test the structure-activity relationship with respect to the capacity to induce SE and auxinic activity in Arabidopsis thaliana. Four analogues induced SE as effectively as 2,4-D and 13 analogues induced SE but were less effective. Based on root growth inhibition and auxin response reporter expression, the 2,4-D analogues were classified into different groups, ranging from very active to not active auxin analogues. A halogen at the 4-position of the aromatic ring was important for auxinic activity, whereas a halogen at the 3-position resulted in reduced activity. Moreover, a small substitution at the carboxylate chain was tolerated, as was extending the carboxylate chain with an even number of carbons. The auxinic activity of most 2,4-D analogues was consistent with their simulated TIR1-Aux/IAA coreceptor binding characteristics. A strong correlation was observed between SE induction efficiency and auxinic activity, which is in line with our observation that 2,4-D-induced SE and stress both require TIR1/AFB auxin co-receptor function. Our data indicate that the stress-related effects triggered by 2,4-D and considered important for SE induction are downstream of auxin signalling.

E- or Z-Isomers Arising from the Geometries of Ligands in the Mercury Complex of 2-(Anthracen-9-ylmethylene)-N-phenylhydrazine Carbothioamide

An anionic mercury(II) complex of 2-(anthracen-9-ylmethylene)-N-phenylhydrazine carbothioamide (HATU) and two isomers of a neutral mercury(II) complex of the anion of the same ligand (ATU) were reported. The anionic complex [Hg(HATU)2Cl2]·CH2Cl2 had a monodentate HATU ligand (a neutral form of the ligand) and chloride ligands. The two conformational isomers were of the neutral mercury(II) complex Hg(ATU)2·2DMF. The two isomers were from the E or Z geometry of the ligands across the conjugated C=N-N=C-N scaffold of the coordinated ligand. The two isomers of the complex were independently prepared and characterized. The spectroscopic properties of the isomers in solution were studied by 1H NMR as well as fluorescence spectroscopy. Facile conversion of the E-isomer to the Z-isomer in solution was observed. Density functional theory (DFT) calculations revealed that the Z-isomer of the complex was stable compared to the E-isomer by an energy of 14.35 kJ/mol; whereas, E isomer of the ligand was more stable than Z isomer by 8.37 KJ/mol. The activation barrier for the conversion of the E-isomer to the Z-isomer of the ligand was 167.37 kJ/mol. The role of the mercury ion in the conversion of the E-form to the Z-form was discussed. The mercury complex [Hg(HATU)2Cl2]·CH2Cl2 had the E-form of the ligand. Distinct photophysical features of these mercury complexes were presented.

Crystal structures of MHC class I complexes reveal the elusive intermediate conformations explored during peptide editing

Studies have suggested that MHC class I (MHC I) molecules fluctuate rapidly between numerous conformational states and these motions support peptide sampling. To date, MHC I intermediates are largely uncharacterized experimentally and remain elusive. Here, we present x-ray crystal structures of HLA-B8 loaded with 20mer peptides that show pronounced distortions at the N-terminus of the groove. Long stretches of N-terminal amino acid residues are missing in the electron density maps creating an open-ended groove. Our structures also reveal highly unusual features in MHC I-peptide interaction at the N-terminus of the groove. Molecular dynamics simulations indicate that the complexes have varying degrees of conformational flexibility in a manner consistent with the structures. We suggest that our structures have captured the remarkable molecular dynamics of MHC I-peptide interaction. The visualization of peptide-dependent conformational motions in MHC I is a major step forward in our conceptual understanding of dynamics in high-affinity peptide selection.

A critical comparison of CH⋯π versus π⋯π interactions in the benzene dimer: obtaining benchmarks at the CCSD(T) level and assessing the accuracy of lower scaling methods

We have established CCSD(T)/CBS (Complete Basis Set) limits for 3 stationary points on the benzene dimer potential energy surface, corresponding to the π⋯π (parallel displaced or PD(C_2h), minimum) and CH⋯π (T-shaped or T(C_2v), transition state) and tilted T-shaped (or TT(C_s), minimum) bonding scenarios considering both the structure and binding energy. The CCSD(T)/CBS binding energies are -2.65 ± 0.02 (PD), -2.74 ± 0.03 (T), and -2.83 ± 0.01 kcal mol^-1 (TT). To this end, the CH⋯π is ∼0.2 kcal mol^-1 stronger than the π⋯π interaction, whereas the tilting of the CH donating benzene molecule with respect to the other benzene is worth 0.1 kcal mol^-1. As previously discussed in the literature, the MP2 level of theory does not provide a close match for either the energy or structure, yet the SCS-MP2 yields structures in excellent agreement with respect to the CCSD(T) result. It is found that the SCS-MI-MP2 also gives optimized structures very close to SCS-MP2 (within ∼0.01 Å of the benchmark). Despite the closer match in structure, the spin-biased MP2 methods (SCS-, SCS-MI-, and SOS-MP2) incorrectly predict the relative stabilities of the isomers. That said, none of the spin biased MP2 methods offers a good compromise between energy and structure for the systems examined. Finally, the CCSD(T)/CBS benchmarks were used to assess the performance of 13 DFT functionals selected from different rungs of Jacob's ladder. Several functionals such as TPSS-D3, B3LYP-D3, B97-D, B97-D3, and B2PLYP-D3 provided a good description of the binding energies for both CH⋯π and π⋯π interactions, yielding values within 6% of the CCSD(T)/CBS benchmark values. Unlike the MP2 methods, these functionals correctly predict the relative stability of the PD(C_2h) and T(C_2v) dimers. Further, we find that there is no systematic improvement as Jacob's ladder is ascended (increased complexity of functional). The best functionals that result in a good compromise between structure and energy accuracy are B97-D3 and B2PLYP-D3 for both the CH⋯π and π⋯π interaction. Despite the impressive performance of these functionals, a challenge that remains is ensuring the transferability of these density functionals in accurately describing the interaction between dimers of larger aromatic molecules, the latter requiring high-level benchmarks for these systems.

Elucidating the Hydrotropic Mechanism of the Antagonistic Salt PPh4Cl

PPh₄Cl is an antagonistic salt that recently showed promise as a hydrotropic agent. Here, we give mechanistic insights into the PPh₄Cl-assisted solubility of a dye molecule using molecular dynamics simulations. Our findings reveal that dye molecules aggregate into a cluster which leads to an accumulation of PPh₄⁺ ions in its vicinity and subsequent exclusion of water molecules from the region. The structural organization is attributed to the preferential interaction of dye molecules and PPh₄Cl. The origin of such preference arises from the difference in π-π and CH-π interaction among the pairs. The hydrodynamic radius of PPh₄Cl indicates a low propensity for cluster formation, which enhances its hydrotropic behavior. The process of dye dissolution is thermodynamically favored and occurs through a cooperative mechanism. Our studies provide molecular insight into experimental observations crucial for the design of novel hydrotropes with enhanced solubilizing properties.

Unusual crystal structures of MHC class I complexes reveal the elusive intermediate conformations explored during peptide editing in antigen presentation

Studies have suggested that MHC class I (MHC I) molecules fluctuate rapidly between conformational states as they sample peptides for potential ligands. To date, MHC I intermediates are largely uncharacterized experimentally and remain elusive. We present x-ray crystal structures of HLA-B8 loaded with 20mer peptides that show significant conformational heterogeneity at the N-terminus of the groove. Long stretches of N-terminal residues were missing in the electron density maps creating an unstructured and widely open-ended groove. Our structures also revealed highly unusual features in MHC I and peptide conformations, and in MHC I-peptide interaction at the N-terminus of the groove. Molecular dynamics simulations showed that the complexes have varying degrees of flexibility in a manner consistent with the structures. We suggest that our structures represent transient substates explored by MHC I molecules during peptide editing. The visualization of peptide-dependent conformational flexibility in MHC I groove is a major step forward in our conceptual understanding of peptide repertoire development in antigen presentation. Our study also raises questions about the role of the N-terminus of the groove in peptide editing.

Piezoresistive Fibers with Large Working Factors for Strain Sensing Applications

Piezoresistive fibers with large working factors remain of great interest for strain sensing applications involving large strains, yet difficult to achieve. Here, we produced strain-sensitive fibers with large working factors by dip-coating nanocomposite piezoresistive inks on surface-modified polyether block amide (PEBA) fibers. Surface modification of neat PEBA fibers was carried out with polydopamine (PDA) while nanocomposite conductive inks consisted of styrene-ethylene-butylene-styrene (SEBS) elastomer and carbon black (CB). As such, the deposition of piezoresistive coatings was enabled through nonconventional hydrogen-bonding interactions. The resultant fibers demonstrated well-defined piezoresistive linear relationships, which increased with CB filler loading in SEBS. In addition, gauge factors decreased with increasing CB mass fractions from ∼15 to ∼7. Furthermore, we used the fatigue theory to predict the endurance limit (C_e) of our fibers toward resistance signal stability. Such a piezoresistive performance allowed us to explore the application of our fibers as strain sensors for monitoring the movement of finger joints.

Noise facilitates entrainment of a population of uncoupled limit cycle oscillators

Many biological oscillators share two properties: they are subject to stochastic fluctuations (noise) and they must reliably adjust their period to changing environmental conditions (entrainment). While noise seems to distort the ability of single oscillators to entrain, in populations of uncoupled oscillators noise allows population-level entrainment for a wider range of input amplitudes and periods. Here, we investigate how this effect depends on the noise intensity and the number of oscillators in the population. We have found that, if a population consists of a sufficient number of oscillators, increasing noise intensity leads to faster entrainment after a phase change of the input signal (jet lag) and increases sensitivity to low-amplitude input signals.

Molecular duplexes featuring NH···N, CH···O and CH···π interactions in solid-state self-assembly of triazine-based compounds

Synthetic supramolecular structures constructed through the cooperative action of numerous non-covalent forces are highly desirable as models to unravel and understand the complexity of systems created in nature via self-assembly. Taking advantage of the low cost of 2,4,6-trichloro-1,3,5-triazine (cyanuric chloride) and the sequential nucleophilic substitution reactions with almost all types of nucleophiles, a series of six structurally related novel s-triazine derivatives 1-6 were synthesized and structurally characterized based on their physical, spectral and crystallographic data. The solid-state structures of all the six compounds showed intriguing and unique molecular duplexes featuring NH···N, CH···O and CH···π interactions. Careful analysis of different geometric parameters of the involved H-bonds indicates that they are linear, significant and are therefore responsible for guiding the three-dimensional structure of these compounds in the solid state. The prevalence of sextuple hydrogen bond array-driven molecular duplexes and the possibility of structural modifications on the s-triazine ring render these novel triazine derivatives 1-6 attractive as a platform to create heteroduplex constructs and their subsequent utility in the field of supramolecular chemistry and crystal engineering.

Miniature β-Hairpin Mimetic by Intramolecular Hydrogen Bond and C-H···π Interactions

Canonically, protein β-hairpin motifs are stabilized by intramolecular hydrogen bonds. Here, we attempt to develop a rational design recipe for a miniature hairpin structure stabilized by hydrogen bonding as well as C-H···π interaction and try to understand how such a stabilization effect varies with different functional groups at each terminus. Database analysis shows that the α-amino acids with an aromatic side chain will not favor that kind of C-H···π stabilized hairpin structure. However, hybrid tripeptides with an N-terminal Boc-Trp-Aib corner residue and C-terminal aromatic ω-amino acids fold into the hairpin conformation with a central β-turn/open-turn that is reinforced by a C-H···π interaction. The CCDC database analysis further confirms that this C-H···π stabilized hairpin motif is general for Boc-protected tripeptides containing Aib in the middle and aromatic functionality at the C-terminus. The different α-amino acids like Leu/Ala/Phe/Pro/Ser at the N-terminus have a minor influence on the C-H···π interaction and stabilities of the folded structures in solid-state. However, the hybrid peptides exhibit different degrees of conformational heterogeneity both in the solid and solution phase, which is common for this kind of flexible small molecule. Conformational heterogeneity in the solution phase including the C-H···π stabilized β-hairpin structures are characterized by the molecular dynamics (MD) simulations explaining their plausible origin at an atomistic level.

Revisiting the concept of peptide bond planarity in an iron-sulfur protein by neutron structure analysis

The planarity of the peptide bond is important for the stability and structure formation of proteins. However, substantial distortion of peptide bonds has been reported in several high-resolution structures and computational analyses. To investigate the peptide bond planarity, including hydrogen atoms, we report a 1.2-Å resolution neutron structure of the oxidized form of high-potential iron-sulfur protein. This high-resolution neutron structure shows that the nucleus positions of the amide protons deviate from the peptide plane and shift toward the acceptors. The planarity of the H─N─C═O plane depends strongly on the pyramidalization of the nitrogen atom. Moreover, the orientation of the amide proton of Cys⁷⁵ is different in the reduced and oxidized states, possibly because of the electron storage capacity of the iron-sulfur cluster.

On the roles of methionine and the importance of its microenvironments in redox metalloproteins

The amino acid residue methionine (Met) is commonly thought of as a ligand in redox metalloproteins, for example in cytochromes c and in blue copper proteins. However, the roles of Met can go beyond a simple ligand. The thioether functional group of Met allows it to be considered as a hydrophobic residue as well as one that is capable of weak dipolar interactions. In addition, the lone pairs on sulphur allow Met to interact with other groups, inluding the aforementioned metal ions. Because of its properties, Met can play diverse roles in metal coordination, fine tuning of redox reactions, or supporting protein structures. These roles are strongly influenced by the nature of the surrounding medium. Herein, we describe several common interactions between Met and surrounding aromatic amino acids and how they affect the physical properties of both copper and iron metalloproteins. While the importance of interactions between Met and other groups is established in biological systems, less is known about their roles in redox metalloproteins and our view is that this is an area that is ready for greater attention.

Carbon-Centered Hydrogen Bonds in Proteins

Hydrogen bonding (H-bonding) without lone pair(s) of electrons and π-electrons is a concept developed 2-3 years ago. H-bonds involving less electronegative tetrahedral carbon are beyond the classical concept of H-bonds. Herein, we present the first report on H-bonds with tetravalent carbons in proteins. A special bonding arrangement is needed to increase the negative charge density around the sp³-hybridized carbon atom. Therefore, less electronegative elements such as As and Mg, when bonded to sp³-C, enable the C-atoms as H-bond acceptors. Careful protein structure analysis aided by several quantum chemical calculations suggests that these H-bonds are weak to moderate in strength. We developed an empirical equation to estimate the C-H···C H-bond energy in proteins from the distances between the C- and H-atoms. In proteins, the binding energies range from -5.4 to -14.0 kJ/mol. The C-H···C H-bonds assist the substrate binding in proteins. We also explored the potential role of these carbon-centered H-bonds in C-H bond activation through σ-bond metathesis. To our surprise, contribution from these H-bonds is almost of similar magnitude as that from C-H···π H-bonds for C-H bond activation.

The synthesis and crystal structures of two μ-oxo iron(III) porphyrin complexes

Two [Formula: see text]-oxo iron(III) porphyrin complexes, [Formula: see text]-oxo-bis[meso-tetra-[Formula: see text],[Formula: see text],[Formula: see text],[Formula: see text]-([Formula: see text]-nicotina-midophenyl)porphinatoiron(III)] ([Fe(TPyPP)]2O) and its copper adduct ([Fe(TPyPP)(CuCl)]2O[CF3SO3][Formula: see text], which are derived from Gunter’s porphyrin model for cytochrome [Formula: see text]oxidase are reported. The two complexes were isolated and characterized by single crystal X-ray diffraction, UV-vis and IR spectroscopies. Crystal structural analysis including comparisons to [Formula: see text]-oxo analogues and Hirshfeld surface calculations revealed several noteworthy structural features [Formula: see text] the various picket orientations which are partially attributed to intra- and/or inter-molecular nonbonded interactions.

Molecular Tailoring Approach for Estimating Individual Intermolecular Interaction Energies in Benzene Clusters

There is no general method available for the estimation of individual intermolecular interaction energies in weakly bound molecular clusters, and such studies are limited only to the dimer. Recently, we proposed a molecular tailoring approach-based method for the estimation of individual O-H···O hydrogen bond energies in water clusters. In the present work, we extend the applicability of this method for estimating the individual intermolecular interaction energies in benzene clusters, which are expected to be small. The basis set superposition error (BSSE)-corrected individual intermolecular interaction energies in linear (LN) benzene clusters, LN-(Bz)_n n = 3-7, were calculated to be in the range from -1.75 to -2.33 kcal/mol with the cooperativity contribution falling between 0.05 and 0.20 kcal/mol, calculated at the MP2.5/aug-cc-pVDZ level of theory. In the case of non-linear (NLN) benzene clusters, NLN-(Bz)_n n = 3-5, the BSSE-corrected individual intermolecular interaction energies exhibit a wider range from -1.16 to -2.55 kcal/mol with cooperativity contribution in the range from 0.02 to -0.61 kcal/mol. The accuracy of these estimated values was validated by adding the sum of interaction energies to the sum of monomer energies. These estimated molecular energies of clusters were compared with their actual calculated values. The small difference (<0.3 kcal/mol) in these two values suggests that our estimated individual intermolecular interaction energies in benzene clusters are quite reliable.

PEGylation Increases the Strength of a Nearby NH-π Hydrogen Bond in the WW Domain

Here we show that an NH-π interaction between a highly conserved Asn and a nearby Trp stabilizes the WW domain of the human protein Pin1. The strength of this NH-π interaction depends on the structure of the arene, with NH-π interactions involving Trp or naphthylalanine being substantially more stabilizing than those involving Tyr or Phe. Calculations suggest arene size and polarizability are key structural determinants of NH-π interaction strength. Methylation or PEGylation of the Asn side-chain amide nitrogen each strengthens the associated NH-π interaction, though likely for different reasons. We hypothesize that methylation introduces steric clashes that destabilize conformations in which the NH-π interaction is not possible, whereas PEGylation strengthens the NH-π interaction via localized desolvation of the protein surface.

Synthesis and Optical Properties of Axially Chiral Bibenzo[b]carbazole Derivatives

Pure organic materials with the circularly polarized luminescence (CPL) property have attracted significant research interest over the past few decades. In this study, a series of axially chiral bibenzo[b]carbazole derivatives were synthesized by adopting palladium- and iridium-catalyzed direct C-H functionalization reactions as the key steps. These compounds exhibited CPL characteristics with considerably large dissymmetry factors up to 2.81 × 10^-2 in the solid state, indicating the formation of well-ordered aggregates.

Photoinduced electron transfer in nano-Saturn complexes of fullerene

The photoinduced electron transfer is studied computationally in several Saturn-shaped inclusion complexes of carbo-aromatic rings and C60 fullerene - C72⊃C60, C96⊃C60, C120⊃C60, and C168⊃C60. Analysis of their structural and electronic properties shows that the charge separation process is efficient in C120⊃C60 and C168⊃C60 where the host molecule resembles the conjugated [24]circulene unit. In contrast, the electron transfer is not feasible in the complexes of the oligophenylene-based rings C72⊃C60 and C96⊃C60.

Dual-Site Binding of Quaternary Ammonium Ions as Internal K+-Ion Channel Blockers: Nonclassical (C-H···O) H Bonding vs Dispersive (C-H···H-C) Interaction

A molecular-level study of the influence of the alkyl chain length of quaternary ammonium ions (QAs) on the blocking action and the mode of binding with the bacterial KcsA K⁺-ion channel is carried out by molecular dynamics (MD) simulations as well as quantum mechanics/molecular mechanics (QM/MM) methods. The present work unveils distinct modes of binding for different QAs, due to differences in size and hydrophobicity. The QAs bind near the channel gate as well as at the central cavity, leading to a possible dual-site blocking action. Small-sized tetraethylammonium (TEA) and tetrabutylammonium (TBA) ions enter inside the channel cavity in the open state of KcsA but bind strongly in the closed state. TEA binds to the polar hydroxyl group of threonine residues situated at the channel gate via nonclassical H-bonding interaction (C-H···O), while TBA binds to a second binding site, the central cavity, with hydrophobic benzyl and sec-butyl side chains of phenylalanine and isoleucine residues via alkyl-π and hydrophobic interactions (C-H···H-C). On the contrary, large tetrahexylammonium (THA) and tetraoctylammonium (TOA) ions bind the hydrophobic side-chain methyl and isopropyl of alanine and valine at the channel gate both in the open and closed states, thereby restricting the free movement of large QAs toward the center of the cavity. However, the binding to the hydrophobic benzyl and sec-butyl side chains of phenylalanine and isoleucine residues in the closed state is thermodynamically preferable. Also, the binding energy is found to increase with an increase in the alkyl chain length from ethyl (-16.4 kcal/mol) to octyl (-65.5 kcal/mol), due to an almost linear increase in dispersive interaction.

Xanthones for melanogenesis inhibition: Molecular docking and QSAR studies to understand their anti-tyrosinase activity

The human skin is constantly exposed to external factors that affect its integrity, UV radiation being one of the main stress factors. The repeated exposure to this radiation leads to increased production of Reactive Oxygen Species (ROS) which activate a series of processes involved in photoaging. Excessive UV exposure also exacerbates melanin production leading to a variety of pigmentation disorders. Xanthones are reported to exhibit properties that prevent deleterious effects of UV exposure and high levels of ROS in the organism, so in this work a wide library of xanthones with different patterns of substitution was synthesized and tested for their inhibitory activity against the skin enzymes tyrosinase, elastase, collagenase and hyaluronidase, many of which were evaluated for the first time. Most of the compounds were tyrosinase inhibitors, with the best one (xanthone 27) presenting an IC₅₀ of 1.9 µM, which is approximately 6 times lower than the IC₅₀ of the positive control kojic acid. Concerning the other enzymes, only one compound presented IC₅₀ lower than 150 µM in elastase inhibition (xanthone 14 = 91.8 µM) and none in collagenase and hyaluronidase inhibition. A QSAR model for tyrosinase inhibitory activity was built using six molecular descriptors, with a partial negative surface area descriptor and the relative number of oxygen atoms being positively contributing to the tyrosinase inhibitory activity. Docking using AutoDock Vina shows that all the tested compounds have more affinity to mushroom tyrosinase than kojic acid. Docking results implied that the tyrosinase inhibitory mechanisms of xanthonic derivatives are attributed to an allosteric interaction. Taken together, these data suggest that xanthones might be useful scaffolds for the development of new and promising candidates for the treatment of pigmentation-related disorders and for skin whitening cosmetic products.