Site-Specific Labeling and 19F NMR Provide Direct Evidence for Dynamic Behavior of the Anthrax Toxin Pore ϕ-Clamp Structure.

The anthrax toxin protective antigen (PA), the membrane binding and pore-forming component of the anthrax toxin, was studied using 19F NMR. We site-specifically labeled PA with p-fluorophenylalanine (pF-Phe) at Phe427, a critically important residue that comprises the ϕ-clamp that is required for translocation of edema factor (EF) and lethal factor (LF) into the host cell cytosol. We utilized 19F NMR to follow low-pH-induced structural changes in the prepore, alone and bound to the N-terminal PA binding domain of LF, LFN. Our studies indicate that pF-Phe427 is dynamic in the prepore state and then becomes more dynamic in the transition to the pore. An increase in dynamic behavior at the ϕ-clamp may provide the necessary room for movement needed in translocating EF and LF into the cell cytosol.

Uncovering a membrane-distal conformation of KRAS available to recruit RAF to the plasma membrane.

The small GTPase KRAS is localized at the plasma membrane where it functions as a molecular switch, coupling extracellular growth factor stimulation to intracellular signaling networks. In this process, KRAS recruits effectors, such as RAF kinase, to the plasma membrane where they are activated by a series of complex molecular steps. Defining the membrane-bound state of KRAS is fundamental to understanding the activation of RAF kinase and in evaluating novel therapeutic opportunities for the inhibition of oncogenic KRAS-mediated signaling. We combined multiple biophysical measurements and computational methodologies to generate a consensus model for authentically processed, membrane-anchored KRAS. In contrast to the two membrane-proximal conformations previously reported, we identify a third significantly populated state using a combination of neutron reflectivity, fast photochemical oxidation of proteins (FPOP), and NMR. In this highly populated state, which we refer to as "membrane-distal" and estimate to comprise ∼90% of the ensemble, the G-domain does not directly contact the membrane but is tethered via its C-terminal hypervariable region and carboxymethylated farnesyl moiety, as shown by FPOP. Subsequent interaction of the RAF1 RAS binding domain with KRAS does not significantly change G-domain configurations on the membrane but affects their relative populations. Overall, our results are consistent with a directional fly-casting mechanism for KRAS, in which the membrane-distal state of the G-domain can effectively recruit RAF kinase from the cytoplasm for activation at the membrane.

Probabilistic identification of saccharide moieties in biomolecules and their protein complexes.

The chemical composition of saccharide complexes underlies their biomedical activities as biomarkers for cardiometabolic disease, various types of cancer, and other conditions. However, because these molecules may undergo major structural modifications, distinguishing between compounds of saccharide and non-saccharide origin becomes a challenging computational problem that hinders the aggregation of information about their bioactive moieties. We have developed an algorithm and software package called "Cheminformatics Tool for Probabilistic Identification of Carbohydrates" (CTPIC) that analyzes the covalent structure of a compound to yield a probabilistic measure for distinguishing saccharides and saccharide-derivatives from non-saccharides. CTPIC analysis of the RCSB Ligand Expo (database of small molecules found to bind proteins in the Protein Data Bank) led to a substantial increase in the number of ligands characterized as saccharides. CTPIC analysis of Protein Data Bank identified 7.7% of the proteins as saccharide-binding. CTPIC is freely available as a webservice at (http://ctpic.nmrfam.wisc.edu).

Fragment screening targeting Ebola virus nucleoprotein C-terminal domain identifies lead candidates.

The Ebola Virus is a causative agent of viral hemorrhagic fever outbreaks and a potential global health risk. The outbreak in West Africa (2013-2016) led to 11,000+ deaths and 30,000+ Ebola infected individuals. The current outbreak in the Democratic Republic of Congo (DRC) with 3000+ confirmed cases and 2000+ deaths attributed to Ebola virus infections provides a reminder that innovative countermeasures are still needed. Ebola virus encodes 7 open reading frames (ORFs). Of these, the nucleocapsid protein (eNP) encoded by the first ORF plays many significant roles, including a role in viral RNA synthesis. Here we describe efforts to target the C-terminal domain of eNP (eNP-CTD) that contains highly conserved residues 641-739 as a pan-Ebola antiviral target. Interactions of eNP-CTD with Ebola Viral Protein 30 (eVP30) and Viral Protein 40 (eVP40) have been shown to be crucial for viral RNA synthesis, virion formation, and virion transport. We used nuclear magnetic response (NMR)-based methods to screened the eNP-CTD against a fragment library. Perturbations of 1D 1H NMR spectra identified of 48 of the 439 compounds screened as potential eNP CTD interactors. Subsequent analysis of these compounds to measure chemical shift perturbations in 2D 1H,15N NMR spectra of 15N-labeled protein identified six with low millimolar affinities. All six perturbed an area consisting mainly of residues at or near the extreme C-terminus that we named "Site 1" while three other sites were perturbed by other compounds. Our findings here demonstrate the potential utility of eNP as a target, several fragment hits, and provide an experimental pipeline to validate viral-viral interactions as potential panfiloviral inhibitor targets.

Templated Collagen "Double Helices" Maintain Their Structure.

The self-assembly of collagen-mimetic peptides (CMPs) that form sticky-ended triple helices has allowed the production of surprisingly stable artificial collagen fibers and hydrogels. Assembly through sticky ends requires the recognition of a single strand by a templated strand dimer. Although CMPs and their triple helices have been studied extensively, the structure of a strand dimer is unknown. Here, we evaluate the physical characteristics of such dimers, using disulfide-templated (PPG)10 dimers as a model. Such "linked-dimers" retain their collagen-like structure even in the absence of a third strand, but only when their strands are capable of adopting a triple-helical fold. The intrinsic collagen-like structure of templated CMP pairs helps to explain the success of sticky-ended CMP association and changes the conception of new synthetic collagen designs.

Tools for Enhanced NMR-Based Metabolomics Analysis.

Metabolomics is the study of profiles of small molecules in biological fluids, cells, or organs. These profiles can be thought of as the "fingerprints" left behind from chemical processes occurring in biological systems. Because of its potential for groundbreaking applications in disease diagnostics, biomarker discovery, and systems biology, metabolomics has emerged as a rapidly growing area of research. Metabolomics investigations often, but not always, involve the identification and quantification of endogenous and exogenous metabolites in biological samples. Software tools and databases play a crucial role in advancing the rigor, robustness, reproducibility, and validation of these studies. Specifically, the establishment of a robust library of spectral signatures with unique compound descriptors and atom identities plays a key role in profiling studies based on data from nuclear magnetic resonance (NMR) spectroscopy. Here, we discuss developments leading to a rigorous basis for unique identification of compounds, reproducible numbering of atoms, the compact representation of NMR spectra of metabolites and small molecules, tools for improved compound identification, quantification and visualization, and approaches toward the goal of rigorous analysis of metabolomics data.

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 .

Automated evaluation of consistency within the PubChem Compound database.

Identification of discrepant data in aggregated databases is a key step in data curation and remediation. We have applied the ALATIS approach, which is based on the international chemical shift identifier (InChI) model, to the full PubChem Compound database to generate unique and reproducible compound and atom identifiers for all entries for which three-dimensional structures were available. This exercise also served to identify entries with discrepancies between structures and chemical formulas or InChI strings. The use of unique compound identifiers and atom nomenclature should support more rigorous links between small-molecule databases including those containing atom-specific information of the type available from crystallography and spectroscopy. The comprehensive results from this analysis are publicly available through our webserver [http://alatis.nmrfam.wisc.edu/].

Applications of Parametrized NMR Spin Systems of Small Molecules.

We have developed technology for producing accurate spectral fingerprints of small molecules through modeling of NMR spin system matrices to encapsulate their chemical shifts and scalar couplings. We describe here how libraries of these spin systems utilizing unique and reproducible atom numbering can be used to improve NMR-based ligand screening and metabolomics studies. We introduce new Web services that facilitate the analysis of NMR spectra of mixtures of small molecules to yield their identification and quantification. The library of parametrized compounds has been expanded to cover simulations of 1H NMR spectra at a variety of magnetic fields of more than 1100 compounds, included are many common metabolites and a library of drug-like molecular fragments used in ligand screening. The compound library and related Web services are freely available from http://gissmo.nmrfam.wisc.edu/ .

Conservation and Divergence in the Candida Species Biofilm Matrix Mannan-Glucan Complex Structure, Function, and Genetic Control.

Candida biofilms resist the effects of available antifungal therapies. Prior studies with Candida albicans biofilms show that an extracellular matrix mannan-glucan complex (MGCx) contributes to antifungal sequestration, leading to drug resistance. Here we implement biochemical, pharmacological, and genetic approaches to explore a similar mechanism of resistance for the three most common clinically encountered non-albicansCandida species (NAC). Our findings reveal that each Candida species biofilm synthesizes a mannan-glucan complex and that the antifungal-protective function of this complex is conserved. Structural similarities extended primarily to the polysaccharide backbone (α-1,6-mannan and ß-1,6-glucan). Surprisingly, biochemical analysis uncovered stark differences in the branching side chains of the MGCx among the species. Consistent with the structural analysis, similarities in the genetic control of MGCx production for each Candida species also appeared limited to the synthesis of the polysaccharide backbone. Each species appears to employ a unique subset of modification enzymes for MGCx synthesis, likely accounting for the observed side chain diversity. Our results argue for the conservation of matrix function among Candida spp. While biogenesis is preserved at the level of the mannan-glucan complex backbone, divergence emerges for construction of branching side chains. Thus, the MGCx backbone represents an ideal drug target for effective pan-Candida species biofilm therapy.IMPORTANCECandida species, the most common fungal pathogens, frequently grow as a biofilm. These adherent communities tolerate extremely high concentrations of antifungal agents, due in large part, to a protective extracellular matrix. The present studies define the structural, functional, and genetic similarities and differences in the biofilm matrix from the four most common Candida species. Each species synthesizes an extracellular mannan-glucan complex (MGCx) which contributes to sequestration of antifungal drug, shielding the fungus from this external assault. Synthesis of a common polysaccharide backbone appears conserved. However, subtle structural differences in the branching side chains likely rely upon unique modification enzymes, which are species specific. Our findings identify MGCx backbone synthesis as a potential pan-Candida biofilm therapeutic target.

Spin System Modeling of Nuclear Magnetic Resonance Spectra for Applications in Metabolomics and Small Molecule Screening.

The exceptionally rich information content of nuclear magnetic resonance (NMR) spectra is routinely used to identify and characterize molecules and molecular interactions in a wide range of applications, including clinical biomarker discovery, drug discovery, environmental chemistry, and metabolomics. The set of peak positions and intensities from a reference NMR spectrum generally serves as the identifying signature for a compound. Reference spectra normally are collected under specific conditions of pH, temperature, and magnetic field strength, because changes in conditions can distort the identifying signatures of compounds. A spin system matrix that parametrizes chemical shifts and coupling constants among spins provides a much richer feature set for a compound than a spectral signature based on peak positions and intensities. Spin system matrices expand the applicability of NMR spectral libraries beyond the specific conditions under which data were collected. In addition to being able to simulate spectra at any field strength, spin parameters can be adjusted to systematically explore alterations in chemical shift patterns due to variations in other experimental conditions, such as compound concentration, pH, or temperature. We present methodology and software for efficient interactive optimization of spin parameters against experimental 1D-1H NMR spectra of small molecules. We have used the software to generate spin system matrices for a set of key mammalian metabolites and are also using the software to parametrize spectra of small molecules used in NMR-based ligand screening. The software, along with optimized spin system matrix data for a growing number of compounds, is available from http://gissmo.nmrfam.wisc.edu/ .

Usb1 controls U6 snRNP assembly through evolutionarily divergent cyclic phosphodiesterase activities.

U6 small nuclear ribonucleoprotein (snRNP) biogenesis is essential for spliceosome assembly, but not well understood. Here, we report structures of the U6 RNA processing enzyme Usb1 from yeast and a substrate analog bound complex from humans. Unlike the human ortholog, we show that yeast Usb1 has cyclic phosphodiesterase activity that leaves a terminal 3' phosphate which prevents overprocessing. Usb1 processing of U6 RNA dramatically alters its affinity for cognate RNA-binding proteins. We reconstitute the post-transcriptional assembly of yeast U6 snRNP in vitro, which occurs through a complex series of handoffs involving 10 proteins (Lhp1, Prp24, Usb1 and Lsm2-8) and anti-cooperative interactions between Prp24 and Lhp1. We propose a model for U6 snRNP assembly that explains how evolutionarily divergent and seemingly antagonistic proteins cooperate to protect and chaperone the nascent snRNA during its journey to the spliceosome.The mechanism of U6 small nuclear ribonucleoprotein (snRNP) biogenesis is not well understood. Here the authors characterize the enzymatic activities and structures of yeast and human U6 RNA processing enzyme Usb1, reconstitute post-transcriptional assembly of yeast U6 snRNP in vitro, and propose a model for U6 snRNP assembly.

Progressive Stereo Locking (PSL): A Residual Dipolar Coupling Based Force Field Method for Determining the Relative Configuration of Natural Products and Other Small Molecules.

Establishing the relative configuration of a bioactive natural product represents the most challenging part in determining its structure. Residual dipolar couplings (RDCs) are sensitive probes of the relative spatial orientation of internuclear vectors. We adapted a force field structure calculation methodology to allow free sampling of both R and S configurations of the stereocenters of interest. The algorithm uses a floating alignment tensor in a simulated annealing protocol to identify the conformations and configurations that best fit experimental RDC and distance restraints (from NOE and J-coupling data). A unique configuration (for rigid molecules) or a very small number of configurations (for less rigid molecules) of the structural models having the lowest chiral angle energies and reasonable magnitudes of the alignment tensor are provided as the best predictions of the unknown configuration. For highly flexible molecules, the progressive locking of their stereocenters into their statistically dominant R or S state dramatically reduces the number of possible relative configurations. The result is verified by checking that the same configuration is obtained by initiating the locking from different regions of the molecule. For all molecules tested having known configurations (with conformations ranging from mostly rigid to highly flexible), the method accurately determined the correct configuration.

Unique identifiers for small molecules enable rigorous labeling of their atoms.

Rigorous characterization of small organic molecules in terms of their structural and biological properties is vital to biomedical research. The three-dimensional structure of a molecule, its 'photo ID', is inefficient for searching and matching tasks. Instead, identifiers play a key role in accessing compound data. Unique and reproducible molecule and atom identifiers are required to ensure the correct cross-referencing of properties associated with compounds archived in databases. The best approach to this requirement is the International Chemical Identifier (InChI). However, the current implementation of InChI fails to provide a complete standard for atom nomenclature, and incorrect use of the InChI standard has resulted in the proliferation of non-unique identifiers. We propose a methodology and associated software tools, named ALATIS, that overcomes these shortcomings. ALATIS is an adaptation of InChI, which operates fully within the InChI convention to provide unique and reproducible molecule and all atom identifiers. ALATIS includes an InChI extension for unique atom labeling of symmetric molecules. ALATIS forms the basis for improving reproducibility and unifying cross-referencing across databases.

Increasing rigor in NMR-based metabolomics through validated and open source tools.

The metabolome, the collection of small molecules associated with an organism, is a growing subject of inquiry, with the data utilized for data-intensive systems biology, disease diagnostics, biomarker discovery, and the broader characterization of small molecules in mixtures. Owing to their close proximity to the functional endpoints that govern an organism's phenotype, metabolites are highly informative about functional states. The field of metabolomics identifies and quantifies endogenous and exogenous metabolites in biological samples. Information acquired from nuclear magnetic spectroscopy (NMR), mass spectrometry (MS), and the published literature, as processed by statistical approaches, are driving increasingly wider applications of metabolomics. This review focuses on the role of databases and software tools in advancing the rigor, robustness, reproducibility, and validation of metabolomics studies.

Integrative NMR for biomolecular research.

NMR spectroscopy is a powerful technique for determining structural and functional features of biomolecules in physiological solution as well as for observing their intermolecular interactions in real-time. However, complex steps associated with its practice have made the approach daunting for non-specialists. We introduce an NMR platform that makes biomolecular NMR spectroscopy much more accessible by integrating tools, databases, web services, and video tutorials that can be launched by simple installation of NMRFAM software packages or using a cross-platform virtual machine that can be run on any standard laptop or desktop computer. The software package can be downloaded freely from the NMRFAM software download page ( http://pine.nmrfam.wisc.edu/download_packages.html ), and detailed instructions are available from the Integrative NMR Video Tutorial page ( http://pine.nmrfam.wisc.edu/integrative.html ).

NMRmix: A Tool for the Optimization of Compound Mixtures in 1D (1)H NMR Ligand Affinity Screens.

NMR ligand affinity screening is a powerful technique that is routinely used in drug discovery or functional genomics to directly detect protein-ligand binding events. Binding events can be identified by monitoring differences in the 1D (1)H NMR spectrum of a compound with and without protein. Although a single NMR spectrum can be collected within a short period (2-10 min per sample), one-by-one screening of a protein against a library of hundreds or thousands of compounds requires a large amount of spectrometer time and a large quantity of protein. Therefore, compounds are usually evaluated in mixtures ranging in size from 3 to 20 compounds to improve the efficiency of these screens in both time and material. Ideally, the NMR signals from individual compounds in the mixture should not overlap so that spectral changes can be associated with a particular compound. We have developed a software tool, NMRmix, to assist in creating ideal mixtures from a large panel of compounds with known chemical shifts. Input to NMRmix consists of an (1)H NMR peak list for each compound, a user-defined overlap threshold, and additional user-defined parameters if default settings are not used. NMRmix utilizes a simulated annealing algorithm to optimize the composition of the mixtures to minimize spectral peak overlaps so that each compound in the mixture is represented by a maximum number of nonoverlapping chemical shifts. A built-in graphical user interface simplifies data import and visual evaluation of the results.

Probabilistic validation of protein NMR chemical shift assignments.

Data validation plays an important role in ensuring the reliability and reproducibility of studies. NMR investigations of the functional properties, dynamics, chemical kinetics, and structures of proteins depend critically on the correctness of chemical shift assignments. We present a novel probabilistic method named ARECA for validating chemical shift assignments that relies on the nuclear Overhauser effect data . ARECA has been evaluated through its application to 26 case studies and has been shown to be complementary to, and usually more reliable than, approaches based on chemical shift databases. ARECA is available online at http://areca.nmrfam.wisc.edu/.

Solution Structural Studies of GTP:Adenosylcobinamide-Phosphateguanylyl Transferase (CobY) from Methanocaldococcus jannaschii.

GTP:adenosylcobinamide-phosphate (AdoCbi-P) guanylyl transferase (CobY) is an enzyme that transfers the GMP moiety of GTP to AdoCbi yielding AdoCbi-GDP in the late steps of the assembly of Ado-cobamides in archaea. The failure of repeated attempts to crystallize ligand-free (apo) CobY prompted us to explore its 3D structure by solution NMR spectroscopy. As reported here, the solution structure has a mixed α/ß fold consisting of seven ß-strands and five α-helices, which is very similar to a Rossmann fold. Titration of apo-CobY with GTP resulted in large changes in amide proton chemical shifts that indicated major structural perturbations upon complex formation. However, the CobY:GTP complex as followed by 1H-15N HSQC spectra was found to be unstable over time: GTP hydrolyzed and the protein converted slowly to a species with an NMR spectrum similar to that of apo-CobY. The variant CobYG153D, whose GTP complex was studied by X-ray crystallography, yielded NMR spectra similar to those of wild-type CobY in both its apo- state and in complex with GTP. The CobYG153D:GTP complex was also found to be unstable over time.