Non-empirical study of the phosphorylation reaction catalyzed by 4-methyl-5-β-hydroxyethylthiazole kinase: relevance of the theory of intermolecular interactions

The subject of this study was an analysis of the role of active site residues in the phosphoryl transfer reaction catalyzed by 4-methyl-5-β-hydroxyethylthiazole kinase (ThiK). The ThiK-catalyzed reaction is of special interest due to the lack of a highly conserved aspartate residue serving as a catalytic base. ONIOM(B3LYP:PM3) models of stationary points along the reaction pathway consisted of reactants, two magnesium ions and several highly conserved ThiK active site residues. The results indicate that an S_N2-like mechanism of ThiK, with γ-phosphate acting as an alcohol-activating base is reasonable. Geometries of substrates, transition state and products were utilized in the non-empirical analysis of the physical nature of catalytic interactions taking place in the ThiK active site. The role of particular residues was investigated in terms of their ability to preferentially stabilize the transition state relative to substrates (differential transition state stabilization, DTSS) or products (differential product stabilization, DPS). It seems that Mg2, Glu126 and Cys198 play a major catalytic role, whereas Mg1 and the same Cys198 are responsible for product release. It is remarkable that no dominant role of an electrostatic term in the interactions involved in catalytic activity is observed for product release. Determination of catalytic fields expressing differential electrostatic potential of the transition state with respect to substrates revealed the optimal electrostatic features of an ideal catalyst for the studied reaction. The predicted catalytic environment is in agreement with experimental data showing increased catalytic activity of ThiK upon mutation of Cys198 to aspartate. Figure Catalytic fields for ThiK-catalyzed reaction juxtaposed with the positions of active site residues of a model system. Magnesium ions are considered part of the transition state/reactants. The surface of constant electronic density is colored according to differential electrostatic potential of transition state with respect to reactants. The sign of the differential potential reflects the electrostatic properties of a complementary molecular environment. Red (green) color denotes regions where a negative (positive) charge would be optimal for catalytic activity     

Dynamic domains and geometrical properties of HIV-1 gp120 during conformational changes induced by CD4 binding

The HIV-1 gp120 exterior envelope glycoprotein undergoes a series of conformational rearrangements while sequentially interacting with the receptor CD4 and coreceptor CCR5 or CXCR4 on the surface of host cells to initiate virus entry. Both the crystal structures of the HIV-1 gp120 core bound by the CD4 and antigen 17b, and the SIV gp120 core pre-bound by the CD4 are known. We have performed dynamic domain studies on the homology models of the CD4-bound and unliganded HIV-1 gp120 with modeled V3 and V4 loops to explore details of conformational changes, hinge axes, and hinge bending regions in the gp120 structures upon CD4 binding. Four dynamic domains were clustered and intricately motional modes for domain pairs were discovered. Together with the detailed comparative analyses of geometrical properties between the unliganded and liganded gp120 models, an induced fit model was proposed to explain events accompanying the CD4 engagement to the gp120, which provided new insight into the dynamics of the molecular induced binding mechanism that complements the molecular dynamics and crystallographic studies.     

Thrombin inhibitors with novel P1 binding pocket functionality: free energy of binding analysis

The high incidence of thrombembolic diseases justifies the development of new antithrombotics. The search for a direct inhibitor has resulted in the synthesis of a considerable number of low molecular weight molecules that inhibit human α-thrombin potently. However, efforts to develop an orally active drug remain in progress as the most active inhibitors with a highly basic P1 moiety exhibit an unsatisfactory bioavailability profile. In our previous work we solved several X-ray structures of human α-thrombin in complexes with (1) novel bicyclic arginine mimetics attached to the glycylproline amide and pyridinone acetamide scaffold and (2) inhibitors with a novel aza scaffold and with charged or neutral P1 moieties. In the present contribution, we correlate the structures of the complex between these inhibitors and the protein with the calculated free energy of binding. The energy of solvation was calculated using the Poisson–Boltzmann approach. In particular, the requirements for successful recognition of an inhibitor at the protein’s active site pocket S1 are discussed. Figure We report here on free energy of binding analysis of thrombin inhibitors with novel aza scaffold and novel bicyclic arginine mimetics in S1 pocket of thrombin     

GST pull-down筛选冠突曲霉MAT的互作蛋白*

目的: 冠突曲霉(Aspergillus cristatus)是一种同宗结合菌,它的产孢受渗透压调控,与构巢曲霉的光调控产孢机制存在较大差异。冠突曲霉的有性生殖主要受MAT1-1-1和MAT1-2-1调控,但MAT基因对该菌有性生殖的调控机制仍不清楚。期望筛选得到冠突曲霉MAT的互作蛋白,为深入研究冠突曲霉有性产孢机制奠定基础。方法: 利用GST pull-down联合液相色谱-串联质谱(LC-MS/MS)技术筛选可能与冠突曲霉MAT1-1-1和MAT1-2-1互作的蛋白,结合ProteinPilot和冠突曲霉基因组注释结果进行互作蛋白的注释及GO分析,其中互作蛋白SI65_00917和SI65_03348利用RT-qPCR探索它们与有性发育的联系,并利用酵母双杂交技术初步验证它们与MAT蛋白的互作关系。结果: 成功构建了GST-MAT1-1-1、GST-MAT1-2-1表达载体,诱导表达纯化出目的诱饵蛋白,分别利用诱饵蛋白捕获冠突曲霉总蛋白中的互作蛋白,经分析、筛选共鉴定出与MAT1-1-1互作的蛋白56个,与MAT1-2-1互作的蛋白413个。GO分析表明,这些蛋白参与翻译调控、代谢过程、蛋白质转运及蛋白结合等生物学过程,具有核苷酸结合活性、催化活性、蛋白结合活性;RT-qPCR结果表明互作蛋白SI65_00917可能与有性发育相关。酵母双杂交结果表明,SI65_00917蛋白具有自激活作用,可能是转录因子;SI65_03348蛋白与MAT1-1-1、MAT1-2-1在酵母中均有互作。结论: MAT通过与其他蛋白直接或间接的相互作用调控其有性发育过程。     

TrioMDR: Detecting SNP interactions in trio families with model-based multifactor dimensionality reduction

Single nucleotide polymorphism (SNP) interactions can explain the missing heritability of common complex diseases. Many interaction detection methods have been proposed in genome-wide association studies, and they can be divided into two types: population-based and family-based. Compared with population-based methods, family-based methods are robust vs. population stratification. Several family-based methods have been proposed, among which Multifactor Dimensionality Reduction (MDR)-based methods are popular and powerful. However, current MDR-based methods suffer from heavy computational burden. Furthermore, they do not allow for main effect adjustment. In this work we develop a two-stage model-based MDR approach (TrioMDR) to detect multi-locus interaction in trio families (i.e., two parents and one affected child). TrioMDR combines the MDR framework with logistic regression models to check interactions, so TrioMDR can adjust main effects. In addition, unlike consuming permutation procedures used in traditional MDR-based methods, TrioMDR utilizes a simple semi-parameter P-values correction procedure to control type I error rate, this procedure only uses a few permutations to achieve the significance of a multi-locus model and significantly speeds up TrioMDR. We performed extensive experiments on simulated data to compare the type I error and power of TrioMDR under different scenarios. The results demonstrate that TrioMDR is fast and more powerful in general than some recently proposed methods for interaction detection in trios. The R codes of TrioMDR are available at: https://github.com/TrioMDR/TrioMDR.     

蛋白质功能预测方法概述

蛋白质是生物体内最必需也是最通用的大分子,对它们功能的认识对于科学领域和农业领域的发展有着至关重要的作用。随着后基因组时代的发展,NCBI数据库中迅速涌现出大量不明结构与功能的蛋白质序列,这些蛋白质序列甚至一跃成了研究的热点。近几十年来蛋白质功能预测的方法不断被完善。由最初的仅基于蛋白质序列或3D结构信息的方法衍生出更多的基于序列相似性、基于结构基序、基于相互作用网络等新方法,这些新型方法采用新的算法、新的研究思路和技术手段,力求得到准确性与普遍性并存,能够被广泛应用的蛋白质功能预测方法。本文综述了近年来蛋白质功能预测的方法,并将这些研究方法分类归纳,各自阐明了每类方法的优缺点。     

Ruminal transport and metabolism of short-chain fatty acids (SCFA) in vitro: effect of SCFA chain length and pH

The unidirectional transport and metabolism of 14C-labeled acetate, propionate and butyrate across the isolated bovine rumen epithelium was measured in vitro by the Ussing chamber technique. There was a significant, but relatively small, net secretion of acetate and propionate, and a large and significant net absorption of butyrate. The results demonstrate that the mucosal-serosal (MS) pathway for short-chain fatty acids (SCFA) is different from the serosal-mucosal (SM) pathway, and that butyrate is treated differently from acetate and propionate by the epithelium. The results support that the main route for epithelial SCFA transport is transcellular. The correlation between SCFA lipophility and the flux rate was positive but weak at both pH 7.3 and 6.0. Decreasing pH increased all SCFA fluxes significantly, but not proportionally to the increase of protonized SCFA in the bathing solution. There was a significant and apparently non-competitive interaction between the transport of acetate, propionate and butyrate. It seems that mediated transport mechanisms must be involved in epithelial SCFA transport in the bovine rumen, but the data do not exclude that passive diffusion could account for a significant part of the flux. The metabolism of SCFA in the Ussing chamber system was considerable, and there was a clear preference for excretion of CO2 from this metabolism to the mucosal side, while side preference for non-CO2 metabolite excretion was not studied. Of the propionate and butyrate transported in the MS direction, 78 and 95% was metabolised, while only 37 and 38% was metabolised in the SM direction (acetate metabolism could not be measured). There was, however, no simple relation between the degree of metabolism and the transport rate or the transport asymmetry of the SCFA.     

Molecular dynamics simulations of 14 HIV protease mutants in complexes with indinavir

The molecular mechanisms of HIV drug resistance were studied using molecular dynamics simulations of HIV-1 protease complexes with the clinical inhibitor indinavir. One nanosecond molecular dynamics simulations were run for solvated complexes of indinavir with wild type protease, a control variant and 12 drug resistant mutants. The quality of the simulations was assessed by comparison with crystallographic and inhibition data. Molecular mechanisms that contribute to drug resistance include structural stability and affinity for inhibitor. The mutants showed a range of structural variation from 70 to 140% of the wild type protease. The protease affinity for indinavir was estimated by calculating the averaged molecular mechanics interaction energy. A correlation coefficient of 0.96 was obtained with observed inhibition constants for wild type and four mutants. Based on this good agreement, the trends in binding were predicted for the other mutants and discussed in relation to the clinical data for indinavir resistance. Figure Poincare map representation for WT protease-indinavir complex. The side chain of Tyr 59 showing the positions of hydrogen atoms.This revised version was published online in October 2004 with corrections to the Graphical Abstract.     

Tumour necrosis factor-alpha interacts with biglycan and decorin

The interaction of beta-arrestin-1 with the somatostatin receptor type 2A (sst2A) was monitored using both biochemical and confocal imaging approaches. We show that, using transient transfection of either beta-arrestin-1 or its dominant negative Delta-arrestin-1 in CHO cells stably transfected with the sst2A, beta-arrestin-1 is colocalized with the receptor in endosomal vesicles after somatostatin-induced sequestration. However, this interaction leads to a role of beta-arrestin-1 in the desensitization of the sst2A rather than in the internalization process of the receptor-ligand complex.