Computational models in the service of X-ray and cryo-electron microscopy structure determination

Structural studies of Parvoviridae capsid assembly and evolution: implications for novel AAV vector design

Adeno-associated virus (AAV) vectors have emerged as powerful tools in gene therapy, potentially treating various genetic disorders. Engineering the AAV capsids through computational methods enables the customization of these vectors to enhance their effectiveness and safety. This engineering allows for the development of gene therapies that are not only more efficient but also personalized to unique genetic profiles. When developing, it is essential to understand the structural biology and the vast techniques used to guide vector designs. This review covers the fundamental biology of the Parvoviridae capsids, focusing on modern structural study techniques, including (a) Cryo-electron microscopy and X-ray Crystallography studies and (b) Comparative analysis of capsid structures across different Parvoviridae species. Along with the structure and evolution of the Parvoviridae capsids, computational methods have provided significant insights into the design of novel AAV vector techniques, which include (a) Structure-guided design of AAV capsids with improved properties, (b) Directed Evolution of AAV capsids for specific applications, and (c) Computational prediction of AAV capsid-receptor interactions. Further discussion addressed the ongoing challenges in the AAV vector design and proposed future directions for exploring enhanced computational tools, such as artificial intelligence/machine learning and deep learning.

Relationship between structural properties and biological activity of (-)-menthol and some menthyl esters

Menthol is a naturally occurring cyclic terpene alcohol and is the major component of peppermint and corn mint essential oils extracted from Mentha piperita L. and Mentha arvensis L.. Menthol and its derivatives are widely used in pharmaceutical, cosmetic and food industries. Among its eight isomers, (-)-menthol is the most effective one in terms of refreshing effect. While the invigorating property of (-)-menthol is generally known, this claim is based on a substantial amount of literature and experience. (-)-Menthol has consistently been reported to possess better cooling and refreshing qualities in comparison to its isomers, making it the preferred choice in a broad range of applications such as personal care products, pharmaceuticals and food additives. Additionally, the (-)-menthol molecular structure allows it to have a tighter fitting with the thermoreceptors in the skin and mucous membranes, and thus to provide a more intense cooling feeling. Thus, although others have similar properties to a degree, (-)-menthol is the best compared to all in its refreshing capacity. This study focuses on menthol and some of its esters, viz. menthyl acetate, propionate, butyrate, valerate and hexanoate, with the purpose of establish a connection between structural, electrostatic and electronic characteristics and biological effects. The mostly favoured interactions of the esters with biotargets were investigated at a molecular level, offering a plausible foundation for their bioactivity elucidation. This study is conducted at a quantum mechanical and molecular docking level. The results may be of possible usefulness in areas of applications, such as pharmacological research and drug.

Deep learning structural insights into heterotrimeric alternatively spliced P2X7 receptors

P2X7 receptors (P2X7Rs) are membrane-bound ATP-gated ion channels that are composed of three subunits. Different subunit structures may be expressed due to alternative splicing of the P2RX7 gene, altering the receptor's function when combined with the wild-type P2X7A subunits. In this study, the application of the deep-learning method, AlphaFold2-Multimer (AF2M), for the generation of trimeric P2X7Rs was validated by comparing an AF2M-generated rat wild-type P2X7A receptor with a structure determined by cryogenic electron microscopy (cryo-EM) (Protein Data Bank Identification: 6U9V). The results suggested AF2M could firstly, accurately predict the structures of P2X7Rs and secondly, accurately identify the highest quality model through the ranking system. Subsequently, AF2M was used to generate models of heterotrimeric alternatively spliced P2X7Rs consisting of one or two wild-type P2X7A subunits in combination with one or two P2X7B, P2X7E, P2X7J, and P2X7L splice variant subunits. The top-ranking models were deemed valid based on AF2M's confidence measures, stability in molecular dynamics simulations, and consistent flexibility of the conserved regions between the models. The structure of the heterotrimeric receptors, which were missing key residues in the ATP binding sites and carboxyl terminal domains (CTDs) compared to the wild-type receptor, help to explain their observed functions. Overall, the models produced in this study (available as supplementary material) unlock the possibility of structure-based studies into the heterotrimeric P2X7Rs.

The power and pitfalls of AlphaFold2 for structure prediction beyond rigid globular proteins

Artificial intelligence-driven advances in protein structure prediction in recent years have raised the question: has the protein structure-prediction problem been solved? Here, with a focus on nonglobular proteins, we highlight the many strengths and potential weaknesses of DeepMind's AlphaFold2 in the context of its biological and therapeutic applications. We summarize the subtleties associated with evaluation of AlphaFold2 model quality and reliability using the predicted local distance difference test (pLDDT) and predicted aligned error (PAE) values. We highlight various classes of proteins that AlphaFold2 can be applied to and the caveats involved. Concrete examples of how AlphaFold2 models can be integrated with experimental data in the form of small-angle X-ray scattering (SAXS), solution NMR, cryo-electron microscopy (cryo-EM) and X-ray diffraction are discussed. Finally, we highlight the need to move beyond structure prediction of rigid, static structural snapshots toward conformational ensembles and alternate biologically relevant states. The overarching theme is that careful consideration is due when using AlphaFold2-generated models to generate testable hypotheses and structural models, rather than treating predicted models as de facto ground truth structures.

Solving protein structures by combining structure prediction, molecular replacement and direct-methods-aided model completion

Highly accurate protein structure prediction can generate accurate models of protein and protein-protein complexes in X-ray crystallography. However, the question of how to make more effective use of predicted models for completing structure analysis, and which strategies should be employed for the more challenging cases such as multi-helical structures, multimeric structures and extremely large structures, both in the model preparation and in the completion steps, remains open for discussion. In this paper, a new strategy is proposed based on the framework of direct methods and dual-space iteration, which can greatly simplify the pre-processing steps of predicted models both in normal and in challenging cases. Following this strategy, full-length models or the conservative structural domains could be used directly as the starting model, and the phase error and the model bias between the starting model and the real structure would be modified in the direct-methods-based dual-space iteration. Many challenging cases (from CASP14) have been tested for the general applicability of this constructive strategy, and almost complete models have been generated with reasonable statistics. The hybrid strategy therefore provides a meaningful scheme for X-ray structure determination using a predicted model as the starting point.

An extensive disulfide bond network prevents tail contraction in Agrobacterium tumefaciens phage Milano

A contractile sheath and rigid tube assembly is a widespread apparatus used by bacteriophages, tailocins, and the bacterial type VI secretion system to penetrate cell membranes. In this mechanism, contraction of an external sheath powers the motion of an inner tube through the membrane. The structure, energetics, and mechanism of the machinery imply rigidity and straightness. The contractile tail of Agrobacterium tumefaciens bacteriophage Milano is flexible and bent to varying degrees, which sets it apart from other contractile tail-like systems. Here, we report structures of the Milano tail including the sheath-tube complex, baseplate, and putative receptor-binding proteins. The flexible-to-rigid transformation of the Milano tail upon contraction can be explained by unique electrostatic properties of the tail tube and sheath. All components of the Milano tail, including sheath subunits, are crosslinked by disulfides, some of which must be reduced for contraction to occur. The putative receptor-binding complex of Milano contains a tailspike, a tail fiber, and at least two small proteins that form a garland around the distal ends of the tailspikes and tail fibers. Despite being flagellotropic, Milano lacks thread-like tail filaments that can wrap around the flagellum, and is thus likely to employ a different binding mechanism.

DEMO-EM2: assembling protein complex structures from cryo-EM maps through intertwined chain and domain fitting

The breakthrough in cryo-electron microscopy (cryo-EM) technology has led to an increasing number of density maps of biological macromolecules. However, constructing accurate protein complex atomic structures from cryo-EM maps remains a challenge. In this study, we extend our previously developed DEMO-EM to present DEMO-EM2, an automated method for constructing protein complex models from cryo-EM maps through an iterative assembly procedure intertwining chain- and domain-level matching and fitting for predicted chain models. The method was carefully evaluated on 27 cryo-electron tomography (cryo-ET) maps and 16 single-particle EM maps, where DEMO-EM2 models achieved an average TM-score of 0.92, outperforming those of state-of-the-art methods. The results demonstrate an efficient method that enables the rapid and reliable solution of challenging cryo-EM structure modeling problems.

DeepMainmast: integrated protocol of protein structure modeling for cryo-EM with deep learning and structure prediction

Three-dimensional structure modeling from maps is an indispensable step for studying proteins and their complexes with cryogenic electron microscopy. Although the resolution of determined cryogenic electron microscopy maps has generally improved, there are still many cases where tracing protein main chains is difficult, even in maps determined at a near-atomic resolution. Here we developed a protein structure modeling method, DeepMainmast, which employs deep learning to capture the local map features of amino acids and atoms to assist main-chain tracing. Moreover, we integrated AlphaFold2 with the de novo density tracing protocol to combine their complementary strengths and achieved even higher accuracy than each method alone. Additionally, the protocol is able to accurately assign the chain identity to the structure models of homo-multimers, which is not a trivial task for existing methods.

Integrated Protocol of Protein Structure Modeling for Cryo-EM with Deep Learning and Structure Prediction

Structure modeling from maps is an indispensable step for studying proteins and their complexes with cryogenic electron microscopy (cryo-EM). Although the resolution of determined cryo-EM maps has generally improved, there are still many cases where tracing protein main-chains is difficult, even in maps determined at a near atomic resolution. Here, we have developed a protein structure modeling method, called DeepMainmast, which employs deep learning to capture the local map features of amino acids and atoms to assist main-chain tracing. Moreover, since Alphafold2 demonstrates high accuracy in protein structure prediction, we have integrated complementary strengths of de novo density tracing using deep learning with Alphafold2's structure modeling to achieve even higher accuracy than each method alone. Additionally, the protocol is able to accurately assign chain identity to the structure models of homo-multimers.

Artificial intelligence and the analysis of cryo-EM data provide structural insight into the molecular mechanisms underlying LN-lamininopathies

Laminins (Lm) are major components of basement membranes (BM), which polymerize to form a planar lattice on cell surface. Genetic alternations of Lm affect their oligomerization patterns and lead to failures in BM assembly manifesting in a group of human disorders collectively defined as Lm N-terminal domain lamininopathies (LN-lamininopathies). We have employed a recently determined cryo-EM structure of the Lm polymer node, the basic repeating unit of the Lm lattice, along with structure prediction and modeling to systematically analyze structures of twenty-three pathogenic Lm polymer nodes implicated in human disease. Our analysis provides the detailed mechanistic explanation how Lm mutations lead to failures in Lm polymerization underlining LN-lamininopathies. We propose the new categorization scheme of LN-lamininopathies based on the insight gained from the structural analysis. Our results can help to facilitate rational drug design aiming in the treatment of Lm deficiencies.

Protein Structure Prediction: Challenges, Advances, and the Shift of Research Paradigms

Protein structure prediction is an interdisciplinary research topic that has attracted researchers from multiple fields, including biochemistry, medicine, physics, mathematics, and computer science. These researchers adopt various research paradigms to attack the same structure prediction problem: biochemists and physicists attempt to reveal the principles governing protein folding; mathematicians, especially statisticians, usually start from assuming a probability distribution of protein structures given a target sequence and then find the most likely structure, while computer scientists formulate protein structure prediction as an optimization problem - finding the structural conformation with the lowest energy or minimizing the difference between predicted structure and native structure. These research paradigms fall into the two statistical modeling cultures proposed by Leo Breiman, namely, data modeling and algorithmic modeling. Recently, we have also witnessed the great success of deep learning in protein structure prediction. In this review, we present a survey of the efforts for protein structure prediction. We compare the research paradigms adopted by researchers from different fields, with an emphasis on the shift of research paradigms in the era of deep learning. In short, the algorithmic modeling techniques, especially deep neural networks, have considerably improved the accuracy of protein structure prediction; however, theories interpreting the neural networks and knowledge on protein folding are still highly desired.

Just Keep Rolling?-An Encompassing Review towards Accelerated Vaccine Product Life Cycles

In light of the recent pandemic, several COVID-19 vaccines were developed, tested and approved in a very short time, a process that otherwise takes many years. Above all, these efforts have also unmistakably revealed the capacity limits and potential for improvement in vaccine production. This review aims to emphasize recent approaches for the targeted rapid adaptation and production of vaccines from an interdisciplinary, multifaceted perspective. Using research from the literature, stakeholder analysis and a value proposition canvas, we reviewed technological innovations on the pharmacological level, formulation, validation and resilient vaccine production to supply bottlenecks and logistic networks. We identified four main drivers to accelerate the vaccine product life cycle: computerized candidate screening, modular production, digitized quality management and a resilient business model with corresponding transparent supply chains. In summary, the results presented here can serve as a guide and implementation tool for flexible, scalable vaccine production to swiftly respond to pandemic situations in the future.

Critical Assessment of Methods for Predicting the 3D Structure of Proteins and Protein Complexes

Advances in a scientific discipline are often measured by small, incremental steps. In this review, we report on two intertwined disciplines in the protein structure prediction field, modeling of single chains and modeling of complexes, that have over decades emulated this pattern, as monitored by the community-wide blind prediction experiments CASP and CAPRI. However, over the past few years, dramatic advances were observed for the accurate prediction of single protein chains, driven by a surge of deep learning methodologies entering the prediction field. We review the mainscientific developments that enabled these recent breakthroughs and feature the important role of blind prediction experiments in building up and nurturing the structure prediction field. We discuss how the new wave of artificial intelligence-based methods is impacting the fields of computational and experimental structural biology and highlight areas in which deep learning methods are likely to lead to future developments, provided that major challenges are overcome.

Structural and enzymatic characterization of the sialidase SiaPG from Porphyromonas gingivalis

The sialidases, which catalyze the hydrolysis of sialic acid from extracellular glycoconjugates, are a group of major virulence factors in various pathogenic bacteria. In Porphyromonas gingivalis, which causes human periodontal disease, sialidase contributes to bacterial pathogenesis via promoting the formation of biofilms and capsules, reducing the ability for macrophage clearance, and providing nutrients for bacterial colonization. Here, the crystal structure of the P. gingivalis sialidase SiaPG is reported at 2.1 Å resolution, revealing an N-terminal carbohydrate-binding domain followed by a canonical C-terminal catalytic domain. Simulation of the product sialic acid in the active-site pocket together with functional analysis enables clear identification of the key residues that are required for substrate binding and catalysis. Moreover, structural comparison with other sialidases reveals distinct features of the active-site pocket which might confer substrate specificity. These findings provide the structural basis for the further design and optimization of effective inhibitors to target SiaPG to fight against P. gingivalis-derived oral diseases.

Function of the bacteriophage P2 baseplate central spike Apex domain in the infection process

The contractile tail of bacteriophage P2 functions to drive the tail tube across the outer membrane of its host bacterium, a prerequisite event for subsequent translocation of phage genomic DNA into the host cell. The tube is equipped with a spike-shaped protein (product of P2 gene V , gpV or Spike) that contains a membrane-attacking Apex domain carrying a centrally positioned Fe ion. The ion is enclosed in a histidine cage that is formed by three symmetry-related copies of a conserved HxH (histidine, any residue, histidine) sequence motif. Here, we used solution biophysics and X-ray crystallography to characterize the structure and properties of Spike mutants in which the Apex domain was either deleted or its histidine cage was either destroyed or replaced with a hydrophobic core. We found that the Apex domain is not required for the folding of full-length gpV or its middle intertwined β-helical domain. Furthermore, despite its high conservation, the Apex domain is dispensable for infection in laboratory conditions. Collectively, our results show that the diameter of the Spike but not the nature of its Apex domain determines the efficiency of infection, which further strengthens the earlier hypothesis of a drill bit-like function of the Spike in host envelope disruption.

Exploring AlphaFold2's Performance on Predicting Amino Acid Side-Chain Conformations and Its Utility in Crystal Structure Determination of B318L Protein

Recent technological breakthroughs in machine-learning-based AlphaFold2 (AF2) are pushing the prediction accuracy of protein structures to an unprecedented level that is on par with experimental structural quality. Despite its outstanding structural modeling capability, further experimental validations and performance assessments of AF2 predictions are still required, thus necessitating the development of integrative structural biology in synergy with both computational and experimental methods. Focusing on the B318L protein that plays an essential role in the African swine fever virus (ASFV) for viral replication, we experimentally demonstrate the high quality of the AF2 predicted model and its practical utility in crystal structural determination. Structural alignment implies that the AF2 model shares nearly the same atomic arrangement as the B318L crystal structure except for some flexible and disordered regions. More importantly, side-chain-based analysis at the individual residue level reveals that AF2's performance is likely dependent on the specific amino acid type and that hydrophobic residues tend to be more accurately predicted by AF2 than hydrophilic residues. Quantitative per-residue RMSD comparisons and further molecular replacement trials suggest that AF2 has a large potential to outperform other computational modeling methods in terms of structural determination. Additionally, it is numerically confirmed that the AF2 model is accurate enough so that it may well potentially withstand experimental data quality to a large extent for structural determination. Finally, an overall structural analysis and molecular docking simulation of the B318L protein are performed. Taken together, our study not only provides new insights into AF2's performance in predicting side-chain conformations but also sheds light upon the significance of AF2 in promoting crystal structural determination, especially when the experimental data quality of the protein crystal is poor.

De novo modelling of HEV replication polyprotein: Five-domain breakdown and involvement of flexibility in functional regulation

Hepatitis E virus (HEV), a major cause of acute viral hepatitis, is a single-stranded, positive-sense RNA virus. As such, it encodes a 1700-residue replication polyprotein pORF1 that directs synthesis of new viral RNA in infected cells. Here we report extensive modeling with AlphaFold2 of the full-length pORF1, and its production by in vitro translation. From this, we give a detailed update on the breakdown into domains of HEV pORF1. We also provide evidence that pORF1's N-terminal domain is likely to oligomerize to form a dodecameric pore, homologously to what has been described for Chikungunya virus. Beyond providing accurate folds for its five domains, our work highlights that there is no canonical protease encoded in pORF1 and that flexibility in several functionally important regions rather than proteolytic processing may serve to regulate HEV RNA synthesis.

The impact of AlphaFold2 on experimental structure solution

AlphaFold2 is a machine-learning based program that predicts a protein structure based on the amino acid sequence. In this article, we report on the current usages of this new tool and give examples from our work in the Coronavirus Structural Task Force. With its unprecedented accuracy, it can be utilized for the design of expression constructs, de novo protein design and the interpretation of Cryo-EM data with an atomic model. However, these methods are limited by their training data and are of limited use to predict conformational variability and fold flexibility; they also lack co-factors, post-translational modifications and multimeric complexes with oligonucleotides. They also are not always perfect in terms of chemical geometry. Nevertheless, machine learning-based fold prediction is a game changer for structural bioinformatics and experimentalists alike, with exciting developments ahead.

Present Impact of AlphaFold2 Revolution on Structural Biology, and an Illustration With the Structure Prediction of the Bacteriophage J-1 Host Adhesion Device

In 2021, the release of AlphaFold2 - the DeepMind's machine-learning protein structure prediction program - revolutionized structural biology. Results of the CASP14 contest were an immense surprise as AlphaFold2 successfully predicted 3D structures of nearly all submitted protein sequences. The AlphaFold2 craze has rapidly spread the life science community since structural biologists as well as untrained biologists have now the possibility to obtain high-confidence protein structures. This revolution is opening new avenues to address challenging biological questions. Moreover, AlphaFold2 is imposing itself as an essential step of any structural biology project, and requires us to revisit our structural biology workflows. On one hand, AlphaFold2 synergizes with experimental methods including X-ray crystallography and cryo-electron microscopy. On the other hand, it is, to date, the only method enabling structural analyses of large and flexible assemblies resistant to experimental approaches. We illustrate this valuable application of AlphaFold2 with the structure prediction of the whole host adhesion device from the Lactobacillus casei bacteriophage J-1. With the ongoing improvement of AlphaFold2 algorithms and notebooks, there is no doubt that AlphaFold2-driven biological stories will increasingly be reported, which questions the future directions of experimental structural biology.

A database of calculated solution parameters for the AlphaFold predicted protein structures

Recent spectacular advances by AI programs in 3D structure predictions from protein sequences have revolutionized the field in terms of accuracy and speed. The resulting “folding frenzy” has already produced predicted protein structure databases for the entire human and other organisms’ proteomes. However, rapidly ascertaining a predicted structure’s reliability based on measured properties in solution should be considered. Shape-sensitive hydrodynamic parameters such as the diffusion and sedimentation coefficients (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D_{t(20,w)}^{0}}$$\end{document}Dt(20,w)0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s_{{\left( {{20},w} \right)}}^{{0}} }$$\end{document}s20,w0) and the intrinsic viscosity ([η]) can provide a rapid assessment of the overall structure likeliness, and SAXS would yield the structure-related pair-wise distance distribution function p(r) vs. r. Using the extensively validated UltraScan SOlution MOdeler (US-SOMO) suite, a database was implemented calculating from AlphaFold structures the corresponding \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D_{t(20,w)}^{0}}$$\end{document}Dt(20,w)0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s_{{\left( {{20},w} \right)}}^{{0}} }$$\end{document}s20,w0, [η], p(r) vs. r, and other parameters. Circular dichroism spectra were computed using the SESCA program. Some of AlphaFold’s drawbacks were mitigated, such as generating whenever possible a protein’s mature form. Others, like the AlphaFold direct applicability to single-chain structures only, the absence of prosthetic groups, or flexibility issues, are discussed. Overall, this implementation of the US-SOMO-AF database should already aid in rapidly evaluating the consistency in solution of a relevant portion of AlphaFold predicted protein structures.

MrParse: finding homologues in the PDB and the EBI AlphaFold database for molecular replacement and more

Crystallographers have an array of search-model options for structure solution by molecular replacement (MR). The well established options of homologous experimental structures and regular secondary-structure elements or motifs are increasingly supplemented by computational modelling. Such modelling may be carried out locally or may use pre-calculated predictions retrieved from databases such as the EBI AlphaFold database. MrParse is a new pipeline to help to streamline the decision process in MR by consolidating bioinformatic predictions in one place. When reflection data are provided, MrParse can rank any experimental homologues found using eLLG, which indicates the likelihood that a given search model will work in MR. Inbuilt displays of predicted secondary structure, coiled-coil and transmembrane regions further inform the choice of MR protocol. MrParse can also identify and rank homologues in the EBI AlphaFold database, a function that will also interest other structural biologists and bioinformaticians.

Real-time structure search and structure classification for AlphaFold protein models

Last year saw a breakthrough in protein structure prediction, where the AlphaFold2 method showed a substantial improvement in the modeling accuracy. Following the software release of AlphaFold2, predicted structures by AlphaFold2 for proteins in 21 species were made publicly available via the AlphaFold Database. Here, to facilitate structural analysis and application of AlphaFold2 models, we provide the infrastructure, 3D-AF-Surfer, which allows real-time structure-based search for the AlphaFold2 models. In 3D-AF-Surfer, structures are represented with 3D Zernike descriptors (3DZD), which is a rotationally invariant, mathematical representation of 3D shapes. We developed a neural network that takes 3DZDs of proteins as input and retrieves proteins of the same fold more accurately than direct comparison of 3DZDs. Using 3D-AF-Surfer, we report structure classifications of AlphaFold2 models and discuss the correlation between confidence levels of AlphaFold2 models and intrinsic disordered regions.

The X-ray crystallography phase problem solved thanks to AlphaFold and RoseTTAFold models: a case-study report

The breakthrough recently made in protein structure prediction by deep-learning programs such as AlphaFold and RoseTTAFold will certainly revolutionize biology over the coming decades. The scientific community is only starting to appreciate the various applications, benefits and limitations of these protein models. Yet, after the first thrills due to this revolution, it is important to evaluate the impact of the proposed models and their overall quality to avoid the misinterpretation or overinterpretation of these models by biologists. One of the first applications of these models is in solving the `phase problem' encountered in X-ray crystallography in calculating electron-density maps from diffraction data. Indeed, the most frequently used technique to derive electron-density maps is molecular replacement. As this technique relies on knowledge of the structure of a protein that shares strong structural similarity with the studied protein, the availability of high-accuracy models is then definitely critical for successful structure solution. After the collection of a 2.45 Å resolution data set, we struggled for two years in trying to solve the crystal structure of a protein involved in the nonsense-mediated mRNA decay pathway, an mRNA quality-control pathway dedicated to the elimination of eukaryotic mRNAs harboring premature stop codons. We used different methods (isomorphous replacement, anomalous diffraction and molecular replacement) to determine this structure, but all failed until we straightforwardly succeeded thanks to both AlphaFold and RoseTTAFold models. Here, we describe how these new models helped us to solve this structure and conclude that in our case the AlphaFold model largely outcompetes the other models. We also discuss the importance of search-model generation for successful molecular replacement.

Application of Homology Modeling by Enhanced Profile-Profile Alignment and Flexible-Fitting Simulation to Cryo-EM Based Structure Determination

Application of cryo-electron microscopy (cryo-EM) is crucially important for ascertaining the atomic structure of large biomolecules such as ribosomes and protein complexes in membranes. Advances in cryo-EM technology and software have made it possible to obtain data with near-atomic resolution, but the method is still often capable of producing only a density map with up to medium resolution, either partially or entirely. Therefore, bridging the gap separating the density map and the atomic model is necessary. Herein, we propose a methodology for constructing atomic structure models based on cryo-EM maps with low-to-medium resolution. The method is a combination of sensitive and accurate homology modeling using our profile-profile alignment method with a flexible-fitting method using molecular dynamics simulation. As described herein, this study used benchmark applications to evaluate the model constructions of human two-pore channel 2 (one target protein in CASP13 with its structure determined using cryo-EM data) and the overall structure of Enterococcus hirae V-ATPase complex.

Implications of AlphaFold2 for crystallographic phasing by molecular replacement

The AlphaFold2 results in the 14th edition of Critical Assessment of Structure Prediction (CASP14) showed that accurate (low root-mean-square deviation) in silico models of protein structure domains are on the horizon, whether or not the protein is related to known structures through high-coverage sequence similarity. As highly accurate models become available, generated by harnessing the power of correlated mutations and deep learning, one of the aspects of structural biology to be impacted will be methods of phasing in crystallography. Here, the data from CASP14 are used to explore the prospects for changes in phasing methods, and in particular to explore the prospects for molecular-replacement phasing using in silico models.

Critical assessment of methods of protein structure prediction (CASP)-Round XIV

Critical assessment of structure prediction (CASP) is a community experiment to advance methods of computing three-dimensional protein structure from amino acid sequence. Core components are rigorous blind testing of methods and evaluation of the results by independent assessors. In the most recent experiment (CASP14), deep-learning methods from one research group consistently delivered computed structures rivaling the corresponding experimental ones in accuracy. In this sense, the results represent a solution to the classical protein-folding problem, at least for single proteins. The models have already been shown to be capable of providing solutions for problematic crystal structures, and there are broad implications for the rest of structural biology. Other research groups also substantially improved performance. Here, we describe these results and outline some of the many implications. Other related areas of CASP, including modeling of protein complexes, structure refinement, estimation of model accuracy, and prediction of inter-residue contacts and distances, are also described.

Cryo-EM targets in CASP14

Structures of seven CASP14 targets were determined using cryo-electron microscopy (cryo-EM) technique with resolution between 2.1 and 3.8 Å. We provide an evaluation of the submitted models versus the experimental data (cryo-EM density maps) and experimental reference structures built into the maps. The accuracy of models is measured in terms of coordinate-to-density and coordinate-to-coordinate fit. A-posteriori refinement of the most accurate models in their corresponding cryo-EM density resulted in structures that are close to the reference structure, including some regions with better fit to the density. Regions that were found to be less "refineable" correlate well with regions of high diversity between the CASP models and low goodness-of-fit to density in the reference structure.

Target highlights in CASP14: Analysis of models by structure providers

The biological and functional significance of selected Critical Assessment of Techniques for Protein Structure Prediction 14 (CASP14) targets are described by the authors of the structures. The authors highlight the most relevant features of the target proteins and discuss how well these features were reproduced in the respective submitted predictions. The overall ability to predict three-dimensional structures of proteins has improved remarkably in CASP14, and many difficult targets were modeled with impressive accuracy. For the first time in the history of CASP, the experimentalists not only highlighted that computational models can accurately reproduce the most critical structural features observed in their targets, but also envisaged that models could serve as a guidance for further studies of biologically-relevant properties of proteins.

Computational models in the service of X-ray and cryo-electron microscopy structure determination

Critical assessment of structure prediction (CASP) conducts community experiments to determine the state of the art in computing protein structure from amino acid sequence. The process relies on the experimental community providing information about not yet public or about to be solved structures, for use as targets. For some targets, the experimental structure is not solved in time for use in CASP. Calculated structure accuracy improved dramatically in this round, implying that models should now be much more useful for resolving many sorts of experimental difficulties. To test this, selected models for seven unsolved targets were provided to the experimental groups. These models were from the AlphaFold2 group, who overall submitted the most accurate predictions in CASP14. Four targets were solved with the aid of the models, and, additionally, the structure of an already solved target was improved. An a posteriori analysis showed that, in some cases, models from other groups would also be effective. This paper provides accounts of the successful application of models to structure determination, including molecular replacement for X-ray crystallography, backbone tracing and sequence positioning in a cryo-electron microscopy structure, and correction of local features. The results suggest that, in future, there will be greatly increased synergy between computational and experimental approaches to structure determination.

Assessing the utility of CASP14 models for molecular replacement

The assessment of CASP models for utility in molecular replacement is a measure of their use in a valuable real‐world application. In CASP7, the metric for molecular replacement assessment involved full likelihood‐based molecular replacement searches; however, this restricted the assessable targets to crystal structures with only one copy of the target in the asymmetric unit, and to those where the search found the correct pose. In CASP10, full molecular replacement searches were replaced by likelihood‐based rigid‐body refinement of models superimposed on the target using the LGA algorithm, with the metric being the refined log‐likelihood‐gain (LLG) score. This enabled multi‐copy targets and very poor models to be evaluated, but a significant further issue remained: the requirement of diffraction data for assessment. We introduce here the relative‐expected‐LLG (reLLG), which is independent of diffraction data. This reLLG is also independent of any crystal form, and can be calculated regardless of the source of the target, be it X‐ray, NMR or cryo‐EM. We calibrate the reLLG against the LLG for targets in CASP14, showing that it is a robust measure of both model and group ranking. Like the LLG, the reLLG shows that accurate coordinate error estimates add substantial value to predicted models. We find that refinement by CASP groups can often convert an inadequate initial model into a successful MR search model. Consistent with findings from others, we show that the AlphaFold2 models are sufficiently good, and reliably so, to surpass other current model generation strategies for attempting molecular replacement phasing.