Displaying publications 1 - 20 of 35 in total

Abstract:
Sort:
  1. Rosli AN, Ahmad MR, Alias Y, Zain SM, Lee VS, Woi PM
    J Mol Model, 2014 Dec;20(12):2533.
    PMID: 25433601 DOI: 10.1007/s00894-014-2533-9
    Design of neutral receptor molecules (ionophores) for beryllium(II) using unsaturated carbonitrile models has been carried out via density functional theory, G3, and G4 calculations. The first part of this work focuses on gas phase binding energies between beryllium(II) and 2-cyano butadiene (2-CN BD), 3-cyano propene (3-CN P), and simpler models with two separate fragments; acrylonitrile and ethylene. Interactions between beryllium(II) and cyano nitrogen and terminal olefin in the models have been examined in terms of geometrical changes, distribution of charge over the entire π-system, and rehybridization of vinyl carbon orbitals. NMR shieldings and vibrational frequencies probed charge centers and strength of interactions. The six-membered cyclic complexes have planar structures with the rehybridized carbon slightly out of plane (16° in 2-CN BD). G3 results show that in 2-CN BD complex participation of vinyl carbon further stabilizes the cyclic adduct by 16.3 kcal mol(-1), whereas, in simpler models, interaction between beryllium(II) and acetonitrile is favorable by 46.4 kcal mol(-1) compared with that of ethylene. The terminal vinyl carbon in 2-CN BD rehybridizes to sp (3) with an increase of 7 % of s character to allow interaction with beryllium(II). G4 calculations show that the Be(II) and 2-CN BD complex is more strongly bound than those with Mg(II) and Ca(II) by 98.5 and 139.2 kcal mol(-1) (-1), respectively. QST2 method shows that the cyclic and acyclic forms of Be(II)-2-CN BD complexes are separated by 12.3 kcal mol(-1) barrier height. Overlap population analysis reveals that Ca(II) can be discriminated based on its tendency to form ionic interaction with the receptor models.
  2. Rosli AN, Abu Bakar MA, Lee VS, Zain SM, Ahmad MR, Abdul Manan NS, et al.
    J Mol Model, 2014 Sep;20(9):2428.
    PMID: 25149440 DOI: 10.1007/s00894-014-2428-9
    In this work, hybrid functional and G4 methods were employed in the rational design of carbonitrile-carboxaldehyde receptor models for cation recognition. Electron-sharing and ionic interactions between the models and the cations were analyzed utilizing the concepts of overlap population, atomic valence, electrostatic potential, and CHELPG charge in order to elucidate the nature of the heteroatom-metal interaction, the N versus O disparity, and the effect of pH. Receptor fragment models from ionomycin were employed to rationalize the selection of receptor models for discriminating group I cations and enhancing the selectivity for Mg(II) rather than Ca(II), and to examine the effects of keto-enol forms and negatively charged sites. The changes in geometries, overlap population, metal valence, and CHELPG charge upon solvation in heptane medium as compared to the gas phase were negligible. The optimized geometries reveal that the interaction between group II cations and the keto, enol, and enolate forms of 2-cyanoethanal causes 12 % bending of the C-C-N angle from linearity. Overlap populations show that the electron-sharing interaction favors group II cations but that the same mechanism allows Li(I) to compete. The total spin of Li(I) is 17 % greater than that of Ca(II), but the G4 binding energies of the two are separated by more than 50 kcal/mol, favoring group II cations, which may eliminate interference from Li(I). 1,2-Dicyanoethylene, which has only one form, shows similar characteristics. CHELPG analysis shows that Mg(II) transfers 25 and 18 % of its positive charge to 2-cyanoethanal enolate and 1,2-dicyanoethylene, respectively. Hydrogen atoms receive most of the positive charge in both receptors, but the N-termini exhibit strikingly different characteristics. Electrostatic potential contour profiles were found to be in good agreement with the atomic charge distributions. The application of uncharged 1,3-dicarbonyl and 2-cyanocarbonyl receptors and a judicious choice of polymeric membrane that suppresses the Hofmeister effect should lead to high selectivity for magnesium, whereas the utilization of multiple negatively charged ionophores should result in selectivity for calcium.
  3. 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.

  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. 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.

  6. Tejo BA, Salleh AB, Pleiss J
    J Mol Model, 2004 Dec;10(5-6):358-66.
    PMID: 15597204
    The effect of organic solvent on the structure and dynamics of proteins was investigated by multiple molecular dynamics simulations (1 ns each) of Candida rugosa lipase in water and in carbon tetrachloride. The choice of solvent had only a minor structural effect. For both solvents the open and the closed conformation of the lipase were near to their experimental X-ray structures (C(alpha) rms deviation 1-1.3 A). However, the solvents had a highly specific effect on the flexibility of solvent-exposed side chains: polar side chains were more flexible in water, but less flexible in organic solvent. In contrast, hydrophobic residues were more flexible in organic solvent, but less flexible in water. As a major effect solvent changed the dynamics of the lid, a mobile element involved in activation of the lipase, which fluctuated as a rigid body about its average position. While in water the deviations were about 1.6 A, organic solvent reduced flexibility to 0.9 A. This increase rigidity was caused by two salt bridges (Lys85-Asp284, Lys75-Asp79) and a stable hydrogen bond (Lys75-Asn 292) in organic solvent. Thus, organic solvents stabilize the lid but render the side chains in the hydrophobic substrate-binding site more mobile. [figure: see text]. Superimposition of open (black, PDB entry 1CRL) and closed (gray, PDB entry 1TRH) conformers of C. rugosa lipase. The mobile lid is indicated.
  7. Al-qattan MN, Mordi MN
    J Mol Model, 2010 May;16(5):1047-58.
    PMID: 19911202 DOI: 10.1007/s00894-009-0618-7
    A molecular docking tool of AutoDock3.05 was evaluated for its ability to reproduce experimentally determined affinities of various sialic acid analogues toward hemagglutinin of influenza A virus. With the exception of those with a C6-modified glycerol side chain, the experimental binding affinities of most sialic acid analogues (C2, C4 and C5-substituted) determined by viral hemadsorption inhibition assay, hemagglutination inhibition assay and nuclear magnetic resonance correlated well with the computationally estimated free energy of binding. Sialic acid analogues with modified glycerol side chains showed only poor correlation between the experimentally determined hemagglutinin inhibitor affinities and AutoDock3.05 scores, suggesting high mobility of the glutamic acid side chain at the glycerol binding pocket, which is difficult to simulate using a flexi-rigid molecular docking approach. In conclusion, except for some glycerol-substituted sialic acid analogues, the results showed the effectiveness of AutoDock3.05 searching and scoring functions in estimating affinities of sialic acid analogues toward influenza A hemagglutinin, making it a reliable tool for screening a database of virtually designed sialic acid analogues for hemagglutinin inhibitors.
  8. Al-qattan MN, Mordi MN
    J Mol Model, 2010 May;16(5):975-91.
    PMID: 19856192 DOI: 10.1007/s00894-009-0606-y
    In this study fragment-based drug design is combined with molecular docking simulation technique, to design databases of virtual sialic acid (SA) analogues with new substitutions at C2, C5 and C6 positions of SA scaffold. Using spaces occupied by C2, C5 and C6 natural moieties of SA when bound to hemagglutinin (HA) crystallographic structure, new fragments that are commercially available were docked independently in all the pockets. The oriented fragments were then connected to the SA scaffold with or without incorporation of linker molecules. The completed analogues were docked to the whole SA binding site to estimate their binding conformations and affinities, generating three databases of HA-bound SA analogues. Selected new analogues showed higher estimated affinities than the natural SA when tested against H3N2, H5N1 and H1N1 subtypes of influenza A. An improvement in the binding energies indicates that fragment-based drug design when combined with molecular docking simulation is capable to produce virtual analogues that can become lead compound candidates for anti-flu drug discovery program.
  9. 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.

  10. 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
  11. 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.
  12. Mosapour Kotena Z, Behjatmanesh-Ardakani R, Hashim R, Manickam Achari V
    J Mol Model, 2013 Feb;19(2):589-99.
    PMID: 22972691 DOI: 10.1007/s00894-012-1576-z
    Density functional theory calculations on two glycosides, namely, n-octyl-β-D-glucopyranoside (C(8)O-β-Glc) and n-octyl-β-D-galactopyranoside (C(8)O-β-Gal) were performed for geometry optimization at the B3LYP/6-31G level. Both molecules are stereoisomers (epimers) differing only in the orientation of the hydroxyl group at the C4 position. Thus it is interesting to investigate electronically the effect of the direction (axial/equatorial) of the hydroxyl group at the C4 position. The structure parameters of X-H∙∙∙Y intramolecular hydrogen bonds were analyzed, while the nature of these bonds and the intramolecular interactions were considered using the atoms in molecules (AIM) approach. Natural bond orbital analysis (NBO) was used to determine bond orders, charge and lone pair electrons on each atom and effective non-bonding interactions. We have also reported electronic energy and dipole moment in gas and solution phases. Further, the electronic properties such as the highest occupied molecular orbital, lowest unoccupied molecular orbital, ionization energy, electron affinity, electronic chemical potential, chemical hardness, softness and electrophilicity index, are also presented here for both C(8)O-β-Glc and C(8)O-β-Gal. These results show that, while C(8)O-β-Glc possess- only one hydrogen bond, C(8)O-β-Gal has two intramolecular hydrogen bonds, which further confirms the anomalous stability of the latter in self-assembly phenomena.
  13. 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.
  14. Behjatmanesh-Ardakani R, Pourroustaei-Ardakani F, Taghdiri M, Kotena ZM
    J Mol Model, 2016 07;22(7):149.
    PMID: 27271162 DOI: 10.1007/s00894-016-3012-2
    This report present the results of natural energy decomposition analysis (NEDA), natural bond orbital (NBO), and quantum theory of atoms in molecules (QTAIM) calculations of three derivatives of biphenyl-1-aza-18-crown-6 ether and their 1:1 complexes with Cd(2+). All calculations used the B3LYP density functional theory in combination with the 6-311G and WTBS basis sets for ligands and Cd(2+) ion, respectively. Ligands 1 and 3 have a single 1-aza-18-crown-6, substituent; ligand 2 has two such substituents. The results show that, in the optimized geometries of the complexes, the distance between N and Cd(2+) is greater than the distance between O and Cd(2+). NBO and QTAIM data confirm these results. There was no stabilization energy or bond critical point for N · · · Cd(2+) in NBO or QTAIM, respectively. Data show that the O · · · Cd(2+) interaction is a kind of closed shell interaction. The trend of the calculated stabilization energy was similar to the experimental data. Different contributions of interaction energies for complex formation were analyzed by NEDA, and the results show that the main component of the interactions is accounted for by polarization.
  15. Parvizpour S, Razmara J, Jomah AF, Shamsir MS, Illias RM
    J Mol Model, 2015 Mar;21(3):63.
    PMID: 25721655 DOI: 10.1007/s00894-015-2617-1
    Here, we present a novel psychrophilic β-glucanase from Glaciozyma antarctica PI12 yeast that has been structurally modeled and analyzed in detail. To our knowledge, this is the first attempt to model a psychrophilic laminarinase from yeast. Because of the low sequence identity (<40%), a threading method was applied to predict a 3D structure of the enzyme using the MODELLER9v12 program. The results of a comparative study using other mesophilic, thermophilic, and hyperthermophilic laminarinases indicated several amino acid substitutions on the surface of psychrophilic laminarinase that totally increased the flexibility of its structure for efficient catalytic reactions at low temperatures. Whereas several structural factors in the overall structure can explain the weak thermal stability, this research suggests that the psychrophilic adaptation and catalytic activity at low temperatures were achieved through existence of longer loops and shorter or broken helices and strands, an increase in the number of aromatic and hydrophobic residues, a reduction in the number of hydrogen bonds and salt bridges, a higher total solvent accessible surface area, and an increase in the exposure of the hydrophobic side chains to the solvent. The results of comparative molecular dynamics simulation and principal component analysis confirmed the above strategies adopted by psychrophilic laminarinase to increase its catalytic efficiency and structural flexibility to be active at cold temperature.
  16. Ramli AN, Azhar MA, Shamsir MS, Rabu A, Murad AM, Mahadi NM, et al.
    J Mol Model, 2013 Aug;19(8):3369-83.
    PMID: 23686283 DOI: 10.1007/s00894-013-1861-5
    A novel α-amylase was isolated successfully from Glaciozyma antarctica PI12 using DNA walking and reverse transcription-polymerase chain reaction (RT-PCR) methods. The structure of this psychrophilic α-amylase (AmyPI12) from G. antarctica PI12 has yet to be studied in detail. A 3D model of AmyPI12 was built using a homology modelling approach to search for a suitable template and to generate an optimum target-template alignment, followed by model building using MODELLER9.9. Analysis of the AmyPI12 model revealed the presence of binding sites for a conserved calcium ion (CaI), non-conserved calcium ions (CaII and CaIII) and a sodium ion (Na). Compared with its template-the thermostable α-amylase from Bacillus stearothermophilus (BSTA)-the binding of CaII, CaIII and Na ions in AmyPI12 was observed to be looser, which suggests that the low stability of AmyPI12 allows the protein to work at different temperature scales. The AmyPI12 amino acid sequence and model were compared with thermophilic α-amylases from Bacillus species that provided the highest structural similarities with AmyPI12. These comparative studies will enable identification of possible determinants of cold adaptation.
  17. 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.

  18. 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.

  19. Ahmadi S, Manickam Achari V, Nguan H, Hashim R
    J Mol Model, 2014 Mar;20(3):2165.
    PMID: 24623320 DOI: 10.1007/s00894-014-2165-0
    Fully atomistic molecular dynamics simulation studies of thermotropic bilayers were performed using a set of glycosides namely n-octyl-β-D-glucopyranoside (β-C8Glc), n-octyl-α-D-glucopyranoside (α-C8Glc), n-octyl-β-D-galactopyranoside (β-C8Gal), and n-octyl-α-D-galactopyranoside (α-C8Gal) to investigate the stereochemical relationship of the epimeric/anomeric quartet liner glycolipids with the same octyl chain group. The results showed that, the anomeric stereochemistry or the axial/equatorial orientation of C1-O1 (α/β) is an important factor controlling the area and d-spacing of glycolipid bilayer systems in the thermotropic phase. The head group tilt angle and the chain ordering properties are affected by the anomeric effect. In addition, the L(C) phase of β-C8Gal, is tilting less compared to those in the fluid L(α). The stereochemistry of the C4-epimeric (axial/equatorial) and anomeric (α/β) centers simultaneously influence the inter-molecular hydrogen bond. Thus, the trend in the values of the hydrogen bond for these glycosides is β-C8Gal > α-C8Glc > β-C8Glc > α-C8Gal. The four bilayer systems showed anomalous diffusion behavior with an observed trend for the diffusion coefficients; and this trend is β-C8Gal > β-C8Glc > α-C8Gal > α-C8Glc. The "bent" configuration of the α-anomer results in an increase of the hydrophobic area, chain vibration and chain disorganization. Since thermal energy is dispensed more entropically for the chain region, the overall molecular diffusion decreases.
  20. Nguan H, Ahmadi S, Hashim R
    J Mol Model, 2012 Dec;18(12):5041-50.
    PMID: 22752540 DOI: 10.1007/s00894-012-1497-x
    A theoretical study of a series of five glucose based glycolipid crown ethers and their complexes with Na(+) and K(+) was performed using the density functional theory with B3LYP/6-31 G* to obtain the optimized geometrical structures and electronic properties. The local nucleophilicity of the five molecules was investigated using Fukui function, while the global nucleophilicity was calculated from the ionization potential and electron affinity. The structures and coordination of the complexes were studied to identify the best match of the glycolipid crown ethers with cations. In general, it was found that the oxygen atoms pairs O2 and O3 (or O4 and O6) on the sugar ring are constrained from moving toward the cation, which results in a weaker O-cation coordination strength for the oxygen pair compared to the other oxygen atoms in the crown ether ring. The thermodynamic properties of the binding of the complexes and the exchange reaction in gas phase were evaluated. The cation selectivity pattern among the five molecules was in good agreement with the experiment.
Related Terms
Filters
Contact Us

Please provide feedback to Administrator (afdal@afpm.org.my)

External Links