Investigation into the site-specific binding interactions between chlorogenic acid and ovalbumin using multi-spectroscopic and in silico simulation studies

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a promising antimicrobial target involved in the folate biosynthesis pathway. Although, the results from crystallographic studies of HPPK have attracted a great interest in the design of novel HPPK inhibitors, the mechanism of action of HPPK due to inhibitor binding remains questionable. Recently, mercaptoguanine derivatives were reported to inhibit the pyrophosphoryl transfer mechanism of Staphylococcus aureus HPPK (SaHPPK). The present study is an attempt to understand the SaHPPK-inhibitors binding mechanism and to highlight the key residues that possibly involve in the complex formation. To decipher these questions, we used the state-of-the-art advanced insilico approach such as molecular docking, molecular dynamics (MD), molecular mechanics-generalized Born surface area approach. Domain cross correlation and principle component analysis were applied to the snapshots obtained from MD revealed that the compounds with high binding affinity stabilize the conformational dynamics of SaHPPK. The binding free energy estimation showed that the van der Waals and electrostatic interactions played a vital role for the binding mechanism. Additionally, the predicted binding free energy was in good agreement with the experimental values (R2 = .78). Moreover, the free energy decomposition on per-residue confirms the key residues that significantly contribute to the complex formation. These results are expected to be useful for rational design of novel SaHPPK inhibitors.

Scaling Behavior and Structure of Denatured Proteins

Major histocompatibility complex class I (MHC-I), a key element of the acquired immune system, plays essential roles in activating CD8+ T cells by recognizing intracellular antigens derived from pa...

Role of decaethoxylated stearylamine in the selective flotation of hornblende and siderite: An experimental and molecular dynamics simulation study

Matrix metalloproteinase-9 (MMP-9) is a zinc and calcium-dependent proteolytic enzyme involved in extracellular matrix degradation. Overexpression of MMP-9 has been confirmed in several disorders, including cancers, Alzheimer's disease, autoimmune diseases, cardiovascular diseases, and dental caries. Therefore, MMP-9 inhibition is recommended as a therapeutic strategy for combating various diseases. Cinnamic acid derivatives have shown therapeutic effects in different cancers, Alzheimer's disease, cardiovascular diseases, and dental caries. A computational drug discovery approach was performed to evaluate the binding affinity of selected cinnamic acid derivatives to the MMP-9 active site. The stability of docked poses for top-ranked compounds was also examined. Twelve herbal cinnamic acid derivatives were tested for possible MMP-9 inhibition using the AutoDock 4.0 tool. The stability of the docked poses for the most potent MMP-9 inhibitors was assessed by molecular dynamics (MD) in 10 nanosecond simulations. Interactions between the best MMP-9 inhibitors in this study and residues incorporated in the MMP-9 active site were studied before and after MD simulations. Cynarin, chlorogenic acid, and rosmarinic acid revealed a considerable binding affinity to the MMP-9 catalytic domain (ΔGbinding < -10 kcal/ mol). The inhibition constant value for cynarin and chlorogenic acid were calculated at the picomolar scale and assigned as the most potent MMP-9 inhibitor from the cinnamic acid derivatives. The root-mean-square deviations for cynarin and chlorogenic acid were below 2 Å in the 10 ns simulation. Cynarin, chlorogenic acid, and rosmarinic acid might be considered drug candidates for MMP-9 inhibition.

The p300/CBP associated factor bromodomain (PCAF Brd) is emerged as one of the promising target proteins for different types of cancers. PCAF is one among the histone acetyltransferase enzymes which involved in the regulation of transcriptase process by modifying the chromatin structure. Anacardic acid, carnosol, garcinol are the experimentally reported inhibitors of PCAF Brd; however, their detailed binding mechanism these inhibitors are not yet known. The intermolecular interaction, binding energy, and the stability of these inhibitors with the active site of PCAF Brd are playing the key role in the binding of these inhibitors with PCAF. The in silico study incorporates the molecular docking and dynamics simulations; these molecular level simulations allow to understand the binding mechanism. In the present study, the induced fit molecular docking and molecular dynamics of anacardic acid, carnosol and garcinol molecules against the PCAF Brd have been performed. The docking score values of these molecules are -5.112 (anacardic acid), -5.141 (carnosol), -5.199 (garcinol) and -3.641 (L45) kcal/mol, respectively. Further, the molecular dynamics simulation was carried out for these docked complexes to understand their conformational their stability and binding energy from the roots means square deviation (RMSD) and root means square of fluctuation (RMSF), and molecular mechanics with the generalized born and surface area solvation (MM/GBSA) binding free energy calculations. The intermolecular interactions and binding free energy values confirm that garcinol forms key interactions and has high binding affinity towards PCAF Brd on compare with the other two inhibitors. Therefore, garcinol may be considered as a potential inhibitor of PCAF Brd.

Effective anti-inflammatory phenolic compounds from dandelion: identification and mechanistic insights using UHPLC-ESI-MS/MS, fluorescence quenching and anisotropy, molecular docking and dynamics simulation

Chlorogenic acid (CGA), an important metabolite in natural plant medicines such as honeysuckle and eucommia, has been shown to have potent antinociceptive effects. Nevertheless, the mechanism by which CGA relieves chronic pain remains unclear. α-amino-3-hydroxy-5-methyl-4-isooxazolpropionic acid receptor (AMPAR) is a major ionotropic glutamate receptor that mediates rapid excitatory synaptic transmission and its glutamate ionotropic receptor AMPA type subunit 1 (GluA1) plays a key role in nociceptive transmission. In this study, we used Western blot, surface plasmon resonance (SPR) assay, and the molecular simulation technologies to investigate the mechanism of interaction between CGA and AMPAR to relieve chronic pain. Our results indicate that the protein expression level of GluA1 showed a dependent decrease as the concentration of CGA increased (0, 50, 100, and 200 μM). The SPR assay demonstrates that CGA can directly bind to GluA1 (KD = 496 μM). Furthermore, CGA forms a stable binding interaction with GluA1, which is validated by molecular dynamics (MD) simulation. The binding free energy between CGA and GluA1 is −39.803 ± 14.772 kJ/mol, where van der Waals interaction and electrostatic interaction are the major contributors to the GluA1–CGA binding, and the key residues are identified (Val-32, Glu-33, Ala-36, Glu-37, Leu-48), which play a crucial role in the binding interaction. This study first reveals the structural basis of the stable interaction between CGA and GluA1 to form a binding complex for the relief of chronic pain. The research provides the structural basis to understand the treatment of chronic pain and is valuable to the design of novel drug molecules in the future.

To obtain a scientific thought and expedition to explore key interactions with Tyr48 in aldose reductase (ALR), combined study of pharmacophore modeling, induced fit docking, and dynamics studies were performed on ALR. A stereo chemically and energetically valid model of ALR-NADP+ complex was developed using homology modeling technique. Statistically a significant five point pharmacophore model was designed on a set of 54 thiazolidinedione derivatives with good external and internal predictive ability. Rigid and induced fit docking protocols were applied on ALR protein for both with and without NADP+ cofactor to identify a suitable binding mode that facilitates the key hydrogen bond interactions with Tyr48. Docking of thiazolidinedione derivatives into ALR-NADP+ complex gave more promising results by reducing false positive binding of inhibitors into the co-factor binding site. Structural changes within Try48 and Asp43 during the binding process in enzyme inhibitor complex were analyzed using molecular dynamics (MD) simulations. The results obtained from dynamic simulations emphasized the role of Tyr48 in maintaining inter or intra molecular hydrogen bond interaction with the protein or inhibitor respectively. New molecules were designed and checked for their binding interactions and showed improved results compared to existing thiazolidinediones derivatives. Hence, these combined protocols will be helpful and cooperative to design and optimize molecules with better inhibitory activity against the biologically active target.

Phenolic environmental endocrine-disrupting chemicals (PEDCs) are persistent EDCs that are widely found in food packaging materials and environmental media and seriously threaten human health and ecological security. Human estrogen-related receptor γ (hERRγ) has been proposed as a mediator for the low-dose effects of many environmental PEDCs; however, the atomic-level descriptions of dynamical structural features and interactions of hERRγ and PEDCs are still unclarified. Herein, how three PEDCs, 4-(1-methylpropyl)phenol (4-sec-butylphenol), 5,6,7,8-tetrahydro-2-naphthol (tetrahydro-2-napthol), and 2,2-bis(4-hydroxy-3,5-dimethoxyphenyl)propane (BP(2,2)(Me)), interact with hERRγ to produce its estrogenic disruption effects was studied. Molecular docking and multiple molecular dynamics (MD) simulations were first conducted to distinguish the detailed interaction pattern of hERRγ with PEDCs. These binding structures revealed that residues around Leu271, Leu309, Leu345, and Phe435 are important when binding with PEDCs. Furthermore, the binding energies of PEDCs with hERRγ were also characterized using the molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) and solvated interaction energy (SIE) methods, and the results showed that the interactions of CH-π, π-π, and hydrogen bonds are the major contributors for hERRγ binding to these three PEDCs. What is striking is that the methoxide groups of BP(2,2)(Me), as hydrophobic groups, can help to reduce the binding energy of PEDCs binding with hERRγ. These results provide important guidance for further understanding the influence of PEDCs on human health problems.

Unraveling the pH-dependent binding mechanism of acrylamide to β-lactoglobulin: Combined insights from multi-spectroscopy, molecular docking, and molecular dynamics simulations.

A comprehensive theoretical study of the hydrogen bonding interactions and microscopic solvation structures of a pyridyl-urea-based hydrogelator in aqueous solution

By performing homology modeling, molecular docking, and molecular dynamics (MD) simulations, we have developed three-dimensional (3D) structural models of the M5 muscarinic acetylcholine receptor (mAChR) and two complexes for M5 mAChR binding with antagonists SVT-40776 and solifenacin in the environment of lipid bilayer and solvent water. According to the simulated results, each of the antagonists is oriented horizontally in the binding pocket formed by transmembrane helices 2, 3, and 5-7. The cationic headgroup of each of the antagonists interacts with a negatively charged residue, Asp110, through electrostatic and hydrogen-bonding interactions. The simulated results also reveal some significant difference between the binding modes of SVT-40776 and solifenacin. In particular, SVT-40776 is persistently hydrogen bonded with the side chain of residue Tyr458, whereas solifenacin cannot form a similar hydrogen bond with residues around its carbonyl group. Such significant difference in the binding structures is consistent with the fact that SVT-40776 has a much higher binding affinity (K(d) = 0.4 nM) to M5 mAChR than that of solifenacin (K(d) = 31 nM) with the same reeptor. The calculated binding free energy change (-2.3 ± 0.3 kcal/mol) from solifenacin to SVT-40776 is in good agreement with the experimentally derived binding free energy change (-2.58 kcal/mol), suggesting that our modeled M5 mAChR structure and its complexes with the antagonists are reliable. The new structural insights obtained from this computational study are expected to stimulate further biochemical and pharmacological studies on the detailed structures of M5 and other subtypes of mAChRs.

This study aimed to identify B-cell epitope candidates using multiple epitope identification software and in silico analysis of the modeled B19 V protein against specific antibodies using molecular docking and dynamics simulation. Materials and Methods : Full-length amino acid sequences of the VP1 protein of B19 V were retrieved from NCBI. A consensus sequence was generated using CLC sequence viewer. Linear B cell epitopes were identified using Bepipred 2.0, ABCpred, and LBTope. The linear epitope was synthesized and validated against B19 V-specific antibodies. A 3D model of the B19 V VP1 consensus protein was generated using the ITASSER server. Discontinuous B cell epitopes were identified using Discotope 2.0 and Ellipro. Molecular docking and molecular dynamics simulation was performed to investigate the interaction between the modeled B19 V protein and specific anti-B19 V antibody. Results : The identified epitope was 100% conserved and similarly identified through ABCpred and LBTope. The HADDOCK score and MDS analysis, such as hydrogen bond interactions and MMPBSA analysis, revealed that the VP1 and mAb H chains formed a significantly stable complex. The MDS demonstrated that the VP1-mAb H chain complexes had lower RMSF values around 130 to 200 residues, a region responsible for the catalytic network for enzyme activity; as a result, the flexibility of the antibody-bound VP1 decreased when compared to Apo-VP1. Conclusion: A viable epitope identified through this process was synthesized and validated using ELISA, which highlighted the role of the epitope identification process in diagnostics. This study also sheds light on the complex interplay between VP1 and the mAb H chain and highlights key binding specificity and stability determinants.

Interaction mechanism between resveratrol and ovalbumin based on fluorescence spectroscopy and molecular dynamic simulation