Homology modeling,molecular dynamic simulation and docking studies of cyclin dependent kinase 1

In order to develop promising cyclin dependent kinase 1 inhibitors, homology modeling, docking and molecular dynamic simulation techniques were applied to get insight into the functional and structural properties of cyclin dependent kinase 1 (CDK1). Since there is no reported CDK1 crystal structural data, the three dimensional structure of CDK1 was constructed based on homology modeling. An extensive dynamic simulation was also performed on a Flavopiridol-CDK1 complex for probing the binding pattern of Flavopiridol in the active site of CDK1. The binding modes of other inhibitors to CDK1 were also proposed by molecular docking. The structural requirement for developing more potent CDK1 inhibitors was obtained by the above-mentioned molecular simulations and pharmacophore modeling.     

Theoretical model of the three-dimensional structure of a disease resistance gene homolog encoding resistance protein in Vigna mungo

Plant disease resistance (R) genes, the key players of innate immunity system in plants encode 'R' proteins. 'R' protein recognizes product of avirulance gene from the pathogen and activate downstream signaling responses leading to disease resistance. No three dimensional (3D) structural information of any 'R' proteins is available as yet. We have reported a 'R' gene homolog, the 'VMYR1', encoding 'R' protein in Vigna mungo. Here, we describe the homology modeling of the 'VMYR1' protein. The model was created by using the 3D structure of an ATP-binding cassette transporter protein from Vibrio cholerae as a template. The strategy for homology modeling was based on the high structural conservation in the superfamily of P-loop containing nucleoside triphosphate hydrolase in which target and template proteins belong. This is the first report of theoretical model structure of any 'R' proteins.     

Substrate recognition by ribosome-inactivating protein studied by molecular modeling and molecular electrostatic potentials

A computer model of dianthin 30, a type 1 ribosome-inactivating protein (RIP), is constructed by homology modeling using two known X-ray structures; a type 1 RIP, pokeweed antiviral protein (PAP), and chain A of a type 2 RIP, ricin. The 3D structure is refined by molecular dynamics and its binding site compared with those of PAP and ricin using molecular electrostatic potential mapping. The differences in the maps obtained clearly show how, despite the similarity of the topology of the binding site, differences in electrostatic potential can account for the experimentally observed differences in substrate recognition and binding. This demonstrates the potential of these techniques for guiding further experimental analyses.     

人白介素15及突变体的分子设计

以人白介素－２晶体结构为模板同源模型建人白介素１５及其两株突变体（Ｎ端缺失４个氨基酸、Ｃ端缺失３个氨基酸）的空间构象。在ＣＶＦＦ力场下,经过分子力学优化、常温分子动力学模拟获得稳定立体结构模型。借助空间构象残基的亲疏水性分析,利用Ｄｅｌｐｈｉ一性分析蛋白表达静电分布,进行从理论上预测人白介素１５及两株突变体生物学功能的性、差异性。结合分子生物实验,在ＰＢＶ２２０载体中克隆表达获得人白介素１５及两载     

Theoretical Model of the Three-dimensional Structure of a Disease Resistance Gene Homolog Encoding Resistance Protein in Vigna mungo

Abstract

Plant disease resistance (R) genes, the key players of innate immunity system in plants encode ‘R’ proteins. ‘R’ protein recognizes product of avirulance gene from the pathogen and activate downstream signaling responses leading to disease resistance. No three dimensional (3D) structural information of any ‘R’ proteins is available as yet. We have reported a ‘R’ gene homolog, the 'VMYR1′, encoding ‘R’ protein in Vigna mungo. Here, we describe the homology modeling of the 'VMYR1′ protein. The model was created by using the 3D structure of an ATP-binding cassette transporter protein from Vibrio cholerae as a template. The strategy for homology modeling was based on the high structural conservation in the superfamily of P-loop containing nucleoside triphosphate hydrolase in which target and template proteins belong. This is the first report of theoretical model structure of any ‘R’ proteins.     

小鼠TAp63γ及两种突变体的原核表达、纯化和DNA结合活性

采用PCR扩增法得到小鼠TAp63γ野生型及两种缺失突变体的cDNA,3种cDNA与表达载体pGEX-2TK重组构建成GST融合表达质粒并转化感受态E.coli BL21 (DE3),经IPTG诱导了小鼠TAp63γ野生型及两种缺失突变体的可溶性表达. 诱导表达的菌液经离心收集菌体、超声破碎及Triton X-100增溶后获得可溶性表达蛋白粗提液. 利用Glutathione Sepharose 4 Fast Flow亲合层析纯化出电泳均一的3种GST融合蛋白. 凝胶滞留分析证实仅野生型小鼠TAp63γ蛋白能特异结合p53靶序列,经序列比对及同源建模分析,表明小鼠TAp63γ DBD结合区的完整性、关键氨基酸的保守性及三维结构的相似性可能是其DNA结合活性所必需的.     

Multibody coarse-grained potentials for native structure recognition and quality assessment of protein models

Multibody potentials have been of much interest recently because they take into account three dimensional interactions related to residue packing and capture the cooperativity of these interactions in protein structures. Our goal was to combine long range multibody potentials and short range potentials to improve recognition of native structure among misfolded decoys. We optimized the weights for four-body nonsequential, four-body sequential, and short range potentials to obtain optimal model ranking results for threading and have compared these data against results obtained with other potentials (26 different coarse-grained potentials from the Potentials 'R'Us web server have been used). Our optimized multibody potentials outperform all other contact potentials in the recognition of the native structure among decoys, both for models from homology template-based modeling and from template-free modeling in CASP8 decoy sets. We have compared the results obtained for this optimized coarse-grained potentials, where each residue is represented by a single point, with results obtained by using the DFIRE potential, which takes into account atomic level information of proteins. We found that for all proteins larger than 80 amino acids our optimized coarse-grained potentials yield results comparable to those obtained with the atomic DFIRE potential.     

Computational protein structure modeling and analysis of UV-B stress protein in Synechocystis PCC 6803

This study focuses on Ultra Violet stress (UVS) gene product which is a UV stress induced protein from cyanobacteria, Synechocystis PCC 6803. Three dimensional structural modeling of target UVS protein was carried out by homology modeling method. 3F2I pdb from Nostoc sp. PCC 7120 was selected as a suitable template protein structure. Ultimately, the detection of active binding regions was carried out for characterization of functional sites in modeled UV-B stress protein. The top five probable ligand binding sites were predicted and the common binding residues between target and template protein was analyzed. It has been validated for the first time that modeled UVS protein structure from Synechocystis PCC 6803 was structurally and functionally similar to well characterized UVS protein of another cyanobacterial species, Nostoc sp PCC 7120 because of having same structural motif and fold with similar protein topology and function. Investigations revealed that UVS protein from Synechocystis sp. might play significant role during ultraviolet resistance. Thus, it could be a potential biological source for remediation for UV induced stress.     

人EEN蛋白的结构预测和功能分析

采用基于神经网络的算法预测了我们自行克隆的新的白血病相关蛋白ＥＥＮ（ｅｘｔｒａ ｅｌｅｖｅｎｎｉｎｅｔｅｅｎ, ＥＥＮ）全长分子的二级结构,结果表明：ＥＥＮ 蛋白可能有三个结构域,Ｎ 端由三段α螺旋和短β折叠组成,中间为四段α螺旋组成的四螺旋结构,Ｃ端为ＳＨ３结构域,类似于在受体酪氨酸激酶信号传导途径中起重要作用的ＳＥＭ－５／ＧＲＢ２ Ｃ端ＳＨ３结构域;利用同源蛋白结构模拟的方法,模拟了ＥＥＮ ＳＨ３结构域的三维结构,结果表明：ＥＥＮ ＳＨ３结构域与ＳＥＭ－５／ＧＲＢ２ ＳＨ３结构域具有相近的结构,构成脯氨酸结合区的氨基酸非常保守．上述结果提示：ＥＥＮ 蛋白可能为新的信号蛋白,可能涉及新的信号传导途径或新的信号传导旁路,ＳＨ３结构域是其功能区域．     

Main structural and functional features of the basic cytosolic bovine 21 kDa protein delineated through hydrophobic cluster analysis and molecular modelling.

A 21 kDa protein purified from bovine brain cytosol was previously described as a hydrophobic ligand binding protein; however, its accurate biological function remained still uncertain. In order to get further information about its potential biological role, an extended prediction of its secondary and three dimensional structures was undertaken. We describe here a process which permitted us to discover a structural homology between the 21 kDa protein and the N-domain of yeast phosphoglycerate kinase (PGK). This process is based on comparing the 21 kDa protein with all the proteins presenting a slight homology, by using the Hydrophobic Cluster Analysis (HCA) method. According to the observed similarity between the N-domain of yeast PGK and the 21 kDa protein, we built a model which was shown to possess a potential binding site for nucleotides. Moreover, the model obtained presents three-dimensional (3D) structure similarity with adenylate kinase. These results suggest two main hypotheses: (i) the 21 kDa protein may belong to the kinase family; (ii) the binding of a nucleotide could imply a modification of the 3D structure of the 21 kDa protein that can promote the transfer of hydrophobic ligands to the plasma membrane. Meanwhile, verification of these hypotheses has been in part performed experimentally: the 21 kDa protein binds MgATP as well as, to a lesser extent, phosphoglycerate.     

Prediction of 3-D structures of fucose-binding proteins and structural analysis of their interaction with ligands

The importance of fucose-binding proteins stems from the presence of fucose as terminal sugars in H and Lewis (a) blood groups. Recently, the structure of a complex between Anguilla anguilla agglutinin (AAA) and alpha-L-fucose has been worked out at 1.9 A resolution. The structure of AAA characterizes the novel fold of an entire lectin family. In the present study, molecular modeling techniques have been used to identify new proteins that can provide a similar fucose binding module in the newly discovered genomic sequences using the above mentioned structural information. We modeled 3-D structures of three such proteins, namely, ebiP5322 protein of Anopheles gambiae, a pentraxin of Xenopus laevis, and the fw gene product of Drosophila melanogaster. alpha-L-fucose was docked in the binding pockets of the modeled structures followed by energy minimization and molecular dynamic runs to obtain the most probable structures of the complexes. Properties of these modeled complexes were studied to examine the nature of physicochemical forces involved in the complex formation and compared with AAA-alpha-L-fucose complex. It was found that ebiP5322 protein of A. gambiae and the pentraxin of X. laevis can provide a fucose-binding fold similar to AAA. We studied structures of four protein-fucose complexes to examine the electrostatic potential surfaces around the binding site and concluded that a highly positive-charged surface was not a necessary condition of fucose-binding.     

Cooperative DNA looping by PRC2 complexes

Polycomb repressive complex 2 (PRC2) is an essential protein complex that silences gene expression via post-translational modifications of chromatin. This paper combined homology modeling, atomistic and coarse-grained molecular dynamics simulations, and single-molecule force spectroscopy experiments to characterize both its full-length structure and PRC2-DNA interactions. Using free energy calculations with a newly parameterized protein-DNA force field, we studied a total of three potential PRC2 conformations and their impact on DNA binding and bending. Consistent with cryo-EM studies, we found that EZH2, a core subunit of PRC2, provides the primary interface for DNA binding, and its curved surface can induce DNA bending. Our simulations also predicted the C2 domain of the SUZ12 subunit to contact DNA. Multiple PRC2 complexes bind with DNA cooperatively via allosteric communication through the DNA, leading to a hairpin-like looped configuration. Single-molecule experiments support PRC2-mediated DNA looping and the role of AEBP2 in regulating such loop formation. The impact of AEBP2 can be partly understood from its association with the C2 domain, blocking C2 from DNA binding. Our study suggests that accessory proteins may regulate the genomic location of PRC2 by interfering with its DNA interactions.     

cDNA Cloning and Molecular Modeling of Procerain B,a Novel Cysteine Endopeptidase Isolated from Calotropis procera

Procerain B, a novel cysteine protease (endopeptidase) isolated from Calotropis procera belongs to Asclepiadaceae family. Purification of the enzyme, biochemical characterization and potential applications are already published by our group. Here, we report cDNA cloning, complete amino acid sequencing and molecular modeling of procerain B. The derived amino acid sequence showed high sequence homology with other papain like plant cysteine proteases of peptidase C1A superfamily. The three dimensional structure of active procerain B was modeled by homology modeling using X-ray crystal structure of actinidin (PDB ID: 3P5U), a cysteine protease from the fruits of Actinidia arguta. The structural aspect of the enzyme is also discussed.     

An actin‐like protein is involved in regulation of mitochondrial and flagellar functions as well as in intramacrophage survival of Leishmania donovani

Actin‐related proteins are ubiquitous actin‐like proteins that show high similarity with actin in terms of their amino acid sequence and three‐dimensional structure. However, in lower eukaryotes, such as trypanosomatids, their functions have not yet been explored. Here, we show that a novel actin‐related protein (ORF LmjF.13.0950) is localized mainly in the Leishmania mitochondrion. We further reveal that depletion of the intracellular levels of this protein leads to an appreciable decrease in the mitochondrial membrane potential as well as in the ATP production, which appears to be accompanied with impairment in the flagellum assembly and motility. Additionally, we report that the mutants so generated fail to survive inside the mouse peritoneal macrophages. These abnormalities are, however, reversed by the episomal gene complementation. Our results, for the first time indicate that apart from their classical roles in the cytoplasm and nucleus, actin‐related proteins may also regulate the mitochondrial function, and in case of Leishmania donovani they may also serve as the essential factor for their survival in the host cells.     

Homology modeling of an immunoglobulin-like domain in the Saccharomyces cerevisiae adhesion protein alpha-agglutinin.

The Saccharomyces cerevisiae adhesion protein alpha-agglutinin is expressed by cells of alpha mating type. On the basis of sequence similarities, alpha-agglutinin has been proposed to contain variable-type immunoglobulin-like (IgV) domains. The low level of sequence similarity to IgV domains of known structure made homology modeling using standard sequence-based alignment algorithms impossible. We have therefore developed a secondary structure-based method that allowed homology modeling of alpha-aggulutinin domain III, the domain most similar to IgV domains. The model was assessed and where necessary refined to accommodate information obtained by biochemical and molecular genetic approaches, including the positions of a disulfide bond, glycosylation sites, and proteolytic sites. The model successfully predicted surface exposure of glycosylation and proteolytic sites, as well as identifying residues essential for binding activity. One side of the domain was predicted to be covered by carbohydrate residues. Surface accessibility and volume packing analyses showed that the regions of the model that have greatest sequence dissimilarity from the IgV consensus sequence are poorly structured in the biophysical sense. Nonetheless, the utility of the model suggests that these alignment and testing techniques should be of general use for building and testing of models of proteins that share limited sequence similarity with known structures.     

Phenotypic effects of Ehlers�CDanlos syndrome-associated mutation on the FnIII domain of tenascin-X

Tenascin‐X (TNX) is an extracellular matrix (ECM) protein and interacts with a wide variety of molecules in the ECM as well as on the membrane. Deficiency of TNX causes a recessive form of Ehlers–Danlos syndrome (EDS) characterized by hyperelastic and fragile skin, easy bruising, and hypermobile joints. Three point mutations in TNX gene were found to be associated with hypermobility type EDS and one of such mutations is the V1195M mutation at the 7th fibronectin Type III domain (TNXfn7). To help elucidate the underlying molecular mechanism connecting this mutation to EDS, here we combined homology modeling, chemical denaturation, single molecule atomic force microscopy, and molecular dynamics (MD) simulation techniques to investigate the phenotypic effects of V1195M on TNXfn7. We found that the V1195M mutation does not alter the three‐dimensional structure of TNXfn7 and had only mild destabilization effects on the thermodynamic and mechanical stability of TNXfn7. However, MD simulations revealed that the mutation V1195M significantly alters the flexibility of the C′E loop of TNXfn7. As loops play important roles in protein–protein and protein–ligand interactions, we hypothesize that the decreased loop flexibility by V1195M mutation may affect the binding of TNX to ECM molecules and thus adversely affect collagen deposition and fibrillogenesis. Our results may provide new insights in understanding the molecular basis for the pathogenesis of V1195M‐resulted EDS.     

Application of molecular modeling to analysis of inhibition of kinesin motor proteins of the BimC subfamily by monastrol and related compounds

Application of molecular modeling approaches has potential to contribute to rational drug design. These approaches may be especially useful when attempting to elucidate the structural features associated with novel drug targets. In this study, molecular docking and molecular dynamics were applied to studies of inhibition of the human motor protein denoted HsEg5 and other homologues in the BimC subfamily. These proteins are essential for mitosis, so compounds that inhibit their activity may have potential as anticancer therapeutics. The discovery of a small-molecule cell-permeable inhibitor, monastrol, has stimulated research in this area. Interestingly, monastrol is reported to inhibit the human and Xenopus forms of Eg5, but not those from Drosophila and Aspergillus. In this study, homology modeling was used to generate models of the Xenopus, Drosophila, and Aspergillus homologues, using the crystal structure of the human protein in complex with monastrol as a template. A series of known inhibitors was docked into each of the homologues, and the differences in binding energies were consistent with reported experimental data. Molecular dynamics revealed significant changes in the structure of the Aspergillus homologue that may contribute to its relative insensitivity to monastrol and related compounds.