MyMedR

Displaying publications 81 - 100 of 226 in total

Abstract:

Sort:

Encapsulation of caffeine into starch matrices: Bitterness evaluation and suppression mechanism

Shao M, Li S, Tan CP, Kraithong S, Gao Q, Fu X, et al.

Int J Biol Macromol, 2021 Mar 15;173:118-127.
PMID: 33444656 DOI: 10.1016/j.ijbiomac.2021.01.043

In this study, caffeine (CA) was encapsulated into food-grade starch matrices, including swelled starch (SS), porous starch (PS), and V-type starch (VS). The bitterness of the microcapsules and suppression mechanisms were investigated using an electronic tongue, molecular dynamics (MD) simulation and the in vitro release kinetics of CA. All the CA-loaded microcapsules showed a lower bitterness intensity than the control. The MD results proved that the weak interactions between starch and CA resulted in a moderate CA release rate for SS-CA microcapsules. The PS-CA microcapsule presented the longest CA release, up to 40 min, whereas the VS-CA microcapsule completely released CA in 9 min. The CA release rate was found to be related to the microcapsule structure and rehydration properties. A low CA bitterness intensity could be attributed to a delay in the CA release rate and resistance to erosion of the microcapsules. The results of this work are valuable for improving starch-based microcapsules (oral-targeted drug-delivery systems) by suppressing the bitterness of alkaloid compounds.

Matched MeSH terms: Molecular Dynamics Simulation
Molecular Dynamics Simulations of Glycyrrhizic Acid Aggregates as Drug-Carriers for Paclitaxel

Hussain M

Curr Drug Deliv, 2019;16(7):618-627.
PMID: 30868954 DOI: 10.2174/1567201816666190313155117

BACKGROUND: Glycyrrhizic acid (GA) is a glycoside that has shown considerable promise as a penetration enhancer and drug carrier to improve the absorption of poorly water-soluble drugs. The aggregation behavior of GA and its ability to form large micelles at higher solution concentrations are thought to contribute to these bioavailability enhancing properties. The oral absorption of Paclitaxel (PTX) for example, an anti-cancer agent which exhibits poor oral bioavailability, has been found to significantly increase in the presence of GA.
METHODS: In an attempt to visualize the aggregation behavior of GA and its subsequent association with PTX, 100 ns molecular dynamics simulation of a 5 mM aqueous solution of GA with 10 molecules of PTX was conducted using GROMACS and an all-atom forcefield.
RESULTS: Aggregation of GA molecules was found to occur quickly at this level of saturation leading to two stable aggregates of 13 and 17 GA molecules with an effective radius of 10.17 nm to 10.92 nm. These aggregates form not in isolation, but together with PTX molecule embedded within the structures, which reduces the number of interactions and hydrogen-bonding with water.
CONCLUSION: GA aggregation occurs around PTX molecules in solution, forming co-joined GA-PTX cluster units at a ratio of 3:1. These clusters remain stable for the remainder of the 100ns simulation and serve to isolate and protect PTX from the aqueous environment.

Matched MeSH terms: Molecular Dynamics Simulation
Insights into the stereospecificity of the d-specific dehalogenase from Rhizobium sp. RC1 toward d- and l-2-chloropropionate

Sudi IY, Hamid AA, Shamsir MS, Jamaluddin H, Wahab RA, Huyop F

Biotechnology, biotechnological equipment, 2014 Jul 04;28(4):608-615.
PMID: 26740767

Halogenated compounds are recalcitrant environmental pollutants prevalent in agricultural fields, waste waters and industrial by-products, but they can be degraded by dehalogenase-containing microbes. Notably, 2-haloalkanoic acid dehalogenases are employed to resolve optically active chloropropionates, as exemplified by the d-specific dehalogenase from Rhizobium sp. RCI (DehD), which acts on d-2-chloropropionate but not on its l-enantiomer. The catalytic residues of this dehalogenase responsible for its affinity toward d-2-chloropropionate have not been experimentally determined, although its three-dimensional crystal structure has been solved. For this study, we performed in silico docking and molecular dynamic simulations of complexes formed by this dehalogenase and d- or l-2-chloropropionate. Arg134 of the enzyme plays the key role in the stereospecific binding and Arg16 is in a position that would allow it to activate a water molecule for hydrolytic attack on the d-2-chloropropionate chiral carbon for release of the halide ion to yield l-2-hydroxypropionate. We propose that within the DehD active site, the NH group of Arg134 can form a hydrogen bond with the carboxylate of d-2-chloropropionate with a strength of ∼4 kcal/mol that may act as an acid-base catalyst, whereas, when l-2-chloropropionate is present, this bond cannot be formed. The significance of the present work is vital for rational design of this dehalogenase in order to confirm the involvement of Arg16 and Arg134 residues implicated in hydrolysis and binding of d-2-chloropropionate in the active site of d-specific dehalogenase from Rhizobium sp. RC1.

Matched MeSH terms: Molecular Dynamics Simulation
Fulltext Structure prediction, molecular dynamics simulation and docking studies of D-specific dehalogenase from Rhizobium sp. RC1

Sudi IY, Wong EL, Joyce-Tan KH, Shamsir MS, Jamaluddin H, Huyop F

Int J Mol Sci, 2012;13(12):15724-54.
PMID: 23443090 DOI: 10.3390/ijms131215724

Currently, there is no three-dimensional structure of D-specific dehalogenase (DehD) in the protein database. We modeled DehD using ab initio technique, performed molecular dynamics (MD) simulation and docking of D-2-chloropropionate (D-2CP), D-2-bromopropionate (D-2BP), monochloroacetate (MCA), monobromoacetate (MBA), 2,2-dichloropropionate (2,2-DCP), d,l-2,3-dichloropropionate (d,l-2,3-DCP), and 3-chloropropionate (3-CP) into the DehD active site. The sequences of DehD and D-2-haloacid dehalogenase (HadD) from Pseudomonas putida AJ1 have 15% sequence similarity. The model had 80% of the amino acid residues in the most favored region when compared to the crystal structure of DehI from Pseudomonas putida PP3. Docking analysis revealed that Arg107, Arg134 and Tyr135 interacted with D-2CP, and Glu20 activated the water molecule for hydrolytic dehalogenation. Single residue substitutions at 25-30 °C showed that polar residues of DehD were stable when substituted with nonpolar residues and showed a decrease in activity within the same temperature range. The molecular dynamics simulation of DehD and its variants showed that in R134A variant, Arg107 interacted with D-2CP, while in Y135A, Gln221 and Arg231 interacted with D-2CP. It is our emphatic belief that the new model will be useful for the rational design of DehDs with enhanced potentials.

Matched MeSH terms: Molecular Dynamics Simulation*
Fulltext An S188V mutation alters substrate specificity of non-stereospecific α-haloalkanoic acid dehalogenase E (DehE)

Hamid AA, Hamid TH, Wahab RA, Omar MS, Huyop F

PLoS One, 2015;10(3):e0121687.
PMID: 25816329 DOI: 10.1371/journal.pone.0121687

The non-stereospecific α-haloalkanoic acid dehalogenase E (DehE) degrades many halogenated compounds but is ineffective against β-halogenated compounds such as 3-chloropropionic acid (3CP). Using molecular dynamics (MD) simulations and site-directed mutagenesis we show here that introducing the mutation S188V into DehE improves substrate specificity towards 3CP. MD simulations showed that residues W34, F37, and S188 of DehE were crucial for substrate binding. DehE showed strong binding ability for D-2-chloropropionic acid (D-2CP) and L-2-chloropropionic acid (L-2CP) but less affinity for 3CP. This reduced affinity was attributed to weak hydrogen bonding between 3CP and residue S188, as the carboxylate of 3CP forms rapidly interconverting hydrogen bonds with the backbone amide and side chain hydroxyl group of S188. By replacing S188 with a valine residue, we reduced the inter-molecular distance and stabilised bonding of the carboxylate of 3CP to hydrogens of the substrate-binding residues. Therefore, the S188V can act on 3CP, although its affinity is less strong than for D-2CP and L-2CP as assessed by Km. This successful alteration of DehE substrate specificity may promote the application of protein engineering strategies to other dehalogenases, thereby generating valuable tools for future bioremediation technologies.

Matched MeSH terms: Molecular Dynamics Simulation
Theoretical analyses on enantiospecificity of L-2-haloacid dehalogenase (DehL) from Rhizobium sp. RC1 towards 2-chloropropionic acid

Adamu A, Abdul Wahab R, Aliyu F, Abdul Razak FI, Mienda BS, Shamsir MS, et al.

J Mol Graph Model, 2019 11;92:131-139.
PMID: 31352207 DOI: 10.1016/j.jmgm.2019.07.012

Dehalogenases continue to garner interest of the scientific community due to their potential applications in bioremediation of halogen-contaminated environment and in synthesis of various industrially relevant products. Example of such enzymes is DehL, an L-2-haloacid dehalogenase (EC 3.8.1.2) from Rhizobium sp. RC1 that catalyses the specific cleavage of halide ion from L-2-halocarboxylic acids to produce the corresponding D-2-hydroxycarboxylic acids. Recently, the catalytic residues of DehL have been identified and its catalytic mechanism has been fully elucidated. However, the enantiospecificity determinants of the enzyme remain unclear. This information alongside a well-defined catalytic mechanism are required for rational engineering of DehL for substrate enantiospecificity. Therefore, using quantum mechanics/molecular mechanics and molecular mechanics Poisson-Boltzmann surface area calculations, the current study theoretically investigated the molecular basis of DehL enantiospecificity. The study found that R51L mutation cancelled out the dehalogenation activity of DehL towards it natural substrate, L-2-chloropropionate. The M48R mutation, however introduced a new activity towards D-2-chloropropionate, conveying the possibility of inverting the enantiospecificity of DehL from L-to d-enantiomer with a minimum of two simultaneous mutations. The findings presented here will play important role in the rational design of DehL dehalogenase for improving substrate utility.

Matched MeSH terms: Molecular Dynamics Simulation
Deciphering the catalytic amino acid residues of l-2-haloacid dehalogenase (DehL) from Rhizobium sp. RC1: An in silico analysis

Adamu A, Wahab RA, Shamsir MS, Aliyu F, Huyop F

Comput Biol Chem, 2017 Oct;70:125-132.
PMID: 28873365 DOI: 10.1016/j.compbiolchem.2017.08.007

The l-2-haloacid dehalogenases (EC 3.8.1.2) specifically cleave carbon-halogen bonds in the L-isomers of halogenated organic acids. These enzymes have potential applications for the bioremediation and synthesis of various industrial products. One such enzyme is DehL, the l-2-haloacid dehalogenase from Rhizobium sp. RC1, which converts the L-isomers of 2-halocarboxylic acids into the corresponding D-hydroxycarboxylic acids. However, its catalytic mechanism has not been delineated, and to enhance its efficiency and utility for environmental and industrial applications, knowledge of its catalytic mechanism, which includes identification of its catalytic residues, is required. Using ab initio fragment molecular orbital calculations, molecular mechanics Poisson-Boltzmann surface area calculations, and classical molecular dynamic simulation of a three-dimensional model of DehL-l-2-chloropropionic acid complex, we predicted the catalytic residues of DehL and propose its catalytic mechanism. We found that when Asp13, Thr17, Met48, Arg51, and His184 were individually replaced with an alanine in silico, a significant decrease in the free energy of binding for the DehL-l-2-chloropropionic acid model complex was seen, indicating the involvement of these residues in catalysis and/or structural integrity of the active site. Furthermore, strong inter-fragment interaction energies calculated for Asp13 and L-2-chloropropionic acid, and for a water molecule and His184, and maintenance of the distances between atoms in the aforementioned pairs during the molecular dynamics run suggest that Asp13 acts as the nucleophile and His184 activates the water involved in DehL catalysis. The results of this study should be important for the rational design of a DehL mutant with improved catalytic efficiency.

Matched MeSH terms: Molecular Dynamics Simulation*
In silico analysis of a putative dehalogenase from the genome of halophilic bacterium Halomonas smyrnensis AAD6T

Oyewusi HA, Akinyede KA, Abdul Wahab R, Huyop F

J Biomol Struct Dyn, 2023 Jan;41(1):319-335.
PMID: 34854349 DOI: 10.1080/07391102.2021.2006085

Microbial-assisted removal of natural or synthetic pollutants is the prevailing green, low-cost technology to treat polluted environments. However, the challenge with enzyme-assisted bioremediation is the laborious nature of dehalogenase-producing microorganisms' bioprospecting. This bottleneck could be circumvented by in-silico analysis of certain microorganisms' whole-genome sequences to predict their protein functions and enzyme versatility for improved biotechnological applications. Herein, this study performed structural analysis on a dehalogenase (DehHsAAD6) from the genome of Halomonas smyrnensis AAD6 by molecular docking and molecular dynamic (MD) simulations. Other bioinformatics tools were also employed to identify substrate preference (haloacids and haloacetates) of the DehHsAAD6. The DehHsAAD6 preferentially degraded haloacids and haloacetates (-3.2-4.8 kcal/mol) and which formed three hydrogen bonds with Tyr12, Lys46, and Asp182. MD simulations data revealed the higher stability of DehHsAAD6-haloacid- (RMSD 0.22-0.3 nm) and DehHsAAD6-haloacetates (RMSF 0.05-0.14 nm) complexes, with the DehHsAAD6-L-2CP complex being the most stable. The detail of molecular docking calculations ranked complexes with the lowest binding free energies as: DehHsAAD6-L-2CP complex (-4.8 kcal/mol) = DehHsAAD6-MCA (-4.8 kcal/mol) < DehHsAAD6-TCA (-4.5 kcal/mol) < DehHsAAD6-2,3-DCP (-4.1 kcal/mol) < DehHsAAD6-D-2CP (-3.9 kcal/mol) < DehHsAAD6-2,2-DCP (-3.5 kcal/mol) < DehHsAAD6-3CP (-3.2 kcal/mol). In a nutshell, the study findings offer valuable perceptions into the elucidation of possible reaction mechanisms of dehalogenases for extended substrate specificity and higher catalytic activity.Communicated by Ramaswamy H. Sarma.

Matched MeSH terms: Molecular Dynamics Simulation*
Comparative modeling and enzymatic affinity of novel haloacid dehalogenase from Bacillus megaterium strain BHS1 isolated from alkaline Blue Lake in Turkey

Wahhab BH, Oyewusi HA, Wahab RA, Mohammad Hood MH, Abdul Hamid AA, Al-Nimer MS, et al.

J Biomol Struct Dyn, 2024;42(3):1429-1442.
PMID: 37038649 DOI: 10.1080/07391102.2023.2199870

This study presents the initial structural model of L-haloacid dehalogenase (DehLBHS1) from Bacillus megaterium BHS1, an alkalotolerant bacterium known for its ability to degrade halogenated environmental pollutants. The model provides insights into the structural features of DehLBHS1 and expands our understanding of the enzymatic mechanisms involved in the degradation of these hazardous pollutants. Key amino acid residues (Arg40, Phe59, Asn118, Asn176, and Trp178) in DehLBHS1 were identified to play critical roles in catalysis and molecular recognition of haloalkanoic acid, essential for efficient binding and transformation of haloalkanoic acid molecules. DehLBHS1 was modeled using I-TASSER, yielding a best TM-score of 0.986 and an RMSD of 0.53 Å. Validation of the model using PROCHECK revealed that 89.2% of the residues were located in the most favored region, providing confidence in its structural accuracy. Molecular docking simulations showed that the non-simulated DehLBHS1 preferred 2,2DCP over other substrates, forming one hydrogen bond with Arg40 and exhibiting a minimum energy of -2.5 kJ/mol. The simulated DehLBHS1 exhibited a minimum energy of -4.3 kJ/mol and formed four hydrogen bonds with Arg40, Asn176, Asp9, and Tyr11, further confirming the preference for 2,2DCP. Molecular dynamics simulations supported this preference, based on various metrics, including RMSD, RMSF, gyration, hydrogen bonding, and molecular distance. MM-PBSA calculations showed that the DehLBHS1-2,2-DCP complex had a markedly lower binding energy (-21.363 ± 1.26 kcal/mol) than the DehLBHS1-3CP complex (-14.327 ± 1.738 kcal/mol). This finding has important implications for the substrate specificity and catalytic function of DehLBHS1, particularly in the bioremediation of 2,2-DCP in contaminated alkaline environments. These results provide a detailed view of the molecular interactions between the enzyme and its substrate and may aid in the development of more efficient biocatalytic strategies for the degradation of halogenated compounds.Communicated by Ramaswamy H. Sarma.

Matched MeSH terms: Molecular Dynamics Simulation
Molecular dynamics simulations suggest changes in electrostatic interactions as a potential mechanism through which serine phosphorylation inhibits DNA Polymerase β's activity

Homouz D, Joyce-Tan KH, Shahir Shamsir M, Moustafa IM, Idriss H

J Mol Graph Model, 2018 01;79:192.
PMID: 29223917 DOI: 10.1016/j.jmgm.2017.11.002

DNA polymerase β is a 39kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31kDa domain responsible for the polymerase activity and an 8kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S44/55 and S55 phosphorylation were less dramatic than S44 and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is likely the major contributor to structural fluctuations that lead to loss of enzymatic activity.

Matched MeSH terms: Molecular Dynamics Simulation
Molecular dynamics simulations suggest changes in electrostatic interactions as a potential mechanism through which serine phosphorylation inhibits DNA polymerase β activity

Homouz D, Joyce-Tan KH, ShahirShamsir M, Moustafa IM, Idriss HT

J Mol Graph Model, 2018 09;84:236-241.
PMID: 30138833 DOI: 10.1016/j.jmgm.2018.08.007

DNA polymerase β is a 39 kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31 kDa domain responsible for the polymerase activity and an 8 kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics simulations. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S55/44 and S55 phosphorylation were less dramatic and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is the major contributor to structural fluctuations that lead to loss of enzymatic activity.

Matched MeSH terms: Molecular Dynamics Simulation*
Structural prediction of a novel laminarinase from the psychrophilic Glaciozyma antarctica PI12 and its temperature adaptation analysis

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.

Matched MeSH terms: Molecular Dynamics Simulation
Structural prediction of a novel chitinase from the psychrophilic Glaciozyma antarctica PI12 and an analysis of its structural properties and function

Ramli AN, Mahadi NM, Shamsir MS, Rabu A, Joyce-Tan KH, Murad AM, et al.

J Comput Aided Mol Des, 2012 Aug;26(8):947-61.
PMID: 22710891 DOI: 10.1007/s10822-012-9585-7

The structure of psychrophilic chitinase (CHI II) from Glaciozyma antarctica PI12 has yet to be studied in detail. Due to its low sequence identity (<30 %), the structural prediction of CHI II is a challenge. A 3D model of CHI II was built by first using a threading approach to search for a suitable template and to generate an optimum target-template alignment, followed by model building using MODELLER9v7. Analysis of the catalytic insertion domain structure in CHI II revealed an increase in the number of aromatic residues and longer loops compared to mesophilic and thermophilic chitinases. A molecular dynamics simulation was used to examine the stability of the CHI II structure at 273, 288 and 300 K. Structural analysis of the substrate-binding cleft revealed a few exposed aromatic residues. Substitutions of certain amino acids in the surface and loop regions of CHI II conferred an increased flexibility to the enzyme, allowing for an adaptation to cold temperatures. A substrate binding comparison of CHI II with the mesophilic chitinase from Coccidioides immitis, 1D2K, suggested that the psychrophilic adaptation and catalytic activity at low temperatures were achieved through a reduction in the number of salt bridges, fewer hydrogen bonds and an increase in the exposure of the hydrophobic side chains to the solvent.

Matched MeSH terms: Molecular Dynamics Simulation*
Computational docking, molecular dynamics simulation and subsite structure analysis of a maltogenic amylase from Bacillus lehensis G1 provide insights into substrate and product specificity

Manas NH, Bakar FD, Illias RM

J Mol Graph Model, 2016 06;67:1-13.
PMID: 27155296 DOI: 10.1016/j.jmgm.2016.04.004

Maltogenic amylase (MAG1) from Bacillus lehensis G1 displayed the highest hydrolysis activity on β-cyclodextrin (β-CD) to produce maltose as a main product and exhibited high transglycosylation activity on malto-oligosaccharides with polymerization degree of three and above. These substrate and product specificities of MAG1 were elucidated from structural point of view in this study. A three-dimensional structure of MAG1 was constructed using homology modeling. Docking of β-CD and malto-oligosaccharides was then performed in the MAG1 active site. An aromatic platform in the active site was identified which is responsible in substrate recognition especially in determining the enzyme's preference toward β-CD. Molecular dynamics (MD) simulation showed MAG1 structure is most stable when docked with β-CD and least stable when docked with maltose. The docking analysis and MD simulation showed that the main subsites for substrate stabilization in the active site are -2, -1, +1 and +2. A bulky residue, Trp359 at the +2 subsite was identified to cause steric interference to the bound linear malto-oligosaccharides thus prevented it to occupy subsite +3, which can only be reached by a highly bent glucose molecule such as β-CD. The resulted modes of binding from docking simulation show a good correlation with the experimentally determined hydrolysis pattern. The subsite structure generated from this study led to a possible mode of action that revealed how maltose was mainly produced during hydrolysis. Furthermore, maltose only occupies subsite +1 and +2, therefore could not be hydrolyzed or transglycosylated by the enzyme. This important knowledge has paved the way for a novel structure-based molecular design for modulation of its catalytic activities.

Matched MeSH terms: Molecular Dynamics Simulation*
Synthesis, characterization, monoamine oxidase inhibition, molecular docking and dynamic simulations of novel 2,1-benzothiazine-2,2-dioxide derivatives

Ahmad S, Zaib S, Jalil S, Shafiq M, Ahmad M, Sultan S, et al.

Bioorg Chem, 2018 10;80:498-510.
PMID: 29996111 DOI: 10.1016/j.bioorg.2018.04.012

In this research work, we report the synthesis and biological evaluation of two new series of 1-benzyl-4-(benzylidenehydrazono)-3,4-dihydro-1H-benzo[c] [1,2]thiazine 2,2-dioxides and 1-benzyl-4-((1-phenylethylidene)hydrazono)-3,4-dihydro-1H-benzo[c][1,2]thiazine 2,2-dioxides. The synthetic plan involves the mesylation of methyl anthranilate with subsequent N-benzylation of the product. The methyl 2-(N-benzylmethylsulfonamido)benzoate was subjected to cyclization reaction in the presence of sodium hydride to obtain 1-benzyl-1H-benzo[c][1,2]thiazin-4(3H)-one 2,2-dioxide which was treated with hydrazine hydrate to get corresponding hydrazone precursor. Finally, the titled compounds were obtained by reaction of hydrazone with various substituted aldehydes and ketones. The synthesized derivatives were subjected to carry out their inhibition activities against monoamine oxidases along with modelling investigations to evaluate their binding interactions and dynamic stability during the docking studies. The inhibition profile of potent compounds was found as competitive for both the isozymes. The compounds were more selective inhibitors of MAO-A as compared to MAO-B. Moreover, drug likeness profile of the derivatives was evaluated to have an additional insight into the physicochemical properties. The molecular dynamic simulations predicted the behaviour of amino acids with the active site residues.

Matched MeSH terms: Molecular Dynamics Simulation
Multi-step molecular docking and dynamics simulation-based screening of large antiviral specific chemical libraries for identification of Nipah virus glycoprotein inhibitors

Kalbhor MS, Bhowmick S, Alanazi AM, Patil PC, Islam MA

Biophys Chem, 2021 03;270:106537.
PMID: 33450550 DOI: 10.1016/j.bpc.2020.106537

Nipah virus (NiV) infections are highly contagious and can cause severe febrile encephalitis. An outbreak of NiV infection has reported high mortality rates in Southeast Asian countries including Bangladesh, East Timor, Malaysia, Papua New Guinea, Vietnam, Cambodia, Indonesia, Madagascar, Philippines, Thailand and India. Considering the high risk for an epidemic outbreak, the World Health Organization (WHO) declared NiV as an emerging priority pathogen. However, there are no effective therapeutics or any FDA approved drugs available for the treatment of this infection. Among the known nine proteins of NiV, glycoprotein plays an important role in initiating the entry of viruses and attaching to the host cell receptors. Herein, three antiviral databases consisting of 79,892 chemical entities have been computationally screened against NiV glycoprotein (NiV-G). Particularly, multi-step molecular docking followed by extensive molecular binding interactions analyses, binding free energy estimation, in silico pharmacokinetics, synthetic accessibility and toxicity profile evaluations have been carried out for initial identification of potential NiV-G inhibitors. Further, molecular dynamics (MD) simulation has been performed to understand the dynamic properties of NiV-G protein-bound with proposed five inhibitors (G1-G5) and their interactions behavior, and any conformational changes in NiV-G protein during simulations. Moreover, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) based binding free energies (∆G) has been calculated from all MD simulation trajectories to understand the energy contribution of each proposed compound in maintaining and stabilizing the complex binding interactions with NiV-G protein. Proposed compounds showed high negative ∆G values ranging from -166.246 to -226.652 kJ/mol indicating a strong affinity towards the NiV-G protein.

Matched MeSH terms: Molecular Dynamics Simulation
Design, synthesis, anti-acetylcholinesterase evaluation and molecular modelling studies of novel coumarin-chalcone hybrids

Hasan AH, Shakya S, Hussain FHS, Murugesan S, Chander S, Pratama MRF, et al.

J Biomol Struct Dyn, 2023;41(21):11450-11462.
PMID: 36591704 DOI: 10.1080/07391102.2022.2162583

The major enzyme responsible for the hydrolytic breakdown of the neurotransmitter acetylcholine (ACh) is acetylcholinesterase (AChE). Acetylcholinesterase inhibitors (AChEIs) are the most prescribed class of medications for the treatment of Alzheimer's disease (AD) and dementia. The limitations of available therapy, like side effects, drug tolerance, and inefficacy in halting disease progression, drive the need for better, more efficacious, and safer drugs. In this study, a series of fourteen novel chalcone-coumarin derivatives (8a-n) were designed, synthesized and characterized by spectral techniques like FT-IR, NMR, and HR-MS. Subsequently, the synthesized compounds were tested for their ability to inhibit acetylcholinesterase (AChE) activity by Ellman's method. All tested compounds showed AChE inhibition with IC50 value ranging from 0.201 ± 0.008 to 1.047 ± 0.043 μM. Hybrid 8d having chloro substitution on ring-B of the chalcone scaffold showed relatively better potency, with IC50 value of 0.201 ± 0.008 μM compared to other members of the series. The reference drug, galantamine, exhibited an IC50 at 1.142 ± 0.027 μM. Computational studies revealed that designed compounds bind to the peripheral anionic site (PAS), the catalytic active site (CAS), and the mid-gorge site of AChE. Putative binding modes, ligand-enzyme interactions, and stability of the best active compound are studied using molecular docking, followed by molecular dynamics (MD) simulations. The cytotoxicity of the synthesised derivatives was determined using the MTT test at three concentrations (100 g/mL, 500 g/mL, and 1 mg/mL). None of the chemicals had a significant effect on the body at the highest dose of 1 mg/mL.Communicated by Ramaswamy H. Sarma.

Matched MeSH terms: Molecular Dynamics Simulation
Fulltext New reaction kinetic model derived from molecular dynamics simulation data on non-linearity of rate coefficients and their relation to activity coefficients

Jesudason, C.G.

ASM Science Journal, 2007;1(1):7-18.
MyJurnal

Molecular dynamics reaction simulation showed that the rate constant is not constant over the concentration profile of reactants and products over a fixed temperature regime, and this variation is expressed in terms of the defined reactivity coefficients. The ratio of these coefficients for the forward and backward reactions were found to equal that of the activity coefficient ratio for the product and reactant species. A theory was developed to explain kinetics in general based on these observations. Several other theorems had first to be developed, most striking of all was the inference that the excess Helmholtz free energy was the thermodynamical function which had a direct relation to these activity factors than the Gibbs free energy. The theory is applied to a class of ionic reactions which could not be rationalized using the standard Bjørn-Bjerrum theory of ionic reactions.

Matched MeSH terms: Molecular Dynamics Simulation
Insight into the molecular mechanism of P-glycoprotein mediated drug toxicity induced by bioflavonoids: an integrated computational approach

Wongrattanakamon P, Lee VS, Nimmanpipug P, Sirithunyalug B, Chansakaow S, Jiranusornkul S

Toxicol. Mech. Methods, 2017 May;27(4):253-271.
PMID: 27996361 DOI: 10.1080/15376516.2016.1273428

In this work, molecular docking, pharmacophore modeling and molecular dynamics (MD) simulation were rendered for the mouse P-glycoprotein (P-gp) (code: 4Q9H) and bioflavonoids; amorphigenin, chrysin, epigallocatechin, formononetin and rotenone including a positive control; verapamil to identify protein-ligand interaction features including binding affinities, interaction characteristics, hot-spot amino acid residues and complex stabilities. These flavonoids occupied the same binding site with high binding affinities and shared the same key residues for their binding interactions and the binding region of the flavonoids was revealed that overlapped the ATP binding region with hydrophobic and hydrophilic interactions suggesting a competitive inhibition mechanism of the compounds. Root mean square deviations (RMSDs) analysis of MD trajectories of the protein-ligand complexes and NBD2 residues, and ligands pointed out these residues were stable throughout the duration of MD simulations. Thus, the applied preliminary structure-based molecular modeling approach of interactions between NBD2 and flavonoids may be gainful to realize the intimate inhibition mechanism of P-gp at NBD2 level and on the basis of the obtained data, it can be concluded that these bioflavonoids have the potential to cause herb-drug interactions or be used as lead molecules for the inhibition of P-gp (as anti-multidrug resistance agents) via the NBD2 blocking mechanism in future.

Matched MeSH terms: Molecular Dynamics Simulation
Molecular Docking and Molecular Dynamics Simulation Studies to Predict Flavonoid Binding on the Surface of DENV2 E Protein

Ismail NA, Jusoh SA

Interdiscip Sci, 2017 Dec;9(4):499-511.
PMID: 26969331 DOI: 10.1007/s12539-016-0157-8

Dengue infections are currently estimated to be 390 million cases annually. Yet, there is no vaccine or specific therapy available. Envelope glycoprotein E (E protein) of DENV mediates viral attachment and entry into the host cells. Several flavonoids have been shown to inhibit HIV-1 and hepatitis C virus entry during the virus-host membrane fusion. In this work, molecular docking method was employed to predict the binding of nine flavonoids (baicalin, baicalein, EGCG, fisetin, glabranine, hyperoside, ladanein, quercetin and flavone) to the soluble ectodomain of DENV type 2 (DENV2) E protein. Interestingly, eight flavonoids were found to dock into the same binding pocket located between the domain I and domain II of different subunits of E protein. Consistent docking results were observed not only for the E protein structures of the DENV2-Thai and DENV2-Malaysia (a homology model) but also for the E protein structures of tick-borne encephalitis virus and Japanese encephalitis virus. In addition, molecular dynamics simulations were performed to further evaluate the interaction profile of the docked E protein-flavonoid complexes. Ile4, Gly5, Asp98, Gly100 and Val151 residues of the DENV2-My E protein that aligned to the same residues in the DENV2-Thai E protein form consistent hydrogen bond interactions with baicalein, quercetin and EGCG during the simulations. This study demonstrates flavonoids potentially form interactions with the E protein of DENV2.

Matched MeSH terms: Molecular Dynamics Simulation*

Filters

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

External Links