AlphaFold Models of Small Proteins Rival the Accuracy of Solution NMR Structures

Differential molecular interactions between iberiotoxin and human SLO3 and SLO1 potassium channels

CONTEXT

SLO1and SLO3 are similar voltage-gated K + channels. However, SLO3 expression is sperm specific and plays an important role in the hyperpolarization of the sperm membrane potential that is crucial for sperm fertilization. This makes SLO3 an excellent molecular target for the development of male contraceptives, and computational methods can facilitate structural insights for this drug development. Here, we evaluated the differential molecular interactions between the human SLO3 (hSLO3) and SLO1 (hSLO1) potassium channels and iberiotoxin (IbTX), a toxin that selectively blocks SLO channels. To do this, molecular docking and dynamics were implemented on the channel-toxin complexes to help elucidate atomistic details of their interaction and binding energy. Our analysis found that IbTX has a similar binding energy to both channels but interacts in a distinct manner with them. Particularly, Trp14 and Arg25 residues of IbTX diverges in their interaction with the residues Val283 and Asn260 residues of hSLO3 and the corresponding residues Tyr359 and Ala336 of hSLO1. Knowledge of key residues in the molecular interface of IbTX blockage can help guide and hasten non-hormonal contraceptive development. Our results encourage the use of toxins as scaffolds for specific SLO3 blockers.

METHODS

Atomistic molecular dynamics were implemented on the channel-toxin complexes. To generate the complexes, IbTX was docked to the channels using HADDOCK. CHARMM-GUI was used to generate simulation systems. GROMACS v2023.1 was used to run the simulations for 500 ns in an NPT ensemble at 297.26 K employing the CHARMM36 force field. Binding energy was evaluated by molecular mechanics generalized born surface area (MM/GBSA) with gmxMMPBGBSA.py.

Leveraging AlphaFold2 and residual dipolar couplings for side-chain methyl group assignment: A case study with S. cerevisiae Xrs2

Side-chain methyl group NMR spectroscopy provides invaluable insights into macromolecular structure, dynamics, and function, particularly for large biomolecular complexes. Accurate assignment of methyl group resonances in two-dimensional spectra is essential for structural and dynamics studies. Traditional methyl group assignment strategies rely on either transferring assignments from backbone resonance data or NOESY data and high-resolution experimental structures; however, these methods are often limited by molecular size or availability of structural information, respectively. Here, we describe the use of AlphaFold2 structural models as a basis for the manual, distance-based assignment of side-chain methyl group resonances in the folded domains of S. cerevisiae Xrs2. While AlphaFold2 models facilitated initial assignments for the methyl resonances, inaccuracies in the side-chain coordinates highlighted the need for improved structural models. By generating >500 ColabFold-derived models and filtering with methyl residual dipolar couplings (RDCs), we identified structural models with superior agreement to experimental data. These refined models enabled additional methyl group assignments while suggesting an iterative approach to simultaneously improve structure prediction and resonance assignment. Our findings outline a workflow that integrates machine learning-based structural predictions with experimental NMR data, offering a pathway for advancing methyl group assignment in systems lacking high-resolution experimental structures.

Deep learning and its applications in nuclear magnetic resonance spectroscopy

Nuclear Magnetic Resonance (NMR), as an advanced technology, has widespread applications in various fields like chemistry, biology, and medicine. However, issues such as long acquisition times for multidimensional spectra and low sensitivity limit the broader application of NMR. Traditional algorithms aim to address these issues but have limitations in speed and accuracy. Deep Learning (DL), a branch of Artificial Intelligence (AI) technology, has shown remarkable success in many fields including NMR. This paper presents an overview of the basics of DL and current applications of DL in NMR, highlights existing challenges, and suggests potential directions for improvement.

Hidden Structural States of Proteins Revealed by Conformer Selection with AlphaFold-NMR

We introduce AlphaFold-NMR, a novel approach to NMR structure determination that reveals previously undetected protein conformational states. Unlike conventional NMR methods that rely on NOE-derived spatial restraints, AlphaFold-NMR combines AI-driven conformational sampling with Bayesian scoring of realistic protein models against NOESY and chemical shift data. This method uncovers alternative conformational states of the enzyme Gaussia luciferase, involving large-scale changes in the lid, binding pockets, and other surface cavities. It also identifies similar yet distinct conformational states of the human tumor suppressor Cyclin-Dependent Kinase 2-Associated Protein 1. These studies demonstrate the potential of AI-based modeling with enhanced sampling to generate diverse structural models followed by conformer selection and validation with experimental data as an alternative to traditional restraint-satisfaction protocols for protein NMR structure determination. The AlphaFold-NMR framework enables discovery of conformational heterogeneity and cryptic pockets that conventional NMR analysis methods do not distinguish, providing new insights into protein structure-function relationships.

Artificial intelligence in drug development: reshaping the therapeutic landscape

Artificial intelligence (AI) is transforming medication research and development, giving clinicians new treatment options. Over the past 30 years, machine learning, deep learning, and neural networks have revolutionized drug design, target identification, and clinical trial predictions. AI has boosted pharmaceutical R&D (research and development) by identifying new therapeutic targets, improving chemical designs, and predicting complicated protein structures. Furthermore, generative AI is accelerating the development and re-engineering of medicinal molecules to cater to both common and rare diseases. Although, to date, no AI-generated medicinal drug has been FDA-approved, HLX-0201 for fragile X syndrome and new molecules for idiopathic pulmonary fibrosis have entered clinical trials. However, AI models are generally considered "black boxes," making their conclusions challenging to understand and limiting the potential due to a lack of model transparency and algorithmic bias. Despite these obstacles, AI-driven drug discovery has substantially reduced development times and costs, expediting the process and financial risks of bringing new medicines to market. In the future, AI is expected to continue to impact pharmaceutical innovation positively, making life-saving drug discoveries faster, more efficient, and more widespread.

Small Disulfide Proteins with Antifungal Impact: NMR Experimental Structures as Compared to Models of Alphafold Versions

In response to the growth of emerging resistance to conventional antifungal drugs, antifungal proteins (AFPs) of filamentous Ascomycetes origin have been discovered in recent years. Understanding the structure of AFPs is crucial for elucidating their antifungal mechanisms and developing new therapeutic agents. While nuclear magnetic resonance (NMR) has proven effective in determining the structures of small proteins, some AFP structures remain unresolved, necessitating the use of alternative prediction methods. Through bioinformatics analysis and heatmaps of amino acid sequence identity and similarity matrix, we categorized AFPs into three major classes and six subcategories, revealing structural and bioactivity differences. We employed AlphaFold (AF) to predict the 3D structures of six different AFPs, with predictions compared to NMR-derived structures. The results demonstrated a high degree of consistency between AF and NMR structures, with AF excelling in structural quality assessment and accurately capturing complex disulfide bond patterns. Both AF2 and AF3 models outperform the NMR model in overall structural quality and coherence, with AF3 showing the best performance. However, the limitations of AF should be considered, including its reduced accuracy in predicting multi-metal ion complexes, suboptimal performance in highly flexible or disordered regions, and its inability to account for multiple conformers, as it generates only a single dominant structure. Moreover, while AF3 accurately predicts all disulfide bond patterns, AF2 falls short in this regard. This study verifies the reliability of AF in the structural prediction of cysteine-rich AFPs while highlighting these constraints, offering important support for the rational design of new protein-based antifungal drugs.

Revealing Missing Protein-Ligand Interactions Using AlphaFold Predictions

Protein-ligand interactions represent an essential step to understand the bases of molecular recognition, an intense field of research in many scientific areas. Structural biology has played a central role in unveiling protein-ligand interactions, but current techniques are still not able to reliably describe the interactions of ligands with highly flexible regions. In this work, we explored the capacity of AlphaFold2 (AF2) to estimate the presence of interactions between ligands and residues belonging to disordered regions. As these interactions are missing in the crystallographic-derived structures, we called them "ghost interactions". Using a set of protein structures experimentally obtained after AF2 was trained, we found that the obtained models are good predictors of regions associated with order-disorder transitions. Additionally, we found that AF2 predicts residues making ghost interactions with ligands, which are mostly buried and show differential evolutionary conservation with the rest of the residues located in the flexible region. Our findings could fuel current areas of research that consider, given their biological relevance and their involvement in diseases, intrinsically disordered proteins as potentially valuable targets for drug development.

Predicting multiple conformations of ligand binding sites in proteins suggests that AlphaFold2 may remember too much

The goal of this paper is predicting the conformational distributions of ligand binding sites using the AlphaFold2 (AF2) protein structure prediction program with stochastic subsampling of the multiple sequence alignment (MSA). We explored the opening of cryptic ligand binding sites in 16 proteins, where the closed and open conformations define the expected extreme points of the conformational variation. Due to the many structures of these proteins in the Protein Data Bank (PDB), we were able to study whether the distribution of X-ray structures affects the distribution of AF2 models. We have found that AF2 generates both a cluster of open and a cluster of closed models for proteins that have comparable numbers of open and closed structures in the PDB and not too many other conformations. This was observed even with default MSA parameters, thus without further subsampling. In contrast, with the exception of a single protein, AF2 did not yield multiple clusters of conformations for proteins that had imbalanced numbers of open and closed structures in the PDB, or had substantial numbers of other structures. Subsampling improved the results only for a single protein, but very shallow MSA led to incorrect structures. The ability of generating both open and closed conformations for six out of the 16 proteins agrees with the success rates of similar studies reported in the literature. However, we showed that this partial success is due to AF2 "remembering" the conformational distributions in the PDB and that the approach fails to predict rarely seen conformations.

Coevolutionary Analysis of the Pfs47-P47Rec Complex: A Bioinformatics Approach

Background

The ability to predict and comprehend molecular interactions offers significant insights into the biological functions of proteins. The interaction between surface protein 47 of Plasmodium falciparum (Pfs47) and receptor of the protein 47 (P47Rec) has attracted increased attention due to their role in parasite evasion of the mosquito immune system and the concept of geographical coevolution between species. The aims of this study were as follows: to apply a bioinformatics approach to investigate the interaction between Pfs47 and P47Rec proteins and to identify the potential binding sites, protein orientations and receptor specificity sites concerning the geographical origins of the vectors and the parasite.

Methods

Public sequences of the pfs47 and p47rec genes were downloaded and subsequently filtered to predict functional and structural annotations of the Pfs47-P47Rec complex. Phylogenetic analyses of both proteins were carried out. In addition, the p47Rec gene was subjected to sequencing and subsequent analysis in 2 distinct Anopheles species collected in Honduras.

Results

The examination of motifs reveals a significant degree of conservation in pfs47, suggesting that Pfs47 might have undergone recent evolutionary development and adaptation. Structural models and docking analyses supported the theory of selectivity of Plasmodium falciparum strains towards their vectors in diverse geographical regions. A detailed description of the putative interaction between the Pfs47-P47Rec complex is shown.

Conclusions

The study identifies coevolutionary patterns between P47Rec and Pfs47 related to the speciation and geographic dispersion of Anopheles species and Plasmodium falciparum, with Pfs47 evolving more recently than P47Rec. This suggests a link between the parasite's adaptability and existing anopheline species across different regions. P47Rec likely has a cytoplasmic localization due to its lack of membrane attachment elements. However, these findings are based on simulations and require validation through methods like cryo-electron microscopy. A significant limitation is the scarcity of sequences in global databases, which restricts precise interaction modelling. Further research with diverse parasite isolates and anopheline species is recommended to enhance understanding of these proteins' structure and interaction.

Insights into AlphaFold's breakthrough in neurodegenerative diseases

Neurodegenerative diseases (ND) are disorders of the central nervous system (CNS) characterized by impairment in neurons' functions, and complete loss, leading to memory loss, and difficulty in learning, language, and movement processes. The most common among these NDs are Alzheimer's disease (AD) and Parkinson's disease (PD), although several other disorders also exist. These are frontotemporal dementia (FTD), amyotrophic lateral syndrome (ALS), Huntington's disease (HD), and others; the major pathological hallmark of NDs is the proteinopathies, either of amyloid-β (Aβ), tauopathies, or synucleinopathies. Aggregation of proteins that do not undergo normal configuration, either due to mutations or through some disturbance in cellular pathway contributes to the diseases. Artificial Intelligence (AI) and deep learning (DL) have proven to be successful in the diagnosis and treatment of various congenital diseases. DL approaches like AlphaFold (AF) are a major leap towards success in CNS disorders. This 3D protein geometry modeling algorithm developed by DeepMind has the potential to revolutionize biology. AF has the potential to predict 3D-protein confirmation at an accuracy level comparable to experimentally predicted one, with the additional advantage of precisely estimating protein interactions. This breakthrough will be beneficial to identify diseases' advancement and the disturbance of signaling pathways stimulating impaired functions of proteins. Though AlphaFold has solved a major problem in structural biology, it cannot predict membrane proteins-a beneficial approach for drug designing.

Integrative Modeling of Protein-Polypeptide Complexes by Bayesian Model Selection using AlphaFold and NMR Chemical Shift Perturbation Data

Protein-polypeptide interactions, including those involving intrinsically-disordered peptides and intrinsically-disordered regions of protein binding partners, are crucial for many biological functions. However, experimental structure determination of protein-peptide complexes can be challenging. Computational methods, while promising, generally require experimental data for validation and refinement. Here we present CSP_Rank, an integrated modeling approach to determine the structures of protein-peptide complexes. This method combines AlphaFold2 (AF2) enhanced sampling methods with a Bayesian conformational selection process based on experimental Nuclear Magnetic Resonance (NMR) Chemical Shift Perturbation (CSP) data and AF2 confidence metrics. Using a curated dataset of 108 protein-peptide complexes from the Biological Magnetic Resonance Data Bank (BMRB), we observe that while AF2 typically yields models with excellent consistency with experimental CSP data, applying enhanced sampling followed by data-guided conformational selection routinely results in ensembles of structures with improved agreement with NMR observables. For two systems, we cross-validate the CSP-selected models using independently acquired nuclear Overhauser effect (NOE) NMR data and demonstrate how CSP and NMR can be combined using our Bayesian framework for model selection. CSP_Rank is a novel method for integrative modeling of protein-peptide complexes and has broad implications for studies of protein-peptide interactions and aiding in understanding their biological functions.

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.

AlphaFold2 as a replacement for solution NMR structure determination of small proteins: Not so fast!

The determination of a protein's structure is often a first step towards the development of a mechanistic understanding of its function. Considerable advances in computational protein structure prediction have been made in recent years, with AlphaFold2 (AF2) emerging as the primary tool used by researchers for this purpose. While AF2 generally predicts accurate structures of folded proteins, we present here a case where AF2 incorrectly predicts the structure of a small, folded and compact protein with high confidence. This protein, pro-interleukin-18 (pro-IL-18), is the precursor of the cytokine IL-18. Interestingly, the structure of pro-IL-18 predicted by AF2 matches that of the mature cytokine, and not the corresponding experimentally determined structure of the pro-form of the protein. Thus, while computational structure prediction holds immense promise for addressing problems in protein biophysics, there is still a need for experimental structure determination, even in the context of small well-folded, globular proteins.

A predicted model-aided reconstruction algorithm for X-ray free-electron laser single-particle imaging

Ultra-intense, ultra-fast X-ray free-electron lasers (XFELs) enable the imaging of single protein molecules under ambient temperature and pressure. A crucial aspect of structure reconstruction involves determining the relative orientations of each diffraction pattern and recovering the missing phase information. In this paper, we introduce a predicted model-aided algorithm for orientation determination and phase retrieval, which has been tested on various simulated datasets and has shown significant improvements in the success rate, accuracy and efficiency of XFEL data reconstruction.

Advancing Molecular Dynamics: Toward Standardization, Integration, and Data Accessibility in Structural Biology

Molecular dynamics (MD) simulations have become a valuable tool in structural biology, offering insights into complex biological systems that are difficult to obtain through experimental techniques alone. The lack of available data sets and structures in most published computational work has limited other researchers' use of these models. In recent years, the emergence of online sharing platforms and MD database initiatives favor the deposition of ensembles and structures to accompany publications, favoring reuse of the data sets. However, the lack of uniform metadata collection, formats, and what data are deposited limits the impact and its use by different communities that are not necessarily experts in MD. This Perspective highlights the need for standardization and better resource sharing for processing and interpreting MD simulation results, akin to efforts in other areas of structural biology. As the field moves forward, we will see an increase in popularity and benefits of MD-based integrative approaches combining experimental data and simulations through probabilistic reasoning, but these too are limited by uniformity in experimental data availability and choices on how the data are modeled that are not trivial to decipher from papers. Other fields have addressed similar challenges comprehensively by establishing task forces with different degrees of success. The large scope and number of communities to represent the breadth of types of MD simulations complicates a parallel approach that would fit all. Thus, each group typically decides what data and which format to upload on servers like Zenodo. Uploading data with FAIR (findable, accessible, interoperable, reusable) principles in mind including optimal metadata collection will make the data more accessible and actionable by the community. Such a wealth of simulation data will foster method development and infrastructure advancements, thus propelling the field forward.

The 100-protein NMR spectra dataset: A resource for biomolecular NMR data analysis

Multidimensional NMR spectra are the basis for studying proteins by NMR spectroscopy and crucial for the development and evaluation of methods for biomolecular NMR data analysis. Nevertheless, in contrast to derived data such as chemical shift assignments in the BMRB and protein structures in the PDB databases, this primary data is in general not publicly archived. To change this unsatisfactory situation, we present a standardized set of solution NMR data comprising 1329 2-4-dimensional NMR spectra and associated reference (chemical shift assignments, structures) and derived (peak lists, restraints for structure calculation, etc.) annotations. With the 100-protein NMR spectra dataset that was originally compiled for the development of the ARTINA deep learning-based spectra analysis method, 100 protein structures can be reproduced from their original experimental data. The 100-protein NMR spectra dataset is expected to help the development of computational methods for NMR spectroscopy, in particular machine learning approaches, and enable consistent and objective comparisons of these methods.

How accurately can one predict drug binding modes using AlphaFold models?

Computational prediction of protein structure has been pursued intensely for decades, motivated largely by the goal of using structural models for drug discovery. Recently developed machine-learning methods such as AlphaFold 2 (AF2) have dramatically improved protein structure prediction, with reported accuracy approaching that of experimentally determined structures. To what extent do these advances translate to an ability to predict more accurately how drugs and drug candidates bind to their target proteins? Here, we carefully examine the utility of AF2 protein structure models for predicting binding poses of drug-like molecules at the largest class of drug targets, the G-protein-coupled receptors. We find that AF2 models capture binding pocket structures much more accurately than traditional homology models, with errors nearly as small as differences between structures of the same protein determined experimentally with different ligands bound. Strikingly, however, the accuracy of ligand-binding poses predicted by computational docking to AF2 models is not significantly higher than when docking to traditional homology models and is much lower than when docking to structures determined experimentally without these ligands bound. These results have important implications for all those who might use predicted protein structures for drug discovery.

Exposing the small protein load of bacterial life

The ever-growing repertoire of genomic techniques continues to expand our understanding of the true diversity and richness of prokaryotic genomes. Riboproteogenomics laid the foundation for dynamic studies of previously overlooked genomic elements. Most strikingly, bacterial genomes were revealed to harbor robust repertoires of small open reading frames (sORFs) encoding a diverse and broadly expressed range of small proteins, or sORF-encoded polypeptides (SEPs). In recent years, continuous efforts led to great improvements in the annotation and characterization of such proteins, yet many challenges remain to fully comprehend the pervasive nature of small proteins and their impact on bacterial biology. In this work, we review the recent developments in the dynamic field of bacterial genome reannotation, catalog the important biological roles carried out by small proteins and identify challenges obstructing the way to full understanding of these elusive proteins.

Blind assessment of monomeric AlphaFold2 protein structure models with experimental NMR data

Recent advances in molecular modeling of protein structures are changing the field of structural biology. AlphaFold-2 (AF2), an AI system developed by DeepMind, Inc., utilizes attention-based deep learning to predict models of protein structures with high accuracy relative to structures determined by X-ray crystallography and cryo-electron microscopy (cryoEM). Comparing AF2 models to structures determined using solution NMR data, both high similarities and distinct differences have been observed. Since AF2 was trained on X-ray crystal and cryoEM structures, we assessed how accurately AF2 can model small, monomeric, solution protein NMR structures which (i) were not used in the AF2 training data set, and (ii) did not have homologous structures in the Protein Data Bank at the time of AF2 training. We identified nine open-source protein NMR data sets for such "blind" targets, including chemical shift, raw NMR FID data, NOESY peak lists, and (for 1 case) 15N-1H residual dipolar coupling data. For these nine small (70-108 residues) monomeric proteins, we generated AF2 prediction models and assessed how well these models fit to these experimental NMR data, using several well-established NMR structure validation tools. In most of these cases, the AF2 models fit the NMR data nearly as well, or sometimes better than, the corresponding NMR structure models previously deposited in the Protein Data Bank. These results provide benchmark NMR data for assessing new NMR data analysis and protein structure prediction methods. They also document the potential for using AF2 as a guiding tool in protein NMR data analysis, and more generally for hypothesis generation in structural biology research.

Solution Structure Ensembles of the Open and Closed Forms of the ∼130 kDa Enzyme I via AlphaFold Modeling, Coarse Grained Simulations, and NMR

Large-scale interdomain rearrangements are essential to protein function, governing the activity of large enzymes and molecular machineries. Yet, obtaining an atomic-resolution understanding of how the relative domain positioning is affected by external stimuli is a hard task in modern structural biology. Here, we show that combining structural modeling by AlphaFold2 with coarse-grained molecular dynamics simulations and NMR residual dipolar coupling data is sufficient to characterize the spatial domain organization of bacterial enzyme I (EI), a ∼130 kDa multidomain oligomeric protein that undergoes large-scale conformational changes during its catalytic cycle. In particular, we solve conformational ensembles for EI at two different experimental temperatures and demonstrate that a lower temperature favors sampling of the catalytically competent closed state of the enzyme. These results suggest a role for conformational entropy in the activation of EI and demonstrate the ability of our protocol to detect and characterize the effect of external stimuli (such as mutations, ligand binding, and post-translational modifications) on the interdomain organization of multidomain proteins. We expect the ensemble refinement protocol described here to be easily transferrable to the investigation of the structure and dynamics of other uncharted multidomain systems and have assembled a Google Colab page (https://potoyangroup.github.io/Seq2Ensemble/) to facilitate implementation of the presented methodology elsewhere.

AlphaFold, allosteric, and orthosteric drug discovery: Ways forward

Drug discovery is arguably a highly challenging and significant interdisciplinary aim. The stunning success of the artificial intelligence-powered AlphaFold, whose latest version is buttressed by an innovative machine-learning approach that integrates physical and biological knowledge about protein structures, raised drug discovery hopes that unsurprisingly, have not come to bear. Even though accurate, the models are rigid, including the drug pockets. AlphaFold's mixed performance poses the question of how its power can be harnessed in drug discovery. Here we discuss possible ways of going forward wielding its strengths, while bearing in mind what AlphaFold can and cannot do. For kinases and receptors, an input enriched in active (ON) state models can better AlphaFold's chance of rational drug design success.

Recent Advances in NMR Protein Structure Prediction with ROSETTA

Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for studying the structure and dynamics of proteins in their native state. For high-resolution NMR structure determination, the collection of a rich restraint dataset is necessary. This can be difficult to achieve for proteins with high molecular weight or a complex architecture. Computational modeling techniques can complement sparse NMR datasets (<1 restraint per residue) with additional structural information to elucidate protein structures in these difficult cases. The Rosetta software for protein structure modeling and design is used by structural biologists for structure determination tasks in which limited experimental data is available. This review gives an overview of the computational protocols available in the Rosetta framework for modeling protein structures from NMR data. We explain the computational algorithms used for the integration of different NMR data types in Rosetta. We also highlight new developments, including modeling tools for data from paramagnetic NMR and hydrogen-deuterium exchange, as well as chemical shifts in CS-Rosetta. Furthermore, strategies are discussed to complement and improve structure predictions made by the current state-of-the-art AlphaFold2 program using NMR-guided Rosetta modeling.

AlphaFold2 and its applications in the fields of biology and medicine

AlphaFold2 (AF2) is an artificial intelligence (AI) system developed by DeepMind that can predict three-dimensional (3D) structures of proteins from amino acid sequences with atomic-level accuracy. Protein structure prediction is one of the most challenging problems in computational biology and chemistry, and has puzzled scientists for 50 years. The advent of AF2 presents an unprecedented progress in protein structure prediction and has attracted much attention. Subsequent release of structures of more than 200 million proteins predicted by AF2 further aroused great enthusiasm in the science community, especially in the fields of biology and medicine. AF2 is thought to have a significant impact on structural biology and research areas that need protein structure information, such as drug discovery, protein design, prediction of protein function, et al. Though the time is not long since AF2 was developed, there are already quite a few application studies of AF2 in the fields of biology and medicine, with many of them having preliminarily proved the potential of AF2. To better understand AF2 and promote its applications, we will in this article summarize the principle and system architecture of AF2 as well as the recipe of its success, and particularly focus on reviewing its applications in the fields of biology and medicine. Limitations of current AF2 prediction will also be discussed.

Ranking Peptide Binders by Affinity with AlphaFold

AlphaFold has revolutionized structural biology by predicting highly accurate structures of proteins and their complexes with peptides and other proteins. However, for protein-peptide systems, we are also interested in identifying the highest affinity binder among a set of candidate peptides. We present a novel competitive binding assay using AlphaFold to predict structures of the receptor in the presence of two peptides. For systems in which the individual structures of the peptides are well predicted, the assay captures the higher affinity binder in the bound state, and the other peptide in the unbound form with statistical significance. We test the application on six protein receptors for which we have experimental binding affinities to several peptides. We find that the assay is best suited for identifying medium to strong peptide binders that adopt stable secondary structures upon binding.

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.

Blind Assessment of Monomeric AlphaFold2 Protein Structure Models with Experimental NMR Data

Recent advances in molecular modeling of protein structures are changing the field of structural biology. AlphaFold-2 (AF2), an AI system developed by DeepMind, Inc., utilizes attention-based deep learning to predict models of protein structures with high accuracy relative to structures determined by X-ray crystallography and cryo-electron microscopy (cryoEM). Comparing AF2 models to structures determined using solution NMR data, both high similarities and distinct differences have been observed. Since AF2 was trained on X-ray crystal and cryoEM structures, we assessed how accurately AF2 can model small, monomeric, solution protein NMR structures which (i) were not used in the AF2 training data set, and (ii) did not have homologous structures in the Protein Data Bank at the time of AF2 training. We identified nine open source protein NMR data sets for such "blind" targets, including chemical shift, raw NMR FID data, NOESY peak lists, and (for 1 case) 15 N- 1 H residual dipolar coupling data. For these nine small (70 - 108 residues) monomeric proteins, we generated AF2 prediction models and assessed how well these models fit to these experimental NMR data, using several well-established NMR structure validation tools. In most of these cases, the AF2 models fit the NMR data nearly as well, or sometimes better than, the corresponding NMR structure models previously deposited in the Protein Data Bank. These results provide benchmark NMR data for assessing new NMR data analysis and protein structure prediction methods. They also document the potential for using AF2 as a guiding tool in protein NMR data analysis, and more generally for hypothesis generation in structural biology research.

Highlights

AF2 models assessed against NMR data for 9 monomeric proteins not used in training.AF2 models fit NMR data almost as well as the experimentally-determined structures. RPF-DP, PSVS , and PDBStat software provide structure quality and RDC assessment. RPF-DP analysis using AF2 models suggests multiple conformational states.