			Home About us Contact

Homology Modeling (homology + modeling)

Distribution by Scientific Domains

Chemistry	73%
Life Sciences	17%
2 Other Domains	10%

Selected Abstracts

Exploring the binding site of the human muscarinic M3 receptor: Homology modeling and docking study

INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY, Issue 8 2007
Liliana Ostopovici
Abstract The human muscarinic M3 receptor (hM3) and its interactions with selective agonists and antagonists were investigated by means of combined homology and docking approach. Also, two pharmacophoric models for the hM3 agonist and antagonist binding sites were proposed. The three-dimensional (3D) structure of hM3 receptor was modeled based on the high-resolution X-ray structure of bovine rhodopsin from the Protein Data Bank (PDB). To validate the reliability of the model obtained, the main chain torsion angles phi (,) and psi (,) were examined in a Ramachandran plot, and all omega angles were measured for peptidic bond planarity. The characteristics of the active site, the position, and the orientation of ligands in situ, as well as the binding modes of the representative agonists and antagonists, were analyzed by applying a molecular docking technique using the AutoDock 3.0.5 program. Specific interactions responsible for recognition of the hM3 receptor, like ionic bond formed between protonated amine of the ligands and the Asp3.6 side chain were identified. Structure,reactivity relationships have been explained by analyzing the 3D structure of the hM3 model and the ligand conformations resulted from molecular docking simulation. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007 [source]

Homology modeling and molecular dynamics simulations of lymphotactin

PROTEIN SCIENCE, Issue 11 2000
Buyong Ma
Abstract We have modeled the structure of human lymphotactin (hLpnt), by homology modeling and molecular dynamics simulations. This chemokine is unique in having a single disulfide bond and a long C-terminal tail. Because other structural classes of chemokines have two pairs of Cys residues, compared to one in Lpnt, and because it has been shown that both disulfide bonds are required for stability and function, the question arises how the Lpnt maintains its structural integrity. The initial structure of hLpnt was constructed by homology modeling. The first 63 residues in the monomer of hLpnt were modeled using the structure of the human CC chemokine, RANTES, whose sequence appeared most similar. The structure of the long C-terminal tail, missing in RANTES, was taken from the human muscle fatty-acid binding protein. In a Protein Data Bank search, this protein was found to contain a sequence that was most homologous to the long tail. Consequently, the modeled hLpnt C-terminal tail consisted of both ,-helical and ,-motifs. The complete model of the hLpnt monomer consisted of two ,-helices located above the five-stranded ,-sheet. Molecular dynamics simulations of the solvated initial model have indicated that the stability of the predicted fold is related to the geometry of Pro78. The five-stranded ,-sheet appeared to be preserved only when Pro78 was modeled in the cis conformation. Simulations were also performed both for the C-terminal truncated forms of the hLpnt that contained one or two (CC chemokine-like) disulfide bonds, and for the chicken Lpnt (cLpnt). Our MD simulations indicated that the turn region (T30-G34) in hLpnt is important for the interactions with the receptor, and that the long C-terminal region stabilizes both the turn (T30-G34) and the five-stranded ,-sheet. The major conclusion from our theoretical studies is that the lack of one disulfide bond and the extension of the C-terminus in hLptn are mutually complementary. It is very likely that removal of two Cys residues sufficiently destabilizes the structure of a chemokine molecule, particularly the core ,-sheet, to abolish its biological function. However, this situation is rectified by the long C-terminal segment. The role of this long region is most likely to stabilize the first ,-turn region and ,-helix H1, explaining how this chemokine can function with a single disulfide bond. [source]

Trifluoroethanol and binding to model membranes stabilize a predicted turn in a peptide corresponding to the first extracellular loop of the angiotensin II AT1A receptor

BIOPOLYMERS, Issue 1 2002
Roberto K. Salinas
Abstract Homology modeling of the angiotensin II AT1A receptor based on rhodopsin,s crystal structure has assigned the 92,100 (YRWPFGNHL) sequence of the receptor to its first extracellular loop. Solution and membrane-bound conformational properties of a peptide containing this sequence (EL1) were examined by CD, fluorescence, and 1H-NMR. CD spectra in aqueous solution revealed an equilibrium between less organized and folded conformers. NMR spectra indicated the coexistence of trans and cis isomers of the Trp3,Pro4 bond. A positive band at 226 nm in the CD spectra suggested aromatic ring stacking, modulated by EL1's ionization degree. CD spectra showed that trifluoroethanol (TFE), or binding to detergent micelles and phospholipid bilayers, shifted the equilibrium toward conformers with higher secondary structure content. Different media gave rise to spectra suggestive of different ,-turns. Chemical shift changes in the NMR spectra corroborated the stabilization of different conformations. Thus, environments of lower polarity or binding to interfaces probably favored the formation of hydrogen bonds, stabilizing ,-turns, predicted for this sequence in the whole receptor. Increases in Trp3 fluorescence intensity and anisotropy, blue shifts of the maximum emission wavelength, and pK changes also evinced the interaction between EL1 and model membranes. Binding was seen to depend on both hydrophobic and electrostatic interactions, as well as lipid phase packing. Studies with water-soluble and membrane-bound fluorescence quenchers demonstrated that Trp3 is located close to the water,membrane interface. The results are discussed with regard to possible implications in receptor folding and function. © 2002 Wiley Periodicals, Inc. Biopolymers 65: 21,31, 2002 [source]

Improvement of low-temperature caseinolytic activity of a thermophilic subtilase by directed evolution and site-directed mutagenesis

BIOTECHNOLOGY & BIOENGINEERING, Issue 5 2009
Chuan-Qi Zhong
Abstract By directed evolution and subsequent site-directed mutagenesis, cold-adapted variants of WF146 protease, a thermophilic subtilase, have been successfully engineered. A four-amino acid substitution variant RTN29 displayed a sixfold increase in caseinolytic activity in the temperature range of 15,25°C, a down-shift of optimum temperature by ,15°C, as well as a decrease in thermostability, indicating it follows the general principle of trade-off between activity and stability. Nevertheless, to some extent RTN29 remained its thermophilic nature, and no loss of activity was observed after heat-treatment at 60°C for 2,h. Notably, RTN29 exhibited a lower hydrolytic activity toward suc-AAPF-pNA, due to an increase in Km and a decrease in kcat, in contrast to other artificially cold-adapted subtilases with increased low-temperature activity toward small synthetic substrates. All mutations (S100P, G108S, D114G, M137T, T153A, and S246N) identified in the cold-adapted variants occurred within or near the substrate-binding region. None of these mutations, however, match the corresponding sites in naturally psychrophilic and other artificially cold-adapted subtilases, implying there are multiple routes to cold adaptation. Homology modeling and structural analysis demonstrated that these mutations led to an increase in mobility of substrate-binding region and a modulation of substrate specificity, which seemed to account for the improvement of the enzyme's catalytic activity toward macromolecular substrates at lower temperatures. Our study may provide valuable information needed to develop enzymes coupling high stability and high low-temperature activity, which are highly desired for industrial use. Biotechnol. Bioeng. 2009; 104: 862,870. © 2009 Wiley Periodicals, Inc. [source]

Characterization of a hemocyte intracellular fatty acid-binding protein from crayfish (Pacifastacus leniusculus) and shrimp (Penaeus monodon)

FEBS JOURNAL, Issue 13 2006
Irene Söderhäll
Intracellular fatty acid-binding proteins (FABPs) are small members of the superfamily of lipid-binding proteins, which occur in invertebrates and vertebrates. Included in this superfamily are the cellular retinoic acid-binding proteins and retinol-binding proteins, which seem to be restricted to vertebrates. Here, we report the cDNA cloning and characterization of two FABPs from hemocytes of the freshwater crayfish Pacifastacus leniusculus and the shrimp Penaeus monodon. In both these proteins, the binding triad residues involved in interaction with ligand carboxylate groups are present. From the sequence and homology modeling, the proteins are probably FABPs and not retinoic acid-binding proteins. The crayfish transcript (plFABP) was detected at high level in hemocytes, hepatopancreas, intestine and ovary and at low level in hematopoietic tissue and testis. Its expression in hematopoietic cells varied depending on the state of the crayfish from which it was isolated. Expression was 10,15 times higher in cultures isolated from crayfish with red colored plasma, in which hemocyte synthesis was high, if retinoic acid was added to the culture medium. In normal colored crayfish, with normal levels of hemocytes, no increase in expression of p1FABP was detected. Two other putative plFABP ligands, stearic acid and oleic acid, did not have any effect on plFABP expression in hematopoietic cells. These results suggest that retinoic acid-dependent signaling may be present in crustaceans. [source]

The structural comparison of the bacterial PepX and human DPP-IV reveals sites for the design of inhibitors of PepX activity

FEBS JOURNAL, Issue 8 2005
Pascal Rigolet
X-prolyl dipeptidyl aminopeptidases (X-PDAP) are enzymes catalysing the release of dipeptides from the amino termini of polypeptides containing a proline or an alanine at the penultimate position. Involved in various mammalian regulation processes, as well as in chronic human diseases, they have been proposed to play a role in pathogenicity for Streptococci. We compared the structure of X-PDAP from Lactococcus lactis (PepX) with its human counterpart DPP-IV. Despite very different overall folds, the residues most implicated for X-PDAP activity are conserved in the same positions and orientations in both enzymes, thus defining a structural signature for the X-PDAP specificity that crosses the species frontiers of evolution. Starting from this observation, we tested some inhibitors of DPP-IV on PepX activity, for which no specific inhibitor is known. We thus found that PepX was highly sensitive to valine-pyrrolidide with a KI of 9.3 µm, close to that reported in DPP-IV inhibition. We finally used the structure of PepX from L. lactis as a template for computer-based homology modeling of PepX from the pathogenic Streptococcus gordonii. Docking simulations of valine-pyrrolidide into the active site of PepX led to the identification of key residues for a rational drug design against PepX from Streptococci. These results could have applications in human health giving new perspectives to the struggle against pathogens. [source]

Complete subunit sequences, structure and evolution of the 6 × 6-mer hemocyanin from the common house centipede, Scutigera coleoptrata

FEBS JOURNAL, Issue 13 2003
Kristina Kusche
Hemocyanins are large oligomeric copper-containing proteins that serve for the transport of oxygen in many arthropod species. While studied in detail in the Chelicerata and Crustacea, hemocyanins had long been considered unnecessary in the Myriapoda. Here we report the complete molecular structure of the hemocyanin from the common house centipede Scutigera coleoptrata (Myriapoda: Chilopoda), as deduced from 2D-gel electrophoresis, MALDI-TOF mass spectrometry, protein and cDNA sequencing, and homology modeling. This is the first myriapod hemocyanin to be fully sequenced, and allows the investigation of hemocyanin structure,function relationship and evolution. S. coleoptrata hemocyanin is a 6 × 6-mer composed of four distinct subunit types that occur in an approximate 2 : 2 : 1 : 1 ratio and are 49.5,55.5% identical. The cDNA of a fifth, highly diverged, putative hemocyanin was identified that is not included in the native 6 × 6-mer hemocyanin. Phylogenetic analyses show that myriapod hemocyanins are monophyletic, but at least three distinct subunit types evolved before the separation of the Chilopoda and Diplopoda more than 420 million years ago. In contrast to the situation in the Crustacea and Chelicerata, the substitution rates among the myriapod hemocyanin subunits are highly variable. Phylogenetic analyses do not support a common clade of Myriapoda and Hexapoda, whereas there is evidence in favor of monophyletic Mandibulata. [source]

Structure,activity relationship of the p55 TNF receptor death domain and its lymphoproliferation mutants

FEBS JOURNAL, Issue 5 2001
Gert De Wilde
Upon stimulation with tumor necrosis factor (TNF), the TNF receptor (TNFR55) mediates a multitude of effects both in normal and in tumor cells. Clustering of the intracellular domain of the receptor, the so-called death domain (DD), is responsible for both the initiation of cell killing and the activation of gene expression. To characterize this domain further, TNFR55 DD was expressed and purified as a thioredoxin fusion protein in Escherichia coli. Circular dichroism, steady-state and time-resolved fluorescence spectroscopy were used to compare TNFR55 DD with DDs of the Fas antigen (Fas), the Fas-associating protein with DD (FADD) and p75 nerve growth factor receptor, for which the 3-dimensional structure are already known. The structural information derived from the measurements strongly suggests that TNFR55 DD adopts a similar fold in solution. This prompted a homology modeling of the TNFR DD 3-D structure using FADD as a template. In vivo studies revealed a difference between the two lymphoproliferation (lpr) mutations. Biophysical techniques were used to analyze the effect of changing Leu351 to Ala and Leu351 to Asn on the global structure and its impact on the overall stability of TNFR55 DD. The results obtained from these experiments in combination with the modeled structure offer an explanation for the in vivo observed difference. [source]

Classic and atypical fibrodysplasia ossificans progressiva (FOP) phenotypes are caused by mutations in the bone morphogenetic protein (BMP) type I receptor ACVR1,

HUMAN MUTATION, Issue 3 2009
Frederick S. Kaplan
Abstract Fibrodysplasia ossificans progressiva (FOP) is an autosomal dominant human disorder of bone formation that causes developmental skeletal defects and extensive debilitating bone formation within soft connective tissues (heterotopic ossification) during childhood. All patients with classic clinical features of FOP (great toe malformations and progressive heterotopic ossification) have previously been found to carry the same heterozygous mutation (c.617G>A; p.R206H) in the glycine and serine residue (GS) activation domain of activin A type I receptor/activin-like kinase 2 (ACVR1/ALK2), a bone morphogenetic protein (BMP) type I receptor. Among patients with FOP-like heterotopic ossification and/or toe malformations, we identified patients with clinical features unusual for FOP. These atypical FOP patients form two classes: FOP-plus (classic defining features of FOP plus one or more atypical features) and FOP variants (major variations in one or both of the two classic defining features of FOP). All patients examined have heterozygous ACVR1 missense mutations in conserved amino acids. While the recurrent c.617G>A; p.R206H mutation was found in all cases of classic FOP and most cases of FOP-plus, novel ACVR1 mutations occur in the FOP variants and two cases of FOP-plus. Protein structure homology modeling predicts that each of the amino acid substitutions activates the ACVR1 protein to enhance receptor signaling. We observed genotype-phenotype correlation between some ACVR1 mutations and the age of onset of heterotopic ossification or on embryonic skeletal development. Hum Mutat 0, 1,12, 2008. © 2008 Wiley-Liss, Inc. [source]

Spectrum of molecular defects and mutation detection rate in patients with severe hemophilia A,

HUMAN MUTATION, Issue 3 2005
Nadja Bogdanova
Abstract Hemophilia A is the most frequently occurring X-linked bleeding disorder, affecting one to two out of 10,000 males worldwide. Various types of mutations in the F8 gene are causative for this condition. It is well known that the most common mutation in severely affected patients is the intron 22 inversion, which accounts for about 45% of cases with F8 residual activity of less than 1%. Therefore, the aim of the present study was to determine the spectrum and distribution of mutations in the F8 gene in a large group of patients with severe hemophilia A who previously tested negative for the common intron 22 inversion. Here we report on a mutation analysis of 86 patients collected under the above-mentioned criterion. The pathogenic molecular defect was identified in all patients, and thus our detection rate was virtually 100%. Thirty-four of the identified mutations are described for the first time. The newly detected amino acid substitutions were scored for potential gross or local conformational changes and influence on molecular stability for every single F8 domain with available structures, using homology modeling. Hum Mutat Res 26(3), 249,254, 2005. © 2005 Wiley-Liss, Inc. [source]

Mechanistic insights into oxidosqualene cyclizations through homology modeling

JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 6 2003
Gasch, Tanja Schulz
Abstract 2,3-Oxidosqualene cyclases (OSC) are key enzymes in sterol biosynthesis. They catalyze the stereoselective cyclization and skeletal rearrangement of (3S)-2,3-oxidosqualene to lanosterol in mammals and fungi and to cycloartenol in algae and higher plants. Sequence information and proposed mechanism of 2,3-oxidosqualene cyclases are closely related to those of squalene-hopene cyclases (SHC), which represent functional analogs of OSCs in bacteria. SHCs catalyze the cationic cyclization cascade converting the linear triterpene squalene to fused ring compounds called hopanoids. High stereoselectivity and precision of the skeletal rearrangements has aroused the interest of researchers for nearly half a century, and valuable data on studying mechanistic details in the complex enzyme-catalyzed cyclization cascade has been collected. Today, interest in cyclases is still unbroken, because OSCs became targets for the development of antifungal and hypocholesterolemic drugs. However, due to the large size and membrane-bound nature of OSCs, three-dimensional structural information is still not available, thus preventing a complete understanding of the atomic details of the catalytic mechanism. In this work, we discuss results gained from homology modeling of human OSC based on structural information of SHC from Alicyclobacillus acidocaldarius and propose a structural model of human OSC. The model is in accordance with previously performed experimental studies with mechanism-based suicide inhibitors and mutagenesis experiments with altered activity and product specificity. Structural insight should strongly stimulate structure-based design of antifungal or cholesterol-lowering drugs. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 741,753, 2003 [source]

Semiautomatic sequence-specific assignment of proteins based on the tertiary structure,The program st2nmr

JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 3 2002
Primo, Pristov
Abstract The sequence-specific assignment of resonances is still the most time-consuming procedure that is necessary as the first step in high-resolution NMR studies of proteins. In many cases a reliable three-dimensional (3D) structure of the protein is available, for example, from X-ray spectroscopy or homology modeling. Here we introduce the st2nmr program that uses the 3D structure and Nuclear Overhauser Effect spectroscopy (NOESY) peak list(s) to evaluate and optimize trial sequence-specific assignments of spin systems derived from correlation spectra to residues of the protein. A distance-dependent target function that scores trial assignments based on the presence of expected NOESY crosspeaks is optimized in a Monte Carlo fashion. The performance of the program st2nmr is tested on real NMR data of an ,-helical (cytochrome c) and ,-sheet (lipocalin) protein using homology models and/or X-ray structures; it succeeded in completely reproducing the correct sequence-specific assignments in most cases using 2D and/or 15N/13C Nuclear Overhauser Effect (NOE) data. Additionally to amino acid residues the program can also handle ligands that are bound to the protein, such as heme, and can be used as a complementary tool to fully automated assignment procedures. © 2002 Wiley Periodicals, Inc. J Comput Chem 23: 335,340, 2002 [source]

Evolution of physics-based methodology for exploring the conformational energy landscape of proteins

JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 1 2002
Harold A. Scheraga
Abstract The evolution of our physics-based computational methods for determining protein conformation without the introduction of secondary-structure predictions, homology modeling, threading, or fragment coupling is described. Initial use of a hard-sphere potential captured much of the structural properties of polypeptide chains, and subsequent more refined force fields, together with efficient methods of global optimization provide indications that progress is being made toward an understanding of the interresidue interactions that underlie protein folding. © 2002 Wiley Periodicals, Inc. J Comput Chem 23: 28,34, 2002 [source]

Analysis of interactions responsible for vasopressin binding to human neurohypophyseal hormone receptors,molecular dynamics study of the activated receptor,vasopressin,G, systems

JOURNAL OF PEPTIDE SCIENCE, Issue 3 2006
Magdalena J., lusarz
Abstract Vasopressin (CYFQNCPRG-NH2, AVP) is a semicyclic endogenous peptide, which exerts a variety of biological effects in mammals. The main physiological roles of AVP are the regulation of water balance and the control of blood pressure and adrenocorticotropin hormone (ACTH) secretion, mediated via three different subtypes of vasopressin receptors: V1a, V1b and V2 receptors (V1aR, V1bR and V2R, respectively). They are the members of the class A, G-protein-coupled receptors (GPCRs). AVP also modulates several behavioral and social functions. In this study, the interactions responsible for AVP binding to vasopressin V1a and V2 receptors versus the closely related oxytocin ([I3,L8]AVP, OT) receptor (OTR) have been investigated. Three-dimensional models of the activated receptors were constructed using multiple sequence alignment, followed by homology modeling using the complex of activated rhodopsin with Gt,C -terminal peptide of transducin MII-Gt(338-350) prototype as a template. AVP was docked into the receptor-G, systems. The three lowest-energy pairs of receptor-AVP-G, (two complexes per each receptor) were selected. The 1-ns unconstrained molecular dynamics (MD) of complexes embedded into the fully hydrated 1-palmitoyl-2-oleoyl- sn -glycero-3-phosphatidylcholine (POPC) lipid bilayer was conducted in the AMBER 7.0 force field. Six relaxed receptor-AVP-G, models were obtained. The residues responsible for AVP binding to vasopressin receptors have been identified and a different mechanism of AVP binding to V2R than to V1aR has been proposed. Copyright © 2005 European Peptide Society and John Wiley & Sons, Ltd. [source]

Molecular models of the procoagulant Factor VIIIa,Factor IXa complex

JOURNAL OF THROMBOSIS AND HAEMOSTASIS, Issue 9 2005
L. AUTIN
Summary.,Background:,Formation of the intrinsic tenase complex is an essential event in the procoagulant reactions that lead to clot formation. The tenase complex is formed when the activated serine protease, Factor IXa (FIXa), and its cofactor Factor VIIIa (FVIIIa) assemble on a phospholipid surface to proteolytically convert the zymogen Factor X (FX) into its active form FXa. The physiological relevance of the tenase complex is evident in hemophilia A or B patients who present with bleeding disorders. Objectives:,The purpose of this study was to establish three-dimensional (3D) models of the FVIIIa,FIXa complex. Methods:,First, we built two new theoretical models of FVIIIa via homology modeling, inter-domain docking and loop simulation algorithms as well as a model for FIXa. This was followed by pseudo-Brownian protein,protein docking in internal coordinates with the ICM (Internal Coordinates Mechanics) program between the two FVIIIa and the FIXa structures. Results:,Ten representative models of this complex are presented based on agreements with known experimental data and according to structural criteria. Conclusions:,These novel 3D models will help guide future site directed mutagenesis aimed at improving the functionality of FVIIIa and/or FIXa and will contribute to a better understanding of the role of this macromolecular complex in the blood coagulation cascade. [source]

Docking and homology modeling explain inhibition of the human vesicular glutamate transporters

PROTEIN SCIENCE, Issue 9 2007
Jonas Almqvist
Abstract As membrane transporter proteins, VGLUT1,3 mediate the uptake of glutamate into synaptic vesicles at presynaptic nerve terminals of excitatory neural cells. This function is crucial for exocytosis and the role of glutamate as the major excitatory neurotransmitter in the central nervous system. The three transporters, sharing 76% amino acid sequence identity in humans, are highly homologous but differ in regional expression in the brain. Although little is known regarding their three-dimensional structures, hydropathy analysis on these proteins predicts 12 transmembrane segments connected by loops, a topology similar to other members in the major facilitator superfamily, where VGLUT1,3 have been phylogenetically classified. In this work, we present a three-dimensional model for the human VGLUT1 protein based on its distant bacterial homolog in the same superfamily, the glycerol-3-phosphate transporter from Escherichia coli. This structural model, stable during molecular dynamics simulations in phospholipid bilayers solvated by water, reveals amino acid residues that face its pore and are likely to affect substrate translocation. Docking of VGLUT1 substrates to this pore localizes two different binding sites, to which inhibitors also bind with an overall trend in binding affinity that is in agreement with previously published experimental data. [source]

Diagnostic cross-linking of paired cysteine pairs demonstrates homologous structures for two chemoreceptor domains with low sequence identity

PROTEIN SCIENCE, Issue 1 2006
Wing-Cheung Lai
Abstract Hundreds of bacterial chemoreceptors from many species have periplasmic, ligand-recognition domains of approximately the same size, but little or no sequence identity. The only structure determined is for the periplasmic domain of chemoreceptor Tar from Salmonella and Escherichia coli. Do sequence-divergent but similarly sized chemoreceptor periplasmic domains have related structures? We addressed this issue for the periplasmic domain of chemoreceptor TrgE from E. coli, which has a low level of sequence similarity to Tar, by combining homology modeling and diagnostic cross-linking between pairs of introduced cysteines. A homology model of the TrgE domain was created using the homodimeric, four-helix bundle structure of the TarS domain from Salmonella. In this model, we chose four pairs of positions at which introduced cysteines would be sufficiently close to form disulfides across each of four different helical interfaces. For each pair we chose a second pair, in which one cysteine of the original pair was shifted by one position around the helix and thus would be less favorably placed for disulfide formation. We created genes coding for proteins containing four such pairs of cysteine pairs and investigated disulfide formation in vivo as well as functional consequences of the substitutions and disulfides between neighboring helices. Results of the experimental tests provided strong support for the accuracy of the model, indicating that the TrgE periplasmic domain is very similar to the TarS domain. Diagnostic cross-linking of paired pairs of introduced cysteines could be applied generally as a stringent test of homology models. [source]

Automatic generation and evaluation of sparse protein signatures for families of protein structural domains

PROTEIN SCIENCE, Issue 1 2005
Matthew J. Blades
Abstract We identified key residues from the structural alignment of families of protein domains from SCOP which we represented in the form of sparse protein signatures. A signature-generating algorithm (SigGen) was developed and used to automatically identify key residues based on several structural and sequence-based criteria. The capacity of the signatures to detect related sequences from the SWISSPROT database was assessed by receiver operator characteristic (ROC) analysis and jack-knife testing. Test signatures for families from each of the main SCOP classes are described in relation to the quality of the structural alignments, the SigGen parameters used, and their diagnostic performance. We show that automatically generated signatures are potently diagnostic for their family (ROC50 scores typically >0.8), consistently outperform random signatures, and can identify sequence relationships in the "twilight zone" of protein sequence similarity (<40%). Signatures based on 15%,30% of alignment positions occurred most frequently among the best-performing signatures. When alignment quality is poor, sparser signatures perform better, whereas signatures generated from higher-quality alignments of fewer structures require more positions to be diagnostic. Our validation of signatures from the Globin family shows that when sequences from the structural alignment are removed and new signatures generated, the omitted sequences are still detected. The positions highlighted by the signature often correspond (alignment specificity >0.7) to the key positions in the original (non-jack-knifed) alignment. We discuss potential applications of sparse signatures in sequence annotation and homology modeling. [source]

Weak alignment offers new NMR opportunities to study protein structure and dynamics

PROTEIN SCIENCE, Issue 1 2003
Ad Bax
Protein solution nuclear magnetic resonance (NMR) can be conducted in a slightly anisotropic environment, where the orientational distribution of the proteins is no longer random. In such an environment, the large one-bond internuclear dipolar interactions no longer average to zero and report on the average orientation of the corresponding vectors relative to the magnetic field. The desired very weak ordering, on the order of 10,3, can be induced conveniently by the use of aqueous nematic liquid crystalline suspensions or by anisotropically compressed hydrogels. The resulting residual dipolar interactions are scaled down by three orders of magnitude relative to their static values, but nevertheless can be measured at high accuracy. They are very precise reporters on the average orientation of bonds relative to the molecular alignment frame, and they can be used in a variety of ways to enrich our understanding of protein structure and function. Applications to date have focused primarily on validation of structures, determined by NMR, X-ray crystallography, or homology modeling, and on refinement of structures determined by conventional NMR approaches. Although de novo structure determination on the basis of dipolar couplings suffers from a severe multiple minimum problem, related to the degeneracy of dipolar coupling relative to inversion of the internuclear vector, a number of approaches can address this problem and potentially can accelerate the NMR structure determination process considerably. In favorable cases, where large numbers of dipolar couplings can be measured, inconsistency between measured values can report on internal motions. [source]

Homology modeling and molecular dynamics simulations of lymphotactin

A series of molecular dynamics and homology modeling computer labs for an undergraduate molecular modeling course

BIOCHEMISTRY AND MOLECULAR BIOLOGY EDUCATION, Issue 4 2010
Donald E. Elmore
Abstract As computational modeling plays an increasingly central role in biochemical research, it is important to provide students with exposure to common modeling methods in their undergraduate curriculum. This article describes a series of computer labs designed to introduce undergraduate students to energy minimization, molecular dynamics simulations, and homology modeling. These labs were created as part of a one-semester course on the molecular modeling of biochemical systems. Students who completed these activities felt that they were an effective component of the course, reporting improved comfort with the conceptual background and practical implementation of the computational methods. Although created as a component of a larger course, these activities could be readily adapted for a variety of other educational contexts. As well, all of these labs utilize software that is freely available in an academic environment and can be run on fairly common computer hardware, making them accessible to teaching environments without extensive computational resources. [source]

Endo- and exo-inulinases: Enzyme-substrate interaction and rational immobilization

BIOTECHNOLOGY PROGRESS, Issue 2 2010
Alessandra Basso
Abstract Three-dimensional models of exoinulinase from Bacillus stearothermophilus and endoinulinase from Aspergillus niger were built up by means of homology modeling. The crystal structure of exoinulinase from Aspergillus awamori was used as a template, which is the sole structure of inulinase resolved so far. Docking and molecular dynamics simulations were performed to investigate the differences between the two inulinases in terms of substrate selectivity. The analysis of the structural differences between the two inulinases provided the basis for the explanation of their different regio-selectivity and for the understanding of enzyme-substrate interactions. Surface analysis was performed to point out structural features that can affect the efficiency of enzymes also after immobilization. The computational analysis of the three-dimensional models proved to be an effective tool for acquiring information and allowed to formulate an optimal immobilized biocatalyst even more active that the native one, thus enabling the full exploitation of the catalytic potential of these enzymes. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2010 [source]