Raman Spectra of Amino Acids and Peptides from Machine Learning Polarizabilities.

Raman spectroscopy is an important tool in the study of vibrational properties and composition of molecules, peptides, and even proteins. Raman spectra can be simulated based on the change of the electronic polarizability with vibrations, which can nowadays be efficiently obtained via machine learning models trained on first-principles data. However, the transferability of the models trained on small molecules to larger structures is unclear, and direct training on large structures is prohibitively expensive. In this work, we first train two machine learning models to predict the polarizabilities of all 20 amino acids. Both models are carefully benchmarked and compared to density functional theory (DFT) calculations, with the neural network method being found to offer better transferability. By combination of machine learning models with classical force field molecular dynamics, Raman spectra of all amino acids are also obtained and investigated, showing good agreement with experiments. The models are further extended to small peptides. We find that adding structures containing peptide bonds to the training set greatly improves predictions, even for peptides not included in training sets.

Substrate specificity and conformational flexibility properties of the Mycobacterium tuberculosis ß-oxidation trifunctional enzyme.

The Mycobacterium tuberculosis trifunctional enzyme (MtTFE) is an α2ß2 tetrameric enzyme. The α-chain harbors the 2E-enoyl-CoA hydratase (ECH) and 3S-hydroxyacyl-CoA dehydrogenase (HAD) activities and the ß-chain provides the 3-ketoacyl-CoA thiolase (KAT) activity. Enzyme kinetic data reported here show that medium and long chain enoyl-CoA molecules are preferred substrates for MtTFE. Modelling studies indicate how the linear medium and long acyl chains of these substrates can bind to each of the active sites. In addition, crystallographic binding studies have identified three new CoA binding sites which are different from the previously known CoA binding sites of the three TFE active sites. Structure comparisons provide new insights into the properties of ECH, HAD and KAT active sites of MtTFE. The interactions of the adenine moiety of CoA with loop-2 of the ECH active site cause a conformational change of this loop by which a competent ECH active site is formed. The NAD+ binding domain (domain C) of the HAD part of MtTFE has only a few interactions with the rest of the complex and adopts a range of open conformations, whereas the A-domain of the ECH part is rigidly fixed with respect to the HAD part. Two loops, the CB1-CA1 region and the catalytic CB4-CB5 loop, near the thiolase active site and the thiolase dimer interface, have high B-factors. Structure comparisons suggest that a competent and stable thiolase dimer is formed only when complexed with the α-chains, highlighting the importance of the assembly for the proper functioning of the complex.

Integrin α11ß1 is a receptor for collagen XIII.

Collagen XIII is a conserved transmembrane collagen mainly expressed in mesenchymal tissues. Previously, we have shown that collagen XIII modulates tissue development and homeostasis. Integrins are a family of receptors that mediate signals from the environment into the cells and vice versa. Integrin α11ß1 is a collagen receptor known to recognize the GFOGER (O=hydroxyproline) sequence in collagens. Interestingly, collagen XIII and integrin α11ß1 both have a role in the regulation of bone homeostasis. To study whether α11ß1 is a receptor for collagen XIII, we utilized C2C12 cells transfected to express α11ß1 as their only collagen receptor. The interaction between collagen XIII and integrin α11ß1 was also confirmed by surface plasmon resonance and pull-down assays. We discovered that integrin α11ß1 mediates cell adhesion to two collagenous motifs, namely GPKGER and GF(S)QGEK, that were shown to act as the recognition sites for the integrin α11-I domain. Furthermore, we studied the in vivo significance of the α11ß1-collagen XIII interaction by crossbreeding α11 null mice (Itga11-/-) with mice overexpressing Col13a1 (Col13a1oe). When we evaluated the bone morphology by microcomputed tomography, Col13a1oe mice had a drastic bone overgrowth followed by severe osteoporosis, whereas the double mutant mouse line showed a much milder bone phenotype. To conclude, our data identifies integrin α11ß1 as a new collagen XIII receptor and demonstrates that this ligand-receptor pair has a role in the maintenance of bone homeostasis.

Structure-based classification of FAD binding sites: A comparative study of structural alignment tools.

A total of six different structural alignment tools (TM-Align, TriangleMatch, CLICK, ProBis, SiteEngine and GA-SI) were assessed for their ability to perform two particular tasks: (i) discriminating FAD (flavin adenine dinucleotide) from non-FAD binding sites, and (ii) performing an all-to-all comparison on a set of 883 FAD binding sites for the purpose of classifying them. For the first task, the consistency of each alignment method was evaluated, showing that every method is able to distinguish FAD and non-FAD binding sites with a high Matthews correlation coefficient. Additionally, GA-SI was found to provide alignments different from those of the other approaches. The results obtained for the second task revealed more significant differences among alignment methods, as reflected in the poor correlation of their results and highlighted clearly by the independent evaluation of the structural superimpositions generated by each method. The classification itself was performed using the combined results of all methods, using the best result found for each comparison of binding sites. A number of different clustering methods (Single-linkage, UPGMA, Complete-linkage, SPICKER and k-Means clustering) were also used. The groups of similar binding sites (proteins) or clusters generated by the best performing method were further analyzed in terms of local sequence identity, local structural similarity and conservation of analogous contacts with the FAD ligands. Each of the clusters was characterized by a unique set of structural features or patterns, demonstrating that the groups generated truly reflect the structural diversity of FAD binding sites. Proteins 2016; 84:1728-1747. © 2016 Wiley Periodicals, Inc.

The structure of the CD3ζζ transmembrane dimer in lipid bilayers.

Virtually every aspect of the human adaptive immune response is controlled by T cells. The T cell receptor (TCR) complex is responsible for the recognition of foreign peptide sequences, forming the initial step in the elimination of germ-infected cells. The recognition leads to an extracellular conformational change that is transmitted intracellularly through the Cluster of Differentiation 3 (CD3) subunits of the TCR-CD3 complex. Here we address the interplay between the disulfide-linked CD3ζζ dimer, an essential signaling component of the TCR-CD3 complex, and its lipidic environment. The disulfide bond formation requires the absolute presence of a nearby conserved aspartic acid, a fact that has mystified the scientific community. We use atomistic simulation methods to demonstrate that the conserved aspartic acid pair of the CD3ζζ dimer leads to a deformation of the membrane. This deformation changes the local environment of the cysteines and promotes disulfide bond formation. We also investigate the role of a conserved Tyr, highlighting its possible role in the interaction with other transmembrane components of the TCR-CD3 complex.

The first identification of carbohydrate binding modules specific to chitosan.

Two carbohydrate binding modules (DD1 and DD2) belonging to CBM32 are located at the C terminus of a chitosanase from Paenibacillus sp. IK-5. We produced three proteins, DD1, DD2, and tandem DD1/DD2 (DD1+DD2), and characterized their binding ability. Transition temperature of thermal unfolding (Tm) of each protein was elevated by the addition of cello-, laminari-, chitin-, or chitosan-hexamer (GlcN)6. The Tm elevation (ΔTm) in DD1 was the highest (10.3 °C) upon the addition of (GlcN)6 and was markedly higher than that in DD2 (1.0 °C). A synergistic effect was observed (ΔTm = 13.6 °C), when (GlcN)6 was added to DD1+DD2. From isothermal titration calorimetry experiments, affinities to DD1 were not clearly dependent upon chain length of (GlcN)n; ΔGr° values were -7.8 (n = 6), -7.6 (n = 5), -7.6 (n = 4), -7.6 (n = 3), and -7.1 (n = 2) kcal/mol, and the value was not obtained for GlcN due to the lowest affinity. DD2 bound (GlcN)n with the lower affinities (ΔGr° = -5.0 (n = 3) ~ -5.2 (n = 6) kcal/mol). Isothermal titration calorimetry profiles obtained for DD1+DD2 exhibited a better fit when the two-site model was used for analysis and provided greater affinities to (GlcN)6 for individual DD1 and DD2 sites (ΔGr° = -8.6 and -6.4 kcal/mol, respectively). From NMR titration experiments, (GlcN)n (n = 2~6) were found to bind to loops extruded from the core ß-sandwich of individual DD1 and DD2, and the interaction sites were similar to each other. Taken together, DD1+DD2 is specific to chitosan, and individual modules synergistically interact with at least two GlcN units, facilitating chitosan hydrolysis.

Interaction of di-N-acetylchitobiosyl moranoline with a family GH19 chitinase from moss, Bryum coronatum.

Tri-N-acetylchitotriosyl moranoline, (GlcNAc)3-M, was previously shown to strongly inhibit lysozyme (Ogata M, Umemoto N, Ohnuma T, Numata T, Suzuki A, Usui T, Fukamizo T. 2013. A novel transition-state analogue for lysozyme, 4-O-ß-tri-Nacetylchitotriosyl moranoline, provided evidence supporting the covalent glycosyl-enzyme intermediate. J Biol Chem. 288:6072-6082). The findings prompted us to examine the interaction of di-N-acetylchitobiosyl moranoline, (GlcNAc)2-M, with a family GH19 chitinase from moss, Bryum coronatum (BcChi19A). Thermal unfolding experiments using BcChi19A and the catalytic acid-deficient mutant (BcChi19A-E61A) revealed that the transition temperature (Tm) was elevated by 4.3 and 5.8°C, respectively, upon the addition of (GlcNAc)2-M, while the chitin dimer, (GlcNAc)2, elevated Tm only by 1.0 and 1.4°C, respectively. By means of isothermal titration calorimetry, binding free energy changes for the interactions of (GlcNAc)3 and (GlcNAc)2-M with BcChi19A-E61A were determined to be -5.2 and -6.6 kcal/mol, respectively, while (GlcNAc)2 was found to interact with BcChi19A-E61A with markedly lower affinity. nuclear magnetic resonance titration experiments using (15)N-labeled BcChi19A and BcChi19A-E61A revealed that both (GlcNAc)2 and (GlcNAc)2-M interact with the region surrounding the catalytic center of the enzyme and that the interaction of (GlcNAc)2-M is markedly stronger than that of (GlcNAc)2 for both enzymes. However, (GlcNAc)2-M was found to moderately inhibit the hydrolytic reaction of chitin oligosaccharides catalyzed by BcChi19A (IC50 = 130-620 µM). A molecular dynamics simulation of BcChi19A in complex with (GlcNAc)2-M revealed that the complex is quite stable and the binding mode does not significantly change during the simulation. The moranoline moiety of (GlcNAc)2-M did not fit into the catalytic cleft (subsite -1) but was rather in contact with subsite +1. This situation may result in the moderate inhibition toward the BcChi19A-catalyzed hydrolysis.

Atomic resolution view into the structure-function relationships of the human myelin peripheral membrane protein P2.

P2 is a fatty acid-binding protein expressed in vertebrate peripheral nerve myelin, where it may function in bilayer stacking and lipid transport. P2 binds to phospholipid membranes through its positively charged surface and a hydrophobic tip, and accommodates fatty acids inside its barrel structure. The structure of human P2 refined at the ultrahigh resolution of 0.93 Šallows detailed structural analyses, including the full organization of an internal hydrogen-bonding network. The orientation of the bound fatty-acid carboxyl group is linked to the protonation states of two coordinating arginine residues. An anion-binding site in the portal region is suggested to be relevant for membrane interactions and conformational changes. When bound to membrane multilayers, P2 has a preferred orientation and is stabilized, and the repeat distance indicates a single layer of P2 between membranes. Simulations show the formation of a double bilayer in the presence of P2, and in cultured cells wild-type P2 induces membrane-domain formation. Here, the most accurate structural and functional view to date on P2, a major component of peripheral nerve myelin, is presented, showing how it can interact with two membranes simultaneously while going through conformational changes at its portal region enabling ligand transfer.

Proteome-wide analysis of lysine acetylation suggests its broad regulatory scope in Saccharomyces cerevisiae.

Post-translational modification of proteins by lysine acetylation plays important regulatory roles in living cells. The budding yeast Saccharomyces cerevisiae is a widely used unicellular eukaryotic model organism in biomedical research. S. cerevisiae contains several evolutionary conserved lysine acetyltransferases and deacetylases. However, only a few dozen acetylation sites in S. cerevisiae are known, presenting a major obstacle for further understanding the regulatory roles of acetylation in this organism. Here we use high resolution mass spectrometry to identify about 4000 lysine acetylation sites in S. cerevisiae. Acetylated proteins are implicated in the regulation of diverse cytoplasmic and nuclear processes including chromatin organization, mitochondrial metabolism, and protein synthesis. Bioinformatic analysis of yeast acetylation sites shows that acetylated lysines are significantly more conserved compared with nonacetylated lysines. A large fraction of the conserved acetylation sites are present on proteins involved in cellular metabolism, protein synthesis, and protein folding. Furthermore, quantification of the Rpd3-regulated acetylation sites identified several previously known, as well as new putative substrates of this deacetylase. Rpd3 deficiency increased acetylation of the SAGA (Spt-Ada-Gcn5-Acetyltransferase) complex subunit Sgf73 on K33. This acetylation site is located within a critical regulatory domain in Sgf73 that interacts with Ubp8 and is involved in the activation of the Ubp8-containing histone H2B deubiquitylase complex. Our data provides the first global survey of acetylation in budding yeast, and suggests a wide-ranging regulatory scope of this modification. The provided dataset may serve as an important resource for the functional analysis of lysine acetylation in eukaryotes.

Theoretical Investigations of a point mutation affecting H5 Hemagglutinin's receptor binding preference.

The avian influenza A H5N1 virus is a subtype of influenza A virus (IAV) that causes a highly infectious and severe respiratory illness in birds and poses significant economic losses in poultry farming. To infect host cell, the virus uses its surface glycoprotein named Hemagglutinin (HA) to recognize and to interact with the host cell receptor containing either α2,6- (SAα2,6 Gal) or α2,3-linked Sialic Acid (SAα2,3 Gal). The H5N1 virus has not yet acquired the capability for efficient human-to-human transmission. However, studies have demonstrated that even a single amino acid substitution in the HA can switch its glycan receptor preference from the avian-type SAα2,3 Gal to the human-type SAα2,6 Gal. The present study aims to explain the underlying mechanism of a mutation (D94N) on the H5 HA that causes the protein to change its glycan receptor-binding preference using molecular dynamics (MD) simulations. Our results reveal that the mutation alters the electrostatic interactions pattern near the HA receptor binding pocket, leading to a reduced stability for the HA-avian-type SAα2,3 Gal complex. On the other hand, the detrimental effect of the mutation D94N is not observed in the HA-human-type SAα2,6 Gal complex due to the glycan's capability to switch its topology.

An In Vivo and In Silico Approach Reveals Possible Sodium Channel Nav1.2 Inhibitors from Ficus religiosa as a Novel Treatment for Epilepsy.

Epilepsy is a neurological disease that affects approximately 50 million people worldwide. Despite an existing abundance of antiepileptic drugs, lifelong disease treatment is often required but could be improved with alternative drugs that have fewer side effects. Given that epileptic seizures stem from abnormal neuronal discharges predominately modulated by the human sodium channel Nav1.2, the quest for novel and potent Nav1.2 blockers holds promise for epilepsy management. Herein, an in vivo approach was used to detect new antiepileptic compounds using the maximum electroshock test on mice. Pre-treatment of mice with extracts from the Ficus religiosa plant ameliorated the tonic hind limb extensor phase of induced convulsions. Subsequently, an in silico approach identified potential Nav1.2 blocking compounds from F. religiosa using a combination of computational techniques, including molecular docking, prime molecular mechanics/generalized Born surface area (MM/GBSA) analysis, and molecular dynamics (MD) simulation studies. The molecular docking and MM/GBSA analysis indicated that out of 82 compounds known to be present in F. religiosa, seven exhibited relatively strong binding affinities to Nav1.2 that ranged from -6.555 to -13.476 kcal/mol; similar or with higher affinity than phenytoin (-6.660 kcal/mol), a known Na+-channel blocking antiepileptic drug. Furthermore, MD simulations revealed that two compounds: 6-C-glucosyl-8-C-arabinosyl apigenin and pelargonidin-3-rhamnoside could form stable complexes with Nav1.2 at 300 K, indicating their potential as lead antiepileptic agents. In summary, the combination of in vivo and in silico approaches supports the potential of F. religiosa phytochemicals as natural antiepileptic therapeutic agents.

An atomistic model for assembly of transmembrane domain of T cell receptor complex.

The T cell receptor (TCR) together with accessory cluster of differentiation 3 (CD3) molecules (TCR-CD3 complex) is a key component in the primary function of T cells. The nature of association of the transmembrane domains is of central importance to the assembly of the complex and is largely unknown. Using multiscale molecular modeling and simulations, we have investigated the structure and assembly of the TCRα-CD3ε-CD3δ transmembrane domains both in membrane and in micelle environments. We demonstrate that in a membrane environment the transmembrane basic residue of the TCR closely interacts with both of the transmembrane acidic residues of the CD3 dimer. In contrast, in a micelle the basic residue interacts with only one of the acidic residues. Simulations of a recent micellar nuclear magnetic resonance structure of the natural killer (NK) cell-activating NKG2C-DAP12-DAP12 trimer in a membrane further indicate that the environment significantly affects the way these trimers associate. Since the currently accepted model for transmembrane association is entirely based on a micellar structure, we propose a revised model for the association of transmembrane domains of the activating immune receptors in a membrane environment.

Periplasmic chitooligosaccharide-binding protein requires a three-domain organization for substrate translocation.

Periplasmic solute-binding proteins (SBPs) specific for chitooligosaccharides, (GlcNAc)n (n = 2, 3, 4, 5 and 6), are involved in the uptake of chitinous nutrients and the negative control of chitin signal transduction in Vibrios. Most translocation processes by SBPs across the inner membrane have been explained thus far by two-domain open/closed mechanism. Here we propose three-domain mechanism of the (GlcNAc)n translocation based on experiments using a recombinant VcCBP, SBP specific for (GlcNAc)n from Vibrio cholerae. X-ray crystal structures of unliganded or (GlcNAc)3-liganded VcCBP solved at 1.2-1.6 Å revealed three distinct domains, the Upper1, Upper2 and Lower domains for this protein. Molecular dynamics simulation indicated that the motions of the three domains are independent and that in the (GlcNAc)3-liganded state the Upper2/Lower interface fluctuated more intensively, compared to the Upper1/Lower interface. The Upper1/Lower interface bound two GlcNAc residues tightly, while the Upper2/Lower interface appeared to loosen and release the bound sugar molecule. The three-domain mechanism proposed here was fully supported by binding data obtained by thermal unfolding experiments and ITC, and may be applicable to other translocation systems involving SBPs belonging to the same cluster.

A class V chitinase from Arabidopsis thaliana: gene responses, enzymatic properties, and crystallographic analysis.

Expression of a class V chitinase gene (At4g19810, AtChiC) in Arabidopsis thaliana was examined by quantitative real-time PCR and by analyzing microarray data available at Genevestigator. The gene expression was induced by the plant stress-related hormones abscisic acid (ABA) and jasmonic acid (JA) and by the stress resulting from the elicitor flagellin, NaCl, and osmosis. The recombinant AtChiC protein was produced in E. coli, purified, and characterized with respect to the structure and function. The recombinant AtChiC hydrolyzed N-acetylglucosamine oligomers producing dimers from the non-reducing end of the substrates. The crystal structure of AtChiC was determined by the molecular replacement method at 2.0 Å resolution. AtChiC was found to adopt an (ß/α)(8) fold with a small insertion domain composed of an α-helix and a five-stranded ß-sheet. From docking simulation of AtChiC with pentameric substrate, the amino acid residues responsible for substrate binding were found to be well conserved when compared with those of the class V chitinase from Nicotiana tabacum (NtChiV). All of the structural and functional properties of AtChiC are quite similar to those obtained for NtChiV, and seem to be common to class V chitinases from higher plants.

A flexible loop controlling the enzymatic activity and specificity in a glycosyl hydrolase family 19 endochitinase from barley seeds (Hordeum vulgare L.).

To examine the role of the loop structure consisting of residues 70-82 (70-82 loop) localized to +3/4 subsite of the substrate binding cleft of a family GH-19 endochitinase from barley seeds, Trp72 and Trp82 were mutated, and the mutated enzymes (W72A, W82A, and W72A/W82A) were characterized. Thermal stability and specific activities toward glycol chitin and chitin hexasaccharide were significantly affected by the individual mutations. When N-acetylglucosamine hexamer was hydrolyzed by the wild type, the beta-anomer of the substrate was preferentially hydrolyzed, producing the trimer predominantly and the dimer and tetramer in lesser amounts. When the mutated enzymes were used instead of the wild type, the enzyme cleavage sites in the hexamer substrate were clearly shifted, and the beta-anomer selectivity was eliminated. The mutation effects on the enzymatic activity and stability were much more substantial in W82A than in W72A, but surprisingly the effects of the W82A/W72A double mutation were intermediate between those of the two single mutations. A molecular dynamics simulation of the wild type and the Trp-mutated enzymes indicated that the 70-82 loop becomes more flexible upon mutation and the flexibility increases in the order of W72A, W72A/W82A and W82A. We conclude that Trp72 interacts with the sugar residue but Trp82 modulates the loop flexibility, which controls the protein stability and enzymatic properties. These tryptophan residues are likely to interact with each other, resulting in the non-additivity of mutational effects.

Inclusion of ionization states of ligands in affinity calculations.

When estimating binding affinities of a ligand, which can exists in multiple forms, for a target molecule, one must consider all possible competing equilibria. Here, a method is presented that estimates the contribution of the protonation equilibria of a ligand in solution to the measured or calculated binding affinity. The method yields a correction to binding constants that are based on the total concentration of inhibitor (the sum of all ionized forms of the inhibitor in solution) to account for the complexed form of the inhibitor only. The method is applied to the calculation of the difference in binding affinity of two inhibitors, 2-phosphoglycolate (PGA) and its phoshonate analog 3-phosphonopropionate (3PP), for the glycolytic enzyme triosephosphate isomerase. Both inhibitors have three titrating sites and exist in solution as a mixture of different forms. In this case the form that actually binds to the enzyme is present at relative low concentrations. The contributions of the alternative forms to the difference in binding energies is estimated by means of molecular dynamics simulations and corrections. The inhibitors undergo a pK(a) shift upon binding that is estimated by ab initio calculations. An interesting finding is that the affinity difference of the two inhibitors is not due to different interactions in the active site of the enzyme, but rather due to the difference in the solvation properties of the inhibitors.

A DFT study on the formation of a phosphohistidine intermediate in prostatic acid phosphatase.

Histidine phosphatases are a class of enzymes that are characterized by the presence of a conserved RHGXRXP motif. This motif contains a catalytic histidine that is being phosphorylated in the course of a dephosphorylation reaction catalyzed by these enzymes. Prostatic acid phosphatase (PAP) is one such enzyme. The dephosphorylation of phosphotyrosine by PAP is a two-step process. The first step involves the transfer of a phosphate group from the substrate to the histidine (His12). The present study reports on the details of the first step of this reaction, which was investigated using a series of quantum chemistry calculations. A number of quantum models were constructed containing various residues that were thought to play a role in the mechanism. In all these models, the transition state displayed an associative character. The transition state is stabilized by three active site arginines (Arg11, Arg15, and Arg79), two of which belong to the aforementioned conserved motif. The work also demonstrated that His12 could act as a nucleophile. The enzyme is further characterized by a His257-Asp258 motif. The role of Asp258 has been elusive. In this work, we propose that Asp258 acts as a proton donor which becomes protonated when the substrate enters the binding pocket. Evidence is also obtained that the transfer of a proton from Asp258 to the leaving group is possibly mediated by a water molecule in the active site. The work also underlines the importance of His257 in lowering the energy barrier for the nucleophilic attack.

Computer modeling of brain tumor growth.

An important objective of brain tumor modeling is to predict the progression of tumors so as to provide guidance about the best possible medical treatment to halt or slow the tumor's growth. Such computer models also provide a deeper insight into the physiology of tumors. In addition, one can study various what-if scenarios, for instance, investigating the response of tumors following the administration of a drug or a variety of drugs. Abrupt changes in growth rate can also be important for surgical decision-making. Despite increased interest in modeling techniques, relatively little progress has been made in improving such technologies. One problem is the limited data available from patients, typically 1 to 3 MRI (magnetic resonance imaging) sessions, from which one has to extrapolate the type of tumor so as to successfully predict its evolution over time. Here, the biological and clinical aspects of tumor growth and treatment with surgery, radiotherapy and drugs are discussed in the light of a patient with a brain tumor showing accelerated growth over time. Then, the contributions of mathematical modeling of tumor growth and effects of treatment are presented. Current tumor growth models can be roughly divided in three main categories, (i) cellular and microscopic models that emphasize isolated cell behavior, (ii) macroscopic models that concentrate on the development of cell density over time, and (iii) hybrid approaches that contain elements of both microscopic and macroscopic models. The mathematical theory that underlies these simulation methods is remarkably similar to the physical theory that forms the basis of protein modeling and molecular mechanics tools. A severe limitation of current models is that they are in fact not patient-specific at all.

Crystal structure of yeast peroxisomal multifunctional enzyme: structural basis for substrate specificity of (3R)-hydroxyacyl-CoA dehydrogenase units.

(3R)-hydroxyacyl-CoA dehydrogenase is part of multifunctional enzyme type 2 (MFE-2) of peroxisomal fatty acid beta-oxidation. The MFE-2 protein from yeasts contains in the same polypeptide chain two dehydrogenases (A and B), which possess difference in substrate specificity. The crystal structure of Candida tropicalis (3R)-hydroxyacyl-CoA dehydrogenase AB heterodimer, consisting of dehydrogenase A and B, determined at the resolution of 2.2A, shows overall similarity with the prototypic counterpart from rat, but also important differences that explain the substrate specificity differences observed. Docking studies suggest that dehydrogenase A binds the hydrophobic fatty acyl chain of a medium-chain-length ((3R)-OH-C10) substrate as bent into the binding pocket, whereas the short-chain substrates are dislocated by two mechanisms: (i) a short-chain-length 3-hydroxyacyl group ((3R)-OH-C4) does not reach the hydrophobic contacts needed for anchoring the substrate into the active site; and (ii) Leu44 in the loop above the NAD(+) cofactor attracts short-chain-length substrates away from the active site. Dehydrogenase B, which can use a (3R)-OH-C4 substrate, has a more shallow binding pocket and the substrate is correctly placed for catalysis. Based on the current structure, and together with the structure of the 2-enoyl-CoA hydratase 2 unit of yeast MFE-2 it becomes obvious that in yeast and mammalian MFE-2s, despite basically identical functional domains, the assembly of these domains into a mature, dimeric multifunctional enzyme is very different.

A previously unobserved conformation for the human Pex5p receptor suggests roles for intrinsic flexibility and rigid domain motions in ligand binding.

BACKGROUND: The C-terminal tetratricopeptide (TPR) repeat domain of Pex5p recognises proteins carrying a peroxisomal targeting signal type 1 (PTS1) tripeptide in their C-terminus. Previously, structural data have been obtained from the TPR domain of Pex5p in both the liganded and unliganded states, indicating a conformational change taking place upon cargo protein binding. Such a conformational change would be expected to play a major role both during PTS1 protein recognition as well as in cargo release into the peroxisomal lumen. However, little information is available on the factors that may regulate such structural changes. RESULTS: We have used a range of biophysical and computational methods to further analyse the conformational flexibility and ligand binding of Pex5p. A new crystal form for the human Pex5p C-terminal domain (Pex5p(C)) was obtained in the presence of Sr2+ ions, and the structure presents a novel conformation, distinct from all previous liganded and apo crystal structures for Pex5p(C). The difference relates to a near-rigid body movement of two halves of the molecule, and this movement is different from that required to reach a ring-like conformation upon PTS1 ligand binding. The bound Sr2+ ion changes the dynamic properties of Pex5p(C) affecting its conformation, possibly by making the Sr2+-binding loop - located near the hinge region for the observed domain motions - more rigid. CONCLUSION: The current data indicate that Pex5p(C) is able to sample a range of conformational states in the absence of bound PTS1 ligand. The domain movements between various apo conformations are distinct from those involved in ligand binding, although the differences between all observed conformations so far can be characterised by the movement of the two halves of Pex5p(C) as near-rigid bodies with respect to each other.