1. Introduction
Diseases are nowadays one of the major problems facing world populations in general and those in sub-Saharan Africa in particular [1] [2] . To cope with diseases, more often chemical entities are used. Thus, several chemical entities have been discovered as drugs and are grouped into different categories based on their actions on the disease. Among them are antibiotics and anti-inflammatories [3] [4] [5] . Unfortunately, these drugs often date back many years, and as a result, resistance has been developed by pathogens that cause diseases. This is currently a major public health concern that needs to be taken into consideration by both governments and researchers [6] [7] . The most promising approach to fighting with these resistances is the more than urgent development of new therapeutic formulas. These new drugs must interact with well-known target proteins as well as the search for new targets that might lead to new successful therapies [8] [9] . Tuberculosis (TB) is a highly contagious disease after COVID-19 caused by Mycobacterium tuberculosis (Mtb). Immunodepressed patients such as the HIV-positive individuals are more susceptible to TB, and its major concern is due to the high resistance of the bacterium to the available drugs such as isoniazid, rifampicin, etc. [10] [11] [12] . The World Health Organization aims to reduce mortality rates and eradicate TB under 2030 by improving access to interventions [13] . TB is a leading cause of death among infectious diseases worldwide, resulting in a significant global public health burden. In 2022, there were 1.30 million deaths worldwide due to TB, with 1.13 million deaths among HIV-negative individuals and approximately 167,000 deaths among HIV-positive individuals. Although there was a 19% reduction in TB deaths between 2015 and 2022, the TB strategy’s goal of a 75% reduction by 2025 was not met [14] . The lungs are affected primarily by Mtb and become active when the immune system weakens. This triggers a fight between the hosts and the pathogen, releasing free radicals (ROS) and survival mechanisms used by Mtb. Patients with TB present more levels of ROS than healthy individuals. Immune response and inflammation keep a significant role in TB disease [15] . TB is a disease that can be treated with antibiotics. Currently active TB is treated with several regimens that can take 4, 6, or 9 months depending on the treatment plan. The available TB treatment regimens include the 4-month Rifapentine-moxifloxacin TB Treatment Regimen and the 6- or 9-month RIPE (rifampicin, isoniazid, pyrazinamide, and ethambutol) TB Treatment Regimen. However, there is a risk of drug resistance [16] [17] [18] . So to combat this, more toxic drugs are available, but they are not easily accessible to those with lower incomes. Additionally, the treatment itself has several unwanted side effects, such as liver toxicity, gastrointestinal issues, skin rashes, eyes problem, and heart complications. It is essential to improve the management of TB patients by reducing these side effects and increasing the effectiveness of the drugs. Despite the increasing worldwide incidence of TB and its alarming threat toward the public health, no novel anti-TB drugs have been introduced into clinical practice over the past decades. The impact of ever-increasing drug resistance, the serious side effects of some current anti-TB drugs, and the lack of efficacy of current treatments in immunodepressed patients combine to make the development of new antimycobacterial agents an urgent priority. One key approach is to combine anti-TB antibiotics with agents that target the antioxidant, which may help to eliminate the disease more quickly. However, research in this area is still in the early stages of evaluation [19] . Previous studies show that the combination of therapeutic supplements such as vitamin C and cysteine has a rapid inhibitory effect on bacteria, preventing their proliferation in the host [20] .
This paves the way for research into natural tropical antioxidants that are accessible to all. Research has been carried out on halimane diterpenoids isolated from plant sources (Marchantiophyta, Magnoliophyta) and animal sources (Porifera, Mollusca, Chordata, Cnidaria) showing their potential to inhibit the virulence factors of the Mtb responsible for this disease [21] . Three labdane-type of diterpenoids (18-nor-labd-13(E)-ene-8α, 15-diol, labd-13(E)-ene-8α, 15-diol and austroinulin) has been extracted from the stem bark of Croton Sylvaticus [22] . From DPPH free radical scavenging, Frum et al. have previously described the relevant antioxidant activity of the extract of the stem bark of Croton Sylvaticus [23] . Research on the ethnobiological activity on plants shows that South African traditional healers used this plant to treat the TB and mental disorder [24] . These findings are a proof of the double competence of these extracts from the stem bark as a tonic and as a therapeutic agent for TB. In addition, the requirements for the incorporation of diastereomers of the same halimane diterpenoids into other halimane diterpenoids in clearly defined ratios blocks the phagolysome maturation and macrophage phagocytosis in human cells [25] . For instance, the compound tuberculosinol can be added in the proportion (1:1) to the mixture of the diastereoisomers isotuberculosinol (R) and isotuberculosinol (S) in turn in the proportion (1:3) [26] . Minor studies on the stereo-clarification of halimane diterpenoids discovered to date have been carried out by Roncero et al. [21] .
Despite these advances, an analysis of the impact of their biological, pharmacokinetic and medicinal activities is lacking. Among 42 halimane diterpenoids compound selected, we have investigated the better anti-TB druglikeness able to block the development and maturation of the Mtb. Consequently, molecular docking and dynamic simulation have been performed. The pharmacokinetic parameters of this elected anti-TB druglikeness estimated have compare to those of three first-line treatments: ethambutol Isoniazid and amikacin.
2. Materials and Methods
2.1. Ligands Preparation and Selection
Forty-two (42) halimane diterpenoids compounds of natural origin, as they occur in African flora. The absolute configurations are presented in Table 1 [21] . They are carried out with under force field in the Open Babel plug-in tool of PyRx 0.8 tool to obtain optimized geometric [27] . After that, they became as
![]()
Table 1. Molecular registration and substituent specification.
ligand and automatically converted to pdbqt format.
2.2. Targets Selection and Preparation
Two validated Mtb drug targets namely Homoserine transacetylase (PDB ID 8F2L) and Mycothiol S-transferase enzyme (PDB ID 8F5V) have been selected according to their potential roles in Mtb maturation besides development [28] [29] . Their 3D co-crystallized structures have been retrieved in PDB code with respectively PDB ID 8F2L and 8F5V. Each of the two selected targets has been obtained by X-ray diffractions with resolution 2.89 AÅ and 1.45 AÅ respectively.
The two receptors were prepared for docking by removing all water molecules, co-crystallized ligands, cofactor ions, as well as charges and hydrogens using the BIOVIA Discovery Studio 2021 client software [30] . Opened in the PyRx virtual screening tool, everyone is treated as a macromolecule structure and saved in pdbqt format [31] .
2.3. Molecular Docking Validation and Study Process
1) Description
In the field of computer-aided drug design, molecular docking is popular technique employed to predict the affinity between the three-dimensional structures of a small molecule (ligand) and a macromolecule (protein). This method uses mathematical algorithms to study the interactions of natural compounds with specific protein targets such as mycobacterial enzymes. It is based on the binding affinity energy of compounds with proteins. The result is the formation of complexes, the most stable of which corresponds to the lowest affinity energy [32] [33] .
2) Validation of Docking Protocol
Separation of the co-crystallized ligand from each protein and repetition of the docking studies in the same region with the re-docked nature ligand is a useful way of validating molecular docking. The docking simulations of the orientation and position of the ligands used have been validated by redocking experiments, through a representation of the valid and reasonable potential binding mode of the inhibitors. The molecular structure of halimane derivatives was optimized using the Gaussian 09 and Gauss View 6.0 suite programs, with density functional theory (DFT) and Becker’s three-parameter calculation method (B3) combined with Lee, Yang, and Parr’s non-local correlation functional (LYP) and a 6-311++G (d, p) basis set [34] [35] [36] [37] [38] . DFT was chosen as it provides reliable results in predicting molecular geometry and properties, such as total energy and molecular shape [37] [39] . The 6-311++G (d, p) basis set was selected due to its ability to account for diffusion, polarization, and angular flexibility in hydrogen-bonded organic systems [40] . All geometric optimizations were followed by determination of the vibrating frequencies. Before carrying out molecular docking experiments on the Mtb enzymes Homoserine transacetylase and Mycothiol S-transferase, a protocol was developed to focus on the active site of these enzymes and take into account surrounding amino acid residues. Cofactors, water molecules and co-crystallized ligands were removed from the enzymes. The co-crystallized ligand was docked into the protein structure to determine the best conformation using the calculated RMSD value [41] . Each amino acid residue interacting with the ligand was realigned using PyMol software [42] .
Virtual screening was then carried out by converting Mtb ligands and enzyme receptors to PDBQT format using PyRx software [43] . Ligands with the best binding energy and lowest RMSD value were selected. Finally, the hydrogen, hydrophobic and protein-ligand interactions of the best ligands were studied using Discovery Studio [30] .
2.4. Pharmacokinetic and Drug-Likeness Study
In the biological system, the pharmacokinetic parameters of candidate therapeutic molecules, such as their ADMET are crucial to the selection of those with drug potential activity. This selection enables the creation and development of new innovative drugs. Availability of a therapeutic compound for administration and absorption in the body is an essential aspect of drug-likeness. It is based on evaluating the most suitable candidates in accordance with the ADMET’s rules. Maestro Schrödinger software and the online resource Swiss-ADME were used to perform the ADMET properties of the compounds [44] [45] [46] .
2.5. DFT Studies
The electrophilic reactivity and stability of the optimized halimane compounds was classified according to their electrophilicity index (ω). The index was calculated using the following expression in Equation (1).
(1)
with the variables μ, η, S, and χ correspond to chemical potential, hardness, global softness, and electronegativity (see Equation (1)). They are determined by expressing the energy of the HOMO and LUMO orbitals using the following formulae [47] .
(2)
(3)
(4)
(5)
2.6. Molecular Dynamics Simulation (MDs)
The MDs is a study conformational changes and the stability of ligand-binding complexes over simulation time. These simulations were conducted using the Desmond simulation package in Schrödinger software, based on two complexes obtained from protein docking [48] . These simulations were performed at standard temperature and pressure, with a 100 ns duration and snapshots recorded every 100 ps. Short-range interactions were analyzed using OPLS4 force field parameters and a cutoff radius of 9.0 Å, while long-range electrostatic interactions were calculated using the Particle-Mesh Ewald approach [49] . The Nose-Hoover and Martyna-Tobias-Klein chain coupling scheme was used for constant temperature and pressure conditions. Water molecules were modeled using a simple point-charge model. This analysis included evaluating factors of root mean square deviation (RMSD) and fluctuation (RMSF), as well as examining ligand contacts with amino acids using the Desmond interaction diagram [50] [51] .
2.7. Interaction Energy Calculations
The unbound interactions between the ligand compound and enzyme residues were analyzed using the Discovery Studio software to calculate the interaction energy values [30] . In addition, the molecular modeling interaction values between the ligand compound and enzyme residues in unbound interactions were determined using the same software such describe in our previous work [52] [53] . This study utilized the interaction energy protocol with the CHARMm force field and dielectric constant to identify electrostatic and van der Waals forces [54] . By analyzing the interaction energy contributions per active site residue, the significance of individual interactions was identified and facilitated comparative analysis. This method allows for the selection of affinity values that predict favorable and unfavorable PEP substitutions.
3. Results and Discussion
3.1. Docking Procedure
3.1.1. Validation Protocol
To create a model able to correlate the structure of candidate molecules with their ability to bind to receptor crystal structures (PDB ID: 8F2L and 8F5V) [55] . For each structure, we started by testing a combination of parameters to best reproduce (with high accuracy) the native ligand conformation, in order to create a reliable protocol. In this context, the ligand was removed and then re-docked according to the rigid docking process with the crystal structure. The RMSD obtained between the native ligand and the realigned ligand in Homoserine transacetylase and Mycothiol S-transferase enzyme was respectively 1.326 Å and 3.404 Å.
3.1.2. Molecular Docking Analysis of the Mtb’s Protein
Now that the submitted protocol has been validated, candidate molecules are subjected to molecular docking to determine the best binding mode for the ligand-Mtb enzyme protein complex. Each ligand that docked into the binding site of the Mtb protein produced different poses. Molecular docking of 42 halimane diterpenoids produced the highest affinity of ligand for binding to each Mtb enzyme (Table 2). Each ligand (target molecule)-protein bond is stable due to this binding affinity. Adjustment of the binding pocket is best when the binding affinity value is negative. For this reason, binding affinity represents a more efficient binding pocket fit [56] .
In view of the higher calculated affinity results in complex Homoserine transacetylase and its top 5 position in complex Mycothiol S-transferase enzyme, (2R,5R,6S)-1,2,3,4,5,6,7,8-octahydro-5-((E)-5-hydroxy-3-methylpent-3-enyl)-1,1,5,6-tetramethylnaphthalen-2-ol (compound 20) can be selected as the most efficient anti-TB halimane diterpenoids in the series. The different energy parameters of the ligand 20-enzyme Mtb bond in each of the crystal structures are presented in Table 2. To analyze the electrostatic effect of the ligand, the ligand of compound 20 and every one of the Mtb enzymes proteins was generated using Discovery Studio Visualizer software [57] .
![]()
Table 2. Energy value of the binding affinity of each enzyme with each compound and ligand.
In Figure 1, 2D graph shows the various molecular interactions between ligand compound 20, the best potential anti-TB drug, and each of the target enzymatic proteins.
1) Interaction with Mycothiol S-transferase enzyme
Figure 1(a) shows three dark green circles, each connected to a dotted line of the same colour, representing the three conventional hydrogen bonds formed between:
a) The oxygen atom of the terminal hydroxide group of the 5-hydroxy-3-methylpent-3-enyl substituent bonded to carbon 5 of the naphthalene ring and an acceptor hydrogen atom of the X (N or O)-H bond of THR 89;
b) The hydrogen atom of the hydroxide group of the same substituent and a donor atom of the X (N or O) -H bond of THR 89;
c) The hydrogen atom of this terminal group of the same substituent bonded to carbon 5 of the naphthalene ring and a donor atom of the X (N or O) -H bond of TYP 91. An analysis of the interaction of ligand 20 with Mtb revealed pi-sigma bonds with the acid TRP 139 (purple color).
Additionally, the dark purple indicates the interaction of the Pi-sigma bond with the amino acid TRP 139, which also forms alkyl and pi-alkyl bonds shown in light purple. The amino acids THR 91, VAL 56 and ILE 59 are also colored in the same way.
2) Interaction with Homoserine transacetylase enzyme
In Figure 1(b), the three conventional hydrogen bonds formed are materialized by:
![]()
![]()
Figure 1. Interaction diagram of (2R,5R,6S)-1,2,3,4,5,6,7,8-octahydro-5-((E)-5-hydroxy-3-methylpent-3-enyl)-1,1,5,6-tetramethylnaphthalen-2-ol with Mycothiol S-transferase enzyme (a) and Homoserine transacetylase enzyme (b).
a) The hydrogen atom of the hydroxide group of the 5-hydroxy-3-methylpent-3-enyl substituent bonded to carbon 5 of the naphthalene ring and a donor atom of the X (N or O)-H bond of PRO 77;
b) The oxygen atom of the hydroxide group bonded to carbon 2 and an acceptor hydrogen atom of the X (N or O)-H bond of GLY 62;
c) The hydrogen atom of the same group bonded to carbon 2 of the naphthalene ring (Figure 2(b)) and a donor atom of the X (N or O) -H bond of THR 61. The final oxygen atom of the substituent 5-hydroxy-3-methylpent-3-enyl connected to carbon 5 of the naphthalene ring and the amino acid TRP 79 also form a carbon-hydrogen bond.
The amino acids contributing to the pi-alkyl and alkyl bonds formation are as for ALA 59: three, and one each for LEU 354, TYR 234 and TRP 79.
3.2. Results of Modelling and Optimization Structures
Evaluation of the chemical library’s molecules by kinetic stability and chemical reactivity requires analysis of the HOMO-LUMO gap presented in Table 3. This descriptor provides a measure of the energies of electronic interactions between orbitals. During the interactions, the frontier orbitals can be electron acceptors for the high altitude electronically unoccupied or electron donors for the low altitude electronically occupied orbitals [58] . A molecule with a high HOMO-LUMO energy gap has low chemical reactivity and high kinetic stability [59] . Compounds considered unstable and highly reactive have a low HOMO-LUMO energy gap.
In Table 3, we show the comparison between the HOMO-LUMO energy gap in gas and water. In gas, the orbitals energies are 6.855 eV. The electronic transition from the HOMO frontier orbital to its LUMO counterpart in this absolute configuration of this compound proves to be therefore very difficult. In aqueous media, there is an average decrease of 1983 eV compared with gaseous media. This decrease in the energies of the frontier orbitals is due to the solvation effect claimed by De Proft and Geerlings [60] [61] .
A more extensive kinetic assessment of stability through the different global reactivity descriptors (electronegativity χ, global hardness (η), global softness (S), and electrophilicity index (ω)) was also performed.
The information on the resistance to electron change by a deformation of the electron cloud caused by a small disturbance from the chemical reaction is found from the hardness (η). The selection of the most kinetically stable compound is made according to the criterion of the highest hardness. The average achieved hardness difference is equal to 1585 eV. Our results prove that solvation in water significantly decreases the hardness of the studied molecules. This reduction in hardness values occasioned by the solvation of neutral molecules has been claimed earlier by De Proft and Geerlings [61] . We have measured the deterioration in binding energy due to maximum electron flow between a donor and an acceptor using the overall electrophilicity index (ω). This descriptor, through the measurement of the capacity of a system to gain an electron, enables the prediction of the chemical reactivity of the compound. In gas phase, the highest electrophilicity index is 1.460 for 2.393 in water. The compounds in aqueous solution show that solvation enhances the reactivity. A correlation between the change in the global reactivity descriptor and the solvation energy observed by Parr et al. in the framework of the reaction field theory has been confirmed [62] .
3.3. Pharmacokinetic and Drug-Likeness Study
Drug-likeness refers to the prediction of whether a specific compound can be developed into a drug. In order to determine drug-likeness, various physiochemical properties of the compound are analyzed using established rules and properties [63] .
In 1997, Lipinski and colleagues established the “rule of five”, a filter commonly used to determine whether a molecule is likely to be absorbed orally. According to this rule, to be considered absorbable, a molecule must not exceed
![]()
Table 3. HOMO-LUMO gap, electronegativity χ, global hardness (η), global softness (S), and electrophilicity index (ω) in eV of some compounds of halimane diterpenoids.
a certain molecular weight (500 Da), octanol/water partition coefficient (5), number of hydrogen bond donors (5), and number of hydrogen bond acceptors (10) [64] . The Egan filter proposes a maximum lipophilicity (5.88) and total polar surface area (131) for optimal oral bioavailability [65] . Veber and colleagues suggest a maximum number of rotating bonds (10) and total polar surface area (140) for good oral bioavailability [66] . Muegge and his team established criteria including molecular weight (200 - 600), lipophilicity (between −2 and −5), total polar surface area (≤150), number of rings (≤7), number of carbons (>4), number of heteroatoms (>1), number of rotatable bonds (≤15), hydrogen bond acceptors (≤10), and hydrogen bond donors (≤5) [67] . The pharmacokinetic study was based on prediction of blood-brain barrier (BBB) crossing with criterion values (between −3.0 and 1.2), prediction of Caco-2 permeability (<25 low and >500 high), percentage oral absorption in humans with values 1, 2, or 3 for low, medium, or high, skin permeability (−8.0 to −1.0) [68] . According to the results of the analysis in Table 4, compound 20 has appropriate physicochemical parameters to be considered as a drug crossing the brain barrier (BBB).
Compounds derived from therapeutic treatment of TB meet the criteria for drug eligibility. Compound 20 also shows good results. With regard to Human Oral Absorption (HOA), which is strongly recommended by the WHO as an absorption criterion, compound 20 has good absorption potential. It is more effective at passively crossing the intestinal permeability barrier than first-line treatments (ethambutol, isoniazid).
3.4. Molecular Dynamic Simulation Results
Molecular dynamics simulation (MDs) is an essential computational tool for understanding the physical and chemical properties of drug candidate compounds. It is seen as a study of relative positional changes of a moving molecule system in a protein over time [69] .
MDs first bring to life biomolecular structures in solution at different timelines, giving a view of natural dynamics. Secondly, it allows us to find values for several molecular properties as close as possible to experimental data. Finally, it explores the different conformations of molecular structures in space [70] .
![]()
Table 4. Assessment of the pharmacokinetic profiles of compound 20 and anti-TB drugs therapy (ethambutol Isoniazid and amikacin), using the recommended screening for 95% of known drugs.
Our study will examine the stability of docking-validated complexes to see how they perform in MDs.
The Root Mean Square Deviation (RMSD) serves as a metric for evaluating the stability of a system. Specifically, it measures the average displacement of a set of atoms relative to a reference frame. For globular and small proteins, the permissible range of structural conformations falls between 1 and 31. The RMSD evolution between the Mtb Homoserine transacetylase and the ligand (compound 20) is depicted in Figure 2(a). Between 0 and 40 nsec, the Lig fit Prot curve is far from the reference protein scaffold, suggesting instability within the complex. However, from 40 to 60 nsec, there is an alignment with the reference protein backbone, indicating structural stability. Beyond 60 nsec, there is an RMSD deviation as a function of time, indicating instability of ligand-protein conformations. In contrast, Figure 3(a) illustrates the stable structural conformation of the ligand and the atoms of the reference protein Mycothiol S-transferase enzyme. Although there is a sudden shift away from each other between 55 and 90 nsec, reflecting instability of the complex, the two structures converge once more in the last 10 nsec. Both complexes exhibit an RMSD of less than 3.0 Å, indicating their stability.
In accordance with the analysis of the protein-RMSF graph, the RMSF value proves to be a useful tool for assessing the fluctuations of residues in molecular dynamics simulations. It is particularly beneficial in evaluating the quality of amino acid binding at a binding site. RMSF values that are high suggest suboptimal binding. Upon examining the protein-RMSF plot depicted in Figure 2(b), it is apparent that the interaction between the compound 20 and the protein’s residue, as indicated by the green lines, is stable, with values below 1.8 Å. Conversely, residue fluctuation values on the protein that are greater than 2 Å suggest a weak interaction between the ligand and the protein. Figure 3(b) also illustrates the stability of the compound 20 and residue, although the RMSF value is 2.2 Å, which is below the cut-off value of 2.5 (vita). However, all protein-mediated fluctuations remain below 2 Å, indicating a strong interaction.
According to Figure 2(c) show the protein amino acids that interact with the ligand. Therefore, we have THR 61, SER 157 and HIS 350 for hydrogen bonds. The amino acids involved in hydrophilic interactions are as follows: ALA 59, LEU 60, ALA 230 TYP 234, TYR 282, LEU 289, LEU 322, PHE 353 and LEU 354. We have selected ALA 59, LEU 60, THR 61, ASP 63, HIS 65, THR 76, PRO 77, SER157, ARG 227 and ASP 351 for protein-ligand interactions with hydrogen in water (water bridges). For Figure 3(c), hydrogen bonds are established between GLN 60, HIS 87, ASP 88, THR 89, TYR 91, TRP 139, ASN 140, SER 151, ASP 154 and ASP 155. In terms of hydrophobic bonds, the ligand interacts with VAL 56, ILE 59, TYR 91, TRP 139, ASN 140, PRO 141 and ILE 152. Atoms of the ligand and the acids HIS 52, GLN 60, ARG 86, HIS 87, ASP 88, THR 89, GLY 90, TYR 91, SER 138, TRP 139, ASN 140, PRO 141, ALA 147, ARG 148, SER 151, ASP 154 and ASP 155 form water bridges. The ligand and ASN 140 bond in polar mode.
![]()
Figure 2. Simulation Interactions graphs of compound 20 (ligand 8F2L (protein)): (a) Protein-Ligand RMSD; (b) Protein RMSF. (c) Protein-Ligand Contacts.
3.5. Interaction Energy Calculations
The study of the non-bonded interactions between each residue of the receptor and the small molecule revealed the individual contribution of each residue to the total interaction energy. This was done by calculating the sum of the van der Waals and electrostatic interaction energies of each residue. This analysis allowed us to understand the role of each residue in binding to the small molecule.
![]()
Figure 3. Simulation Interactions graphs of compound 20 (ligand) −8F5V (protein): (a) Protein-Ligand RMSD. (b) Protein RMSF. (c) Protein-Ligand Contacts.
Several residues showed lower interaction energy values ranging from 0 to −12.0 kcal/mol. The results allow us to rank our residues in descending order according to their protein-ligand interaction stability: THR 89, GLY 90, ARG 148, VAL 56, GLN 60. This compensated for the perturbing contribution energies of residues such as ILE 59, HIS 87 and ASP 88 whose total energy was greater than.
![]()
Figure 4. Interaction energies (electrostatic, van der waals, total) of some residues of Mycothiol S-transferase enzyme and ligand.
4. Conclusion
Evaluation of the molecules in our chemical library by molecular docking shows that (2R,5R,6S)-1,2,3,4,5,6,7,8-octahydro-5-((E)-5-hydroxy-3-methylpent-3-enyl)-1,1,5,6-tetramethylnaphthalene-2-ol is the best ligand for binding affinity with the Mtb Mycothiol S-transferase and Mtb Homoserine transacetylase enzymes, respectively. Molecular dynamics simulations reveal that the interaction between the crystalline structure of the Mtb Mycothiol S-transferase enzyme and is the most effective over time. Analysis of kinetic stability and chemical reactivity through the HOMO-LUMO gap study shows good stability of the compound with excellent stability of the series in both media. For pharmacokinetic requirements, our molecule meets all the criteria. Further studies can be carried out, such as in vitro and in vivo evaluations are recommended.
Acknowledgements
The author is immensely grateful to Prof. Dr. Fidele Ntie-Kang for his technical assistance during this work.
Author Contributions
Conceptualization, B.D.B. and B.M.D.; methodology, E.A.L.G., B.D.B., B.M.D., P.A.O., and M.M.L.L.; writing-original draft preparation: E.A.L.G. and B.D.B.; writing-review and editing: B.M.D. L.C.O.O. and M.M.L.L. All authors have read and agreed to the published version of the manuscript.
Statement of Usage of Artificial Intelligence
Not applicable.
Data Availability Statement
Not applicable.
Funding Information
Not applicable.