An evolutionary treasure: unification of a broad set of amidohydrolases related to urease

The recent determination of the three-dimensional structure of urease revealed striking similarities of enzyme architecture to adenosine deaminase and phosphotriesterase, evidence of a distant evolutionary relationship that had gone undetected by one-dimensional sequence comparisons. Here, based on an analysis of conservation patterns in three dimensions, we report the discovery of the same active-site architecture in an even larger set of enzymes involved primarily in nucleotide metabolism. As a consequence, we predict the three-dimensional fold and details of the active site architecture for dihydroorotases, allantoinases, hydantoinases, AMP-, adenine and cytosine deaminases, imidazolonepropionase, aryldialkylphosphatase, chlorohydrolases, formylmethanofuran dehydrogenases, and proteins involved in animal neuronal development. Two member families are common to archaea, eubacteria, and eukaryota. Thirteen other functions supported by the same structural motif and conserved chemical mechanism apparently represent later adaptations for different substrate specificities in different cellular contexts. © 1997 Wiley-Liss Inc.     

Analysis of six protein structures predicted by comparative modeling techniques

The protein structures of six comparative modeling targets were predicted in a procedure that relied on improved energy minimization, without empirical rules, to position all new atoms. The structures of human nucleoside diphosphate kinase NM23-H2, HPr from Mycoplasma capricolum, 2Fe-2S ferredoxin from Haloarcula marismortui, eosinophil-derived neurotoxin (EDN), mouse cellular retinoic acid protein I (CRABP1), and P450eryf were predicted with root mean square deviations on Cα atoms of 0.69, 0.73, 1.11, 1.48, 1.69, and 1.73 Å, respectively, compared to the target crystal structures. These differences increased as the sequence similarity between the target and parent proteins decreased from about 60 to 20% identity. More residues were predicted than form the common region shared by the two crystal structures. In most cases insertions or deletions between the target and the related protein of known structure were not correctly positioned. One two residue insertion in CRABP1 was predicted in the correct conformation, while a nine residue insertion in EDN was predicted in the correct spatial region, although not in the correct conformation. The positions of common cofactors and their binding sites were predicted correctly, even when overall sequence similarity was low. © 1995 Wiley-Liss, Inc.     

Effect of natural polymorphism on structure and function of the Yersinia pestis outer membrane porin F (OmpF protein): a computational study

The Yersinia pestis outer membrane porin F (OmpF) is a transmembrane protein located in the outer membrane of this Gram-negative bacterium which is the causative agent of plague, where it plays a significant role in controlling the selective permeability of the membrane. The amino acid sequences of OmpF proteins from 48 Y. pestis strains representing all currently available phylogenetic groups of this Gram-negative bacterium were recently deduced. Comparison of these amino acid sequences revealed that the OmpF can be present in four isoforms, the pestis-pestis type, and the pestis-microtus types I, II, and III. OmpF of the most recent pestis-pestis type has an alanine residue at the position 148, where all the pestis-microtus types have threonine there (T148A polymorphism). The variability of different pestis-microtus types is caused by an additional polymorphism at the 193rd position, where the OmpFs of the pestis-microtus type II and type III have isoleucine-glycine (IG⁺₁₉₃) or isoleucine-glycine-isoleucine-glycine (IGIG⁺₁₉₃) insertions, respectively (IG⁺₁₉₃ and IGIG⁺₁₉₃ polymorphism). To investigate potential effects of these sequence polymorphisms on the structural properties of the OmpF protein, we conducted multi-level computational analysis of its isoforms. Analysis of the I-TASSER-generated 3D-models revealed that the Yersinia OmpF is very similar to other non-specific enterobacterial porins. The T148A polymorphism affected a loop located in the external vestibule of the OmpF channel, whereas IG⁺₁₉₃ and IGIG⁺₁₉₃ polymorphisms affected one of its β-strands. Our analysis also suggested that polymorphism has moderate effect on the predicted local intrinsic disorder predisposition of OmpF, but might have some functional implementations.     

胰岛素受体底物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与胰岛素信号传导途径中其他信号分子的相互作用机制尚不明确。     

A novel parameterization scheme for energy equations and its use to calculate the structure of protein molecules

A novel scheme for the parameterization of a type of “potential energy” function for protein molecules is introduced. The function is parameterized based on the known conformations of previously determined protein structures and their sequence similarity to a molecule whose conformation is to be calculated. Once parameterized, minima of the potential energy function can be located using a version of simulated annealing which has been previously shown to locate global and near-global minima with the given functional form. As a test problem, the potential was parameterized based on the known structures of the rubredoxins from Desulfovibrio vulgaris, Desulfovibrio desulfuricans, and Clostridium pasteurianum, which vary from 45 to 54 amino acids in length, and the sequence alignments of these molecules with the rubredoxin sequence from Desulfovibrio gigas. Since the Desulfovibrio gigas rubredeoxin conformation has also been determined, it is possible to check the accuracy of the results. Ten simulated-annealing runs from random starting conformations were performed. Seven of the 10 resultant conformations have an all-C_α rms deviation from the crystallographically determined conformation of less than 1.7 Å. For five of the structures, the rms deviation is less than 0.8 Å. Four of the structures have conformations which are virtually identical to each other except for the position of the carboxy-terminal residue. This is also the conformation which is achieved if the determined crystal structure is minimized with the same potential. The all-C_α rms difference between the crystal and minimized crystal structures is 0.6 Å. It is further observed that the “energies” of the structures according to the potential function exhibit a strong correlation with rms deviation from the native structure. The conformations of the individual model structures and the computational aspects of the modeling procedure are discussed. © 1993 Wiley-Liss, Inc.     

Discovery of novel inhibitors of Mycobacterium tuberculosis MurG: homology modelling,structure based pharmacophore,molecular docking,and molecular dynamics simulations

MurG (Rv2153c) is a key player in the biosynthesis of the peptidoglycan layer in Mycobacterium tuberculosis (Mtb). This work is an attempt to highlight the structural and functional relationship of Mtb MurG, the three-dimensional (3D) structure of protein was constructed by homology modelling using Discovery Studio 3.5 software. The quality and consistency of generated model was assessed by PROCHECK, ProSA and ERRAT. Later, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with substrate Uridine-diphosphate-N-acetylglucosamine (UD1) facilitated us to employ structure-based virtual screening approach to obtain new hits from Asinex database using energy-optimized pharmacophore modelling (e-pharmacophore). The pharmacophore model was validated using enrichment calculations, and finally, validated model was employed for high-throughput virtual screening and molecular docking to identify novel Mtb MurG inhibitors. This study led to the identification of 10 potential compounds with good fitness, docking score, which make important interactions with the protein active site. The 25 ns MD simulations of three potential lead compounds with protein confirmed that the structure was stable and make several non-bonding interactions with amino acids, such as Leu290, Met310 and Asn167. Hence, we concluded that the identified compounds may act as new leads for the design of Mtb MurG inhibitors.     

Refinement of protein structures by iterative comparative modeling and CryoEM density fitting

We developed a method for structure characterization of assembly components by iterative comparative protein structure modeling and fitting into cryo-electron microscopy (cryoEM) density maps. Specifically, we calculate a comparative model of a given component by considering many alternative alignments between the target sequence and a related template structure while optimizing the fit of a model into the corresponding density map. The method relies on the previously developed Moulder protocol that iterates over alignment, model building, and model assessment. The protocol was benchmarked using 20 varied target-template pairs of known structures with less than 30% sequence identity and corresponding simulated density maps at resolutions from 5A to 25A. Relative to the models based on the best existing sequence profile alignment methods, the percentage of C(alpha) atoms that are within 5A of the corresponding C(alpha) atoms in the superposed native structure increases on average from 52% to 66%, which is half-way between the starting models and the models from the best possible alignments (82%). The test also reveals that despite the improvements in the accuracy of the fitness function, this function is still the bottleneck in reducing the remaining errors. To demonstrate the usefulness of the protocol, we applied it to the upper domain of the P8 capsid protein of rice dwarf virus that has been studied by cryoEM at 6.8A. The C(alpha) root-mean-square deviation of the model based on the remotely related template, bluetongue virus VP7, improved from 8.7A to 6.0A, while the best possible model has a C(alpha) RMSD value of 5.3A. Moreover, the resulting model fits better into the cryoEM density map than the initial template structure. The method is being implemented in our program MODELLER for protein structure modeling by satisfaction of spatial restraints and will be applicable to the rapidly increasing number of cryoEM density maps of macromolecular assemblies.