Omit Maps Research Articles

Direct Observation of the Protonation States in the Mutant Green Fluorescent Protein.

Neutron crystallography has been used to elucidate the protonation states for the enhanced green fluorescent protein, which has revolutionized imaging technologies. The structure has a deprotonated hydroxyl group in the fluorescent chromophore. Also, the protonation states of His148 and Thr203, as well as the orientation of a critical water molecule in direct contact with the chromophore, could be determined. The results demonstrate that the deprotonated hydroxyl group in the chromophore and the nitrogen atom ND1 in His148 are charged negatively and positively, respectively, forming an ion pair. The position of the two deuterium atoms in the critical water molecule appears to be displaced slightly toward the acceptor oxygen atoms according to their omit maps. This displacement implies the formation of an intriguing electrostatic potential realized inside of the protein. Our findings provide new insights into future protein design strategies along with developments in quantum chemical calculations.

X-ray crystallographic studies on the hydrogen isotope effects of green fluorescent protein at sub-ångström resolutions.

Hydrogen atoms are critical to the nature and properties of proteins, and thus deuteration has the potential to influence protein function. In fact, it has been reported that some deuterated proteins show different physical and chemical properties to their protiated counterparts. Consequently, it is important to investigate protonation states around the active site when using deuterated proteins. Here, hydrogen isotope effects on the S65T/F99S/M153T/V163A variant of green fluorescent protein (GFP), in which the deprotonated B form isdominant at pH 8.5, were investigated. The pH/pD dependence of the absorption and fluorescence spectra indicates that the protonation state of the chromophore is the same in protiated GFP in H2O and protiated GFP in D2O at pH/pD 8.5, while the pKa of the chromophore became higher in D2O. Indeed, X-ray crystallographic analyses at sub-ångström resolution revealed no apparent changes in the protonation state of the chromophore between the two samples. However, detailed comparisons of the hydrogen OMIT maps revealed that the protonation state of His148 in the vicinity of the chromophore differed between the two samples. This indicates that protonation states around the active site should be carefully adjusted to be the same as those of the protiated protein when neutron crystallographic analyses of proteins are performed.

Molstack: A platform for interactive presentations of electron density and cryo-EM maps and their interpretations.

In the Special Issue on Tools for Protein Science in 2018, we presented Molstack: a concept of a cloud-based platform for sharing electron density maps and their interpretations. Molstack is a web platform that allows the interactive visualization of density maps through the simultaneous presentation of multiple datasets and models in a way that allows for easy pairwise comparison. We anticipated that the users of this conceptually simple platform would find many different uses for their projects, and we were not mistaken. We have observed researchers use Molstack to present experimental evidence for their models in the form of electron density maps, omit maps, and anomalous difference density maps. Users also employed Molstack to present alternative interpretations of densities, including rerefinements and speculative interpretations. While we anticipated these types of projects to be the main use cases, we were pleased to see Molstack used to display superpositions of different models, as a tool for story-driven presentations, and for collaboration as well. Here, we present developments in the platform that were driven by user feedback, highlight several cases that used Molstack to enhance the publication, and discuss possible directions for the platform.

Polder maps: improving OMIT maps by excluding bulk solvent.

The crystallographic maps that are routinely used during the structure-solution workflow are almost always model-biased because model information is used for their calculation. As these maps are also used to validate the atomic models that result from model building and refinement, this constitutes an immediate problem: anything added to the model will manifest itself in the map and thus hinder the validation. OMIT maps are a common tool to verify the presence of atoms in the model. The simplest way to compute an OMIT map is to exclude the atoms in question from the structure, update the corresponding structure factors and compute a residual map. It is then expected that if these atoms are present in the crystal structure, the electron density for the omitted atoms will be seen as positive features in this map. This, however, is complicated by the flat bulk-solvent model which is almost universally used in modern crystallographic refinement programs. This model postulates constant electron density at any voxel of the unit-cell volume that is not occupied by the atomic model. Consequently, if the density arising from the omitted atoms is weak then the bulk-solvent model may obscure it further. A possible solution to this problem is to prevent bulk solvent from entering the selected OMIT regions, which may improve the interpretative power of residual maps. This approach is called a polder (OMIT) map. Polder OMIT maps can be particularly useful for displaying weak densities of ligands, solvent molecules, side chains, alternative conformations and residues both in terminal regions and in loops. The tools described in this manuscript have been implemented and are available in PHENIX.

Volumic omit maps in ab initio dual-space phasing.

Alternating-projection-type dual-space algorithms have a clear construction, but are susceptible to stagnation and, thus, inefficient for solving the phase problem ab initio. To improve this behaviour new omit maps are introduced, which are real-space perturbations applied periodically during the iteration process. The omit maps are called volumic, because they delete some predetermined subvolume of the unit cell without searching for atomic regions or analysing the electron density in any other way. The basic algorithms of positivity, histogram matching and low-density elimination are tested by their solution statistics. It is concluded that, while all these algorithms based on weak constraints are practically useless in their pure forms, appropriate volumic omit maps can transform them to practically useful methods. In addition, the efficiency of the already useful reflector-type charge-flipping algorithm can be further improved. Itis important that these results are obtained by using non-sharpened structure factors and without any weighting scheme or reciprocal-space perturbation. The mathematical background of volumic omit maps and their expected applications are also discussed.

A Threonine Stabilizes the NiC and NiR Catalytic Intermediates of [NiFe]-hydrogenase

The heterodimeric [NiFe] hydrogenase from Desulfovibrio fructosovorans catalyzes the reversible oxidation of H2 into protons and electrons. The catalytic intermediates have been attributed to forms of the active site (NiSI, NiR, and NiC) detected using spectroscopic methods under potentiometric but non-catalytic conditions. Here, we produced variants by replacing the conserved Thr-18 residue in the small subunit with Ser, Val, Gln, Gly, or Asp, and we analyzed the effects of these mutations on the kinetic (H2 oxidation, H2 production, and H/D exchange), spectroscopic (IR, EPR), and structural properties of the enzyme. The mutations disrupt the H-bond network in the crystals and have a strong effect on H2 oxidation and H2 production turnover rates. However, the absence of correlation between activity and rate of H/D exchange in the series of variants suggests that the alcoholic group of Thr-18 is not necessarily a proton relay. Instead, the correlation between H2 oxidation and production activity and the detection of the NiC species in reduced samples confirms that NiC is a catalytic intermediate and suggests that Thr-18 is important to stabilize the local protein structure of the active site ensuring fast NiSI-NiC-NiR interconversions during H2 oxidation/production.

Structure, Activity, and Substrate Selectivity of the Orf6 Thioesterase from Photobacterium profundum

Thioesterase activity is typically required for the release of products from polyketide synthase enzymes, but no such enzyme has been characterized in deep-sea bacteria associated with the production of polyunsaturated fatty acids. In this work, we have expressed and purified the Orf6 thioesterase from Photobacterium profundum. Enzyme assays revealed that Orf6 has a higher specific activity toward long-chain fatty acyl-CoA substrates (palmitoyl-CoA and eicosapentaenoyl-CoA) than toward short-chain or aromatic acyl-CoA substrates. We determined a high resolution (1.05 Å) structure of Orf6 that reveals a hotdog hydrolase fold arranged as a dimer of dimers. The putative active site of this structure is occupied by additional electron density not accounted for by the protein sequence, consistent with the presence of an elongated compound. A second crystal structure (1.40 Å) was obtained from a crystal that was grown in the presence of Mg(2+), which reveals the presence of a binding site for divalent cations at a crystal contact. The Mg(2+)-bound structure shows localized conformational changes (root mean square deviation of 1.63 Å), and its active site is unoccupied, suggesting a mechanism to open the active site for substrate entry or product release. These findings reveal a new thioesterase enzyme with a preference for long-chain CoA substrates in a deep-sea bacterium whose potential range of applications includes bioremediation and the production of biofuels.

Structural Basis of the Interaction of MbtH-like Proteins, Putative Regulators of Nonribosomal Peptide Biosynthesis, with Adenylating Enzymes

The biosynthesis of nonribosomally formed peptides (NRPs), which include important antibiotics such as vancomycin, requires the activation of amino acids through adenylate formation. The biosynthetic gene clusters of NRPs frequently contain genes for small, so-called MbtH-like proteins. Recently, it was discovered that these MbtH-like proteins are required for some of the adenylation reactions in NRP biosynthesis, but the mechanism of their interaction with the adenylating enzymes has remained unknown. In this study, we determined the structure of SlgN1, a 3-methylaspartate-adenylating enzyme involved in the biosynthesis of the hybrid polyketide/NRP antibiotic streptolydigin. SlgN1 contains an MbtH-like domain at its N terminus, and our analysis defines the parameters required for an interaction between MbtH-like domains and an adenylating enzyme. Highly conserved tryptophan residues of the MbtH-like domain critically contribute to this interaction. Trp-25 and Trp-35 form a cleft on the surface of the MbtH-like domain, which accommodates the alanine side chain of Ala-433 of the adenylating domain. Mutation of Ala-433 to glutamate abolished the activity of SlgN1. Mutation of Ser-23 of the MbtH-like domain to tyrosine resulted in strongly reduced activity. However, the activity of this S23Y mutant could be completely restored by addition of the intact MbtH-like protein CloY from another organism. This suggests that the interface found in the structure of SlgN1 is the genuine interface between MbtH-like proteins and adenylating enzymes.

Structural Basis for dsRNA Recognition, Filament Formation, and Antiviral Signal Activation by MDA5

MDA5, a viral double-stranded RNA (dsRNA) receptor, shares sequence similarity and signaling pathways with RIG-I yet plays essential functions in antiviral immunity through distinct specificity for viral RNA. Revealing the molecular basis for the functional divergence, we report here the crystal structure ofMDA5 bound to dsRNA, which shows how, using the same domain architecture, MDA5 recognizes the internal duplex structure, whereas RIG-I recognizes the terminus of dsRNA. We further show that MDA5 uses direct protein-protein contacts to stack along dsRNA in a head-to-tail arrangement, and that the signaling domain (tandem CARD), which decorates the outside of the core MDA5 filament, also has an intrinsic propensity to oligomerize into an elongated structure that activates the signaling adaptor, MAVS. These data support a model in which MDA5 uses long dsRNA as a signaling platform to cooperatively assemble the core filament, which in turn promotes stochastic assembly of the tandem CARD oligomers for signaling.

Crystal Structure of a Voltage-gated K+ Channel Pore Module in a Closed State in Lipid Membranes

Voltage-gated K(+) channels underlie the electrical excitability of cells. Each subunit of the functional tetramer consists of the tandem fusion of two modules, an N-terminal voltage-sensor and a C-terminal pore. To investigate how sensor coupling to the pore generates voltage-dependent channel opening, we solved the crystal structure and characterized the function of a voltage-gated K(+) channel pore in a lipid membrane. The structure of a functional channel in a membrane environment at 3.1 Å resolution establishes an unprecedented connection between channel structure and function. The structure is unique in delineating an ion-occupied ready to conduct selectivity filter, a confined aqueous cavity, and a closed activation gate, embodying a dynamic entity trapped in an unstable closed state.

Substrate binds in the S1 site of the F253A mutant of LeuT, a neurotransmitter sodium symporter homologue

LeuT serves as the model protein for understanding the relationships between structure, mechanism and pharmacology in neurotransmitter sodium symporters (NSSs). At the present time, however, there is a vigorous debate over whether there is a single high-affinity substrate site (S1) located at the original, crystallographically determined substrate site or whether there are two high-affinity substrates sites, one at the primary or S1 site and the other at a second site (S2) located at the base of the extracellular vestibule. In an effort to address the controversy over the number of high-affinity substrate sites in LeuT, one group studied the F253A mutant of LeuT and asserted that in this mutant substrate binds exclusively to the S2 site and that 1 mM clomipramine entirely ablates substrate binding to the S2 site. Here we study the binding of substrate to the F253A mutant of LeuT using ligand binding and X-ray crystallographic methods. Both experimental methods unambiguously show that substrate binds to the S1 site of the F253A mutant and that binding is retained in the presence of 1 mM clomipramine. These studies, in combination with previous work, are consistent with a mechanism for LeuT that involves a single high-affinity substrate binding site.

Structural studies of the effect that dimethyl sulfoxide (DMSO) has on cisplatin and carboplatin binding to histidine in a protein

The anticancer complexes cisplatin and carboplatin target the DNA major groove, forming intrastrand and interstrand cross-links between guanine bases through their N7 atoms, causing distortion of the DNA helix and apoptotic cell death. A major side effect of these drugs is toxicity, which is caused via binding to many proteins in the body. A range of crystallographic studies have been carried out involving the cocrystallization of hen egg-white lysozyme (HEWL) as a test protein with cisplatin and carboplatin in aqueous and dimethyl sulfoxide (DMSO) conditions. Different cryoprotectants, glycerol and Paratone, were used for each of the cisplatin and carboplatin cocrystallization cases, while silicone oil was used for studies involving N-acetylglucosamine (NAG). Both cisplatin and carboplatin do not bind to HEWL in aqueous media on the timescales of the conditions used here, but upon addition of DMSO two molecules of cisplatin or carboplatin bind either side of His15, which is the only His residue in lysozyme and is assumed to be an imidazolyl anion or a chemical resonance moiety, i.e. both imidazole N atoms are chemically reactive. To identify the platinum-peak positions in the 'with DMSO conditions', anomalous scattering maps were calculated as a cross-check with the F(o) - F(c) OMIT maps. Platinum-occupancy σ values were established using three different software programs in each case. The use of EVAL15 to process all of the diffraction data sets provided a consistent platform for a large ensemble of data sets for the various protein and platinum-compound model refinements with REFMAC5 and then SHELXTL. Overall, this extensive set of crystallization and cryoprotectant conditions allowed a systematic evaluation of cisplatin and carboplatin binding to lysozyme as a test protein via detailed X-ray crystal structure characterizations. DMSO is used as a super-solvent for drug delivery as it is deemed to cause no effect upon drug binding. However, these results show that addition of DMSO causes the platinum anticancer drugs to bind to HEWL. This effect should be considered in toxicity assessments of these drugs and perhaps more widely.

Crystal Structure of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-associated Csn2 Protein Revealed Ca2+-dependent Double-stranded DNA Binding Activity

Clustered regularly interspaced short palindromic repeats (CRISPR) and their associated protein genes (cas genes) are widespread in bacteria and archaea. They form a line of RNA-based immunity to eradicate invading bacteriophages and malicious plasmids. A key molecular event during this process is the acquisition of new spacers into the CRISPR loci to guide the selective degradation of the matching foreign genetic elements. Csn2 is a Nmeni subtype-specific cas gene required for new spacer acquisition. Here we characterize the Enterococcus faecalis Csn2 protein as a double-stranded (ds-) DNA-binding protein and report its 2.7 Å tetrameric ring structure. The inner circle of the Csn2 tetrameric ring is ∼26 Å wide and populated with conserved lysine residues poised for nonspecific interactions with ds-DNA. Each Csn2 protomer contains an α/β domain and an α-helical domain; significant hinge motion was observed between these two domains. Ca(2+) was located at strategic positions in the oligomerization interface. We further showed that removal of Ca(2+) ions altered the oligomerization state of Csn2, which in turn severely decreased its affinity for ds-DNA. In summary, our results provided the first insight into the function of the Csn2 protein in CRISPR adaptation by revealing that it is a ds-DNA-binding protein functioning at the quaternary structure level and regulated by Ca(2+) ions.

The Structural Basis of Endocannabinoid Oxygenation by Cyclooxygenase-2

The cyclooxygenases (COX-1 and COX-2) oxygenate arachidonic acid (AA) in the committed step of prostaglandin biogenesis. Substitutions of I434V, H513R, and I523V constitute the only differences in residues lining the cyclooxygenase channel between COX-1 and COX-2. These changes create a hydrophobic pocket in COX-2, with Arg-513 located at the base of the pocket, which has been exploited in the design of COX-2-selective inhibitors. Previous studies have shown that COX-2, but not COX-1, can oxygenate endocannabinoid substrates, including 2-arachidonoyl glycerol (2-AG). To investigate the isoform-specific structural basis of endocannabinoid binding to COX-2, we determined the crystal structure of the 2-AG isomer 1-arachidonoyl glycerol (1-AG) in complex with wild type and R513H murine (mu) COX-2 to 2.2 and 2.35 Å, respectively, and R513H muCOX-2 in complex with AA to 2.45 Å resolution. The 2,3-dihydroxypropyl moiety of 1-AG binds near the opening of the cyclooxygenase channel in the space vacated by the movement of the Leu-531 side chain, validating our previous hypothesis implicating the flexibility of the Leu-531 side chain as a determinant for the ability of COX-2 to oxygenate endocannabinoid substrates. Functional analyses carried out to compliment our structural findings indicated that Y355F and R513H muCOX-2 constructs had no effect on the oxygenation of 1-AG and 2-AG, whereas substitutions that resulted in a shortened side chain for Leu-531 had only modest effects. Both AA and 1-AG bind to R513H muCOX-2 in conformations similar to those observed in the co-crystal structures of these substrates with wild type enzyme.

Product-assisted Catalysis as the Basis of the Reaction Specificity of Threonine Synthase

Threonine synthase (TS), which is a pyridoxal 5'-phosphate (PLP)-dependent enzyme, catalyzes the elimination of the γ-phosphate group from O-phospho-L-homoserine (OPHS) and the subsequent addition of water at Cβ to form L-threonine. The catalytic course of TS is the most complex among the PLP enzymes, and it is an intriguing problem how the elementary steps are controlled in TS to carry out selective reactions. When L-vinylglycine was added to Thermus thermophilus HB8 TS in the presence of phosphate, L-threonine was formed with k(cat) and reaction specificity comparable with those when OPHS was used as the substrate. However, in the absence of phosphate or when sulfate was used in place of phosphate, only the side reaction product, α-ketobutyrate, was formed. Global analysis of the spectral changes in the reaction of TS with L-threonine showed that compared with the more acidic sulfate ion, the phosphate ion decreased the energy levels of the transition states of the addition of water at the Cβ of the PLP-α-aminocrotonate aldimine (AC) and the transaldimination to form L-threonine. The x-ray crystallographic analysis of TS complexed with an analog for AC gave a distinct electron density assigned to the phosphate ion derived from the solvent near the Cβ of the analog. These results indicated that the phosphate ion released from OPHS by γ-elimination acts as the base catalyst for the addition of water at Cβ of AC, thereby providing the basis of the reaction specificity. The phosphate ion is also considered to accelerate the protonation/deprotonation at Cγ.

Structural Basis of Drug Binding to CYP46A1, an Enzyme That Controls Cholesterol Turnover in the Brain

Cytochrome P450 46A1 (CYP46A1) initiates the major pathway of cholesterol elimination from the brain and thereby controls cholesterol turnover in this organ. We determined x-ray crystal structures of CYP46A1 in complex with four structurally distinct pharmaceuticals; antidepressant tranylcypromine (2.15 Å), anticonvulsant thioperamide (1.65 Å), antifungal voriconazole (2.35 Å), and antifungal clotrimazole (2.50 Å). All four drugs are nitrogen-containing compounds that have nanomolar affinity for CYP46A1 in vitro yet differ in size, shape, hydrophobicity, and type of the nitrogen ligand. Structures of the co-complexes demonstrate that each drug binds in a single orientation to the active site with tranylcypromine, thioperamide, and voriconazole coordinating the heme iron via their nitrogen atoms and clotrimazole being at a 4 Å distance from the heme iron. We show here that clotrimazole is also a substrate for CYP46A1. High affinity for CYP46A1 is determined by a set of specific interactions, some of which were further investigated by solution studies using structural analogs of the drugs and the T306A CYP46A1 mutant. Collectively, our results reveal how diverse inhibitors can be accommodated in the CYP46A1 active site and provide an explanation for the observed differences in the drug-induced spectral response. Co-complexes with tranylcypromine, thioperamide, and voriconazole represent the first structural characterization of the drug binding to a P450 enzyme.

SET7/9 Catalytic Mutants Reveal the Role of Active Site Water Molecules in Lysine Multiple Methylation

SET domain lysine methyltransferases (KMTs) methylate specific lysine residues in histone and non-histone substrates. These enzymes also display product specificity by catalyzing distinct degrees of methylation of the lysine ε-amino group. To elucidate the molecular mechanism underlying this specificity, we have characterized the Y245A and Y305F mutants of the human KMT SET7/9 (also known as KMT7) that alter its product specificity from a monomethyltransferase to a di- and a trimethyltransferase, respectively. Crystal structures of these mutants in complex with peptides bearing unmodified, mono-, di-, and trimethylated lysines illustrate the roles of active site water molecules in aligning the lysine ε-amino group for methyl transfer with S-adenosylmethionine. Displacement or dissociation of these solvent molecules enlarges the diameter of the active site, accommodating the increasing size of the methylated ε-amino group during successive methyl transfer reactions. Together, these results furnish new insights into the roles of active site water molecules in modulating lysine multiple methylation by SET domain KMTs and provide the first molecular snapshots of the mono-, di-, and trimethyl transfer reactions catalyzed by these enzymes.

Metal Ion Roles and the Movement of Hydrogen during Reaction Catalyzed by D-Xylose Isomerase: A Joint X-Ray and Neutron Diffraction Study

Conversion of aldo to keto sugars by the metalloenzyme D-xylose isomerase (XI) is a multistep reaction that involves hydrogen transfer. We have determined the structure of this enzyme by neutron diffraction in order to locate H atoms (or their isotope D). Two studies are presented, one of XI containing cadmium and cyclic D-glucose (before sugar ring opening has occurred), and the other containing nickel and linear D-glucose (after ring opening has occurred but before isomerization). Previously we reported the neutron structures of ligand-free enzyme and enzyme with bound product. The data show that His54 is doubly protonated on the ring N in all four structures. Lys289 is neutral before ring opening and gains a proton after this; the catalytic metal-bound water is deprotonated to hydroxyl during isomerization and O5 is deprotonated. These results lead to new suggestions as to how changes might take place over the course of the reaction.

Direct Determination of Protonation States of Histidine Residues in a 2 Å Neutron Structure of Deoxy-Human Normal Adult Hemoglobin and Implications for the Bohr Effect

We have investigated the protonation states of histidine residues (potential Bohr groups) in the deoxy form (T state) of human hemoglobin by direct determination of hydrogen (deuterium) positions with the neutron protein crystallography technique. The reversible binding of protons is key to the allosteric regulation of human hemoglobin. The protonation states of 35 of the 38 His residues were directly determined from neutron scattering omit maps, with 3 of the remaining residues being disordered. Protonation states of 5 equivalent His residues—αHis20, αHis50, αHis89, βHis143, and βHis146—differ between the symmetry-related globin subunits. The distal His residues, αHis58 and βHis63, are protonated in the α1β1 heterodimer and are neutral in α2β2. Buried residue αHis103 is found to be protonated in both subunits. These distal and buried residues have the potential to act as Bohr groups. The observed protonation states of His residues are compared to changes in their p K a values during the transition from the T to the R state and the results provide some new insights into our understanding of the molecular mechanism of the Bohr effect.

Stability of ligand-binding domain dimer assembly controls kainate receptor desensitization

AMPA and kainate receptors mediate fast synaptic transmission. AMPA receptor ligand-binding domains form dimers, which are key functional units controlling ion-channel activation and desensitization. Dimer stability is inversely related to the rate and extent of desensitization. Kainate and AMPA receptors share common structural elements, but functional measurements suggest that subunit assembly and gating differs between these subtypes. To investigate this, we constructed a library of GluR6 kainate receptor mutants and directly measured changes in kainate receptor dimer stability by analytical ultracentrifugation, which, combined with electrophysiological experiments, revealed an inverse correlation between dimer stability and the rate of desensitization. We solved crystal structures for a series of five GluR6 mutants, to understand the molecular mechanisms for dimer stabilization. We demonstrate that the desensitized state of kainate receptors acts as a deep energy well offsetting the stabilizing effects of dimer interface mutants, and that the deactivation of kainate receptor responses is dominated by entry into desensitized states. Our results show how neurotransmitter receptors with similar structures and gating mechanisms can exhibit strikingly different functional properties.