Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Carousel with three slides shown at a time. Use the Previous and Next buttons to navigate three slides at a time, or the slide dot buttons at the end to jump three slides at a time. Plastic Trestle Table

Reetuparna Acharya, Shinu Chacko, … Shakti Prasad Pattanayak

Firoz A. Dain Md Opo, Mohammed M. Rahman, Abdullah M. Asiri

Deepa Mandal, Paritosh Patel, … Tithi Parija

Dharmendra K. Yadav, Surendra Kumar, … Mi-hyun Kim

Kiran Bharat Lokhande, Shuchi Nagar & K. Venkateswara Swamy

Juan Marcelo Carpio Arevalo & Juliana Carolina Amorim

Ninad V. Puranik, Pratibha Srivastava, … Jayanthi Sivaraman

Rajkumar Veligeti, Rajesh Bagepalli Madhu, … Grigory Vasilievich Zyryanov

Scientific Reports volume  12, Article number: 22446 (2022 ) Cite this article

Cladosporium spp. have been reported for their great diversity of secondary metabolites which represent as a prominent base material for verifying the biological activities. Several bioactive compounds which have antimicrobial, cytotoxic, quorum sensing inhibitory and phytotoxic activities have been isolated from Cladosporium species. Most of them are still needed to be explored for their anticancer properties. Therefore, the present study is focused on screening and identifying the bioactive compounds of Cladosporium spp. for their anticancer activity via the integrated approaches of Molecular Docking (MD), Molecular Dynamics Simulation (MDS) and Density Functional Theory (DFT) studies. A total of 123 bioactive compounds of Cladosporium spp. were explored for their binding affinity with the selected breast cancer drug target receptor such as estrogen receptor alpha (PDB:6CBZ). The Molecular Docking studies revealed that amongst the bioactive compounds screened, Altertoxin X and Cladosporol H showed a good binding affinity of − 10.5 kcal/mol and − 10.3 kcal/mol, respectively, with the estrogen receptor alpha when compared to the reference compound (17\(\upbeta\) -Estradiol: − 10.2 kcal/mol). The MDS study indicated the stable binding patterns and conformation of the estrogen receptor alpha-Altertoxin X complex in a stimulating environment. In addition, in silico absorption, distribution, metabolism, excretion and toxicity (ADMET) study suggested that Altertoxin X has a good oral bioavailability with a high LD\(_{50}\) value of 2.375 mol/kg and did not cause any hepatotoxicity and skin sensitization. In summary, the integrated approaches revealed that Altertoxin X possesses a promising anticancer activity and could serve as a new therapeutic drug for breast cancer treatment.

Presently, breast cancer is one of the most common cancers in women in the world in terms of 2.26 million new cases and 0.685 million deaths occurring worldwide in 20201,2. The global breast cancer mortality rate is alarmingly high due to delayed diagnosis or treatment availability3. Approximately 80% of breast cancers in women are from hormone-dependent estrogen receptor positive breast cancer. Since the estrogen receptor alpha is mainly responsible for the initiation and progression of breast cancer, it has emerged as the single most important target for treating the disease4,5. Therefore, there is a lot of emphasis on developing new effective therapeutic drugs which specifically bind to estrogen receptor alpha and prevent the progression of hormone-dependent breast cancer. To date, many ligands such as estradiol, genistein, daidzein, biochanin A, 6,3’,4’-trihydroxyflavone, phloretin, ellagic acid, (-)-epigallocatechin-3-gallate, ursolic acid, kaempferol, naringenin, toxaphene, chlordane, thiadiazole acrylamide (XCT790), diethyl stilbestrol and deketene curcumin derivatives have been successfully identified as antagonists of estrogen receptor alpha by the virtual screening of structural model44,6,7. Even though plants have long been utilized to treat several diseases, including cancer, the bioactive compounds produced by microbes have yet to be investigated. Although plants have long been utilized to treat several diseases, including cancer, the bioactive compounds produced by microbes have yet to be investigated.

The availability of many microbial-originated antibiotics in the market8 has gained a great interest in the discovery of anticancer drug lead compounds from the unexplored microbial sources, especially fungi. Fungi are well-known for their ubiquitous occurrence as they are one of the most abundant sources of chemodiversity9,10. The genus Cladosporium is considered as a rich source of diverse and bioactive natural compounds that are reported for their biological properties including antimicrobial, cytotoxic, quorum sensing inhibitory and phytotoxic activities belonging to various classes of secondary metabolites such as naphtalenones, alkaloids, flavonoids, lactones, benzofluoranthenones, macrolides, coumarins, perylenequinones, azaphilones, isocoumarins, sterols and others11. Taxol extracted from Cladosporium oxysporum has been extensively reported for its anticancer as well as antibacterial properties12. However, most of the compounds in Cladosporium spp. have not yet been exploited for their anticancer properties.

The use of bioinformatics tools to find drugs that could be used to combat several diseases is commonplace nowadays. Structure-based drug design enables virtual screening of compounds for possible therapeutic candidates by analyzing their binding affinities with protein receptors before in vitro and in vivo experiments13,14. The integrated computational approaches of MD, MDS and DFT calculation are used to analyze their performance during interaction with protein receptors. Therefore, the present study was focused on the exploration of Cladosporium spp. bioactive compounds for their anticancer properties through the integrated approaches of Molecular Docking, Molecular Dynamics Simulation and DFT studies as well as absorption, distribution, metabolism, excretion and toxicity (ADMET) prediction calculations. The estrogen receptor alpha, a key transcription factor in breast cancer was used as the breast cancer drug target.

A total of 123 bioactive compounds from Cladosporium spp. reported from a recently published literature11 and the reference compound (17\(\upbeta\) -estradiol) were selected for molecular docking and simulation studies against the estrogen receptor alpha (PDB: 6CBZ) protein receptor responsible for breast cancer. Since the estrogen receptor alpha is an important biological target mediating 17\(\upbeta\) -estradiol driven breast cancer development, 17\(\upbeta\) -estradiol was used as the reference compound in the study. The 3D and 2D structures of the compounds were downloaded from the online PubChem database (https://pubchem.ncbi.nlm.nih.gov/) in structure data file (SDF) format. Some of them were created using MarvinSketch (version 18.30) ChemAxon chemical drawing tool (https://chemaxon.com/products/marvin/). The 2D structures were converted into 3D coordinates and geometries which are then converted into protein data bank (PDB) format using the open-source chemical toolbox, Open Babel15. Before molecular docking, the geometries of the ligand PDB files were optimized in the PRODRG server (https://davapc1.bioch.dundee.ac.uk/cgi-bin/prodrg/submit.html). These optimized structures were used as the ligand molecules for the molecular docking studies.

The three-dimensional structure of the estrogen receptor alpha ligand-binding domain Y537S mutant in complex with estradiol and GRIP peptide (PDB: 6CBZ) was retrieved from the PDB database (https://www.rcsb.org/) in PDB format. The protein receptor was prepared for Molecular Docking by removing ligand and water molecules attached to it using Discovery Studio Visualizer (version 20.1.0.19295) software to avoid interfering with the docking study. The energy minimization studies were performed with empirical force fields using Swiss-Pdb Viewer software (version 4.10) to generate the protein structure’s lower energy conformations, which infers a greater stable conformation. The process optimizes conformational errors in the structure’s geometry during the protein structure modeling. In addition, the steepest descent algorithm was employed using the GROMOS 96 force field for the geometry optimization process16. The PDBsum database (https://www.ebi.ac.uk/thornton-srv/databases/pdbsum/Generate.html) was used to validate the protein structure with the Ramachandran plot17, which suggested that the majority (96.8%) of the amino acids residues were found within the most favored regions of the protein used.

The AutoDock Vina implicated in PyRx Virtual Screening software (version 0.8) was used to execute the molecular docking study by considering the estrogen receptor alpha (PDB: 6CBZ) protein and bioactive compounds as macromolecule and ligand molecules, respectively. AutoDock Vina was used for fixing the cubical grid box size at 60 \(\times\) 60 \(\times\) 60 with 0.375 Å around the active sites of the protein. The 100 independent docking runs were performed for each compound. The most favorable binding pose was selected based on the lowest free binding energy (kcal/mol). The molecular interaction between the amino acid residues of protein receptors and ligand molecules was studied by visualizing the docking result using Discovery Studio Visualizer software. The accuracy of the docking protocol was validated through re-docking (self-docking) of the compounds with the protein used during the study.

The molecular dynamics simulations for the protein alone and protein-ligand complex were carried out using the GROMACS software (version 5.1.4) with the AMBER force field, modified Berendsen thermostat and LINear Constraint Solver (LINCS) constraint algorithm. The docked complex was cleaned and optimized first, followed by the orientation of hydrogen bond network systems. The parameter files were generated for the ligand using SwissParam web server. The TIP3P solvation model was used to solvate the cubic simulation system with water and NaCl counter-ions were added to neutralize the charge of the system. The steepest (gradient) descent algorithm was used to perform the initial energy minimization of the system with 5000 steps. The system was then equilibrated using NVT (or canonical) ensemble followed by NPT (isothermal-isobaric) ensemble for 100 ps. The protonation states of key amino acid residues such as Histidine (His) in 6CBZ protein as well as selected ligand molecules (Altertoxin X, Cladosporol H and 17\(\upbeta\) -estradiol) were analyzed at neutral pH using the PROPKA web server, which provided the pKa value around His residues18. Finally, the molecular dynamics simulations were performed for 100 ns with an integration time of 0.002 ps under the physiological parameters (temperature: 300 K and pressure: 1 bar) of simulation systems. The simulation trajectories were exploited to calculate the root mean square deviation (RMSD), root mean square fluctuation (RMSF), the radius of gyration (Rg), solvent accessible surface area (SASA) and the number of hydrogen bonds.

The outcome of the MD simulation run for the target protein complexed with Altertoxin X, Cladosporol H and 17\(\upbeta\) -Estradiol was subjected to binding free energy calculations using the Molecular MechanicsPoisson-Boltzmann Surface Area (MM-PBSA) technique. It was another application of molecular dynamics simulations and thermodynamics for determining the extent of ligand binding with protein. The gmx\(\_\) MMPBSA program with MMPBSA Stat.py script, which utilizes the GROMACS 2018.1 trajectories as input, was used to determine the binding free energy for each ligand-protein combination. The gmx\(\_\) MMPBSA program used three components to calculate the binding free energy: molecular mechanical energy, polar and apolar solvation energies, and molecular mechanical energy. The calculations were based on MD trajectories of last 100 ns, which compute \(\Delta\) G with dt 1000 frames. It was evaluated using molecular mechanical energy, polar and apolar solvation energies. The equations (1) and (2) used to calculate the free binding energy are given below.

where G\(_{Binding}\) : Binding free energy, G\(_{Complex}\) : Total free energy of the protein-ligand complex, G\(_{Protein}\) and G\(_{Ligand}\) : Total free energies of the isolated protein and ligand in solvent, respectively, \(\Delta\) G: Standard free energy, \(\Delta\) E\(_{MM}\) : Average molecular mechanics potential energy in vacuum, \(\Delta\) G\(_{Solvation}\) : Solvation energy, \(\Delta\) E: Total energy of bonded as well as non-bonded interactions, \(\Delta\) S: Change in entropy of the system upon ligand binding; T. Temperature in Kelvin.

The molecular energy, electronic density, and frontier orbital energies, chemical reactivity descriptors and of the studied Altertoxin X, Cladosporol H and 17\(\upbeta\) -Estradiol molecular systems were determined using the Kohn-Sham (KS) approach19,20 while making use of the Conceptual DFT (CDFT) methodology21,22,23,24. Many different conformers of the studied compounds were determined using MarvinView 17.15 from ChemAxon [http://www.chemaxon.com] through the consideration of the MMFF94 force field to perform Molecular Mechanics calculations25,26,27,28,29. This was followed by a geometry optimization and frequency calculation by means of the Density Functional Tight Binding (DFTBA) methodology30 and a later geometry reoptimization, frequency analysis and calculation of the electronic properties and the chemical reactivity descriptors by means of the MN12SX/Def2TZVP/H2O model chemistry31,32,33 on their optimized molecular structures. The charge of the molecule was taken as equal to zero while the radical anion and cation have been considered in the doublet spin state. This determination was performed with the aid of the Gaussian 16 software30 and the SMD solvation model34 and owing to the fact that the mentioned model chemistry has been previously proved as verifying the ’Koopmans in DFT’ (KID) procedure35,36,37,38, This last step was also required for the verification of the absence of imaginary frequencies as a check for the stability of the optimized structure as being a minimum in the energy landscape.

The physicochemical and ADMET properties of the selected potential compounds were predicted by utilizing the pkCSM platferm39. These parameters are related to the absorption, distribution, metabolism, excretion and toxicity of drugs.

Absorption parameters: water solubility in a buffer system, Caco2 cell permeability, intestinal (human) absorption, P-glycoprotein inhibition and skin permeability.

Distribution parameters: Lipinski’s rule, blood-brain barrier (BBB) and central nervous system (CNS) permeability.

Metabolism parameters: Cytochrome P450 (CYP)1A2, CYP2C19, CYP2C9, CYP2D6 and CYP3A4 inhibition, CYP2D6 and CYP3A4 substrate.

Excretion parameters: total renal clearance and renal OCT2 substrate.

Toxicity parameters: AMES test, oral rat acute and chronic toxicity, skin sensitization, hepatotoxicity, Tetrahymena pyriformis toxicity and Minnow toxicity.

Extensive studies have demonstrated the anticancer properties of the compounds found in many plants, but microbial compounds’ use for cancer treatment are still yet to be explored. In the present study, the binding potential of bioactive compounds from Cladosporium spp. was explored with the estrogen receptor alpha through the integrated in silico approaches to identify inhibitory effect. The estrogen receptor alpha is a critical transcription factor for human breast cancer development, which accounts for approximately 70% of breast cancer40. It significantly contributes to the progression and inhibition of breast cancer in women. Due to its over-expression, the estrogen receptor alpha could help us to find the new effective ways to treat the breast cancer. Therefore, the estrogen receptor alpha was used as a promising protein target to determine the drug candidate for breast cancer therapy.

The outcome of the molecular docking analysis suggested that six bioactive compounds (viz., Cladosporol H, Cladosporol J, Altertoxin VIII, Altertoxin IX, Altertoxin X and Altertoxin XI) were found to obtain the higher binding affinities ranging from − 9.9 to − 10.5 kcal/mol (Table 1). The interaction of compounds with the protein target (PDB: 6CBZ) resulted in varied binding potential. Among the evaluated compounds, the best binding energy results were noticed with Altertoxin X (− 10.5 kcal/mol) followed by Cladosporol H (− 10.3 kcal/mol), while − 10.2 kcal/mol was noticed with the reference compound 17\(\upbeta\) -Estradiol. These results indicated that the binding affinity was increased by modifying (replacing) the hydroxyl (-OH) group with a long carbon chain/aromatic ring molecule41. The best binding pose clearly suggested that the potential bioactive compounds as drugs candidates bind within the estrogen receptor alpha (Fig. 1). The 2D schematics indicated that amino acids played a significant role in the pattern of interactions between the protein and the ligands, significantly contributing to the total energy of the interaction.

Besides, the hydrogen bonding significantly enhanced the binding energy and binding mode, which is important for influencing the ligand binding specificity with the receptor, drug design in chemical and biological processes, molecular recognition, and biological activity42. The best docked complex of Altertoxin X with the estrogen receptor alpha protein with binding energy − 10.5 kcal/mol was firmly bound through conventional hydrogen bonds with the residues of LEU387, ARG394 and HIS524 (Fig. 2). It has been reported that the molecular docking was utilized for virtual screening of the compounds produced in fungal endophytes (Chaetomium sp.) against the human estrogen receptor alpha (PDB: 1G50) for searching anti-breast cancer agents43. The observed binding energies for virtual active compounds ranged from − 9.2 to -4.9 kcal/mol with the corresponding active amino acid residues (such as ARG394, GLU323, GLU353, LEU345, LEU346, LYS449, PRO324, PRO325 and TRP393).

The findings of Ervina et al44 have noticed that Quercetin 3-O-(2”,6”-digalloyl)-\(\upbeta\) -D-galactopyranoside isolated from Melia azearach leaves bioactive ethyl acetate fraction offered binding energy of − 9.9 kcal/mol during molecular docking with estrogen receptor (PDB: 3ERT) which was lower that its agonist 17\(\upbeta\) -Estradiol (− 9.4 kcal/mol). The better potency of Quercetin 3-O-(2”,6”-digalloyl)-\(\upbeta\) -D-galactopyranoside to bind the estrogen receptor alpha indicatred its higher binding affinity to the active sites of 3ERT. In addition, Muhammad et al45 have reported that among the 10 selective bioactive Curcumin derivatives, Salicydenecurcumin, 4-benzylidene Curcumin, and difluorinated Curcumin were found to show the lowest binding energy of -8.6 kcal/mol, -8.8 kcal/mol, and − 9.0 kcal/mol, respectively with the same protein 3ERT. Therefore, the MD simulation studies were conducted with Altertoxin X and Cladosporol H based on their better binding potential along with the reference compound.

Docking model of the potential bioactive compounds of Cladosporium spp. in the active binding site of the estrogen receptor alpha (PDB: 6CBZ) protein.

Non-bonding interaction of the potential bioactive compounds extracted from Cladosporium spp. against estrogen receptor alpha (PDB: 6CBZ) protein.

Among the different amino acids, only His amino acid residues are known to play a physiological role at neutral pH due to different protonated structures based on the pKa value around them46. At neutral pH, the His residues in proteins can attain three different protonation states: (1) Hid protonation has a hydrogen atom at the \(\updelta\) -site of its imidazole ring; (2) Hie protonation has a hydrogen atom at the \(\upepsilon\) -site of the imidazole ring, and (3) positively charged Hip protonation has hydrogen atoms at both the \(\updelta\) - and \(\upepsilon\) -sites. Therefore, the His residues with pKa value > 6 have the Hip protonation, while those residues with pKa value <6 have the Hid or Hie protonation18 . During the determination of the protonation state of the protein selected (6CBZ), it was noted that the protein consists of 24 His residues with varying pKa values (Supplementary Table 1) wherein it was noted that only 12 His residues had pKa values < 6 and designated to have Hie protonation. The His residues with only Hie protonation were considered because practically it is impossible to consider all the combinations of Hid/Hie protonation.

In addition, the protein binding potential of the molecules can be encouraged or discouraged by changing hydrogen-bond, electrostatic, and Van der Waals interactions. During protonation, the proton (or hydrogen cation, H\(^+\) ) is added to an atom, molecule, or ion to form a conjugate acid, thereby enhancing its binding affinity with the protein receptor. The study results showed that Altertoxin X, Cladosporol H and 17\(\upbeta\) -Estradiol at the neutral pH were found to be in protonated state (pKa > 7), thereby improving their binding ability.

The best docked bioactive compounds, Altertoxin X and Cladosporol H, along with a reference compound (17\(\upbeta\) -Estradiol) in complex with the receptor (estrogen receptor alpha; PDB: 6CBZ), were further subjected to the molecular dynamics simulation to study the interaction and stability of the protein-ligand complex in an aqueous system for a simulation time of 100 ns. The simulation study provided the geometric properties of the protein-ligand complex that included the analysis of RMSD (Root Mean Square Deviation), RMSF (Root Mean Square Fluctuation), Rg (Radius of Gyration), SASA (Solvent Accessible Surface Area) and a number of hydrogen bonds maintained throughout the simulation time and variation in receptor alone and the three complexes (Figure 3)47. These calculations were analyzed to understand the simulation systems’ structural variations and stability. The RMSD of atomic positions is used to calculate the average distance between the atoms of superimposed protein and ligand structures over time48. The RMSD plot of the estrogen receptor alpha protein alone showed that the estrogen receptor alpha protein alone reached the equilibrium approximately at 80–100 ns time and the remaining showed the stable trajectory simulation with minimal deviation in 0.35-0.47 nm RMSD range, whereas the RMSD plot of estrogen receptor alpha-Altertoxin X complex reached the equilibrium at 40–100 ns and 0.15-0.21 nm. The structural flexibility of protein is reserved when the estrogen receptor alpha protein is in free form in its complex with Altertoxin X compared to Cladosporol H and 17\(\upbeta\) -Estradiol. The ligand molecule, Altertoxin X bound to the estrogen receptor alpha protein, reached equilibrium after initial fluctuations, showing that this molecule has fewer deviations than the other molecules.

The RMSF concentrates on the protein structural regions that differ the most/least from the mean. Further, by calculating root mean square distances with respect to the central axis of rotation. The RMSF plots were predicted with fluctuations only at the terminal ends and loop regions, indicating the stable interactions between the complexes. The protein complexed with Altertoxin X exhibited fewer fluctuations, indicating it to be a stable molecule. The Rg is the root mean square distance between each protein atom and its center of mass in a system49 . The Rg plot analyzed the capability, shape and folding during every time step of the trajectory throughout the simulation. The estrogen receptor alpha protein complexed with Altertoxin X exhibited a similar pattern of Rg values with a deviation of 1.3–1.7 nm. The solvent-accessible surface area of the simulation complex was explored to understand better the complex’s changes in surface area, where a higher SASA indicates the extension of surface volumes and a lower SASA indicates the truncated nature of the complex50. The SASA calculates the surface area of the hydrophobic core generated by protein-ligand interactions. Consistent SASA values were observed in the estrogen receptor alpha protein-Altertoxin X complex compared to other complexes.

In addition, the hydrogen bond must be assessed in the biological system to determine the bonding and structural changes in the complex. The stability of the complex is defined by the hydrogen bond of the simulation system, where the entire complex exhibits a stable trend51. The hydrogen bonds during molecular docking are evaluated across the simulation duration. All intermolecular hydrogen bonds between estrogen receptor alpha (protein) and Altertoxin X (ligand) were solely examined during the analysis and shown accordingly. The plot showed that the number of hydrogen bonds formed during simulation runs was consistent with the molecular docking study, with only a few bonds being broken and repaired simultaneously compared to Cladosporol H and 17\(\upbeta\) -Estradiol. The MD simulation verified the stability of docked complexes Salicylidenecurcumin and Curcumin difluorinated within the interaction cavity of the estrogen receptor alpha 3ERT in humans and confirmed the steadiness of the complexes over the simulated trajectories at 120 ns time scale45.

The molecular dynamics simulation trajectories from 100 ns simulation time represents the (A) RMSD plot, (B) RMSF plot, (C) RG plot, (D) SASA plot, and (E) Hydrogen bond plot (Black: 6CBZ alone, Red: 6CBZ-Altertoxin X complex, Blue: 6CBZ-Cladosporol H complex, and Green: 6CBZ-17\(\upbeta\) -Estradiol complex.

Various energy metrics such as Van der Waals, electrostatic, polar solvation, SASA, and binding energies were utilized to measure the extent of ligand-target protein binding interactions during MD simulations. In this study, the protein-ligand complex was majorly formed using the Van der Waals energy followed by the binding energy, SASA energy, and electrostatic energy. Polar solvation energy was predicted with no contribution to the protein-ligand complex formation, as the values appeared positive. The 6CBZ complexed with Altertoxin X, Cladosporol H, and 17\(\upbeta\) -Estradiol (reference compound) were considered for binding energy calculation studies. In addition, the protein-ligand complex standard deviations were calculated. A lower standard deviation means the data values are closer to the mean (or expected value), whereas a high standard deviation means the data values are spread out over a wider range. However, there was no high standard deviation in the Altertoxin X- receptor complex compared to Cladosporol H and 17\(\upbeta\) -Estradiol, which has high standard deviation values, which indicates that Altertoxin X binds to the protein with high binding affinity and stable interaction. The binding free energy calculations of the protein-ligand complex have been represented in Table 2:

Oral bioavailability is one of the most important properties in drug design. A higher score reduces the amount of an administered drug necessary to achieve the desired pharmacological effect, thereby lessening the risk of side effects and toxicity52. The 90% of orally bioactive drugs that have progressed to the second stage of the clinical trial are associated with the four simple physicochemical properties ranges the molecular weight \(\le\) 500, logP \(\le\) 5, number of hydrogen bond acceptors \(\le\) 10, and number of hydrogen bond donors \(\le\) 553. The excellent physicochemical properties of Altertoxin X (Table 3) are responsible for its interaction with the amino acid residues of estrogen receptor alpha protein by forming hydrogen bonds. In addition, the orally bioactive drugs that are passively transported via the transcellular route should not exceed their polar surface area of about 120 Å2 and should be tailored to less than 60–70 Å2 for good brain penetration of drugs54. The polar surface area \(\le\) 140 Å2 predicted to have a high probability of good oral bioavailability. Altertoxin X was found to have good oral bioavailability. The prediction of ADMET properties is critical in drug discovery and development because these properties are responsible for around 60% of all drug failures in clinical trials. The ADMET is used early in the drug development process to eliminate compounds with poor ADMET properties from the pipeline, reducing the research and development cost. In addition, the ADMET analysis of the selected compounds corroborates with the MD simulation studies’ findings, thereby supporting the utilization of the same for further in vitro and in vivo studies (Table 4 and Supplementary Table 2).

The absorption of drugs from an orally administered solution depends on the factors including water solubility, Caco-2 cell membrane permeability, human intestinal absorption, skin permeability threshold, and substrate or inhibitor of P-glycoprotein. High water solubility is one of the useful factors for delivering a sufficient quantity of active ingredients in the small volume of such pharmaceutical dosage55. According to log S scale, the water solubility of the drug molecules is considered as poorly soluble if the value is < − 10 mol/L; moderately soluble if < − 6 mol/L; soluble if < − 4 mol/L; very soluble if < − 2 mol/L, and highly soluble if < 0 mol/L. From the results, it was observed that Altertoxin X tested was moderately soluble in water. The apparent permeability coefficient (Papp) value of drug molecules measures the rate at which they can cross the area of Caco-2 cell monolayer. The high Caco-2 permeability is translated into the predicted log Papp value of \(\ge\) 0.90 cm/s in the pkCSM predictive model. The predicted Caco-2 permeability value (log Papp 1.612 cm/s) of Altertoxin X indicated a low Caco-2 permeability. The human intestinal absorption value greater than 90% for Altertoxin X indicated its excellent absorption. The skin permeability (Kp) measures the rate at which the drug molecules penetrate the stratum corneum. The Kp value is widely used to quantitatively describe the transport of drug molecules into the outermost layer of epidermal skin (or stratum corneum) and indicates the significance of skin absorption39. The drug molecules will easily penetrate the skin if their log Kp value is > − 2.5 cm/h. Altertoxin X also has the predicted log Kp value of − 2.738 cm/h (\(\le\) − 2.5 cm/h) indicating its poor skin permeability. P-glycoprotein is a member of the ATP-binding cassette (ABC) transporter family that actively transports various compounds out of the cells. Here, Altertoxin X was predicted to be transported across the cell membrane through ABC transporter.

The steady-state volume of distribution (VDss) is one of the pharmacokinetic parameters that representing the propensity of drugs to either remain in blood plasma or redistribute to another tissues56. The high VDss value (> 0.45 L/kg) indicates the propensity of drugs to leave the plasma and enter the other tissues of the body, while the low VDss value (< − 0.15 L/kg) indicates the propensity of drugs to remain in the plasma. A higher drug dose is required to achieve in a given plasma concentration if it has a high VDss value due to more distribution of drug to other tissues. Conversely, a lower drug dose is required to achieve in a given plasma concentration if it has a low VDss value due to less drug distribution to other tissues. The compound Altertoxin X used in the present study has a relatively low VDss value (0.025 L/kg). The prediction of blood-brain barrier (BBB) permeability is an important factor for regulating the transportation of drugs from and to the central nervous system (CNS)57. The log BB is the most common numerical value describing the BBB permeability. In the qualitative model, the drug is considered BBB permeable if the log BB value \(\ge\) 0.3 and non-permeable if the log BB value \(\le\) − 0.3. Altertoxin X could be easily passed the BBB because it has the log BB value of 0.331. In addition, the drugs that have log PS value \(\ge\) − 2 are considered to penetrate the CNS, while those with log PS value \(\le\) − 3 have difficulty in penetrating the CNS. It can be noticed that Altertoxin X with the log PS value − 1.839 was believed to penetrate the CNS.

The inhibition of major human cytochrome P450 (CYPs) monooxygenase enzymes is involved in the metabolism of drugs in order to prevent undesirable adverse effects. It was observed that Altertoxin X was found to inhibit all the enzymes except CYP2D6 and CYP3A4; thus it could not be metabolized by the CYPs monooxygenase enzymes in the body. The clearance of drugs quantitatively describes the volume of plasma from which they would be removed per unit. The pkCSM pharmacokinetics model predicts the given compounds’ total clearance (log mL/min/kg). The higher total clearance value of the compounds indicates their faster excretion processes. The predicted excretion rate of Altertoxin X was 0.003 mL/min/kg. The positive AMES toxicity test suggests that Altertoxin X was mutagenic. Altertoxin X has high LD\(_{50}\) value (2.375 mol/kg), which indicates that it was lethal to the human body only at extremely high dose. Moreover, it did not cause any hepatotoxicity and skin sensitization.

Conceptual Density Functional Theory (CDFT) is the branch of DFT that deals with the chemical reactivity of atoms and molecules and their interacting behavior. As such, it offers a large number of tools for the study the interactions between chemical systems. Thus, it has great importance in the process of drug design and discovery, mainly through virtual screening. It is expressed in the form of descriptors that can help to get a glimpse of how drugs interact with the receptors in the body. Indeed, this is not an easy task because those descriptors cannot be always related to the physiological chemistry of the human body. However, the estimation of the CDFT descriptors can help to get a qualitative idea of how these interactions proceed.

The structures of the selected ligands, Altertoxin X and Cladosporol H, together with the reference compound, 17\(\upbeta\) -Estradiol, have been optimized by following the methodology presented in the Materials and Methods section. We have resorted to three different functionals: B3LYP58,59,60, PBE061 and MN12SX31, being in all cases, Def2TZVP32,33 the chosen basis set, and \({\hbox {H}_{2}\hbox {O}}\) the solvent simulated with the SMD model34. The required frequencies calculations for the verification of the absence of imaginary frequencies as a check for the stability of the optimized structures as being a minimum in the energy landscape was performed for all the systems.

The quality of the chosen density functionals may be realized by comparing their results with those from high-level computations or from experimental values. Nevertheless, this comparison is not always computationally practicable because of the large size of the molecules or the lack of experimental results for the chemical methods being explored. Our research group has developed a methodology known as KID35,36,37,38, in order to evaluate a particular density functional with regard to its internal coherence. Within the Generalized Kohn-Sham (GKS) version of DFT, some relationships exist between the KID methodology and the Ionization Energy Theorem, which is a corollary of Janak theorem62,63. This is done by connecting \(\epsilon _H\) to -I and \(\epsilon _L\) to -A, through \(J_{I} = {\epsilon _H + E_{gs}(N-1) - E_{gs}(N)}\) , \(J_{A} = \epsilon _L + E_{gs}(N) - E_{gs}(N+1)\) and \(J_{HL} = \sqrt{{J_{I}}^2 + {J_{A}}^2}\) , being \(\epsilon _H\) and \(\epsilon _L\) the frontier orbital energies of the studied molecules. Another KID descriptor \(\Delta\) SL related to the difference in energies between the SOMO and the LUMO of the neutral system has been devised to aid in the verification of the accuracy of the methodology. A general Global KID Descriptor has been defined as \({\mathrm{GKD}} = \sqrt{{J_{I}}^2 + {J_{A}}^2 + {J_{HL}}^2 + {\Delta {SL}}^2}\) , whose value must be zero for the exact density functional meaning that it verifies the Ionization Energy theorem. The results are presented in Table 5:

It can be concluded from an evaluation of the results in Table 5 that the MN12SX density functional is the best for verifying the Ionization Energy theorem while the B3LYP and PBE0 density functionals will render values of the calculated properties with large errors. Thus, the calculated results for the Global Reactivity Descriptors derived from CDFT (including the Nucleophilicity N)21,22,23,24 for the Altertoxin X, Cladosporol H, and 17\(\upbeta\) -Estradiol acquired utilizing the in-house CDFT tool software based on the MN12SX/Def2TZVP/H2O geometry optimizations and frequency calculations are displayed in Table 6:

The electronegativity \(\chi\) and global hardness \(\eta\) are absolute chemical reactivity parameters for which no experimental analogue exists. Indeed, the observed vertical ionization energy (I) and vertical electron affinity (A) can be used to approximate them, but these values are unknown for the molecular system under investigation. However, in comparing Altertoxin X and Cladosporol H and the results from Table 1, it can be said that they correlate with their global hardness values. The electrophilicity index \(\omega\) is a balance between an electrophile’s proclivity for acquiring more electron density and its reluctance to exchanging electron density with its surroundings64. An electrophilicity scale was established by Domingo et al65,66,67, with larger than 1.5 eV for strong electrophiles, between 0.8 and 1.5 eV for medium electrophiles, and smaller than 0.8 eV for the marginal ones65,66,67, On the basis of Table 5, Altertoxin X and Cladosporol H may be classified as strong electrophiles while 17\(\upbeta\) -Estradiol is a medium electrophile. In the same way, Domingo and colleagues68,69 presented a Nucleophilicity index N based on the HOMO energy calculated using the KS technique with an arbitrary shift of the origin, using the molecule of tetracyanoethylene (TCE) as a reference. They were able to classify organic molecules as strong nucleophiles with N > 3.0 eV, moderate nucleophiles with 2.0 < N < 3.0 eV, and marginal nucleophiles with N < 2.0 eV. By re-examining Table 6, it is clear that Altrotoxin X may be classified as a strong nucleophile while Cladosporol H will act as a moderate nucleophile. Although the difference between the binding energies for Altrotoxin X and Cladosporol H is small, perhaps the better behavior of the first could be attributed to its observed nucleophilicity. The same conclusions may obtained by relating the binding energies to the inverse of the \(\omega ^{-}\) , \(\omega ^{+}\) and \(\Delta \omega \pm\) . These results could also be considered as complement of a recent study70, where Altertoxin X was identified as new perylenequinone found in the marine derived fungus Cladosporium spp. possessing quorum sensing inhibitory activity against Chromobacterium violaceum with a minimum inhibitory concentration (MIC) of 20 \(\upmu\) g/well, and whose structure was elucidated through quantum ECD calculations.

Since Alexander Fleming discovered Penicillin, the microbes have also been considered a rich source of most bioactive compounds. In the present study, the integrated approaches of Molecular Docking and MD simulation studies, as well as DFT calculations, were utilized to explore the structural insights into possible binding modes of drug-like bioactive compounds of Cladosporium spp. against the estrogen receptor alpha. Among the compounds screened, Altertoxin X and Cladosporol H were found to have better binding affinity than the reference standard (17\(\upbeta\) -Estradiol). The MD simulation study for up to 100 ns confirmed that the complexes of estrogen receptor alpha with Altertoxin X and Cladosporol H were highly stable in the biological system. The ADMET results revealed that Altertoxin X has a good oral bioavailability with a high LD\(_{50}\) value (2.375 mol/kg) and did not cause hepatotoxicity or skin sensitization. From the present study, it can be concluded that Altertoxin X and Cladosporol H, which possess the most significant inhibitory potential against the estrogen receptor alpha than the reference compound (17\(\upbeta\) -Estradiol), were reported as the antagonists of estrogen receptor alpha and as the promising therapeutic drug candidates for the treatment of breast cancer. Therefore, they can be considered candidates for pre-clinical and clinical trials as the future of clinical therapy. In a future study, the in vitro and in vivo studies will also be conducted to further validate the potential of Altertoxin X and Cladosporol H as clinical drugs for breast cancer treatment before the clinical trial.

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Sung, H. et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71, 209–249. https://doi.org/10.3322/caac.21660 (2021).

Ferlay, J. et al. Cancer statistics for the year 2020: an overview. Int. J. Cancer 149, 778–789. https://doi.org/10.1002/ijc.33588 (2021).

Pucci, C., Martinelli, C. & Ciofani, G. Innovative approaches for cancer treatment: Current perspectives and new challenges. Ecancermedicalscience 13, 961. https://doi.org/10.3332/ecancer.2019.961 (2019).

Ariazi, E., Ariazi, J., Cordera, F. & Jordan, V. Estrogen receptors as therapeutic targets in breast cancer. Curr. Top. Med. Chem. 6, 181–202. https://doi.org/10.2174/156802606776173483 (2006).

Sharma, D., Kumar, S. & Narasimhan, B. Estrogen alpha receptor antagonists for the treatment of breast cancer: A review. Chem. Central J.https://doi.org/10.1186/s13065-018-0472-8 (2018).

Pang, X. et al. Identification of estrogen receptor \(\upalpha\) antagonists from natural products via in vitro and in silico approaches. Oxid. Med. Cell. Longev. 1–11, 2018. https://doi.org/10.1155/2018/6040149 (2018).

Shah, V., Bhaliya, J. & Patel, G. M. In silico docking and ADME study of Deketene Curcumin Derivatives (DKC) as an aromatase inhibitor or antagonist to the Estrogen-alpha positive receptor (Er\(\upalpha\) +): potent application of breast cancer. Struct. Chem. 33, 571–600. https://doi.org/10.1007/s11224-021-01871-2 (2022).

Peláez, F. The historical delivery of antibiotics from microbial natural products—Can history repeat?. Biochem. Pharmacol. 71, 981–990. https://doi.org/10.1016/j.bcp.2005.10.010 (2006).

Greco, C., Keller, N. P. & Rokas, A. Unearthing fungal chemodiversity and prospects for drug discovery. Curr. Opin. Microbiol. 51, 22–29. https://doi.org/10.1016/j.mib.2019.03.002 (2019).

Staszczak, M. Fungal secondary metabolites as inhibitors of the ubiquitin-proteasome system. Int. J. Mol. Sci. 22, 13309. https://doi.org/10.3390/ijms222413309 (2021).

Salvatore, M. M., Andolfi, A. & Nicoletti, R. The genus cladosporium: A rich source of diverse and bioactive natural compounds. Molecules 26, 3959. https://doi.org/10.3390/molecules26133959 (2021).

Raj, K. G., Manikandan, R., Arulvasu, C. & Pandi, M. Anti-proliferative effect of fungal taxol extracted from cladosporium oxysporum against human pathogenic bacteria and human colon cancer cell line HCT 15. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 138, 667–674. https://doi.org/10.1016/j.saa.2014.11.036 (2015).

Article  ADS  CAS  Google Scholar 

Lionta, E., Spyrou, G., Vassilatis, D. & Cournia, Z. Structure-based virtual screening for drug discovery: Principles, applications and recent advances. Curr. Top. Med. Chem. 14, 1923–1938. https://doi.org/10.2174/1568026614666140929124445 (2014).

Murali, M. et al. Repositioning therapeutics for SARS-CoV-2: Virtual screening of plant-based anti-HIV compounds as possible inhibitors against COVID-19 viral RdRp. Curr. Pharm. Des. 28, 969–980. https://doi.org/10.2174/1381612828666220428120939 (2022).

O’Boyle, N. M. et al. Open babel: An open chemical toolbox. J. Cheminform. 3, 33. https://doi.org/10.1186/1758-2946-3-33 (2011).

Schwede, T. SWISS-MODEL: An automated protein homology-modeling server. Nucleic Acids Res. 31, 3381–3385. https://doi.org/10.1093/nar/gkg520 (2003).

Laskowski , RA , Jabłońska , J. , Pravda , L. , Vařekova , RS & Thornton , JM PDBsum: Structural summaries of PDB entries .Protein Sci.27 , 129–1https://doi.org/10.1002/pro.3289 (2017).

Kasumi, T., Araki, K., Mizushima, T., Kobayashi, H. & Kurita, N. A proposal of potent inhibitor for cancer metastasis blocking the pocket of urokinase receptor: ab initio molecular simulations. J. Phys. Conf. Ser. 433, 012034. https://doi.org/10.1088/1742-6596/433/1/012034 (2013).

Lewars, E. Computational Chemistry—Introduction to the Theory and Applications of Molecular and Quantum Mechanics (Kluwer Academic Publishers, Dordrecht, 2003).

Cramer, C. Essentials of Computational Chemistry—Theories and Models 2nd edn. (Wiley, 2004).

Parr, R. & Yang, W. Density-Functional Theory of Atoms and Molecules (Oxford University Press, 1989).

Geerlings, P., De Proft, F. & Langenaeker, W. Conceptual density functional theory. Chem. Rev. 103, 1793–1873. https://doi.org/10.1021/cr990029p (2003).

Geerlings, P. et al. Conceptual density functional theory: Status, prospects, issues. Theor. Chem. Accounts 139, 36. https://doi.org/10.1007/s00214-020-2546-7 (2020).

Chattaraj, P. K. (ed.) Chemical Reactivity Theory—A Density Functional View (CRC Press, 2009).

Halgren, T. A. Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J. Comput. Chem. 17, 490–519. https://doi.org/10.1002/(SICI)1096-987X(199604)17:5/6<490::AID-JCC1>3.0.CO;2-P (1996).

Halgren, T. A. Merck molecular force field. II. MMFF94 van der Waals and electrostatic parameters for intermolecular interactions. J. Comput. Chem. 17, 520–552 (1996). https://doi.org/10.1002/(SICI)1096-987X(199604)17:5/6\(<\) 520::AID-JCC2\(>\) 3.0.CO;2-W

Halgren, T. A. MMFF VI. MMFF94s option for energy minimization studies. J. Comput. Chem. 20, 720–729. https://doi.org/10.1002/(SICI)1096-987X(199905)20:7<720::AID-JCC7>3.0.CO;2-X (1999).

Halgren, T. A. & Nachbar, R. B. Merck molecular force field. IV. Conformational energies and geometries for MMFF94. J. Comput. Chem. 17, 587–615. https://doi.org/10.1002/(SICI)1096-987X(199604)17:5/6<587::AID-JCC4>3.0.CO;2-Q (1996).

Halgren, T. A., Merck molecular force field. V. Extension of MMFF94 using experimental data, additional computational data, and empirical rules. J. Comput. Chem., 17, 616–641. https://doi.org/10.1002/(SICI)1096-987X(199604)17:5/6%3C616::AID-JCC5%3E3.0.CO;2-X (1996).

Frisch, M. J. et al. Gaussian 16 Revision C.01. Gaussian Inc. (2016).

Peverati, R. & Truhlar, D. G. Screened-exchange density functionals with broad accuracy for chemistry and solid-state physics. Phys. Chem. Chem. Phys. 14, 16187–16191. https://doi.org/10.1039/C2CP42576A (2012).

Weigend, F. & Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 7, 3297–3305. https://doi.org/10.1039/B508541A (2005).

Weigend, F. Accurate Coulomb-fitting basis sets for H to R. Phys. Chem. Chem. Phys. 8, 1057–1065. https://doi.org/10.1039/B515623H (2006).

Marenich, A., Cramer, C. & Truhlar, D. Universal solvation model based on solute electron density and a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 113, 6378–6396. https://doi.org/10.1021/jp810292n (2009).

Flores-Holguín, N., Frau, J. & Glossman-Mitnik, D. A fast and simple evaluation of the chemical reactivity properties of the Pristinamycin family of antimicrobial peptides. Chem. Phys. Lett. 739, 137021. https://doi.org/10.1016/j.cplett.2019.137021 (2020).

Flores-Holguín, N., Frau, J. & Glossman-Mitnik, D. Conceptual DFT-based computational peptidology of marine natural compounds: Discodermins A-H. Molecules 25, 4158. https://doi.org/10.3390/molecules25184158 (2020).

Flores-Holguín, N., Frau, J. & Glossman-Mitnik, D. Virtual screening of marine natural compounds by means of chemoinformatics and CDFT-based computational peptidology. Mar. Drugs 18, 478. https://doi.org/10.3390/md18090478 (2020).

Flores-Holguín, N., Frau, J. & Glossman-Mitnik, D. Conceptual DFT as a helpful chemoinformatics tool for the study of the Clavanin family of antimicrobial marine peptides. In Density Functional Theory (eds De Lazaro, S. R. et al.), chap. 3, 57–67 (IntechOpen, London, 2021).

Pires, D. E. V., Blundell, T. L. & Ascher, D. B. pkCSM: Predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures. J. Med. Chem. 58, 4066–4072. https://doi.org/10.1021/acs.jmedchem.5b00104 (2015).

Liu, Y., Ma, H. & Yao, J. ER\(\upalpha\) , A key target for cancer therapy: A review. Onco. Targets. Ther. 13, 2183–2191. https://doi.org/10.2147/ott.s236532 (2020).

Amin, M. R. et al. Synthesis, antimicrobial, anticancer, PASS, molecular docking, molecular dynamic simulations & pharmacokinetic predictions of some methyl \(\upbeta\) -D-galactopyranoside analogs. Molecules 26, 7016. https://doi.org/10.3390/molecules26227016 (2021).

Schmidtke, P., Luque, F. J., Murray, J. B. & Barril, X. Shielded hydrogen bonds as structural determinants of binding kinetics: application in drug design. J. Am. Chem. Soc. 133, 18903–18910. https://doi.org/10.1021/ja207494u (2011).

Hariono, M. & Rollando, R. Molecular docking of compounds from Chaetomium sp. against human estrogen receptor alpha in searching anti breast cancer. J. Pharm. Sci. Commun. 13, 35–43. https://doi.org/10.24071/jpsc.131128 (2016).

Ervina, M. et al. In silico estrogen receptor alpha antagonist studies and toxicity prediction of melia azedarach leaves bioactive ethyl acetate fraction. J. Adv. Pharm. Technol. Res. 12, 241–246. https://doi.org/10.4103/japtr.JAPTR_198_21 (2021).

Muhammad, S. et al. Virtual screening of potential inhibitor against breast cancer-causing estrogen receptor alpha (ER\(\upalpha\) ): Molecular docking and dynamic simulations. Mol. Simul. 48, 1163–1174. https://doi.org/10.1080/08927022.2022.2072840 (2022).

Jiao, D. & Rempe, S. B. Combined density functional theory (DFT) and continuum calculations of pKa in carbonic anhydrase. Biochemistry 51, 5979–5989. https://doi.org/10.1021/bi201771q (2012).

Pradeep, S. et al. Synthesis, computational pharmacokinetics report, conceptual DFT-based calculations and anti-acetylcholinesterase activity of hydroxyapatite nanoparticles derived from Acorus calamus plant extract. Front. Chem.https://doi.org/10.3389/fchem.2021.741037 (2021).

Prasad, S. K. et al. Evaluation of Annona muricata acetogenins as potential anti-SARS-CoV-2 agents through computational approaches. Front. Chem.https://doi.org/10.3389/fchem.2020.624716 (2021).

Lobanov, M. Y., Bogatyreva, N. S. & Galzitskaya, O. V. Radius of gyration as an indicator of protein structure compactness. Mol. Biol. 42, 623–628. https://doi.org/10.1134/s0026893308040195 (2008).

Chen, H. & Panagiotopoulos, A. Z. Molecular modeling of surfactant micellization using solvent-accessible surface area. Langmuir 35, 2443–2450. https://doi.org/10.1021/acs.langmuir.8b03440 (2019).

Chikalov, I., Yao, P., Moshkov, M. & Latombe, J.-C. Learning probabilistic models of hydrogen bond stability from molecular dynamics simulation trajectories. BMC Bioinform.https://doi.org/10.1186/1471-2105-12-s1-s34 (2011).

Aungst, B. J. Optimizing oral bioavailability in drug discovery: An overview of design and testing strategies and formulation options. J. Pharm. Sci. 106, 921–929. https://doi.org/10.1016/j.xphs.2016.12.002 (2017).

Lipinski, C. A. Lead- and drug-like compounds: The rule-of-five revolution. Drug Discov. Today Technol. 1, 337–341. https://doi.org/10.1016/j.ddtec.2004.11.007 (2004).

Kelder, J., Grootenhuis, P. D. J., Bayada, D. M., Delbressine, L. P. C. & Ploemen, J.-P. Polar molecular surface as a dominating determinant for oral absorption and brain penetration of drugs. Pharm. Res. 16, 1514–1519. https://doi.org/10.1023/a:1015040217741 (1999).

Savjani, K. T., Gajjar, A. K. & Savjani, J. K. Drug solubility: Importance and enhancement techniques. ISRN Pharm. 2012, 195727. https://doi.org/10.5402/2012/195727 (2012).

Smith, D. A., Beaumont, K., Maurer, T. S. & Di, L. Volume of distribution in drug design. J. Med. Chem. 58, 5691–5698. https://doi.org/10.1021/acs.jmedchem.5b00201 (2015).

Muehlbacher, M., Spitzer, G. M., Liedl, K. R. & Kornhuber, J. Qualitative prediction of blood-brain barrier permeability on a large and refined dataset. J. Comput. Aided Mol. Des. 25, 1095–1106. https://doi.org/10.1007/s10822-011-9478-1 (2011).

Article  ADS  CAS  Google Scholar 

Becke, A. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98, 5648–5652. https://doi.org/10.1063/1.464913 (1993).

Article  ADS  CAS  Google Scholar 

Lee, C., Yang, W. & Parr, R. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 37, 785–789. https://doi.org/10.1103/physrevb.37.785 (1988).

Article  ADS  CAS  Google Scholar 

Vosko, S., Wilk, L. & Nusair, M. Accurate spin-dependent electron liquid correlation energies for local spin density calculations: A critical analysis. Can. J. Phys. 58, 1200–1211. https://doi.org/10.1139/p80-159 (1980).

Article  ADS  CAS  Google Scholar 

Adamo, C. & Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. Chem. Phys. 110, 6158–6169. https://doi.org/10.1063/1.478522 (1999).

Article  ADS  CAS  Google Scholar 

Janak, J. Proof that \(\partial E / \partial n_{i} = \epsilon\) in density functional theory. Phys. Rev. B 18, 7165–7168. https://doi.org/10.1103/PhysRevB.18.7165 (1978).

Article  ADS  CAS  Google Scholar 

Kar, R., Song, J.-W. & Hirao, K. Long-range corrected functionals satisfy Koopmans’ theorem: Calculation of correlation and relaxation energies. J. Comput. Chem. 34, 958–964. https://doi.org/10.1002/jcc.23222 (2013).

Domingo, L. R., Ríos-Gutiérrez, M. & Pérez, P. Applications of the conceptual density functional theory indices to organic chemistry reactivity. Molecules 21, 748. https://doi.org/10.3390/molecules21060748 (2016).

Domingo, L. R., Aurell, M., Pérez, P. & Contreras, R. Quantitative characterization of the global electrophilicity power of common diene/dienophile pairs in diels-alder reactions. Tetrahedron 58, 4417–4423. https://doi.org/10.1016/s0040-4020(02)00410-6 (2002).

Domingo, L. R. & Sáez, J. A. Understanding the mechanism of polar diels-alder reactions. Org. Biomol. Chem. 7, 3576–3583. https://doi.org/10.1039/B909611F (2009).

Pérez, P., Domingo, L. R., Aurell, M. J. & Contreras, R. Quantitative characterization of the global electrophilicity pattern of some reagents involved in 1,3-dipolar cycloaddition reactions. Tetrahedron 59, 3117–3125. https://doi.org/10.1016/S0040-4020(03)00374-0 (2003).

Domingo, L. R., Chamorro, E. & Perez, P. Understanding the reactivity of captodative ethylenes in polar cycloaddition reactions. A theoretical study. J. Org. Chem. 73, 4615–4624. https://doi.org/10.1021/jo800572a (2008).

Jaramillo, P., Domingo, L. R., Chamorro, E. & Pérez, P. A further exploration of a nucleophilicity index based on the gas-phase ionization potentials. J. Mol. Struct. (Thoechem.) 865, 68–72. https://doi.org/10.1016/j.theochem.2008.06.022 (2008).

Zhang, F. et al. Altertoxins with quorum sensing inhibitory activities from the marine-derived fungus Cladosporium sp. KFD33. Mar. Drugs 18, 67. https://doi.org/10.3390/md18010067 (2020).

Authors are grateful to Department of Studies in Botany and Microbiology, University of Mysore and JSSAHER for analytical facilities. SPK thankfully acknowledges the Director, Amrita Vishwa Vidyapeetham, Mysuru Campus for laboratory facilities. The authors extend their appreciation to the Deanship of Scientific research at King Khalid University for funding this work through large groups (Project under Grant Number R.G.P. 2/78/43). N.F.H. and D.G.M. are researchers of CIMAV and CONACYT and want to thank both institutions for partial support.

These authors contributed equally: Satish Anandan, Hittanahallikoppal Gajendramurthy Gowtham, CS Shivakumar, Anjana Thampy, Sudarshana Brijesh Singh, Mahadevamurthy Murali, Chandan Shivamallu, Sushma Pradeep, Natrajamurthy Shilpa, Ali A. Shati, Mohammad Y. Alfaifi, Serag Eldin I. Elbehairi, Joaquín Ortega-Castro, Juan Frau, Norma Flores-Holguín, Shiva Prasad Kollur and Daniel Glossman-Mitnik.

Department of Clinical Nutrition and Dietetics, Sri Devaraj Urs Academy of Higher Education and Research, Kolar, Karnataka, 563101, India

Satish Anandan, CS Shivakumara & Anjana Thampy

Department of PG Studies in Biotechnology, Government Science College (Autonomous), Nrupathunga Road, Bangalore, 560001, India

Department of Studies in Botany, University of Mysore, Manasagangotri, Mysore, Karnataka, 570006, India

Sudarshana Brijesh Singh & Mahadevamurthy Murali

Department of Biotechnology and Bioinformatics, JSS Academy of Higher Education and Research, Mysore, Karnataka, 570015, India

Chandan Shivamallu & Sushma Pradeep

Department of Studies in Microbiology, University of Mysore, Manasagangotri, Mysore, Karnataka, 570006, India

Biology Department, Faculty of Sciences, King Khalid University, Abha, Saudi Arabia

Ali A. Shati & Serag Eldin I. Elbehairi

Cell Culture Lab, Egyptian Organization for Biological Products and Vaccines (VACSERA Holding Company), 51 Wezaret El-Zeera St., Agouza, Giza, Egypt

Mohammad Y. Alfaifi & Serag Eldin I. Elbehairi

Department of Chemistry, University of the Balearic Islands, 07122, Palma de Mallorca, Spain

Joaquin Ortega-Castro & Juan Frau

NANOCOSMOS Virtual Laboratory, Department of Environment and Energy, Center for Research in Advanced Materials, 31136, Chihuahua, Chih, Mexico

Norma Flores-Holguín & Daniel Glossman-Mitnik

Department of Clinical Sciences, College of Veterinary Medicine, Kansas State University, Manhattan, Kansas, 66506-5606, USA

Midwest Veterinary Services, Inc., Oakland, NE, 68045, USA

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

S.A., H.G.G. and C.S.S.: Conceptualization; A.T. and S.B.S.: Methodology; C.S. and S.P.: Writing and Reviewing; N.S., M.M., A.A.S., M.Y.A., and S.E.I.E.: Visualization and Validation; J.O.C., J.F., N.F.H. and D.G.M.: Software; M.M., S.P.K., C.S. and D.G.M.: Supervision and reviewing the final manuscript.

Correspondence to Mahadevamurthy Murali, Chandan Shivamallu, Shiva Prasad Kollur or Daniel Glossman-Mitnik.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Anandan, S., Gowtham, H.G., Shivakumara, C.S. et al. Integrated approach for studying bioactive compounds from Cladosporium spp. against estrogen receptor alpha as breast cancer drug target. Sci Rep 12, 22446 (2022). https://doi.org/10.1038/s41598-022-22038-x

DOI: https://doi.org/10.1038/s41598-022-22038-x

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Plastic Banquet Tables Sign up for the Nature Briefing: Cancer newsletter — what matters in cancer research, free to your inbox weekly.