Sh3glb1/Bif-1 and mitophagy

Evasion of apoptosis, which enables cells to survive and proliferate under metabolic stress, is one of the hallmarks of cancer. We have recently reported that SH3GLB1/Bif-1 functions as a haploinsufficient tumor suppressor to prevent the acquisition of apoptosis resistance and malignant transformation during Myc-driven lymphomagenesis. SH3GLB1 is a membrane curvature-inducing protein that interacts with BECN1 though UVRAG and regulates the post-Golgi trafficking of membrane-integrated ATG9A for autophagy. At the premalignant stage, allelic loss of Sh3glb1 enhances Myc-induced chromosomal instability and results in the upregulation of anti-apoptotic proteins, including MCL1 and BCL2L1. Notably, we found that Sh3glb1 haploinsufficiency increases mitochondrial mass in overproliferated prelymphomatous Eμ-Myc cells. Moreover, loss of Sh3glb1 suppresses autophagy-dependent mitochondrial clearance (mitophagy) in PARK2/Parkin-expressing mouse embryonic fibroblasts (MEFs) treated with the mitochondrial uncoupler CCCP. Interestingly, PARK2-expressing Sh3glb1-deficient cells accumulate ER-associated immature autophagosome-like structures after treatment with CCCP. Taken together, we propose a model of mitophagy in which SH3GLB1 together with the class III phosphatidylinositol 3-kinase complex II (PIK3C3CII) (PIK3R4-PIK3C3-BECN1-UVRAG) regulates the trafficking of ATG9A-containing Golgi-derived membranes (A9⁺GDMs) to damaged mitochondria for autophagosome formation to counteract oncogene-driven tumorigenesis.     

Biosensors for the Detection of Organophosphorous Pesticides

Cytokinins are master regulators of plant growth and development. They are involved in the regulation of many important physiological and metabolic processes. Recent progress in cytokinin research at the molecular level, including identification of related genes and cytokinin receptors, plus elucidation of signal transduction, has greatly increased our understanding of cytokinin actions. Although still in its infant stage, molecular breeding of crops with altered cytokinin metabolism, when combined with the transgenic approach, has shown very promising potential for application to agriculture. In this review we briefly introduce recent progress in cytokinin molecular biology, discuss applications of cytokinin genetic engineering to agriculture, and present implications and future research directions.     

Ultrafast molecular dynamics of biofuel extraction for microalgae and bacteria milking: blocking membrane folding pathways to damaged lipid-bilayer conformations with nanomicelles

Cell milking is a 100% renewable green energy for CO₂ by extraction of biofuels inside the cytosol of photosynthetic micro-organisms as microalgae and bacteria. The cells are exposed to a hydrophobic solvent forming holes and cracks through their membranes from which the biofuels can leak out. In protein folding, the goal would be to find pathways to the unique functional protein conformer. However, in the lipid-bilayer interaction with the solvent for milking, the objective is to block the pathways for damaged membrane conformations of low free energy with undesired nanostructures, using the solvent properties, as shown with an ab initio structural bioinformatic model. Statistical thermodynamics is used to compute the free energy (including entropy) from the molecular dynamics trajectory of the biomolecular system with many conformational changes. This model can be extended to the general problem of biomolecules folding as for proteins and nucleic acids. Using an adaptation of the Einstein diffusion law, the conformational change dynamics of the lipid bilayer depends on the two diffusion coefficients of the solvent: D₁ before the irreversible folding transition time and the much smaller D₂ thereafter. In contrast to the n-hexane and n-heptane hydrocarbons of smaller size, the residual D₂?=?4.7?×?10^?7?cm²/s of the n-decane solvent, with the highest partition coefficient among the three extractors, is the only to present a D₂ value that is significantly below the critical threshold of 10^?6?cm²/s. Therefore, the membrane would resist to long hydrocarbons and the exposed cells would remain viable for milking.     

3D-QSAR,molecular docking and molecular dynamics studies of a series of RORγt inhibitors

The discovery of clinically relevant inhibitors of retinoic acid receptor-related orphan receptor-gamma-t (RORγt) for autoimmune diseases therapy has proven to be a challenging task. In the present work, to find out the structural features required for the inhibitory activity, we show for the first time a three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations for a series of novel thiazole/thiophene ketone amides with inhibitory activity at the RORγt receptor. The optimum CoMFA and CoMSIA models, derived from ligand-based superimposition I, exhibit leave-one-out cross-validated correlation coefficient (R²_cv) of .859 and .805, respectively. Furthermore, the external predictive abilities of the models were evaluated by a test set, producing the predicted correlation coefficient (R²_pred) of .7317 and .7097, respectively. In addition, molecular docking analysis was applied to explore the binding modes between the inhibitors and the receptor. MD simulation and MM/PBSA method were also employed to study the stability and rationality of the derived conformations, and the binding free energies in detail. The QSAR models and the results of molecular docking, MD simulation, binding free energies corroborate well with each other and further provide insights regarding the development of novel RORγt inhibitors with better activity.     

Molecular dynamics simulation studies of GSK-3β ATP competitive inhibitors: understanding the factors contributing to selectivity

Glycogen synthase kinase-3 is a constitutively acting, multifunctional serine threonine kinase, the role of which has been implicated in several physiological pathways and has emerged as a promising target for the treatment of type-II diabetes and Alzheimer’s disease. In order to provide a detailed understanding of the origin of selectivity determinants of ATP competitive inhibitors, molecular dynamics simulations in combination with MM-PBSA binding energy calculations were performed using crystal structures of GSK-3β and CDK-2 in complex with 12 ATP competitive inhibitors. Analysis of energy contributions indicate that electrostatic interaction energy dictates the selectivity of ATP competitive inhibitors against CDK-2. Key interactions as well as residues that potentially make a major contribution to the binding free energy were identified at the ATP binding site. This analysis stresses the need for the inhibitors to interact with Lys85, Thr138, and Arg141 in the binding site of GSK-3β to show selectivity. The residue-wise energy decomposition analysis further suggested the additional role of Gln185 in determining the selectivity of maleimides. The results obtained in this study can be utilized to design new selective GSK-3 ATP competitive inhibitors.     

Computational DNA binding studies of (–)-epigallocatechin-3-gallate

The catechin family of molecules that are present in the leaves of green tea has been under investigation since the antioxidant and anti-inflammatory properties of tea were discovered. Among multiple proposed therapeutic targets of these molecules, the direct interaction with nucleic acids has been proposed and experimentally observed but without clear knowledge about the potential binding modes between these ligands and DNA. One of these catechin structures, (–)-epigallocatechin gallate (EGCG), has three aromatic rings that could interact with double-stranded DNA via terminal base-pair stacking, intercalation, or through groove binding. Using enhanced sampling techniques and molecular dynamics simulations, we have found a stable complex between the EGCG ligand and DNA through intercalation of the trihydroxybenzoate aromatic ring and an ApC step. Moreover, we have calculated the absorption spectra of four possible binding modes and compared these to absorption profiles reported in the literature, and explored the possible DNA sequence preference for the EGCG ligand to bind. Our results suggest that an intercalative mode of interaction through the major groove is possible between the EGCG ligands and DNA with apparently very little DNA sequence selectivity.     

Effect of Solute Adsorption Properties on Its Separation from Supercritical Carbon Dioxide with a Thin Porous Silica Membrane

A thin porous silica membrane (average pore size of 3.3 mm) was prepared by the sol–gel method and used to separate the solute from supercritical carbon dioxide. The characteristics of solute permeation were investigated in respect of the adsorption properties of the solute, the desorption rate of the solute from the membrane being measured and the potential energy of solute near the silica surface being calculated by the molecular modeling technique. It was found that caffeine was strongly adsorbed to the surface and then slowly desorbed to form an adsorption layer, making the pores narrower and causing a molecular-sieving effect. Therefore, the rejection value was positive. On the other hand, the rejection value of n-octanoic acid, which was well adsorbed and rapidly desorbed, was negative. It is presumed that the molecules filled the pores due to their potential energy and were then forced to flow through the pores by the transmembrane pressure.     

Forced Unfolding of Apocytochrome <Emphasis Type="Italic">b</Emphasis><Subscript>5</Subscript> by Steered Molecular Dynamics Simulation

Apocytochrome b5 (apocyt b5), a small b-type cytochrome with heme prosthetic group removal, has been subjected to steered molecular dynamics (SMD) simulations for investigating the consequences of mechanical force-induced unfolding. Both constant velocity (0.5 and 1.0 A/ps) and constant force (500, 750 and 1000 pN) stretching have been employed to model forced unfolding of apocyt b5. The results of SMD simulations elucidate that apocyt b5 is protected against external stress mainly through the interstrand hydrogen bonding between its beta1-beta2 and beta2-beta3 strands, highlighting the importance of hydrophobic core 2 in stabilization of apocyt b5. The existence of intermediate states manifested by current simulations in the forced unfolding pathway of apocyt b5 is different from the observations in pervious thermal or chemical unfolding studies in the absence of force. The present study could thus provide insights into the relationship between the two cooperative functional modules of apocyt b5 and also guide the rational molecular design of heme proteins.     

The mechanism of high affinity pentasaccharide binding to antithrombin,insights from Gaussian accelerated molecular dynamics simulations

Abstract

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from ‘static’ structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a ‘non-activated’ AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. ‘Generalized correlation’ matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.

Communicated by Ramaswamy H. Sarma     

植物病媒昆虫的翅型分化

翅多型现象是昆虫非遗传多型性的一种表现,包括不具飞行能力的短翅型或无翅型,以及可以进行长距离迁飞的长翅型或有翅型。翅多型现象常发生在可以携带病原并将其传播给植物宿主的媒介昆虫中,对植物病害的时空分布与暴发有重要影响。本文从翅型分化的遗传规律、诱导因素、分子机制和伴随翅型分化的其他生理表现4个方面,对植物病原主要传播媒介蚜虫和飞虱的翅型分化研究进行综述和梳理。昆虫翅型分化的诱导因素主要包括温度、湿度和光周期等非生物因素以及虫口密度、宿主营养、病毒等生物因素;而其内在的分子机制大多是通过胰岛素/胰岛素样生长因子信号(IIS)通路、c-Jun氨基末端激酶(c-Jun NH₂-terminal kinase, JNK)信号通路、Wingless和嗅觉受体SaveOrco等调控。翅型分化的同时伴随着生理状态的变化,表现为短翅型具有更强的繁殖能力和长翅型含有更丰富的飞行肌结构成分。目前,昆虫翅型分化的研究尚不够完善,有许多需要解答的问题,如找到胰岛素/胰岛素样生长因子信号通路中真正发挥功能的靶基因,JNK如何调控翅型分化以及虫媒病毒影响媒介昆虫翅型的分子机理。本综述可为控制虫...