key: cord-0921682-nbwlscxu
authors: Flower, Thomas G.; Hurley, James H.
title: Crystallographic molecular replacement using an in silico‐generated search model of SARS‐CoV‐2 ORF8
date: 2021-03-04
journal: Protein Sci
DOI: 10.1002/pro.4050
sha: 709f9238714ad3f7c5053466aef544514ebb0a86
doc_id: 921682
cord_uid: nbwlscxu

The majority of crystal structures are determined by the method of molecular replacement (MR). The range of application of MR is limited mainly by the need for an accurate search model. In most cases, pre‐existing experimentally determined structures are used as search models. In favorable cases, ab initio predicted structures have yielded search models adequate for MR. The ORF8 protein of SARS‐CoV‐2 represents a challenging case for MR using an ab initio prediction because ORF8 has an all β‐sheet fold and few orthologs. We previously determined experimentally the structure of ORF8 using the single anomalous dispersion (SAD) phasing method, having been unable to find an MR solution to the crystallographic phase problem. Following a report of an accurate prediction of the ORF8 structure, we assessed whether the predicted model would have succeeded as an MR search model. A phase problem solution was found, and the resulting structure was refined, yielding structural parameters equivalent to the original experimental solution.

Molecular replacement (MR) is an in silico technique that provides phase information required to solve macromolecular crystal structures. MR relies on the existence of an experimentally determined structure, known as the "search model," which is similar to the target. Here we show that an ab initio, SARS CoV-2 ORF8 protein model, generated by Google DeepMind's AlphaFold team, is sufficiently accurate to provide a phase solution by MR, bypassing the need for an experimentally determined search model.

Two key pieces of information are required to determine a macromolecular structure from a crystallographic diffraction experiment, namely the amplitudes and phases of the diffracted waves. The amplitudes of the diffracted waves are calculated directly from the measured intensities of the scattered waves, while the phase information is lost. 1 The crystallographic phase problem remains a substantial bottleneck in macromolecular structure determination. In the majority of cases, the phase problem is overcome using an in silico method known as molecular replacement (MR). [2] [3] [4] [5] In MR, a related structure, known as the search model, is used to provide initial phaseestimates for the target structure. While the advent of MR has rendered phase determination near-trivial in many cases, it was, until recently, limited to circumstances where structures of homologous proteins already existed.

In principle, advances in in silico protein structure prediction could provide MR search models sufficiently similar to the target structures so as to produce a phase solution. [6] [7] [8] Such an advance could, in theory, forever bypass the need for experimentally derived phases. The use of ab initio models as MR search models has seen some success, principally with α helix-rich proteins. 9 Accurate ab initio prediction of β-rich folds has long proven to be exceptionally challenging 10 due, in part, to the high proportion of non-local interactions between β-strands. To our knowledge only one β-rich crystal structure has previously been phased using an ab initio generated search model. 11 This approach relied on a highly parallelized trial and error approach where hundreds of input MR ensemble search models were tested. 10 Recent advances in ab initio modeling have come from the long-accepted principle that evolutionary covariance of residues can aid in inter-residue contact prediction. This relies on the fact that deleterious point mutations are often paired with compensatory ones during evolutionary development. Multiple sequence alignments (MSA) of many related protein sequences can be used to identify these correlated mutations, with several systems applying neural networks to do so. [12] [13] [14] [15] [16] Many groups, including Google DeepMind's "AlphaFold," have taken this principle one step further by predicting specified distances between residue pairs which provide more information about the structure than contact predictions alone. [16] [17] [18] [19] [20] In the case of AlphaFold, these inter-residue distances, along with backbone torsion angles, are predicted in a first step using convolutional neural networks. 16 They are then provided as a target in gradient descent algorithm, which aims to bring the 3D structure as close as possible to these predicted distance and torsion angles. These principles were applied to great effect during the 14th biennial Critical Assessment of protein Structure Prediction (CASP) competition, where the AlphaFold2 team contributed a wealth of ab initio predictions, the majority of which were considered highly similar to experimentally determined protein structures with a median Global Distance Test (GDT) score of 92.4 out of 100. 21, 22 Our group recently determined the crystal structure of ORF8 at 2.04 Å, 23 a SARS-CoV-2 protein which has been implicated in immune evasion. The structure revealed that the protein is composed entirely of β-strands and unstructured regions. While a full native dataset was collected during the early stages of the project, a lack of suitable search model made phase determination by MR unfeasible. Initial phases were obtained by anomalous dispersion.

Here we assessed whether a template-free, ab initio protein model, generated by the AlphaFold2 group, was of sufficient quality to phase the native ORF8 dataset by MR. No truncation of the model was required nor was there any need to provide an ensemble of search models. Not only is this approach likely to prove useful for future structural determination campaigns where a homologous structure is not available but could aid in the determination of preexisting "unsolvable" datasets.

Superposition of the ab initio AlphaFold2 ORF8 model (CASP14 model ID: T1064TS427_1-D1) with the previously deposited crystal structure protomer (PDBID: 7JTL) showed they are highly similar, exhibiting all atom RMSDs of 1.4 Å and 1.6 Å with chains A and B respectively ( Figure 1 ). This suggests the potential of the AlphaFold2 model as an appropriate MR search model candidate.

To test this, MR was performed using native ORF8 structure factor amplitudes (PDBID: 7JTL) and the unaltered AlphaFold2 ORF8 prediction as a search model. A single MR solution was identified with two copies in the asymmetric unit. The Log Likelihood Gain (LLG) was 167 with values of 120 or greater indicating that a correct solution has almost certainly been found. 24 The MR solution places the search model at the correct position within the crystal lattice ( Figure 2 ). F I G U R E 1 Superposition of experimentally determined CoV-2 ORF8 structure with AlphaFold2 prediction. Chains A and B of the ORF8 crystal structure (PDBID: 7JTL) are superposed and colored yellow and orange respectively. The AlphaFold2 prediction is colored cyan. N and C termini are labeled accordingly. The ORF8-specific region is highlighted in pink A 2Fo-Fc electron density map was generated using phases from the placed but un-modified AlphaFold2 model and the native structure factor amplitudes ( Figure 3 ). The density was of sufficient quality to allow the majority of sidechains and unmodeled main chain regions to be built unambiguously. The majority of sidechain positions in the starting model proved to be remarkably close to the final structure (Figure 4a,b) . Other map features were present that were in poor agreement with the input model (Figure 3b ,c-f). These differences exclude the possibility that the density is dominated by input model-based phase-bias. Such differences include incorrect side-chain positioning/orientation, minor main chain deviations and the presence of unbuilt/unmodeled regions (Figures 3 and 4) . Of the side-chains that were included in the AlphaFold2 model, 18.5% adopted rotamer conformations that were not consistent with the experimentally determined structures (PDBIDs: 7JTL, 7JX6). 23 This value increased to 26.5% when unmodeled regions were considered as having incorrect rotamer assignment. Iterative rounds of manual model building and refinement produced a final structure determination with global quality metrics on par with those of the previously deposited structure ( Table 1) .

The region corresponding to residues 62-74 was absent from the AlphaFold2 prediction. This missing region is F I G U R E 2 Position of two placed CoV-2 ORF8 AlphaFold2 search models following MR. MR was performed using native ORF8 diffraction data. The first and second placed copies of the search model are colored cyan and purple respectively. The positions of chains A and B of the experimentally-phased, previously determined ORF8 crystal structure (PDBID: 7JTL) asymmetric unit are shown for reference and are colored yellow and orange respectively F I G U R E 3 Representative regions of 2Fo-Fc electron density following MR with AlphaFold2 CoV-2 ORF8 search model. Panel (a) shows a region of density that is in good agreement with the unmodified AlphaFold2 search model. Panel (b) provides an example where there is an obvious discrepancy between the experimental data and the search model, suggesting that the position and orientation of the Arg48 side-chain should be modified. Map is contoured at 1.2 σ and represented as a grey mesh located within a~35 amino acid insertion known as the "ORF8-specific region" which has been proposed to be involved in the formation of higher-order assemblies ( Figure 1 ). 23 This majority of this region is well ordered in both our structure (PDBID: 7JTL) and a subsequently determined crystal structure produced by another group (PDBID: 7JX6). Despite being ordered, this region is largely devoid of secondary structure with each individual chain adopting a unique conformation in each of the two molecules within the asymmetric unit, suggesting that it is dynamic and context-dependent in nature.

MR was systematically attempted with the nine nextbest ranked ORF8 predictions that were submitted to CASP14 (Table 2 ). Other than the AlphaFold2 prediction, none of the top 10 ranking models produced a convincing, single, MR solution. Top Phaser LLG scores ranged from 43.0 to 45.7 and visual comparison with the experimentally phased model (PDBID: 7JTL) confirmed that they were incorrectly placed, with the resulting electron densities dominated by input model phase bias.

Immunoglobulin domains from the Drosophila neural receptor Dscam1 (PDBID: 4X83) and Murine Natural Killer Cell Receptor 2B4 (PDBID: 2PTU) 25,26 were retrospectively identified as the most-similar existing PDB entries to chain A (RMSD: 3.4 Å, sequence identity: 11%) and chain B (RMSD: 3.7 Å, sequence identity 9%) of SARS CoV-2 ORF8 (PDBID: 7JTL) respectively (Table 3) . These models also failed to produce convincing MR phase-solutions, with top LLG scores of 43.4 and 44.1 for Dscam1 and Natural Killer Cell Receptor 2B4 respectively. Visual comparison of these top-ranking solutions with the completed ORF8 structure (PDBID: 7JTL) confirmed that a correct solution had not been found. and AlphaFold prediction are colored yellow and cyan respectively. 2Fo-Fc electron density was calculated using the previously deposited ORF8 coordinates and structure factor amplitudes (PDBID: 7JTL), contoured at 1.2 σ and represented as a grey mesh. Notable residues are labeled accordingly

The success of the AlphaFold2 predictions at CASP14 has led to considerable excitement and discussion of the relative roles of ab initio prediction and experimental structure determination in the future. 22 Our own assessment found that the AlphaFold2 predicted structure of SARS-CoV-2 ORF8 was of sufficient quality to yield a correct MR solution. This represents both a stringent test of accuracy and bodes well for the practical utility of AlphaFold2 predictions in MR. It will clearly be important for the AlphaFold2 method, and the computing power needed to support it, to become available to the structural biology community.

Despite the impressive accuracy of the fold assignment and backbone structure in the AlphaFold2 prediction, critical details were missing and could only be resolved by crystallographic structure solution. Side-chain conformations were often inaccurate in the predicted structure. Moreover, dimerization of ORF8 is thought to be important for function, yet the AlphaFold2 prediction only provided the structure of the monomer. Even before AlphaFold2, high quality ab initio structure predictions were accurate enough to reduce, if not eliminate, the need for experimental determinations of the overall fold of single domain proteins. 27 In our view, the real utility of the "accuracy revolution" in structure prediction will be to increase synergy with experimentation. Not only will this be invaluable for crystallographic MR, as assessed here, but in the interpretation of density in cryo-EM as well.

The AlphaFold2 ORF8 prediction was obtained from the CASP14 website (www.predictioncenter.org; model ID: T1064TS427_1-D1) and prepared for Molecular Replacement using the Phenix software suite. 28 Molecular replacement was performed using the program Phaser with default search parameters. The prepared AlphaFold2 ORF8 prediction was provided as a search model and native ORF8 structure factors as experimental data. The search was limited to the known space group P4 1 2 1 2. A single solution was obtained with two copies in the asymmetric unit, an LLG value of 167 and TFZ score of 14.7 indicating that the correct solution had been found. Identical search parameters were used for all other search models.

Iterative sounds of manual model building and refinement were performed using Coot 29 and Phenix Refine 28 respectively (for statistics see Table 1 ). Figures were produced using the program PyMOL (https://pymol.org/2/).

Rotamer analysis was performed using the Structure Comparison Tool within the Phenix Suite which assigns rotamer IDs according to the Ultimate Rotamer Library designation. 30 An AlphaFold2 rotamer was considered incorrect when its assigned Ultimate Rotamer Library ID failed to match those of any of the equivalent residues within the experimentally determined structures (all chains from PDBIDs: 7JTL, 7JX6). 

The phase problem

The detection of sub-units within the crystallographic asymmetric unit

Molecular replacement-Historical background

Molecular replacement then and now

Molecular replacement: Tricks and treats

Evaluating the usefulness of protein structure models for molecular replacement

Advances in rosetta structure prediction for difficult molecular-replacement problems

Rosetta structure prediction as a tool for solving difficult molecular replacement problems

Molecular replacement using structure predictions from databases

Ample: A cluster-and-truncate approach to solve the crystal structures of small proteins using rapidly computed ab initio models

Crystal structure of human chondroadherin: Solving a difficult molecular-replacement problem using de novo models

High precision in protein contact prediction using fully convolutional neural networks and minimal sequence features

Combining coevolution methods for accurate prediction of contacts and long range hydrogen bonding in proteins

Improved contact predictions using the recognition of protein like contact patterns

Accurate de novo prediction of protein contact map by ultra-deep learning model

Improved protein structure prediction using potentials from deep learning

Analysis of distance-based protein structure prediction by deep learning in casp13

Prediction of interresidue contacts with deepmetapsicov in casp13

A position-specific distance-dependent statistical potential for protein structure and functional study

Global fold determination from a small number of distance restraints

Protein structure prediction using multiple deep neural networks in the 13th critical assessment of protein structure prediction (casp13)

It will change everything': Deepmind's ai makes gigantic leap in solving protein structures

Structure of sars-cov-2 orf8, a rapidly evolving immune evasion protein

Phaser crystallographic software

Structural basis of dscam1 homodimerization: Insights into context constraint for protein recognition

Structure of natural killer receptor 2b4 bound to cd48 reveals basis for heterophilic recognition in signaling lymphocyte activation molecule family

High-resolution comparative modeling with rosettacm

Real-space refinement in phenix for cryo-em and crystallography

Features and development of coot

Molprobity's ultimate rotamer-library distributions for model validation

Crystallographic molecular replacement using an in silico-generated search model of SARS-CoV-2 ORF8