FUGUE: sequence-structure homology recognition using environment-specific substitution tables and structure-dependent gap penalties.

FUGUE, a program for recognizing distant homologues by sequence-structure comparison (http://www-cryst.bioc.cam.ac.uk/fugue/), has three key features. (1) Improved environment-specific substitution tables. Substitutions of an amino acid in a protein structure are constrained by its local structural environment, which can be defined in terms of secondary structure, solvent accessibility, and hydrogen bonding status. The environment-specific substitution tables have been derived from structural alignments in the HOMSTRAD database (http://www-cryst.bioc. cam.ac.uk/homstrad/). (2) Automatic selection of alignment algorithm with detailed structure-dependent gap penalties. FUGUE uses the global-local algorithm to align a sequence-structure pair when they greatly differ in length and uses the global algorithm in other cases. The gap penalty at each position of the structure is determined according to its solvent accessibility, its position relative to the secondary structure elements (SSEs) and the conservation of the SSEs. (3) Combined information from both multiple sequences and multiple structures. FUGUE is designed to align multiple sequences against multiple structures to enrich the conservation/variation information. We demonstrate that the combination of these three key features implemented in FUGUE improves both homology recognition performance and alignment accuracy.     

SCORE: predicting the core of protein models

MOTIVATION: The prediction of the regions of homology models that can be 'restrained by' or 'copied from' the basis structures is a vital step in correct model generation, because these regions are the models most accurate part. However, there is no ideal method for the identification of their limits. In most algorithms their length depends on the number of family members and definitions of secondary structure. RESULTS: The algorithm SCORE steps away from the conventional definitions of the core to identify from large numbers of basis structures those regions that can be considered structurally related to a target sequence. The use of phi, psi constraints to accurately pinpoint the regions that are conserved across a family and environmentally constrained substitution tables to extend these regions allows SCORE to rapidly (generally in under 1 s, an order of magnitude faster than methods such as MODELLER) identify and build the core of homology models from the alignments of the target sequence to the basis structures. The SCORE algorithm was used to build 114 model cores. In only two cases was the core size less than 50% of the structure and all the cores built had an RMSD of 3.7 A or less to the target structure.     

Incorporating knowledge-based biases into an energy-based side-chain modeling method: application to comparative modeling of protein structure

The performance of the self-consistent mean field theory (SCMFT) method for side-chain modeling, employing rotamer energies calculated with the flexible rotamer model (FRM), is evaluated in the context of comparative modeling of protein structure. Predictions were carried out on a test set of 56 model backbones of varying accuracy, to allow side-chain prediction accuracy to be analyzed as a function of backbone accuracy. A progressive decrease in the accuracy of prediction was observed as backbone accuracy decreased. However, even for very low backbone accuracy, prediction was substantially higher than random, indicating that the FRM can, in part, compensate for the errors in the modeled tertiary environment. It was also investigated whether the introduction in the FRM-SCMFT method of knowledge-based biases, derived from a backbone-dependent rotamer library, could enhance its performance. A bias derived from the backbone-dependent rotamer conformations alone did not improve prediction accuracy. However, a bias derived from the backbone-dependent rotamer probabilities improved prediction accuracy considerably. This bias was incorporated through two different strategies. In one (the indirect strategy), rotamer probabilities were used to reject unlikely rotamers a priori, thus restricting prediction by FRM-SCMFT to a subset containing only the most probable rotamers in the library. In the other (the direct strategy), rotamer energies were transformed into pseudo-energies that were added to the average potential energies of the respective rotamers, thereby creating hybrid energy-based/knowledge-based average rotamer energies, which were used by the FRM-SCMFT method for prediction. For all degrees of backbone accuracy, an optimal strength of the knowledge-based bias existed for both strategies for which predictions were more accurate than pure energy-based predictions, and also than pure knowledge-based predictions. Hybrid knowledge-based/energy-based methods were obtained from both strategies and compared with the SCWRL method, a hybrid method based on the same backbone-dependent rotamer library. The accuracy of the indirect method was approximately the same as that of the SCWRL method, but that of the direct method was significantly higher.     

Evolution of floral meristem identity genes. Analysis of Lolium temulentum genes related to APETALA1 and LEAFY of Arabidopsis

Flowering (inflorescence formation) of the grass Lolium temulentum is strictly regulated, occurring rapidly on exposure to a single long day (LD). During floral induction, L. temulentum differs significantly from dicot species such as Arabidopsis in the expression, at the shoot apex, of two APETALA1 (AP1)-like genes, LtMADS1 and LtMADS2, and of L. temulentum LEAFY (LtLFY). As shown by in situ hybridization, LtMADS1 and LtMADS2 are expressed in the vegetative shoot apical meristem, but expression increases strongly within 30 h of LD floral induction. Later in floral development, LtMADS1 and LtMADS2 are expressed within spikelet and floret meristems and in the glume and lemma primordia. It is interesting that LtLFY is detected quite late (about 12 d after LD induction) within the spikelet meristems, glumes, and lemma primordia. These patterns contrast with Arabidopsis, where LFY and AP1 are consecutively activated early during flower formation. LtMADS2, when expressed in transgenic Arabidopsis plants under the control of the AP1 promoter, could partially complement the organ number defect of the severe ap1-15 mutant allele, confirming a close relationship between LtMADS2 and AP1.     

X-ray structure of human stromelysin catalytic domain complexed with nonpeptide inhibitors: implications for inhibitor selectivity.

Effective inhibitors of matrix metalloproteinases (MMPs), a family of connective tissue-degrading enzymes, could be useful for the treatment of diseases such as cancer, multiple sclerosis, and arthritis. Many of the known MMP inhibitors are derived from peptide substrates, with high potency in vitro but little selectivity among MMPs and poor bioavailability. We have discovered nonpeptidic MMP inhibitors with improved properties, and report here the crystal structures of human stromelysin-1 catalytic domain (SCD) complexed with four of these inhibitors. The structures were determined and refined at resolutions ranging from 1.64 to 2.0 A. Each inhibitor binds in the active site of SCD such that a bulky diphenyl piperidine moiety penetrates a deep, predominantly hydrophobic S'1 pocket. The active site structure of the SCD is similar in all four inhibitor complexes, but differs substantially from the peptide hydroxamate complex, which has a smaller side chain bound in the S'1 pocket. The largest differences occur in the loop forming the "top" of this pocket. The occupation of these nonpeptidic inhibitors in the S'1 pocket provides a structural basis to explain their selectivity among MMPs. An analysis of the unique binding mode predicts structural modifications to design improved MMP inhibitors.     

Analysis of interactive packing of secondary structural elements in alpha/beta units in proteins

An alpha-helix and a beta-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical alpha/beta units in protein structures, defined at < or = 2.5 A resolution, shows that the interaxial distance between the alpha-helix and the beta-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The beta-sheet unit (beta u), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded beta-strands, packing with an alpha-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in alpha/beta units using the interacting residue information of equivalent alpha/beta units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.     

蛋白质序列和结构的保守性与其功能的关系

以蛋白质酪氨酸磷酸酶为模型,通过序列和结构的比较,对蛋白质序列、结构的保守性与其和进化的关系问题进行了研究。结果显示,虽然在酪氨酸磷酸酶超家族分子的序列中,仅有３个与其催化功能密切相关的残基是高度保守的,但它们功能结构域的核心拓扑结构却明显类似,其中存在着βａβ和βａβａ。２个保守的结构单元;此外,它们活性位点的拓扑结构也极其相似,因此,在保持蛋白南功能方面具有重要作用的残基是高度保守的,而具维持     

The insulin receptor: from protein sequence to structure

Sequences of the insulin receptor (IR), the type-I insulin-like growth-factor receptor (IGFR) and the insulin-receptor-related receptor show that they belong to a homologous family but, until recently, have given few clues about their structures. Three repeats of fibronectin type III have been identified close to the membrane. Although the N-terminal L1, Cys-rich and L2 domains of the IGFR have been identified from their sequences and their structures determined by X-ray crystallography, little is known of their relative positions in the complete receptor in vivo. Here, we ask what can be learnt further from the analysis of sequences, about the structure, organization and function of the extracellular regions of the IR family.     

Homozygosity mapping of the Achromatopsia locus in the Pingelapese.

Achromatopsia, or total color blindness (also referred to as "rod monochromacy"), is a severe retinal disorder characterized clinically by an inability to distinguish colors, impaired visual acuity in daylight, photophobia, and nystagmus. Inherited as an autosomal recessive trait, achromatopsia is rare in the general population (1:20,000-1:50,000). Among the Pingelapese people of the Eastern Caroline Islands, however, the disorder occurs at an extremely high frequency, as recounted in Oliver Sacks's popular book The Island of the Colorblind: 4%-10% of this island population have the disorder and approximately 30% carry the gene. This extraordinary enrichment of the disease allele most likely resulted from a sharp reduction in population in the late 18th century, in the aftermath of a typhoon and subsequent geographic and cultural isolation. To obtain insights into the genetic basis of achromatopsia, as well as into the genetic history of this region of Micronesia, a genomewide search for linkage was performed in three Pingelapese kindreds with achromatopsia. A two-step search was used with a DNA pooling strategy, followed by genotyping of individual family members. Genetic markers that displayed a shift toward homozygosity in the affected DNA pool were used to genotype individual members of the kindreds, and an achromatopsia locus was identified on 8q21-q22. A maximal multipoint LOD score of 9.5 was observed with marker D8S1707. Homozygosity was seen for three adjacent markers (D8S275, D8S1119, and D8S1707), whereas recombination was observed with the flanking markers D8S1757 and D8S270, defining the outer boundaries of the disease-gene locus that spans a distance of <6.5cM.