Topology Files Research Articles

Structure-based Mechanism of CMP-2-keto-3-deoxymanno-octulonic Acid Synthetase

The enzyme CMP-Kdo synthetase (KdsB) catalyzes the addition of 2-keto-3-deoxymanno-octulonic acid (Kdo) to CTP to form CMP-Kdo, a key reaction in the biosynthesis of lipopolysaccharide. The reaction catalyzed by KdsB and the related CMP-acylneuraminate synthase is unique among the sugar-activating enzymes in that the respective sugars are directly coupled to a cytosine monophosphate. Using inhibition studies, in combination with isothermal calorimetry, we show the substrate analogue 2beta-deoxy-Kdo to be a potent competitive inhibitor. The ligand-free Escherichia coli KdsB and ternary complex KdsB-CTP-2beta-deoxy-Kdo crystal structures reveal that Kdo binding leads to active site closure and repositioning of the CTP phosphates and associated Mg(2+) ion (Mg-B). Both ligands occupy conformations compatible with an S(n)2-type attack on the alpha-phosphate by the Kdo 2-hydroxyl group. Based on strong similarity with DNA/RNA polymerases, both in terms of overall chemistry catalyzed as well as active site configuration, we postulate a second Mg(2+) ion (Mg-A) is bound by the catalytically competent KdsB-CTP-Kdo ternary complex. Modeling of this complex reveals the Mg-A coordinated to the conserved Asp(100) and Asp(235) in addition to the CTP alpha-phosphate and both the Kdo carboxylic and 2-hydroxyl groups. EPR measurements on the Mn(2+)-substituted ternary complex support this model. We propose the KdsB/CNS sugar-activating enzymes catalyze the formation of activated sugars, such as the abundant CMP-5-N-acetylneuraminic acid, by recruitment of two Mg(2+) to the active site. Although each metal ion assists in correct positioning of the substrates and activation of the alpha-phosphate, Mg-A is responsible for activation of the sugar-hydroxyl group.

CHAMBER: Comprehensive support for CHARMM force fields within the AMBER software

AbstractThe similarity of the AMBER force field's energy functional form with that of the CHARMM force field, gives the potential for direct translation of common bonding and nonbonding terms, along with their parameters, present in CHARMM topology and parameter files, with the intent of evaluation within the AMBER software; specifically the SANDER and PMEMD dynamics engines. To this extent, we have created a tool, CHAMBER, which can take a CHARMM protein structure file (PSF), coordinate file (COR) and associated forcefield files, and convert these to an AMBER topology file (prmtop) and associated coordinate file (inpcrd). CHAMBER opens a conversion route which enables the simulation of CHARMM parameterized models using AMBER's PMEMD engine; thus providing improved serial efficiency as well as parallel efficiency over large numbers of CPUs. Significant effort has been expended in ensuring a true representation of the CHARMM force field in AMBER providing energies and forces that are the same to the limits of machine precision. This software will be released in the upcoming version 1.3 of the free AMBERTools suite.1 © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009

Synthesis and Structural Characterization of Sialic Acid−Glutamic Acid Hybrid Foldamers as Conformational Surrogates of α-2,8-Linked Polysialic Acid

Surface expression of alpha-(2,8)-linked polymers of sialic acid in adult tissues has been correlated with metastasis of several human cancers. One approach to chemotherapeutic intervention against the spread of these cancers involves the development of immunogenic molecules that elicit an antibody response against alpha-(2,8)-linked polysialic acids. Naturally occurring polysialic acids are not viable candidates because they are present during embryonic development and are recognized as self by the immune system. These natural polymers also have poor pharmacokinetic properties because they are readily degraded by neuraminidase enzymes. We have been interested in developing structural surrogates of polysialic acids in an effort to overcome these limitations. Reported herein are microwave-assisted solid-phase peptide syntheses and structural characterization studies of a series of alpha/delta hybrid peptides derived from Fmoc-Neu2en and Fmoc-Glu(OtBu)-OH. Conformational experiments including circular dichroism, NH/ND exchange, and ROESY in aqueous solution were performed to study the secondary structures of these hybrid foldamers. ROESY data were analyzed with the assistance of XPLOR-NIH that was modified to include parameter and topology files to accommodate unnatural amino acids and the delta amide linkages. The results indicate that stable secondary structure is dependent upon both the amino acid sequence and the configuration of Glu. The most stable foldamer was composed of a total of 6 residues beginning with L-Glu at the carboxy terminus and alternating Neu2en and L-Glu residues. In water, this foldamer adopts a right-handed helical conformation with 3.7 residues per turn, 7.4 A pitch, 5.8 A diameter, and a length of 18.5 A, which is stabilized by both classical C=O...H-N backbone interactions and by pyranose ring O and L-Glu HN H-bonding. These structural features orient the L-Glu carboxylates along the helical backbone with a periodicity that matches the carboxylate positions along the reported G2(+) left-handed helix of alpha-(2,8)-polysialic acid. However, the charge density of the foldamer is one-half that of the natural polymer. These findings provide a fundamental understanding of the factors that influence stable secondary structure in hybrid Neu2en/Glu systems, and the tools we have developed establish a viable platform for the rational design of alpha-(2,8)-polysialic acid surrogates.

PURY: a database of geometric restraints of hetero compounds for refinement in complexes with macromolecular structures

The number and variety of macromolecular structures in complex with `hetero' ligands is growing. The need for rapid delivery of correct geometric parameters for their refinement, which is often crucial for understanding the biological relevance of the structure, is growing correspondingly. The current standard for describing protein structures is the Engh–Huber parameter set. It is an expert data set resulting from selection and analysis of the crystal structures gathered in the Cambridge Structural Database (CSD). Clearly, such a manual approach cannot be applied to the vast and ever‐growing number of chemical compounds. Therefore, a database, named PURY, of geometric parameters of chemical compounds has been developed, together with a server that accesses it. PURY is a compilation of the whole CSD. It contains lists of atom classes and bonds connecting them, as well as angle, chirality, planarity and conformation parameters. The current compilation is based on CSD 5.28 and contains 1978 atom classes and 32 702 bonding, 237 068 angle, 201 860 dihedral and 64 193 improper geometric restraints. Analysis has confirmed that the restraints from the PURY database are suitable for use in macromolecular crystal structure refinement and should be of value to the crystallographic community. The database can be accessed through the web server http://pury.ijs.si/, which creates topology and parameter files from deposited coordinates in suitable forms for the refinement programs MAIN, CNS and REFMAC. In the near future, the server will move to the CSD website http://pury.ccdc.cam.ac.uk/.

Molecular Basis of Prodrug Activation by Human Valacyclovirase, an α-Amino Acid Ester Hydrolase

Chemical modification to improve biopharmaceutical properties, especially oral absorption and bioavailability, is a common strategy employed by pharmaceutical chemists. The approach often employs a simple structural modification and utilizes ubiquitous endogenous esterases as activation enzymes, although such enzymes are often unidentified. This report describes the crystal structure and specificity of a novel activating enzyme for valacyclovir and valganciclovir. Our structural insights show that human valacyclovirase has a unique binding mode and specificity for amino acid esters. Biochemical data demonstrate that the enzyme hydrolyzes esters of alpha-amino acids exclusively and displays a broad specificity spectrum for the aminoacyl moiety similar to tricorn-interacting aminopeptidase F1. Crystal structures of the enzyme, two mechanistic mutants, and a complex with a product analogue, when combined with biochemical analysis, reveal the key determinants for substrate recognition; that is, a flexible and mostly hydrophobic acyl pocket, a localized negative electrostatic potential, a large open leaving group-accommodating groove, and a pivotal acidic residue, Asp-123, after the nucleophile Ser-122. This is the first time that a residue immediately after the nucleophile has been found to have its side chain directed into the substrate binding pocket and play an essential role in substrate discrimination in serine hydrolases. These results as well as a phylogenetic analysis establish that the enzyme functions as a specific alpha-amino acid ester hydrolase. Valacyclovirase is a valuable target for amino acid ester prodrug-based oral drug delivery enhancement strategies.

MKTOP: a program for automatic construction of molecular topologies

Molecular mechanics methods are widely-used for simulations of biomolecules. A molecular topology file containing all atom types, which depend on the chemical environment, must be constructed. MKTOP is a free-software capable of automatic atom type recognition and construction of molecular topologies for GROMACS, being a useful tool for its users.

Stereospecific Proton Transfer by a Mobile Catalyst in Mammalian Fructose-1,6-bisphosphate Aldolase

Class I fructose-1,6-bisphosphate aldolases catalyze the interconversion between the enamine and iminium covalent enzymatic intermediates by stereospecific exchange of the pro(S) proton of the dihydroxyacetone-phosphate C3 carbon, an obligatory reaction step during substrate cleavage. To investigate the mechanism of stereospecific proton exchange, high resolution crystal structures of native and a mutant Lys(146) --> Met aldolase were solved in complex with dihydroxyacetone phosphate. The structural analysis revealed trapping of the enamine intermediate at Lys(229) in native aldolase. Mutation of conserved active site residue Lys(146) to Met drastically decreased activity and enabled trapping of the putative iminium intermediate in the crystal structure showing active site attachment by C-terminal residues 360-363. Attachment positions the conserved C-terminal Tyr(363) hydroxyl within 2.9A of the C3 carbon in the iminium in an orientation consistent with incipient re face proton transfer. We propose a catalytic mechanism by which the mobile C-terminal Tyr(363) is activated by the iminium phosphate via a structurally conserved water molecule to yield a transient phenate, whose developing negative charge is stabilized by a Lys(146) positive charge, and which abstracts the C3 pro(S) proton forming the enamine. An identical C-terminal binding mode observed in the presence of phosphate in the native structure corroborates Tyr(363) interaction with Lys(146) and is consistent with transient C terminus binding in the enamine. The absence of charge stabilization and of a mobile C-terminal catalyst explains the extraordinary stability of enamine intermediates in transaldolases.

H++: a server for estimating pKas and adding missing hydrogens to macromolecules

The structure and function of macromolecules depend critically on the ionization (protonation) states of their acidic and basic groups. A number of existing practical methods predict protonation equilibrium pK constants of macromolecules based upon their atomic resolution Protein Data Bank (PDB) structures; the calculations are often performed within the framework of the continuum electrostatics model. Unfortunately, these methodologies are complex, involve multiple steps and require considerable investment of effort. Our web server provides access to a tool that automates this process, allowing both experts and novices to quickly obtain estimates of pKs as well as other related characteristics of biomolecules such as isoelectric points, titration curves and energies of protonation microstates. Protons are added to the input structure according to the calculated ionization states of its titratable groups at the user-specified pH; the output is in the PQR (PDB + charges + radii) format. In addition, corresponding coordinate and topology files are generated in the format supported by the molecular modeling package AMBER. The server is intended for a broad community of biochemists, molecular modelers, structural biologists and drug designers; it can also be used as an educational tool in biochemistry courses.

Probing the Role of Asp-120(81) of Metallo-β-lactamase (IMP-1) by Site-directed Mutagenesis, Kinetic Studies, and X-ray Crystallography

Metallo-beta-lactamase IMP-1 is a di-Zn(II) metalloenzyme that efficiently hydrolyzes beta-lactam antibiotics. Wild-type (WT) IMP-1 has a conserved Asp-120(81) in the active site, which plays an important role in catalysis. To probe the catalytic role of Asp-120(81) in IMP-1, the IMP-1 mutants, D120(81)A and D120(81)E, were prepared by site-directed mutagenesis, and various kinetics studies were conducted. The IMP-1 mutants exhibited 10(2)-10(4)-fold drops in k(cat) values compared with WT despite the fact that they contained two Zn(II) ions in the active site. To evaluate the acid-base characteristics of Asp-120(81), the pH dependence for hydrolysis was examined by stopped-flow studies. No observable pK(a) values between pH 5 and 9 were found for WT and D120(81)A. The rapid mixing of equimolar amounts of nitrocefin and all enzymes failed to result in the detection of an anion intermediate of nitrocefin at 650 nm. These results suggest that Asp-120(81) of IMP-1 is not a factor in decreasing the pK(a) for the water bridging two Zn(II) ions and is not a proton donor to the anionic intermediate. In the case of D120(81)E, the nitrocefin hydrolysis product, which shows a maximum absorption at 460 nm, was bound to D120(81)E in the protonated form. The three-dimensional structures of D120(81)A and D120(81)E were also determined at 2.0 and 3.0 A resolutions, respectively. In the case of D120(81)E, the Zn-Zn distance was increased by 0.3 A compared with WT, due to the change in the coordination mode of Glu-120(81)OE1 and the positional shift in the conserved His-263(197) at the active site.

Binding of Two Flaviolin Substrate Molecules, Oxidative Coupling, and Crystal Structure of Streptomyces coelicolor A3(2) Cytochrome P450 158A2

Cytochrome P450 158A2 (CYP158A2) is encoded within a three-gene operon (sco1206-sco1208) in the prototypic soil bacterium Streptomyces coelicolor A3(2). This operon is widely conserved among streptomycetes. CYP158A2 has been suggested to produce polymers of flaviolin, a pigment that may protect microbes from UV radiation, in combination with the adjacent rppA gene, which encodes the type III polyketide synthase, 1,3,6,8-tetrahydroxynaphthalene synthase. Following cloning, expression, and purification of this cytochrome P450, we have shown that it can produce dimer and trimer products from the substrate flaviolin and that the structures of two of the dimeric products were established using mass spectrometry and multiple NMR methods. A comparison of the x-ray structures of ligand-free (1.75 angstroms) and flaviolin-bound (1.62 angstroms) forms of CYP158A2 demonstrates a major conformational change upon ligand binding that closes the entry into the active site, partly due to repositioning of the F and G helices. Particularly interesting is the presence of two molecules of flaviolin in the closed active site. The flaviolin molecules form a quasi-planar three-molecule stack including the heme of CYP158A2, suggesting that oxidative C-C coupling of these phenolic molecules leads to the production of flaviolin dimers.

Structural Basis for the Insensitivity of a Serine Enzyme (Palmitoyl-Protein Thioesterase) to Phenylmethylsulfonyl Fluoride

Palmitoyl-protein thioesterase-1 (PPT1) is a newly described lysosomal enzyme that hydrolyzes long chain fatty acids from lipid-modified cysteine residues in proteins. Deficiency in this enzyme results in a severe neurodegenerative storage disorder, infantile neuronal ceroid lipofuscinosis. Although the primary structure of PPT1 contains a serine lipase consensus sequence, the enzyme is insensitive to commonly used serine-modifying reagents phenylmethylsulfonyl fluoride (PMSF) and diisopropylfluorophosphate. In the current paper, we show that the active site serine in PPT1 is modified by a substrate analog of PMSF, hexadecylsulfonylfluoride (HDSF) in a specific and site-directed manner. The apparent K(i) of the inhibition was 125 micrometer (in the presence of 1.5 mm Triton X-100), and the catalytic rate constant for sulfonylation (k(2)) was 3.3/min, a value similar to previously described sulfonylation reactions. PPT1 was crystallized after inactivation with HDSF, and the structure of the inactive form was determined to 2.4 A resolution. The hexadecylsulfonyl was found to modify serine 115 and to snake through a narrow hydrophobic channel that would not accommodate an aromatic sulfonyl fluoride. Therefore, the geometry of the active site accounts for the reactivity of PPT1 with HDSF but not PMSF. These observations suggest a structural explanation as to why certain serine lipases are resistant to modification by commonly used serine-modifying reagents.

29] Protein data bank archives of three-dimensional macromolecular structures

Publisher Summary This chapter discusses the Protein Data Bank (PDB) archives of three-dimensional (3D) macromolecular structures. Several pieces of information related to an entry are archived by the PDB. In addition to the coordinate entry file, the PDB stores files related to the experiment such as structure factors, nuclear Overhauser effect (NOE) restraints, and lists of chemical shifts. Also archived are auxiliary files used in structure analysis and refinement such as X-PLOR parameter and topology files. Currently, the archives are managed as a set of individual files, and each entry may have several associated files. The PDB is in the process of building a relational database, 3Dbase that will replace the current data management and access system. A description of 3DBase, including an outline of the way users can access its contents is provided in the chapter. Coordinate entries in the PDB are stored in separate files, each of which reports the results of an experiment or analysis that elucidates the structure of proteins, nucleic acids, polysaccharides, and other biological macromolecules. Although most of the data are generated from single crystal X-ray diffraction studies, a growing number of PDB entries are from nuclear magnetic resonance (NMR) studies.

PRODRG, a program for generating molecular topologies and unique molecular descriptors from coordinates of small molecules.

A software package is described that operates on small molecules observed in the PDB collection of protein structures. Molecular topology files for many molecular modeling programs can be generated automatically. The three-dimensional coordinates of small molecules can be converted to molecular descriptor strings that encode them uniquely in order to enable small-molecule recognition, despite high variability in atom and molecule nomenclature. From this descriptor a plausible 3D structure can be regenerated using energy minimisation. Alternatively, an ensemble of structures can be generated using a distance-geometry-based algorithm.