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

Recent advances in molecular modeling using deep learning can revolutionize our understanding of dynamic protein structures. NMR is particularly well-suited for determining dynamic features of biomolecular structures. The conventional process for determining biomolecular structures from experimental NMR data involves its representation as conformation-dependent restraints, followed by generation of structural models guided by these spatial restraints. Here we describe an alternative approach: generating a distribution of realistic protein conformational models using artificial intelligence-(AI-) based methods and then selecting the sets of conformers that best explain the experimental data. We applied this conformational selection approach to redetermine the solution NMR structure of the enzyme Gaussia luciferase. First, we generated a diverse set of conformer models using AlphaFold2 (AF2) with an enhanced sampling protocol. The models that best-fit NOESY and chemical shift data were then selected with a Bayesian scoring metric. The resulting models include features of both the published NMR structure and the standard AF2 model generated without enhanced sampling. This "AlphaFold-NMR" protocol also generated an alternative "open" conformational state that fits nearly as well to the overall NMR data but accounts for some NOESY data that is not consistent with first "closed" conformational state; while other NOESY data consistent with this second state are not consistent with the first conformational state. The structure of this "open" structural state differs from that of the "closed" state primarily by the position of a thumb-shaped loop between α-helices H5 and H6, revealing a cryptic surface pocket. These alternative conformational states of Gluc are supported by "double recall" analysis of NOESY data and AF2 models. Additional structural states are also indicated by backbone chemical shift data indicating partially-disordered conformations for the C-terminal segment. Considered as a multistate ensemble, these multiple states of Gluc together fit the NOESY and chemical shift data better than the "restraint-based" NMR structure and provide novel insights into its structure-dynamic-function relationships. This study demonstrates the potential of AI-based modeling with enhanced sampling to generate conformational ensembles followed by conformer selection with experimental data as an alternative to conventional restraint satisfaction protocols for protein NMR structure determination.

A Computational Pipeline for Accurate Prioritization of Protein-Protein Binding Candidates in High-Throughput Protein Libraries.

Identifying the interactome for a protein of interest is challenging due to the large number of possible binders. High-throughput experimental approaches narrow down possible binding partners but often include false positives. Furthermore, they provide no information about what the binding region is (e.g., the binding epitope). We introduce a novel computational pipeline based on an AlphaFold2 (AF) Competitive Binding Assay (AF-CBA) to identify proteins that bind a target of interest from a pull-down experiment and the binding epitope. Our focus is on proteins that bind the Extraterminal (ET) domain of Bromo and Extraterminal domain (BET) proteins, but we also introduce nine additional systems to show transferability to other peptide-protein systems. We describe a series of limitations to the methodology based on intrinsic deficiencies of AF and AF-CBA to help users identify scenarios where the approach will be most useful. Given the method's speed and accuracy, we anticipate its broad applicability to identify binding epitope regions among potential partners, setting the stage for experimental verification.

Restraint validation of biomolecular structures determined by NMR in the Protein Data Bank.

Biomolecular structure analysis from experimental NMR studies generally relies on restraints derived from a combination of experimental and knowledge-based data. A challenge for the structural biology community has been a lack of standards for representing these restraints, preventing the establishment of uniform methods of model-vs-data structure validation against restraints and limiting interoperability between restraint-based structure modeling programs. The NEF and NMR-STAR formats provide a standardized approach for representing commonly used NMR restraints. Using these restraint formats, a standardized validation system for assessing structural models of biopolymers against restraints has been developed and implemented in the wwPDB OneDep data deposition-validation-biocuration system. The resulting wwPDB restraint violation report provides a model vs. data assessment of biomolecule structures determined using distance and dihedral restraints, with extensions to other restraint types currently being implemented. These tools are useful for assessing NMR models, as well as for assessing biomolecular structure predictions based on distance restraints.

Restraint Validation of Biomolecular Structures Determined by NMR in the Protein Data Bank.

Biomolecular structure analysis from experimental NMR studies generally relies on restraints derived from a combination of experimental and knowledge-based data. A challenge for the structural biology community has been a lack of standards for representing these restraints, preventing the establishment of uniform methods of model-vs-data structure validation against restraints and limiting interoperability between restraint-based structure modeling programs. The NMR exchange (NEF) and NMR-STAR formats provide a standardized approach for representing commonly used NMR restraints. Using these restraint formats, a standardized validation system for assessing structural models of biopolymers against restraints has been developed and implemented in the wwPDB OneDep data deposition-validation-biocuration system. The resulting wwPDB Restraint Violation Report provides a model vs. data assessment of biomolecule structures determined using distance and dihedral restraints, with extensions to other restraint types currently being implemented. These tools are useful for assessing NMR models, as well as for assessing biomolecular structure predictions based on distance restraints.

An isothermal calorimetry assay for determining steady state kinetic and enzyme inhibition parameters for SARS-CoV-2 3CL-protease.

This manuscript describes the application of Isothermal Titration Calorimetry (ITC) to characterize the kinetics of 3CL pro from the Severe Acute Respiratory Syndrome CoronaVirus-2 (SARS-CoV-2) and its inhibition by Ensitrelvir, a known non-covalent inhibitor. 3CL pro is the main protease that plays a crucial role of producing the whole array of proteins necessary for the viral infection that caused the spread of COVID-19, responsible for millions of deaths worldwide as well as global economic and healthcare crises in recent years. The proposed calorimetric method proved to have several advantages over the two types of enzymatic assays so far applied to this system, namely Förster Resonance Energy Transfer (FRET) and Liquid Chromatography-Mass Spectrometry (LC-MS). The developed ITC-based assay provided a rapid response to 3CL pro activity, which was used to directly derive the kinetic enzymatic constants K M and k cat reliably and reproducibly, as well as their temperature dependence, from which the activation energy of the reaction was obtained for the first time. The assay further revealed the existence of two modes of inhibition of 3CL pro by Ensitrelvir, namely a competitive mode as previously inferred by crystallography as well as an unprecedented uncompetitive mode, further yielding the respective inhibition constants with high precision. The calorimetric method described in this paper is thus proposed to be generally and widely used in the discovery and development of drugs targeting 3CL pro .

The NS1 protein of influenza B virus binds 5'-triphosphorylated dsRNA to suppress RIG-I activation and the host antiviral response.

Influenza A and B viruses overcome the host antiviral response to cause a contagious and often severe human respiratory disease. Here, integrative structural biology and biochemistry studies on non-structural protein 1 of influenza B virus (NS1B) reveal a previously unrecognized viral mechanism for innate immune evasion. Conserved basic groups of its C-terminal domain (NS1B-CTD) bind 5'triphosphorylated double-stranded RNA (5'-ppp-dsRNA), the primary pathogen-associated feature that activates the host retinoic acid-inducible gene I protein (RIG-I) to initiate interferon synthesis and the cellular antiviral response. Like RIG-I, NS1B-CTD preferentially binds blunt-end 5'ppp-dsRNA. NS1B-CTD also competes with RIG-I for binding 5'ppp-dsRNA, and thus suppresses activation of RIG-I's ATPase activity. Although the NS1B N-terminal domain also binds dsRNA, it utilizes a different binding mode and lacks 5'ppp-dsRNA end preferences. In cells infected with wild-type influenza B virus, RIG-I activation is inhibited. In contrast, RIG-I activation and the resulting phosphorylation of transcription factor IRF-3 are not inhibited in cells infected with a mutant virus encoding NS1B with a R208A substitution it its CTD that eliminates its 5'ppp-dsRNA binding activity. These results reveal a novel mechanism in which NS1B binds 5'ppp-dsRNA to inhibit the RIG-I antiviral response during influenza B virus infection, and open the door to new avenues for antiviral drug discovery.

Sifting Through the Noise: A Computational Pipeline for Accurate Prioritization of Protein-Protein Binding Candidates in High-Throughput Protein Libraries.

Identifying the interactome for a protein of interest is challenging due to the large number of possible binders. High-throughput experimental approaches narrow down possible binding partners, but often include false positives. Furthermore, they provide no information about what the binding region is (e.g. the binding epitope). We introduce a novel computational pipeline based on an AlphaFold2 (AF) Competition Assay (AF-CBA) to identify proteins that bind a target of interest from a pull-down experiment, along with the binding epitope. Our focus is on proteins that bind the Extraterminal (ET) domain of Bromo and Extraterminal domain (BET) proteins, but we also introduce nine additional systems to show transferability to other peptide-protein systems. We describe a series of limitations to the methodology based on intrinsic deficiencies to AF and AF-CBA, to help users identify scenarios where the approach will be most useful. Given the speed and accuracy of the methodology, we expect it to be generally applicable to facilitate target selection for experimental verification starting from high-throughput protein libraries.

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.

Breaking the conformational ensemble barrier: Ensemble structure modeling challenges in CASP15.

For the first time, the 2022 CASP (Critical Assessment of Structure Prediction) community experiment included a section on computing multiple conformations for protein and RNA structures. There was full or partial success in reproducing the ensembles for four of the nine targets, an encouraging result. For protein structures, enhanced sampling with variations of the AlphaFold2 deep learning method was by far the most effective approach. One substantial conformational change caused by a single mutation across a complex interface was accurately reproduced. In two other assembly modeling cases, methods succeeded in sampling conformations near to the experimental ones even though environmental factors were not included in the calculations. An experimentally derived flexibility ensemble allowed a single accurate RNA structure model to be identified. Difficulties included how to handle sparse or low-resolution experimental data and the current lack of effective methods for modeling RNA/protein complexes. However, these and other obstacles appear addressable.

Representing structures of the multiple conformational states of proteins.

Biomolecules exhibit dynamic behavior that single-state models of their structures cannot fully capture. We review some recent advances for investigating multiple conformations of biomolecules, including experimental methods, molecular dynamics simulations, and machine learning. We also address the challenges associated with representing single- and multiple-state models in data archives, with a particular focus on NMR structures. Establishing standardized representations and annotations will facilitate effective communication and understanding of these complex models to the broader scientific community.

Protocol for production and purification of SARS-CoV-2 3CLpro.

3CLpro protease from SARS-CoV-2 is a primary target for COVID-19 antiviral drug development. Here, we present a protocol for 3CLpro production in Escherichia coli. We describe steps to purify 3CLpro, expressed as a fusion with the Saccharomyces cerevisiae SUMO protein, with yields up to 120 mg L-1 following cleavage. The protocol also provides isotope-enriched samples suitable for nuclear magnetic resonance (NMR) studies. We also present methods to characterize 3CLpro by mass spectrometry, X-ray crystallography, heteronuclear NMR, and a Förster-resonance-energy-transfer-based enzyme assay. For complete details on the use and execution of this protocol, please refer to Bafna et al.1.

Cell-permeable chameleonic peptides: Exploiting conformational dynamics in de novo cyclic peptide design.

Membrane-traversing peptides offer opportunities for targeting intracellular proteins and oral delivery. Despite progress in understanding the mechanisms underlying membrane traversal in natural cell-permeable peptides, there are still several challenges to designing membrane-traversing peptides with diverse shapes and sizes. Conformational flexibility appears to be a key determinant of membrane permeability of large macrocycles. We review recent developments in the design and validation of chameleonic cyclic peptides, which can switch between alternative conformations to enable improved permeability through cell membranes, while still maintaining reasonable solubility and exposed polar functional groups for target protein binding. Finally, we discuss the principles, strategies, and practical considerations for rational design, discovery, and validation of permeable chameleonic peptides.

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.

Structure Determination of Challenging Protein-Peptide Complexes Combining NMR Chemical Shift Data and Molecular Dynamics Simulations.

Intrinsically disordered regions of proteins often mediate important protein-protein interactions. However, the folding-upon-binding nature of many polypeptide-protein interactions limits the ability of modeling tools to predict the three-dimensional structures of such complexes. To address this problem, we have taken a tandem approach combining NMR chemical shift data and molecular simulations to determine the structures of peptide-protein complexes. Here, we use the MELD (Modeling Employing Limited Data) technique applied to polypeptide complexes formed with the extraterminal domain (ET) of bromo and extraterminal domain (BET) proteins, which exhibit a high degree of binding plasticity. This system is particularly challenging as the binding process includes allosteric changes across the ET receptor upon binding, and the polypeptide binding partners can adopt different conformations (e.g., helices and hairpins) in the complex. In a blind study, the new approach successfully modeled bound-state conformations and binding poses, using only protein receptor backbone chemical shift data, in excellent agreement with experimentally determined structures for moderately tight (Kd ∼100 nM) binders. The hybrid MELD + NMR approach required additional peptide ligand chemical shift data for weaker (Kd ∼250 µM) peptide binding partners. AlphaFold also successfully predicts the structures of some of these peptide-protein complexes. However, whereas AlphaFold can provide qualitative peptide rankings, MELD can directly estimate relative binding affinities. The hybrid MELD + NMR approach offers a powerful new tool for structural analysis of protein-polypeptide complexes involving disorder-to-order transitions upon complex formation, which are not successfully modeled with most other complex prediction methods, providing both the 3D structures of peptide-protein complexes and their relative binding affinities.

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.

Structural similarities between SARS-CoV2 3CLpro and other viral proteases suggest potential lead molecules for developing broad spectrum antivirals.

Considering the significant impact of the recent COVID-19 outbreak, development of broad-spectrum antivirals is a high priority goal to prevent future global pandemics. Antiviral development processes generally emphasize targeting a specific protein from a particular virus. However, some antiviral agents developed for specific viral protein targets may exhibit broad spectrum antiviral activity, or at least provide useful lead molecules for broad spectrum drug development. There is significant potential for repurposing a wide range of existing viral protease inhibitors to inhibit the SARS-CoV2 3C-like protease (3CLpro). If effective even as relatively weak inhibitors of 3CLpro, these molecules can provide a diverse and novel set of scaffolds for new drug discovery campaigns. In this study, we compared the sequence- and structure-based similarity of SARS-CoV2 3CLpro with proteases from other viruses, and identified 22 proteases with similar active-site structures. This structural similarity, characterized by secondary-structure topology diagrams, is evolutionarily divergent within taxonomically related viruses, but appears to result from evolutionary convergence of protease enzymes between virus families. Inhibitors of these proteases that are structurally similar to the SARS-CoV2 3CLpro protease were identified and assessed as potential inhibitors of SARS-CoV2 3CLpro protease by virtual docking. Several of these molecules have docking scores that are significantly better than known SARS-CoV2 3CLpro inhibitors, suggesting that these molecules are also potential inhibitors of the SARS-CoV2 3CLpro protease. Some have been previously reported to inhibit SARS-CoV2 3CLpro. The results also suggest that established inhibitors of SARS-CoV2 3CLpro may be considered as potential inhibitors of other viral 3C-like proteases.

Accurate de novo design of membrane-traversing macrocycles.

We use computational design coupled with experimental characterization to systematically investigate the design principles for macrocycle membrane permeability and oral bioavailability. We designed 184 6-12 residue macrocycles with a wide range of predicted structures containing noncanonical backbone modifications and experimentally determined structures of 35; 29 are very close to the computational models. With such control, we show that membrane permeability can be systematically achieved by ensuring all amide (NH) groups are engaged in internal hydrogen bonding interactions. 84 designs over the 6-12 residue size range cross membranes with an apparent permeability greater than 1 × 10-6 cm/s. Designs with exposed NH groups can be made membrane permeable through the design of an alternative isoenergetic fully hydrogen-bonded state favored in the lipid membrane. The ability to robustly design membrane-permeable and orally bioavailable peptides with high structural accuracy should contribute to the next generation of designed macrocycle therapeutics.