X-ray Structure Of Protein Research Articles

Multilevel superposition for deciphering the conformational variability of protein ensembles.

The dynamics and variability of protein conformations are directly linked to their functions. Many comparative studies of X-ray protein structures have been conducted to elucidate the relevant conformational changes, dynamics and heterogeneity. The rapid increase in the number of experimentally determined structures has made comparison an effective tool for investigating protein structures. For example, it is now possible to compare structural ensembles formed by enzyme species, variants or the type of ligands bound to them. In this study, the author developed a multilevel model for estimating two covariance matrices that represent inter- and intra-ensemble variability in the Cartesian coordinate space. Principal component analysis using the two estimated covariance matrices identified the inter-/intra-enzyme variabilities, which seemed to be important for the enzyme functions, with the illustrative examples of cytochrome P450 family 2 enzymes and class A $\beta$-lactamases. In P450, in which each enzyme has its own active site of a distinct size, an active-site motion shared universally between the enzymes was captured as the first principal mode of the intra-enzyme covariance matrix. In this case, the method was useful for understanding the conformational variability after adjusting for the differences between enzyme sizes. The developed method is advantageous in small ensemble-size problems and hence promising for use in comparative studies on experimentally determined structures where ensemble sizes are smaller than those generated, for example, by molecular dynamics simulations.

Universal pictures: A lithophane codex helps teenagers with blindness visualize nanoscopic systems.

People with blindness have limited access to the high-resolution graphical data and imagery of science. Here, a lithophane codex is reported. Its pages display tactile and optical readouts for universal visualization of data by persons with or without eyesight. Prototype codices illustrated microscopy of butterfly chitin-from N-acetylglucosamine monomer to fibril, scale, and whole insect-and were given to high schoolers from the Texas School for the Blind and Visually Impaired. Lithophane graphics of Fischer-Spier esterification reactions and electron micrographs of biological cells were also 3D-printed, along with x-ray structures of proteins (as millimeter-scale 3D models). Students with blindness could visualize (describe, recall, distinguish) these systems-for the first time-at the same resolution as sighted peers (average accuracy = 88%). Tactile visualization occurred alongside laboratory training, synthesis, and mentoring by chemists with blindness, resulting in increased student interest and sense of belonging in science.

X-ray structures of Enterobacter cloacae allose-binding protein in complexes with monosaccharides demonstrate its unique recognition mechanism for high affinity to allose

d-Allose is an aldohexose of the C3-epimer of d-glucose, existing in very small amounts in nature, called a rare sugar. The operon responsible for d-allose metabolism, the allose operon, was found in several bacteria, which consists of seven genes: alsR, alsB, alsA, alsC, alsE, alsK, and rpiB. To understand the biological implication of the allose operon utilizing a rare sugar of d-allose as a carbon source, it is important to clarify whether the allose operon functions specifically for d-allose or also functions for other ligands. It was proposed that the allose operon can function for d-ribose, which is essential as a component of nucleotides and abundant in nature. Allose-binding protein, AlsB, coded in the allose operon, is thought to capture a ligand outside the cell, and is expected to show high affinity for the specific ligand. X-ray structure determinations of Enterobacter cloacae AlsB (EtcAlsB) in ligand-free form, and in complexes with d-allose, d-ribose, and d-allulose, and measurements of the thermal parameters of the complex formation using an isothermal titration calorimeter were performed. The results demonstrated that EtcAlsB has a unique recognition mechanism for high affinity to d-allose by changing its conformation from an open to a closed form depending on d-allose-binding, and that the binding of d-ribose to EtcAlsB could not induce a completely closed form but an intermediate form, explaining the low affinity for d-ribose.

Predicting Mutation-Induced Changes in the Electronic Properties of Photosynthetic Proteins from First Principles: The Fenna-Matthews-Olson Complex Example.

Multiscale molecular modeling is utilized to predict optical absorption and circular dichroism spectra of two single-point mutants of the Fenna-Matthews-Olson photosynthetic pigment-protein complex. The modeling approach combines classical molecular dynamics simulations with structural refinement of photosynthetic pigments and calculations of their excited states in a polarizable protein environment. The only experimental input to the modeling protocol is the X-ray structure of the wild-type protein. The first-principles modeling reproduces changes in the experimental optical spectra of the considered mutants, Y16F and Q198V. Interestingly, the Q198V mutation has a negligible effect on the electronic properties of the targeted bacteriochlorophyll a pigment. Instead, the electronic properties of several other pigments respond to this mutation. The molecular modeling demonstrates that a single-point mutation can induce long-range effects on the protein structure, while extensive structural changes near a pigment do not necessarily lead to significant changes in the electronic properties of that pigment.

Modeling with Alternate Locations in X-ray Protein Structures.

In many molecular modeling applications, the standard procedure is still to handle proteins as single, rigid structures. While the importance of conformational flexibility is widely known, handling it remains challenging. Even the crystal structure of a protein usually contains variability exemplified in alternate side chain orientations or backbone segments. This conformational variability is encoded in PDB structure files by so-called alternate locations (AltLocs). Most modeling approaches either ignore AltLocs or resolve them with simple heuristics early on during structure import. We analyzed the occurrence and usage of AltLocs in the PDB and developed an algorithm to automatically handle AltLocs in PDB files enabling all structure-based methods using rigid structures to take the alternative protein conformations described by AltLocs into consideration. A respective software tool named AltLocEnumerator can be used as a structure preprocessor to easily exploit AltLocs. While the amount of data makes it difficult to show impact on a statistical level, handling AltLocs has a substantial impact on a case-by-case basis. We believe that the inspection and consideration of AltLocs is a very valuable approach in many modeling scenarios.

LSSmScarlet2 and LSSmScarlet3, Chemically Stable Genetically Encoded Red Fluorescent Proteins with a Large Stokes' Shift.

Red fluorescent proteins with a large Stokes’ shift (LSSRFPs) are genetically encoded and efficiently excited by 488 nm light, allowing simultaneous dual-color one- and two-photon fluorescence imaging and fluorescence correlation spectroscopy in combination with green fluorescent proteins FPs. Recently, based on the conventional bright mScarlet RFP, we developed the LSSRFP LSSmScarlet. LSSmScarlet is characterized by two pKa values at pH values of 1.9 and 5.8. In this study, we developed improved versions of LSSmScarlet, named LSSmScarlet2 and LSSmScarlet3, which are characterized by a Stokes’ shift of 128 nm and extreme pH stability with a single pKa value of 2.2. LSSmScarlet2 and LSSmScarlet3 had 1.8-fold faster and 3-fold slower maturation than LSSmScarlet, respectively. In addition, both LSSRFPs were 1.5- to 1.6-fold more photostable and more chemically resistant to denaturation by guanidinium chloride and guanidinium thiocyanate. We also compared the susceptibility of the LSSmScarlet2, LSSmScarlet3, and other LSSRFPs to the reagents used for whole-mount imaging, expansion microscopy, and immunostaining techniques. Due to higher pH stability and faster maturation, the LSSmScarlet3-LAMP3 fusion was 2.2-fold brighter than LSSmScarlet-LAMP3 in lysosomes of mammalian cells. The LSSmScarlet3-hLAMP2A fusion was similar in brightness to LSSmScarlet-hLAMP2A in lysosomes. We successfully applied the monomeric LSSmScarlet2 and LSSmScarlet3 proteins for confocal imaging of structural proteins in live mammalian cells. We also solved the X-ray structure of the LSSmScarlet2 protein at a resolution of 1.41 Å. Site-directed mutagenesis of the LSSmScarlet2 protein demonstrated the key role of the T74 residue in improving the pH and chemical stability of the LSSmScarlet2 protein.

In-silico Validation of the Proposed Treatment Strategy of Periodontitis.

The present study aims to assess a proposed treatment approach or therapy for periodontitis by using the in-silico technique. The proposed treatment strategy offers a singular vehicular system consisting of minocycline (antibiotic), celecoxib (selective COX-II inhibitor), doxycycline hyclate (matrix metalloproteinase inhibitor), and hydroxyapatite (osteogenic agent). Molecular docking studies of drugs were performed using Maestro version 9.4 software Schrödinger, and 3-Dimensional Crystallographic X-ray protein structures of targeted proteins were downloaded from RCSB protein data bank in .pdb file format. These agents were docked, and their affinities towards the receptors/protein/enzyme were calculated. Furthermore, their affinities were compared with the standard drug. The study suggests that minocycline and metronidazole possess equal affinity towards the RGPB and Inlj protein of P.gingivalis. Celecoxib, a well-known inhibitor of the COX-II enzyme, showed very high affinity. Selective inhibitor of MMP-8 possessed higher affinity than doxycycline, whereas CMT-3 showed equal affinity as doxycycline for MMP-13. Similarly, hydroxyapatite and simvastatin also showed a comparatively similar affinity for osteopontin receptor. Based upon molecular docking results, it can be concluded that the proposed treatment strategy would be a suitable approach for periodontitis and all the selected therapeutic agents have potential similar to the standard drugs, thereby constituting a reliable system for periodontitis.

Unmasking allosteric binding sites: Novel targets for GPCR drug discovery

ABSTRACT Introduction Unexpected non-apparent and hidden allosteric binding sites are non-classical and non-apparent allosteric centers in 3-D X-ray protein structures until orthosteric or allosteric ligands bind to them. The orthosteric center of one protomer that modulates binding centers of the other protomers within an oligomer is also an unexpected allosteric site. Furthermore, another partner protein can also produce these effects, acting as an unexpected allosteric modulator. Areas covered This review summarizes both classical and non-classical allosterism. The authors focus on G protein-coupled receptor (GPCR) oligomers as a paradigm of allosteric molecules. Moreover, they show several examples of unexpected allosteric sites such as hidden allosteric sites in a protomer that appear after the interaction with other molecules and the allosterism exerted between orthosteric sites within GPCR oligomer, emphasizing on the allosteric modulations that can occur between binding sites. Expert opinion The study of these new non-classical allosteric sites will expand the diversity of allosteric control on the function of orthosteric sites within proteins, whether GPCRs or other receptors, enzymes or transporters. Moreover, the design of new drugs targeting these hidden allosteric sites or already known orthosteric sites acting as allosteric sites in protein homo- or hetero-oligomers will increase the therapeutic potential of allosterism.

Mapping of Protein Binding Sites using clustering algorithms - Development of a pharmacophore based drug discovery tool

Discovering new hit small molecules binding to a specific protein binding site can be a difficult task. In support of existing procedures, a proof of concept methodology has been developed to process fragment flooded X-ray protein structures using the K-means clustering algorithm in order to derive pharmacophore models of the binding site. The novel method includes the implementation of several K-means initialisation methods in serial and parallel versions. Furthermore, required parameter optimisations for two initialisation methods was achieved, which was necessary to determine their validity and performance. A graph theory algorithm was adapted to compare the clustering-derived pharmacophores with X-ray ligand structure-derived pharmacophores to confirm that they mapped to each other. Initial proof of concept method validation was demonstrated using the Androgen Receptor (AR).

Solid state NMR chemical shift assignment of the non-structural single-stranded DNA binding protein gVp from fd bacteriophage.

The non-structural gene V protein (pV, gVp) from fd virus is a non-specific single-stranded DNA binding protein. The role of gVp is to sequester the single-stranded DNA thus reducing the generation of the replicative DNA form and leading to the formation of progeny phage. In this study, we assigned the 13C and 15N resonances of the crystalline unbound protein by magic-angle spinning solid-state NMR. The secondary structure predicted by the NMR shifts is in excellent agreement with the X-ray structure of the same 87-residue protein.

New Insight into the Effects of Various Parameters on the Crystallization of Ribulose-1,5-Bisphosphate Carboxylase/Oxygenase (RuBisCO) from Alcaligenes eutrophus

Crystallization remains a bottleneck for determining the three-dimensional X-ray structure of proteins. Many parameters influence the complexity of protein crystallization. Therefore, it is not easy to systematically examine all of these parameters individually during crystallization because of a limited quantity of purified protein. We studied several factors that influence crystallization including protein concentration, pH, temperature, age, volume of crystallization, inhibitors, metal ions, seeding, and precipitating agents on RuBisCO samples from Alcaligenes eutrophus which are not only freshly purified, but are also dissolved both individually and in combination from microcrystals and precipitated droplets of recycled RuBisCO. Single-, twin-, and/or microcrystals are dependent upon the concentration of RuBisCO by both RuBisCO samples. The morphology, either orthorhombic- or monoclinic-space group, depends upon pH. Furthermore, ammonium sulfate((NH4)2SO4) concentration at 20 °C (22% saturated) and/or at 4 °C (28% saturated) affected the crystallization of RuBisCO differently from one another. Finally, the age of RuBisCO also affected more uniformity and forming sharp edge during crystallization. Unexpected surprising monoclinic RuBisCO crystals were grown from dissolved microcrystals and precipitated droplets recycled RuBisCO samples. This quaternary RuBisCO single crystal, which contained Mg2+ and HCO3 for an activated ternary complex and is inhibited with a transition substrate analogue, CABP (2-carboxyarabinitol-1,5-bisphosphate)−, diffracts better than 2.2 Å. It is different from Hansen S. et al. reported RuBisCO crystals which were grown ab initio in absence of Mg2+, HCO3− and CABP, a structure which was determined at 2.7 Å resolution.

Divalent-Metal-Ion Selectivity of the CRISPR-Cas System-Associated Cas1 Protein: Insights from Classical Molecular Dynamics Simulations and Electronic Structure Calculations.

CRISPR-associated protein 1 (Cas1) is a universally conserved essential metalloenzyme of the clustered regularly interspaced short palindromic repeat (CRISPR) immune system of prokaryotes (bacteria, archaea) that can cut and integrate a part of viral DNA to its host genome with the help of other proteins. The integrated DNA acts as a memory of viral infection, which can be transcribed to RNA and stop future infection by recognition (based on the RNA/DNA complementarity principle) followed by protein-mediated degradation of the viral DNA. It has been proposed that the presence of a single manganese (Mn2+) ion in a conserved divalent-metal-ion binding pocket (key residues: E190, H254, D265, D268) of Cas1 is crucial for its function. Cas1-mediated DNA degradation was proposed to be hindered by metal substitution, metal chelation, or mutation of the binding pocket residues. Cas1 is active toward dsDNA degradation with both Mn2+ and Mg2+. X-ray structures of Cas1 revealed an intricate atomic interaction network of the divalent-metal-ion binding pocket and opened up the possibility of modeling related metal ions (viz., Mg2+, Ca2+) in the binding pocket of wild-type (WT) and mutated Cas1 proteins for computational analysis, which includes (1) quantitative estimation of the energetics of the divalent-metal-ion preference and (2) exploring the structural and dynamical aspects of the protein in response to divalent-metal-ion substitution or amino acid mutation. Using the X-ray structure of the Cas1 protein from Pseudomonas aeruginosa as a template (PDB 3GOD), we performed (∼2.23 μs) classical molecular dynamics (MD) simulations to compare structural and dynamical differences between Mg2+- and Ca2+-bound binding pockets of wild-type (WT) and mutant (E190A, H254A, D265A, D268A) Cas1. Furthermore, reduced binding pocket models were generated from X-ray and molecular dynamics (MD) trajectories, and the resulting structures were subjected to quantum chemical calculations. Results suggest that Cas1 prefers Mg2+ binding relative to Ca2+ and the preference is the strongest for WT and the weakest for the D268A mutant. Quantum chemical calculations indicate that Mn2+ is the most preferred relative to both Mg2+ and Ca2+ in the wild-type and mutant Cas1. Substitution of Mg2+ by Ca2+ does not alter the interaction network between Cas1 and the divalent metal ion but increases the wetness of the binding pocket by introducing a single water molecule in the first coordination shell of the latter. The strength of metal-ion preference (Mg2+ versus Ca2+) seems to be dependent on the solvent accessibility of the divalent-metal-ion binding pocket, strongest for wild-type Cas1 (in which the metal-ion binding pocket is dry, which includes two water molecules) and the weakest for the D268A mutant (in which the metal-ion binding pocket is wet, which includes four water molecules).

A computational methodology to diagnose sequence-variant dynamic perturbations by comparing atomic protein structures.

MotivationThe objective is to diagnose dynamics perturbations caused by amino acid mutations as prerequisite to assess protein functional health or drug failure, simply using network models of protein X-ray structures.ResultsWe find that the differences in the allocation of the atomic interactions of each amino acid to 1D, 2D, 3D, 4D structural levels between variants structurally robust, recover experimental dynamic perturbations. The allocation measure validated on two B-pentamers variants of AB5 toxins having 17 mutations, also distinguishes dynamic perturbations of pathogenic and non-pathogenic Transthyretin single-mutants. Finally, the main proteases of the coronaviruses SARS-CoV and SARS-CoV-2 exhibit changes in the allocation measure, raising the possibility of drug failure despite the main proteases structural similarity.AvailabilityThe Python code used for the production of the results is available at github.com/lorpac/protein_partitioning_atomic_contacts. The authors will run the analysis on any PDB structures of protein variants upon request.Supplementary informationSupplementary data are available at Bioinformatics online.

NMR structural analysis of the yeast cytochrome c oxidase subunit Cox13 and its interaction with ATP

BackgroundMitochondrial respiration is organized in a series of enzyme complexes in turn forming dynamic supercomplexes. In Saccharomyces cerevisiae (baker’s yeast), Cox13 (CoxVIa in mammals) is a conserved peripheral subunit of Complex IV (cytochrome c oxidase, CytcO), localized at the interface of dimeric bovine CytcO, which has been implicated in the regulation of the complex.ResultsHere, we report the solution NMR structure of Cox13, which forms a dimer in detergent micelles. Each Cox13 monomer has three short helices (SH), corresponding to disordered regions in X-ray or cryo-EM structures of homologous proteins. Dimer formation is mainly induced by hydrophobic interactions between the transmembrane (TM) helix of each monomer. Furthermore, an analysis of chemical shift changes upon addition of ATP revealed that ATP binds at a conserved region of the C terminus with considerable conformational flexibility.ConclusionsTogether with functional analysis of purified CytcO, we suggest that this ATP interaction is inhibitory of catalytic activity. Our results shed light on the structural flexibility of an important subunit of yeast CytcO and provide structure-based insight into how ATP could regulate mitochondrial respiration.

The Reaction Cycle of Dasscarboxylate Transporters

The divalent anion sodium symporter (DASS) family of plasma membrane transporters import di- and tricarboxylates essential to metabolic signaling and cellular regulation. While DASS proteins transport substrates via an elevator mechanism, previous structures were of only a single DASS cotransporter protein in a substrate-bound, inward-facing state. We have determined cryo-EM and X-ray structures of multiple DASS proteins in several different states, including three hitherto unseen states. Comparison of these structures reveal how the transport domain translates and rotates within the framework of the scaffold domain through the transport cycle. Additionally, we propose that DASS transporters ensure substrate coupling by a charge-compensation mechanism, and by structural changes upon substrate release.

Conservation of binding properties in protein models

We study the models submitted to round 12 of the Critical Assessment of protein Structure Prediction (CASP) experiment to assess how well the binding properties are conserved when the X-ray structures of the target proteins are replaced by their models. To explore small molecule binding we generate distributions of molecular probes – which are fragment-sized organic molecules of varying size, shape, and polarity – around the protein, and count the number of interactions between each residue and the probes, resulting in a vector of interactions we call a binding fingerprint. The similarity between two fingerprints, one for the X-ray structure and the other for a model of the protein, is determined by calculating the correlation coefficient between the two vectors. The resulting correlation coefficients are shown to correlate with global measures of accuracy established in CASP, and the relationship yields an accuracy threshold that has to be reached for meaningful binding surface conservation. The clusters formed by the probe molecules reliably predict binding hot spots and ligand binding sites in both X-ray structures and reasonably accurate models of the target, but ensembles of models may be needed for assessing the availability of proper binding pockets. We explored ligand docking to the few targets that had bound ligands in the X-ray structure. More targets were available to assess the ability of the models to reproduce protein–protein interactions by docking both the X-ray structures and models to their interaction partners in complexes. It was shown that this application is more difficult than finding small ligand binding sites, and the success rates heavily depend on the local structure in the potential interface. In particular, predicted conformations of flexible loops are frequently incorrect in otherwise highly accurate models, and may prevent predicting correct protein–protein interactions.

Structural characterization of HIV-1 matrix mutants implicated in envelope incorporation

During the late phase of HIV-1 infection, viral Gag polyproteins are targeted to the plasma membrane (PM) for assembly. Gag localization at the PM is a prerequisite for the incorporation of the envelope protein (Env) into budding particles. Gag assembly and Env incorporation are mediated by the N-terminal myristoylated matrix (MA) domain of Gag. Nonconservative mutations in the trimer interface of MA (A45E, T70R, and L75G) were found to impair Env incorporation and infectivity, leading to the hypothesis that MA trimerization is an obligatory step for Env incorporation. Conversely, Env incorporation can be rescued by a compensatory mutation in the MA trimer interface (Q63R). The impact of these MA mutations on the structure and trimerization properties of MA is not known. In this study, we employed NMR spectroscopy, X-ray crystallography, and sedimentation techniques to characterize the structure and trimerization properties of HIV-1 MA A45E, Q63R, T70R, and L75G mutant proteins. NMR data revealed that these point mutations did not alter the overall structure and folding of MA but caused minor structural perturbations in the trimer interface. Analytical ultracentrifugation data indicated that mutations had a minimal effect on the MA monomer–trimer equilibrium. The high-resolution X-ray structure of the unmyristoylated MA Q63R protein revealed hydrogen bonding between the side chains of adjacent Arg-63 and Ser-67 on neighboring MA molecules, providing the first structural evidence for an additional intermolecular interaction in the trimer interface. These findings advance our knowledge of the interplay of MA trimerization and Env incorporation into HIV-1 particles.

BANΔIT: B'-Factor Analysis for Drug Design and Structural Biology.

The analysis of B‐factor profiles from X‐ray protein structures can be utilized for structure‐based drug design since protein mobility changes have been associated with the quality of protein‐ligand interactions. With the BANΔIT (B’‐factor analysis and ΔB’ interpretation toolkit), we have developed a JavaScript‐based browser application that provides a graphical user interface for the normalization and analysis of B’‐factor profiles. To emphasize the usability for rational drug design applications, we have analyzed a selection of crystallographic protein‐ligand complexes and have given exemplary conclusions for further drug optimization including the development of a B’‐factor‐supported pharmacophore model for SARS CoV‐2 main protease inhibitors. BANΔIT is available online at https://bandit.uni‐mainz.de. The source code can be downloaded from https://github.com/FBarthels/BANDIT.

A3 adenosine receptor activation mechanisms: molecular dynamics analysis of inactive, active, and fully active states.

We investigated the Gi-coupled A3 adenosine receptor (A3AR) activation mechanism by running 7.2µs of molecular dynamics (MD) simulations. Based on homology to G protein-coupled receptor (GPCR) structures, three constitutively active mutant (CAM) and the wild-type (WT) A3ARs in the apo form were modeled. Conformational signatures associated with three different receptor states (inactive R, activeR*, and bound to Gi protein mimic) were predicted by analyzing and comparing the CAMs with WT receptor and by considering site-directed mutagenesis data. Detected signatures that were correlated with receptor state included: Persistent salt-bridges involving key charged residues for activation (including a novel, putative ionic lock), rotameric state of conserved W6.48, and Na+ ions and water molecules present. Active-coupled state signatures similar to the X-ray structures of β2 adrenergic receptor-Gs protein and A2AAR-mini-Gs and the recently solved cryo-EM A1AR-Gi complexes were found. Our MD analysis suggests that constitutive activation might arise from the D1073.49-R1083.50 ionic lock destabilization in R and the D1073.49-R1113.53 ionic lock stabilization in R* that presumably lowers the energy barrier associated with an R to R* transition. This study provides new opportunities to understand the underlying interactions of different receptor states of other Gi protein-coupled GPCRs.

Massive in Silico Study of Noble Gas Binding to the Structural Proteome.

Noble gases are chemically inert, and it was therefore thought they would have little effect on biology. Paradoxically, it was found that they do exhibit a wide range of biological effects, many of which are target-specific and potentially useful and some of which have been demonstrated in vivo. The underlying mechanisms by which useful pharmacology, such as tissue and neuroprotection, anti-addiction effects, and analgesia, is elicited are relatively unexplored. Experiments to probe the interactions of noble gases with specific proteins are more difficult with gases than those with other chemicals. It is clearly impractical to conduct the large number of gas-protein experiments required to gain a complete picture of noble gas biology. Given the simplicity of atoms as ligands, in silico methods provide an opportunity to gain insight into which noble gas-protein interactions are worthy of further experimental or advanced computational investigation. Our previous validation studies showed that in silico methods can accurately predict experimentally determined noble gas binding sites in X-ray structures of proteins. Here, we summarize the largest reported in silico reverse docking study involving 127 854 protein structures and the five nonradioactive noble gases. We describe how these computational screening methods are implemented, summarize the main types of interactions that occur between noble gases and target proteins, describe how the massive data set that this study generated can be analyzed (freely available at group18.csiro.au), and provide the NDMA receptor as an example of how these data can be used to understand the molecular pharmacology underlying the biology of the noble gases. We encourage chemical biologists to access the data and use them to expand the knowledge base of noble gas pharmacology, and to use this information, together with more efficient delivery systems, to develop "atomic drugs" that can fully exploit their considerable and relatively unexplored potential in medicine.