Displaying publications 1 - 20 of 34 in total

Abstract:
Sort:
  1. 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.
  2. 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.
  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. Chaudhry AR, Ahmed R, Irfan A, Mohamad M, Muhammad S, Ul Haq B, et al.
    J Mol Model, 2016 Oct;22(10):248.
    PMID: 27683259
    Some important optoelectronic properties of naphtho[2,1-b:6,5-b']difuran (DPNDF) and its two derivatives have been limelighted by applying the density functional theory (DFT). Due to their low cost, high stability and earth abundance, the DPNDF and its derivatives are considered as potential organic semiconductor materials for their optoelectronics applications. Highly proficient derivatives are obtained systematically by attaching -CN (cyanide) to DPNDF at various sites. Our calculations indicate that DPNDF has a wide and direct band gap with an energy gap of 3.157 eV. Whereas the band gaps of its derivatives are found to be decreased by 88 meV for derivative "a" and 300 meV for derivative "b" as a consequence of p orbitals present in C and N atoms in derivative structures. The narrowing of the energy gap and density of states for the derivatives of DPNDF in the present investigation suggest that energy gap can be engineered for desirable optoelectronic applications via derivatives designing. Furthermore, their obtained results for optical parameters such as the dielectric and conductivity functions, reflectivity, refractive index, and the extinction coefficients endorses their aptness for optoelectronic applications. Graphical Abstract Real part of dielectric function for derivative "b".
  5. Chaudhry AR, Ahmed R, Irfan A, Muhammad S, Shaari A, Al-Sehemi AG
    J Mol Model, 2014 Dec;20(12):2547.
    PMID: 25503703 DOI: 10.1007/s00894-014-2547-3
    We have investigated computationally the effects of π-conjugation extension on naphtha[2,1-b:6,5-b'] difuran (DPNDF); where we increase the number of fused NDF (central core) and furan rings in the parent molecule. The molecular structures of all analogues have been optimized at the ground (S0) and first excited (S1) states using density functional theory (DFT) and time-dependent density functional theory (TD-DFT), respectively. Then highest occupied molecular orbitals (HOMOs), the lowest unoccupied molecular orbitals (LUMOs), photophysical properties, adiabatic/vertical electron affinities (EAa)/(EAv), adiabatic/vertical ionization potentials (IPa)/(IPv), and hole/electron reorganization energies λh/λe have been investigated. The effect of NDF and furan rings on structural and electro-optical properties has also been studied. Our calculated reorganization energies of 1a, 1b, and 2c reveal them, materials with balanced hole/electron charge transport, whereas 2a and 2b are good hole-transport materials. By increasing the number of furan rings; the photostability was augmented in 2a, 2b, and 2c.
  6. Chaudhry AR, Ahmed R, Irfan A, Shaari A, Isa AR, Muhammad S, et al.
    J Mol Model, 2015 Aug;21(8):199.
    PMID: 26177706 DOI: 10.1007/s00894-015-2743-9
    The present study spotlights the designing of new derivatives of 2,7-bis (4-octylphenyl) naphtho [2,1-b:6,5-b'] difuran (C8-DPNDF) by substituting the alkyl groups (methyl, ethyl, propyl, butyl, pentyl, hexyl, and heptyl groups) at para position. Density functional theory (DFT) and time-dependent density functional theory (TD-DFT) methods are employed to optimize the molecular structures in ground and first excited states, respectively. Several electro-optical properties including hole/electron reorganization energies (λh/λe), electron affinities (EAs), ionization potentials (IPs), molecular electrostatic potentials (MEP), and frontier molecular orbitals (FMOs) have been evaluated. Furthermore their transfer integrals and intrinsic mobilities values have also been calculated. From this study, it is found that hole mobility of octyl containing derivative is raised to 4.69 cm(2) V(-1) s(-1). Moreover with attaching octyl group, hole transfer integral values have also been enhanced in newly designed derivatives. The balanced hole and electron reorganization energies, and improved transfer integrals lead to enhanced mobility in derivatives with octyl group, highlighting them as an efficient hole transfer material. Unlike the other electro-optical properties, the intrinsic hole mobility has increased because of transfer integral values of octyl containing derivative C8-DPNDF due to the dense and close crystal packing of C8-DPNDF. However, photostability of furan-based materials has not changed by increasing length of extended alkyl chain. Thus our present investigation highlights the importance of alkyl auxiliary groups that are often neglected/replaced with simple methyl group to save computation costs. Graphical Abstract The hole and electron reorganization energies of naphtho[2,1-b:6,5-b']difuran derivatives.
  7. 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.
  8. Woi PM, Bakar MA, Rosli AN, Lee VS, Ahmad MR, Zain S, et al.
    J Mol Model, 2014 May;20(5):2219.
    PMID: 24770548 DOI: 10.1007/s00894-014-2219-3
    DFT and G4 results reveal that cations display the following trends in imparting its positive charge to acrylonitrile; H⁺ > Li⁺ > Na⁺ > K⁺ for group I and Be²⁺ > Mg²⁺ > Ca²⁺ for group II. Solvation by water molecules and interaction with cation make the cyano bond more polarized and exhibits ketene-imine character. Bond order in nitrile-cation complexes has been predicted based on the s character of the covalent bond orbitals. Mulliken, CHELPG, and NPA charges are in good agreement in predicting positive charge buildup and GIAO nuclear deshileding on C1. G4 enthalpies show that Mg²⁺ is more strongly bound to acrylonitrile than to acetonitrile by 3 kcal mol⁻¹, and the proton affinity of the former is higher by 0.8 kcal mol⁻¹. G4 enthalpies of reductions support prior experimental observation that metalated conjugated nitriles show enhanced reactivity toward weak nucleophiles to afford Michael addition products.
  9. Lim YY, Lim TS, Choong YS
    J Mol Model, 2019 Sep 05;25(10):301.
    PMID: 31486892 DOI: 10.1007/s00894-019-4192-3
    The sigma-E transcription factor (σETF) can be found in most of the bacteria cells including Bacillus thuringiensis. However, the cellular regulatory mechanisms of these transcription factors in the mass production of δ-endotoxins during sporulation stage are yet to be revealed. In addition, the recognition of DNA towards σETF DNA binding motifs that led to the transcription activities is also being poorly studied. Therefore, this work studied the possible DNA binding motifs of σETF by utilising in silico approaches. The structure of σETF was first built via three different computational methods. A cognate DNA sequence was then docked to the predicted σETF DNA-binding motifs. The binding free energy calculated using molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) for triplicate 50 ns simulation of σETF-DNA complex revealed favourable binding energy of DNA to σETF (average ∆Gbind = -34.57 kcal/mol) mainly driven by non-polar interactions. This study revealed that σETF LYS131, ARG133, PHE138, TRP146, ARG222, LYS225 and ARG226 are most likely the key residues upon the binding and recognition of DNA prior to transcription actives. Since determination of genome-regulating protein which recognises specific DNA sequence is important to discriminate between the proteins preferences for different genes, this study might provide some understanding on the possible σETF-DNA recognition prior to transcription initiated for the δ-endotoxins production.
  10. Yap BK, Buckle MJ, Doughty SW
    J Mol Model, 2012 Aug;18(8):3639-55.
    PMID: 22354276 DOI: 10.1007/s00894-012-1368-5
    5-HT(1A) serotonin and D1 dopamine receptor agonists have been postulated to be able to improve negative and cognitive impairment symptoms of schizophrenia, while partial agonists and antagonists of the D2 and 5-HT(2A) receptors have been reported to be effective in reducing positive symptoms. There is therefore a need for well-defined homology models for the design of more selective antipsychotic agents, since no three-dimensional (3D) crystal structures of these receptors are currently available. In this study, homology models were built based on the high-resolution crystal structure of the β(2)-adrenergic receptor (2RH1) and further refined via molecular dynamics simulations in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer system with the GROMOS96 53A6 united atom force field. Docking evaluations with representative agonists and antagonists using AutoDock 4.2 revealed binding modes in agreement with experimentally determined site-directed mutagenesis data and significant correlations between the computed and experimental pK (i) values. The models are also able to distinguish between antipsychotic agents with different selectivities and binding affinities for the four receptors, as well as to differentiate active compounds from decoys. Hence, these human 5-HT(1A), 5-HT(2A), D1 and D2 receptor homology models are capable of predicting the activities of novel ligands, and can be used as 3D templates for antipsychotic drug design and discovery.
  11. Yahya R, Karjiban RA, Basri M, Rahman MB, Girardi M
    J Mol Model, 2014 Nov;20(11):2512.
    PMID: 25381172 DOI: 10.1007/s00894-014-2512-1
    Nonionic surfactants such as the Brij® series are important in the preparation of transdermal drug nanodelivery products using nanoemulsions because of their low toxicity and low irritancy. Here, Monte Carlo (MC) simulation was used to examine the physical behavior of the model deterministic system by using sampling procedures. Metropolis MC simulations were run on three mixtures of two different nonionic surfactants, Brij92 and Brij96, with different compositions in aqueous solution. The system was simulated in the canonical ensemble with constant temperature, volume and number of molecules. Hence, the acceptance ratio for single atom moves of the mixed surfactants increased as the concentration of surfactants increased from 0.494 to 0.591. The lowest total energy for the mixed surfactant systems was -99,039 kcal mol(-1) due to the interaction between all molecules in the system simulated. The physicochemical properties of models such as the radius of gyration and radial distribution function, were also determined. These observations indicate that the behavior and physicochemical of mixed surfactant and PKOEs nanoemulsion systems were described adequately during the simulation.
  12. Mohamad M, Ahmed R, Shaari A, Goumri-Said S
    J Mol Model, 2015 Feb;21(2):27.
    PMID: 25631921 DOI: 10.1007/s00894-015-2582-8
    Escalating demand for sustainable energy resources, because of the rapid exhaustion of conventional energy resources as well as to maintain the environmental level of carbon dioxide (CO2) to avoid its adverse effect on the climate, has led to the exploitation of photovoltaic technology manifold more than ever. In this regard organic materials have attracted great attention on account of demonstrating their potential to harvest solar energy at an affordable rate for photovoltaic technology. 2-vinyl-4,5-dicyanoimidazole (vinazene) is considered as a suitable material over the fullerenes for photovoltaic applications because of its particular chemical and physical nature. In the present study, DFT approaches are employed to provide an exposition of optoelectronic properties of vinazene molecule and molecular crystal. To gain insight into its properties, different forms of exchange correlation energy functional/potential such as LDA, GGA, BLYP, and BL3YP are used. Calculated electronic structure of vinazene molecule has been displayed via HOMO-LUMO isosurfaces, whereas electronic structure of the vinazene molecular crystal, via electronic band structure, is presented. The calculated electronic and optical properties were analyzed and compared as well. Our results endorse vinazene as a suitable material for organic photovoltaic applications.
  13. 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.
  14. Pang WC, Ramli ANM, Hamid AAA
    J Mol Model, 2020 May 16;26(6):142.
    PMID: 32417971 DOI: 10.1007/s00894-020-04398-1
    Fruit bromelain is a cysteine protease accumulated in pineapple fruits. This proteolytic enzyme has received high demand for industrial and therapeutic applications. In this study, fruit bromelain sequences QIM61759, QIM61760 and QIM61761 were retrieved from the National Center for Biotechnology Information (NCBI) Genbank Database. The tertiary structure of fruit bromelain QIM61759, QIM61760 and QIM61761 was generated by using MODELLER. The result revealed that the local stereochemical quality of the generated models was improved by using multiple templates during modelling process. Moreover, by comparing with the available papain model, structural analysis provides an insight on how pro-peptide functions as a scaffold in fruit bromelain folding and contributing to inactivation of mature protein. The structural analysis also disclosed the similarities and differences between these models. Lastly, thermal stability of fruit bromelain was studied. Molecular dynamics simulation of fruit bromelain structures at several selected temperatures demonstrated how fruit bromelain responds to elevation of temperature.
  15. 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.

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

  20. 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.
Related Terms
Filters
Contact Us

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

External Links