PINE-SPARKY.2 for automated NMR-based protein structure research

Predicting multiple conformations via sequence clustering and AlphaFold2

AlphaFold2 (ref. 1) has revolutionized structural biology by accurately predicting single structures of proteins. However, a protein's biological function often depends on multiple conformational substates2, and disease-causing point mutations often cause population changes within these substates3,4. We demonstrate that clustering a multiple-sequence alignment by sequence similarity enables AlphaFold2 to sample alternative states of known metamorphic proteins with high confidence. Using this method, named AF-Cluster, we investigated the evolutionary distribution of predicted structures for the metamorphic protein KaiB5 and found that predictions of both conformations were distributed in clusters across the KaiB family. We used nuclear magnetic resonance spectroscopy to confirm an AF-Cluster prediction: a cyanobacteria KaiB variant is stabilized in the opposite state compared with the more widely studied variant. To test AF-Cluster's sensitivity to point mutations, we designed and experimentally verified a set of three mutations predicted to flip KaiB from Rhodobacter sphaeroides from the ground to the fold-switched state. Finally, screening for alternative states in protein families without known fold switching identified a putative alternative state for the oxidoreductase Mpt53 in Mycobacterium tuberculosis. Further development of such bioinformatic methods in tandem with experiments will probably have a considerable impact on predicting protein energy landscapes, essential for illuminating biological function.

Backbone NMR assignment of the yeast expressed Fab fragment of the NISTmAb reference antibody

The monoclonal antibody (mAb) protein class has become a primary therapeutic platform for the production of new life saving drug products. MAbs are comprised of two domains: the antigen-binding fragment (Fab) and crystallizable fragment (Fc). Despite the success in the clinic, NMR assignments of the complete Fab domain have been elusive, in part due to problems in production of properly folded, triply-labeled ²H,¹³C,¹⁵N Fab domain. Here, we report the successful recombinant expression of a triply-labeled Fab domain, derived from the standard IgG1κ known as NISTmAb, in yeast. Using the ²H,¹³C,¹⁵N Fab domain, we assigned 94% of the ¹H, ¹³C, and ¹⁵N backbone atoms.

A Novel Mechanism of Coactivator Recruitment by the Nurr1 Nuclear Receptor

Nuclear receptors constitute one of the largest families of transcription factors that regulate genes in metazoans in response to small molecule ligands. Many receptors harbor two transactivation domains, one at each end of the protein sequence. Whereas the molecular mechanisms of transactivation mediated by the ligand-binding domain at the C-terminus of the protein are generally well established, the mechanism involving the N-terminal domain called activation function 1 (AF1) has remained elusive. Previous studies implicated the AF1 domain as a significant contributor towards the overall transcriptional activity of the NR4A family of nuclear receptors and suggested that the steroid receptor coactivators (SRCs) play an important role in this process. Here we show that a short segment within the AF1 domain of the NR4A receptor Nurr1 can directly engage with the SRC1 PAS-B domain. We also show that this segment forms a helix upon binding to a largely hydrophobic groove on PAS-B, overlapping with the surface engaged by the STAT6 transcription factor, suggesting that this mode of recruitment could be shared by diverse transcription factors including other nuclear receptors.

ssPINE: Probabilistic Algorithm for Automated Chemical Shift Assignment of Solid-State NMR Data from Complex Protein Systems

The heightened dipolar interactions in solids render solid-state NMR (ssNMR) spectra more difficult to interpret than solution NMR spectra. On the other hand, ssNMR does not suffer from severe molecular weight limitations like solution NMR. In recent years, ssNMR has undergone rapid technological developments that have enabled structure-function studies of increasingly larger biomolecules, including membrane proteins. Current methodology includes stable isotope labeling schemes, non-uniform sampling with spectral reconstruction, faster magic angle spinning, and innovative pulse sequences that capture different types of interactions among spins. However, computational tools for the analysis of complex ssNMR data from membrane proteins and other challenging protein systems have lagged behind those for solution NMR. Before a structure can be determined, thousands of signals from individual types of multidimensional ssNMR spectra of samples, which may have differing isotopic composition, must be recognized, correlated, categorized, and eventually assigned to atoms in the chemical structure. To address these tedious steps, we have developed an automated algorithm for ssNMR spectra called "ssPINE". The ssPINE software accepts the sequence of the protein plus peak lists from a variety of ssNMR experiments as inputs and offers automated backbone and side-chain assignments. The alpha version of ssPINE, which we describe here, is freely available through a web submission form.

Solution structure and dynamics of the mitochondrial-targeted GTPase-activating protein (GAP) VopE by an integrated NMR/SAXS approach

The bacterial pathogen Vibrio cholerae use a type III secretion system to inject effector proteins into a host cell. Recently, a putative Toxic GTPase Activating Protein (ToxGAP) called Vibrio outer protein E (VopE) was identified as a T3SS substrate and virulence factor that affected host mitochondrial dynamics and immune response. However, biophysical and structural characterization has been absent. Here, we describe solution NMR structure of the putative GTPase-activating protein (GAP) domain (73-204) of VopE. Using size exclusion chromatography coupled with small-angle x-ray scattering and residual dipolar coupling data, we restrained the MD process to efficiently determine the overall fold and improve the quality of the output calculated structures. Comparing the structure of VopE with other ToxGAP's revealed a similar overall fold with several features unique to VopE. Specifically, the "Bulge 1," α1 helix, and noteworthy "backside linker" elements on the N-terminus are dissimilar to the other ToxGAP's. By using NMR relaxation dispersion experiments, we demonstrate that these regions undergo motions on a > 6 s-1 timescale. Based on the disposition of these mobile regions relative to the putative catalytic arginine residue, we hypothesize that the protein may undergo structural changes to bind cognate GTPases.

Backbone and sidechain 1H, 15N and 13C resonance assignments of the free and RNA-bound tandem zinc finger domain of the tristetraprolin family member from Selaginella moellendorffii

Members of the tristetraprolin (TTP) family of RNA binding proteins (RBPs) regulate the metabolism of a variety of mRNA targets. In mammals, these proteins modulate many physiological processes, including immune cell activation, hematopoiesis, and embryonic development. Regulation of mRNA stability by these proteins requires that the tandem zinc finger (TZF) domain binds initially and directly to target mRNAs, ultimately leading to their deadenylation and decay. Proteins of this type throughout eukarya possess a highly conserved TZF domain, suggesting that they are all capable of high-affinity RNA binding. However, the mechanism of TTP-mediated mRNA decay is largely undefined. Given the vital role that these TTP family proteins play in maintaining RNA homeostasis throughout eukaryotes, we focused here on the first, key step in this process: recognition and binding of the TZF domain to target RNA. For these studies, we chose a primitive plant, the spikemoss Selaginella moellendorffii, which last shared a common ancestor with humans more than a billion years ago. Here we report the near complete backbone and side chain resonance assignments of the spikemoss TZF domain, including: (1) the assignment of the RNA-TZF domain complex, representing one of only two data sets currently available for the entire TTP family of proteins; and (2) the first NMR resonance assignments of the entire TZF domain, in the RNA-free form. This work will serve as the basis for further NMR structural investigations aimed at gaining insights into the process of RNA recognition and the mechanisms of TTP-mediated mRNA decay.

Cuprous Ions can Disrupt Structure and Functions of the RING Finger Domain of RNF11

Copper is an essential element that plays crucial roles in a variety of biological processes, while excessive copper is harmful to cells. RNF11 is a RING finger protein associated with...

POKY software tools encapsulating assignment strategies for solution and solid-state protein NMR data

•

New tools support efficient analysis of solution and solid-state NMR spectra of proteins.

•

POKY integrates a powerful suite of software packages for automated assignments.

•

The Versatile Assigner module validates assignments through probabilistic analysis.

•

The operation of these tools is supported by on-line guidance.

•

The performance of these tools is evaluated in reference to competing software.

NMR spectroscopy provides structural and functional information about biomolecules and their complexes. The complexity of these systems can make the NMR data difficult to interpret, particularly for newer users of NMR technology, who may have limited understanding of the tools available and how they are used. To alleviate this problem, we have created software based on standardized workflows for both solution and solid-state NMR spectroscopy of proteins. These tools assist with manual and automated peak picking and with chemical shift assignment and validation. They provide users with an optimized path through spectral analysis that can help them perform the necessary tasks more efficiently.

Toho-1 β-lactamase: backbone chemical shift assignments and changes in dynamics upon binding with avibactam

Backbone chemical shift assignments for the Toho-1 β-lactamase (263 amino acids, 28.9 kDa) are reported based on triple resonance solution-state NMR experiments performed on a uniformly 2H,13C,15N-labeled sample. These assignments allow for subsequent site-specific characterization at the chemical, structural, and dynamical levels. At the chemical level, titration with the non-β-lactam β-lactamase inhibitor avibactam is found to give chemical shift perturbations indicative of tight covalent binding that allow for mapping of the inhibitor binding site. At the structural level, protein secondary structure is predicted based on the backbone chemical shifts and protein residue sequence using TALOS-N and found to agree well with structural characterization from X-ray crystallography. At the dynamical level, model-free analysis of 15N relaxation data at a single field of 16.4 T reveals well-ordered structures for the ligand-free and avibactam-bound enzymes with generalized order parameters of ~ 0.85. Complementary relaxation dispersion experiments indicate that there is an escalation in motions on the millisecond timescale in the vicinity of the active site upon substrate binding. The combination of high rigidity on short timescales and active site flexibility on longer timescales is consistent with hypotheses for achieving both high catalytic efficiency and broad substrate specificity: the induced active site dynamics allows variously sized substrates to be accommodated and increases the probability that the optimal conformation for catalysis will be sampled.

CHESPA/CHESCA-SPARKY: automated NMR data analysis plugins for SPARKY to map protein allostery

MOTIVATION

Correlated Nuclear Magnetic Resonance (NMR) chemical shift changes identified through the CHEmical Shift Projection Analysis (CHESPA) and CHEmical Shift Covariance Analysis (CHESCA) reveal pathways of allosteric transitions in biological macromolecules. To address the need for an automated platform that implements CHESPA and CHESCA and integrates them with other NMR analysis software packages, we introduce here integrated plugins for NMRFAM-SPARKY that implement the seamless detection and visualization of allosteric networks.

AVAILABILITY AND IMPLEMENTATION

CHESCA-SPARKY and CHESPA-SPARKY are available in the latest version of NMRFAM-SPARKY from the National Magnetic Resonance Facility at Madison (http://pine.nmrfam.wisc.edu/download_packages.html), the NMRbox Project (https://nmrbox.org) and to subscribers to the SBGrid (https://sbgrid.org). The assigned spectra involved in this study and tutorial videos using this dataset are available at https://sites.google.com/view/chescachespa-sparky.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics Online.

Resonance assignments and secondary structure of thermophile single-stranded DNA binding protein from Sulfolobus solfataricus at 323K

Single-stranded DNA (ssDNA)-binding proteins (SSBs) are essential for DNA replication, recombination, and repair processes in all organisms. Sulfolobus solfataricus (S. solfataricus), a hyperthermophilic species, overexpresses its SSB (S. solfataricus SSB (SsoSSB)) to protect ssDNA during DNA metabolisms. Even though the crystal structure of apo SsoSSB and its ssDNA-bound solution structure have been reported at room temperature, structural information at high temperature is not yet available. To find out how SsoSSB maintains its structure and ssDNA binding affinity at high temperatures, we performed multidimensional NMR experiments for SsoSSB at 323K. In this study, we present the backbone and side chain chemical shifts and predict the secondary structure of SsoSSB from the chemical shifts. We found that SsoSSB is ordered, even at high temperatures, and has the same fold at high temperature as at room temperature. Our data will help improve structural analyses and our understanding of the features of thermophilic proteins.

1H, 13C, and 15N backbone and side chain chemical shift assignments of the SARS-CoV-2 non-structural protein 7

The SARS-CoV-2 genome encodes for approximately 30 proteins. Within the international project covid19-nmr, we distribute the spectroscopic analysis of the viral proteins and RNA. Here, we report NMR chemical shift assignments for the protein nsp7. The 83 amino acid nsp7 protein is an essential cofactor in the RNA-dependent RNA polymerase. The polymerase activity and processivity of nsp12 are greatly enhanced by binding 1 copy of nsp7 and 2 copies of nsp8 to form a 160 kD complex. A separate hexadecameric complex of nsp7 and nsp8 (8 copies of each) forms a large ring-like structure. Thus, nsp7 is an important component of several large protein complexes that are required for replication of the large and complex coronavirus genome. We here report the near-complete NMR backbone and sidechain resonance assignment (¹H,¹³C,¹⁵N) of isolated nsp7 from SARS-CoV-2 in solution. Further, we derive the secondary structure and compare it to the previously reported assignments and structure of the SARS-CoV nsp7.

POKY: a software suite for multidimensional NMR and 3D structure calculation of biomolecules

SUMMARY

The need for an efficient and cost-effective method is compelling in biomolecular NMR. To tackle this problem, we have developed the Poky suite, the revolutionized platform with boundless possibilities for advancing research and technology development in signal detection, resonance assignment, structure calculation and relaxation studies with the help of many automation and user interface tools. This software is extensible and scalable by scripting and batching as well as providing modern graphical user interfaces and a diverse range of modules right out of the box.

AVAILABILITY AND IMPLEMENTATION

Poky is freely available to non-commercial users at https://poky.clas.ucdenver.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Solution NMR Determination of the CDHR3 Rhinovirus-C Binding Domain, EC1

Cadherin Related Family Member 3 (CDHR3) is the identified and required cellular receptor for all virus isolates in the rhinovirus-C species (RV-C). Cryo-EM determinations recently resolved the atomic structure of RV-C15a, and subsequently, a complex of this virus bound to CDHR3 extracellular domain 1 (EC1), the N-terminal portion of this receptor responsible for virus interactions. The EC1 binds to a hypervariable sequence footprint on the virus surface, near the 3-fold axis of icosahedral symmetry. The key contacts involve discontinuous residues from 3 viral proteins, VP1, VP2 and VP3. That single cryo-EM EC1 structure, however, could not resolve whether the virus-receptor interface was structurally adaptable to accommodate multiple virus sequences. We now report the solution NMR determination of CDHR3 EC1, showing that this protein, in fact, is mostly inflexible, particularly in the virus-binding face. The new, higher resolution dataset identifies 3 cis-Pro residues in important loop regions, where they can influence both rigidity and overall protein conformation. The data also provide clarification about the residues involved in essential calcium ion binding, and a potential CDHR3 surface groove feature that may be involved in native protein interactions with cellular partners.

Solution structures of the Shewanella woodyi H-NOX protein in the presence and absence of soluble guanylyl cyclase stimulator IWP-051

Heme-nitric oxide/oxygen binding (H-NOX) domains bind gaseous ligands for signal transduction in organisms spanning prokaryotic and eukaryotic kingdoms. In the bioluminescent marine bacterium Shewanella woodyi (Sw), H-NOX proteins regulate quorum sensing and biofilm formation. In higher animals, soluble guanylyl cyclase (sGC) binds nitric oxide with an H-NOX domain to induce cyclase activity and regulate vascular tone, wound healing and memory formation. sGC also binds stimulator compounds targeting cardiovascular disease. The molecular details of stimulator binding to sGC remain obscure but involve a binding pocket near an interface between H-NOX and coiled-coil domains. Here, we report the full NMR structure for CO-ligated Sw H-NOX in the presence and absence of stimulator compound IWP-051, and its backbone dynamics. Nonplanar heme geometry was retained using a semi-empirical quantum potential energy approach. Although IWP-051 binding is weak, a single binding conformation was found at the interface of the two H-NOX subdomains, near but not overlapping with sites identified in sGC. Binding leads to rotation of the subdomains and closure of the binding pocket. Backbone dynamics are similar across both domains except for two helix-connecting loops, which display increased dynamics that are further enhanced by compound binding. Structure-based sequence analyses indicate high sequence diversity in the binding pocket, but the pocket itself appears conserved among H-NOX proteins. The largest dynamical loop lies at the interface between Sw H-NOX and its binding partner as well as in the interface with the coiled coil in sGC, suggesting a critical role for the loop in signal transduction.

Backbone assignments and conformational dynamics in the S. typhimurium tryptophan synthase α-subunit from solution-state NMR

Backbone assignments for the isolated α-subunit of Salmonella typhimurium tryptophan synthase (TS) are reported based on triple resonance solution-state NMR experiments on a uniformly 2H,13C,15N-labeled sample. From the backbone chemical shifts, secondary structure and random coil index order parameters (RCI-S2) are predicted. Titration with the 3-indole-D-glycerol 3'-phosphate analog, N-(4'-trifluoromethoxybenzenesulfonyl)-2-aminoethyl phosphate (F9), leads to chemical shift perturbations indicative of conformational changes from which an estimate of the dissociation constant is obtained. Comparisons of the backbone chemical-shifts, RCI-S2 values, and site-specific relaxation times with and without F9 reveal allosteric changes including modulation in secondary structures and loop rigidity induced upon ligand binding. A comparison is made to the X-ray crystal structure of the α-subunit in the full TS αββα bi-enzyme complex and to two new X-ray crystal structures of the isolated TS α-subunit reported in this work.

Solution structure of human myeloid-derived growth factor suggests a conserved function in the endoplasmic reticulum

Human myeloid-derived growth factor (hMYDGF) is a 142-residue protein with a C-terminal endoplasmic reticulum (ER) retention sequence (ERS). Extracellular MYDGF mediates cardiac repair in mice after anoxic injury. Although homologs of hMYDGF are found in eukaryotes as distant as protozoans, its structure and function are unknown. Here we present the NMR solution structure of hMYDGF, which consists of a short α-helix and ten β-strands distributed in three β-sheets. Conserved residues map to the unstructured ERS, loops on the face opposite the ERS, and the surface of a cavity underneath the conserved loops. The only protein or portion of a protein known to have a similar fold is the base domain of VNN1. We suggest, in analogy to the tethering of the VNN1 nitrilase domain to the plasma membrane via its base domain, that MYDGF complexed to the KDEL receptor binds cargo via its conserved residues for transport to the ER.

Backbone resonance assignments and secondary structure of Ebola nucleoprotein 600-739 construct

Ebola viral infections have resulted in several deadly epidemics in recent years in West and Central Africa. Because only one of the seven proteins encoded by the viral genome possesses enzymatic activity, disruption of protein-protein interactions is a promising route for antiviral drug development. We carried out a screening campaign to identify small, drug-like compounds that bind to the C-terminal region of the multifunctional Ebola nucleoprotein (eNP) with the objective of discovering ones that disrupt its binding to other Ebola proteins or to the single-stranded RNA genome. In the course of this effort we assigned the backbone ¹H, ¹⁵N, and ¹³C resonances of residues 600‒739, the region that contains the critical eVP30 binding region 600‒615 targeted by host factors, and used the assigned chemical shifts to predict secondary structural features and peptide dynamics. This work supports and extends the previous X-ray crystal structures and NMR studies of residues 641‒739. We found that the 600‒739 domain consists of separate regions that are largely disordered and ordered.

I-PINE web server: an integrative probabilistic NMR assignment system for proteins

Various methods for understanding the structural and dynamic properties of proteins rely on the analysis of their NMR chemical shifts. These methods require the initial assignment of NMR signals to particular atoms in the sequence of the protein, a step that can be very time-consuming. The probabilistic interaction network of evidence (PINE) algorithm for automated assignment of backbone and side chain chemical shifts utilizes a Bayesian probabilistic network model that analyzes sequence data and peak lists from multiple NMR experiments. PINE, which is one of the most popular and reliable automated chemical shift assignment algorithms, has been available to the protein NMR community for longer than a decade. We announce here a new web server version of PINE, called Integrative PINE (I-PINE), which supports more types of NMR experiments than PINE (including three-dimensional nuclear Overhauser enhancement and four-dimensional J-coupling experiments) along with more comprehensive visualization of chemical shift based analysis of protein structure and dynamics. The I-PINE server is freely accessible at http://i-pine.nmrfam.wisc.edu . Help pages and tutorial including browser capability are available at: http://i-pine.nmrfam.wisc.edu/instruction.html . Sample data that can be used for testing the web server are available at: http://i-pine.nmrfam.wisc.edu/examples.html .

15N, 13C and 1H resonance assignments of FKBP12 proteins from the pathogenic fungi Mucor circinelloides and Aspergillus fumigatus

Invasive fungal infections are a leading cause of death in immunocompromised patients and remain difficult to treat since fungal pathogens, like mammals, are eukaryotes and share many orthologous proteins. As a result, current antifungal drugs have limited clinical value, are sometimes toxic, can adversely affect human reaction pathways and are increasingly ineffective due to emerging resistance. One potential antifungal drug, FK506, establishes a ternary complex between the phosphatase, calcineurin, and the 12-kDa peptidyl-prolyl isomerase FK506-binding protein, FKBP12. It has been well established that calcineurin, highly conserved from yeast to mammals, is necessary for invasive fungal disease and is inhibited when in complex with FK506/FKBP12. Unfortunately, FK506 is also immunosuppressive in humans, precluding its usage as an antifungal drug, especially in immunocompromised patients. Whereas the homology between human and fungal calcineurin proteins is > 80%, the human and fungal FKBP12s share 48-58% sequence identity, making them more amenable candidates for drug targeting efforts. Here we report the backbone and sidechain NMR assignments of recombinant FKBP12 proteins from the pathogenic fungi Mucor circinelloides and Aspergillus fumigatus in the apo form and compare these to the backbone assignments of the FK506 bound form. In addition, we report the backbone assignments of the apo and FK506 bound forms of the Homo sapiens FKBP12 protein for evaluation against the fungal forms. These data are the first steps towards defining, at a residue specific level, the impacts of FK506 binding to fungal and mammalian FKBP12 proteins. Our data highlight differences between the human and fungal FKBP12s that could lead to the design of more selective anti-fungal drugs.

Applications of In-Cell NMR in Structural Biology and Drug Discovery

In-cell nuclear magnetic resonance (NMR) is a method to provide the structural information of a target at an atomic level under physiological conditions and a full view of the conformational changes of a protein caused by ligand binding, post-translational modifications or protein⁻protein interactions in living cells. Previous in-cell NMR studies have focused on proteins that were overexpressed in bacterial cells and isotopically labeled proteins injected into oocytes of Xenopus laevis or delivered into human cells. Applications of in-cell NMR in probing protein modifications, conformational changes and ligand bindings have been carried out in mammalian cells by monitoring isotopically labeled proteins overexpressed in living cells. The available protocols and successful examples encourage wide applications of this technique in different fields such as drug discovery. Despite the challenges in this method, progress has been made in recent years. In this review, applications of in-cell NMR are summarized. The successful applications of this method in mammalian and bacterial cells make it feasible to play important roles in drug discovery, especially in the step of target engagement.

Non-Uniform and Absolute Minimal Sampling for High-Throughput Multidimensional NMR Applications

Many biomolecular NMR applications can benefit from the faster acquisition of multidimensional NMR data with high resolution and their automated analysis and interpretation. In recent years, a number of non-uniform sampling (NUS) approaches have been introduced for the reconstruction of multidimensional NMR spectra, such as compressed sensing, thereby bypassing traditional Fourier-transform processing. Such approaches are applicable to both biomacromolecules and small molecules and their complex mixtures and can be combined with homonuclear decoupling (pure shift) and covariance processing. For homonuclear 2D TOCSY experiments, absolute minimal sampling (AMS) permits the drastic shortening of measurement times necessary for high-throughput applications for identification and quantification of components in complex biological mixtures in the field of metabolomics. Such TOCSY spectra can be comprehensively represented by graphic theoretical maximal cliques for the identification of entire spin systems and their subsequent query against NMR databases. Integration of these methods in webservers permits the rapid and reliable identification of mixture components. Recent progress is reviewed in this Minireview.