Distributed computing for macromolecular crystallography

Emerging Topics in Protein Crystallography

Protein crystallography is the discipline concerned with the determination of the three-dimensional structure of biological macromolecules in a crystalline state [...].

Diverse crystalline protein scaffolds through metal-dependent polymorphism

As protein crystals are increasingly finding diverse applications as scaffolds, controlled crystal polymorphism presents a facile strategy to form crystalline assemblies with controllable porosity with minimal to no protein engineering. Polymorphs of consensus tetratricopeptide repeat proteins with varying porosity were obtained through co-crystallization with metal salts, exploiting the innate metal ion geometric requirements. A single structurally exposed negative amino acid cluster was responsible for metal coordination, despite the abundance of negatively charged residues. Density functional theory calculations showed that while most of the crystals were the most thermodynamically stable assemblies, some were kinetically trapped states. Thus, crystalline porosity diversity is achieved and controlled with metal coordination, opening a new scope in the application of proteins as biocompatible protein-metal-organic frameworks (POFs). In addition, metal-dependent polymorphic crystals allow direct comparison of metal coordination preferences.

Structural and biochemical characterisation of the N-carbamoyl-β-alanine amidohydrolase from Rhizobium radiobacter MDC 8606

N-carbamoyl-β-alanine amidohydrolase (CβAA) constitutes one of the most important groups of industrially relevant enzymes used in the production of optically pure amino acids and derivatives. In this study, a CβAA-encoding gene from Rhizobium radiobacter strain MDC 8606 was cloned and overexpressed in Escherichia coli. The purified recombinant enzyme (RrCβAA) showed a specific activity of 14 U·mg-1 using N-carbamoyl-β-alanine as a substrate with an optimum activity at 55 °C and pH 8.0. In this work, we report also the first prokaryotic CβAA structure at a resolution of 2.0 Å. A discontinuous catalytic domain and a dimerisation domain attached through a flexible hinge region at the domain interface have been revealed. We identify key ligand binding residues, including a conserved glutamic acid (Glu131), histidine (H385) and arginine (Arg291). Our results allowed us to explain the preference of the enzyme for linear carbamoyl substrates, as large and branched carbamoyl substrates cannot fit in the active site of the enzyme. This work envisages the use of RrCβAA from R. radiobacter MDC 8606 for the industrial production of L-α-, L-β- and L-γ-amino acids. The structural analysis provides new insights on enzyme-substrate interaction, which shed light on engineering of CβAAs for high catalytic activity and broad substrate specificity.

Structure-Based Lead Optimization of Enterovirus D68 2A Protease Inhibitors

Enterovirus D68 (EV-D68) virus is a nonpolio enterovirus that typically causes respiratory illness and, in severe cases, can lead to paralysis and death in children. There is currently no vaccine or antiviral for EV-D68. We previously discovered the viral 2A protease (2Apro) as a viable antiviral drug target and identified telaprevir as a 2Apro inhibitor. 2Apro is a viral cysteine protease that cleaves the viral VP1-2A polyprotein junction. In this study, we report the X-ray crystal structures of EV-D68 2Apro, wild-type, and the C107A mutant and the structure-based lead optimization of telaprevir. Guided by the X-ray crystal structure, we predicted the binding pose of telaprevir in 2Apro using molecular dynamics simulations. We then utilized this model to inform structure-based optimization of the telaprevir's reactive warhead and P1-P4 substitutions. These efforts led to the discovery of 2Apro inhibitors with improved antiviral activity than telaprevir. These compounds represent promising lead compounds for further development as EV-D68 antivirals.

XDSGUI: a graphical user interface for XDS, SHELX and ARCIMBOLDO

XDSGUI is a lightweight graphical user interface (GUI) for the XDS, SHELX and ARCIMBOLDO program packages that serves both novice and experienced users in obtaining optimal processing and phasing results for X-ray, neutron and electron diffraction data. The design of the program enables data processing and phasing without command line usage, and supports advanced command flows in a simple user-modifiable and user-extensible way. The GUI supplies graphical information based on the tabular log output of the programs, which is more intuitive, comprehensible and efficient than text output can be.

Predicted models and CCP4

In late 2020, the results of CASP14, the 14th event in a series of competitions to assess the latest developments in computational protein structure-prediction methodology, revealed the giant leap forward that had been made by Google's Deepmind in tackling the prediction problem. The level of accuracy in their predictions was the first instance of a competitor achieving a global distance test score of better than 90 across all categories of difficulty. This achievement represents both a challenge and an opportunity for the field of experimental structural biology. For structure determination by macromolecular X-ray crystallography, access to highly accurate structure predictions is of great benefit, particularly when it comes to solving the phase problem. Here, details of new utilities and enhanced applications in the CCP4 suite, designed to allow users to exploit predicted models in determining macromolecular structures from X-ray diffraction data, are presented. The focus is mainly on applications that can be used to solve the phase problem through molecular replacement.

BtuB TonB-dependent transporters and BtuG surface lipoproteins form stable complexes for vitamin B12 uptake in gut Bacteroides

Vitamin B12 (cobalamin) is required for most human gut microbes, many of which are dependent on scavenging to obtain this vitamin. Since bacterial densities in the gut are extremely high, competition for this keystone micronutrient is severe. Contrasting with Enterobacteria, members of the dominant genus Bacteroides often encode several BtuB vitamin B12 outer membrane transporters together with a conserved array of surface-exposed B12-binding lipoproteins. Here we show that the BtuB transporters from Bacteroides thetaiotaomicron form stable, pedal bin-like complexes with surface-exposed BtuG lipoprotein lids, which bind B12 with high affinities. Closing of the BtuG lid following B12 capture causes destabilisation of the bound B12 by a conserved BtuB extracellular loop, causing translocation of the vitamin to BtuB and subsequent transport. We propose that TonB-dependent, lipoprotein-assisted small molecule uptake is a general feature of Bacteroides spp. that is important for the success of this genus in colonising the human gut.

Outer membrane utilisomes mediate glycan uptake in gut Bacteroidetes

Bacteroidetes are abundant members of the human microbiota, utilizing a myriad of diet- and host-derived glycans in the distal gut¹. Glycan uptake across the bacterial outer membrane of these bacteria is mediated by SusCD protein complexes, comprising a membrane-embedded barrel and a lipoprotein lid, which is thought to open and close to facilitate substrate binding and transport. However, surface-exposed glycan-binding proteins and glycoside hydrolases also play critical roles in the capture, processing and transport of large glycan chains. The interactions between these components in the outer membrane are poorly understood, despite being crucial for nutrient acquisition by our colonic microbiota. Here we show that for both the levan and dextran utilization systems of Bacteroides thetaiotaomicron, the additional outer membrane components assemble on the core SusCD transporter, forming stable glycan-utilizing machines that we term utilisomes. Single-particle cryogenic electron microscopy structures in the absence and presence of substrate reveal concerted conformational changes that demonstrate the mechanism of substrate capture, and rationalize the role of each component in the utilisome.

Important Structural Features of Thiolate-Rich Four-Helix Bundles for Cu(I) Uptake and Removal

A family of bacterial copper storage proteins (the Csps) possess thiolate-lined four-helix bundles whose cores can be filled with Cu(I) ions. The majority of Csps are cytosolic (Csp3s), and in vitro studies carried out to date indicate that the Csp3s from Methylosinus trichosporium OB3b (MtCsp3), Bacillus subtilis (BsCsp3), and Streptomyces lividans (SlCsp3) are alike. Bioinformatics have highlighted homologues with potentially different Cu(I)-binding properties from these characterized "classical" Csp3s. Determination herein of the crystal structure of the protein (RkCsp3) from the methanotroph Methylocystis sp. strain Rockwell with Cu(I) bound identifies this as the first studied example of a new subgroup of Csp3s. The most significant structural difference from classical Csp3s is the presence of only two Cu(I) sites at the mouth of the bundle via which Cu(I) ions enter and leave. This is due to the absence of three Cys residues and a His-containing motif, which allow classical Csp3s to bind five to six Cu(I) ions in this region. Regardless, RkCsp3 exhibits rapid Cu(I) binding and the fastest measured Cu(I) removal rate for a Csp3 when using high-affinity ligands as surrogate partners. New experiments on classical Csp3s demonstrate that their His-containing motif is not essential for fast Cu(I) uptake and removal. Other structural features that could be important for these functionally relevant in vitro properties are discussed.

An ancient metalloenzyme evolves through metal preference modulation

Evolution creates functional diversity of proteins, the essential building blocks of all biological systems. However, studies of natural proteins sampled across the tree of life and evaluated in a single experimental system are lacking. Almost half of enzymes require metals, and metalloproteins tend to optimally utilize the physicochemical properties of a specific metal co-factor. Life must adapt to changes in metal bioavailability, including those during the transition from anoxic to oxic Earth or pathogens' exposure to nutritional immunity. These changes can challenge the ability of metalloenzymes to maintain activity, presumptively driving their evolution. Here we studied metal-preference evolution within the natural diversity of the iron/manganese superoxide dismutase (SodFM) family of reactive oxygen species scavengers. We identified and experimentally verified residues with conserved roles in determining metal preference that, when combined with an understanding of the protein's evolutionary history, improved prediction of metal utilization across the five SodFM subfamilies defined herein. By combining phylogenetics, biochemistry and structural biology, we demonstrate that SodFM metal utilization can be evolutionarily fine tuned by sliding along a scale between perfect manganese and iron specificities. Over the history of life, SodFM metal preference has been modulated multiple independent times within different evolutionary and ecological contexts, and can be changed within short evolutionary timeframes.

Structure and Dynamics of Three Escherichia coli NfsB Nitro-Reductase Mutants Selected for Enhanced Activity with the Cancer Prodrug CB1954

Escherichia coli NfsB has been studied extensively for its potential for cancer gene therapy by reducing the prodrug CB1954 to a cytotoxic derivative. We have previously made several mutants with enhanced activity for the prodrug and characterised their activity in vitro and in vivo. Here, we determine the X-ray structure of our most active triple and double mutants to date, T41Q/N71S/F124T and T41L/N71S. The two mutant proteins have lower redox potentials than wild-type NfsB, and the mutations have lowered activity with NADH so that, in contrast to the wild-type enzyme, the reduction of the enzyme by NADH, rather than the reaction with CB1954, has a slower maximum rate. The structure of the triple mutant shows the interaction between Q41 and T124, explaining the synergy between these two mutations. Based on these structures, we selected mutants with even higher activity. The most active one contains T41Q/N71S/F124T/M127V, in which the additional M127V mutation enlarges a small channel to the active site. Molecular dynamics simulations show that the mutations or reduction of the FMN cofactors of the protein has little effect on its dynamics and that the largest backbone fluctuations occur at residues that flank the active site, contributing towards its broad substrate range.

Structural and Biochemical Characterization of Silver/Copper Binding by Dendrorhynchus zhejiangensis Ferritin

Ferritin with a highly symmetrical cage-like structure is not only key in the reversible storage of iron in efficient ferroxidase activity; it also provides unique coordination environments for the conjugation of heavy metal ions other than those associated with iron. However, research regarding the effect of these bound heavy metal ions on ferritin is scarce. In the present study, we prepared a marine invertebrate ferritin from Dendrorhynchus zhejiangensis (DzFer) and found that it could withstand extreme pH fluctuation. We then demonstrated its capacity to interact with Ag⁺ or Cu²⁺ ions using various biochemical and spectroscopic methods and X-ray crystallography. Structural and biochemical analyses revealed that both Ag⁺ and Cu²⁺ were able to bind to the DzFer cage via metal-coordination bonds and that their binding sites were mainly located inside the three-fold channel of DzFer. Furthermore, Ag⁺ was shown to have a higher selectivity for sulfur-containing amino acid residues and appeared to bind preferentially at the ferroxidase site of DzFer as compared with Cu²⁺. Thus, it is far more likely to inhibit the ferroxidase activity of DzFer. The results provide new insights into the effect of heavy metal ions on the iron-binding capacity of a marine invertebrate ferritin.

A Cold-Active Flavin-Dependent Monooxygenase from Janthinobacterium svalbardensis Unlocks Applications of Baeyer–Villiger Monooxygenases at Low Temperature

Cold-active enzymes maintain a large part of their optimal activity at low temperatures. Therefore, they can be used to avoid side reactions and preserve heat-sensitive compounds. Baeyer-Villiger monooxygenases (BVMO) utilize molecular oxygen as a co-substrate to catalyze reactions widely employed for steroid, agrochemical, antibiotic, and pheromone production. Oxygen has been described as the rate-limiting factor for some BVMO applications, thereby hindering their efficient utilization. Considering that oxygen solubility in water increases by 40% when the temperature is decreased from 30 to 10 °C, we set out to identify and characterize a cold-active BVMO. Using genome mining in the Antarctic organism Janthinobacterium svalbardensis, a cold-active type II flavin-dependent monooxygenase (FMO) was discovered. The enzyme shows promiscuity toward NADH and NADPH and high activity between 5 and 25 °C. The enzyme catalyzes the monooxygenation and sulfoxidation of a wide range of ketones and thioesters. The high enantioselectivity in the oxidation of norcamphor (eeS = 56%, eeP > 99%, E > 200) demonstrates that the generally higher flexibility observed in the active sites of cold-active enzymes, which compensates for the lower motion at cold temperatures, does not necessarily reduce the selectivity of these enzymes. To gain a better understanding of the unique mechanistic features of type II FMOs, we determined the structure of the dimeric enzyme at 2.5 Å resolution. While the unusual N-terminal domain has been related to the catalytic properties of type II FMOs, the structure shows a SnoaL-like N-terminal domain that is not interacting directly with the active site. The active site of the enzyme is accessible only through a tunnel, with Tyr-458, Asp-217, and His-216 as catalytic residues, a combination not observed before in FMOs and BVMOs.

Structural insights into 3Fe-4S ferredoxins diversity in M. tuberculosis highlighted by a first redox complex with P450

Ferredoxins are small iron-sulfur proteins and key players in essential metabolic pathways. Among all types, 3Fe-4S ferredoxins are less studied mostly due to anaerobic requirements. Their complexes with cytochrome P450 redox partners have not been structurally characterized. In the present work, we solved the structures of both 3Fe-4S ferredoxins from M. tuberculosis-Fdx alone and the fusion FdxE-CYP143. Our SPR analysis demonstrated a high-affinity binding of FdxE to CYP143. According to SAXS data, the same complex is present in solution. The structure reveals extended multipoint interactions and the shape/charge complementarity of redox partners. Furthermore, FdxE binding induced conformational changes in CYP143 as evident from the solved CYP143 structure alone. The comparison of FdxE-CYP143 and modeled Fdx-CYP51 complexes further revealed the specificity of ferredoxins. Our results illuminate the diversity of electron transfer complexes for the production of different secondary metabolites.

The Crystal Structure of Mouse Ces2c, a Potential Ortholog of Human CES2, Shows Structural Similarities in Substrate Regulation and Product Release to Human CES1

Members of the carboxylesterase 2 (Ces2/CES2) family have been studied intensively with respect to their hydrolytic function on (pro)drugs, whereas their physiological role in lipid and energy metabolism has been realized only within the last few years. Humans have one CES2 gene which is highly expressed in liver, intestine, and kidney. Interestingly, eight homologous Ces2 (Ces2a to Ces2h) genes exist in mice and the individual roles of the corresponding proteins are incompletely understood. Mouse Ces2c (mCes2c) is suggested as potential ortholog of human CES2. Therefore, we aimed at its structural and biophysical characterization. Here, we present the first crystal structure of mCes2c to 2.12 Å resolution. The overall structure of mCes2c resembles that of the human CES1 (hCES1). The core domain adopts an α/β hydrolase-fold with S230, E347, and H459 forming a catalytic triad. Access to the active site is restricted by the cap, the flexible lid, and the regulatory domain. The conserved gate (M417) and switch (F418) residues might have a function in product release similar as suggested for hCES1. Biophysical characterization confirms that mCes2c is a monomer in solution. Thus, this study broadens our understanding of the mammalian carboxylesterase family and assists in delineating the similarities and differences of the different family members.

The 3D-structure, kinetics and dynamics of the E. coli nitroreductase NfsA with NADP+ provide glimpses of its catalytic mechanism

Nitroreductases activate nitroaromatic antibiotics and cancer prodrugs to cytotoxic hydroxylamines and reduce quinones to quinols. Using steady-state and stopped-flow kinetics, we show that the Escherichia coli nitroreductase NfsA is 20-50 fold more active with NADPH than with NADH and that product release may be rate-limiting. The crystal structure of NfsA with NADP⁺ shows that a mobile loop forms a phosphate-binding pocket. The nicotinamide ring and nicotinamide ribose are mobile, as confirmed in molecular dynamics (MD) simulations. We present a model of NADPH bound to NfsA. Only one NADP⁺ is seen bound to the NfsA dimers, and MD simulations show that binding of a second NADP(H) cofactor is unfavourable, suggesting that NfsA and other members of this protein superfamily may have a half-of-sites mechanism.

ModelCraft: an advanced automated model-building pipeline using Buccaneer

Interactive model building can be a difficult and time-consuming step in the structure-solution process. Automated model-building programs such as Buccaneer often make it quicker and easier by completing most of the model in advance. However, they may fail to do so with low-resolution data or a poor initial model or map. The Buccaneer pipeline is a relatively simple program that iterates Buccaneer with REFMAC to refine the model and update the map. A new pipeline called ModelCraft has been developed that expands on this to include shift-field refinement, machine-learned pruning of incorrect residues, classical density modification, addition of water and dummy atoms, building of nucleic acids and final rebuilding of side chains. Testing was performed on 1180 structures solved by experimental phasing, 1338 structures solved by molecular replacement using homologues and 2030 structures solved by molecular replacement using predicted AlphaFold models. Compared with the previous Buccaneer pipeline, ModelCraft increased the mean completeness of the protein models in the experimental phasing cases from 91% to 95%, the molecular-replacement cases from 50% to 78% and the AlphaFold cases from 82% to 91%.

CCP4 Cloud for structure determination and project management in macromolecular crystallography

Nowadays, progress in the determination of three-dimensional macromolecular structures from diffraction images is achieved partly at the cost of increasing data volumes. This is due to the deployment of modern high-speed, high-resolution detectors, the increased complexity and variety of crystallographic software, the use of extensive databases and high-performance computing. This limits what can be accomplished with personal, offline, computing equipment in terms of both productivity and maintainability. There is also an issue of long-term data maintenance and availability of structure-solution projects as the links between experimental observations and the final results deposited in the PDB. In this article, CCP4 Cloud, a new front-end of the CCP4 software suite, is presented which mitigates these effects by providing an online, cloud-based environment for crystallographic computation. CCP4 Cloud was developed for the efficient delivery of computing power, database services and seamless integration with web resources. It provides a rich graphical user interface that allows project sharing and long-term storage for structure-solution projects, and can be linked to data-producing facilities. The system is distributed with the CCP4 software suite version 7.1 and higher, and an online publicly available instance of CCP4 Cloud is provided by CCP4.

De novo determination of mosquitocidal Cry11Aa and Cry11Ba structures from naturally-occurring nanocrystals

Cry11Aa and Cry11Ba are the two most potent toxins produced by mosquitocidal Bacillus thuringiensis subsp. israelensis and jegathesan, respectively. The toxins naturally crystallize within the host; however, the crystals are too small for structure determination at synchrotron sources. Therefore, we applied serial femtosecond crystallography at X-ray free electron lasers to in vivo-grown nanocrystals of these toxins. The structure of Cry11Aa was determined de novo using the single-wavelength anomalous dispersion method, which in turn enabled the determination of the Cry11Ba structure by molecular replacement. The two structures reveal a new pattern for in vivo crystallization of Cry toxins, whereby each of their three domains packs with a symmetrically identical domain, and a cleavable crystal packing motif is located within the protoxin rather than at the termini. The diversity of in vivo crystallization patterns suggests explanations for their varied levels of toxicity and rational approaches to improve these toxins for mosquito control.

The integration of AlphaFold-predicted and crystal structures of human trans-3-hydroxy-l-proline dehydratase reveals a regulatory catalytic mechanism

Computational methods for protein structure prediction have made significant strides forward, as evidenced by the last development of the neural network AlphaFold, which outperformed the CASP14 competitors by consistently predicting the structure of target proteins. Here we show an integrated structural investigation that combines the AlphaFold and crystal structures of human trans-3-Hydroxy-l-proline dehydratase, an enzyme involved in hydroxyproline catabolism and whose structure had never been reported before, identifying a structural element, absent in the AlphaFold model but present in the crystal structure, that was subsequently proved to be functionally relevant. Although the AlphaFold model lacked information on protein oligomerization, the native dimer was reconstructed using template-based and ab initio computational approaches. Moreover, molecular phasing of the diffraction data using the AlphaFold model resulted in dimer reconstruction and straightforward structure solution. Our work adds to the integration of AlphaFold with experimental structural and functional data for protein analysis, crystallographic phasing and structure solution.

Identification, structure determination and analysis of Mycobacterium smegmatis acyl-carrier protein synthase (AcpS) crystallized serendipitously

The unintended crystallization of proteins which generally originate from the expression host instead of the target recombinant proteins is periodically reported. Despite the massive technological advances in the field, assigning a structural model to the corresponding diffraction data is not a trivial task. Here, the structure of acyl-carrier protein synthase (AcpS) from Mycobacterium smegmatis (msAcpS), which crystallized inadvertently in an experimental setup to grow crystals of a Mycobacterium tuberculosis protein using M. smegmatis as an expression system, is reported. After numerous unsuccessful attempts to solve the structure of the target protein by the molecular-replacement method no convincing solutions were obtained, indicating that the diffraction data may correspond to a crystal of an artifactual protein, which was finally identified by the Sequence-Independent Molecular replacement Based on Available Databases (SIMBAD) server. The msAcpS structure was solved at 2.27 Å resolution and structural analysis showed an overall conserved fold. msAcpS formed a trimeric structure similar to those of other reported structures of AcpS from various organisms; however, the residues involved in trimer formation are not strictly conserved. An unrelated metal ion (Ni2+), which was possibly incorporated during protein purification, was observed in the proximity of His49 and His116. Structural and sequence differences were observed in the loop connecting the α3 and α4 helices that is responsible for the open and closed conformations of the enzyme. Moreover, the structural analysis of msAcpS augments the current understanding of this enzyme, which plays a crucial role in the functional activation of acyl-carrier proteins in the fatty-acid biosynthesis pathway.

Characterization of inositol lipid metabolism in gut-associated Bacteroidetes

Inositol lipids are ubiquitous in eukaryotes and have finely tuned roles in cellular signalling and membrane homoeostasis. In Bacteria, however, inositol lipid production is relatively rare. Recently, the prominent human gut bacterium Bacteroides thetaiotaomicron (BT) was reported to produce inositol lipids and sphingolipids, but the pathways remain ambiguous and their prevalence unclear. Here, using genomic and biochemical approaches, we investigated the gene cluster for inositol lipid synthesis in BT using a previously undescribed strain with inducible control of sphingolipid synthesis. We characterized the biosynthetic pathway from myo-inositol-phosphate (MIP) synthesis to phosphoinositol dihydroceramide, determined the crystal structure of the recombinant BT MIP synthase enzyme and identified the phosphatase responsible for the conversion of bacterially-derived phosphatidylinositol phosphate (PIP-DAG) to phosphatidylinositol (PI-DAG). In vitro, loss of inositol lipid production altered BT capsule expression and antimicrobial peptide resistance. In vivo, loss of inositol lipids decreased bacterial fitness in a gnotobiotic mouse model. We identified a second putative, previously undescribed pathway for bacterial PI-DAG synthesis without a PIP-DAG intermediate, common in Prevotella. Our results indicate that inositol sphingolipid production is widespread in host-associated Bacteroidetes and has implications for symbiosis.

The pathways responsible for inositol lipid production in human gut Bacteroides are characterized and these lipids are important for capsule expression and antimicrobial peptide resistance in vitro and colonization in vivo.

Optimal structure determination from sub-optimal diffraction data

Herein we present the newest version of the HKL-3000 system that integrates data collection, data reduction, phasing, model building, refinement, and validation. The system significantly accelerates the process of structure determination and has proven its high value for the determination of very high-quality structures. The heuristic for choosing the best approach for every step of structure determination for various quality samples and diffraction data has been optimized. The latest modifications increase the likelihood of a successful structure determination with challenging data. The HKL-3000 is a successor of HKL and HKL-2000 programs. The use of the HKL family of programs has been reported for over 73,000 PDB deposits, that is, almost 50% of macromolecular structures determined with X-ray diffraction.

findMySequence: a neural-network-based approach for identification of unknown proteins in X-ray crystallography and cryo-EM

Although experimental protein-structure determination usually targets known proteins, chains of unknown sequence are often encountered. They can be purified from natural sources, appear as an unexpected fragment of a well characterized protein or appear as a contaminant. Regardless of the source of the problem, the unknown protein always requires characterization. Here, an automated pipeline is presented for the identification of protein sequences from cryo-EM reconstructions and crystallographic data. The method's application to characterize the crystal structure of an unknown protein purified from a snake venom is presented. It is also shown that the approach can be successfully applied to the identification of protein sequences and validation of sequence assignments in cryo-EM protein structures.

Modelling covalent linkages in CCP4

In this contribution, the current protocols for modelling covalent linkages within the CCP4 suite are considered. The mechanism used for modelling covalent linkages is reviewed: the use of dictionaries for describing changes to stereochemistry as a result of the covalent linkage and the application of link-annotation records to structural models to ensure the correct treatment of individual instances of covalent linkages. Previously, linkage descriptions were lacking in quality compared with those of contemporary component dictionaries. Consequently, AceDRG has been adapted for the generation of link dictionaries of the same quality as for individual components. The approach adopted by AceDRG for the generation of link dictionaries is outlined, which includes associated modifications to the linked components. A number of tools to facilitate the practical modelling of covalent linkages available within the CCP4 suite are described, including a new restraint-dictionary accumulator, the Make Covalent Link tool and AceDRG interface in Coot, the 3D graphical editor JLigand and the mechanisms for dealing with covalent linkages in the CCP4i2 and CCP4 Cloud environments. These integrated solutions streamline and ease the covalent-linkage modelling workflow, seamlessly transferring relevant information between programs. Current recommended practice is elucidated by means of instructive practical examples. By summarizing the different approaches to modelling linkages that are available within the CCP4 suite, limitations and potential pitfalls that may be encountered are highlighted in order to raise awareness, with the intention of improving the quality of future modelled covalent linkages in macromolecular complexes.

Advancements in macromolecular crystallography: from past to present

Protein Crystallography or Macromolecular Crystallography (MX) started as a new discipline of science with the pioneering work on the determination of the protein crystal structures by John Kendrew in 1958 and Max Perutz in 1960. The incredible achievements in MX are attributed to the development of advanced tools, methodologies, and automation in every aspect of the structure determination process, which have reduced the time required for solving protein structures from years to a few days, as evident from the tens of thousands of crystal structures of macromolecules available in PDB. The advent of brilliant synchrotron sources, fast detectors, and novel sample delivery methods has shifted the paradigm from static structures to understanding the dynamic picture of macromolecules; further propelled by X-ray Free Electron Lasers (XFELs) that explore the femtosecond regime. The revival of the Laue diffraction has also enabled the understanding of macromolecules through time-resolved crystallography. In this review, we present some of the astonishing method-related and technological advancements that have contributed to the progress of MX. Even with the rapid evolution of several methods for structure determination, the developments in MX will keep this technique relevant and it will continue to play a pivotal role in gaining unprecedented atomic-level details as well as revealing the dynamics of biological macromolecules. With many exciting developments awaiting in the upcoming years, MX has the potential to contribute significantly to the growth of modern biology by unraveling the mechanisms of complex biological processes as well as impacting the area of drug designing.

Crystal and solution structures reveal oligomerization of individual capsid homology domains of Drosophila Arc

Synaptic plasticity is vital for brain function and memory formation. One of the key proteins in long-term synaptic plasticity and memory is the activity-regulated cytoskeleton-associated protein (Arc). Mammalian Arc forms virus-like capsid structures in a process requiring the N-terminal domain and contains two C-terminal lobes that are structural homologues to retroviral capsids. Drosophila has two isoforms of Arc, dArc1 and dArc2, with low sequence similarity to mammalian Arc, but lacking a large N-terminal domain. Both dArc isoforms are related to the Ty3/gypsy retrotransposon capsid, consisting of N- and C-terminal lobes. Structures of dArc1, as well as capsids formed by both dArc isoforms, have been recently determined. We carried out structural characterization of the four individual dArc lobe domains. As opposed to the corresponding mammalian Arc lobe domains, which are monomeric, the dArc lobes were all oligomeric in solution, indicating a strong propensity for homophilic interactions. A truncated N-lobe from dArc2 formed a domain-swapped dimer in the crystal structure, resulting in a novel dimer interaction that could be relevant for capsid assembly or other dArc functions. This domain-swapped structure resembles the dimeric protein C of flavivirus capsids, as well as the structure of histones dimers, domain-swapped transcription factors, and membrane-interacting BAK domains. The strong oligomerization properties of the isolated dArc lobe domains explain the ability of dArc to form capsids in the absence of any large N-terminal domain, in contrast to the mammalian protein.

Specificity of AMPylation of the human chaperone BiP is mediated by TPR motifs of FICD

To adapt to fluctuating protein folding loads in the endoplasmic reticulum (ER), the Hsp70 chaperone BiP is reversibly modified with adenosine monophosphate (AMP) by the ER-resident Fic-enzyme FICD/HYPE. The structural basis for BiP binding and AMPylation by FICD has remained elusive due to the transient nature of the enzyme-substrate-complex. Here, we use thiol-reactive derivatives of the cosubstrate adenosine triphosphate (ATP) to covalently stabilize the transient FICD:BiP complex and determine its crystal structure. The complex reveals that the TPR-motifs of FICD bind specifically to the conserved hydrophobic linker of BiP and thus mediate specificity for the domain-docked conformation of BiP. Furthermore, we show that both AMPylation and deAMPylation of BiP are not directly regulated by the presence of unfolded proteins. Together, combining chemical biology, crystallography and biochemistry, our study provides structural insights into a key regulatory mechanism that safeguards ER homeostasis.

Structural and functional study of Legionella pneumophila effector RavA

Legionella pneumophila is an intracellular pathogen that causes Legionnaire's disease in humans. This bacterium can be found in freshwater environments as a free-living organism, but it is also an intracellular parasite of protozoa. Human infection occurs when inhaled aerosolized pathogen comes into contact with the alveolar mucosa and replicates in alveolar macrophages. Legionella enters the host cell by phagocytosis and redirects the Legionella-containing phagosomes from the phagocytic maturation pathway. These nascent phagosomes fuse with ER-derived secretory vesicles and membranes forming the Legionella-containing vacuole. Legionella subverts many host cellular processes by secreting over 300 effector proteins into the host cell via the Dot/Icm type IV secretion system. The cellular function for many Dot/Icm effectors is still unknown. Here, we present a structural and functional study of L. pneumophila effector RavA (Lpg0008). Structural analysis revealed that the RavA consists of four ~85 residue long α-helical domains with similar folds, which show only a low level of structural similarity to other protein domains. The ~90 residues long C-terminal segment is predicted to be natively unfolded. We show that during L. pneumophila infection of human cells, RavA localizes to the Golgi apparatus and to the plasma membrane. The same localization is observed when RavA is expressed in human cells. The localization signal resides within the C-terminal sequence C⁴⁰⁹ WTSFCGLF⁴¹⁷ . Yeast-two-hybrid screen using RavA as bait identified RAB11A as a potential binding partner. RavA is present in L. pneumophila strains but only distant homologs are found in other Legionella species, where the number of repeats varies.

Energetic and structural effects of the Tanford transition on ligand recognition of bovine β-lactoglobulin

Bovine β-lactoglobulin, an abundant protein in whey, is a promising nanocarrier for peroral administration of drug-like hydrophobic molecules, a process that involves transit through the different acidic conditions of the human digestive tract. Among the several pH-induced conformational rearrangements that this lipocalin undergoes, the Tanford transition is particularly relevant. This transition, which occurs with a midpoint around neutral pH, involves a conformational change of the E-F loop that regulates accessibility to the primary binding site. The effect of this transition on the ligand binding properties of this protein has scarcely been explored. In this study, we carried out an energetic and structural characterization of β-lactoglobulin molecular recognition at pH values above and below the zone in which the Tanford transition occurs. The combined analysis of crystallographic, calorimetric, and molecular dynamics data sheds new light on the interplay between self-association, ligand binding, and the Tanford pre- and post-transition conformational states, revealing novel aspects underlying the molecular recognition mechanism of this enigmatic lipocalin.

Phasertng: directed acyclic graphs for crystallographic phasing

Crystallographic phasing strategies increasingly require the exploration and ranking of many hypotheses about the number, types and positions of atoms, molecules and/or molecular fragments in the unit cell, each with only a small chance of being correct. Accelerating this move has been improvements in phasing methods, which are now able to extract phase information from the placement of very small fragments of structure, from weak experimental phasing signal or from combinations of molecular replacement and experimental phasing information. Describing phasing in terms of a directed acyclic graph allows graph-management software to track and manage the path to structure solution. The crystallographic software supporting the graph data structure must be strictly modular so that nodes in the graph are efficiently generated by the encapsulated functionality. To this end, the development of new software, Phasertng, which uses directed acyclic graphs natively for input/output, has been initiated. In Phasertng, the codebase of Phaser has been rebuilt, with an emphasis on modularity, on scripting, on speed and on continuing algorithm development. As a first application of phasertng, its advantages are demonstrated in the context of phasertng.xtricorder, a tool to analyse and triage merged data in preparation for molecular replacement or experimental phasing. The description of the phasing strategy with directed acyclic graphs is a generalization that extends beyond the functionality of Phasertng, as it can incorporate results from bioinformatics and other crystallographic tools, and will facilitate multifaceted search strategies, dynamic ranking of alternative search pathways and the exploitation of machine learning to further improve phasing strategies.

Crystal structure of the SPRY domain-containing protein 7 reveals unique structural features

The SPRY/B30.2 domain is one of the most abundant protein domains found in eukaryotes. Vast majority of the SPRY domain-containing proteins are multi-domain proteins. The SPRY domain-containing protein 7 (SPRY7, also named C13orf1, and named chronic lymphocytic leukemia deletion region gene 6 protein, CCLD6, encoded by the spryd7 gene) is the smallest SPRY domain protein in human that does not contain other accessory domains. Here we have determined the crystal structure of human SPRY7 at a resolution of 1.62 Å and found that SPRY7 has some unique structural features that are not present in other previously reported SRPY domain structures. Overall, SPRY7 may represent an evolutionary early version of the SPRY domain, and subsequent loop insertions and expansions, residue substitutions, as well as domain combinations have rendered the SPRY domain versatile binding specificities and broad biological functions. These results serve as a useful basis for a profound characterization of the molecular interactions of SPRY7.

Predicting protein model correctness in Coot using machine learning

Manually identifying and correcting errors in protein models can be a slow process, but improvements in validation tools and automated model-building software can contribute to reducing this burden. This article presents a new correctness score that is produced by combining multiple sources of information using a neural network. The residues in 639 automatically built models were marked as correct or incorrect by comparing them with the coordinates deposited in the PDB. A number of features were also calculated for each residue using Coot, including map-to-model correlation, density values, B factors, clashes, Ramachandran scores, rotamer scores and resolution. Two neural networks were created using these features as inputs: one to predict the correctness of main-chain atoms and the other for side chains. The 639 structures were split into 511 that were used to train the neural networks and 128 that were used to test performance. The predicted correctness scores could correctly categorize 92.3% of the main-chain atoms and 87.6% of the side chains. A Coot ML Correctness script was written to display the scores in a graphical user interface as well as for the automatic pruning of chains, residues and side chains with low scores. The automatic pruning function was added to the CCP4i2 Buccaneer automated model-building pipeline, leading to significant improvements, especially for high-resolution structures.

The use of local structural similarity of distant homologues for crystallographic model building from a molecular-replacement solution

The performance of automated protein model building usually decreases with resolution, mainly owing to the lower information content of the experimental data. This calls for a more elaborate use of the available structural information about macromolecules. Here, a new method is presented that uses structural homologues to improve the quality of protein models automatically constructed using ARP/wARP. The method uses local structural similarity between deposited models and the model being built, and results in longer main-chain fragments that in turn can be more reliably docked to the protein sequence. The application of the homology-based model extension method to the example of a CFA synthase at 2.7 Å resolution resulted in a more complete model with almost all of the residues correctly built and docked to the sequence. The method was also evaluated on 1493 molecular-replacement solutions at a resolution of 4.0 Å and better that were submitted to the ARP/wARP web service for model building. A significant improvement in the completeness and sequence coverage of the built models has been observed.

Using Phaser and ensembles to improve the performance of SIMBAD

The conventional approach to search-model identification in molecular replacement (MR) is to screen a database of known structures using the target sequence. However, this strategy is not always effective, for example when the relationship between sequence and structural similarity fails or when the crystal contents are not those expected. An alternative approach is to identify suitable search models directly from the experimental data. SIMBAD is a sequence-independent MR pipeline that uses either a crystal lattice search or MR functions to directly locate suitable search models from databases. The previous version of SIMBAD used the fast AMoRe rotation-function search. Here, a new version of SIMBAD which makes use of Phaser and its likelihood scoring to improve the sensitivity of the pipeline is presented. It is shown that the additional compute time potentially required by the more sophisticated scoring is counterbalanced by the greater sensitivity, allowing more cases to trigger early-termination criteria, rather than running to completion. Using Phaser solved 17 out of 25 test cases in comparison to the ten solved with AMoRe, and it is shown that use of ensemble search models produces additional performance benefits.

Facilities that make the PDB data collection more powerful

We describe a series of databases and tools that directly or indirectly support biomedical research on macromolecules, with focus on their applicability in protein structure bioinformatics research. DSSP, that determines secondary structures of proteins, has been updated to work well with extremely large structures in multiple formats. The PDBREPORT database that lists anomalies in protein structures has been remade to remove many small problems. These reports are now available as PDF-formatted files with a computer-readable summary. The VASE software has been added to analyze and visualize HSSP multiple sequence alignments for protein structures. The Lists collection of databases has been extended with a series of databases, most noticeably with a database that gives each protein structure a grade for usefulness in protein structure bioinformatics projects. The PDB-REDO collection of reanalyzed and re-refined protein structures that were solved by X-ray crystallography has been improved by dealing better with sugar residues and with hydrogen bonds, and adding many missing surface loops. All academic software underlying these protein structure bioinformatics applications and databases are now publicly accessible, either directly from the authors or from the GitHub software repository.

Elucidation of a sialic acid metabolism pathway in mucus-foraging Ruminococcus gnavus unravels mechanisms of bacterial adaptation to the gut

Sialic acid (Neu5Ac) is commonly found in terminal location of colonic mucins glycans where it is a much-coveted nutrient for gut bacteria including Ruminococcus gnavus. R. gnavus is part of the healthy gut microbiota in humans but shows a disproportionate representation in diseases. There is therefore a need in understanding the molecular mechanisms underpinning its adaptation to the gut. Previous in vitro work demonstrated that R. gnavus mucin glycan-foraging strategy is strain-dependent and associated with the expression of an intramolecular trans-sialidase releasing 2,7-anhydro-Neu5Ac instead of Neu5Ac from mucins. Here, we have unravelled the metabolism pathway of 2,7-anhydro-Neu5Ac in R. gnavus which is underpinned by the exquisite specificity of the sialic transporter for 2,7-anhydro-Neu5Ac, and by the action of an oxidoreductase converting 2,7-anhydro-Neu5Ac into Neu5Ac which then becomes substrate of a Neu5Ac-specific aldolase. Having generated a R. gnavus nan cluster deletion mutant that lost the ability to grow on sialylated substrates, we showed that in gnotobiotic mice colonised with R. gnavus wild-type and mutant strains, the fitness of the nan mutant was significantly impaired with a reduced ability to colonise the mucus layer. Overall, our study revealed a unique sialic acid pathway in bacteria, with significant implications for the spatial adaptation of mucin-foraging gut symbionts in health and disease.

Structural Evidence for an Octameric Ring Arrangement of SARM1

SARM1 induces axonal degeneration in response to various insults and is therefore considered an attractive drug target for the treatment of neuro-degenerative diseases as well as for brain and spinal cord injuries. SARM1 activity depends on the integrity of the protein's SAM domains, as well as on the enzymatic conversion of NAD+ to ADPR (ADP Ribose) products by the SARM1's TIR domain. Therefore, inhibition of either SAM or TIR functions may constitute an effective therapeutic strategy. However, there is currently no SARM1-directed therapeutic approach available because of an insufficient structural and mechanistic understanding of this protein. In this study, we found that SARM1 assembles into an octameric ring. This arrangement was not described before in other SAM proteins, but is reminiscent of the apoptosome and inflammasome-well-known apoptotic ring-like oligomers. We show that both SARM1 and the isolated tandem SAM^1-2 domains form octamers in solution, and electron microscopy analysis reveals an octameric ring of SARM1. We determined the crystal structure of SAM^1-2 and found that it also forms a closed octameric ring in the crystal lattice. The SAM^1-2 ring interactions are mediated by complementing "lock and key" hydrophobic grooves and inserts and electrostatic charges between the neighboring protomers. We have mutated several interacting SAM^1-2 interfaces and measured how these mutations affect SARM1 apoptotic activity in cultured cells, and in this way identified critical oligomerization sites that facilitate cell death. These results highlight the importance of oligomerization for SARM1 function and reveal critical epitopes for future targeted drug development.

Announcing mandatory submission of PDBx/mmCIF format files for crystallographic depositions to the Protein Data Bank (PDB)

This letter announces that PDBx/mmCIF format files will become mandatory for crystallographic depositions to the Protein Data Bank (PDB).

Structure of mammalian plasma fetuin-B and its mechanism of selective metallopeptidase inhibition

Mammalian fetuin-A and fetuin-B are abundant serum proteins with pleiotropic functions. Fetuin-B is a highly selective and potent inhibitor of metallo-peptidases (MPs) of the astacin family, which includes ovastacin in mammals. By inhibiting ovastacin, fetuin-B is essential for female fertility. The crystal structure of fetuin-B was determined unbound and in complex with archetypal astacin, and it was found that the inhibitor has tandem cystatin-type modules (CY1 and CY2). They are connected by an exposed linker with a rigid, disulfide-linked 'CPDCP-trunk', and are followed by a C-terminal region (CTR) with little regular secondary structure. The CPDCP-trunk and a hairpin of CY2 form a bipartite wedge, which slots into the active-site cleft of the MP. These elements occupy the nonprimed and primed sides of the cleft, respectively, but spare the specificity pocket so that the inhibitor is not cleaved. The aspartate in the trunk blocks the catalytic zinc of astacin, while the CY2 hairpin binds through a QWVXGP motif. The CY1 module assists in structural integrity and the CTR is not involved in inhibition, as verified by in vitro studies using a cohort of mutants and variants. Overall, the inhibition conforms to a novel 'raised-elephant-trunk' mechanism for MPs, which is reminiscent of single-domain cystatins that target cysteine peptidases. Over 200 sequences from vertebrates have been annotated as fetuin-B, underpinning its ubiquity and physiological relevance; accordingly, sequences with conserved CPDCP- and QWVXGP-derived motifs have been found from mammals to cartilaginous fishes. Thus, the raised-elephant-trunk mechanism is likely to be generally valid for the inhibition of astacins by orthologs of fetuin-B.

SIMBAD: a sequence-independent molecular-replacement pipeline

The conventional approach to finding structurally similar search models for use in molecular replacement (MR) is to use the sequence of the target to search against those of a set of known structures. Sequence similarity often correlates with structure similarity. Given sufficient similarity, a known structure correctly positioned in the target cell by the MR process can provide an approximation to the unknown phases of the target. An alternative approach to identifying homologous structures suitable for MR is to exploit the measured data directly, comparing the lattice parameters or the experimentally derived structure-factor amplitudes with those of known structures. Here, SIMBAD, a new sequence-independent MR pipeline which implements these approaches, is presented. SIMBAD can identify cases of contaminant crystallization and other mishaps such as mistaken identity (swapped crystallization trays), as well as solving unsequenced targets and providing a brute-force approach where sequence-dependent search-model identification may be nontrivial, for example because of conformational diversity among identifiable homologues. The program implements a three-step pipeline to efficiently identify a suitable search model in a database of known structures. The first step performs a lattice-parameter search against the entire Protein Data Bank (PDB), rapidly determining whether or not a homologue exists in the same crystal form. The second step is designed to screen the target data for the presence of a crystallized contaminant, a not uncommon occurrence in macromolecular crystallography. Solving structures with MR in such cases can remain problematic for many years, since the search models, which are assumed to be similar to the structure of interest, are not necessarily related to the structures that have actually crystallized. To cater for this eventuality, SIMBAD rapidly screens the data against a database of known contaminant structures. Where the first two steps fail to yield a solution, a final step in SIMBAD can be invoked to perform a brute-force search of a nonredundant PDB database provided by the MoRDa MR software. Through early-access usage of SIMBAD, this approach has solved novel cases that have otherwise proved difficult to solve.

Recent developments in MrBUMP: better search-model preparation, graphical interaction with search models, and solution improvement and assessment

Increasing sophistication in molecular-replacement (MR) software and the rapid expansion of the PDB in recent years have allowed the technique to become the dominant method for determining the phases of a target structure in macromolecular X-ray crystallography. In addition, improvements in bioinformatic techniques for finding suitable homologous structures for use as MR search models, combined with developments in refinement and model-building techniques, have pushed the applicability of MR to lower sequence identities and made weak MR solutions more amenable to refinement and improvement. MrBUMP is a CCP4 pipeline which automates all stages of the MR procedure. Its scope covers everything from the sourcing and preparation of suitable search models right through to rebuilding of the positioned search model. Recent improvements to the pipeline include the adoption of more sensitive bioinformatic tools for sourcing search models, enhanced model-preparation techniques including better ensembling of homologues, and the use of phase improvement and model building on the resulting solution. The pipeline has also been deployed as an online service through CCP4 online, which allows its users to exploit large bioinformatic databases and coarse-grained parallelism to speed up the determination of a possible solution. Finally, the molecular-graphics application CCP4mg has been combined with MrBUMP to provide an interactive visual aid to the user during the process of selecting and manipulating search models for use in MR. Here, these developments in MrBUMP are described with a case study to explore how some of the enhancements to the pipeline and to CCP4mg can help to solve a difficult case.