Homology modeling and consensus protein disorder prediction of human filamin

Filamins are dimeric actin-binding proteins participating in the organization of the actin-based cytoskeleton. Their modular domain organization is made up of an N-terminal actin-binding domain composed of two CH domains followed by flexible rod regions that consist of 24 Ig-like domains. Homology modeling was used to model human filamin using Modeller 9v5. The resulting model assessed by Verify 3D and PROCHECK showed that the final model is reliable. The conformational disorder prediction of human filamin residues were also mapped on the validated structure of human filamin. Prediction of protein disorder in filamin structures will help structural biologists to find suitable targets to be analyzed and for understanding protein function.     

米根霉富马酸酶三维结构的同源建模与分析

利用Modeller7v7软件对米根霉（Rhizopus oryzoe）富马酸酶（fumarase）进行了三级结构的同源建模并对结果的空间和能量上的合理性进行了验证,进一步对酶的结构域和催化活性位点进行了研究。结果表明：富马酸酶由三个结构域组成,中心区域为一个由五个几乎平行的α螺旋组成的独特的束型结构,其催化活性位点是由三个亚基上的氨基酸相互靠近共同组成的。为以后有针对性的进行富马酸酶的定点突变提高富马酸产量提供分子水平上的理论指导。     

Homology modeling threedimensional structure of AnxB1 and reducing its immunogenicity by sequence-deleted mutagenesis

AnxB1,a novel annexin previously isolated from Cysticercus cellulose,shows high thrombi affinity and anticoagulant activity in vivo.In order to investigate the relationship between structure and biological function,a predicted three-dimensional(3D)model of AnxB1 was generated by homology modeling.This model contains four homologous internal-domains and the Cα trace of domain Ⅰ,Ⅱ and IV shows high similarity.Based on the structure characterization,four sequence-deleted mutants were constructed and expressed as GST fusion proteins in E.coli.Two of the mutants,GST-M3 and GST-M4 reserved high anticoagulant activity(p＜0.01 vs.GST).Furthermore,compared with the wild type GST-AnxB1,the immunogenicity of GST-M3 and GST-M4 was reduced significantly(p＜0.01)and the molecular weight was lowered to 27 kD and 34 kD,respectively.These observations laid a solid foundation for further study on developing new thrombolytic agents with higher efficiency and lower side effect.     

Homology modeling three-dimensional structure of AnxB1 and reducing its immunogenicity by sequence-deleted mutagenesis

We previously isolated and functionally identified a new annexin from Cysticercus cellulosae, which belongs to a novel annexin subfamily and is desig-nated as AnxB1[1]. In the presence of Ca2+, AnxB1 can bind to negatively charged phospholipids like phos-phatidylserine (PS) with high affinity, whereas hardly associates with phosphatidylcholine (PC) even at high Ca2+ concentration. PS is the major component of the platelet membrane and plays a key role in the process of coagulation cascade.…     

Homology modeling provides insights into the binding mode of the PAAD/DAPIN/pyrin domain,a fourth member of the CARD/DD/DED domain family

The PAAD/DAPIN/pyrin domain is the fourth member of the death domain superfamily, but unlike other members of this family, it is involved not only in apoptosis but also in innate immunity and several other processes. We have identified 40 PAAD domain-containing proteins by extensively searching the genomes of higher eukaryotes and viruses. Phylogenetic analyses suggest that there are five categories of PAAD domains that correlate with the domain architecture of the entire proteins. Homology models built on CARD and DD structures identified functionally important residues by studying conservation patterns on the surface of the models. Surface maps of each subfamily show different distributions of these residues, suggesting that domains from different subfamilies do not interact with each other, forming independent regulatory networks. Helix3 of PAAD is predicted to be critical for dimerization. Multiple alignment analysis and modeling suggest that it may be partly disordered, following a new paradigm for interaction proteins that are stabilized by protein-protein interactions.     

Influence of human tryptophan hydroxylase 2 N- and C-terminus on enzymatic activity and oligomerization

Tryptophan hydroxylase (TPH) catalyses the first and rate limiting step in the biosynthesis of the neurotransmitter serotonin. There are two TPH isoenzymes in humans, encoded by two different genes: TPH1 and the recently described TPH2. We have expressed both human enzymes and various deletion mutants of TPH2 (DeltaN44, DeltaC17, DeltaC19, DeltaC51) in COS7 cells. TPH1 and 2 displayed different kinetic properties with a lower K(m) value of TPH1. Removal of 44 amino acids from the N-terminus of TPH2 resulted in a 3-4-fold increased V(max), which indicates a strong inhibitory function of this part on the enzymes activity. TPH1 and 2 were able to form homooligomers and also heterooligomers with each other. The different deletion mutants (DeltaC17, DeltaC19 and DeltaC51), which lack the putative C-terminal leucine zipper tetramerization domain, existed as monomeric enzymes. While short deletions (DeltaC17 and DeltaC19) hardly changed V(max) values, the DeltaC51 mutant lost 99% of TPH activity. These data identify a region between the C-terminal oligomerization domain and the catalytic domain, which is indispensable for TPH2 activity.     

Homology modeling of human histone deacetylase 10 and design of potential selective inhibitors

Abstract

Histone deacetylases (HDACs) are implicated in the pathology of various cancers, and their pharmacological blockade has proven to be promising in reversing the malignant phenotypes. However, lack of crystal structures of some of the human HDAC isoforms (e.g., HDAC10) hinders the design of the isoform-selective inhibitor. Here, the recently solved X-ray crystal structure of Danio rerio (zebrafish) HDAC10 (Protein Data Bank (PDB) ID; 5TD7, released on 24 May 2017) was retrieved from the PDB and used as a template structure to model the three-dimensional structure of human HDAC10. The overall quality of the best model (M0017) was assessed by computing its z-score—a measure of the deviation of the total energy of the structure with respect to an energy distribution derived from random conformations and by docking of known HDAC10 inhibitors to its catalytic cavity. Furthermore, to identify potential HDAC10-selective inhibitor ligand-based virtual screening was carried out against the ZINC database. The free modeled structure of HDAC10 and its complexes with quisinostat and the highest-ranked compound ZINC19749069 were submitted to molecular dynamics simulation. The comparative analysis of root-mean-squared deviation, root-mean-squared fluctuation, radius of gyration (Rg), and potential energy of these systems showed that HDAC10-ZINC19749069 complex remained the most stable over time. Thus, M0017 could be potentially used for structure-based inhibitor against HDAC10, and ZINC19749069 may provide a scaffold for further optimization.

Communicated by Ramaswamy H. Sarma     

Homology modeling and molecular dynamics simulations of the alpha1 glycine receptor reveals different states of the channel

Homology modeling is used to build initial models of the transmembrane domain of the human alpha1 glycine receptor (GlyR) based on the most recently published refined structure of nAChR (PDB ID: 2BG9). Six preliminary GlyR models are constructed using two different approaches. In one approach, five different homopentamers are built by symmetric assembly of alpha1 GlyR subunits using only one of the five unique chains of nAChR as a template. In a second approach, each nAChR subunit serves as a template for an alpha1 GlyR subunit. All six initial GlyR constructs are then embedded into a hydrated POPC lipid bilayer and subjected to molecular dynamics simulation for at least six nanoseconds. Each model is stable throughout the simulation, and the final models fall into three distinct categories. Homopentameric GlyR bundles using a single alpha nAChR subunit as a template appear to be in an open conformation. Under an applied external potential, permeation of Cl(-) ions is observed within several ns in a channel built on an alpha chain. Model channels built on non-alpha chains have a constriction either near the intracellular mouth or more centrally located in the pore domain, both of which may be narrow enough to close the channel and whose locations correspond to putative gates observed in nicotinicoid receptors. The differences between these three general models suggest that channel closure may be effected by either rotation or tangential tilting of TM2.     

Walking through the protein sequence space: towards new generation of the homology modeling

A new method is proposed to reveal apparent evolutionary relationships between protein fragments with similar 3D structures by finding "intermediate" sequences in the proteomic database. Instead of looking for homologies and intermediates for a whole protein domain, we build a chain of intermediate short sequences, which allows one to link similar structural modules of proteins belonging to the same or different families. Several such chains of intermediates can be combined into an evolutionary tree of structural protein modules. All calculations were made for protein fragments of 20 aa residues. Three evolutionary trees for different module structures are described. The aim of the paper is to introduce the new method and to demonstrate its potential for protein structural predictions. The approach also opens new perspectives for protein evolution studies.     

Homology modeling of human serum paraoxonase1 and its molecular interaction studies with aspirin and cefazolin

Human serum paraoxonase1 (HuPON1) belongs to the family of A-esterases (EC.3.1.8.1). It is associated with HDL particle and prevents atherosclerosis by cleaving lipid hydroperoxides and other proatherogenic molecules of oxidized low density lipoproteins (LDL). Since the precise structure of HuPON1 is not yet available, the structure-function relationship between HuPON1 and activators/inhibitors is still unknown. Therefore, a theoretical model of HuPON1 was generated using homology modelling and precise molecular interactions of an activator aspirin and an inhibitor cefazolin with PON1 were studied using Autodock software. The ligand binding residues were found to be similar to the predicted active site residues. Both cefazolin and aspirin were found to dock in the vicinity of the predicted active sites of PON1; cefazolin bound at residues N166, S193 and Y71, while aspirin at residues N309, I310 and L311. Binding region in the PON1 by prediction (3D2GO server) and docking studies provide useful insight into mechanism of substrate and inhibitor binding to the enzyme active site.     

Homology modeling and active-site residues probing of the thermophilic Alicyclobacillus acidocaldarius esterase 2.

The moderate thermophilic eubacterium Alicyclobacillus (formerly Bacillus) acidocaldarius expresses a thermostable carboxylesterase (esterase 2) belonging to the hormone-sensitive lipase (HSL)-like group of the esterase/lipase family. Based on secondary structures predictions and a secondary structure-driven multiple sequence alignment with remote homologous protein of known three-dimensional (3D) structure, we previously hypothesized for this enzyme the alpha/beta-hydrolase fold typical of several lipases and esterases and identified Ser155, Asp252, and His282 as the putative members of the catalytic triad. In this paper we report the construction of a 3D model for this enzyme based on the structure of mouse acetylcholinesterase complexed with fasciculin. The model reveals the topological organization of the fold corroborating our predictions. As regarding the active-site residues, Ser155, Asp252, and His282 are located close to each other at hydrogen bond distances. Their catalytic role was here probed by biochemical and mutagenic studies. Moreover, on the basis of the secondary structure-driven multiple sequence alignment and the 3D structural model, a residue supposed important for catalysis, Gly84, was mutated to Ser. The activity of the mutated enzyme was drastically reduced. We propose that Gly84 is part of a putative "oxyanion hole" involved in the stabilization of the transition state similar to the C group of the esterase/lipase family.     

Homology modeling of Mycoplasma pneumoniae enolase and its molecular interaction with human plasminogen

Alpha (α)-enolase (e), a glycolytic enzyme, has an alternative role as a surface receptor of several bacteria mediating plasminogen (pg) binding. It is also recognized as a virulence factor of some pathogenic bacteria facilitating plasminogen activation and host cell invasion. A mycoplasmal α-enolase is also a plasminogen binding protein. Molecular interactions of enolase from Mycoplasma pneumoniae with host plasminogen would be useful for exploring the pathogen-host interaction. In an attempt to identify plasminogen binding sites of M. pneumoniae enolase, homology modeling and docking studies were conducted to obtain modeled structures of the M. pneumoniae enolase-plasminogen complex. The refined model was validated further by standard methods. Molecular docking revealed hydrogen bonding of eLys70-pgTyr50, eAsn165-pgThr66, eAla168-pgGlu21, eAsp17-pgLys70, and eAsn213-pgPro68/pgAsn69. Substantial decreases in accessible surface area (ASA) were observed and in concurrence with hydrogen bond pattern. These findings provide a detailed prediction of key residues that interact at the protein-protein interface. Our theoretical prediction is consistent with known biochemical data. The predicted interaction complex can be of great assistance in understanding structural insights, which is necessary to pathogen and host-component interaction. The ability of M. pneumoniae enolase to bind plasminogen may be indicative of an important role in invasion of this pathogen to host.     

Homology modeling of human Toll‐like receptors TLR7, 8, and 9 ligand‐binding domains

Toll‐like receptors (TLRs) play a key role in the innate immune system. The TLR7, 8, and 9 compose a family of intracellularly localized TLRs that signal in response to pathogen‐derived nucleic acids. So far, there are no crystallographic structures for TLR7, 8, and 9. For this reason, their ligand‐binding mechanisms are poorly understood. To enable first predictions of the receptor–ligand interaction sites, we developed three‐dimensional structures for the leucine‐rich repeat ectodomains of human TLR7, 8, and 9 based on homology modeling. To achieve a high sequence similarity between targets and templates, structural segments from all known TLR ectodomain structures (human TLR1/2/3/4 and mouse TLR3/4) were used as candidate templates for the modeling. The resulting models support previously reported essential ligand‐binding residues. They also provide a basis to identify three potential receptor dimerization mechanisms. Additionally, potential ligand‐binding residues are identified using combined procedures. We suggest further investigations of these residues through mutation experiments. Our modeling approach can be extended to other members of the TLR family or other repetitive proteins.     

肺炎链球菌双组份系统中的组氨酸激酶(YycG)的同源模建与分析

对肺炎链球菌双组份系统中的组氨酸激酶YycG进行同源模建, 并分析其与底物ADP的相互作用, 为寻找特异性的激酶抑制剂提供了理论依据。采用同源模建的方法构建YycG蛋白的三维结构, 并用ProCheck、Profile_3D软件对此结构模型的合理性进行验证; 用Autodock4.0软件将结构模型与ADP进行自动对接, 分析二者之间的相互作用。序列比对结果显示肺炎链球菌YycG蛋白与Thermotoga maritima X-ray晶体结构序列的同一性达33%; YycG模建后的结构与模板能很好的叠合; 在活性口袋处的保守的氨基酸残基Asn145、Asn149、Lys152以及口袋内部的疏水残基在结合、水解底物ADP的过程中发挥重要作用。组氨酸激酶YycG的模建合理, 该结构模型可作为设计抗菌药的研究起点。     

胰岛素受体底物PH结构域相互作用蛋白结构和功能研究进展

PHIP是一种与胰腺β细胞中胰岛素受体底物（IRS）的PH结构域相互作用的蛋白。根据小鼠PHIP（mPHIP）mRNA翻译的不同起始位点,除全长的PHIP1外,mPHIP基因还编码其他3种不同变异体。在胰岛素诱导的信号途径中,主要分布于细胞核的PHIP1和IRS-1的PH结构域相互作用,介导IRS蛋白酪氨酸的磷酸化。IRS-2和PHIP1的共表达能诱导IRS在细胞膜上的定位,促进葡萄糖转运蛋白4（GLUT4）向细胞质膜的转移。PHIP1的表达能提高β-细胞内细胞周期蛋白D2的表达,促进β细胞的生长。PHIP1的表达活化蛋白激酶B（PKB）,活化的PKB能明显抑制β细胞的凋亡。PHIP与胰岛素信号传导途径中其他信号分子的相互作用机制尚不明确。     

Homology modeling and assigned functional annotation of an uncharacterized antitoxin protein from Streptomyces xinghaiensis

Streptomyces xinghaiensis is a Gram-positive, aerobic and non-motile bacterium. The bacterial genome is known. Therefore, it is ofinterest to study the uncharacterized proteins in the genome. An uncharacterized protein (gi|518540893|86 residues) in thegenome was selected for a comprehensive computational sequence-structure-function analysis using available data and tools. Subcellularlocalization of the targeted protein with conserved residues and assigned secondary structures is documented. Sequencehomology search against the protein data bank (PDB) and non-redundant GenBank proteins using BLASTp showed differenthomologous proteins with known antitoxin function. A homology model of the target protein was developed using a knowntemplate (PDB ID: 3CTO:A) with 62% sequence similarity in HHpred after assessment using programs PROCHECK and QMEAN6.The predicted active site using CASTp is analyzed for assigned anti-toxin function. This information finds specific utility inannotating the said uncharacterized protein in the bacterial genome.     

Homology modeling of histidine-containing phosphocarrier protein and eosinophil-derived neurotoxin: Construction of models and comparison with experiment

Homology modeling methods have been used to construct models of two proteins—the histidine-containing phosphocarrier protein (HPr) from Mycoplasma capricolum and human eosinophil-derived neurotoxin (EDN). Comparison of the models with the subsequently determined X-ray crystal structures indicates that the core regions of both proteins are reasonably well reproduced, although the template structures are closer to the X-ray structures in these regions—possible enhancements are discussed. The conformations of most of the side chains in the core of HPr are well reproduced in the modeled structure. As expected, the conformations of surface side chains in this protein differ significantly from the X-ray structure. The loop regions of EDN were incorrectly modeled—reasons for this and possible enhancements are discussed. © 1995 Wiley-Liss, Inc.     

S-adenosylhomocysteine hydrolase from the archaeon Pyrococcus furiosus: biochemical characterization and analysis of protein structure by comparative molecular modeling

S-adenosylhomocysteine hydrolase (AdoHcyHD) is an ubiquitous enzyme that catalyzes the breakdown of S-adenosylhomocysteine, a powerful inhibitor of most transmethylation reactions, to adenosine and L-homocysteine. AdoHcyHD from the hyperthermophilic archaeon Pyrococcus furiosus (PfAdoHcyHD) was cloned, expressed in Escherichia coli, and purified. The enzyme is thermoactive with an optimum temperature of 95 degrees C, and thermostable retaining 100% residual activity after 1 h at 90 degrees C and showing an apparent melting temperature of 98 degrees C. The enzyme is a homotetramer of 190 kDa and contains four cysteine residues per subunit. Thiol groups are not involved in the catalytic process whereas disulfide bond(s) could be present since incubation with 0.8 M dithiothreitol reduces enzyme activity. Multiple sequence alignment of hyperthermophilic AdoHcyHD reveals the presence of two cysteine residues in the N-terminus of the enzyme conserved only in members of Pyrococcus species, and shows that hyperthermophilic AdoHcyHD lack eight C-terminal residues, thought to be important for structural and functional properties of the eukaryotic enzyme. The homology-modeled structure of PfAdoHcyHD shows that Trp220, Tyr181, Tyr184, and Leu185 of each subunit and Ile244 from a different subunit form a network of hydrophobic and aromatic interactions in the central channel formed at the subunits interface. These contacts partially replace the interactions of the C-terminal tail of the eukaryotic enzyme required for tetramer stability. Moreover, Cys221 and Lys245 substitute for Thr430 and Lys426, respectively, of the human enzyme in NAD-binding. Interestingly, all these residues are fairly well conserved in hyperthermophilic AdoHcyHDs but not in mesophilic ones, thus suggesting a common adaptation mechanism at high temperatures.