Automated electron-density sampling reveals widespread conformational polymorphism in proteins

Identifying and avoiding radiation damage in macromolecular crystallography

Radiation damage remains one of the major impediments to accurate structure solution in macromolecular crystallography. The artefacts of radiation damage can manifest as structural changes that result in incorrect biological interpretations being drawn from a model, they can reduce the resolution to which data can be collected and they can even prevent structure solution entirely. In this article, we discuss how to identify and mitigate against the effects of radiation damage at each stage in the macromolecular crystal structure-solution pipeline.

Alternate conformations found in protein structures implies biological functions: A case study using cyclophilin A

Protein dynamics linked to numerous biomolecular functions, such as ligand binding, allosteric regulation, and catalysis, must be better understood at the atomic level. Reactive atoms of key residues drive a repertoire of biomolecular functions by flipping between alternate conformations or conformational substates, seldom found in protein structures. Probing such sparsely sampled alternate conformations would provide mechanistic insight into many biological functions. We are therefore interested in evaluating the instance of amino acids adopted alternate conformations, either in backbone or side-chain atoms or in both. Accordingly, over 70000 protein structures appear to contain alternate conformations only 'A' and 'B' for any atom, particularly the instance of amino acids that adopted alternate conformations are more for Arg, Cys, Met, and Ser than others. The resulting protein structure analysis depicts that amino acids with alternate conformations are mainly found in the helical and β-regions and are often seen in high-resolution X-ray crystal structures. Furthermore, a case study on human cyclophilin A (CypA) was performed to explain the pre-existing intrinsic dynamics of catalytically critical residues from the CypA and how such intrinsic dynamics perturbed upon Ser99Thr mutation using molecular dynamics simulations on the ns-μs timescale. Simulation results demonstrated that the Ser99Thr mutation had impaired the alternate conformations or the catalytically productive micro-environment of Phe113, mimicking the experimentally observed perturbation captured by X-ray crystallography. In brief, a deeper comprehension of alternate conformations adopted by the amino acids may shed light on the interplay between protein structure, dynamics, and function.

FLEXR GUI: a graphical user interface for multi-conformer modeling of proteins

Proteins are well known 'shapeshifters' which change conformation to function. In crystallography, multiple conformational states are often present within the crystal and the resulting electron-density map. Yet, explicitly incorporating alternative states into models to disentangle multi-conformer ensembles is challenging. We previously reported the tool FLEXR, which, within a few minutes, automatically separates conformational signal from noise and builds the corresponding, often missing, structural features into a multi-conformer model. To make the method widely accessible for routine multi-conformer building as part of the computational toolkit for macromolecular crystallography, we present a graphical user interface (GUI) for FLEXR, designed as a plugin for Coot 1. The GUI implementation seamlessly connects FLEXR models with the existing suite of validation and modeling tools available in Coot. We envision that FLEXR will aid crystallographers by increasing access to a multi-conformer modeling method that will ultimately lead to a better representation of protein conformational heterogeneity in the Protein Data Bank. In turn, deeper insights into the protein conformational landscape may inform biology or provide new opportunities for ligand design. The code is open source and freely available on GitHub at https://github.com/TheFischerLab/FLEXR-GUI.

Structural characterization of ligand binding and pH-specific enzymatic activity of mouse Acidic Mammalian Chitinase

Chitin is an abundant biopolymer and pathogen-associated molecular pattern that stimulates a host innate immune response. Mammals express chitin-binding and chitin-degrading proteins to remove chitin from the body. One of these proteins, Acidic Mammalian Chitinase (AMCase), is an enzyme known for its ability to function under acidic conditions in the stomach but is also active in tissues with more neutral pHs, such as the lung. Here, we used a combination of biochemical, structural, and computational modeling approaches to examine how the mouse homolog (mAMCase) can act in both acidic and neutral environments. We measured kinetic properties of mAMCase activity across a broad pH range, quantifying its unusual dual activity optima at pH 2 and 7. We also solved high resolution crystal structures of mAMCase in complex with oligomeric GlcNAcn, the building block of chitin, where we identified extensive conformational ligand heterogeneity. Leveraging these data, we conducted molecular dynamics simulations that suggest how a key catalytic residue could be protonated via distinct mechanisms in each of the two environmental pH ranges. These results integrate structural, biochemical, and computational approaches to deliver a more complete understanding of the catalytic mechanism governing mAMCase activity at different pH. Engineering proteins with tunable pH optima may provide new opportunities to develop improved enzyme variants, including AMCase, for therapeutic purposes in chitin degradation.

Pushed to extremes: distinct effects of high temperature versus pressure on the structure of STEP

Protein function hinges on small shifts of three-dimensional structure. Elevating temperature or pressure may provide experimentally accessible insights into such shifts, but the effects of these distinct perturbations on protein structures have not been compared in atomic detail. To quantitatively explore these two axes, we report the first pair of structures at physiological temperature versus. high pressure for the same protein, STEP (PTPN5). We show that these perturbations have distinct and surprising effects on protein volume, patterns of ordered solvent, and local backbone and side-chain conformations. This includes interactions between key catalytic loops only at physiological temperature, and a distinct conformational ensemble for another active-site loop only at high pressure. Strikingly, in torsional space, physiological temperature shifts STEP toward previously reported active-like states, while high pressure shifts it toward a previously uncharted region. Altogether, our work indicates that temperature and pressure are complementary, powerful, fundamental macromolecular perturbations.

Ligand 1 H NMR Chemical Shifts as Accurate Reporters for Protein-Ligand Binding Interfaces in Solution

The availability of high-resolution 3D structural information is crucial for investigating guest-host systems across a wide range of fields. In the context of drug discovery, the information is routinely used to establish and validate structure-activity relationships, grow initial hits from screening campaigns, and to guide molecular docking. For the generation of protein-ligand complex structural information, X-ray crystallography is the experimental method of choice, however, with limited information on protein flexibility. An experimentally verified structural model of the binding interface in the native solution-state would support medicinal chemists in their molecular design decisions. Here we demonstrate that protein-bound ligand 1 H NMR chemical shifts are highly sensitive and accurate probes for the immediate chemical environment of protein-ligand interfaces. By comparing the experimental ligand 1 H chemical shift values with those computed from the X-ray structure using quantum mechanics methodology, we identify significant disagreements for parts of the ligand between the two experimental techniques. We show that quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) ensembles can be used to refine initial X-ray co-crystal structures resulting in a better agreement with experimental 1 H ligand chemical shift values. Overall, our findings highlight the usefulness of ligand 1 H NMR chemical shift information in combination with a QM/MM MD workflow for generating protein-ligand ensembles that accurately reproduce solution structural data.

High-resolution double vision of the allosteric phosphatase PTP1B

Protein tyrosine phosphatase 1B (PTP1B) plays important roles in cellular homeostasis and is a highly validated therapeutic target for multiple human ailments, including diabetes, obesity and breast cancer. However, much remains to be learned about how conformational changes may convey information through the structure of PTP1B to enable allosteric regulation by ligands or functional responses to mutations. High-resolution X-ray crystallography can offer unique windows into protein conformational ensembles, but comparison of even high-resolution structures is often complicated by differences between data sets, including non-isomorphism. Here, the highest resolution crystal structure of apo wild-type (WT) PTP1B to date is presented out of a total of ∼350 PTP1B structures in the PDB. This structure is in a crystal form that is rare for PTP1B, with two unique copies of the protein that exhibit distinct patterns of conformational heterogeneity, allowing a controlled comparison of local disorder across the two chains within the same asymmetric unit. The conformational differences between these chains are interrogated in the apo structure and between several recently reported high-resolution ligand-bound structures. Electron-density maps in a high-resolution structure of a recently reported activating double mutant are also examined, and unmodeled alternate conformations in the mutant structure are discovered that coincide with regions of enhanced conformational heterogeneity in the new WT structure. These results validate the notion that these mutations operate by enhancing local dynamics, and suggest a latent susceptibility to such changes in the WT enzyme. Together, these new data and analysis provide a detailed view of the conformational ensemble of PTP1B and highlight the utility of high-resolution crystallography for elucidating conformational heterogeneity with potential relevance for function.

Combining temperature perturbations with X-ray crystallography to study dynamic macromolecules: A thorough discussion of experimental methods

Temperature is an important state variable that governs the behavior of microscopic systems, yet crystallographers rarely exploit temperature changes to study the structure and dynamics of biological macromolecules. In fact, approximately 90% of crystal structures in the Protein Data Bank were determined under cryogenic conditions, because sample cryocooling makes crystals robust to X-ray radiation damage and facilitates data collection. On the other hand, cryocooling can introduce artifacts into macromolecular structures, and can suppress conformational dynamics that are critical for function. Fortunately, recent advances in X-ray detector technology, X-ray sources, and computational data processing algorithms make non-cryogenic X-ray crystallography easier and more broadly applicable than ever before. Without the reliance on cryocooling, high-resolution crystallography can be combined with various temperature perturbations to gain deep insight into the conformational landscapes of macromolecules. This Chapter reviews the historical reasons for the prevalence of cryocooling in macromolecular crystallography, and discusses its potential drawbacks. Next, the Chapter summarizes technological developments and methodologies that facilitate non-cryogenic crystallography experiments. Finally, the chapter discusses the theoretical underpinnings and practical aspects of multi-temperature and temperature-jump crystallography experiments, which are powerful tools for understanding the relationship between the structure, dynamics, and function of proteins and other biological macromolecules.

MD simulations of macromolecular crystals: Implications for the analysis of Bragg and diffuse scattering

Some of our most detailed information about structure and dynamics of macromolecules comes from X-ray-diffraction studies in crystalline environments. More than 170,000 atomic models have been deposited in the Protein Data Bank, and the number of observations (typically of intensities of Bragg diffraction peaks) is generally quite large, when compared to other experimental methods. Nevertheless, the general agreement between calculated and observed intensities is far outside the experimental precision, and the majority of scattered photons fall between the sharp Bragg peaks, and are rarely taken into account. This chapter considers how molecular dynamics simulations can be used to explore the connections between microscopic behavior in a crystalline lattice and observed scattering intensities, and point the way to new atomic models that could more faithfully recapitulate Bragg intensities and extract useful information from the diffuse scattering that lies between those peaks.

Modeling diffuse scattering with simple, physically interpretable models

Diffuse scattering has long been proposed to probe protein dynamics relevant for biological function, and more recently, as a tool to aid structure determination. Despite recent advances in measuring and modeling this signal, the field has not been able to routinely use experimental diffuse scattering for either application. A persistent challenge has been to devise models that are sophisticated enough to robustly reproduce experimental diffuse features but remain readily interpretable from the standpoint of structural biology. This chapter presents eryx, a suite of computational tools to evaluate the primary models of disorder that have been used to analyze protein diffuse scattering. By facilitating comparative modeling, eryx aims to provide insights into the physical origins of this signal and help identify the sources of disorder that are critical for reproducing experimental features. This framework also lays the groundwork for the development of more advanced models that integrate different types of disorder without loss of interpretability.

Refinement of multiconformer ensemble models from multi-temperature X-ray diffraction data

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363 K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

Refinement of Multiconformer Ensemble Models from Multi-temperature X-ray Diffraction Data

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

Pushed to extremes: distinct effects of high temperature vs. pressure on the structure of an atypical phosphatase

Protein function hinges on small shifts of three-dimensional structure. Elevating temperature or pressure may provide experimentally accessible insights into such shifts, but the effects of these distinct perturbations on protein structures have not been compared in atomic detail. To quantitatively explore these two axes, we report the first pair of structures at physiological temperature vs. high pressure for the same protein, STEP (PTPN5). We show that these perturbations have distinct and surprising effects on protein volume, patterns of ordered solvent, and local backbone and side-chain conformations. This includes novel interactions between key catalytic loops only at physiological temperature, and a distinct conformational ensemble for another active-site loop only at high pressure. Strikingly, in torsional space, physiological temperature shifts STEP toward previously reported active-like states, while high pressure shifts it toward a previously uncharted region. Together, our work argues that temperature and pressure are complementary, powerful, fundamental macromolecular perturbations.

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.

Obtaining anomalous and ensemble information from protein crystals from 220 K up to physiological temperatures

X-ray crystallography has been invaluable in delivering structural information about proteins. Previously, an approach has been developed that allows high-quality X-ray diffraction data to be obtained from protein crystals at and above room temperature. Here, this previous work is built on and extended by showing that high-quality anomalous signal can be obtained from single protein crystals using diffraction data collected at 220 K up to physiological temperatures. The anomalous signal can be used to directly determine the structure of a protein, i.e. to phase the data, as is routinely performed under cryoconditions. This ability is demonstrated by obtaining diffraction data from model lysozyme, thaumatin and proteinase K crystals, the anomalous signal from which allowed their structures to be solved experimentally at 7.1 keV X-ray energy and at room temperature with relatively low data redundancy. It is also demonstrated that the anomalous signal from diffraction data obtained at 310 K (37°C) can be used to solve the structure of proteinase K and to identify ordered ions. The method provides useful anomalous signal at temperatures down to 220 K, resulting in an extended crystal lifetime and increased data redundancy. Finally, we show that useful anomalous signal can be obtained at room temperature using X-rays of 12 keV energy as typically used for routine data collection, allowing this type of experiment to be carried out at widely accessible synchrotron beamline energies and enabling the simultaneous extraction of high-resolution data and anomalous signal. With the recent emphasis on obtaining conformational ensemble information for proteins, the high resolution of the data allows such ensembles to be built, while the anomalous signal allows the structure to be experimentally solved, ions to be identified, and water molecules and ions to be differentiated. Because bound metal-, phosphorus- and sulfur-containing ions all have anomalous signal, obtaining anomalous signal across temperatures and up to physiological temperatures will provide a more complete description of protein conformational ensembles, function and energetics.

Room-temperature serial synchrotron crystallography of the human phosphatase PTP1B

Room-temperature X-ray crystallography provides unique insights into protein conformational heterogeneity, but obtaining sufficiently large protein crystals is a common hurdle. Serial synchrotron crystallography (SSX) helps to address this hurdle by allowing the use of many medium- to small-sized crystals. Here, a recently introduced serial sample-support chip system has been used to obtain the first SSX structure of a human phosphatase, specifically protein tyrosine phosphatase 1B (PTP1B) in the unliganded (apo) state. In previous apo room-temperature structures, the active site and allosteric sites adopted alternate conformations, including open and closed conformations of the active-site WPD loop and of a distal allosteric site. By contrast, in our SSX structure the active site is best fitted with a single conformation, but the distal allosteric site is best fitted with alternate conformations. This observation argues for additional nuance in interpreting the nature of allosteric coupling in this protein. Overall, our results illustrate the promise of serial methods for room-temperature crystallography, as well as future avant-garde crystallography experiments, for PTP1B and other proteins.

Determining biomolecular structures near room temperature using X-ray crystallography: concepts, methods and future optimization

For roughly two decades, cryocrystallography has been the overwhelmingly dominant method for determining high-resolution biomolecular structures. Competition from single-particle cryo-electron microscopy and micro-electron diffraction, increased interest in functionally relevant information that may be missing or corrupted in structures determined at cryogenic temperature, and interest in time-resolved studies of the biomolecular response to chemical and optical stimuli have driven renewed interest in data collection at room temperature and, more generally, at temperatures from the protein-solvent glass transition near 200 K to ∼350 K. Fischer has recently reviewed practical methods for room-temperature data collection and analysis [Fischer (2021), Q. Rev. Biophys. 54, e1]. Here, the key advantages and physical principles of, and methods for, crystallographic data collection at noncryogenic temperatures and some factors relevant to interpreting the resulting data are discussed. For room-temperature data collection to realize its potential within the structural biology toolkit, streamlined and standardized methods for delivering crystals prepared in the home laboratory to the synchrotron and for automated handling and data collection, similar to those for cryocrystallography, should be implemented.

Computational analysis of protein conformational heterogeneity

In this paper, we applied the molecular dynamics (MD) simulations and used thermolysin as the system to study the overall protein dynamics and side chain dihedral angles across the Arrhenius break. Simulations were performed at a high temperature of 36 °C which is above the previously observed Arrhenius break, and the lower temperature of 20 °C which is below the Arrhenius break. We observed different protein dynamics and conformational heterogeneity of side chain dihedral angles of thermolysin at the two temperatures. Our results indicated that certain regions of thermolysin have a higher level of fluctuation at lower temperature. A temperature dependent dihedral angles were also observed at the two temperatures. The changes observed in the protein indicated key areas of temperature sensitivity within the protein.Communicated by Ramaswamy H. Sarma.

Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery

EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 μM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 μM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.

Exploiting Sequence-Dependent Rotamer Information in Global Optimization of Proteins

Rotamers, namely amino acid side chain conformations common to many different peptides, can be compiled into libraries. These rotamer libraries are used in protein modeling, where the limited conformational space occupied by amino acid side chains is exploited. Here, we construct a sequence-dependent rotamer library from simulations of all possible tripeptides, which provides rotameric states dependent on adjacent amino acids. We observe significant sensitivity of rotamer populations to sequence and find that the library is successful in locating side chain conformations present in crystal structures. The library is designed for applications with basin-hopping global optimization, where we use it to propose moves in conformational space. The addition of rotamer moves significantly increases the efficiency of protein structure prediction within this framework, and we determine parameters to optimize efficiency.

Water Networks Repopulate Protein-Ligand Interfaces with Temperature

High-resolution crystal structures highlight the importance of water networks in protein-ligand interactions. However, as these are typically determined at cryogenic temperature, resulting insights may be structurally precise but not biologically accurate. By collecting 10 matched room-temperature and cryogenic datasets of the biomedical target Hsp90α, we identified changes in water networks that impact protein conformations at the ligand binding interface. Water repositioning with temperature repopulates protein ensembles and ligand interactions. We introduce Flipper conformational barcodes to identify temperature-sensitive regions in electron density maps. This revealed that temperature-responsive states coincide with ligand-responsive regions and capture unique binding signatures that disappear upon cryo-cooling. Our results have implications for discovering Hsp90 selective ligands, and, more generally, for the utility of hidden protein and water conformations in drug discovery.

Evaluating the impact of X-ray damage on conformational heterogeneity in room-temperature (277 K) and cryo-cooled protein crystals

Cryo-cooling has been nearly universally adopted to mitigate X-ray damage and facilitate crystal handling in protein X-ray crystallography. However, cryo X-ray crystallographic data provide an incomplete window into the ensemble of conformations that is at the heart of protein function and energetics. Room-temperature (RT) X-ray crystallography provides accurate ensemble information, and recent developments allow conformational heterogeneity (the experimental manifestation of ensembles) to be extracted from single-crystal data. Nevertheless, high sensitivity to X-ray damage at RT raises concerns about data reliability. To systematically address this critical issue, increasingly X-ray-damaged high-resolution data sets (1.02–1.52 Å resolution) were obtained from single proteinase K, thaumatin and lysozyme crystals at RT (277 K). In each case a modest increase in conformational heterogeneity with X-ray damage was observed. Merging data with different extents of damage (as is typically carried out) had negligible effects on conformational heterogeneity until the overall diffraction intensity decayed to ∼70% of its initial value. These effects were compared with X-ray damage effects in cryo-cooled crystals by carrying out an analogous analysis of increasingly damaged proteinase K cryo data sets (0.9–1.16 Å resolution). X-ray damage-associated heterogeneity changes were found that were not observed at RT. This property renders it difficult to distinguish real from artefactual conformations and to determine the conformational response to changes in temperature. The ability to acquire reliable heterogeneity information from single crystals at RT, together with recent advances in RT data collection at accessible synchrotron beamlines, provides a strong motivation for the widespread adoption of RT X-ray crystallography to obtain conformational ensemble information.

Xtrapol8 enables automatic elucidation of low-occupancy intermediate-states in crystallographic studies

Unstable states studied in kinetic, time-resolved and ligand-based crystallography are often characterized by a low occupancy, which hinders structure determination by conventional methods. To automatically extract structural information pertaining to these states, we developed Xtrapol8, a program which (i) applies various flavors of Bayesian-statistics weighting to generate the most informative Fourier difference maps; (ii) determines the occupancy of the intermediate states by use of methods hitherto not available; (iii) calculates extrapolated structure factors using the various proposed formalisms while handling the issue of negative structure factor amplitudes, and (iv) refines the corresponding structures in real and reciprocal-space. The use of Xtrapol8 could accelerate data processing in kinetic and time-resolved crystallographic studies, and as well foster the identification of drug-targetable states in ligand-based crystallography.

The software Xtrapol8 extracts low occupancy states and can thereby accelerate and facilitate data processing in kinetic, timeresolved and ligand-based crystallographic studies.

Ligand binding remodels protein side-chain conformational heterogeneity

While protein conformational heterogeneity plays an important role in many aspects of biological function, including ligand binding, its impact has been difficult to quantify. Macromolecular X-ray diffraction is commonly interpreted with a static structure, but it can provide information on both the anharmonic and harmonic contributions to conformational heterogeneity. Here, through multiconformer modeling of time- and space-averaged electron density, we measure conformational heterogeneity of 743 stringently matched pairs of crystallographic datasets that reflect unbound/apo and ligand-bound/holo states. When comparing the conformational heterogeneity of side chains, we observe that when binding site residues become more rigid upon ligand binding, distant residues tend to become more flexible, especially in non-solvent-exposed regions. Among ligand properties, we observe increased protein flexibility as the number of hydrogen bonds decreases and relative hydrophobicity increases. Across a series of 13 inhibitor-bound structures of CDK2, we find that conformational heterogeneity is correlated with inhibitor features and identify how conformational changes propagate differences in conformational heterogeneity away from the binding site. Collectively, our findings agree with models emerging from nuclear magnetic resonance studies suggesting that residual side-chain entropy can modulate affinity and point to the need to integrate both static conformational changes and conformational heterogeneity in models of ligand binding.

Revealing A-T and G-C Hoogsteen base pairs in stressed protein-bound duplex DNA

Watson–Crick base pairs (bps) are the fundamental unit of genetic information and the building blocks of the DNA double helix. However, A-T and G-C can also form alternative ‘Hoogsteen’ bps, expanding the functional complexity of DNA. We developed ‘Hoog-finder’, which uses structural fingerprints to rapidly screen Hoogsteen bps, which may have been mismodeled as Watson–Crick in crystal structures of protein–DNA complexes. We uncovered 17 Hoogsteen bps, 7 of which were in complex with 6 proteins never before shown to bind Hoogsteen bps. The Hoogsteen bps occur near mismatches, nicks and lesions and some appear to participate in recognition and damage repair. Our results suggest a potentially broad role for Hoogsteen bps in stressed regions of the genome and call for a community-wide effort to identify these bps in current and future crystal structures of DNA and its complexes.

Binding specificity and function of the SWI/SNF subunit SMARCA4 bromodomain interaction with acetylated histone H3K14

Bromodomains (BD) are conserved reader modules that bind acetylated lysine residues on histones. Although much has been learned regarding the in vitro properties of these domains, less is known about their function within chromatin complexes. SWI/SNF chromatin-remodeling complexes modulate transcription and contribute to DNA damage repair. Mutations in SWI/SNF subunits have been implicated in many cancers. Here we demonstrate that the BD of Caenorhabditis elegans SMARCA4/BRG1, a core SWI/SNF subunit, recognizes acetylated lysine 14 of histone H3 (H3K14ac), similar to its Homo sapiens ortholog. We identify the interactions of SMARCA4 with the acetylated histone peptide from a 1.29 Å-resolution crystal structure of the CeSMARCA4 BD-H3K14ac complex. Significantly, most of the SMARCA4 BD residues in contact with the histone peptide are conserved with other proteins containing family VIII bromodomains. Based on the premise that binding specificity is conserved among bromodomain orthologs, we propose that loop residues outside of the binding pocket position contact residues to recognize the H3K14ac sequence. CRISPR-Cas9-mediated mutations in the SMARCA4 BD that abolish H3K14ac binding in vitro had little or no effect on C. elegans viability or physiological function in vivo. However, combining SMARCA4 BD mutations with knockdown of the SWI/SNF accessory subunit PBRM-1 resulted in severe developmental defects in animals. In conclusion, we demonstrated an essential function for the SWI/SNF bromodomain in vivo and detected potential redundancy in epigenetic readers in regulating chromatin remodeling. These findings have implications for the development of small-molecule BD inhibitors to treat cancers and other diseases.

Temperature artifacts in protein structures bias ligand-binding predictions

X-ray crystallography is the gold standard to resolve conformational ensembles that are significant for protein function, ligand discovery, and computational methods development. However, relevant conformational states may be missed at common cryogenic (cryo) data-collection temperatures but can be populated at room temperature. To assess the impact of temperature on making structural and computational discoveries, we systematically investigated protein conformational changes in response to temperature and ligand binding in a structural and computational workhorse, the T4 lysozyme L99A cavity. Despite decades of work on this protein, shifting to RT reveals new global and local structural changes. These include uncovering an apo helix conformation that is hidden at cryo but relevant for ligand binding, and altered side chain and ligand conformations. To evaluate the impact of temperature-induced protein and ligand changes on the utility of structural information in computation, we evaluated how temperature can mislead computational methods that employ cryo structures for validation. We find that when comparing simulated structures just to experimental cryo structures, hidden successes and failures often go unnoticed. When using structural information in ligand binding predictions, both coarse docking and rigorous binding free energy calculations are influenced by temperature effects. The trend that cryo artifacts limit the utility of structures for computation holds across five distinct protein classes. Our results suggest caution when consulting cryogenic structural data alone, as temperature artifacts can conceal errors and prevent successful computational predictions, which can mislead the development and application of computational methods in discovering bioactive molecules.

A Minireview on Temperature Dependent Protein Conformational Sampling

In this minireview we discuss the role of the more subtle conformational change-protein conformational sampling and connect it to the classic relationship of protein structure and function. The theory of pre-existing functional states of protein are discussed in context of alternate protein conformational sampling. Last, we discuss how temperature, ligand binding and mutations affect the protein conformational sampling mode which is linked to the protein function regulation. The review includes several protein systems that showed temperature dependent protein conformational sampling. We also specifically included two enzyme systems, thermophilic alcohol dehydrogenase (ht-ADH) and thermolysin which we previously studied when discussing temperature dependent protein conformational sampling.

Substrate-engaged type III secretion system structures reveal gating mechanism for unfolded protein translocation

Many bacterial pathogens rely on virulent type III secretion systems (T3SSs) or injectisomes to translocate effector proteins in order to establish infection. The central component of the injectisome is the needle complex which assembles a continuous conduit crossing the bacterial envelope and the host cell membrane to mediate effector protein translocation. However, the molecular principles underlying type III secretion remain elusive. Here, we report a structure of an active Salmonella enterica serovar Typhimurium needle complex engaged with the effector protein SptP in two functional states, revealing the complete 800Å-long secretion conduit and unraveling the critical role of the export apparatus (EA) subcomplex in type III secretion. Unfolded substrates enter the EA through a hydrophilic constriction formed by SpaQ proteins, which enables side chain-independent substrate transport. Above, a methionine gasket formed by SpaP proteins functions as a gate that dilates to accommodate substrates while preventing leaky pore formation. Following gate penetration, a moveable SpaR loop first folds up to then support substrate transport. Together, these findings establish the molecular basis for substrate translocation through T3SSs and improve our understanding of bacterial pathogenicity and motility.

Macromolecular room temperature crystallography

X-ray crystallography enables detailed structural studies of proteins to understand and modulate their function. Conducting crystallographic experiments at cryogenic temperatures has practical benefits but potentially limits the identification of functionally important alternative protein conformations that can be revealed only at room temperature (RT). This review discusses practical aspects of preparing, acquiring, and analyzing X-ray crystallography data at RT to demystify preconceived impracticalities that freeze progress of routine RT data collection at synchrotron sources. Examples are presented as conceptual and experimental templates to enable the design of RT-inspired studies; they illustrate the diversity and utility of gaining novel insights into protein conformational landscapes. An integrative view of protein conformational dynamics enables opportunities to advance basic and biomedical research.

qFit 3: Protein and ligand multiconformer modeling for X-ray crystallographic and single-particle cryo-EM density maps

New X-ray crystallography and cryo-electron microscopy (cryo-EM) approaches yield vast amounts of structural data from dynamic proteins and their complexes. Modeling the full conformational ensemble can provide important biological insights, but identifying and modeling an internally consistent set of alternate conformations remains a formidable challenge. qFit efficiently automates this process by generating a parsimonious multiconformer model. We refactored qFit from a distributed application into software that runs efficiently on a small server, desktop, or laptop. We describe the new qFit 3 software and provide some examples. qFit 3 is open-source under the MIT license, and is available at https://github.com/ExcitedStates/qfit-3.0.

Instrumentation and experimental procedures for robust collection of X-ray diffraction data from protein crystals across physiological temperatures

Traditional X-ray diffraction data collected at cryo-temperatures have delivered invaluable insights into the three-dimensional structures of proteins, providing the backbone of structure-function studies. While cryo-cooling mitigates radiation damage, cryo-temperatures can alter protein conformational ensembles and solvent structure. Furthermore, conformational ensembles underlie protein function and energetics, and recent advances in room-temperature X-ray crystallography have delivered conformational heterogeneity information that can be directly related to biological function. Given this capability, the next challenge is to develop a robust and broadly applicable method to collect single-crystal X-ray diffraction data at and above room temperature. This challenge is addressed herein. The approach described provides complete diffraction data sets with total collection times as short as ∼5 s from single protein crystals, dramatically increasing the quantity of data that can be collected within allocated synchrotron beam time. Its applicability was demonstrated by collecting 1.09-1.54 Å resolution data over a temperature range of 293-363 K for proteinase K, thaumatin and lysozyme crystals at BL14-1 at the Stanford Synchrotron Radiation Lightsource. The analyses presented here indicate that the diffraction data are of high quality and do not suffer from excessive dehydration or radiation damage.

Discovery of allosteric binding sites by crystallographic fragment screening

Understanding allosteric regulation of proteins is fundamental to our study of protein structure and function. Moreover, allosteric binding pockets have become a major target of drug discovery efforts in recent years. However, even though the function of almost every protein can be influenced by allostery, it remains a challenge to discover, rationalise and validate putative allosteric binding pockets. This review examines how the discovery and analysis of putative allosteric binding sites have been influenced by the availability of centralised facilities for crystallographic fragment screening, along with newly developed computational methods for modelling low occupancy features. We discuss the experimental parameters required for success, and how new methods could influence the field in the future. Finally, we reflect on the general problem of how to translate these findings into actual ligand development programs.

Conformational heterogeneity of Savinase from NMR, HDX-MS and X-ray diffraction analysis

Background

Several examples have emerged of enzymes where slow conformational changes are of key importance for function and where low populated conformations in the resting enzyme resemble the conformations of intermediate states in the catalytic process. Previous work on the subtilisin protease, Savinase, from Bacillus lentus by NMR spectroscopy suggested that this enzyme undergoes slow conformational dynamics around the substrate binding site. However, the functional importance of such dynamics is unknown.

Methods

Here we have probed the conformational heterogeneity in Savinase by following the temperature dependent chemical shift changes. In addition, we have measured changes in the local stability of the enzyme when the inhibitor phenylmethylsulfonyl fluoride is bound using hydrogen-deuterium exchange mass spectrometry (HDX-MS). Finally, we have used X-ray crystallography to compare electron densities collected at cryogenic and ambient temperatures and searched for possible low populated alternative conformations in the crystals.

Results

The NMR temperature titration shows that Savinase is most flexible around the active site, but no distinct alternative states could be identified. The HDX shows that modification of Savinase with inhibitor has very little impact on the stability of hydrogen bonds and solvent accessibility of the backbone. The most pronounced structural heterogeneities detected in the diffraction data are limited to alternative side-chain rotamers and a short peptide segment that has an alternative main-chain conformation in the crystal at cryo conditions. Collectively, our data show that there is very little structural heterogeneity in the resting state of Savinase and hence that Savinase does not rely on conformational selection to drive the catalytic process.

Advances in methods for atomic resolution macromolecular structure determination

Recent technical advances have dramatically increased the power and scope of structural biology. New developments in high-resolution cryo-electron microscopy, serial X-ray crystallography, and electron diffraction have been especially transformative. Here we highlight some of the latest advances and current challenges at the frontiers of atomic resolution methods for elucidating the structures and dynamical properties of macromolecules and their complexes.

Ectopic suicide inhibition of thioredoxin glutathione reductase

Selective suicide inhibitors represent a seductively attractive approach for inactivation of therapeutically relevant enzymes since they are generally devoid of off-target toxicity in vivo. While most suicide inhibitors are converted to reactive species at enzyme active sites, theoretically bioactivation can also occur in ectopic (secondary) sites that have no known function. Here, we report an example of such an "ectopic suicide inhibition", an unprecedented bioactivation mechanism of a suicide inhibitor carried out by a non-catalytic site of thioredoxin glutathione reductase (TGR). TGR is a promising drug target to treat schistosomiasis, a devastating human parasitic disease. Utilizing hits selected from a high throughput screening campaign, time-resolved X-ray crystallography, molecular dynamics, mass spectrometry, molecular modeling, protein mutagenesis and functional studies, we find that 2-naphtholmethylamino derivatives bound to this novel ectopic site of Schistosoma mansoni (Sm)TGR are transformed to covalent modifiers and react with its mobile selenocysteine-containing C-terminal arm. In particular, one 2-naphtholmethylamino compound is able to specifically induce the pro-oxidant activity in the inhibited enzyme. Since some 2-naphtholmethylamino analogues show worm killing activity and the ectopic site is not conserved in human orthologues, a general approach to development of novel and selective anti-parasitic therapeutics against schistosoma is proposed.

Conformational sampling and kinetics changes across a non-Arrhenius break point in the enzyme thermolysin

Numerous studies have suggested a significant role that protein dynamics play in optimizing enzyme catalysis, and changes in conformational sampling offer a window to explore this role. Thermolysin from Bacillus thermoproteolyticus rokko, which is a heat-stable zinc metalloproteinase, serves here as a model system to study changes of protein function and conformational sampling across a temperature range of 16-36 °C. The temperature dependence of kinetics of thermolysin showed a biphasic transition at 26 °C that points to potential conformational and dynamic differences across this temperature. The non-Arrhenius behavior observed resembled results from previous studies of a thermophilic alcohol dehydrogenase enzyme, which also indicated a biphasic transition at ambient temperatures. To explore the non-Arrhenius behavior of thermolysin, room temperature crystallography was applied to characterize structural changes in a temperature range across the biphasic transition temperature. The alternate conformation of side chain fitting to electron density of a group of residues showed a higher variability in the temperature range from 26 to 29 °C, which indicated a change in conformational sampling that correlated with the non-Arrhenius break point.

Computational Analysis of Energy Landscapes Reveals Dynamic Features That Contribute to Binding of Inhibitors to CFTR-Associated Ligand

The CFTR-associated ligand PDZ domain (CALP) binds to the cystic fibrosis transmembrane conductance regulator (CFTR) and mediates lysosomal degradation of mature CFTR. Inhibition of this interaction has been explored as a therapeutic avenue for cystic fibrosis. Previously, we reported the ensemble-based computational design of a novel peptide inhibitor of CALP, which resulted in the most binding-efficient inhibitor to date. This inhibitor, kCAL01, was designed using osprey and evinced significant biological activity in in vitro cell-based assays. Here, we report a crystal structure of kCAL01 bound to CALP and compare structural features against iCAL36, a previously developed inhibitor of CALP. We compute side-chain energy landscapes for each structure to not only enable approximation of binding thermodynamics but also reveal ensemble features that contribute to the comparatively efficient binding of kCAL01. Finally, we compare the previously reported design ensemble for kCAL01 vs the new crystal structure and show that, despite small differences between the design model and crystal structure, significant biophysical features that enhance inhibitor binding are captured in the design ensemble. This suggests not only that ensemble-based design captured thermodynamically significant features observed in vitro, but also that a design eschewing ensembles would miss the kCAL01 sequence entirely.

Solvent flows, conformation changes and lattice reordering in a cold protein crystal

When protein crystals are abruptly cooled, the unit-cell, protein and solvent-cavity volumes all contract, but the volume of bulk-like internal solvent may expand. Outflow of this solvent from the unit cell and its accumulation in defective interior crystal regions has been suggested as one cause of the large increase in crystal mosaicity on cooling. It is shown that when apoferritin crystals are abruptly cooled to temperatures between 220 and 260 K, the unit cell contracts, solvent is pushed out and the mosaicity grows. On temperature-dependent timescales of 10 to 200 s, the unit-cell and solvent-cavity volume then expand, solvent flows back in, and the mosaicity and B factor both drop. Expansion and reordering at fixed low temperature are associated with small-amplitude but large-scale changes in the conformation and packing of apoferritin. These results demonstrate that increases in mosaicity on cooling arise due to solvent flows out of or into the unit cell and to incomplete, arrested relaxation of protein conformation. They indicate a critical role for time in variable-temperature crystallographic studies, and the feasibility of probing interactions and cooperative conformational changes that underlie cold denaturation in the presence of liquid solvent at temperatures down to ∼200 K.

A comparison of gas stream cooling and plunge cooling of macromolecular crystals

Cryocooling for macromolecular crystallography is usually performed via plunging the crystal into a liquid cryogen or placing the crystal in a cold gas stream. These two approaches are compared here for the case of nitro-gen cooling. The results show that gas stream cooling, which typically cools the crystal more slowly, yields lower mosaicity and, in some cases, a stronger anomalous signal relative to rapid plunge cooling. During plunging, moving the crystal slowly through the cold gas layer above the liquid surface can produce mosaicity similar to gas stream cooling. Annealing plunge cooled crystals by warming and recooling in the gas stream allows the mosaicity and anomalous signal to recover. For tetragonal thermolysin, the observed effects are less pronounced when the cryosolvent has smaller thermal contraction, under which conditions the protein structures from plunge cooled and gas stream cooled crystals are very similar. Finally, this work also demonstrates that the resolution dependence of the reflecting range is correlated with the cooling method, suggesting it may be a useful tool for discerning whether crystals are cooled too rapidly. The results support previous studies suggesting that slower cooling methods are less deleterious to crystal order, as long as ice formation is prevented and dehydration is limited.

ACMS: a database of alternate conformations found in the atoms of main and side chains of protein structures

Proteins are usually dynamic biological macromolecules, thereby exhibiting a large number of conformational ensembles which influence the association with their neighbours and interacting partners. Most of the side-chain atoms and a few main-chain atoms of the high-resolution crystal structures deposited in the Protein Data Bank adopt alternate conformations. This kind of conformational behaviour prompted the authors to explore the relationship, if any, between the alternate conformations and the function of the protein molecule. Thus, a knowledge base of the alternate conformations of the main- and side-chain atoms of protein structures has been developed. It provides a detailed description of the alternate conformations of various residues for more than 60 000 high-resolution crystal structures. The proposed knowledge base is very user friendly and has various flexible options. The knowledge base will be updated periodically and can be accessed at http://iris.physics.iisc.ac.in/acms.

Statistical Framework for Uncertainty Quantification in Computational Molecular Modeling

As computational modeling, simulation, and predictions are becoming integral parts of biomedical pipelines, it behooves us to emphasize the reliability of the computational protocol. For any reported quantity of interest (QOI), one must also compute and report a measure of the uncertainty or error associated with the QOI. This is especially important in molecular modeling, since in most practical applications the inputs to the computational protocol are often noisy, incomplete, or low-resolution. Unfortunately, currently available modeling tools do not account for uncertainties and their effect on the final QOIs with sufficient rigor. We have developed a statistical framework that expresses the uncertainty of the QOI as the probability that the reported value deviates from the true value by more than some user-defined threshold. First, we provide a theoretical approach where this probability can be bounded using Azuma-Hoeffding like inequalities. Second, we approximate this probability empirically by sampling the space of uncertainties of the input and provide applications of our framework to bound uncertainties of several QOIs commonly used in molecular modeling. Finally, we also present several visualization techniques to effectively and quantitavely visualize the uncertainties: in the input, final QOIs, and also intermediate states.

Dynamic ion pair behavior stabilizes single α-helices in proteins

Ion pairs are key stabilizing interactions between oppositely charged amino acid side chains in proteins. They are often depicted as single conformer salt bridges (hydrogen-bonded ion pairs) in crystal structures, but it is unclear how dynamic they are in solution. Ion pairs are thought to be particularly important in stabilizing single α-helix (SAH) domains in solution. These highly stable domains are rich in charged residues (such as Arg, Lys, and Glu) with potential ion pairs across adjacent turns of the helix. They provide a good model system to investigate how ion pairs can contribute to protein stability. Using NMR spectroscopy, small-angle X-ray light scattering (SAXS), and molecular dynamics simulations, we provide here experimental evidence that ion pairs exist in a SAH in murine myosin 7a (residues 858-935), but that they are not fixed or long lasting. In silico modeling revealed that the ion pairs within this α-helix exhibit dynamic behavior, rapidly forming and breaking and alternating between different partner residues. The low-energy helical state was compatible with a great variety of ion pair combinations. Flexible ion pair formation utilizing a subset of those available at any one time avoided the entropic penalty of fixing side chain conformations, which likely contributed to helix stability overall. These results indicate the dynamic nature of ion pairs in SAHs. More broadly, thermodynamic stability in other proteins is likely to benefit from the dynamic behavior of multi-option solvent-exposed ion pairs.

Journey to the center of the protein: allostery from multitemperature multiconformer X-ray crystallography

Proteins inherently fluctuate between conformations to perform functions in the cell. For example, they sample product-binding, transition-state-stabilizing and product-release states during catalysis, and they integrate signals from remote regions of the structure for allosteric regulation. However, there is a lack of understanding of how these dynamic processes occur at the basic atomic level. This gap can be at least partially addressed by combining variable-temperature (instead of traditional cryogenic temperature) X-ray crystallography with algorithms for modeling alternative conformations based on electron-density maps, in an approach called multitemperature multiconformer X-ray crystallography (MMX). Here, the use of MMX to reveal alternative conformations at different sites in a protein structure and to estimate the degree of energetic coupling between them is discussed. These insights can suggest testable hypotheses about allosteric mechanisms. Temperature is an easily manipulated experimental parameter, so the MMX approach is widely applicable to any protein that yields well diffracting crystals. Moreover, the general principles of MMX are extensible to other perturbations such as pH, pressure, ligand concentration etc. Future work will explore strategies for leveraging X-ray data across such perturbation series to more quantitatively measure how different parts of a protein structure are coupled to each other, and the consequences thereof for allostery and other aspects of protein function.

qFit-ligand Reveals Widespread Conformational Heterogeneity of Drug-Like Molecules in X-Ray Electron Density Maps

Proteins and ligands sample a conformational ensemble that governs molecular recognition, activity, and dissociation. In structure-based drug design, access to this conformational ensemble is critical to understand the balance between entropy and enthalpy in lead optimization. However, ligand conformational heterogeneity is currently severely underreported in crystal structures in the Protein Data Bank, owing in part to a lack of automated and unbiased procedures to model an ensemble of protein-ligand states into X-ray data. Here, we designed a computational method, qFit-ligand, to automatically resolve conformationally averaged ligand heterogeneity in crystal structures, and applied it to a large set of protein receptor-ligand complexes. In an analysis of the cancer related BRD4 domain, we found that up to 29% of protein crystal structures bound with drug-like molecules present evidence of unmodeled, averaged, relatively isoenergetic conformations in ligand-receptor interactions. In many retrospective cases, these alternate conformations were adventitiously exploited to guide compound design, resulting in improved potency or selectivity. Combining qFit-ligand with high-throughput screening or multitemperature crystallography could therefore augment the structure-based drug design toolbox.

A Multi-model Approach to Assessing Local and Global Cryo-EM Map Quality

There does not currently exist a standardized indicator of how well cryo-EM-derived models represent the density from which they were generated. We present a straightforward methodology that utilizes freely available tools to generate a suite of independent models and to evaluate their convergence in an EM density. These analyses provide both a quantitative and qualitative assessment of the precision of the models and their representation of the density, respectively, while concurrently providing a platform for assessing both global and local EM map quality. We further use standardized datasets to provide an expected deviation within a suite of models refined against EM maps reported to be at 5 Å resolution or better. Associating multiple atomic models with a deposited EM map provides a rapid and accessible reporter of convergence, a strong indicator of highly resolved molecular detail, and is an important step toward an FSC-independent assessment of map and model quality.