Exercise 6: Homology Modeling

Total Page:16

File Type:pdf, Size:1020Kb

Exercise 6: Homology Modeling Workshop in Computational Structural Biology (81813), Spring 2019 Exercise 6: Homology Modeling Homology modeling, sometimes referred to as comparative modeling, is a commonly-used alternative to ab-initio modeling (which requires zero knowledge about the native structure, but also overlooks much of our understanding of structural diversity). Homology modeling leverages the massive amount of structural knowledge available today to guide our search of the native structure. Practically, it assumes that sequence similarity implies structural similarity, i.e. assumes a protein homologous to our query protein cannot be very different from it. Therefore, we can roughly "thread" the query sequence onto the template (homologue) structure, and exploit existing knowledge to guide our modeling. Contents Homology modeling with Rosetta 1 Overview 1 Template Selection 1 Loop Modeling 2 Running the full simulation 4 Bonus section! Homology modeling with automated/pre-calculated servers 6 MODBASE 6 I-Tasser 6 Homology modeling with Rosetta Overview In this section, we will model the structure of the Erbin PDZ domain, using homologue template structures, and analyze the quality of the models. The structure of this protein has been solved (PDB ID 1MFG), so we will be able to compare prediction with experiment. As taught in class, classical homology modeling in Rosetta is composed of four phases: 1. Template selection 2. Template-query alignment 3. Rebuilding loop regions 4. Full-atom refinement of the model Our model protein is the Erbin PDZ domain: Erbin contains a class-I PDZ domain that binds to the C-terminal region of the receptor tyrosine kinase ErbB2. PDB ID 1MFG is the crystal structure of the human Erbin PDZ bound to the peptide EYLGLDVPV corresponding to the C-terminal residues of human ErbB2. ErbB2 is implicated in both signaling and oncogenicity. Template Selection The three-dimensional organization of the template determines to a great extent the quality of our result model. That is why it is such a critical step in homology modeling, and why sequence alignment and homology detection methods are so important in this context. Usually, successful modeling (i.e. generation of models within 3-5 Å RMSD of the experimentally solved structure) requires 20-30% sequence identity, and/or 30-50% sequence similarity to the template structure, given a reasonable alignment with a limited number of gaps. 1 Workshop in Computational Structural Biology (81813), Spring 2019 Question 1: Given two homologous structures and their sequence alignment, what structural regions are likely to be well and poorly aligned? Explain briefly. In the real world scenario, we have a sequence for which we want to know the structure. Here, we will first get the sequence from the structure that is already solved... The following commands extract chain A of 1MFG and renumber it to start from residue 1: $> prody fetch 1mfg $> prody select "protein and chain A" 1mfg -o tmp $> sequentialPdbResSeq.pl -pdbfile tmp.pdb > 1mfgA.pdb Next, extract the sequence of 1mfg chain A from the pdb file: $> getFastaFromCoords.pl -pdbfile 1mfgA.pdb -chain A > 1mfgA.fasta HHPred is a popular server for finding template models and for creating homology models. Go to https://toolkit.tuebingen.mpg.de/#/tools/hhpred. In the "Input" section, press "Upload file" and select your FASTA file (or paste the sequence to the box above). Press "Submit job". While the job is running, answer the following question: Question 2: What features (e.g., biological, structural), other than sequence alignments, do you think could be used for choosing a template structure for homology modeling? Try to think of at least two features. The results list the best template candidates from the Protein Data Bank. Press the "Help" button (top right corner) for detailed info. We'll use as template the 4th structure in the list, 2Q3G:A. You can see the alignment between 2Q3G and 1MFG, including some useful indicators, such as secondary structure predictions for 1MFG (E = sheet, H = helix, C = coiled loop), amino-acids sequence and the consensus sequence for multiple-sequence alignment to 1MFG. Question 3: What is the sequence identity between 2Q3G and 1MFG? Is 2Q3G expected to be a good template for this modeling task? Once we decide on the template we'll be using, we're left with some technical tasks of acquiring and cleaning the template structure, and generating a simple alignment file of both sequences, to be used in the next phases. Rosetta uses this alignment to determine which parts of the target can be modeled based on the template (i.e. copied from it), and which require explicit modeling. The alignment file can be created by simply copying the sequence lines from HHPred. We provide you with this file, along with others, for the purpose of this exercise. The last two steps of homology modeling (loop modeling, refinement) are also useful in other contexts, e.g. ab-initio modeling or interface design. In the ab-initio exercise, we've practiced refinement (via the fast relax protocol), and in this exercise we will practice loop modeling as an independent module. Loop Modeling Loop modeling is a common sub-problem of homology modeling and is often one of the toughest parts of the modeling, because loops have conformational diversity and inherent flexibility. They play an important role in how proteins interact with one another and with other molecules. Question 4: i. Loop regions are sometimes missing from crystal structures - why do you think this might happen? ii. Explain the relation between gaps in the sequence alignment of homologue proteins and loop regions. 2 Workshop in Computational Structural Biology (81813), Spring 2019 iii. Loops are known to be important for protein function. Explain how they contribute to function, and provide an example. Examine our Erbin2 structure in PyMOL. Pay special attention to the loop found at residues 10-17 - this is the loop we are about to model (note that in real life scenarios, we would model the loop regions after copying template segments from aligned regions, making this loop modeling task slightly easier than in full homology modeling, having exact "take-off" sites). Question 5: How flexible do you think this region is? Do you think it would be easy to model this loop? Do you think it's too long/short? Create a directory for this loop modeling session, placing the structure there: $> mkdir loopmodel $> cp 1mfgA.pdb loopmodel $> cd loopmodel We will now define which loops we want to model, using a loop file with a very simple format. Make a file (called loopfile, for example) that contains the following line: LOOP 10 17 This instructs Rosetta to model the loop between residue 10 to 17 using random positions within the loop as cutpoints (note that if a specific position is desired as the cutting point of the loop, a fourth column stating that position can be added). There are two main algorithms used in Rosetta to model loops, namely Kinematic loop closure (KIC) and coordinate cyclic descent (CCD). Here we will use KIC 1 . Run the loop modeling protocol using the loop file and input structure. loopmodel.linuxgccrelease \ -database $ROSETTA_DB \ -nstruct 15 \ -s 1mfgA.pdb \ -native 1mfgA.pdb \ -out:prefix myloop \ -loops:loop_file loopfile \ # define the loops here -loops:extended false \ # use the input loop as a starting conformation -loops:remodel perturb_kic \ # choose modeling protocol -out:file:fullatom \ > log Note For kinematic loop closure you don’t need fragment files. If you choose to use quick_ccd, replace the modeling protocol option with -loops:remodel quick_ccd and add -loops:frag_files aa1mfgA09_05.200_v1_3.gz aa1mfgA03_05.200_v1_3.gz none (which you can obtain from the Robetta server ). Let's examine how well we modeled the loop: Inspect the models in Pymol. Question 6: Is there any loop conformation that looks similar to the native? How many? 3 Workshop in Computational Structural Biology (81813), Spring 2019 In this run, you have used the -loops:extended false flag, which instructs Rosetta to use the input loop conformation as the starting structure for the loop modeling protocol. We ran the protocol, once with this flag set to true and once to false, only with -nstruct 5000. A pse with the lowest energy decoys is provided under $CSBW_HOME/resources/ex6/lowestEnergyLoops.pse. The score files for the two runs are also provided (their names are extendedstart_score.sc and startconf_score.sc). Question 7: a. How do results differ between the two runs, both with respect to similarity to native and with respect to scores? b. Here we use the information from the crystal structure when trying to model its loop. When starting from a homolog template, when would it be recommended to use its starting conformation and when would it not be advisable? If you recall, the 1MFG pdb is a complex of a PDZ domain binding a peptide. Let's add the peptide to the mix. Create a new folder - loopWithPep. Copy loopfile to the new directory and create there a clean copy of the complex: $> prody fetch 1mfg $> prody select "protein" 1mfg -o 1mfg Go to the new directory and run the same command as before to model the loop (except use 1mfg.pdb instead of 1mfgA.pdb both as the starting structure -s and as the native). In the meanwhile, inspect the complex (1mfg.pdb) in PyMOL. Question 8: (a) Why should the peptide influence the results of the loop modeling? (b) How many loops are close to native for this run? Why is that? We provide loopWithPep.pse and loopWithPep_score.sc (found in $CSBW_HOME/resources/ex6) with additional 10 conformations, the lowest-energy structures created with the same protocol only with -nstruct 5000, just in case you didn’t sample enough.
Recommended publications
  • NIH Public Access Author Manuscript Proteins
    NIH Public Access Author Manuscript Proteins. Author manuscript; available in PMC 2015 February 01. NIH-PA Author ManuscriptPublished NIH-PA Author Manuscript in final edited NIH-PA Author Manuscript form as: Proteins. 2014 February ; 82(0 2): 208–218. doi:10.1002/prot.24374. One contact for every twelve residues allows robust and accurate topology-level protein structure modeling David E. Kim, Frank DiMaio, Ray Yu-Ruei Wang, Yifan Song, and David Baker* Department of Biochemistry, University of Washington, Seattle 98195, Washington Abstract A number of methods have been described for identifying pairs of contacting residues in protein three-dimensional structures, but it is unclear how many contacts are required for accurate structure modeling. The CASP10 assisted contact experiment provided a blind test of contact guided protein structure modeling. We describe the models generated for these contact guided prediction challenges using the Rosetta structure modeling methodology. For nearly all cases, the submitted models had the correct overall topology, and in some cases, they had near atomic-level accuracy; for example the model of the 384 residue homo-oligomeric tetramer (Tc680o) had only 2.9 Å root-mean-square deviation (RMSD) from the crystal structure. Our results suggest that experimental and bioinformatic methods for obtaining contact information may need to generate only one correct contact for every 12 residues in the protein to allow accurate topology level modeling. Keywords protein structure prediction; rosetta; comparative modeling; homology modeling; ab initio prediction; contact prediction INTRODUCTION Predicting the three-dimensional structure of a protein given just the amino acid sequence with atomic-level accuracy has been limited to small (<100 residues), single domain proteins.
    [Show full text]
  • Homology Modeling and Analysis of Structure Predictions of the Bovine Rhinitis B Virus RNA Dependent RNA Polymerase (Rdrp)
    Int. J. Mol. Sci. 2012, 13, 8998-9013; doi:10.3390/ijms13078998 OPEN ACCESS International Journal of Molecular Sciences ISSN 1422-0067 www.mdpi.com/journal/ijms Article Homology Modeling and Analysis of Structure Predictions of the Bovine Rhinitis B Virus RNA Dependent RNA Polymerase (RdRp) Devendra K. Rai and Elizabeth Rieder * Foreign Animal Disease Research Unit, United States Department of Agriculture, Agricultural Research Service, Plum Island Animal Disease Center, Greenport, NY 11944, USA; E-Mail: [email protected] * Author to whom correspondence should be addressed; E-Mail: [email protected]; Tel.: +1-631-323-3177; Fax: +1-631-323-3006. Received: 3 May 2012; in revised form: 3 July 2012 / Accepted: 11 July 2012 / Published: 19 July 2012 Abstract: Bovine Rhinitis B Virus (BRBV) is a picornavirus responsible for mild respiratory infection of cattle. It is probably the least characterized among the aphthoviruses. BRBV is the closest relative known to Foot and Mouth Disease virus (FMDV) with a ~43% identical polyprotein sequence and as much as 67% identical sequence for the RNA dependent RNA polymerase (RdRp), which is also known as 3D polymerase (3Dpol). In the present study we carried out phylogenetic analysis, structure based sequence alignment and prediction of three-dimensional structure of BRBV 3Dpol using a combination of different computational tools. Model structures of BRBV 3Dpol were verified for their stereochemical quality and accuracy. The BRBV 3Dpol structure predicted by SWISS-MODEL exhibited highest scores in terms of stereochemical quality and accuracy, which were in the range of 2Å resolution crystal structures. The active site, nucleic acid binding site and overall structure were observed to be in agreement with the crystal structure of unliganded as well as template/primer (T/P), nucleotide tri-phosphate (NTP) and pyrophosphate (PPi) bound FMDV 3Dpol (PDB, 1U09 and 2E9Z).
    [Show full text]
  • Comparative Protein Structure Modeling of Genes and Genomes
    P1: FPX/FOZ/fop/fok P2: FHN/FDR/fgi QC: FhN/fgm T1: FhN January 12, 2001 16:34 Annual Reviews AR098-11 Annu. Rev. Biophys. Biomol. Struct. 2000. 29:291–325 Copyright c 2000 by Annual Reviews. All rights reserved COMPARATIVE PROTEIN STRUCTURE MODELING OF GENES AND GENOMES Marc A. Mart´ı-Renom, Ashley C. Stuart, Andras´ Fiser, Roberto Sanchez,´ Francisco Melo, and Andrej Sˇali Laboratories of Molecular Biophysics, Pels Family Center for Biochemistry and Structural Biology, Rockefeller University, 1230 York Ave, New York, NY 10021; e-mail: [email protected] Key Words protein structure prediction, fold assignment, alignment, homology modeling, model evaluation, fully automated modeling, structural genomics ■ Abstract Comparative modeling predicts the three-dimensional structure of a given protein sequence (target) based primarily on its alignment to one or more pro- teins of known structure (templates). The prediction process consists of fold assign- ment, target–template alignment, model building, and model evaluation. The number of protein sequences that can be modeled and the accuracy of the predictions are in- creasing steadily because of the growth in the number of known protein structures and because of the improvements in the modeling software. Further advances are nec- essary in recognizing weak sequence–structure similarities, aligning sequences with structures, modeling of rigid body shifts, distortions, loops and side chains, as well as detecting errors in a model. Despite these problems, it is currently possible to model with useful accuracy significant parts of approximately one third of all known protein sequences. The use of individual comparative models in biology is already rewarding and increasingly widespread.
    [Show full text]
  • FORCE FIELDS for PROTEIN SIMULATIONS by JAY W. PONDER
    FORCE FIELDS FOR PROTEIN SIMULATIONS By JAY W. PONDER* AND DAVIDA. CASEt *Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, 51. Louis, Missouri 63110, and tDepartment of Molecular Biology, The Scripps Research Institute, La Jolla, California 92037 I. Introduction. ...... .... ... .. ... .... .. .. ........ .. .... .... ........ ........ ..... .... 27 II. Protein Force Fields, 1980 to the Present.............................................. 30 A. The Am.ber Force Fields.............................................................. 30 B. The CHARMM Force Fields ..., ......... 35 C. The OPLS Force Fields............................................................... 38 D. Other Protein Force Fields ....... 39 E. Comparisons Am.ong Protein Force Fields ,... 41 III. Beyond Fixed Atomic Point-Charge Electrostatics.................................... 45 A. Limitations of Fixed Atomic Point-Charges ........ 46 B. Flexible Models for Static Charge Distributions.................................. 48 C. Including Environmental Effects via Polarization................................ 50 D. Consistent Treatment of Electrostatics............................................. 52 E. Current Status of Polarizable Force Fields........................................ 57 IV. Modeling the Solvent Environment .... 62 A. Explicit Water Models ....... 62 B. Continuum Solvent Models.......................................................... 64 C. Molecular Dynamics Simulations with the Generalized Born Model........
    [Show full text]
  • Structural Bioinformatics
    Structural Bioinformatics Sanne Abeln K. Anton Feenstra Centre for Integrative Bioinformatics (IBIVU), and Department of Computer Science, Vrije Universiteit, De Boelelaan 1081A, 1081 HV Amsterdam, Netherlands December 4, 2017 arXiv:1712.00425v1 [q-bio.BM] 1 Dec 2017 2 Abstract This chapter deals with approaches for protein three-dimensional structure prediction, starting out from a single input sequence with unknown struc- ture, the `query' or `target' sequence. Both template based and template free modelling techniques are treated, and how resulting structural models may be selected and refined. We give a concrete flowchart for how to de- cide which modelling strategy is best suited in particular circumstances, and which steps need to be taken in each strategy. Notably, the ability to locate a suitable structural template by homology or fold recognition is crucial; without this models will be of low quality at best. With a template avail- able, the quality of the query-template alignment crucially determines the model quality. We also discuss how other, courser, experimental data may be incorporated in the modelling process to alleviate the problem of missing template structures. Finally, we discuss measures to predict the quality of models generated. Structural Bioinformatics c Abeln & Feenstra, 2014-2017 Contents 7 Strategies for protein structure model generation5 7.1 Template based protein structure modelling . .7 7.1.1 Homology based Template Finding . .7 7.1.2 Fold recognition . .8 7.1.3 Generating the target-template alignment . 10 7.1.4 Generating a model . 11 7.1.5 Loop or missing substructure modelling . 12 7.2 Template-free protein structure modelling .
    [Show full text]
  • Ten Quick Tips for Homology Modeling of High-Resolution Protein 3D Structures
    PLOS COMPUTATIONAL BIOLOGY EDUCATION Ten quick tips for homology modeling of high- resolution protein 3D structures 1,2 1,2 1,2 Yazan HaddadID , Vojtech AdamID , Zbynek HegerID * 1 Department of Chemistry and Biochemistry, Mendel University in Brno, Brno, Czech Republic, 2 Central European Institute of Technology, Brno University of Technology, Brno, Czech Republic * [email protected] Author summary The purpose of this quick guide is to help new modelers who have little or no background in comparative modeling yet are keen to produce high-resolution protein 3D structures a1111111111 for their study by following systematic good modeling practices, using affordable personal a1111111111 computers or online computational resources. Through the available experimental 3D- a1111111111 structure repositories, the modeler should be able to access and use the atomic coordinates a1111111111 for building homology models. We also aim to provide the modeler with a rationale a1111111111 behind making a simple list of atomic coordinates suitable for computational analysis abiding to principles of physics (e.g., molecular mechanics). Keeping that objective in mind, these quick tips cover the process of homology modeling and some postmodeling computations such as molecular docking and molecular dynamics (MD). A brief section OPEN ACCESS was left for modeling nonprotein molecules, and a short case study of homology modeling Citation: Haddad Y, Adam V, Heger Z (2020) Ten is discussed. quick tips for homology modeling of high- resolution protein 3D structures. PLoS Comput Biol 16(4): e1007449. https://doi.org/10.1371/ journal.pcbi.1007449 Introduction Editor: Francis Ouellette, University of Toronto, Protein 3D-structure folding from a simple sequence of amino acids was seen as a very difficult CANADA problem in the past.
    [Show full text]
  • Methods for the Refinement of Protein Structure 3D Models
    International Journal of Molecular Sciences Review Methods for the Refinement of Protein Structure 3D Models Recep Adiyaman and Liam James McGuffin * School of Biological Sciences, University of Reading, Reading RG6 6AS, UK; [email protected] * Correspondence: l.j.mcguffi[email protected]; Tel.: +44-0-118-378-6332 Received: 2 April 2019; Accepted: 7 May 2019; Published: 1 May 2019 Abstract: The refinement of predicted 3D protein models is crucial in bringing them closer towards experimental accuracy for further computational studies. Refinement approaches can be divided into two main stages: The sampling and scoring stages. Sampling strategies, such as the popular Molecular Dynamics (MD)-based protocols, aim to generate improved 3D models. However, generating 3D models that are closer to the native structure than the initial model remains challenging, as structural deviations from the native basin can be encountered due to force-field inaccuracies. Therefore, different restraint strategies have been applied in order to avoid deviations away from the native structure. For example, the accurate prediction of local errors and/or contacts in the initial models can be used to guide restraints. MD-based protocols, using physics-based force fields and smart restraints, have made significant progress towards a more consistent refinement of 3D models. The scoring stage, including energy functions and Model Quality Assessment Programs (MQAPs) are also used to discriminate near-native conformations from non-native conformations. Nevertheless, there are often very small differences among generated 3D models in refinement pipelines, which makes model discrimination and selection problematic. For this reason, the identification of the most native-like conformations remains a major challenge.
    [Show full text]
  • Advances in Rosetta Protein Structure Prediction on Massively Parallel Systems
    UC San Diego UC San Diego Previously Published Works Title Advances in Rosetta protein structure prediction on massively parallel systems Permalink https://escholarship.org/uc/item/87g6q6bw Journal IBM Journal of Research and Development, 52(1) ISSN 0018-8646 Authors Raman, S. Baker, D. Qian, B. et al. Publication Date 2008 Peer reviewed eScholarship.org Powered by the California Digital Library University of California Advances in Rosetta protein S. Raman B. Qian structure prediction on D. Baker massively parallel systems R. C. Walker One of the key challenges in computational biology is prediction of three-dimensional protein structures from amino-acid sequences. For most proteins, the ‘‘native state’’ lies at the bottom of a free- energy landscape. Protein structure prediction involves varying the degrees of freedom of the protein in a constrained manner until it approaches its native state. In the Rosetta protein structure prediction protocols, a large number of independent folding trajectories are simulated, and several lowest-energy results are likely to be close to the native state. The availability of hundred-teraflop, and shortly, petaflop, computing resources is revolutionizing the approaches available for protein structure prediction. Here, we discuss issues involved in utilizing such machines efficiently with the Rosetta code, including an overview of recent results of the Critical Assessment of Techniques for Protein Structure Prediction 7 (CASP7) in which the computationally demanding structure-refinement process was run on 16 racks of the IBM Blue Gene/Le system at the IBM T. J. Watson Research Center. We highlight recent advances in high-performance computing and discuss future development paths that make use of the next-generation petascale (.1012 floating-point operations per second) machines.
    [Show full text]
  • Homology Modeling
    25 HOMOLOGY MODELING Elmar Krieger, Sander B. Nabuurs, and Gert Vriend The ultimate goal of protein modeling is to predict a structure from its sequence with an accuracy that is comparable to the best results achieved experimentally. This would allow users to safely use rapidly generated in silico protein models in all the contexts where today only experimental structures provide a solid basis: structure-based drug design, analysis of protein function, interactions, antigenic behavior, and rational design of proteins with increased stability or novel functions. In addition, protein modeling is the only way to obtain structural information if experimental techniques fail. Many proteins are simply too large for NMR analysis and cannot be crystallized for X-ray diffraction. Among the three major approaches to three-dimensional (3D) structure prediction described in this and the following two chapters, homology modeling is the easiest one. It is based on two major observations: 1. The structure of a protein is uniquely determined by its amino acid sequence (Epstain, Goldberger, and Anfinsen, 1963). Knowing the sequence should, at least in theory, suffice to obtain the structure. 2. During evolution, the structure is more stable and changes much slower than the associated sequence, so that similar sequences adopt practically identical structures, and distantly related sequences still fold into similar structures. This relationship was first identified by Chothia and Lesk (1986) and later quantified by Sander and Schneider (1991). Thanks to the exponential growth of the Protein Data Bank (PDB), Rost (1999) could recently derive a precise limit for this rule, shown in Figure 25.1.
    [Show full text]
  • Homology Modeling and Optimized Expression of Truncated IK Protein, Tik, As an Anti-Inflammatory Peptide
    molecules Article Homology Modeling and Optimized Expression of Truncated IK Protein, tIK, as an Anti-Inflammatory Peptide Yuyoung Song, Minseon Kim and Yongae Kim * Department of Chemistry, Hankuk University of Foreign Studies, Yongin 17035, Korea; [email protected] (Y.S.); [email protected] (M.K.) * Correspondence: [email protected]; Tel.: +82-2-2173-8705; Fax: +82-31-330-4566 Academic Editors: Manuela Pintado, Ezequiel Coscueta and María Emilia Brassesco Received: 11 August 2020; Accepted: 21 September 2020; Published: 23 September 2020 Abstract: Rheumatoid arthritis, caused by abnormalities in the autoimmune system, affects about 1% of the population. Rheumatoid arthritis does not yet have a proper treatment, and current treatment has various side effects. Therefore, there is a need for a therapeutic agent that can effectively treat rheumatoid arthritis without side effects. Recently, research on pharmaceutical drugs based on peptides has been actively conducted to reduce negative effects. Because peptide drugs are bio-friendly and bio-specific, they are characterized by no side effects. Truncated-IK (tIK) protein, a fragment of IK protein, has anti-inflammatory effects, including anti-rheumatoid arthritis activity. This study focused on the fact that tIK protein phosphorylates the interleukin 10 receptor. Through homology modeling with interleukin 10, short tIK epitopes were proposed to find the essential region of the sequence for anti-inflammatory activity. TH17 differentiation experiments were also performed with the proposed epitope. A peptide composed of 18 amino acids with an anti-inflammatory effect was named tIK-18mer. Additionally, a tIK 9-mer and a 14-mer were also found. The procedure for the experimental expression of the proposed tIK series (9-mer, 14-mer, and 18-mer) using bacterial strain is discussed.
    [Show full text]
  • Homology Modeling and Docking Studies of a 9-Fatty Acid Desaturase
    Homology modeling and docking studies of a 19-fatty acid desaturase from a Cold-tolerant Pseudomonas sp. AMS8 Lawal Garba1,2,3, Mohamad Ariff Mohamad Yussoff4, Khairul Bariyyah Abd Halim4, Siti Nor Hasmah Ishak1, Mohd Shukuri Mohamad Ali1,5, Siti Nurbaya Oslan1,5 and Raja Noor Zaliha Raja Abd. Rahman1,2 1 Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia 2 Department of Microbiology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia 3 Department of Microbiology, Faculty of Science, Gombe State University, Gombe, Gombe State, Nigeria 4 Department of Biotechnology, Kulliyyah of Science, International Islamic University Malaysia, Kuantan, Pahang Darul Makmur, Malaysia 5 Department of Biochemistry, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang, Selangor, Malaysia ABSTRACT Membrane-bound fatty acid desaturases perform oxygenated desaturation reactions to insert double bonds within fatty acyl chains in regioselective and stereoselective manners. The 19-fatty acid desaturase strictly creates the first double bond between C9 and 10 positions of most saturated substrates. As the three-dimensional structures of the bacterial membrane fatty acid desaturases are not available, relevant information about the enzymes are derived from their amino acid sequences, site-directed mutagenesis and domain swapping in similar membrane-bound desaturases. The cold-tolerant Pseudomonas sp. AMS8 was found to produce high amount of monounsaturated fatty acids at low temperature. Subsequently, an active 19-fatty acid desaturase was isolated and functionally expressed in Escherichia coli. In this paper we report homology modeling and docking studies of a 19-fatty acid desaturase from a Cold- Submitted 8 May 2017 tolerant Pseudomonas sp.
    [Show full text]
  • A New Computer Model for Evaluating the Selective Binding Affinity Of
    pharmaceuticals Article A New Computer Model for Evaluating the Selective Binding Affinity of Phenylalkylamines to T-Type Ca2+ Channels You Lu 1 and Ming Li 2,* 1 Center for Aging, School of Medicine, Tulane University, New Orleans, LA 70112, USA; [email protected] 2 Department of Physiology, School of Medicine, Tulane University, New Orleans, LA 70112, USA * Correspondence: [email protected]; Tel.: +1-504-988-8207 Abstract: To establish a computer model for evaluating the binding affinity of phenylalkylamines (PAAs) to T-type Ca2+ channels (TCCs), we created new homology models for both TCCs and a L-type calcium channel (LCC). We found that PAAs have a high affinity for domains I and IV of TCCs and a low affinity for domains III and IV of the LCC. Therefore, they should be considered as favorable candidates for TCC blockers. The new homology models were validated with some commonly recognized TCC blockers that are well characterized. Additionally, examples of the TCC blockers created were also evaluated using these models. Keywords: T-type calcium channel blocker; homology modeling; computer-aid drug design; virtual drug screening; L-type calcium channel 1. Introduction Citation: Lu, Y.; Li, M. A New As the only type of voltage-gated Ca2+ channels that are activated at or near resting Computer Model for Evaluating the membrane potentials, T-type Ca2+ channels (TCCs) play an important role in regulat- Selective Binding Affinity of 2+ ing [Ca ]i homeostasis in a variety of tissues, including pancreatic β-cells and tumor Phenylalkylamines to T-Type Ca2+ cells [1–4]. Therefore, TCC antagonists could be potentially useful for the treatment of Channels.
    [Show full text]