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

Crystal Structure of Caryolan-1-ol Synthase, a Sesquiterpene Synthase Catalyzing an Initial Anti-Markovnikov Cyclization Reaction

In a continuing effort to understand reaction mechanisms of terpene synthases catalyzing initial anti-Markovnikov cyclization reactions, we solved the X-ray crystal structure of (+)-caryolan-1-ol synthase (CS) from Streptomyces griseus , with and without an inactive analog of the FPP substrate, 2-fluorofarnesyl diphosphate (2FFPP), bound in the active site of the enzyme. The CS-2FFPP complex was solved to 2.65 Å resolution and showed the ligand in a linear, elongated orientation, incapable of undergoing the initial cyclization event to form a bond between carbons C1 and C11. Intriguingly, the apo CS structure (2.2 Å) also had electron density in the active site, in this case density that was well fit with a curled-up tetraethylene glycol molecule presumably recruited from the crystallization medium. The density was also well fit by a molecule of farnesene suggesting that the structure may mimic an intermediate along the reaction coordinate. The curled-up conformation of tetraethylene glycol was accompanied by dramatic rotamer shifts among active-site residues. Most notably, W56 was observed to undergo a 90° rotation between the 2FFPP complex and apo-enzyme structures, suggesting that it contributes to steric interactions that help curl the tetraethylene glycol molecule in the active site, and by extension perhaps also a derivative of the FPP substrate in the normal course of the cyclization reaction. In support of this proposal, the CS W56L variant lost the ability to cyclize the FPP substrate and produced only the linear terpene products farnesol and α- and β-farnesene.

Modes and model building in SHELXE

Density modification is a standard step to provide a route for routine structure solution by any experimental phasing method, with single-wavelength or multi-wavelength anomalous diffraction being the most popular methods, as well as to extend fragments or incomplete models into a full solution. The effect of density modification on the starting maps from either source is illustrated in the case of SHELXE. The different modes in which the program can run are reviewed; these include less well known uses such as reading external phase values and weights or phase distributions encoded in Hendrickson-Lattman coefficients. Typically in SHELXE, initial phases are calculated from experimental data, from a partial model or map, or from a combination of both sources. The initial phase set is improved and extended by density modification and, if the resolution of the data and the type of structure permits, polyalanine tracing. As a feature to systematically eliminate model bias from phases derived from predicted models, the trace can be set to exclude the area occupied by the starting model. The trace now includes an extension into the gamma position or hydrophobic and aromatic side chains if a sequence is provided, which is performed in every tracing cycle. Once a correlation coefficient of over 30% between the structure factors calculated from such a trace and the native data indicates that the structure has been solved, the sequence is docked in all model-building cycles and side chains are fitted if the map supports it. The extensions to the tracing algorithm brought in to provide a complete model are discussed. The improvement in phasing performance is assessed using a set of tests.

Archaeal GPN-loop GTPases involve a lock-switch-rock mechanism for GTP hydrolysis

IMPORTANCE

GPN-loop GTPases have been found to be crucial for eukaryotic RNA polymerase II assembly and nuclear trafficking. Despite their ubiquitous occurrence in eukaryotes and archaea, the mechanism by which these GTPases mediate their function is unknown. Our study on an archaeal representative from Sulfolobus acidocaldarius showed that these dimeric GTPases undergo large-scale conformational changes upon GTP hydrolysis, which can be summarized as a lock-switch-rock mechanism. The observed requirement of SaGPN for motility appears to be due to its large footprint on the archaeal proteome.

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.

Structure of the photoreceptor synaptic assembly of the extracellular matrix protein pikachurin with the orphan receptor GPR179

Precise synapse formation is essential for normal functioning of the nervous system. Retinal photoreceptors establish selective contacts with bipolar cells, aligning the neurotransmitter release apparatus with postsynaptic signaling cascades. This involves transsynaptic assembly between the dystroglycan-dystrophin complex on the photoreceptor and the orphan receptor GPR179 on the bipolar cell, which is mediated by the extracellular matrix protein pikachurin (also known as EGFLAM). This complex plays a critical role in the synaptic organization of photoreceptors and signal transmission, and mutations affecting its components cause blinding disorders in humans. Here, we investigated the structural organization and molecular mechanisms by which pikachurin orchestrates transsynaptic assembly and solved structures of the human pikachurin domains by x-ray crystallography and of the GPR179-pikachurin complex by single-particle, cryo-electron microscopy. The structures reveal molecular recognition principles of pikachurin by the Cache domains of GPR179 and show how the interaction is involved in the transsynaptic alignment of the signaling machinery. Together, these data provide a structural basis for understanding the synaptic organization of photoreceptors and ocular pathology.

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 staphylococcal inhibitory protein SPIN binds to human myeloperoxidase with picomolar affinity but only dampens halide oxidation

The heme enzyme myeloperoxidase (MPO) is one of the key players in the neutrophil-mediated killing of invading pathogens as part of the innate immune system. MPO generates antimicrobial oxidants which indiscriminately and effectively kill phagocytosed pathogens. Staphylococcus aureus however is able to escape this fate, in part by secreting a small protein called SPIN (Staphylococcal Peroxidase Inhibitor), which specifically targets and inhibits MPO in a structurally complex manner. Here we present the first crystal structures of the complex of SPIN-aureus and a truncated version (SPIN-truncated) with mature dimeric leukocyte MPO. We unravel the contributions of the two domains to the kinetics and thermodynamics of SPIN-aureus binding to MPO by using a broad array of complementary biochemical and biophysical methods. The C-terminal "recognition" domain is shown to mediate specific binding to MPO, while interaction of the N-terminal "inhibitory" domain is guided mainly by hydrophobic effects and thus is less sequence-dependent. We found that inhibition of MPO is achieved by reducing substrate migration, but SPIN-aureus cannot completely block MPO activity. Its' effectiveness is inversely related to substrate size, with no discernible dependence on other factors. Thus, SPIN-aureus is an extremely high-affinity inhibitor and highly efficient for substrates larger than halogens. As aberrant MPO activity is implicated in a number of chronic inflammatory diseases, SPIN-aureus is the first promising protein inhibitor for specific inhibition of human MPO.

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%.

Structural basis of synergistic neutralization of Crimean-Congo hemorrhagic fever virus by human antibodies

Crimean-Congo hemorrhagic fever virus (CCHFV) is the most widespread tick-borne zoonotic virus, with a 30% case fatality rate in humans. Structural information is lacking in regard to the CCHFV membrane fusion glycoprotein Gc—the main target of the host neutralizing antibody response—as well as antibody–mediated neutralization mechanisms. We describe the structure of prefusion Gc bound to the antigen-binding fragments (Fabs) of two neutralizing antibodies that display synergy when combined, as well as the structure of trimeric, postfusion Gc. The structures show the two Fabs acting in concert to block membrane fusion, with one targeting the fusion loops and the other blocking Gc trimer formation. The structures also revealed the neutralization mechanism of previously reported antibodies against CCHFV, providing the molecular underpinnings essential for developing CCHFV–specific medical countermeasures for epidemic preparedness.

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.

Towards the automatic crystal structure solution of nucleic acids: automated model building using the new CAB program

CAB, a recently described automated model-building (AMB) program, has been modified to work effectively with nucleic acids. To this end, several new algorithms have been introduced and the libraries have been updated. To reduce the input average phase error, ligand heavy atoms are now located before starting the CAB interpretation of the electron-density maps. Furthermore, alternative approaches are used depending on whether the ligands belong to the target or to the model chain used in the molecular-replacement step. Robust criteria are then applied to decide whether the AMB model is acceptable or whether it must be modified to fit prior information on the target structure. In the latter case, the model chains are rearranged to fit prior information on the target chains. Here, the performance of the new AMB program CAB applied to various nucleic acid structures is discussed. Other well documented programs such as Nautilus, ARP/wARP and phenix.autobuild were also applied and the experimental results are described.

The accuracy of protein models automatically built into cryo-EM maps with ARP/wARP

A new module of the ARP/wARP suite for automated model building into cryo-EM maps is presented.

Recent developments in cryogenic electron microscopy (cryo-EM) have enabled structural studies of large macromolecular complexes at resolutions previously only attainable using macromolecular crystallography. Although a number of methods can already assist in de novo building of models into high-resolution cryo-EM maps, automated and reliable map interpretation remains a challenge. Presented here is a systematic study of the accuracy of models built into cryo-EM maps using ARP/wARP. It is demonstrated that the local resolution is a good indicator of map interpretability, and for the majority of the test cases ARP/wARP correctly builds 90% of main-chain fragments in regions where the local resolution is 4.0 Å or better. It is also demonstrated that the coordinate accuracy for models built into cryo-EM maps is comparable to that of X-ray crystallographic models at similar local cryo-EM and crystallographic resolutions. The model accuracy also correlates with the refined atomic displacement parameters.

Pairwise running of automated crystallographic model-building pipelines

For the last two decades, researchers have worked independently to automate protein model building, and four widely used software pipelines have been developed for this purpose: ARP/wARP, Buccaneer, Phenix AutoBuild and SHELXE. Here, the usefulness of combining these pipelines to improve the built protein structures by running them in pairwise combinations is examined. The results show that integrating these pipelines can lead to significant improvements in structure completeness and Rfree. In particular, running Phenix AutoBuild after Buccaneer improved structure completeness for 29% and 75% of the data sets that were examined at the original resolution and at a simulated lower resolution, respectively, compared with running Phenix AutoBuild on its own. In contrast, Phenix AutoBuild alone produced better structure completeness than the two pipelines combined for only 7% and 3% of these data sets.

ALEPH: a network-oriented approach for the generation of fragment-based libraries and for structure interpretation

The analysis of large structural databases reveals general features and relationships among proteins, providing useful insight. A different approach is required to characterize ubiquitous secondary-structure elements, where flexibility is essential in order to capture small local differences. The ALEPH software is optimized for the analysis and the extraction of small protein folds by relying on their geometry rather than on their sequence. The annotation of the structural variability of a given fold provides valuable information for fragment-based molecular-replacement methods, in which testing alternative model hypotheses can succeed in solving difficult structures when no homology models are available or are successful. ARCIMBOLDO_BORGES combines the use of composite secondary-structure elements as a search model with density modification and tracing to reveal the rest of the structure when both steps are successful. This phasing method relies on general fold libraries describing variations around a given pattern of β-sheets and helices extracted using ALEPH. The program introduces characteristic vectors defined from the main-chain atoms as a way to describe the geometrical properties of the structure. ALEPH encodes structural properties in a graph network, the exploration of which allows secondary-structure annotation, decomposition of a structure into small compact folds, generation of libraries of models representing a variation of a given fold and finally superposition of these folds onto a target structure. These functions are available through a graphical interface designed to interactively show the results of structure manipulation, annotation, fold decomposition, clustering and library generation. ALEPH can produce pictures of the graphs, structures and folds for publication purposes.

SEQUENCE SLIDER: expanding polyalanine fragments for phasing with multiple side-chain hypotheses

Fragment-based molecular-replacement methods can solve a macromolecular structure quasi-ab initio. ARCIMBOLDO, using a common secondary-structure or tertiary-structure template or a library of folds, locates these with Phaser and reveals the rest of the structure by density modification and autotracing in SHELXE. The latter stage is challenging when dealing with diffraction data at lower resolution, low solvent content, high β-sheet composition or situations in which the initial fragments represent a low fraction of the total scattering or where their accuracy is low. SEQUENCE SLIDER aims to overcome these complications by extending the initial polyalanine fragment with side chains in a multisolution framework. Its use is illustrated on test cases and previously unknown structures. The selection and order of fragments to be extended follows the decrease in log-likelihood gain (LLG) calculated with Phaser upon the omission of each single fragment. When the starting substructure is derived from a remote homolog, sequence assignment to fragments is restricted by the original alignment. Otherwise, the secondary-structure prediction is matched to that found in fragments and traces. Sequence hypotheses are trialled in a brute-force approach through side-chain building and refinement. Scoring the refined models through their LLG in Phaser may allow discrimination of the correct sequence or filter the best partial structures for further density modification and autotracing. The default limits for the number of models to pursue are hardware dependent. In its most economic implementation, suitable for a single laptop, the main-chain trace is extended as polyserine rather than trialling models with different sequence assignments, which requires a grid or multicore machine. SEQUENCE SLIDER has been instrumental in solving two novel structures: that of MltC from 2.7 Å resolution data and that of a pneumococcal lipoprotein with 638 residues and 35% solvent content.