Ligand efficacy modulates conformational dynamics of the µ-opioid receptor.

The µ-opioid receptor (µOR) is an important target for pain management1 and molecular understanding of drug action on µOR will facilitate the development of better therapeutics. Here we show, using double electron-electron resonance and single-molecule fluorescence resonance energy transfer, how ligand-specific conformational changes of µOR translate into a broad range of intrinsic efficacies at the transducer level. We identify several conformations of the cytoplasmic face of the receptor that interconvert on different timescales, including a pre-activated conformation that is capable of G-protein binding, and a fully activated conformation that markedly reduces GDP affinity within the ternary complex. Interaction of ß-arrestin-1 with the µOR core binding site appears less specific and occurs with much lower affinity than binding of Gi.

KARATE: PKA-induced KRAS4B-RHOA-mTORC2 supercomplex phosphorylates AKT in insulin signaling and glucose homeostasis.

AKT is a serine/threonine kinase that plays an important role in metabolism, cell growth, and cytoskeletal dynamics. AKT is activated by two kinases, PDK1 and mTORC2. Although the regulation of PDK1 is well understood, the mechanism that controls mTORC2 is unknown. Here, by investigating insulin receptor signaling in human cells and biochemical reconstitution, we found that insulin induces the activation of mTORC2 toward AKT by assembling a supercomplex with KRAS4B and RHOA GTPases, termed KARATE (KRAS4B-RHOA-mTORC2 Ensemble). Insulin-induced KARATE assembly is controlled via phosphorylation of GTP-bound KRAS4B at S181 and GDP-bound RHOA at S188 by protein kinase A. By developing a KARATE inhibitor, we demonstrate that KRAS4B-RHOA interaction drives KARATE formation. In adipocytes, KARATE controls insulin-dependent translocation of the glucose transporter GLUT4 to the plasma membrane for glucose uptake. Thus, our work reveals a fundamental mechanism that activates mTORC2 toward AKT in insulin-regulated glucose homeostasis.

Unlocking translational machinery for antitubercular drug development.

Targeting translational factor proteins (TFPs) presents significant promise for the development of innovative antitubercular drugs. Previous insights from antibiotic binding mechanisms and recently solved 3D crystal structures of Mycobacterium tuberculosis (Mtb) elongation factor thermo unstable-GDP (EF-Tu-GDP), elongation factor thermo stable-EF-Tu (EF-Ts-EF-Tu), and elongation factor G-GDP (EF-G-GDP) have opened up new avenues for the design and development of potent antituberculosis (anti-TB) therapies.

Structure-function relationships of the G domain, a canonical switch motif.

GTP-binding (G) proteins constitute a class of P-loop (phosphate-binding loop) proteins that work as molecular switches between the GDP-bound OFF and the GTP-bound ON state. The common principle is the 160-180-residue G domain with an α,ß topology that is responsible for nucleotide-dependent conformational changes and drives many biological functions. Although the G domain uses a universally conserved switching mechanism, its structure, function, and GTPase reaction are modified for many different pathways and processes.

Cryo-EM of elongating ribosome with EF-Tu•GTP elucidates tRNA proofreading.

Ribosomes accurately decode mRNA by proofreading each aminoacyl-tRNA that is delivered by the elongation factor EF-Tu1. To understand the molecular mechanism of this proofreading step it is necessary to visualize GTP-catalysed elongation, which has remained a challenge2-4. Here we use time-resolved cryogenic electron microscopy to reveal 33 ribosomal states after the delivery of aminoacyl-tRNA by EF-Tu•GTP. Instead of locking cognate tRNA upon initial recognition, the ribosomal decoding centre dynamically monitors codon-anticodon interactions before and after GTP hydrolysis. GTP hydrolysis enables the GTPase domain of EF-Tu to extend away, releasing EF-Tu from tRNA. The 30S subunit then locks cognate tRNA in the decoding centre and rotates, enabling the tRNA to bypass 50S protrusions during accommodation into the peptidyl transferase centre. By contrast, the decoding centre fails to lock near-cognate tRNA, enabling the dissociation of near-cognate tRNA both during initial selection (before GTP hydrolysis) and proofreading (after GTP hydrolysis). These findings reveal structural similarity between ribosomes in initial selection states5,6 and in proofreading states, which together govern the efficient rejection of incorrect tRNA.

Unveiling the catalytic mechanism of GTP hydrolysis in microtubules.

Microtubules (MTs) are large cytoskeletal polymers, composed of αß-tubulin heterodimers, capable of stochastically converting from polymerizing to depolymerizing states and vice versa. Depolymerization is coupled with hydrolysis of guanosine triphosphate (GTP) within ß-tubulin. Hydrolysis is favored in the MT lattice compared to a free heterodimer with an experimentally observed rate increase of 500- to 700-fold, corresponding to an energetic barrier lowering of 3.8 to 4.0 kcal/mol. Mutagenesis studies have implicated α-tubulin residues, α:E254 and α:D251, as catalytic residues completing the ß-tubulin active site of the lower heterodimer in the MT lattice. The mechanism for GTP hydrolysis in the free heterodimer, however, is not understood. Additionally, there has been debate concerning whether the GTP-state lattice is expanded or compacted relative to the GDP state and whether a "compacted" GDP-state lattice is required for hydrolysis. In this work, extensive quantum mechanics/molecular mechanics simulations with transition-tempered metadynamics free-energy sampling of compacted and expanded interdimer complexes, as well as a free heterodimer, have been carried out to provide clear insight into the GTP hydrolysis mechanism. α:E254 was found to be the catalytic residue in a compacted lattice, while in the expanded lattice, disruption of a key salt bridge interaction renders α:E254 less effective. The simulations reveal a barrier decrease of 3.8 ± 0.5 kcal/mol for the compacted lattice compared to a free heterodimer, in good agreement with experimental kinetic measurements. Additionally, the expanded lattice barrier was found to be 6.3 ± 0.5 kcal/mol higher than compacted, demonstrating that GTP hydrolysis is variable with lattice state and slower at the MT tip.

Structural basis of TRAPPIII-mediated Rab1 activation.

The GTPase Rab1 is a master regulator of the early secretory pathway and is critical for autophagy. Rab1 activation is controlled by its guanine nucleotide exchange factor, the multisubunit TRAPPIII complex. Here, we report the 3.7 Å cryo-EM structure of the Saccharomyces cerevisiae TRAPPIII complex bound to its substrate Rab1/Ypt1. The structure reveals the binding site for the Rab1/Ypt1 hypervariable domain, leading to a model for how the complex interacts with membranes during the activation reaction. We determined that stable membrane binding by the TRAPPIII complex is required for robust activation of Rab1/Ypt1 in vitro and in vivo, and is mediated by a conserved amphipathic α-helix within the regulatory Trs85 subunit. Our results show that the Trs85 subunit serves as a membrane anchor, via its amphipathic helix, for the entire TRAPPIII complex. These findings provide a structural understanding of Rab activation on organelle and vesicle membranes.

Cryo-EM structure of metazoan TRAPPIII, the multi-subunit complex that activates the GTPase Rab1.

The TRAPP complexes are nucleotide exchange factors that play essential roles in membrane traffic and autophagy. TRAPPII activates Rab11, and TRAPPIII activates Rab1, with the two complexes sharing a core of small subunits that affect nucleotide exchange but being distinguished by specific large subunits that are essential for activity in vivo. Crystal structures of core subunits have revealed the mechanism of Rab activation, but how the core and the large subunits assemble to form the complexes is unknown. We report a cryo-EM structure of the entire Drosophila TRAPPIII complex. The TRAPPIII-specific subunits TRAPPC8 and TRAPPC11 hold the catalytic core like a pair of tongs, with TRAPPC12 and TRAPPC13 positioned at the joint between them. TRAPPC2 and TRAPPC2L link the core to the two large arms, with the interfaces containing residues affected by disease-causing mutations. The TRAPPC8 arm is positioned such that it would contact Rab1 that is bound to the core, indicating how the arm could determine the specificity of the complex. A lower resolution structure of TRAPPII shows a similar architecture and suggests that the TRAPP complexes evolved from a single ur-TRAPP.

Effects of bound nucleotides on the secondary structure, thermal stability, and phosphorylation of Rab3A.

Rab3A is a member of the Rab GTPase family involved in synaptic vesicle trafficking. Recent evidence has demonstrated that Rab3A is phosphorylated by leucine-rich repeat kinase 2 (LRRK2) that is implicated in both familial and sporadic forms of Parkinson's disease (PD), and an abnormal increase in Rab3A phosphorylation has been proposed as a cause of PD. Despite the potential importance of Rab3A in PD pathogenesis, its structural information is limited and the effects of bound nucleotides on its biophysical and biochemical properties remain unclear. Here, we show that GDP-bound Rab3A is preferentially phosphorylated by LRRK2 compared with GTP-bound Rab3A. The secondary structure of Rab3A, measured by circular dichroism (CD) spectroscopy, revealed that Rab3A is resistant to heat-induced denaturation at pH 7.4 or 9.0 regardless of the nucleotides bound. In contrast, Rab3A underwent heat-induced denaturation at pH 5.0 at a lower temperature in its GDP-bound form than in its GTP-bound form. The unfolding temperature of Rab3A was studied by differential scanning fluorimetry, which showed a significantly higher unfolding temperature in GTP-bound Rab3A than in GDP-bound Rab3A, with the highest at pH 7.4. These results suggest that Rab3A has unusual thermal stability under physiologically relevant conditions and that bound nucleotides influence both thermal stability and phosphorylation by LRRK2.

Structural analysis of a trimeric assembly of the mitochondrial dynamin-like GTPase Mgm1.

The fusion of inner mitochondrial membranes requires dynamin-like GTPases, Mgm1 in yeast and OPA1 in mammals, but how they mediate membrane fusion is poorly understood. Here, we determined the crystal structure of Saccharomyces cerevisiae short Mgm1 (s-Mgm1) in complex with GDP. It revealed an N-terminal GTPase (G) domain followed by two helix bundles (HB1 and HB2) and a unique C-terminal lipid-interacting stalk (LIS). Dimers can form through antiparallel HB interactions. Head-to-tail trimers are built by intermolecular interactions between the G domain and HB2-LIS. Biochemical and in vivo analyses support the idea that the assembly interfaces observed here are native and critical for Mgm1 function. We also found that s-Mgm1 interacts with negatively charged lipids via both the G domain and LIS. Based on these observations, we propose that membrane targeting via the G domain and LIS facilitates the in cis assembly of Mgm1, potentially generating a highly curved membrane tip to allow inner membrane fusion.

Structural insights into the binding of nanobody Rh57 to active RhoA-GTP.

RhoA protein is a small GTPase that acts as a molecular switch. When bound to guanosine triphosphate (GTP), RhoA can activate several key signal pathways. Recently, nanobody Rh57 specific binding with GTP bound active RhoA was discovered and developed as a BRET biosensor without cytotoxicity. To further clarify the nanobody Rh57's mechanism of action, we co-expressed, purified, and crystallized the RhoA-Rh57 nanobody complex and solved the structure by X-ray diffraction with a resolution of 2.76 Å. The structure showed that the interaction is mainly through hydrogen bonds, salt bridges, aromatic-aromatic interactions, and hydrophobic interactions. The involved regions include CDR3 and non-hypervariable loop of Rh57, and the SWI switch loops of RhoA, respectively. The different SWI conformation of inactivated RhoA-GDP prevented the Rh57's binding. The possible explanation of Rh57 as a non-cytotoxic BRET intracellular tracer is that Rh57's binding did not overlap with downstream PRK1 and thus did not interfere with the downstream signaling pathway. Our research provides an in-depth understanding of how nanobodies recognize activated RhoA-GTP while not binding inactivated RhoA-GDP. This structural information may also provide critical information for further optimization of relevant nanobodies.

Markov State Models and Molecular Dynamics Simulations Reveal the Conformational Transition of the Intrinsically Disordered Hypervariable Region of K-Ras4B to the Ordered Conformation.

K-Ras4B, the most frequently mutated Ras isoform in human tumors, plays a vital part in cell growth, differentiation, and survival. Its tail, the C-terminal hypervariable region (HVR), is involved in anchoring K-Ras4B at the cellular plasma membrane and in isoform-specific protein-protein interactions and signaling. In the inactive guanosine diphosphate-bound state, the intrinsically disordered HVR interacts with the catalytic domain at the effector-binding region, rendering K-Ras4B in its autoinhibited state. Activation releases the HVR from the catalytic domain, with its ensemble favoring an ordered α-helical structure. The large-scale conformational transition of the HVR from the intrinsically disordered to the ordered conformation remains poorly understood. Here, we deploy a computational scheme that integrates a transition path-generation algorithm, extensive molecular dynamics simulation, and Markov state model analysis to investigate the conformational landscape of the HVR transition pathway. Our findings reveal a stepwise pathway for the HVR transition and uncover several key conformational substates along the transition pathway. Importantly, key interactions between the HVR and the catalytic domain are unraveled, highlighting the pathogenesis of K-Ras4B mild mutations in several congenital developmental anomaly syndromes. Together, these findings provide a deeper understanding of the HVR transition mechanism and the regulation of K-Ras4B activity at an atomic level.

Conformational switch that induces GDP release from Gi.

Heterotrimeric guanine nucleotide-binding proteins (G proteins) are composed of α, ß, and γ subunits. Gα switches between guanosine diphosphate (GDP)-bound inactive and guanosine triphosphate (GTP)-bound active states, and Gßγ interacts with the GDP-bound state. The GDP-binding regions are composed of two sites: the phosphate-binding and guanine-binding regions. The turnover of GDP and GTP is induced by guanine nucleotide-exchange factors (GEFs), including G protein-coupled receptors (GPCRs), Ric8A, and GIV/Girdin. However, the key structural factors for stabilizing the GDP-bound state of G proteins and the direct structural event for GDP release remain unclear. In this study, we investigated structural factors affecting GDP release by introducing point mutations in selected, conserved residues in Gαi3. We examined the effects of these mutations on the GDP/GTP turnover rate and the overall conformation of Gαi3 as well as the binding free energy between Gαi3 and GDP. We found that dynamic changes in the phosphate-binding regions are an immediate factor for the release of GDP.

Structural basis of substrate recognition and catalysis by fucosyltransferase 8.

Fucosylation of the innermost GlcNAc of N-glycans by fucosyltransferase 8 (FUT8) is an important step in the maturation of complex and hybrid N-glycans. This simple modification can dramatically affect the activities and half-lives of glycoproteins, effects that are relevant to understanding the invasiveness of some cancers, development of mAb therapeutics, and the etiology of a congenital glycosylation disorder. The acceptor substrate preferences of FUT8 are well-characterized and provide a framework for understanding N-glycan maturation in the Golgi; however, the structural basis of these substrate preferences and the mechanism through which catalysis is achieved remain unknown. Here we describe several structures of mouse and human FUT8 in the apo state and in complex with GDP, a mimic of the donor substrate, and with a glycopeptide acceptor substrate at 1.80-2.50 Å resolution. These structures provide insights into a unique conformational change associated with donor substrate binding, common strategies employed by fucosyltransferases to coordinate GDP, features that define acceptor substrate preferences, and a likely mechanism for enzyme catalysis. Together with molecular dynamics simulations, the structures also revealed how FUT8 dimerization plays an important role in defining the acceptor substrate-binding site. Collectively, this information significantly builds on our understanding of the core fucosylation process.

Ras-guanine nucleotide complexes: A UV spectral deconvolution method to analyze protein concentration, nucleotide stoichiometry, and purity.

The many members of the Ras superfamily are small GTPases that serve as molecular switches. These proteins bind the guanine nucleotides GTP and GDP with picomolar affinities, thereby stabilizing on- and off-signaling states, respectively. Quantitative in vitro Ras studies require accurate determination of total protein, its fractional occupancy with guanine nucleotide, and spectroscopic purity. Yet the high nucleotide affinity of Ras and the overlapping UV spectra of the protein and bound nucleotide make such determinations challenging. Here we describe a generalizable UV spectral deconvolution method to analyze the total protein concentration, total nucleotide stoichiometry, and purity of Ras complexes.

Synergistic role of nucleotides and lipids for the self-assembly of Shs1 septin oligomers.

Budding yeast septins are essential for cell division and polarity. Septins assemble as palindromic linear octameric complexes. The function and ultra-structural organization of septins are finely governed by their molecular polymorphism. In particular, in budding yeast, the end subunit can stand either as Shs1 or Cdc11. We have dissected, here, for the first time, the behavior of the Shs1 protomer bound to membranes at nanometer resolution, in complex with the other septins. Using electron microscopy, we have shown that on membranes, Shs1 protomers self-assemble into rings, bundles, filaments or two-dimensional gauzes. Using a set of specific mutants we have demonstrated a synergistic role of both nucleotides and lipids for the organization and oligomerization of budding yeast septins. Besides, cryo-electron tomography assays show that vesicles are deformed by the interaction between Shs1 oligomers and lipids. The Shs1-Shs1 interface is stabilized by the presence of phosphoinositides, allowing the visualization of micrometric long filaments formed by Shs1 protomers. In addition, molecular modeling experiments have revealed a potential molecular mechanism regarding the selectivity of septin subunits for phosphoinositide lipids.

Functional Characterization of Cj1427, a Unique Ping-Pong Dehydrogenase Responsible for the Oxidation of GDP-d-glycero-α-d-manno-heptose in Campylobacter jejuni.

The capsular polysaccharides (CPS) of Campylobacter jejuni contain multiple heptose residues with variable stereochemical arrangements at C3-C6. The immediate precursor to all of these possible variations is currently believed to be GDP-d-glycero-α-d-manno-heptose. Oxidation of this substrate at C4 enables subsequent epimerization reactions at C3-C5 that can be coupled to the dehydration/reduction at C5/C6. However, the enzyme responsible for the critical oxidation of C4 from GDP-d-glycero-α-d-manno-heptose has remained elusive. The enzyme Cj1427 from C. jejuni NCTC 11168 was shown to catalyze the oxidation of GDP-d-glycero-α-d-manno-heptose to GDP-d-glycero-4-keto-α-d-lyxo-heptose in the presence of α-ketoglutarate using mass spectrometry and nuclear magnetic resonance spectroscopy. At pH 7.4, the apparent kcat is 0.6 s-1, with a value of kcat/Km of 1.0 × 104 M-1 s-1 for GDP-d-glycero-α-d-manno-heptose. α-Ketoglutarate is required to recycle the tightly bound NADH nucleotide in the active site of Cj1427, which does not dissociate from the enzyme during catalysis.

Parkinson's disease-associated mutations in the GTPase domain of LRRK2 impair its nucleotide-dependent conformational dynamics.

Mutation in leucine-rich repeat kinase 2 (LRRK2) is a common cause of familial Parkinson's disease (PD). Recently, we showed that a disease-associated mutation R1441H rendered the GTPase domain of LRRK2 catalytically less active and thereby trapping it in a more persistently "on" conformation. However, the mechanism involved and characteristics of this on conformation remained unknown. Here, we report that the Ras of complex protein (ROC) domain of LRRK2 exists in a dynamic dimer-monomer equilibrium that is oppositely driven by GDP and GTP binding. We also observed that the PD-associated mutations at residue 1441 impair this dynamic and shift the conformation of ROC to a GTP-bound-like monomeric conformation. Moreover, we show that residue Arg-1441 is critical for regulating the conformational dynamics of ROC. In summary, our results reveal that the PD-associated substitutions at Arg-1441 of LRRK2 alter monomer-dimer dynamics and thereby trap its GTPase domain in an activated state.

Recurrent De Novo Mutations Disturbing the GTP/GDP Binding Pocket of RAB11B Cause Intellectual Disability and a Distinctive Brain Phenotype.

The Rab GTPase family comprises ∼70 GTP-binding proteins, functioning in vesicle formation, transport and fusion. They are activated by a conformational change induced by GTP-binding, allowing interactions with downstream effectors. Here, we report five individuals with two recurrent de novo missense mutations in RAB11B; c.64G>A; p.Val22Met in three individuals and c.202G>A; p.Ala68Thr in two individuals. An overlapping neurodevelopmental phenotype, including severe intellectual disability with absent speech, epilepsy, and hypotonia was observed in all affected individuals. Additionally, visual problems, musculoskeletal abnormalities, and microcephaly were present in the majority of cases. Re-evaluation of brain MRI images of four individuals showed a shared distinct brain phenotype, consisting of abnormal white matter (severely decreased volume and abnormal signal), thin corpus callosum, cerebellar vermis hypoplasia, optic nerve hypoplasia and mild ventriculomegaly. To compare the effects of both variants with known inactive GDP- and active GTP-bound RAB11B mutants, we modeled the variants on the three-dimensional protein structure and performed subcellular localization studies. We predicted that both variants alter the GTP/GDP binding pocket and show that they both have localization patterns similar to inactive RAB11B. Evaluation of their influence on the affinity of RAB11B to a series of binary interactors, both effectors and guanine nucleotide exchange factors (GEFs), showed induction of RAB11B binding to the GEF SH3BP5, again similar to inactive RAB11B. In conclusion, we report two recurrent dominant mutations in RAB11B leading to a neurodevelopmental syndrome, likely caused by altered GDP/GTP binding that inactivate the protein and induce GEF binding and protein mislocalization.

Recognizing five molecular ligand-binding sites with similar chemical structure.

Accurate identification of ligand-binding sites and discovering the protein-ligand interaction mechanism are important for understanding proteins' functions and designing new drugs. Meanwhile, accurate computational prediction and mechanism research are two grand challenges in proteomics. In this article, ligand-binding residues of five ligands (ATP, ADP, GTP, GDP, and NAD) are predicted as a group, due to their similar chemical structures and close biological function relations. The data set of binding sites by five ligands (ATP, ADP, GTP, GDP, and NAD) are collated from Biolip database. Then, five features, containing increment of diversity value, matrix scoring value, auto-covariance, secondary structure information, and surface accessibility information are used in binding site predictions. The support vector machine (SVM) model is used with the five features to predict ligand-binding sites. Finally, prediction results are tested by fivefold cross validation. Accuracy (Acc) of five ligands (ATP, ADP, GTP, GDP, and NAD) achieves 77.4%, 71.2%, 82.1%, 82.9%, and 85.3%, respectively; and Matthew correlation coefficient (MCC) of the above five ligands achieves 0.549, 0.424, 0.643, 0.659, and 0.702, respectively. The research result shows that for ligands with similar chemical structures, microenvironment of their binding sites and their sensitivities to features are similar, while, differences of their ligand-binding properties exist at the same time. © 2019 Wiley Periodicals, Inc.