Point mutations in Arf1 reveal cooperative effects of the N-terminal extension and myristate for GTPase-activating protein catalytic activity.

The ADP-ribosylation factors (Arfs) constitute a family of small GTPases within the Ras superfamily, with a distinguishing structural feature of a hypervariable N-terminal extension of the G domain modified with myristate. Arf proteins, including Arf1, have roles in membrane trafficking and cytoskeletal dynamics. While screening for Arf1:small molecule co-crystals, we serendipitously solved the crystal structure of the non-myristoylated engineered mutation [L8K]Arf1 in complex with a GDP analogue. Like wild-type (WT) non-myristoylated Arf1•GDP, we observed that [L8K]Arf1 exhibited an N-terminal helix that occludes the hydrophobic cavity that is occupied by the myristoyl group in the GDP-bound state of the native protein. However, the helices were offset from one another due to the L8K mutation, with a significant change in position of the hinge region connecting the N-terminus to the G domain. Hypothesizing that the observed effects on behavior of the N-terminus affects interaction with regulatory proteins, we mutated two hydrophobic residues to examine the role of the N-terminal extension for interaction with guanine nucleotide exchange factors (GEFs) and GTPase Activating Proteins (GAPs. Different than previous studies, all mutations were examined in the context of myristoylated Arf. Mutations had little or no effect on spontaneous or GEF-catalyzed guanine nucleotide exchange but did affect interaction with GAPs. [F13A]myrArf1 was less than 1/2500, 1/1500, and 1/200 efficient as substrate for the GAPs ASAP1, ARAP1 and AGAP1; however, [L8A/F13A]myrArf1 was similar to WT myrArf1. Using molecular dynamics simulations, the effect of the mutations on forming alpha helices adjacent to a membrane surface was examined, yet no differences were detected. The results indicate that lipid modifications of GTPases and consequent anchoring to a membrane influences protein function beyond simple membrane localization. Hypothetical mechanisms are discussed.

Myr-Arf1 conformational flexibility at the membrane surface sheds light on the interactions with ArfGAP ASAP1.

ADP-ribosylation factor 1 (Arf1) interacts with multiple cellular partners and membranes to regulate intracellular traffic, organelle structure and actin dynamics. Defining the dynamic conformational landscape of Arf1 in its active form, when bound to the membrane, is of high functional relevance and key to understanding how Arf1 can alter diverse cellular processes. Through concerted application of nuclear magnetic resonance (NMR), neutron reflectometry (NR) and molecular dynamics (MD) simulations, we show that, while Arf1 is anchored to the membrane through its N-terminal myristoylated amphipathic helix, the G domain explores a large conformational space, existing in a dynamic equilibrium between membrane-associated and membrane-distal conformations. These configurational dynamics expose different interfaces for interaction with effectors. Interaction with the Pleckstrin homology domain of ASAP1, an Arf-GTPase activating protein (ArfGAP), restricts motions of the G domain to lock it in what seems to be a conformation exposing functionally relevant regions.

Exploring Sulfur Sites in Proteins via Triple-Resonance 1H-Detected 77Se NMR.

NMR spectroscopy has been applied to virtually all sites within proteins and biomolecules; however, the observation of sulfur sites remains very challenging. Recent studies have examined 77Se as a replacement for sulfur and applied 77Se NMR in both the solution and solid states. As a spin-1/2 nuclide, 77Se is attractive as a probe of sulfur sites, and it has a very large chemical shift range (due to a large chemical shift anisotropy), which makes it potentially very sensitive to structural and/or binding interactions as well as dynamics. Despite being a spin-1/2 nuclide, there have been rather limited studies of 77Se, and the ability to use 1H-indirect detection has been sparse. Some examples exist, but in the absence of a directly bonded, nonexchangeable 1H, these have been largely limited to smaller molecules. We develop and illustrate approaches using double-labeling of 13C and 77Se in proteins that enable more sensitive triple-resonance schemes via multistep coherence transfers and 1H-detection. These methods require specialized hardware and decoupling schemes, which we developed and will be discussed.

The small molecule inhibitor NAV-2729 has a complex target profile including multiple ADP-ribosylation factor regulatory proteins.

The ADP-ribosylation factor (Arf) GTPases and their regulatory proteins are implicated in cancer progression. NAV-2729 was previously identified as a specific inhibitor of Arf6 that reduced progression of uveal melanoma in an orthotopic xenograft. Here, our goal was to assess the inhibitory effects of NAV-2729 on the proliferation of additional cell types. We found NAV-2729 inhibited proliferation of multiple cell lines, but Arf6 expression did not correlate with NAV-2729 sensitivity, and knockdown of Arf6 affected neither cell viability nor sensitivity to NAV-2729. Furthermore, binding to native Arf6 was not detected; however, we determined that NAV-2729 inhibited both Arf exchange factors and Arf GTPase-activating proteins. ASAP1, a GTPase-activating protein linked to cancer progression, was further investigated. We demonstrated that NAV-2729 bound to the PH domain of ASAP1 and changed ASAP1 cellular distribution. However, ASAP1 knockdown did not fully recapitulate the cytoskeletal effects of NAV-2729 nor affect cell proliferation. Finally, our screens identified 48 other possible targets of NAV-2729. These results illustrate the complexities of defining targets of small molecules and identify NAV-2729 as a model PH domain-binding inhibitor.

A simple protocol for the production of highly deuterated proteins for biophysical studies.

Highly deuterated protein samples expand the biophysics and biological tool kit by providing, among other qualities, contrast matching in neutron diffraction experiments and reduction of dipolar spin interactions from normally protonated proteins in magnetic resonance studies, impacting both electron paramagnetic resonance and NMR spectroscopy. In NMR applications, deuteration is often combined with other isotopic labeling patterns to expand the range of conventional NMR spectroscopy research in both solution and solid-state conditions. However, preparation of deuterated proteins is challenging. We present here a simple, effective, and user-friendly protocol to produce highly deuterated proteins in Escherichia coli cells. The protocol utilizes the common shaker flask growth method and the well-known pET system (which provides expression control via the T7 promotor) for large-scale recombinant protein expression. One liter expression typically yields 5 to 50 mg of highly deuterated protein. Our data demonstrate that the optimized procedure produces a comparable quantity of protein in deuterium (2H2O) oxide M9 medium compared with that in 1H2O M9 medium. The protocol will enable a broader utilization of deuterated proteins in a number of biophysical techniques.

Insights into the Cross Talk between Effector and Allosteric Lobes of KRAS from Methyl Conformational Dynamics.

KRAS is the most frequently mutated RAS protein in cancer patients, and it is estimated that about 20% of the cancer patients in the United States carried mutant RAS proteins. To accelerate therapeutic development, structures and dynamics of RAS proteins had been extensively studied by various biophysical techniques for decades. Although 31P NMR studies revealed population equilibrium of the two major states in the active GMPPNP-bound form, more complex conformational dynamics in RAS proteins and oncogenic mutants subtly modulate the interactions with their downstream effectors. We established a set of customized NMR relaxation dispersion techniques to efficiently and systematically examine the ms-µs conformational dynamics of RAS proteins. This method allowed us to observe varying synchronized motions that connect the effector and allosteric lobes in KRAS. We demonstrated the role of conformational dynamics of KRAS in controlling its interaction with the Ras-binding domain of the downstream effector RAF1, the first kinase in the MAPK pathway. This allows one to explain, as well as to predict, the altered binding affinities of various KRAS mutants, which was neither previously reported nor apparent from the structural perspective.

nightshift: A Python program for plotting simulated NMR spectra from assigned chemical shifts from the Biological Magnetic Resonance Data Bank.

Nuclear magnetic resonance (NMR) provides site specific information on local environments through chemical shifts. NMR is widely used in the study of proteins, ranging from determination of three-dimensional (3D) structures to characterizing dynamics and binding of small molecules and other proteins or ligands. Assigned chemical shift data for the atoms within proteins is a treasure trove of information that can facilitate a broad range of biochemical and biophysical studies. The Biological Magnetic Resonance Data Bank (BMRB) is a publicly accessible database that contains a large number of assigned chemical shifts; however, translating this wealth of knowledge into a practical application is not straightforward. Herein we present nightshift: a Python command line utility and library for plotting simulated two-dimensional (2D) and 3D NMR spectra from assigned chemical shifts in the BMRB. This tool allows users to simulate routinely collected amide and methyl fingerprint spectra, backbone triple-resonance assignment spectra, and user-defined custom correlations, including ones that do not necessarily correspond to published experiments. This tool enables experienced NMR spectroscopists, those learning the craft, and interested scientists seeking to utilize NMR the ability to preview or examine a wide range of spectra for proteins whose assignments are deposited in the BMRB, irrespective of whether those experiments have been executed or reported. The tool applies equally to folded and intrinsically disordered proteins, limited only by the existence of a BMRB deposition. The features of nightshift are described along with applications that illustrate the ease with which complicated correlation spectra and binding events can be simulated.

A structurally conserved site in AUP1 binds the E2 enzyme UBE2G2 and is essential for ER-associated degradation.

Endoplasmic reticulum-associated degradation (ERAD) is a protein quality control pathway of fundamental importance to cellular homeostasis. Although multiple ERAD pathways exist for targeting topologically distinct substrates, all pathways require substrate ubiquitination. Here, we characterize a key role for the UBE2G2 Binding Region (G2BR) of the ERAD accessory protein ancient ubiquitous protein 1 (AUP1) in ERAD pathways. This 27-amino acid (aa) region of AUP1 binds with high specificity and low nanomolar affinity to the backside of the ERAD ubiquitin-conjugating enzyme (E2) UBE2G2. The structure of the AUP1 G2BR (G2BRAUP1) in complex with UBE2G2 reveals an interface that includes a network of salt bridges, hydrogen bonds, and hydrophobic interactions essential for AUP1 function in cells. The G2BRAUP1 shares significant structural conservation with the G2BR found in the E3 ubiquitin ligase gp78 and in vitro can similarly allosterically activate ubiquitination in conjunction with ERAD E3s. In cells, AUP1 is uniquely required to maintain normal levels of UBE2G2; this is due to G2BRAUP1 binding to the E2 and preventing its rapid degradation. In addition, the G2BRAUP1 is required for both ER membrane recruitment of UBE2G2 and for its activation at the ER membrane. Thus, by binding to the backside of a critical ERAD E2, G2BRAUP1 plays multiple critical roles in ERAD.

Theoretical classification of exchange geometries from the perspective of NMR relaxation dispersion.

NMR relaxation dispersion experiments have been widely applied to probe important conformational exchange of macro-molecules in many biological systems. The current improvements in computational techniques as well as the theoretical breakthroughs make the quantitative data analysis of complex exchange models possible. However, the topology of a given exchange model is also one of the main factors affecting the solution of Bloch-McConnell equation. The lack of a theoretical analysis of the exchange topologies at n-site exchange hinders further progress of such data analysis. Here, using graph theory, we reveal the topological complexity of n-site exchange and present all exchange models when n is less than 6. Furthermore, we introduce an alternative way, using machine learning, to select an exchange model based on a set of relaxation dispersion data without fitting them with every individual exchange model.

Combined multi-band decoupling in biomolecular NMR spectroscopy.

Multi-resonance NMR experiments are powerful analytical and structural tools. Their conceptualization assumes that RF fields may be combined independently to manipulate spin interactions. However, practical implementation can compromise performance. One limitation is the generation of combination bands when two or more RF fields are applied simultaneously within the NMR probe. The combination bands can lead to significant interference with the detection circuitry. A facile approach to combined multi-band decoupling can resolve these problems and increase sensitivity two-fold (or more), by time sharing the application of the individual frequencies rather than time sharing decoupling and data acquisition.

Protein Dynamics revealed by NMR Relaxation Methods.

Structural biology often focuses primarily on three-dimensional structures of biological macromolecules, deposited in the Protein Data Bank (PDB). This resource is a remarkable entity for the world-wide scientific and medical communities, as well as the general public, as it is a growing translation into three-dimensional space of the vast information in genomic databases, e.g. GENBANK. There is, however, significantly more to understanding biological function than the three-dimensional coordinate space for ground-state structures of biomolecules. The vast array of biomolecules experiences natural dynamics, interconversion between multiple conformational states, and molecular recognition and allosteric events that play out on timescales ranging from picoseconds to seconds. This wide range of timescales demands ingenious and sophisticated experimental tools to sample and interpret these motions, thus enabling clearer insight into functional annotation of the PDB. NMR spectroscopy is unique in its ability to sample this range of timescales at atomic resolution and in physiologically relevant conditions using spin relaxation methods. The field is constantly expanding to provide new creative experiments, to yield more detailed coverage of timescales, and to broaden the power of interpretation and analysis methods. This review highlights the current state of the methodology and examines the extension of analysis tools for more complex experiments and dynamic models. The future for understanding protein dynamics is bright, and these extended tools bring greater compatibility with developments in computational molecular dynamics, all of which will further our understanding of biological molecular functions. These facets place NMR as a key component in integrated structural biology.

Membrane surface recognition by the ASAP1 PH domain and consequences for interactions with the small GTPase Arf1.

Adenosine diphosphate-ribosylation factor (Arf) guanosine triphosphatase-activating proteins (GAPs) are enzymes that need to bind to membranes to catalyze the hydrolysis of guanosine triphosphate (GTP) bound to the small GTP-binding protein Arf. Binding of the pleckstrin homology (PH) domain of the ArfGAP With SH3 domain, ankyrin repeat and PH domain 1 (ASAP1) to membranes containing phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2] is key for maximum GTP hydrolysis but not fully understood. By combining nuclear magnetic resonance, neutron reflectometry, and molecular dynamics simulation, we show that binding of multiple PI(4,5)P2 molecules to the ASAP1 PH domain (i) triggers a functionally relevant allosteric conformational switch and (ii) maintains the PH domain in a well-defined orientation, allowing critical contacts with an Arf1 mimic to occur. Our model provides a framework to understand how binding of the ASAP1 PH domain to PI(4,5)P2 at the membrane may play a role in the regulation of ASAP1.

Achieving pure spin effects by artifact suppression in methyl adiabatic relaxation experiments.

Recent methyl adiabatic relaxation dispersion experiments provide examination of conformational dynamics across a very wide timescale (102-105 s-1) and, particularly, provide insight into the hydrophobic core of proteins and allosteric effects associated with modulators. The experiments require efficient decoupling of 1H and 13C spin interactions, and some artifacts have been discovered, which are associated with the design of the proton decoupling scheme. The experimental data suggest that the original design is valid; however, pulse sequences with either no proton decoupling or proton decoupling with imperfect pulses can potentially exhibit complications in the experiments. Here, we demonstrate that pulse imperfections in the proton decoupling scheme can be dramatically alleviated by using a single composite π pulse and provide pure single-exponential relaxation data. It allows the opportunity to access high-quality methyl adiabatic relaxation dispersion data by removing the cross-correlation between dipole-dipole interaction and chemical shift anisotropy. The resulting high-quality data is illustrated with the binding of an allosteric modulator (G2BR) to the ubiquitin conjugating enzyme Ube2g2.

Optimization of sortase A ligation for flexible engineering of complex protein systems.

Engineering and bioconjugation of proteins is a critically valuable tool that can facilitate a wide range of biophysical and structural studies. The ability to orthogonally tag or label a domain within a multidomain protein may be complicated by undesirable side reactions to noninvolved domains. Furthermore, the advantages of segmental (or domain-specific) isotopic labeling for NMR, or deuteration for neutron scattering or diffraction, can be realized by an efficient ligation procedure. Common methods-expressed protein ligation, protein trans-splicing, and native chemical ligation-each have specific limitations. Here, we evaluated the use of different variants of Staphylococcus aureus sortase A for a range of ligation reactions and demonstrate that conditions can readily be optimized to yield high efficiency (i.e. completeness of ligation), ease of purification, and functionality in detergents. These properties may enable joining of single domains into multidomain proteins, lipidation to mimic posttranslational modifications, and formation of cyclic proteins to aid in the development of nanodisc membrane mimetics. We anticipate that the method for ligating separate domains into a single functional multidomain protein reported here may enable many applications in structural biology.

Interaction of the N terminus of ADP-ribosylation factor with the PH domain of the GTPase-activating protein ASAP1 requires phosphatidylinositol 4,5-bisphosphate.

Arf GAP with Src homology 3 domain, ankyrin repeat, and pleckstrin homology (PH) domain 1 (ASAP1) is a multidomain GTPase-activating protein (GAP) for ADP-ribosylation factor (ARF)-type GTPases. ASAP1 affects integrin adhesions, the actin cytoskeleton, and invasion and metastasis of cancer cells. ASAP1's cellular function depends on its highly-regulated and robust ARF GAP activity, requiring both the PH and the ARF GAP domains of ASAP1, and is modulated by phosphatidylinositol 4,5-bisphosphate (PIP2). The mechanistic basis of PIP2-stimulated GAP activity is incompletely understood. Here, we investigated whether PIP2 controls binding of the N-terminal extension of ARF1 to ASAP1's PH domain and thereby regulates its GAP activity. Using [Δ17]ARF1, lacking the N terminus, we found that PIP2 has little effect on ASAP1's activity. A soluble PIP2 analog, dioctanoyl-PIP2 (diC8PIP2), stimulated GAP activity on an N terminus-containing variant, [L8K]ARF1, but only marginally affected activity on [Δ17]ARF1. A peptide comprising residues 2-17 of ARF1 ([2-17]ARF1) inhibited GAP activity, and PIP2-dependently bound to a protein containing the PH domain and a 17-amino acid-long interdomain linker immediately N-terminal to the first ß-strand of the PH domain. Point mutations in either the linker or the C-terminal α-helix of the PH domain decreased [2-17]ARF1 binding and GAP activity. Mutations that reduced ARF1 N-terminal binding to the PH domain also reduced the effect of ASAP1 on cellular actin remodeling. Mutations in the ARF N terminus that reduced binding also reduced GAP activity. We conclude that PIP2 regulates binding of ASAP1's PH domain to the ARF1 N terminus, which may partially regulate GAP activity.

Probing the Broad Time Scale and Heterogeneous Conformational Dynamics in the Catalytic Core of the Arf-GAP ASAP1 via Methyl Adiabatic Relaxation Dispersion.

Methyl-TROSY is one of the most powerful NMR spectroscopic tools for studying structures and conformational dynamics of large protein complexes in solution. In studying conformational dynamics, side chains usually display heterogeneous dynamics, including collective and local motions, that can be difficult to detect and analyze by conventional relaxation dispersion (RD) approaches. The combination of NH-based heteronuclear adiabatic relaxation dispersion (HARD) experiments and a geometric approximation (geoHARD) has been shown to have several advantages over conventional RD in revealing conformational dynamics over a broad time scale. Here, we demonstrate a new technique that has been developed to detect both heterogeneous and wide time scale conformational dynamics in the hydrophobic interior of large macromolecules utilizing methyl-geoHARD. It is shown that methyl-geoHARD will be feasible at ultrahigh magnetic fields (>1 GHz), when this technology becomes available. For the ZA domain of Arf-GAP ASAP1, with a global correlational time of 24 ns at 15 °C, a wide range of conformational dynamics (exhibiting chemical exchange rates (kex) between 102 and 105 s-1) are observed in the methyl groups of isoleucine, leucine, and valine. The dynamics include collective and independent local motions. Furthermore, portions of the collective motions have been confirmed by single-quantum Carr-Purcell-Meiboom-Gill (SQ-CPMG) RD experiments; however, motions outside of the detectable CPMG window (400-8000 s-1) cannot be accurately determined by SQ-CPMG experiments. The methyl-geoHARD experiment allows the dissection of heterogeneous conformational dynamics and pinpoints important motions that, potentially, can be correlated with important biological functions and recognition.

Functional Expression and Characterization of Human Myristoylated-Arf1 in Nanodisc Membrane Mimetics.

Lipidated small GTP-binding proteins of the Arf family interact with multiple cellular partners and with membranes to regulate intracellular traffic and organelle structure. Here, we focus on the ADP-ribosylation factor 1 (Arf1), which interacts with numerous proteins in the Arf pathway, such as the ArfGAP ASAP1 that is highly expressed and activated in several cancer cell lines and associated with enhanced migration, invasiveness, and poor prognosis. Understanding the molecular and mechanistic details of Arf1 regulation at the membrane via structural and biophysical studies requires large quantities of fully functional protein bound to lipid bilayers. Here, we report on the production of a functional human Arf1 membrane platform on nanodiscs for biophysical studies. Large scale bacterial production of highly pure, N-myristoylated human Arf1 has been achieved, including complex isotopic labeling for nuclear magnetic resonance (NMR) studies, and the myr-Arf1 can be readily assembled in small nanoscale lipid bilayers (nanodiscs, NDs). It is determined that myr-Arf1 requires a minimum binding surface in the NDs of ∼20 lipids. Fluorescence and NMR were used to establish nucleotide exchange and ArfGAP-stimulated GTP hydrolysis at the membrane, indicating that phophoinositide stimulation of the activity of the ArfGAP ASAP1 is ≥2000-fold. Differences in nonhydrolyzable GTP analogues are observed, and GMPPCP is found to be the most stable. Combined, these observations establish a functional environment for biophysical studies of Arf1 effectors and interactions at the membrane.

Conformational Dynamics and Allostery in E2:E3 Interactions Drive Ubiquitination: gp78 and Ube2g2.

Conformational dynamics plays a fundamental role in molecular recognition and activity in enzymes. The ubiquitin-conjugating enzyme (E2) Ube2g2 functions with the ubiquitin ligase (E3) gp78 to assemble poly-ubiquitin chains on target substrates. Two domains in gp78, RING and G2BR, bind to two distant regions of Ube2g2, and activate it for ubiquitin (Ub) transfer. G2BR increases the affinity between the RING and Ube2g2 by 50-fold, while the RING catalyzes the transfer of Ub from the Ube2g2∼Ub conjugate. How G2BR and RING activate Ube2g2 is unclear. In this work, conformational dynamics in Ube2g2 revealed a clear correlation of binding G2BR and RING with the sequential progression toward Ub transfer. The interrelationship of the existence and exchange between ground and excited states leads to a dynamic energy landscape model, in which redistribution of populations contributes to allostery and activation. These findings provide insight into gp78's modulation of conformational exchange in Ube2g2 to stimulate ubiquitination.

Application of geometric approximation to the CPMG experiment: Two- and three-site exchange.

The Carr-Purcell-Meiboom-Gill (CPMG) experiment is one of the most classical and well-known relaxation dispersion experiments in NMR spectroscopy, and it has been successfully applied to characterize biologically relevant conformational dynamics in many cases. Although the data analysis of the CPMG experiment for the 2-site exchange model can be facilitated by analytical solutions, the data analysis in a more complex exchange model generally requires computationally-intensive numerical analysis. Recently, a powerful computational strategy, geometric approximation, has been proposed to provide approximate numerical solutions for the adiabatic relaxation dispersion experiments where analytical solutions are neither available nor feasible. Here, we demonstrate the general potential of geometric approximation by providing a data analysis solution of the CPMG experiment for both the traditional 2-site model and a linear 3-site exchange model. The approximate numerical solution deviates less than 0.5% from the numerical solution on average, and the new approach is computationally 60,000-fold more efficient than the numerical approach. Moreover, we find that accurate dynamic parameters can be determined in most cases, and, for a range of experimental conditions, the relaxation can be assumed to follow mono-exponential decay. The method is general and applicable to any CPMG RD experiment (e.g. N, C', Cα, Hα, etc.) The approach forms a foundation of building solution surfaces to analyze the CPMG experiment for different models of 3-site exchange. Thus, the geometric approximation is a general strategy to analyze relaxation dispersion data in any system (biological or chemical) if the appropriate library can be built in a physically meaningful domain.