Identification of ions in experimental electrostatic potential maps

Exploiting the full potential of cryo-EM maps

The Coulomb potential maps generated by electron microscopy (EM) experiments contain not only information about the position but also about the charge state of the atom. This feature of EM maps allows the identification of specific ions and the protonation state of amino acid side chains in the sample. Here, we summarize qualitative observations of charges in EM maps, discuss the difficulties in interpreting the charge in Coulomb potential maps with respect to distinguishing it from radiation damage, and outline considerations to implement the correct charge in fitting algorithms.

Identification of potential inhibitors against E.coli via novel approaches based on deep learning and quantum mechanics-based atomistic investigations

Currently, drug resistance to commercially available antibiotics is imparting negative consequences to global health, and the development of novel antibiotics in a timely manner is a prime need of the hour. In the current study, an e-pharmacophore model was built using the 3D structure of DNA gyrase in complex with a standard inhibitor. The generated model was subjected to a pharmacophore based virtual screening against 45,257,086 molecules having 223,460,579 conformers available in MCULE database. Pharmacophore based screening retrieved eight molecules as top hit based on pharmacophoric features in comparison to standard inhibitors. Afterward, all eight compounds were subjected molecular docking based on deep learning algorithm. The molecular docking revealed that compound MCULE-6042843173 and MCULE-2362244223 had significant binding orientation inside active pocket of targeted protein with binding affinity of -9.52 and -9.24 kcal/mol respectively. In addition, density functional theory studies (DFT) were performed to evaluate quantum mechanics of top ranked compounds which were investigated through quantum mechanics (QM) computations which strongly assisted the findings of other in-silico investigations. Consequently, the MCULE-6042843173 and MCULE-2362244223 were subjected to MD simulation studies for evaluation of stability, hydrogen bond analysis, van der Waals interactions, and the contact profile of compounds with targeted amino acid residues. Findings of current study suggested MCULE-6042843173 and MCULE-2362244223 as potential and novel inhibitor of DNA Gyrase enzyme.

Structural insights into the mechanism of the sodium/iodide symporter

The sodium/iodide symporter (NIS) is the essential plasma membrane protein that mediates active iodide (I-) transport into the thyroid gland, the first step in the biosynthesis of the thyroid hormones-the master regulators of intermediary metabolism. NIS couples the inward translocation of I- against its electrochemical gradient to the inward transport of Na+ down its electrochemical gradient1,2. For nearly 50 years before its molecular identification3, NIS was the molecule at the centre of the single most effective internal radiation cancer therapy: radioiodide (131I-) treatment for thyroid cancer2. Mutations in NIS cause congenital hypothyroidism, which must be treated immediately after birth to prevent stunted growth and cognitive deficiency2. Here we report three structures of rat NIS, determined by single-particle cryo-electron microscopy: one with no substrates bound; one with two Na+ and one I- bound; and one with one Na+ and the oxyanion perrhenate bound. Structural analyses, functional characterization and computational studies show the substrate-binding sites and key residues for transport activity. Our results yield insights into how NIS selects, couples and translocates anions-thereby establishing a framework for understanding NIS function-and how it transports different substrates with different stoichiometries and releases substrates from its substrate-binding cavity into the cytosol.

How to correct relative voxel scale factors for calculations of vector-difference Fourier maps in cryo-EM

The atomic coordinates derived from cryo-electron microscopy (cryo-EM) maps can be inaccurate when the voxel scaling factors are not properly calibrated. Here, we describe a method for correcting relative voxel scaling factors between pairs of cryo-EM maps for the same or similar structures that are expanded or contracted relative to each other. We find that the correction of scaling factors reduces the amplitude differences of Fourier-inverted structure factors from voxel-rescaled maps by up to 20-30%, as shown by two cryo-EM maps of the SARS-CoV-2 spike protein measured at pH 4.0 and pH 8.0. This allows for the calculation of the difference map after properly scaling, revealing differences between the two structures for individual amino acid residues. Unexpectedly, the analysis uncovers two previously overlooked differences of amino acid residues in structures and their local structural changes. Furthermore, we demonstrate the method as applied to two cryo-EM maps of monomeric apo-photosystem II from the cyanobacteria Synechocystis sp. PCC 6803 and Thermosynechococcus elongatus. The resulting difference maps reveal many changes in the peripheral transmembrane PsbX subunit between the two species.

Structural Insights into Binding of Remdesivir Triphosphate within the Replication-Transcription Complex of SARS-CoV-2

Remdesivir is an adenosine analogue that has a cyano substitution in the C1' position of the ribosyl moiety and a modified base structure to stabilize the linkage of the base to the C1' atom with its strong electron-withdrawing cyano group. Within the replication-transcription complex (RTC) of SARS-CoV-2, the RNA-dependent RNA polymerase nsp12 selects remdesivir monophosphate (RMP) over adenosine monophosphate (AMP) for nucleotide incorporation but noticeably slows primer extension after the added RMP of the RNA duplex product is translocated by three base pairs. Cryo-EM structures have been determined for the RTC with RMP at the nucleotide-insertion (i) site or at the i + 1, i + 2, or i + 3 sites after product translocation to provide a structural basis for a delayed-inhibition mechanism by remdesivir. In this study, we applied molecular dynamics (MD) simulations to extend the resolution of structures to the measurable maximum that is intrinsically limited by MD properties of these complexes. Our MD simulations provide (i) a structural basis for nucleotide selectivity of the incoming substrates of remdesivir triphosphate over adenosine triphosphate and of ribonucleotide over deoxyribonucleotide, (ii) new detailed information on hydrogen atoms involved in H-bonding interactions between the enzyme and remdesivir, and (iii) direct information on the catalytically active complex that is not easily captured by experimental methods. Our improved resolution of interatomic interactions at the nucleotide-binding pocket between remedesivir and the polymerase could help to design a new class of anti-SARS-CoV-2 inhibitors.

Gastric proton pump with two occluded K+ engineered with sodium pump-mimetic mutations

The gastric H⁺,K⁺-ATPase mediates electroneutral exchange of 1H⁺/1K⁺ per ATP hydrolysed across the membrane. Previous structural analysis of the K⁺-occluded E2-P transition state of H⁺,K⁺-ATPase showed a single bound K⁺ at cation-binding site II, in marked contrast to the two K⁺ ions occluded at sites I and II of the closely-related Na⁺,K⁺-ATPase which mediates electrogenic 3Na⁺/2K⁺ translocation across the membrane. The molecular basis of the different K⁺ stoichiometry between these K⁺-counter-transporting pumps is elusive. We show a series of crystal structures and a cryo-EM structure of H⁺,K⁺-ATPase mutants with changes in the vicinity of site I, based on the structure of the sodium pump. Our step-wise and tailored construction of the mutants finally gave a two-K⁺ bound H⁺,K⁺-ATPase, achieved by five mutations, including amino acids directly coordinating K⁺ (Lys791Ser, Glu820Asp), indirectly contributing to cation-binding site formation (Tyr340Asn, Glu936Val), and allosterically stabilizing K⁺-occluded conformation (Tyr799Trp). This quintuple mutant in the K⁺-occluded E2-P state unambiguously shows two separate densities at the cation-binding site in its 2.6 Å resolution cryo-EM structure. These results offer new insights into how two closely-related cation pumps specify the number of K⁺ accommodated at their cation-binding site.

The gastric H⁺,K⁺-ATPase is a proton pump that creates the acidic environment of the stomach lumen, maintaining high proton gradient across the gastric mucosa cell membrane. Here, structural analysis of rationally designed H⁺,K⁺-ATPase mutants provides insight into this and other P-type ATPases cation binding stoichiometry and mechanisms.

Catalytic machinery of methane oxidation in particulate methane monooxygenase (pMMO)

In this focused review, we portray the recently reported 2.5 Å cyro-EM structure of the particulate methane monooxygenase (pMMO) from M. capsulatus (Bath). The structure of the functional holo-pMMO near atomic resolution has uncovered the sites of the copper cofactors including the location of the active site in the enzyme. The three coppers seen in the original X-ray crystal structures of the enzyme are now augmented by additional coppers in the transmembrane domain as well as in the water-exposed C-terminal subdomain of the PmoB subunit. The cryo-EM structure offers the first glimpse of the catalytic machinery capable of methane oxidation with high selectivity and efficiency. The findings are entirely consistent with the biochemical and biophysical findings previously reported in the literature, including the chemistry of hydrocarbon hydroxylation, regeneration of the catalyst for multiple turnovers, and the mechanism of aborting non-productive cycles to ensure kinetic competence.

Copper Centers in the Cryo-EM Structure of Particulate Methane Monooxygenase Reveal the Catalytic Machinery of Methane Oxidation

The particulate methane monooxygenase (pMMO) is the first enzyme in the C1 metabolic pathway in methanotrophic bacteria. As this enzyme converts methane into methanol efficiently near room temperature, it has become the paradigm for developing an understanding of this difficult C1 chemistry. pMMO is a membrane-bound protein with three subunits (PmoB, PmoA, and PmoC) and 12-14 coppers distributed among different sites. X-ray crystal structures that have revealed only three mononuclear coppers at three sites have neither disclosed the location of the active site nor the catalytic mechanism of the enzyme. Here we report a cyro-EM structure of holo-pMMO from Methylococcus capsulatus (Bath) at 2.5 Å, and develop quantitative electrostatic-potential profiling to scrutinize the nonprotein densities for signatures of the copper cofactors. Our results confirm a mononuclear Cu^I at the A site, resolve two Cu^Is at the B site, and uncover additional Cu^I clusters at the PmoA/PmoC interface within the membrane (D site) and in the water-exposed C-terminal subdomain of the PmoB (E clusters). These findings complete the minimal set of copper factors required for catalytic turnover of pMMO, offering a glimpse of the catalytic machinery for methane oxidation according to the chemical principles underlying the mechanism proposed earlier.

Structural analyses of an RNA stability element interacting with poly(A)

Cis-acting RNA elements are crucial for the regulation of polyadenylated RNA stability. The element for nuclear expression (ENE) contains a U-rich internal loop flanked by short helices. An ENE stabilizes RNA by sequestering the poly(A) tail via formation of a triplex structure that inhibits a rapid deadenylation-dependent decay pathway. Structure-based bioinformatic studies identified numerous ENE-like elements in evolutionarily diverse genomes, including a subclass containing two ENE motifs separated by a short double-helical region (double ENEs [dENEs]). Here, the structure of a dENE derived from a rice transposable element (TWIFB1) before and after poly(A) binding (∼24 kDa and ∼33 kDa, respectively) is investigated. We combine biochemical structure probing, small angle X-ray scattering (SAXS), and cryo-electron microscopy (cryo-EM) to investigate the dENE structure and its local and global structural changes upon poly(A) binding. Our data reveal 1) the directionality of poly(A) binding to the dENE, and 2) that the dENE-poly(A) interaction involves a motif that protects the 3'-most seven adenylates of the poly(A). Furthermore, we demonstrate that the dENE does not undergo a dramatic global conformational change upon poly(A) binding. These findings are consistent with the recently solved crystal structure of a dENE+poly(A) complex [S.-F. Torabi et al., Science 371, eabe6523 (2021)]. Identification of additional modes of poly(A)-RNA interaction opens new venues for better understanding of poly(A) tail biology.

Identification of Mg2+ ions next to nucleotides in cryo-EM maps using electrostatic potential maps

Cryo electron microscopy (cryo-EM) can produce maps of macromolecules that have resolutions that are sufficiently high that structural details such as chemical modifications, water molecules and bound metal ions can be discerned. However, those accustomed to interpreting the electron-density maps of macromolecules produced by X-ray crystallography need to be careful when assigning features such as these in cryo-EM maps because cations, for example, interact far more strongly with electrons than they do with X-rays. Using simulated electrostatic potential (ESP) maps as a tool led us to re-examine a recent cryo-EM map of the human ribosome, and we realized that some of the ESP peaks originally identified as novel groups covalently bonded to the N7, O6 or O4 atoms of several guanines, adenines or uridines, respectively, in this structure are likely to instead represent Mg2+ ions coordinated to these atoms, which provide only partial charge compensation compared with Mg2+ ions located next to phosphate groups. In addition, direct evidence is provided for a variation in the level of 2'-O ribose methylation of nucleotides in the human ribosome. ESP maps can thus help in identifying ions next to nucleotide bases, i.e. at positions that can be difficult to address in cryo-EM maps due to charge effects, which are specifically encountered in cryo-EM. This work is particularly relevant to nucleoprotein complexes and shows that it is important to consider charge effects when interpreting cryo-EM maps, thus opening possibilities for localizing charges in structures that may be relevant for enzymatic mechanisms and drug interactions.

Changes in Membrane Protein Structural Biology

Membrane proteins are essential components of many biochemical processes and are important pharmaceutical targets. Membrane protein structural biology provides the molecular rationale for these biochemical process as well as being a highly useful tool for drug discovery. Unfortunately, membrane protein structural biology is a difficult area of study due to low protein yields and high levels of instability especially when membrane proteins are removed from their native environments. Despite this instability, membrane protein structural biology has made great leaps over the last fifteen years. Today, the landscape is almost unrecognisable. The numbers of available atomic resolution structures have increased 10-fold though advances in crystallography and more recently by cryo-electron microscopy. These advances in structural biology were achieved through the efforts of many researchers around the world as well as initiatives such as the Membrane Protein Laboratory (MPL) at Diamond Light Source. The MPL has helped, provided access to and contributed to advances in protein production, sample preparation and data collection. Together, these advances have enabled higher resolution structures, from less material, at a greater rate, from a more diverse range of membrane protein targets. Despite this success, significant challenges remain. Here, we review the progress made and highlight current and future challenges that will be overcome.

Identification of a Na+-Binding Site near the Oxygen-Evolving Complex of Spinach Photosystem II

The oxygen-evolving complex (OEC) of photosystem II (PSII) is an oxomanganese cluster composed of four redox-active Mn ions and one redox-inactive Ca²⁺ ion, with two nearby bound Cl^- ions. Sodium is a common counterion of both chloride and hydroxide anions, and a sodium-specific binding site has not been identified near the OEC. Here, we find that the oxygen-evolution activity of spinach PSII increases with Na⁺ concentration, particularly at high pH. A Na⁺-specific binding site next to the OEC, becomes available after deprotonation of the D1-H337 amino acid residue, is suggested by the analysis of two recently published PSII cryo-electron microscopy maps in combination with quantum mechanical calculations and multiconformation continuum electrostatics simulations.

Cryo-EM structures reveal intricate Fe-S cluster arrangement and charging in Rhodobacter capsulatus formate dehydrogenase

Metal-containing formate dehydrogenases (FDH) catalyse the reversible oxidation of formate to carbon dioxide at their molybdenum or tungsten active site. They display a diverse subunit and cofactor composition, but structural information on these enzymes is limited. Here we report the cryo-electron microscopic structures of the soluble Rhodobacter capsulatus FDH (RcFDH) as isolated and in the presence of reduced nicotinamide adenine dinucleotide (NADH). RcFDH assembles into a 360 kDa dimer of heterotetramers revealing a putative interconnection of electron pathway chains. In the presence of NADH, the RcFDH structure shows charging of cofactors, indicative of an increased electron load.

Rhodobacter capsulatus NAD⁺ dependent formate dehydrogenase (RcFDH) is a molybdoenzyme that catalyses the reversible oxidation of formate to carbon dioxide, and is of interest for biotechnological applications. Here the authors present the cryo-EM structures of RcFDH as isolated from R. capsulatus and in the reduced state with bound NADH, and discuss the enzyme mechanism.

MicroED methodology and development

Microcrystal electron diffraction, or MicroED, is a method that is capable of determining structure from very small and thin 3D crystals using a transmission electron microscope. MicroED has been successfully used on microcrystalline samples, including proteins, peptides, and small organic molecules, in many cases to very high resolutions. In this work, the MicroED workflow will be briefly described and areas of future method development will be highlighted. These areas include improvements in sample preparation, data collection, and structure determination.

Late steps in bacterial translation initiation visualized using time-resolved cryo-EM

The initiation of bacterial translation involves the tightly regulated joining of the 50S ribosomal subunit to an initiator transfer RNA (fMet-tRNAfMet)-containing 30S ribosomal initiation complex to form a 70S initiation complex, which subsequently matures into a 70S elongation-competent complex. Rapid and accurate formation of the 70S initiation complex is promoted by initiation factors, which must dissociate from the 30S initiation complex before the resulting 70S elongation-competent complex can begin the elongation of translation1. Although comparisons of the structures of the 30S2-5 and 70S4,6-8 initiation complexes have revealed that the ribosome, initiation factors and fMet-tRNAfMet can acquire different conformations in these complexes, the timing of conformational changes during formation of the 70S initiation complex, the structures of any intermediates formed during these rearrangements, and the contributions that these dynamics might make to the mechanism and regulation of initiation remain unknown. Moreover, the absence of a structure of the 70S elongation-competent complex formed via an initiation-factor-catalysed reaction has precluded an understanding of the rearrangements to the ribosome, initiation factors and fMet-tRNAfMet that occur during maturation of a 70S initiation complex into a 70S elongation-competent complex. Here, using time-resolved cryogenic electron microscopy9, we report the near-atomic-resolution view of how a time-ordered series of conformational changes drive and regulate subunit joining, initiation factor dissociation and fMet-tRNAfMet positioning during formation of the 70S elongation-competent complex. Our results demonstrate the power of time-resolved cryogenic electron microscopy to determine how a time-ordered series of conformational changes contribute to the mechanism and regulation of one of the most fundamental processes in biology.

Reduced Occupancy of the Oxygen-Evolving Complex of Photosystem II Detected in Cryo-Electron Microscopy Maps

Computational simulations of electrostatic potentials (ESPs), based on atomistic models and independent atomic scattering factors, have remained challenging when applied to the oxygen-evolving complex (OEC) of photosystem II (PSII). Here, we overcome that challenge by using an ESP function obtained with density functional theory and atomic coordinates for the OEC of PSII obtained by optimization of the dark-adapted S₁ state. We find that the ESP is much higher for the OEC than for the nearby reference side chain of amino acid residue D1-H190. In contrast, experimental ESP maps recently published for two PSII-light-harvesting complex II super-complexes show that the ESP of the OEC is approximately half the value of the D1-H190 side chain. The apparent disparity is attributed to a reduced 31-38% occupancy of the OEC, likely associated with its reduction by electron scattering.