Recent advances in molecular biomarkers for diabetes mellitus: a systematic review

Context: Diabetes is a growing global metabolic epidemic. Current research is focussing on exploring how the biological processes and clinical outcomes of diabetes are related and developing novel biomarkers to measure these relationships, as this can subsequently improve diagnostic, therapeutic and management capacity.

Objective: The objective of this study is to identify the most recent advances in molecular biomarkers of diabetes and directions that warrant further research.

Methods: Using a systematic search strategy, the MEDLINE, CINAHL and OVID MEDLINE databases were canvassed for articles that investigated molecular biomarkers for diabetes. Initial selections were made based on article title, whilst final inclusion was informed by a critical appraisal of the full text of each article.

Results: The systematic search returned 246 records, of which 113 were unique. Following screening, 29 records were included in the final review. Three main research strategies (the development of novel technologies, broad biomarker panels, and targeted approaches) identified a number of potential biomarkers for diabetes including miR-126, C-reactive protein, 2-aminoadipic acid and betatrophin.

Conclusion: The most promising research avenue identified is the detection and quantification of micro RNA. Further, the utilisation of functionalised electrodes as a means to detect biomarker compounds also warrants attention.     

High throughput screening,docking, and molecular dynamics studies to identify potential inhibitors of human calcium/calmodulin-dependent protein kinase IV

Calcium/calmodulin-dependent protein kinase IV (CAMKIV) is associated with many diseases including cancer and neurodegenerative disorders and thus being considered as a potential drug target. Here, we have employed the knowledge of three-dimensional structure of CAMKIV to identify new inhibitors for possible therapeutic intervention. We have employed virtual high throughput screening of 12,500 natural compounds of Zinc database to screen the best possible inhibitors of CAMKIV. Subsequently, 40 compounds which showed significant docking scores (?11.6 to ?10.0?kcal/mol) were selected and further filtered through Lipinski rule and drug likeness parameter to get best inhibitors of CAMKIV. Docking results are indicating that ligands are binding to the hydrophobic cavity of the kinase domain of CAMKIV and forming a significant number of non-covalent interactions. Four compounds, ZINC02098378, ZINC12866674, ZINC04293413, and ZINC13403020, showing excellent binding affinity and drug likeness were subjected to molecular dynamics simulation to evaluate their mechanism of interaction and stability of protein-ligand complex. Our observations clearly suggesting that these selected ligands may be further employed for therapeutic intervention to address CAMKIV associated diseases.

Communicated by Ramaswamy H. Sarma     

3D-QSAR-aided design of potent c-Met inhibitors using molecular dynamics simulation and binding free energy calculation

Abstract

Mesenchymal-epithelial transition factor (c-Met) is a member of receptor tyrosine kinase. It involves in various cellular signaling pathways which includes proliferation, motility, migration, and invasion. Over-expression of c-Met has been reported in various cancers. Hence, it is an ideal therapeutic target for cancer. The main objective of the study is to identify crucial residues involved in the inhibition of c-Met kinase and to design a series of potent imidazo [4,5-b] pyrazine derivatives as c-Met inhibitors. Docking was used to identify important active site residues involved in the inhibition of c-Met kinase which was further validated by 100 ns of molecular dynamics simulation and free energy calculation using molecular mechanics generalized born surface area. Furthermore, binding energy decomposition identified that residues Tyr1230, Met1211, Asp1222, Tyr1159, Met1160, Val1092, Ala1108, and Leu1157 contributed favorably to the binding stability of compound 32. Receptor-guided Comparative Molecular Field Analysis (CoMFA) (q² = 0.751, NOC = 6, r² = 0.933) and Comparative Molecular Similarity Indices Analysis (COMSIA) (q² = 0.744, NOC = 6, r² = 0.950) models were generated based on the docked conformation of the most active compound 32. The robustness of these models was tested using various validation techniques and found to be predictive. The results of CoMFA and CoMSIA contour maps exposed the regions favorable to enhance the activity. Based on this information, 27 novel c-Met inhibitors were designed. These designed compounds exhibited potent activity than the most active compound of the existing dataset.

Communicated by Ramaswamy H. Sarma     

Molecular mechanism study of several inhibitors binding to BRD9 bromodomain based on molecular simulations

Bromodomain-containing protein 9 (BRD9) has been employed as a potential target for anticancer drugs in recent years. In this work, molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and per residue energy decomposition approaches were performed to elucidate the different binding modes between four pyridinone-like scaffold inhibitors and BRD9 bromodomain. Analysis results indicate that non-polar contribution mainly deriving from van der Waals energy is a critical impact on binding affinity of inhibitors against BRD9. Some key residues Phe44, Phe47, Val49, and Ile53 (at ZA loop) enhance the binding energy of inhibitors in BRD9 by means of providing hydrophobic interactions. Moreover, it is observed that BRD9 is anchored by the formation of a stable hydrogen bond between the carbonyl of the inhibitors and the residue Asn100 (at BC loop), and a strong π–π stacking interaction formed between the residue Tyr106 (at BC loop) and the inhibitors. The existence of dimethoxyphenyl structure and the aromatic ring merged to pyridinone scaffold are useful to enhance the BRD9 binding affinity. These findings should guide the rational design of more prospective inhibitors targeting BRD9.

Communicated by Ramaswamy H. Sarma     

Structural insights into the inhibitor binding and new inhibitor design to Polo-like kinase-1 Polo-box domain using computational studies

Polo box domain (PBD) from Polo-Like Kinase-1 (PLK-1) a cell cycle regulator is one of the important non-kinase targets implicated in various cancers. The crystal structure of PLK-1 PBD bound to phosphopeptide inhibitor is available and acylthiourea derivatives have been reported as potent PBD inhibitors. In this work, structure and ligand-based pharmacophore methods have been used to identify new PBD inhibitors. The binding of acylthiourea analogs and new inhibitors to PBD were assessed using molecular docking and molecular dynamics simulations to understand their binding interactions, investigate the complex stability and reveal the molecular basis for inhibition. This study provides the binding free energies and residue-wise contributions to decipher the essential interactions in the protein-inhibitor complementarity for complex formation and the design of new PBD inhibitors with better binding.

Communicated by Ramaswamy H. Sarma     

Extrapolation of phenolic compounds as multi-target agents against cancer and inflammation

Natural products acquire massive structural and chemical diversity, which cannot be coordinated by any synthetic libraries for small molecules and they are continuing to inspire novel discoveries in health sciences. We have performed the computational calculations for geometry optimization and prediction of electronic and structural properties of some plant phenolic compounds through Gaussian 09 program. Energies of molecular orbitals were computed, to mimic out the stabilities arising from charge delocalization and intramolecular interactions. This process indicated the eventual charge transfer within the molecules. The molecular docking and ADMET properties of these compounds with a novel anticancer (HER2) and anti-inflammatory (COX-2) targets revealed that two molecules were capable of inhibiting both the targets, and could be used as multi target inhibitors. Furthermore, molecular dynamics simulation studies were performed to elucidate the binding mechanism and the comparison of inhibitor’s binding mode with diverse biological activities as anticancer and anti-inflammatory agents. A high-quality association was reported among quantum chemical, ADMET, docking, dynamics and MMGBSA results.

Communicated By Ramaswamy H. Sarma     

Synthesis,biological evaluation and molecular modeling studies on novel quinonoid inhibitors of CDC25 phosphatases

The cell division cycle 25 phosphatases (CDC25A, B, and C; E.C. 3.1.3.48) are key regulator of the cell cycle in human cells. Their aberrant expression has been associated with the insurgence and development of various types of cancer, and with a poor clinical prognosis. Therefore, CDC25 phosphatases are a valuable target for the development of small molecule inhibitors of therapeutic relevance. Here, we used an integrated strategy mixing organic chemistry with biological investigation and molecular modeling to study novel quinonoid derivatives as CDC25 inhibitors. The most promising molecules proved to inhibit CDC25 isoforms at single digit micromolar concentration, becoming valuable tools in chemical biology investigations and profitable leads for further optimization.

     

Molecular dynamics and combined docking studies for the identification of Zaire ebola virus inhibitors

Ebola virus (EBOV) is a lethal human pathogen with a risk of global spread of its zoonotic infections, and Ebolavirus Zaire specifically has the highest fatality rate amongst other species. There is a need for continuous effort towards having therapies, as a single licensed treatment to neutralize the EBOV is yet to come into reality. This present study virtually screened the MCULE database containing almost 36 million compounds against the structure of a Zaire Ebola viral protein (VP) 35 and a consensus scoring of both MCULE and CLCDDW docking programs remarked five compounds as potential hits. These compounds, with binding energies ranging from –7.9 to –8.9?kcal/mol, were assessed for predictions of their physicochemical and bioactivity properties, as well as absorption, distribution, metabolism, excretion, and toxicity (ADMET) criteria. The results of the 50?ns molecular dynamics simulations showed the presence of dynamic stability between ligand and protein complexes, and the structures remained significantly unchanged at the ligand-binding site throughout the simulation period. Both docking analysis and molecular dynamics simulation studies suggested strong binding affinity towards the receptor cavity and these selected compounds as potential inhibitors against the Zaire Ebola VP 35. With respect to inhibition constant values, bioavailability radar and other physicochemical properties, compound A (MCULE-1018045960-0-1) appeared to be the most promising hit compound. However, the ligand efficiency and ligand efficiency scale need improvement during optimization, and also validation via in vitro and in vivo studies are necessary to finally make a lead compound in treating Ebola virus diseases.

Communicated by Ramaswamy H. Sarma     

Ser1333 phosphorylation indicates ROCKI activation

Background

Two isoforms of Rho-associated protein kinase (ROCK), ROCKI and ROCKII, play a pivotal role in regulation of cytoskeleton and are involved in multiple cellular processes in mammalian cells. Knockout mice experiments have indicated that the functions of ROCKI and II are probably non-redundant in physiology. However, it is difficult to differentiate the activation status of ROCKI and ROCKII in biological samples. Previously, we have identified phosphorylation site of ROCKII at Ser1366 residue sensitive to ROCK inhibition. We further investigated the activity-dependent phosphorylation site in ROCKI to establish the reagents that can be used to detect their individual activation.

Results

The phosphorylation site of ROCKI sensitive to its inhibition was identified to be the Ser1333 residue. The ROCKI pSer1333-specific antibody does not cross-react with phosphorylated ROCKII. The extent of S1333 phosphorylation of ROCKI correlates with myosin II light chain phosphorylation in cells in response to RhoA stimulation.

Conclusions

Active ROCKI is phosphorylated at Ser1333 site. Antibodies that recognize phospho-Ser1333 of ROCKI and phospho-S1366 residues of ROCKII offer a means to discriminate their individual active status in cells and tissues.     

Exploring binding modes of the selected inhibitors to phosphodiesterase delta by all-atom molecular dynamics simulations and free energy calculations

It’s favorable to alter KRas mutation’s location to endomembrane by interfering the binding of PDEδ (the prenyl-binding protein phosphodiesterase delta) to KRas. In the present work, the binding of four inhibitors (Deltarasin, allyl analogue, pyrazolopyridazinone derivative, and Deltazinone 1) to PDEδ is investigated with all-atom Molecular Dynamic (MD) simulations. The binding free energy calculation results reveal that van der Waals (VDW) energy provides the major force for affinity binding. Moreover, the binding energy decomposition indicates that residues R61 and I129 provide important contributions to binding energies in all systems. The conserved hydrogen bonds play crucial roles in anchoring the inhibitors to the exact site for binding. The results for conformational analysis of PDEδ/free and PDEδ/inhibitors systems show that the structures are more stable after the inhibitors’ binding to the PDEδ. It is also found that the most unstable system among four complexes is PDEδ/pyrazolopyridazinone derivative system whose α3-helix formed by the residues P113-Q116 disappears. This study may provide valuable information for the design of high potency PDEδ inhibitors.

Communicated by Ramaswamy H. Sarma     

Analysis of sulfates on low molecular weight heparin using mass spectrometry: structural characterization of enoxaparin

Introduction: Structural characterization of low molecular weight heparin (LMWH) is critical to meet biosimilarity standards. In this context, the review focuses on structural analysis of labile sulfates attached to the side-groups of LMWH using mass spectrometry. A comprehensive review of this topic will help readers to identify key strategies for tackling the problem related to sulfate loss. At the same time, various mass spectrometry techniques are presented to facilitate compositional analysis of LMWH, mainly enoxaparin.

Areas covered: This review summarizes findings on mass spectrometry application for LMWH, including modulation of sulfates, using enzymology and sample preparation approaches. Furthermore, popular open-source software packages for automated spectral data interpretation are also discussed. Successful use of LC/MS can decipher structural composition for LMWH and help evaluate their sameness or biosimilarity with the innovator molecule. Overall, the literature has been searched using PubMed by typing various search queries such as ‘enoxaparin’, ‘mass spectrometry’, ‘low molecular weight heparin’, ‘structural characterization’, etc.

Expert commentary: This section highlights clinically relevant areas that need improvement to achieve satisfactory commercialization of LMWHs. It also primarily emphasizes the advancements in instrumentation related to mass spectrometry, and discusses building automated software for data interpretation and analysis.     

Understanding the binding of quinoline amines with human serum albumin by spectroscopic and induced fit docking methods

Seven new quinoline-based bioorganic compounds were prepared by solvent-free synthesis and characterized using spectral techniques. The binding of these compounds with human serum albumin (HSA) was investigated by multi-spectroscopic methods. The quenching of Trp fluorescence upon addition of these compounds to HSA confirmed their significant binding. The quenching analysis at three different temperatures revealed that the complex formation is static and the reaction is entropy driven, spontaneous, and exothermic. Hydrogen bonds and van der Waals forces mainly contributed in the interactions as confirmed by the negative ΔH and ΔS values as well as molecular docking. The results from the circular dichroism (CD) spectroscopy indicated the minimal conformational changes of the protein upon binding with these quinoline compounds. The specific binding site and mode of interactions with HSA were also modeled using induced fit molecular docking procedure and their binding site was found to be in the interface of domains II and III, which is similar to the binding of the drug iodipamide with serum albumin.

Communicated by Ramaswamy H. Sarma     

Understanding the molecular consequences of inherited muscular dystrophies: advancements through proteomic experimentation

Introduction: Proteomic techniques offer insights into the molecular perturbations occurring in muscular-dystrophies (MD). Revisiting published datasets can highlight conserved downstream molecular alterations, which may be worth re-assessing to determine whether their experimental manipulation is capable of modulating disease severity.

Areas covered: Here, we review the MD literature, highlighting conserved molecular insights warranting mechanistic investigation for therapeutic potential. We also describe a workflow currently proving effective for efficient identification of biomarkers & therapeutic targets in other neurodegenerative conditions, upon which future MD proteomic investigations could be modelled.

Expert commentary: Studying disease models can be useful for identifying biomarkers and model specific degenerative cascades, but rarely offer translatable mechanistic insights into disease pathology. Conversely, direct analysis of human samples undergoing degeneration presents challenges derived from complex chronic degenerative molecular processes. This requires a carefully planed & reproducible experimental paradigm accounting for patient selection through to grouping by disease severity and ending with proteomic data filtering and processing.     

Screening for the selective inhibitors of MMP-9 from natural products based on pharmacophore modeling and molecular docking in combination with bioassay experiment,hybrid QM/MM calculation,and MD simulation

Matrix metalloproteinase-9 (MMP-9) has been considered as an attractive target involving cancer therapy. In this study, the 3D QSAR pharmacophore model of MMP-9 inhibitors is built, and its reliability is subsequently validated based on different methods. The built pharmacophore model consists of the four chemical features, including two hydrogen bond acceptors (HBA), one hydrophobic (HY), and one ring aromatic (RA). Among them, both HY and RA are found to be especially important features because they involve the interactions of inhibitors with the S1′ pocket of MMP-9, which determines the selectivity of MMP-9 inhibitors. By combining pharmacophore model with molecular docking, the virtual screening is carried out to identify the selective MMP-9 inhibitors from natural products. The four potential selective MMP-9 inhibitors of natural products are found. One of them was used to carry out the bioassay experiment inhibiting MMP-9, and the estimated IC₅₀ value of only 26.94 µM clearly shows its strongly inhibitory activity; besides, both the hybrid quantum mechanics/molecular mechanics (QM/MM) calculation and the molecular dynamics simulation are performed to examine the reliability regarding the binding mode of this inhibitor with MMP-9 active sites predicted by molecular docking. All the screened four natural products are found to well bind with the MMP-9 active sites by different kinds of interactions. Finally, the ADMET properties of screened four natural products are assessed. These screened MMP-9 inhibitors of natural products could be used as the lead compounds to perform structural modifications and optimizations in the future work.

Communicated by Ramaswamy H. Sarma     

Pharmacophore modeling,multiple docking,and molecular dynamics studies on Wee1 kinase inhibitors

Wee1-like protein kinase (Wee1) is a tyrosine kinase that regulates the G₂ checkpoint and prevents entry into mitosis in response to DNA damage. Based on a series of signaling pathways initiated by Wee1, Wee1 has been recognized as a potential target for cancer therapy. To discover potent Wee1 inhibitors with novel scaffolds, ligand-based pharmacophore model has been built based on 101 known Wee1 inhibitors. Then the best pharmacophore model, AADRRR.340, with good partial least square (PLS) statistics (R²?=?0.9212, Q²?=?0.7457), was selected and validated. The validated model was used as a three-dimensional (3D) search query for databases virtual screening. The filtered molecules were further analyzed and refined by Lipinski’s rule of 5, multiple docking procedures (high throughput virtual screening (HTVS), standard precision (SP), genetic optimization for ligand docking (GOLD), extra precision (XP), and unique quantum polarized ligand docking (QPLD)); absorption, distribution, metabolism, excretion, and toxicity (ADMET) screening; and the Prime/molecular mechanics generalized born surface area (MM-GBSA) method binding free energy calculations. Eight leads were identified as potential Wee1 inhibitors, and a 50?ns molecular dynamics (MD) simulation was carried out for top four inhibitors to predict the stability of ligand–protein complex. Molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) based on MD simulation and the energy contribution per residue to the binding energy were calculated. In the end, three hits with good stabilization and affinity to protein were identified.

Communicated by Ramaswamy H. Sarma     

Human p38α mitogen-activated protein kinase in the Asp168-Phe169-Gly170-in (DFG-in) state can bind allosteric inhibitor Doramapimod

Doramapimod (BIRB-796) is widely recognized as one of the most potent and selective type II inhibitors of human p38α mitogen-activated protein kinase (MAPK); however, the understanding of its binding mechanism remains incomplete. Previous studies indicated high affinity of the ligand to a so-called allosteric pocket revealed only in the ‘out’ state of the DFG motif (i.e. Asp168-Phe169-Gly170) when Phe169 becomes fully exposed to the solvent. The possibility of alternative binding in the DFG-in state was hypothesized, but the molecular mechanism was not known. Methods of bioinformatics, docking and long-time scale classical and accelerated molecular dynamics have been applied to study the interaction of Doramapimod with the human p38α MAPK. It was shown that Doramapimod can bind to the protein even when the Phe169 is fully buried inside the allosteric pocket and the kinase activation loop is in the DFG-in state. Orientation of the inhibitor in such a complex is significantly different from that in the known crystallographic complex formed by the kinase in the DFG-out state; however, the Doramapimod’s binding is followed by the ligand-induced conformational changes, which finally improve accommodation of the inhibitor. Molecular modelling has confirmed that Doramapimod combines the features of type I and II inhibitors of p38α MAPK, i.e. can directly and indirectly compete with the ATP binding. It can be concluded that optimization of the initial binding in the DFG-in state and the final accommodation in the DFG-out state should be both considered at designing novel efficient type II inhibitors of MAPK and homologous proteins.

Communicated by Ramaswamy H. Sarma     

Cellular and molecular mechanisms of acute exposure to sulfur mustard: a systematic review

Context: Mustard gas (e.g. sulfur mustard (SM)) has been used as a chemical agent in several battles and is still a potential worldwide menace. Besides local absorption, particularly in the skin, eyes and lungs, systemic spread of the agent also has detrimental effects on gonads, bone marrow and nervous system. Moreover, chronic exposure of SM to respiratory system causes death. Inducing oxidative stress, and disturbing DNA and tissue repair systems, inflammation and cell death signaling pathways have been introduced as molecular mechanisms of the injury.

Methods: In this systematic review, more than 1200 (2000–2014) articles focusing on gross or molecular pathological reports in the acute phase of the respiratory injury after SM exposure were reviewed, followed by two different layers of gross and molecular pathological data (clinic and laboratory) integrated together in a spatio-temporal order. Role of epithelial, neutrophil and macrophage cells and three signaling pathways of inflammation, oxidative stress and cell death are covered in details.

Results and conclusion: Our results propose a critical role of interleukin-17 producing cells in acute and chronic inflammatory responses.