Structure and functions of titin,a giant protein of skeletal and cardiac muscle: Evidence and suppositions

The present-day data on the structural and functional properties of titin in skeletal and cardiac muscle are reviewed.     

New potential inhibitors of mTOR: a computational investigation integrating molecular docking,virtual screening and molecular dynamics simulation

The mTOR (mammalian or mechanistic Target Of Rapamycin), a complex metabolic pathway that involves multiple steps and regulators, is a major human metabolic pathway responsible for cell growth control in response to multiple factors and that is dysregulated in various types of cancer. The classical inhibition of the mTOR pathway is performed by rapamycin and its analogs (rapalogs). Considering that rapamycin binds to an allosteric site and performs a crucial role in the inhibition of the mTOR complex without causing the deleterious side effects common to ATP-competitive inhibitors, we employ ligand-based drug design strategies, such as virtual screening methodology, computational determination of ADME/Tox properties of selected molecules, and molecular dynamics in order to select molecules with the potential to become non-ATP-competitive inhibitors of the mTOR enzymatic complex. Our findings suggest five novel potential mTOR inhibitors, with similar or better properties than the classic inhibitor complex, rapamycin.     

Apta-nanosensor preparation and in vitro assay for rapid Diazinon detection using a computational molecular approach

Aptamers (ss-DNA or ss-RNA), also known as artificial antibodies, have been selected in vitro median to bind target molecules with high affinity and selectivity. Diazinon is one of the most widely used organophosphorus insecticides in developing and underdeveloped countries as insecticide and acaricide. Diazinon is readily absorbed from the gastrointestinal system and rapidly distributed throughout the body. Thus, the design of clinical and laboratory diagnostics using nanobiosensors is necessary. A computational approach allows us to screen or rank receptor structure and predict interaction outcomes with a deeper understanding, and it is much more cost effective than laboratory attempts. In this research, the best sequence (high affinity bind Diazinon-ssDNA) was ranked among 12 aptamers isolated from SELEX experimentation. Docking results, as the first virtual screening stage and static technique, selected frequent conformation of each aptamer. Then, the quantity and quality of aptamer–Diazinon interaction were simulated using molecular dynamics as a mobility technique. RMSD, RMSF, radius of gyration, and the number of hydrogen bonds formed between Diazinon–aptamer were monitored to assess the quantity and quality of interactions. G-quadruplex DNA aptamer (DF20) showed to be a reliable candidate for Diazinon biosensing. The apta-nanosensor designed using simulation results allowed with linearity detection in the range of .141–.65 nM and a LOD of 17.903 nM, and it was validated using a computational molecular approach.     

Developing imidazole analogues as potential inhibitor for Leishmania donovani trypanothione reductase: virtual screening,molecular docking,dynamics and ADMET approach

Visceral leishmaniasis (VL) affects Indian subcontinent, African and South American continent, and it covers 70 countries worldwide. Visceral form of leishmaniasis is caused by Leishmania donovani in Indian subcontinent which is lethal if left untreated. Extensive resistance to antileishmanial drugs such as sodium stibogluconate, pentamidine and miltefosine and their decreased efficacy has been reported in the endemic region. Amphotericin B drug has shown good antileishmanial activity with significant toxicity, but its cost of treatment has limited the outreach of this treatment to affected people living in endemic zone. So, there is an urgent need to identify new antileishmanial drugs with excellent activity and minimal toxicity issues. Trypanothione reductase, a component of antioxidant system, is necessary for parasite growth and survival to raise infection. To develop potential inhibitor, we docked nine hundred and eighty-four 5-nitroimidazole analogues along with clomipramine which is a well-known inhibitor for TR. Total one hundred and forty-seven 5-nitroimidazole analogues with better docking score than clomipramine were chosen for ADMET and QikProp studies. Among these imidazole analogues, total twenty-four imidazole analogues and clomipramine were chosen on the basis of their ADMET, QikProp, and prime MM-GBSA study. Later on, two analogues with best MM-GBSA dG bind were undergone molecular dynamic simulation to ensure protein–ligand interactions. Using above approach, we confirm that ethyl 2-acetyl-5-[4-butyl-2-(3-hydroxypentyl)-5-nitro-1H-imidazol-1-yl]pent-2-enoate can be a drug candidate against L. donovani for the treatment of VL in the Indian subcontinent.     

The effect of fucoidan on tyrosinase: computational molecular dynamics integrating inhibition kinetics

Fucoidan is a complex sulfated polysaccharide extracted from brown seaweed and has a wide variety of biological activities. In this study, we investigated the inhibitory effect of fucoidan on tyrosinase via a combination of inhibition kinetics and computational simulations. Fucoidan reversibly inhibited tyrosinase in a mixed-type manner. Time-interval kinetics showed that the inhibition was processed as first order with biphasic processes. For further insight, we simulated dockings with various sizes of molecular models (monomer to decamer) of fucoidan and showed that the best binding energy change results were obtained from the pentamer (?1.89?kcal/mol) and the hexamer (?1.97?kcal/mol) models of AutoDock Vina. The molecular dynamics simulation confirmed the binding mechanisms between tyrosinase and fucoidan and suggested that fucoidan mostly interacts with several residues including copper ions located in the active site. Our study suggests that fucoidan might be a potential natural antipigment agent.     

Drug discovery for DNA break repair system by screening from TCM database and molecular dynamics approach

In search for novel anti-tumour agents targeting non-homologous end joining pathway, we conducted a virtual screening using traditional Chinese medicine (TCM) on Ku heterodimer which plays an important role on DNA break repair system. Docking results gave glycyrrhizic acid and macedonoside C as potential TCM candidates. Both glycyrrhizic acid and macedonoside C were docked onto the Ku heterodimer with high-predicted binding affinities as evidenced in their Dock Score. These two ligands also interact with key residues which have been previously shown to influence Ku DNA-binding affinity. Both compounds show consistent hydrogen bonding interactions with key residues throughout the 10- ns dynamics simulation run. Furthermore, glycyrrhizic acid and macedonoside C are able to form additional hydrogen bond interactions with positively charged Ku heterodimer surface. Such interactions strengthen the binding interaction of the top two TCM compounds with Ku heterodimer. Glycyrrhizic acid and macedonoside C are products of licorice (Glycyrrhiza glabra), which is a commonly used herb in TCM formulations. Overall, we report glycyrrhizic acid and macedonoside C as potential anti-tumour agents.     

Atomic insight into designed carbamate-based derivatives as acetylcholine esterase (AChE) inhibitors: a computational study by multiple molecular docking and molecular dynamics simulation

Over 100 variants have been designed and studied, using multiple docking methods such as Autodock Vina, ArgusLab, Molegro Virtual Docker, and Hex-Cuda, to study the effect of alteration in the structure of carbamate-based acetylcholyne esterase (AChE) inhibitors. Sixteen selected systems were then subjected to 14 ns molecular dynamics (MD) simulations. Results from all the docking methods are in agreement. Variants that involved biphenyl substituents possess the most negative binding energies in the ?37.64 to ?39.31 kJ mol^?1 range due to their π–π interactions with AChE aromatic residues. The root mean square deviation values showed that all of these components achieved equilibration after 6 ns. Gyration radius (R_g) and solvent accessibility surface area were calculated to further investigate the AChE conformational changes in the presence of these components. MD simulation results suggested that these components might interact with AChE, possibly with no major changes in AChE secondary and tertiary structures.     

Predicting substrate selectivity between UGT1A9 and UGT1A10 using molecular modelling and molecular dynamics approach

Uridine 5′-diphospho-glucuronosyltransferase-1A9 (UGT1A9) expressed in the liver, shows good sequence identity with UGT1A10, expressed in the intestine. Both uridine 5′-diphospho-glucuronosyltransferase (UGT) isoforms show comprehensive overlapping substrate selectivity but there are differences in stereoselectivity, regiospecificity and rate of glucuronidation of the substrates. Multiple sequence alignment analyses of UGT1A9 and UGT1A10 showed that 13% of the residues in N-terminal domain (NTD) are non-identical between them. Herein, authors attempted homology modelling of UGT1A9 and UGT1A10 and validation using software tools and reported mutagenic studies. A molecular docking study of the known substrates is performed on UGT1A9 and UGT1A10 homology models. The non-identical N-terminal residues ranging from 111 to 117 in UGT1A9 and UGT1A10 were identified to play a central role in the substrate selectivity. However, substrate binding is performed by Ser111, Gly115 and Leu117 in UGT1A10 and Gly111, Asp115 and Phe117 in UGT1A9. This study reports new residues in NTD, showing interaction with uridine 5′-diphospho-glucuronic acid which binds with C-terminal domain. Further, molecular dynamics simulations were carried out to study the role of non-identical residues in substrate identification. The study demonstrates the folding of the UGT enzyme, particularly, helix-loop-helix transition and movement of Nα3-2 helix, in response to substrate and co-substrate binding.     

An integrated approach of network-based systems biology,molecular docking,and molecular dynamics approach to unravel the role of existing antiviral molecules against AIDS-associated cancer

A serious challenge in cancer treatment is to reposition the activity of various already known drug candidates against cancer. There is a need to rewrite and systematically analyze the detailed mechanistic aspect of cellular networks to gain insight into the novel role played by various molecules. Most Human Immunodeficiency Virus infection-associated cancers are caused by oncogenic viruses like Human Papilloma Viruses and Epstein–Bar Virus. As the onset of AIDS-associated cancers marks the severity of AIDS, there might be possible interconnections between the targets and mechanism of both the diseases. We have explored the possibility of certain antiviral compounds to act against major AIDS-associated cancers: Kaposi’s Sarcoma, Non-Hodgkin Lymphoma, and Cervical Cancer with the help of systems pharmacology approach that includes screening for targets and molecules through the construction of a series of drug–target and drug–target–diseases network. Two molecules (Calanolide A and Chaetochromin B) and the target “HRAS” were finally screened with the help of molecular docking and molecular dynamics simulation. The results provide novel antiviral molecules against HRAS target to treat AIDS defining cancers and an insight for understanding the pharmacological, therapeutic aspects of similar unexplored molecules against various cancers.     

Identification of InhA inhibitors: A combination of virtual screening,molecular dynamics simulations and quantum chemical studies

In the present work, multiple pharmacophore-based virtual screening of the SPECS natural product database was carried out to identify novel inhibitors of the validated biological target, InhA. The pharmacophore models were built from the five different groups of the co-crystallized ligands present within the active site. The generated models with the same features from each group were pooled and subjected to the test set validation, receiver–operator characteristic analysis and Güner–Henry studies. A set of five hypotheses with sensitivity > 0.5, specificity > 0.5, area under curve (AUC) > 0.7, and goodness of hit score > 0.7 were retrieved and exploited for the virtual screening. The common hits (87 molecules) obtained from these hypotheses were processed via drug-likeness filters. The filtered molecules (27 molecules) were compared for the binding modes and the top scored molecules (12 molecules) along with the reference (triclosan (TCL), docking score = ?11.65 kcal/mol) were rescored and reprioritized via molecular mechanics-generalized Born surface area approach. Eventually, the stability of reprioritized (10 molecules) docked complexes was scrutinized via molecular dynamics simulations. Moreover, the quantum chemical studies of the dynamically stable compounds (9 molecules) were performed to understand structural features essential for the activity. Overall, the protocol resulted in the recognition of nine lead compounds that can be targeted against InhA.     

Capecitabine as a minor groove binder of DNA: molecular docking,molecular dynamics,and multi-spectroscopic studies

The interaction mechanism and binding mode of capecitabine with ctDNA was extensively investigated using docking and molecular dynamics simulations, fluorescence and circular dichroism (CD) spectroscopy, DNA thermal denaturation studies, and viscosity measurements. The possible binding mode and acting forces on the combination between capecitabine and DNA had been predicted through molecular simulation. Results indicated that capecitabine could relatively locate stably in the G-C base-pairs-rich DNA minor groove by hydrogen bond and several weaker nonbonding forces. Fluorescence spectroscopy and fluorescence lifetime measurements confirmed that the quenching was static caused by ground state complex formation. This phenomenon indicated the formation of a complex between capecitabine and ctDNA. Fluorescence data showed that the binding constants of the complex were approximately 2 × 10⁴ M^?1. Calculated thermodynamic parameters suggested that hydrogen bond was the main force during binding, which were consistent with theoretical results. Moreover, CD spectroscopy, DNA melting studies, and viscosity measurements corroborated a groove binding mode of capecitabine with ctDNA. This binding had no effect on B-DNA conformation.     

Mechanical unfolding of a titin Ig domain: structure of transition state revealed by combining atomic force microscopy,protein engineering and molecular dynamics simulations

Titin I27 shows a high resistance to unfolding when subject to external force. To investigate the molecular basis of this mechanical stability, protein engineering Phi-value analysis has been combined with atomic force microscopy to investigate the structure of the barrier to forced unfolding. The results indicate that the transition state for forced unfolding is significantly structured, since highly destabilising mutations in the core do not affect the force required to unfold the protein. As has been shown before, mechanical strength lies in the region of the A' and G-strands but, contrary to previous suggestions, the results indicate clearly that side-chain interactions play a significant role in maintaining mechanical stability. Since Phi-values calculated from molecular dynamics simulations are the same as those determined experimentally, we can, with confidence, use the molecular dynamics simulations to analyse the structure of the transition state in detail, and are able to show loss of interactions between the A' and G-strands with associated A-B and E-F loops in the transition state. The key event is not a simple case of loss of hydrogen bonding interactions between the A' and G-strands alone. Comparison with Phi-values from traditional folding studies shows differences between the force and "no-force" transition states but, nevertheless, the region important for kinetic stability is the same in both cases. This explains the correspondence between hierarchy of kinetic stability (measured in stopped-flow denaturant studies) and mechanical strength in these titin domains.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Introducing a simple model system for binding studies of known and novel inhibitors of AMPK: a therapeutic target for prostate cancer

Prostate cancer (PC) is one of the leading cancers in men, raising a serious health issue worldwide. Due to lack of suitable biomarker, their inhibitors and the platform for testing those inhibitors result in poor prognosis of PC. AMP-activated protein kinase (AMPK) is a highly conserved protein kinase found in eukaryotes that is involved in growth and development, and also acts as a therapeutic target for PC. The aim of the present study is to identify novel potent inhibitors of AMPK and propose a simple cellular model system for understanding its biology. Structural modelling and MD simulations were performed to construct and refine the 3D models of Dictyostelium and human AMPK. Binding mechanisms of different drug compounds were studied by performing molecular docking, molecular dynamics and MM-PBSA methods. Two novel drugs were isolated having higher binding affinity over the known drugs and hydrophobic forces that played a key role during protein–ligand interactions. The study also explored the simple cellular model system for drug screening and understanding the biology of a therapeutic target by performing in vitro experiments.     

Mechanical unfolding of a titin Ig domain: structure of unfolding intermediate revealed by combining AFM,molecular dynamics simulations,NMR and protein engineering

The mechanical unfolding of an immunoglobulin domain from the human muscle protein titin (TI I27) has been shown to proceed via a metastable intermediate in which the A-strand is detached. The structure and properties of this intermediate are characterised in this study. A conservative destabilising mutation in the A-strand has no effect on the unfolding force, nor the dependence of the unfolding force on the pulling speed, indicating that the unfolding forces measured in an AFM experiment are those required for the unfolding of the intermediate and not the native state. A mutant of TI I27 with the A-strand deleted (TI I27-A) is studied by NMR and standard biophysical techniques, combined with protein engineering. Molecular dynamics simulations show TI I27-A to be a good model for the intermediate. It has a structure very similar to the native state, and is surprisingly stable. Comparison with a Phi-value analysis of the unfolding pathway clearly shows that the protein unfolds by a different pathway under an applied force than on addition of denaturant.     

Identification of potent inhibitors against snake venom metalloproteinase (SVMP) using molecular docking and molecular dynamics studies

Snake venom metalloproteinase (SVMP) (Echis coloratus (Carpet viper) is a multifunctional enzyme that is involved in producing several symptoms that follow a snakebite, such as severe local hemorrhage, nervous system effects and tissue necrosis. Because the three-dimensional (3D) structure of SVMP is not known, models were constructed, and the best model was selected based on its stereo-chemical quality. The stability of the modeled protein was analyzed through molecular dynamics (MD) simulation studies. Structure-based virtual screening was performed, and 15 potential molecules with the highest binding energies were selected. Further analysis was carried out with induced fit docking, Prime/MM–GBSA (ΔG_Bind calculations), quantum-polarized ligand docking, and density functional theory calculations. Further, the stability of the lead molecules in the SVMP-active site was examined using MD simulation. The results showed that the selected lead molecules were highly stable in the active site of SVMP. Hence, these molecules could potentially be selective inhibitors of SVMP. These lead molecules can be experimentally validated, and their backbone structural scaffold could serve as building blocks in designing drug-like molecules for snake antivenom.