Search

Rapid protein-structure determination relies greatly on software that can automatically build a protein model into an experimental electron-density map. In favorable circumstances, various software systems are capable of building over 90% of the final model. However, completeness falls off rapidly with the resolution of the diffraction data. Manual completion of these partial models is usually feasible, but is time-consuming and prone to subjective interpretation. Except for the N- and C-termini of the chain, the end points of each missing fragment are known from the initial model. Hence, fitting fragments reduces to an inverse-kinematics problem. A method has been developed that combines fast inverse-kinematics algorithms with a real-space torsion-angle refinement procedure in a two-stage approach to fit missing main-chain fragments into the electron density between two anchor points. The first stage samples a large number of closing conformations, guided by the electron density. These candidates are ranked according to density fit. In a subsequent refinement stage, optimization steps are projected onto a carefully chosen subspace of conformation space to preserve rigid geometry and closure. Experimental results show that fitted fragments are in excellent agreement with the final refined structure for lengths of up to 12-15 residues in areas of weak or ambiguous electron density, even at medium to low resolution.

Functional mechanisms of biomolecules often manifest themselves precisely in transient conformational substates. Researchers have long sought to structurally characterize dynamic processes in non-coding RNA, combining experimental data with computer algorithms. However, adequate exploration of conformational space for these highly dynamic molecules, starting from static crystal structures, remains challenging. Here, we report a new conformational sampling procedure, KGSrna, which can efficiently probe the native ensemble of RNA molecules in solution. We found that KGSrna ensembles accurately represent the conformational landscapes of 3D RNA encoded by NMR proton chemical shifts. KGSrna resolves motionally averaged NMR data into structural contributions; when coupled with residual dipolar coupling data, a KGSrna ensemble revealed a previously uncharacterized transient excited state of the HIV-1 trans-activation response element stem–loop. Ensemble-based interpretations of averaged data can aid in formulating and testing dynamic, motion-based hypotheses of functional mechanisms in RNAs with broad implications for RNA engineering and therapeutic intervention.

Protein function often depends on the exchange between conformational substates. Allosteric ligand binding or distal mutations can stabilize specific active site conformations and consequently alter protein function. In addition to comparing independently determined X-ray crystal structures, alternative conformations observed at low levels of electron density have the potential to provide mechanistic insights into conformational dynamics. Here, we report a new multi-conformer contact network algorithm (CONTACT) that identifies networks of conformationally heterogeneous residues directly from high-resolution X-ray crystallography data. Contact networks in Escherichia coli dihydrofolate reductase (ecDHFR) predict the long-range pattern of NMR chemical shift perturbations of an allosteric mutation. A comparison of contact networks in wild type and mutant ecDHFR suggests how mutations that alter optimized networks of coordinated motions can impair catalytic function. Thus, CONTACT-guided mutagenesis will allow the structure-dynamics-function relationship to be exploited in protein engineering and design.

Modern protein crystal structures are based nearly exclusively on X-ray data collected at cryogenic temperatures (generally 100 K). The cooling process is thought to introduce little bias in the functional interpretation of structural results, because cryogenic temperatures minimally perturb the overall protein backbone fold. In contrast, here we show that flash cooling biases previously hidden structural ensembles in protein crystals. By analyzing available data for 30 different proteins using new computational tools for electron-density sampling, model refinement, and molecular packing analysis, we found that crystal cryocooling remodels the conformational distributions of more than 35% of side chains and eliminates packing defects necessary for functional motions. In the signaling switch protein, H-Ras, an allosteric network consistent with fluctuations detected in solution by NMR was uncovered in the room-temperature, but not the cryogenic, electron-density maps. These results expose a bias in structural databases toward smaller, overpacked, and unrealistically unique models. Monitoring room-temperature conformational ensembles by X-ray crystallography can reveal motions crucial for catalysis, ligand binding, and allosteric regulation.

Proteins must transition between multiple conformations to accomplish their functions in the cell, yet the dominant three-dimensional perspective involves a static crystal structure representing the most populated single conformation. Recent advances have revealed previously "hidden" alternative conformations with the potential to bridge dynamic descriptions from NMR spectroscopy with the static descriptions from X-ray crystallography. However, it remains challenging to determine which of these alternative conformations are collectively responsive to perturbations. As a result, the fundamental interplay between conformational heterogeneity and biological functions such as catalysis and allostery remains poorly understood. To overcome these obstacles, we develop new algorithms to identify hidden alternative conformations and test whether they form cooperative interaction networks across a set of benchmark proteins. We also combine these computational tools with novel X-ray crystallography experiments to characterize the inherent conformational heterogeneity of the diabetes therapeutic target PTP1B and identify new susceptible allosteric sites. Our approach will open new avenues for the design of novel therapeutics and dynamically enabled protein functions.View Large Image | View Hi-Res Image | Download PowerPoint Slide

The crystal structure of an essential bacterial protein, YeaZ, from T. maritima identifies an interface that potentially mediates protein–protein interaction., YeaZ is involved in a protein network that is essential for bacteria. The crystal structure of YeaZ from Thermotoga maritima was determined to 2.5 Å resolution. Although this protein belongs to a family of ancient actin-like ATPases, it appears that it has lost the ability to bind ATP since it lacks some key structural features that are important for interaction with ATP. A conserved surface was identified, supporting its role in the formation of protein complexes.

The crystal structure of an RmlC-type cupin with zinc acetate bound at the putative active site reveals significant differences from a previous structure without any bound ligand. The functional implications of the ligand-induced conformational changes are discussed., In the plant pathogen Xanthomonas campestris pv. campestris, the product of the tcmJ gene, XcTcmJ, encodes a protein belonging to the RmlC family of cupins. XcTcmJ was crystallized in a monoclinic space group (C2) in the presence of zinc acetate and the structure was determined to 1.6 Å resolution. Previously, the apo structure has been reported in the absence of any bound metal ion [Chin et al. (2006 ▶), Proteins, 65, 1046–1050]. The most significant difference between the apo structure and the structure of XcTcmJ described here is a reorganization of the binding site for zinc acetate, which was most likely acquired from the crystallization solution. This site is located in the conserved metal ion-binding domain at the putative active site of XcTcmJ. In addition, an acetate was also bound within coordination distance of the zinc. In order to accommodate this binding, rearrangement of a conserved histidine ligand is required as well as several nearby residues within and around the putative active site. These observations indicate that binding of zinc serves a functional role in this cupin protein.

The crystal structure of the highly specific γ-d-glutamyl-l-diamino acid endopeptidase YkfC from Bacillus cereus in complex with l-Ala-γ-d-Glu reveals the structural basis for the substrate specificity of NlpC/P60-family cysteine peptidases., Dipeptidyl-peptidase VI from Bacillus sphaericus and YkfC from Bacillus subtilis have both previously been characterized as highly specific γ-d-glutamyl-l-­diamino acid endopeptidases. The crystal structure of a YkfC ortholog from Bacillus cereus (BcYkfC) at 1.8 Å resolution revealed that it contains two N-terminal bacterial SH3 (SH3b) domains in addition to the C-terminal catalytic NlpC/P60 domain that is ubiquitous in the very large family of cell-wall-related cysteine peptidases. A bound reaction product (l-Ala-γ-d-Glu) enabled the identification of conserved sequence and structural signatures for recognition of l-Ala and γ-d-Glu and, therefore, provides a clear framework for understanding the substrate specificity observed in dipeptidyl-peptidase VI, YkfC and other NlpC/P60 domains in general. The first SH3b domain plays an important role in defining substrate specificity by contributing to the formation of the active site, such that only murein peptides with a free N-terminal alanine are allowed. A conserved tyrosine in the SH3b domain of the YkfC subfamily is correlated with the presence of a conserved acidic residue in the NlpC/P60 domain and both residues interact with the free amine group of the alanine. This structural feature allows the definition of a subfamily of NlpC/P60 enzymes with the same N-terminal substrate requirements, including a previously characterized cyanobacterial l-­alanine-γ-d-glutamate endopeptidase that contains the two key components (an NlpC/P60 domain attached to an SH3b domain) for assembly of a YkfC-like active site.

The crystal structures of two orthologous proteins from different Shewanella species have uncovered a resemblance to CRAL-TRIO carrier proteins, which suggest that they function as transporters of small nonpolar molecules. One protein adopts an open conformation, while the other adopts a closed structure that may act as a conformational switch in the transport of ligands at the membrane surface., The crystal structures of the proteins encoded by the YP_749275.1 and YP_001095227.1 genes from Shewanella frigidimarina and S. loihica, respectively, have been determined at 1.8 and 2.25 Å resolution, respectively. These proteins are members of a novel family of bacterial proteins that adopt the α/β SpoIIAA-like fold found in STAS and CRAL-TRIO domains. Despite sharing 54% sequence identity, these two proteins adopt distinct conformations arising from different dispositions of their α2 and α3 helices. In the ‘open’ conformation (YP_001095227.1), these helices are 15 Å apart, leading to the creation of a deep nonpolar cavity. In the ‘closed’ structure (YP_749275.1), the helices partially unfold and rearrange, occluding the cavity and decreasing the solvent-exposed hydrophobic surface. These two complementary structures are reminiscent of the conformational switch in CRAL-TRIO carriers of hydrophobic compounds. It is suggested that both proteins may associate with the lipid bilayer in their ‘open’ monomeric state by inserting their amphiphilic helices, α2 and α3, into the lipid bilayer. These bacterial proteins may function as carriers of nonpolar substances or as interfacially activated enzymes.

The crystal structure of the first representative of the Pfam PF07336 (DUF1470) family reveals a two-domain organization that contains a new fold, termed the ABATE domain, at the N-terminus and a treble-clef zinc finger that is likely to bind DNA at the C-terminus., The crystal structure of Jann_2411 from Jannaschia sp. strain CCS1, a member of the Pfam PF07336 family classified as a domain of unknown function (DUF1470), was solved to a resolution of 1.45 Å by multiple-wavelength anomalous dispersion (MAD). This protein is the first structural representative of the DUF1470 Pfam family. Structural analysis revealed a two-domain organization, with the N-terminal domain presenting a new fold called the ABATE domain that may bind an as yet unknown ligand. The C-terminal domain forms a treble-clef zinc finger that is likely to be involved in DNA binding. Analysis of the Jann_2411 protein and the broader ABATE-domain family suggests a role as stress-induced transcriptional regulators.

The crystal structure of tryptophanyl-tRNA synthetase from T. maritima unexpectedly revealed an iron–sulfur cluster bound to the tRNA anticodon-binding region., A novel aminoacyl-tRNA synthetase that contains an iron–sulfur cluster in the tRNA anticodon-binding region and efficiently charges tRNA with tryptophan has been found in Thermotoga maritima. The crystal structure of TmTrpRS (tryptophanyl-tRNA synthetase; TrpRS; EC 6.1.1.2) reveals an iron–sulfur [4Fe–­4S] cluster bound to the tRNA anticodon-binding (TAB) domain and an l-­tryptophan ligand in the active site. None of the other T. maritima aminoacyl-tRNA synthetases (AARSs) contain this [4Fe–4S] cluster-binding motif (C-x 22-C-x 6-C-x 2-C). It is speculated that the iron–sulfur cluster contributes to the stability of TmTrpRS and could play a role in the recognition of the anticodon.

The crystal structure of the BVU2987 gene product from B. vulgatus (UniProt A6L4L1) reveals that members of the new Pfam family PF11396 (domain of unknown function; DUF2874) are similar to β-lactamase inhibitor protein and YpmB., Proteins that contain the DUF2874 domain constitute a new Pfam family PF11396. Members of this family have predominantly been identified in microbes found in the human gut and oral cavity. The crystal structure of one member of this family, BVU2987 from Bacteroides vulgatus, has been determined, revealing a β-lactamase inhibitor protein-like structure with a tandem repeat of domains. Sequence analysis and structural comparisons reveal that BVU2987 and other DUF2874 proteins are related to β-lactamase inhibitor protein, PepSY and SmpA_OmlA proteins and hence are likely to function as inhibitory proteins.

NE1406, the first structural representative of PF09410, reveals a lipocalin-like fold with features that suggest involvement in lipid metabolism. In addition, NE1406 provides potential structural templates for two other protein families (PF07143 and PF08622)., The first structural representative of the domain of unknown function DUF2006 family, also known as Pfam family PF09410, comprises a lipocalin-like fold with domain duplication. The finding of the calycin signature in the N-­terminal domain, combined with remote sequence similarity to two other protein families (PF07143 and PF08622) implicated in isoprenoid metabolism and the oxidative stress response, support an involvement in lipid metabolism. Clusters of conserved residues that interact with ligand mimetics suggest that the binding and regulation sites map to the N-terminal domain and to the interdomain interface, respectively.

The crystal structure of the first representative of DUF364 family reveals a combination of enolase N-terminal-like and C-terminal Rossmann-like folds. Analysis of the interdomain cleft combined with sequence and genome context conservation among homologs, suggests a unique catalytic site likely involved in the synthesis of a flavin or pterin derivative., The crystal structure of Dhaf4260 from Desulfitobacterium hafniense DCB-2 was determined by single-wavelength anomalous diffraction (SAD) to a resolution of 2.01 Å using the semi-automated high-throughput pipeline of the Joint Center for Structural Genomics (JCSG) as part of the NIGMS Protein Structure Initiative (PSI). This protein structure is the first representative of the PF04016 (DUF364) Pfam family and reveals a novel combination of two well known domains (an enolase N-terminal-like fold followed by a Rossmann-like domain). Structural and bioinformatic analyses reveal partial similarities to Rossmann-like methyltransferases, with residues from the enolase-like fold combining to form a unique active site that is likely to be involved in the condensation or hydrolysis of molecules implicated in the synthesis of flavins, pterins or other siderophores. The genome context of Dhaf4260 and homologs additionally supports a role in heavy-metal chelation.

Structures of the first representatives of PF06684 (DUF1185) reveal a Bacillus chorismate mutase-like fold with a potential role in amino-acid synthesis., The crystal structures of BB2672 and SPO0826 were determined to resolutions of 1.7 and 2.1 Å by single-wavelength anomalous dispersion and multiple-wavelength anomalous dispersion, respectively, using the semi-automated high-throughput pipeline of the Joint Center for Structural Genomics (JCSG) as part of the NIGMS Protein Structure Initiative (PSI). These proteins are the first structural representatives of the PF06684 (DUF1185) Pfam family. Structural analysis revealed that both structures adopt a variant of the Bacillus chorismate mutase fold (BCM). The biological unit of both proteins is a hexamer and analysis of homologs indicates that the oligomer interface residues are highly conserved. The conformation of the critical regions for oligomerization appears to be dependent on pH or salt concentration, suggesting that this protein might be subject to environmental regulation. Structural similarities to BCM and genome-context analysis suggest a function in amino-acid synthesis.

The crystal structure of a putative NTP pyrophosphohydrolase, YP_001813558.1 from E. sibiricum, reveals a novel segment-swapped linked-dimer assembly., The crystal structure of a putative NTPase, YP_001813558.1 from Exiguo­bacterium sibiricum 255-15 (PF09934, DUF2166) was determined to 1.78 Å resolution. YP_001813558.1 and its homologs (dimeric dUTPases, MazG proteins and HisE-encoded phosphoribosyl ATP pyrophosphohydrolases) form a superfamily of all-α-helical NTP pyrophosphatases. In dimeric dUTPase-like proteins, a central four-helix bundle forms the active site. However, in YP_001813558.1, an unexpected intertwined swapping of two of the helices that compose the conserved helix bundle results in a ‘linked dimer’ that has not previously been observed for this family. Interestingly, despite this novel mode of dimerization, the metal-binding site for divalent cations, such as magnesium, that are essential for NTPase activity is still conserved. Furthermore, the active-site residues that are involved in sugar binding of the NTPs are also conserved when compared with other α-helical NTPases, but those that recognize the nucleotide bases are not conserved, suggesting a different substrate specificity.

The native state of a protein is regarded to be an ensemble of conformers, which allows association with binding partners. While some of this structural heterogeneity is retained upon crystallization, reliably extracting heterogeneous features from diffraction data has remained a challenge. In this study, a new algorithm for the automatic modelling of discrete hetero­geneity is presented. At high resolution, the authors’ single multi-conformer model, with correlated structural features to represent heterogeneity, shows improved agreement with the diffraction data compared with a single-conformer model. The model appears to be representative of the set of structures present in the crystal. In contrast, below 2 A resolution representing ambiguous electron density by correlated multi-conformers in a single model does not yield better agreement with the experimental data. Consistent with previous studies, this suggests that variability in multi-conformer models at lower resolution levels reflects uncertainty more than co­ordinated motion.

Crystal structures of two homologous peptidases from cyanobacteria Anabaena variabilis and Nostoc punctiforme at 1.05 A and 1.60 A resolution represent the first structures of a large class of cell-wall, cysteine peptidases that contain an N-terminal bacterial SH3-like domain (SH3b) and a C-terminal NlpC/P60 cysteine peptidase domain. The NlpC/P60 domain is a primitive, papain-like peptidase in the CA clan of cysteine peptidases with a Cys126/His176/His188 catalytic triad and a conserved catalytic core. We deduced from structure and sequence analysis, and then experimentally, that that these two proteins act as γ-D-glutamyl-L-diamino acid endopeptidases (EC 3.4.22.-). The active site is located near the interface between the SH3b and NlpC/P60 domains, where the SH3b domain may help define substrate specificity, instead of functioning as a targeting domain, so that only muropeptides with an N-terminal L-alanine can bind to the active site.

ECX21941 represents a very large family (over 600 members) of novel, ocean metagenome–specific proteins identified by clustering of the dataset from the Global Ocean Sampling expedition. The crystal structure of ECX21941 reveals unexpected similarity to Sm/LSm proteins, which are important RNA-binding proteins, despite no detectable sequence similarity. The ECX21941 protein assembles as a homopentamer in solution and in the crystal structure when expressed in Escherichia coli and represents the first pentameric structure for this Sm/LSm family of proteins, although the actual oligomeric form in vivo is currently not known. The genomic neighborhood analysis of ECX21941 and its homologs combined with sequence similarity searches suggest a cyanophage origin for this protein. The specific functions of members of this family are unknown, but our structure analysis of ECX21941 indicates nucleic acid-binding capabilities and suggests a role in RNA and/or DNA processing.

SsgA-like proteins (SALPs) are a family of homologous cell division-related proteins that occur exclusively in morphologically complex actinomycetes. We show that SsgB, a subfamily of SALPs, is the archetypal SALP that is functionally conserved in all sporulating actinomycetes. Sporulation-specific cell division of Streptomyces coelicolor ssgB mutants is restored by introduction of distant ssgB orthologues from other actinomycetes. Interestingly, the number of septa (and spores) of the complemented null mutants is dictated by the specific ssgB orthologue that is expressed. The crystal structure of the SsgB from Thermobifida fusca was determined at 2.6 Å resolution and represents the first structure for this family. The structure revealed similarities to a class of eukaryotic “whirly” single-stranded DNA/RNA-binding proteins. However, the electro-negative surface of the SALPs suggests that neither SsgB nor any of the other SALPs are likely to interact with nucleotide substrates. Instead, we show that a conserved hydrophobic surface is likely to be important for SALP function and suggest that proteins are the likely binding partners.

Regulatory inactivation of DnaA is dependent on Hda (homologous to DnaA), a protein homologous to the AAA+ (ATPases associated with diverse cellular activities) ATPase region of the replication initiator DnaA. When bound to the sliding clamp loaded onto duplex DNA, Hda can stimulate the transformation of active DnaA-ATP into inactive DnaA-ADP. The crystal structure of Hda from Shewanella amazonensis SB2B at 1.75 A resolution reveals that Hda resembles typical AAA+ ATPases. The arrangement of the two subdomains in Hda (residues 1-174 and 175-241) differs dramatically from that of DnaA. A CDP molecule anchors the Hda domains in a conformation that promotes dimer formation. The Hda dimer adopts a novel oligomeric assembly for AAA+ proteins in which the arginine finger, crucial for ATP hydrolysis, is fully exposed and available to hydrolyze DnaA-ATP through a typical AAA+ type of mechanism. The sliding clamp binding motifs at the N-terminus of each Hda monomer are partially buried and combine to form an antiparallel beta-sheet at the dimer interface. The inaccessibility of the clamp binding motifs in the CDP-bound structure of Hda suggests that conformational changes are required for Hda to form a functional complex with the clamp. Thus, the CDP-bound Hda dimer likely represents an inactive form of Hda.

BtDyP from Bacteroides thetaiotaomicron (strain VPI-5482) and TyrA from Shewanella oneidensis are dye-decolorizing peroxidases (DyPs), members of a new family of heme-dependent peroxidases recently identified in fungi and bacteria. Here, we report the crystal structures of BtDyP and TyrA at 1.6 and 2.7 Angstroms, respectively. BtDyP assembles into a hexamer, while TyrA assembles into a dimer; the dimerization interface is conserved between the two proteins. Each monomer exhibits a two-domain, {alpha}+{beta} ferredoxin-like fold. A site for heme binding was identified computationally, and modeling of a heme into the proposed active site allowed for identification of residues likely to be functionally important. Structural and sequence comparisons with other DyPs demonstrate a conservation of putative heme-binding residues, including an absolutely conserved histidine. Isothermal titration calorimetry experiments confirm heme binding, but with a stoichiometry of 0.3:1 (heme:protein).

Two efficient and complementary sampling algorithms are presented to explore the space of closed clash-free conformations of a flexible protein loop. The "seed sampling" algorithm samples conformations broadly distributed over this space, while the "deformation sampling" algorithm uses these conformations as starting points to explore more finely selected regions of the space. Computational results are shown for loops ranging from 5 to 25 residues. The algorithms are implemented in a toolkit, LoopTK, available at https://simtk.org/home/looptk.

Proteins must move between different conformations of their native ensemble to perform their functions. Crystal structures obtained from high-resolution X-ray diffraction data reflect this heterogeneity as a spatial and temporal conformational average. Although movement between natively populated alternative conformations can be critical for characterizing molecular mechanisms, it is challenging to identify these conformations within electron density maps. Alternative side chain conformations are generally well separated into distinct rotameric conformations, but alternative backbone conformations can overlap at several atomic positions. Our model building program qFit uses mixed integer quadratic programming (MIQP) to evaluate an extremely large number of combinations of sidechain conformers and backbone fragments to locally explain the electron density. Here, we describe two major modeling enhancements to qFit: peptide flips and alternative glycine conformations. We find that peptide flips fall into four stereotypical clusters and are enriched in glycine residues at the n+1 position. The potential for insights uncovered by new peptide flips and glycine conformations is exemplified by HIV protease, where different inhibitors are associated with peptide flips in the “flap” regions adjacent to the inhibitor binding site. Our results paint a picture of peptide flips as conformational switches, often enabled by glycine flexibility, that result in dramatic local rearrangements. Our results furthermore demonstrate the power of large-scale computational analysis to provide new insights into conformational heterogeneity. Overall, improved modeling of backbone heterogeneity with high-resolution X-ray data will connect dynamics to the structure-function relationship and help drive new design strategies for inhibitors of biomedically important systems., Author Summary Describing the multiple conformations of proteins is important for understanding the relationship between molecular flexibility and function. However, most methods for interpreting data from X-ray crystallography focus on building a single structure of the protein, which limits the potential for biological insights. Here we introduce an improved algorithm for using crystallographic data to model these multiple conformations that addresses two previously overlooked types of protein backbone flexibility: peptide flips and glycine movements. The method successfully models known examples of these types of multiple conformations, and also identifies new cases that were previously unrecognized but are well supported by the experimental data. For example, we discover glycine-driven peptide flips in the inhibitor-gating “flaps” of the drug target HIV protease that were not modeled in the original structures. Automatically modeling “hidden” multiple conformations of proteins using our algorithm may help drive biomedically relevant insights in structural biology pertaining to, e.g., drug discovery for HIV–1 protease and other therapeutic targets.

Proteins perform their function or interact with partners by exchanging between conformational substates on a wide range of spatiotemporal scales. Structurally characterizing these exchanges is challenging, both experimentally and computationally. Large, diffusional motions are often on timescales that are difficult to access with molecular dynamics simulations, especially for large proteins and their complexes. The low frequency modes of normal mode analysis (NMA) report on molecular fluctuations associated with biological activity. However, NMA is limited to a second order expansion about a minimum of the potential energy function, which limits opportunities to observe diffusional motions. By contrast, kino-geometric conformational sampling (KGS) permits large perturbations while maintaining the exact geometry of explicit conformational constraints, such as hydrogen bonds. Here, we extend KGS and show that a conformational ensemble of the α subunit Gαs of heterotrimeric stimulatory protein Gs exhibits structural features implicated in its activation pathway. Activation of protein Gs by G protein-coupled receptors (GPCRs) is associated with GDP release and large conformational changes of its α-helical domain. Our method reveals a coupled α-helical domain opening motion while, simultaneously, Gαs helix α5 samples an activated conformation. These motions are moderated in the activated state. The motion centers on a dynamic hub near the nucleotide-binding site of Gαs, and radiates to helix α4. We find that comparative NMA-based ensembles underestimate the amplitudes of the motion. Additionally, the ensembles fall short in predicting the accepted direction of the full activation pathway. Taken together, our findings suggest that nullspace sampling with explicit, holonomic constraints yields ensembles that illuminate molecular mechanisms involved in GDP release and protein Gs activation, and further establish conformational coupling between key structural elements of Gαs., Author Summary Multi-cellular physiology is an emergent property, which depends critically on inter-cellular signaling pathways. Transmembrane G protein-coupled receptors (GPCRs) mediate a large variety of physiological events throughout the body, such as vision or cardiovascular regulation. It is thus no surprise that GPCRs are targeted by more than one third of all FDA-approved drugs. Molecules such as hormones and neurotransmitters transmit messages to cells via GPCRs complexed to cytosolic heterotrimeric G proteins. G proteins, upon activation, interact with other molecules to trigger a cellular response. Despite an increasing amount of structural data, the precise conformational dynamics and activation mechanism of G proteins remain poorly understood. The size of the multi-protein complexes and the time scales at which conformational changes occur hinder adequate sampling of the conformational landscape with molecular dynamics simulations. Here, we extend and use an efficient, robotics-inspired conformational sampling procedure to probe the conformational landscape of protein G during activation. Our procedure reveals coupled, molecular mechanisms of the activation pathway, which are absent in a comparative analysis with normal modes. Our exciting results can ultimately lead to modulation of biological activity by drug design or fine-tuning of conformational heterogeneity.

Glutathione S-transferases (GSTs) comprise a diverse superfamily of enzymes found in organisms from all kingdoms of life. GSTs are involved in diverse processes, notably small-molecule biosynthesis or detoxification, and are frequently also used in protein engineering studies or as biotechnology tools. Here, we report the high-resolution X-ray structure of Atu5508 from the pathogenic soil bacterium Agrobacterium tumefaciens (atGST1). Through use of comparative sequence and structural analysis of the GST superfamily, we identified local sequence and structural signatures, which allowed us to distinguish between different GST classes. This approach enables GST classification based on structure, without requiring additional biochemical or immunological data. Consequently, analysis of the atGST1 crystal structure suggests a new GST class, distinct from previously characterized GSTs, which would make it an attractive target for further biochemical studies.

Michael DiDonato, S. Sri Krishna, Robert Schwarzenbacher, Daniel McMullan, Sanjay Agarwalla, Scott M. Brittain, Mitchell D. Miller, Polat Abdubek, Eileen Ambing, Herbert L. Axelrod, JaumeM. Canaves, Hsiu-Ju Chiu, Ashley M. Deacon, Lian Duan, Marc-Andre Elsliger, Adam Godzik, Slawomir K. Grzechnik, Joanna Hale, Eric Hampton, Justin Haugen, Lukasz Jaroszewski, Kevin K. Jin, Heath E. Klock, Mark W. Knuth, Eric Koesema, Andreas Kreusch, Peter Kuhn, Scott A. Lesley, Inna Levin, Andrew T. Morse, Edward Nigoghossian, Linda Okach, Silvya Oommachen, Jessica Paulsen, Kevin Quijano, Ron Reyes, Christopher L. Rife, Glen Spraggon, Raymond C. Stevens, Henry van den Bedem, AprilfawnWhite, Guenter Wolf, Qingping Xu, Keith O. Hodgson, John Wooley, and Ian A. Wilson* The Joint Center for Structural Genomics Stanford Synchrotron Radiation Laboratory, Stanford University, Menlo Park, California The University of California, San Diego, La Jolla, California The Genomics Institute of the Novartis Research Foundation, San Diego, California The Scripps Research Institute, La Jolla, California