Displaying publications 21 - 39 of 39 in total

Abstract:
Sort:
  1. Kwiatkowski M, Hameed BH
    J Mol Model, 2017 Aug;23(8):238.
    PMID: 28735499 DOI: 10.1007/s00894-017-3402-0
    The paper presents the results of research devoted to reliability evaluation of the analysis of results of the porous structure of activated carbons based on incomplete nitrogen adsorption isotherms using the BET, t-plot, and NLDFT methods, as well as the LBET method comprising the unique numerical fast multivariant procedure of adsorption system identification. The research involved the application of the nitrogen adsorption isotherms obtained for five samples of activated carbons produced from waste materials of organic origin by way of chemical activation with potassium hydroxide, sodium hydroxide, and potassium carbonate with the use of microwave heating. The analyses performed pointed to a good correlation between the results obtained using the BET, t-plot, NLDFT, and LBET methods. Moreover, the parameters of the porous structure determined using these methods based on incomplete adsorption isotherms of nitrogen are in fact as reliable as these methods allow.
  2. Alipour E, Alimohammady F, Yumashev A, Maseleno A
    J Mol Model, 2019 Dec 13;26(1):7.
    PMID: 31834504 DOI: 10.1007/s00894-019-4267-1
    Today, drug delivery systems based on nanostructures have become the most efficient to be studied. Recent studies revealed that the fullerenes can be used as drug carriers and transport drugs in a target cell. The aim of the present work is to study the interaction of C60 fullerene containing porphyrin-like transition metal-N4 clusters (TMN4C55, TM = Fe, Co, and Ni) with a non-steroidal anti-inflammatory drug (ibuprofen (Ibp)) by employing the method of the density functional theory. Results showed that the C60 fullerene with TMN4 clusters could significantly enhance the tendency of C60 for adsorption of ibuprofen drug. Also, our ultraviolet-visible results show that the electronic spectra of Ibp/TMN4C55 complexes exhibit a blue shift toward lower wavelengths (higher energies). It was found that the NiN4C55 fullerene had high chemical reactivity, which was important for binding of the drug onto the carrier surface. In order to gain insight into the binding features of Ibp/TMN4C55 complexes, the atoms in molecules analysis was also performed. Our results exhibit the electrostatic features of the Ibp/TMN4C55 bonding. Consequently, this study demonstrated that the TMN4C55 fullerenes could be used as potential carriers for delivery of Ibp drug in the nanomedicine domain. Graphical Abstract The TMN4C55 (TM=Fe, Co, and Ni) fullerenes could be used as potential carriers for delivery of ibuprofen drug in the nanomedicine domain.
  3. Mosapour Kotena Z, Razi M, Ahmadi S
    J Mol Model, 2021 Oct 09;27(11):315.
    PMID: 34625848 DOI: 10.1007/s00894-021-04916-9
    Rare sugars are monosaccharides with tremendous potential for applications in pharmaceutical, cosmetics, nutraceutical, and flavors industries. The four rare sugars, including gulose, allose, altrose, and talose, are stereoisomers that are different in the hydroxyl group orientation (axial or equatorial) on the C2-4 atoms. The basis sets effect in evaluation of the possibility intramolecular hydrogen bonding (H-bonds) in the selected rare sugars was studied from 6-31G* to 6-311 ++ G(d,p) basis sets using DFT, AIM, and NBO methods. The results show that the selected rare sugars are more stable at 6-311 ++ G(d,p) basis sets compared to 6-31G* because their electronic energies were reduced between 158 and 164 (kcal.mol-1). The overall effect of basis set enhancement is to decrease H-bond energies in the range of  1.25 to 2.51 (kcal.mol-1) and stabilization energies between 2 and 5 (kcal.mol-1) in the selected rare sugars at the DFT level of theory. The intramolecular H-bond distances, H-bond energies obtained from the AIM analysis, and also the second-order stabilization energies obtained from the NBO analysis were fluctuated largely depending on the basis set. In summary, it was found that the use of 6-311 ++ G(d,p) basis set to be more efficient results in rare sugars geometry than the 6-31G* basis set.
  4. Nwaefulu ON, Al-Shar'i NA, Owolabi JO, Sagineedu SR, Woei LC, Wai LK, et al.
    J Mol Model, 2022 Oct 04;28(11):340.
    PMID: 36194315 DOI: 10.1007/s00894-022-05326-1
    Cancer is imposing a global health burden because of the steady increase in new cases. Moreover, current anticancer therapeutics are associated with many drawbacks, mainly the emergence of resistance and the severe adverse effects. Therefore, there is a continuous need for developing new anticancer agents with novel mechanisms of action and lower side effects. Natural products have been a rich source of anticancer medication. Cycleanine, a natural product, was reported to exert an antiproliferative effect on ovarian cancer cells by causing apoptosis through activation of caspases 3/7 and cleavage of poly (ADP-ribose) polymerase to form poly (ADP-ribose) polymerase-1 (PARP1). It is well-established that PARP1 is associated with carcinogenesis, and different PARP1 inhibitors are approved as anticancer drugs. In this study, the cytotoxic activity of cycleanine was computationally investigated to determine whether it is a PARP1 inhibitor or a caspase activator. Molecular docking and molecular dynamics (MD) simulations were utilized for this purpose. The results showed that cycleanine has a good binding affinity to PARP1; moreover, MD simulation showed that it forms a stable complex with the enzyme. Consequently, the results showed that cycleanine is a potential inhibitor of the PARP1 enzyme.
  5. Aissaoui T, Benguerba Y, AlOmar MK, AlNashef IM
    J Mol Model, 2017 Sep 14;23(10):277.
    PMID: 28913646 DOI: 10.1007/s00894-017-3450-5
    Recently, there has been significant interest in the possibility of using deep eutectic solvents (DESs) as novel green media and alternatives to conventional solvents and ionic liquids (ILs) in many applications. Due to their attractive properties, such as their biodegradability, low cost, easy preparation, and nontoxicity, DESs appear to be very promising solvents for use in the field of green chemistry. This computational study investigated six glycerol-based DESs: DES1(glycerol:methyl triphenyl phosphonium bromide), DES2(glycerol:benzyl triphenyl phosphonium chloride), DES3(glycerol:allyl triphenyl phosphonium bromide), DES4(glycerol:choline chloride), DES5(glycerol:N,N-diethylethanolammonium chloride), and DES6(glycerol:tetra-n-butylammonium bromide). The chemical structures and combination mechanisms as well as the sigma profiles and sigma potentials of the studied DESs were explored in detail. Moreover, density, viscosity, vapor pressure, and IR analytical data were predicted and compared with the corresponding experimental values reported in the literature for these DESs. To achieve these goals, the conductor-like screening model for realistic solvents (COSMO-RS) and the Amsterdam Density Functional (ADF) software package were used. The predicted results were found to be in good agreement with the corresponding experimental values reported in the literature. Further theoretical investigations are needed to confirm the experimental results-regarding both properties and applications-reported for these DESs.
  6. Mustafa SF, Maarof H, Ahmed R, Abdallah HH
    J Mol Model, 2016 Dec;22(12):290.
    PMID: 27866329
    Diffusion at the atomic or molecular level is a source of many physical, chemical, and biological processes taking place in plentiful materials. This work is an endeavor toward investigating the diffusional behavior of two different type of guests, hexadecane-1,16-diol and hexadecane enclathration in urea tunnel architecture, whereby the correlation of the diffusion mechanism with the guest's structural and conformational properties is explored. To carry out this study, molecular dynamics simulation approach is adopted. It is found that hexadecane-1,16-diol exhibit slower diffusion with an average diffusion coefficient value [Formula: see text], where hexadecane diffuse more rapidly with an average diffusion coefficient value [Formula: see text]. It is also observed that the structural properties influence the guest's travel distance and torsion angle distribution of the trans and gauche conformational proportion. Furthermore, the observed high energy barrier accounted for hexadecane-1,16-diol and low energy barrier for hexadecane along urea tunnel systems was analyzed. The comparison of our obtained results are in close agreement with the available experimental measurements, i.e., gauche proportion properties between two different guest molecules correlate well with Raman spectroscopy investigation on α,ω-dihalogenoalkane/urea inclusion compounds. Our calculations also successfully endorse the structure-property relation between the two systems.
  7. Jusoh N, Zainal H, Abdul Hamid AA, Bunnori NM, Abd Halim KB, Abd Hamid S
    J Mol Model, 2018 Mar 15;24(4):93.
    PMID: 29546582 DOI: 10.1007/s00894-018-3619-6
    Recent outbreaks of highly pathogenic influenza strains have highlighted the need to develop new anti-influenza drugs. Here, we report an in silico study of carvone derivatives to analyze their binding modes with neuraminidase (NA) active sites. Two proposed carvone analogues, CV(A) and CV(B), with 36 designed ligands were predicted to inhibit NA (PDB ID: 3TI6) using molecular docking. The design is based on structural resemblance with the commercial inhibitor, oseltamivir (OTV), ligand polarity, and amino acid residues in the NA active sites. Docking simulations revealed that ligand A18 has the lowest energy binding (∆Gbind) value of -8.30 kcal mol-1, comparable to OTV with ∆Gbind of -8.72 kcal mol-1. A18 formed seven hydrogen bonds (H-bonds) at residues Arg292, Arg371, Asp151, Trp178, Glu227, and Tyr406, while eight H-bonds were formed by OTV with amino acids Arg118, Arg292, Arg371, Glu119, Asp151, and Arg152. Molecular dynamics (MD) simulation was conducted to compare the stability between ligand A18 and OTV with NA. Our simulation study showed that the A18-NA complex is as stable as the OTV-NA complex during the MD simulation of 50 ns through the analysis of RMSD, RMSF, total energy, hydrogen bonding, and MM/PBSA free energy calculations.
  8. Gagui S, Meradji H, Ghemid S, Naeem S, Haq BU, Ahmed R, et al.
    J Mol Model, 2022 Dec 21;29(1):14.
    PMID: 36542150 DOI: 10.1007/s00894-022-05412-4
    Understanding the physical properties of a material is crucial to know its applicability for practical applications. In this study, we investigate the phase stability, elastic, electronic, thermal, and optical properties of the ternary alloying of the scandium and yttrium nitrides (Sc1-xYxN) for different compositions. To do so, we apply a "density functional theory (DFT)" based scheme of calculations named as "full potential (FP) linearized (L) augmented plane wave plus local orbitals (APW + lo) method" realized in the WIEN2k computational package. At first, the phase stability of the investigated compositions of the mentioned alloy is determined. The analysis of our calculations shows that Sc1-xYxN alloy is stable in rock salt crystal structure for all investigated compositions. Next to that, the elastic properties of the rock-salt phase of the studied ternary alloy Sc1-xYxN at all above said compositions were done at the level of "Wu-Cohen generalized gradient approximation (Wu-GGA)" within DFT. However, Trans-Blaha (TB) approximation of the "modified Becke-Johson (mBJ)" potential is also used in combination with Wu-GGA where the thermal properties are calculated at the level of the "quasi-harmonic Debye model." The obtained results for the absorption coefficients, and optical bandgap, represent that the title alloy may be a suitable candidate for the applications in optoelectronic devices.
  9. Kueh R, Rahman NA, Merican AF
    J Mol Model, 2003 Apr;9(2):88-98.
    PMID: 12707802
    The arginine repressor (ArgR) of Escherichia coli binds to six L-arginine molecules that act as its co-repressor in order to bind to DNA. The binding of L-arginine molecules as well as its structural analogues is compared by means of computational docking. A grid-based energy evaluation method combined with a Monte Carlo simulated annealing process was used in the automated docking. For all ligands, the docking procedure proposed more than one binding site in the C-terminal domain of ArgR (ArgRc). Interaction patterns of ArgRc with L-arginine were also observed for L-canavanine and L-citrulline. L-lysine and L-homoarginine, on the other hand, were shown to bind poorly at the binding site. Figure A general overview of the sites found from docking the various ligands into ArgRc ( grey ribbons). Red coloured sticks: residues in binding site H that was selected for docking
  10. Oumelaz F, Nemiri O, Boumaza A, Meradji H, Ghemid S, Khenata R, et al.
    J Mol Model, 2023 Mar 31;29(4):124.
    PMID: 37000284 DOI: 10.1007/s00894-023-05497-5
    CONTEXT: In this work, a comprehensive study concerning the physical properties of ternary alloys system (AlP1-xBix) at different concentrations is presented. The obtained results from our first-principle calculations are compared with previously reported studies in the literature and discussed in detail. Our computed results are found in a nice agreement where available with earlier reported results. Electronic band structures at the above-mentioned concentrations are also determined. Likewise, the impact of the varying temperature and pressure on Debye temperature, heat capacity, and entropy is analyzed as well. Furthermore, elastic constants and related elastic moduli results are also computed. Our results show that alloys are stable and found to be in brittle nature. This is the first quantitative study related to ternary alloys (AlP1-xBix) at mentioned concentrations. We soon expect the experimental confirmation of our predictions.

    METHOD: The calculations are performed, at concentrations x=0.0, 0.25, 0.5, 0.75, and 1.0 by using the "full potential (FP) linearized (L) augmented plane wave plus local orbital (APW+lo) method framed within density functional theory (DFT)" as recognized in the "WIEN2k computational code". The "quasi-harmonic Debye model" approach is employed to determine the thermal properties of the title alloys.

  11. Mohamed EAR, Abdel-Rahman IM, Zaki MEA, Al-Khdhairawi A, Abdelhamid MM, Alqaisi AM, et al.
    J Mol Model, 2023 Feb 20;29(3):70.
    PMID: 36808314 DOI: 10.1007/s00894-023-05457-z
    BACKGROUND: In November 2021, variant B.1.1.529 of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified by the World Health Organization (WHO) and designated Omicron. Omicron is characterized by a high number of mutations, thirty-two in total, making it more transmissible than the original virus. More than half of those mutations were found in the receptor-binding domain (RBD) that directly interacts with human angiotensin-converting enzyme 2 (ACE2). This study aimed to discover potent drugs against Omicron, which were previously repurposed for coronavirus disease 2019 (COVID-19). All repurposed anti-COVID-19 drugs were compiled from previous studies and tested against the RBD of SARS-CoV-2 Omicron.

    METHODS: As a preliminary step, a molecular docking study was performed to investigate the potency of seventy-one compounds from four classes of inhibitors. The molecular characteristics of the best-performing five compounds were predicted by estimating the drug-likeness and drug score. Molecular dynamics simulations (MD) over 100 ns were performed to inspect the relative stability of the best compound within the Omicron receptor-binding site.

    RESULTS: The current findings point out the crucial roles of Q493R, G496S, Q498R, N501Y, and Y505H in the RBD region of SARS-CoV-2 Omicron. Raltegravir, hesperidin, pyronaridine, and difloxacin achieved the highest drug scores compared with the other compounds in the four classes, with values of 81%, 57%, 18%, and 71%, respectively. The calculated results showed that raltegravir and hesperidin had high binding affinities and stabilities to Omicron with ΔGbinding of - 75.7304 ± 0.98324 and - 42.693536 ± 0.979056 kJ/mol, respectively. Further clinical studies should be performed for the two best compounds from this study.

  12. Al-Qattan MNM, Mordi MN
    J Mol Model, 2023 Aug 16;29(9):281.
    PMID: 37584781 DOI: 10.1007/s00894-023-05650-0
    CONTEXT: Modulation of disease progression is frequently started by identifying biochemical pathway catalyzed by biomolecule that is prone to inhibition by small molecular weight ligands. Such ligands (leads) can be obtained from natural resources or synthetic libraries. However, de novo design based on fragments assembly and optimization is showing increasing success. Plasmodium falciparum parasite depends on glutathione-S-transferase (PfGST) in buffering oxidative heme as an approach to resist some antimalarials. Therefore, PfGST is considered an attractive target for drug development. In this research, fragment-based approaches were used to design molecules that can fit to glutathione (GSH) binding site (G-site) of PfGST.

    METHODS: The involved approaches build molecules from fragments that are either isosteric to GSH sub-moieties (ligand-based) or successfully docked to GSH binding sub-pockets (structure-based). Compared to reference GST inhibitor of S-hexyl GSH, ligands with improved rigidity, synthetic accessibility, and affinity to receptor were successfully designed. The method involves joining fragments to create ligands. The ligands were then explored using molecular docking, Cartesian coordinate's optimization, and simplified free energy determination as well as MD simulation and MMPBSA calculations. Several tools were used which include OPENEYE toolkit, Open Babel, Autodock Vina, Gromacs, and SwissParam server, and molecular mechanics force field of MMFF94 for optimization and CHARMM27 for MD simulation. In addition, in-house scripts written in Matlab were used to control fragments connection and automation of the tools.

  13. Zhou C, Yan L, Xu J, Hamezah HS, Wang T, Du F, et al.
    J Mol Model, 2024 Feb 13;30(3):68.
    PMID: 38347278 DOI: 10.1007/s00894-024-05875-7
    CONTEXT: Adipose triglyceride lipase (ATGL), a key enzyme responsible for lipolysis, catalyzes the first step of lipolysis and converts triglycerides to diacylglycerols and free fatty acids (FFA). Our previous work suggested that phillyrin treatment improves insulin resistance in HFD-fed mice, which was associated with ATGL inhibition. In this study, using docking simulation, we explored the binding pose of phillyrin and atglistatin (a mouse ATGL inhibitor) to ATGL in mouse. From the docking results, the interactions with Ser47 and Asp166 were speculated to have caused phillyrin to inhibit ATGL in mice. Further, molecular dynamics simulation of 100 ns and MM-GBSA were conducted for the protein-ligand complex, which indicated that the system was stable and that phillyrin displayed a better affinity to ATGL than did atglistatin throughout the simulation period. Moreover, the results of pharmacological validation were consistent with those of the in silico simulations. In summary, our study illustrates the potential of molecular docking to accurately predict the binding protein produced by AlphaFold and suggests that phillyrin is a potential small molecule that targets and inhibits ATGL enzymatic activity.

    METHODS: The ATGL-predicted protein structure, verified by PROCHECK, was determined using AlphaFold. Molecular docking, molecular dynamics simulation, and prime molecular mechanic-generalized born surface area were performed using LigPrep, Desmond, and prime MM-GBSA modules of Schrödinger software release 2021-2, respectively. For pharmacological validation, immunoblotting was performed to assess ATGL protein expression. The fluorescence intensity and glycerol concentration were quantified to evaluate the efficiency of phillyrin in inhibiting ATGL.

  14. Al-Maari MA, Hizaddin HF, Salleh MZM, Hayyan A
    J Mol Model, 2024 Feb 17;30(3):73.
    PMID: 38368310 DOI: 10.1007/s00894-024-05876-6
    CONTEXT: One of the prevalent methods for evaluating separation performance is to predict the interactions of solvent and solute molecules. The infinite dilution activity coefficient, Gibbs free energy, and sigma profiles provided insights into the solubilization of a solute and revealed the intensity of the solution's molecular interactions. The effective thermodynamic tools (infinite dilution activity coefficient, Gibbs free energy) were evaluated for predicting the efficiency of 18 polar and non-polar organic solvents in rubber seed oil (RSO) extraction. An infinite dilution activity coefficient was computed to evaluate the solubility of the rubber seed oil model compound (linoleic acid) in the organic solvents. Gibbs free energy was evaluated to show the energy change associated with the molecules mixing process and forecast the miscibility of linoleic acid molecules in the solvents. Moreover, the study examined the sigma profiles and sigma surfaces of organic solvents and linoleic acid to acquire a deeper insight into their similarities and how they interact molecularly. According to the computational prediction and experimental verification, the thermodynamic properties of Gibbs free energy and activity coefficient proved to be highly effective tools for screening polar and moderately polar solvents, predicting the molecular interactions with solute. Whereas the sigma profile and sigma surface were found to be the most efficient tools for evaluating the efficacy of non-polar solvents. Solvents with moderate polarity, such as tetrahydrofuran and diethyl ether, as well as non-polar solvents like pentane, heptane, and n-hexane, proved to be effective and favorable for oil extraction, resulting in the highest oil yields of approximately 27.0%. Overall, the COSMO-RS method demonstrates its utility in estimating the solubility of RSO in organic solvents, enabling early identification of the most effective solvent.

    METHODS: The initial geometry optimization of every component was conducted through density functional theory (DFT) using TmoleX software. A single-point density functional theory (DFT) computation using Becke Perdew 86 (BP86) and the Triple-Zeta Valence Potential (TZVPD) was performed to produce.cosmo files. COSMO-RS calculations were performed by applying the parameterization file BP_TZVPD_FINE_19.ctd using COSMOthermX software. The practical extraction of oil from plant seeds was performed using a sonicator bath to verify the accuracy of the COSMO-RS predictions.

  15. Maulana AF, Maksum IP, Sriwidodo S, Rukayadi Y
    J Mol Model, 2024 Apr 18;30(5):136.
    PMID: 38634946 DOI: 10.1007/s00894-024-05934-z
    CONTEXT: Further understanding of the molecular mechanisms is necessary since it is important for designing new drugs. This study aimed to understand the molecular mechanisms involved in the design of drugs that are inhibitors of the α-glucosidase enzyme. This research aims to gain further understanding of the molecular mechanisms underlying antidiabetic drug design. The molecular docking process yielded 4 compounds with the best affinity energy, including γ-Mangostin, 1,6-dimethyl-ester-3-isomangostin, 1,3,6-trimethyl-ester-α-mangostin, and 3,6,7-trimethyl-ester-γ-mangostin. Free energy calculation with molecular mechanics with generalized born and surface area solvation indicated that the 3,6,7-trimethyl-γ-mangostin had a better free energy value compared to acarbose and simulated maltose together with 3,6,7-trimethyl-γ-mangostin compound. Based on the analysis of electrostatic, van der Waals, and intermolecular hydrogen interactions, 3,6,7-trimethyl-γ-mangostin adopts a noncompetitive inhibition mechanism, whereas acarbose adopts a competitive inhibition mechanism. Consequently, 3,6,7-trimethyl-ester-γ-mangostin, which is a derivative of γ-mangostin, can provide better activity in silico with molecular docking approaches and molecular dynamics simulations.

    METHOD: This research commenced with retrieving protein structures from the RCSB database, generating the formation of ligands using the ChemDraw Professional software, conducting molecular docking with the Autodock Vina software, and performing molecular dynamics simulations using the Amber software, along with the evaluation of RMSD values and intermolecular hydrogen bonds. Free energy, electrostatic interactions, and Van der Waals interaction were calculated using MM/GBSA. Acarbose, used as a positive control, and maltose are simulated together with test compound that has the best free energy. The forcefields used for molecular dynamics simulations are ff19SB, gaff2, and tip3p.

  16. Maksum IP, Rustaman R, Deawati Y, Rukayadi Y, Utami AR, Nafisa ZK
    J Mol Model, 2024 Jul 09;30(8):260.
    PMID: 38981921 DOI: 10.1007/s00894-024-06060-6
    CONTEXT: Diabetes mellitus (DM) is a metabolic disorder disease that causes hyperglycemia conditions and associated with various chronic complications leading to mortality. Due to high toxicity of conventional diabetic drugs, the exploration of natural compounds as alternative diabetes treatments has been widely carried out. Previous in silico studies have highlighted berberine, a natural compound, as a promising alternative in antidiabetic therapy, potentially acting through various pathways, including the inhibition of the FOXO1 transcription factor in the gluconeogenesis pathway. However, the specific mechanism by which berberine interacts with FOXO1 remains unclear, and research in this area is relatively limited. Therefore, this study aims to determine the stability of berberine structure with FOXO1 based on RMSD, RMSF, binding energy, and trajectory analysis to determine the potential of berberine to inhibit the gluconeogenesis pathway. This research was conducted by in silico method with molecular docking using AutoDock4.2 and molecular dynamics study using Amber20, then visualized by VMD.

    METHODS: Docking between ligand and FOXO1 receptor was carried out with Autodock4.2. For molecular dynamics simulations, the force fields of DNA.OL15, protein.ff14SB, gaff2, and tip3p were used.

  17. Ng S, Xie Y, Schaefer HF
    J Mol Model, 2024 Sep 14;30(10):333.
    PMID: 39276285 DOI: 10.1007/s00894-024-06135-4
    CONTEXT: 2,2,2-Trifluoroethanol (TFE) is known as a membrane mimetic solvent. The IR spectrum, 1H NMR spectrum, 13C NMR spin‒lattice relaxation times (T1), and nuclear Overhauser effect (NOE) data are consistent with extensive hydrogen bonding in TFE, but do not lead to structural features of the hydrogen bonding. Hence, DFT computations were carried out. The results predict the existence of a set of H-bonded dimers and trimers. The bond lengths and dihedral angles in these complexes are obtained, together with their dissociation energies. Computations were also performed for the geometry of the two conformers of the isolated monomer. The structure of one of the dimers consists of a 7-member cyclic fragment with a free CF3CH2 side chain. One set of the trimer structures involves the OH of a third monomer H-bonding to one of the F atoms in the CF3 group of the side chain of this dimer, thereby creating three trimer isomers. A fourth trimer cluster is formed from three monomers in which three OH∙∙∙O bonds create a cyclic fragment with three CF3CH2 side chains. The high dissociation energy (with respect to three monomers) indicates the high stability of the trimer complexes. The structural features of the trimer complexes resemble the structure of a conventional liquid crystal molecule and are postulated to resemble the latter in properties and function in solution, but at a much shorter timescale because of the noncovalent bonding. This hydrogen bonding phenomenon of TFE may be related to its function as a membrane memetic solvent.

    METHODS: Initially, IR and NMR spectroscopic methods were used. Standard procedures were followed. For the computations, a hybrid DFT method with empirical dispersion, ωB97X-D, was used. The basis set, 6-311++G**, is of triple-ζ quality, in which polarization functions and diffuse functions were added for all atoms.

  18. Yau MQ, Liew CWY, Toh JH, Loo JSE
    J Mol Model, 2024 Oct 31;30(11):390.
    PMID: 39480515 DOI: 10.1007/s00894-024-06189-4
    CONTEXT: The substantial increase in the number of active and inactive-state CB1 receptor experimental structures has provided opportunities for CB1 drug discovery using various structure-based drug design methods, including the popular end-point methods for predicting binding free energies-Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA). In this study, we have therefore evaluated the performance of MM/PBSA and MM/GBSA in calculating binding free energies for CB1 receptor. Additionally, with both MM/PBSA and MM/GBSA being known for their highly individualized performance, we have evaluated the effects of various simulation parameters including the use of energy minimized structures, choice of solute dielectric constant, inclusion of entropy, and the effects of the five GB models. Generally, MM/GBSA provided higher correlations than MM/PBSA (rMM/GBSA = 0.433 - 0.652 vs. rMM/PBSA = 0.100 - 0.486) regardless of the simulation parameters, while also offering faster calculations. Improved correlations were observed with the use of molecular dynamics ensembles compared with energy minimized structures and larger solute dielectric constants. Incorporation of entropic terms led to unfavorable results for both MM/PBSA and MM/GBSA for a majority of the dataset, while the evaluation of the various GB models exerted a varying effect on both the datasets. The findings obtained in this study demonstrate the utility of MM/PBSA and MM/GBSA in predicting binding free energies for the CB1 receptor, hence providing a useful benchmark for their applicability in the endocannabinoid system as well as other G protein-coupled receptors.

    METHODS: The study utilized the docked dataset (Induced Fit Docking with Glide XP scoring function) from Loo et al., consisting of 46 ligands-23 agonists and 23 antagonists. The equilibrated structures from Loo et al. were subjected to 30 ns production simulations using GROMACS 2018 at 300 K and 1 atm with the velocity rescaling thermostat and the Parinello-Rahman barostat. AMBER ff99SB*-ILDN was used for the proteins, General Amber Force Field (GAFF) was used for the ligands, and Slipids parameters were used for lipids. MM/PBSA and MM/GBSA binding free energies were then calculated using gmx_MMPBSA. The solute dielectric constant was varied between 1, 2, and 4 to study the effect of different solute dielectric constants on the performance of MM/PB(GB)SA. The effect of entropy on MM/PB(GB)SA binding free energies was evaluated using the interaction entropy module implemented in gmx_MMPBSA. Five GB models, GBHCT, GBOBC1, GBOBC2, GBNeck, and GBNeck2, were evaluated to study the effect of the choice of GB models in the performance of MM/GBSA. Pearson correlation coefficients were used to measure the correlation between experimental and predicted binding free energies.

  19. Borikhonov B, Berdimurodov E, Kholikov T, Nik WBW, Katin KP, Demir M, et al.
    J Mol Model, 2024 Oct 02;30(11):359.
    PMID: 39356293 DOI: 10.1007/s00894-024-06157-y
    CONTEXT: This study addresses the development of sustainable pyridinium ionic liquids (ILs) because of their potential applications in agriculture and pharmaceuticals. Pyridinium-based ILs are known for their low melting points, high thermal stability, and moderate solvation properties. We synthesized three novel pyridinium-based ILs: 1-(2-(isopentyloxy)-2-oxoethyl)pyridin-1-ium chloride, 1-(2-(hexyloxy)-2-oxoethyl)pyridin-1-ium chloride, and 1-(2-(benzyloxy)-2-oxoethyl)pyridin-1-ium chloride. The biological activities of these compounds were evaluated through plant growth promotion, herbicidal, and insecticidal assays. Our results show that the benzyloxy derivative significantly enhances wheat and cucumber growth, whereas the isopentyloxy compound has potent herbicidal effects. Computational methods, including DFT calculations and molecular docking, were applied to understand the structure‒activity relationships (SARs) and mechanisms of action.

    METHODS: The computational techniques involved dispersion-corrected density functional theory (DFT) with the B3LYP functional and the 6-311G** basis set. Grimme's D3 corrections were included to account for dispersion interactions. The calculations were performed via GAMESS-US software. Quantum descriptors of reactivity, such as ionization potential, electron affinity, chemical potential, and electrophilicity index, were derived from the HOMO and LUMO energies. Molecular docking studies were conducted via the CB-Dock server via AutoDock Vina software to predict binding affinities to cancer-related proteins. Petra/Osiris/Molinspiration (POM) analysis was used to predict the drug likeness and other pharmaceutical properties of the synthesized ILs.

Related Terms
Filters
Contact Us

Please provide feedback to Administrator ([email protected])

External Links