Catalytic site mutations confer multiple states of G protein activation.

Heterotrimeric guanine nucleotide-binding proteins (G proteins) that function as molecular switches for cellular growth and metabolism are activated by GTP and inactivated by GTP hydrolysis. In uveal melanoma, a conserved glutamine residue critical for GTP hydrolysis in the G protein α subunit is often mutated in Gαq or Gα11 to either leucine or proline. In contrast, other glutamine mutations or mutations in other Gα subtypes are rare. To uncover the mechanism of the genetic selection and the functional role of this glutamine residue, we analyzed all possible substitutions of this residue in multiple Gα isoforms. Through cell-based measurements of activity, we showed that some mutants were further activated and inactivated by G protein-coupled receptors. Through biochemical, molecular dynamics, and nuclear magnetic resonance-based structural studies, we showed that the Gα mutants were functionally distinct and conformationally diverse, despite their shared inability to hydrolyze GTP. Thus, the catalytic glutamine residue contributes to functions beyond GTP hydrolysis, and these functions include subtype-specific, allosteric modulation of receptor-mediated subunit dissociation. We conclude that G proteins do not function as simple on-off switches. Rather, signaling emerges from an ensemble of active states, a subset of which are favored in disease and may be uniquely responsive to receptor-directed ligands.

Dynamic spatiotemporal determinants modulate GPCR:G protein coupling selectivity and promiscuity.

Recent studies have shown that G protein coupled receptors (GPCRs) show selective and promiscuous coupling to different Gα protein subfamilies and yet the mechanisms of the range of coupling preferences remain unclear. Here, we use Molecular Dynamics (MD) simulations on ten GPCR:G protein complexes and show that the location (spatial) and duration (temporal) of intermolecular contacts at the GPCR:Gα protein interface play a critical role in how GPCRs selectively interact with G proteins. We identify that some GPCR:G protein interface contacts are common across Gα subfamilies and others specific to Gα subfamilies. Using large scale data analysis techniques on the MD simulation snapshots we derive a spatio-temporal code for contacts that confer G protein selective coupling and validated these contacts using G protein activation BRET assays. Our results demonstrate that promiscuous GPCRs show persistent sampling of the common contacts more than G protein specific contacts. These findings suggest that GPCRs maintain contact with G proteins through a common central interface, while the selectivity comes from G protein specific contacts at the periphery of the interface.

Sequence coevolution and structure stabilization modulate olfactory receptor expression.

Olfactory receptors (ORs) belong to class A G-protein coupled receptors (GPCRs) and are activated by a variety of odorants. To date, there is no three-dimensional structure of an OR available. One of the major bottlenecks in obtaining purified protein for structural studies of ORs is their poor expression in heterologous cells. To design mutants that enhance expression and thereby enable protein purification, we first identified computable physical properties that recapitulate OR and class A GPCR expression and further conducted an iterative computational prediction-experimental test cycle and generated human OR mutants that express as high as biogenic amine receptors for which structures have been solved. In the process of developing the computational method to recapitulate the expression of ORs in membranes, we identified properties, such as amino acid sequence coevolution, and the strength of the interactions between intracellular loop 1 (ICL1) and the helix 8 region of ORs, to enhance their heterologous expression. We identified mutations that are directly located in these regions as well as other mutations not located in these regions but allosterically strengthen the ICL1-helix 8 enhance expression. These mutants also showed functional responses to known odorants. This method to enhance heterologous expression of mammalian ORs will facilitate high-throughput "deorphanization" of ORs, and enable OR purification for biochemical and structural studies to understand odorant-OR interactions.

A universal allosteric mechanism for G protein activation.

G proteins play a central role in signal transduction and pharmacology. Signaling is initiated by cell-surface receptors, which promote guanosine triphosphate (GTP) binding and dissociation of Gα from the Gßγ subunits. Structural studies have revealed the molecular basis of subunit association with receptors, RGS proteins, and downstream effectors. In contrast, the mechanism of subunit dissociation is poorly understood. We use cell signaling assays, molecular dynamics (MD) simulations, and biochemistry and structural analyses to identify a conserved network of amino acids that dictates subunit release. In the presence of the terminal phosphate of GTP, a glycine forms a polar network with an arginine and glutamate, putting torsional strain on the subunit binding interface. This "G-R-E motif" secures GTP and, through an allosteric link, discharges the Gßγ dimer. Replacement of network residues prevents subunit dissociation regardless of agonist or GTP binding. These findings reveal the molecular basis of the final committed step of G protein activation.

How Do Branched Detergents Stabilize GPCRs in Micelles?

The structural and functional properties of G protein-coupled receptors (GPCRs) are often studied in a detergent micellar environment, but many GPCRs tend to denature or aggregate in short alkyl chain detergents. In our previous work [Lee, S., et al. (2016) J. Am. Chem. Soc. 138, 15425-15433], we showed that GPCRs in alkyl glucosides were highly dynamic, resulting in the penetration of detergent molecules between transmembrane α-helices, which is the initial step in receptor denaturation. Although this was not observed for GPCRs in dodecyl maltoside (DDM, also known as lauryl maltoside), even this detergent is not mild enough to preserve the integrity of many GPCRs during purification. Lauryl maltose neopentylglycol (LMNG) detergents have been found to have significant advantages for purifying GPCRs in a native state as they impart more stability to the receptor than DDM. To gain insights into how they stabilize GPCRs, we used atomistic molecular dynamics simulations of wild type adenosine A2A receptor (WT-A2AR), thermostabilized A2AR (tA2AR), and wild type ß2-adrenoceptor (ß2AR) in a variety of detergents (LMNG, DMNG, OGNG, and DDM). Analysis of molecular dynamics simulations of tA2AR in LMNG, DMNG, and OGNG showed that this series of detergents exhibited behavior very similar to that of an analogous series of detergents DDM, DM, and OG in our previous study. However, there was a striking difference upon comparison of the behavior of LMNG to that of DDM. LMNG showed considerably less motion than DDM, which resulted in the enhanced density of the aliphatic chains around the hydrophobic regions of the receptor and considerably more hydrogen bond formation between the head groups. This contributed to enhanced interaction energies between both detergent molecules and between the receptor and detergent, explaining the enhanced stability of GPCRs purified in this detergent. Branched detergents occlude between transmembrane helices and reduce their flexibility. Our results provide a rational foundation to develop detergent variants for stabilizing membrane proteins.

Structural instability and divergence from conserved residues underlie intracellular retention of mammalian odorant receptors.

Mammalian odorant receptors are a diverse and rapidly evolving set of G protein-coupled receptors expressed in olfactory cilia membranes. Most odorant receptors show little to no cell surface expression in nonolfactory cells due to endoplasmic reticulum retention, which has slowed down biochemical studies. Here we provide evidence that structural instability and divergence from conserved residues of individual odorant receptors underlie intracellular retention using a combination of large-scale screening of odorant receptors cell surface expression in heterologous cells, point mutations, structural modeling, and machine learning techniques. We demonstrate the importance of conserved residues by synthesizing consensus odorant receptors that show high levels of cell surface expression similar to conventional G protein-coupled receptors. Furthermore, we associate in silico structural instability with poor cell surface expression using molecular dynamics simulations. We propose an enhanced evolutionary capacitance of olfactory sensory neurons that enable the functional expression of odorant receptors with cryptic mutations.

Machine Learning for Prioritization of Thermostabilizing Mutations for G-Protein Coupled Receptors.

Although the three-dimensional structures of G-protein coupled receptors (GPCRs), the largest superfamily of drug targets, have enabled structure-based drug design, there are no structures available for 87% of GPCRs. This is due to the stiff challenge in purifying the inherently flexible GPCRs. Identifying thermostabilized mutant GPCRs via systematic alanine scanning mutations has been a successful strategy in stabilizing GPCRs, but it remains a daunting task for each GPCR. We developed a computational method that combines sequence-, structure-, and dynamics-based molecular properties of GPCRs that recapitulate GPCR stability, with four different machine learning methods to predict thermostable mutations ahead of experiments. This method has been trained on thermostability data for 1231 mutants, the largest publicly available data set. A blind prediction for thermostable mutations of the complement factor C5a receptor 1 retrieved 36% of the thermostable mutants in the top 50 prioritized mutants compared to 3% in the first 50 attempts using systematic alanine scanning.

Prediction of Conformation Specific Thermostabilizing Mutations for Class A G Protein-Coupled Receptors.

G protein-coupled receptors (GPCRs) are highly flexible and prone to denaturation during protein extraction in detergents and purification. This poses a huge challenge to purify a conformationally homogeneous solution of GPCRs. Thermostabilizing mutations have been used widely to purify and obtain crystal structures of several GPCRs. However, identifying thermostabilizing mutations for GPCRs remains a tedious and expensive task as they are not transferable even among closely related GPCRs. Additionally, the mutations stabilizing one conformational state of a GPCR do not always stabilize other conformational state(s) of the same GPCR. Previously we developed a computational method, LiticonDesign, for rapid prediction of thermostabilizing mutations for a specific GPCR conformation. In this study, we have used LiticonDesign to predict thermostabilizing mutations for the agonist bound active-intermediate state of the human adenosine receptor (A2AR) using the structure of the inactive state of the same GPCR and vice versa. Our study shows that the thermostable mutation predictions using LiticonDesign, for an active-intermediate state of a GPCR (A2AR in our case), requires a homology model that is derived from an active/active-intermediate state GPCR structure as a template. Similarly, the homology models derived from inactive state GPCR conformations are better in predicting the thermostable mutations for the inactive state of A2AR. Overall, LiticonDesign method is not only efficient in predicting thermostabilizing mutations for a given GPCR sequence but also can recover conformation specific mutations for a state of interest, if a suitable starting structure of desired conformation is chosen.

Variable G protein determinants of GPCR coupling selectivity.

G protein-coupled receptors (GPCRs) activate four families of heterotrimeric G proteins, and individual receptors must select a subset of G proteins to produce appropriate cellular responses. Although the precise mechanisms of coupling selectivity are uncertain, the Gα subunit C terminus is widely believed to be the primary determinant recognized by cognate receptors. Here, we directly assess coupling between 14 representative GPCRs and 16 Gα subunits, including one wild-type Gα subunit from each of the four families and 12 chimeras with exchanged C termini. We use a sensitive bioluminescence resonance energy transfer (BRET) assay that provides control over both ligand and nucleotide binding, and allows direct comparison across G protein families. We find that the Gs- and Gq-coupled receptors we studied are relatively promiscuous and always couple to some extent to Gi1 heterotrimers. In contrast, Gi-coupled receptors are more selective. Our results with Gα subunit chimeras show that the Gα C terminus is important for coupling selectivity, but no more so than the Gα subunit core. The relative importance of the Gα subunit core and C terminus is highly variable and, for some receptors, the Gα core is more important for selective coupling than the C terminus. Our results suggest general rules for GPCR-G protein coupling and demonstrate that the critical G protein determinants of selectivity vary widely, even for different receptors that couple to the same G protein.

Azadirachtin inhibits amyloid formation, disaggregates pre-formed fibrils and protects pancreatic ß-cells from human islet amyloid polypeptide/amylin-induced cytotoxicity.

The human islet amyloid polypeptide (hIAPP) or amylin is the major constituent of amyloidogenic aggregates found in pancreatic islets of type 2 diabetic patients that have been associated with ß-cell dysfunction and/or death associated with type 2 diabetes mellitus (T2DM). Therefore, developing and/or identifying inhibitors of hIAPP aggregation pathway and/or compound that can mediate disaggregation of preformed aggregates holds promise as a medical intervention for T2DM management. In the current study, the anti-amyloidogenic potential of Azadirachtin (AZD)-a secondary metabolite isolated from traditional medicinal plant Neem (Azadirachta indica)-was investigated by using a combination of biophysical and cellular assays. Our results indicate that AZD supplementation not only inhibits hIAPP aggregation but also disaggregates pre-existing hIAPP fibrils by forming amorphous aggregates that are non-toxic to pancreatic ß-cells. Furthermore, AZD supplementation in pancreatic ß-cells (INS-1E) resulted in inhibition of oxidative stress; along with restoration of the DNA damage, lipid peroxidation and the associated membrane damage, endoplasmic reticulum stress and mitochondrial membrane potential. AZD treatment also restored glucose-stimulated insulin secretion from pancreatic islets exposed to hIAPP. All-atom molecular dynamics simulation studies on full-length hIAPP pentamer with AZD suggested that AZD interacted with four possible binding sites in the amyloidogenic region of hIAPP. In summary, our results suggest AZD to be a promising candidate for combating T2DM and related amyloidogenic disorders.

Salt Induced Structural Collapse, Swelling, and Signature of Aggregation of Two ssDNA Strands: Insights from Molecular Dynamics Simulation.

Molecular dynamics simulations elucidate the structural collapse shown by two ssDNAs of the same base sequence in the presence of either Na+ or Mg2+, starting from in vivo ionic concentration to higher concentrations. Initially, an increase in ion concentration facilitates the structural distortion of individual ssDNA and helps to bring them close, and for this, Mg2+ is better than Na+. However, further addition of ions leads to structural reswelling of the DNA strands and inhibits their proximity. The structural changes are found to be guided by the strong interaction of the cations with the phosphinyl oxygen (pn_O). Additionally, a significant difference has been noticed in the interaction of the cations with phosphoester oxygen (pe_O) depending on the nature of the ion. The sequential and nonsequential base-pair stacking is one of the major factors in the structural collapse of individual ssDNA. Overall, the present investigation highlights some of the important aspects of aggregation of two ssDNA with the same base sequence at varying cationic concentration.

Engineering Salt Bridge Networks between Transmembrane Helices Confers Thermostability in G-Protein-Coupled Receptors.

Introduction of specific point mutations has been an effective strategy in enhancing the thermostability of G-protein-coupled receptors (GPCRs). Our previous work showed that a specific residue position on transmembrane helix 6 (TM6) in class A GPCRs consistently yields thermostable mutants. The crystal structure of human chemokine receptor CCR5 also showed increased thermostability upon mutation of two positions, A233D6.33 and K303E7.59. With the goal of testing the transferability of these two thermostabilizing mutations in other chemokine receptors, we tested the mutations A237D6.33 and R307E7.59 in human CCR3 for thermostability and aggregation properties in detergent solution. Interestingly, the double mutant exhibited a 6-10-fold decrease in the aggregation propensity of the wild-type protein. This is in stark contrast to the two single mutants whose aggregation properties resemble the wild type (WT). Moreover, unlike in CCR5, the two single mutants separately showed no increase in thermostability compared to the wild-type CCR3, while the double-mutant A237D6.33/R307E7.59 confers an increase of 2.6 °C in the melting temperature compared to the WT. Extensive all-atom molecular dynamics (MD) simulations in detergent micelles show that a salt bridge network between transmembrane helices TM3, TM6, and TM7 that is absent in the two single mutants confers stability in the double mutant. The free energy surface of the double mutant shows conformational homogeneity compared to the single mutants. An annular n-dodecyl maltoside detergent layer packs tighter to the hydrophobic surface of the double-mutant CCR3 compared to the single mutants providing additional stability. The purification of other C-C chemokine receptors lacking such stabilizing residues may benefit from the incorporation of these two point mutations.

Can an ammonium-based room temperature ionic liquid counteract the urea-induced denaturation of a small peptide?

The folding/unfolding equilibrium of proteins in aqueous medium can be altered by adding small organic molecules generally termed as co-solvents. Denaturants such as urea are instrumental in the unfolding of proteins while protecting osmolytes favour the folded ensemble. Recently, room temperature ionic liquids (ILs) have been shown to counteract the deleterious effect of urea on proteins. In this paper, using atomistic molecular dynamics we show that a ternary mixture containing a particular ammonium-based IL, triethylammonium acetate (TEAA), and urea (in 1 : 5 molar ratio) helps a small 15-residue S-peptide analogue regain most of its native structure, whereas a binary aqueous mixture containing a large amount of urea alone completely distorts it. Our simulations show that the denaturant urea directly interacts with the peptide backbone in the binary mixture while for the ternary mixture both urea as well as the IL are preferentially excluded from the peptide surface.

Spontaneous Unzipping of Xylonucleic Acid Assisted by a Single-Walled Carbon Nanotube: A Computational Study.

Xenonucleic acids are synthetic nucleic acid analogues that are potential candidates for antisense or antigene therapy owing to their higher thermal and enzymatic stability compared to that of naturally occurring ones at physiological conditions. We investigate the binding and unzipping of xylonucleic acid (XNA) containing xylose (a stereoisomer of ribose) sugar in its backbone assisted by a single walled carbon nanotube (SWCNT) using extensive atomistic molecular dynamics simulations. Our simulations confirm XNA to undergo faster unzipping compared to a double stranded RNA with the same nucleobase sequence which is presumably due to the near orthogonal base pairing arrangement of the constituent nucleobases of XNA at physiologically relevant conditions (in terms of temperature and salt concentration). The interaction between XNA and SWCNT mimics that of a small interfering RNA (siRNA) and RNA-induced silencing complex (RISC) during a typical RNA induced gene silencing process where the duplex chain unwinding of the siRNA is of primordial importance. Our study may find relevance in designing a more efficient and safer delivery platform for xenonucleic acids by grafting these RNA analogues to SWCNT into an infected target cell. This unveils promising applications of XNA in the field of gene delivery for antisense therapies.

Probing the Salt Concentration Dependent Nucelobase Distribution in a Single-Stranded DNA-Single-Walled Carbon Nanotube Hybrid with Molecular Dynamics.

The hybrids of single-walled carbon nanotube (SWCNT) and single stranded DNA (ssDNA) are novel nanoscale materials having remarkable applications in nanotechnology. The absorption of nucleobases on the surface of a SWCNT depends strongly on the ionic strength of the medium. In this paper, using atomistic molecular dynamics we have shown that at low salt concentration ssDNA wraps on the surface of SWCNT through hydrophobic π-π stacking between the DNA bases and the sp(2)-hybridized carbon atoms of the carbon nanotube. At high salt concentration, however, the DNA molecule adopts a partially folded structure and the ssDNA-SWCNT wrapping gets weakened significantly due to the self-stacking of the DNA bases. Our study can find relevance in CNT mediated gene delivery processes where subsequent unwrapping of the gene from its carrier is anticipated across the cell membrane regulated by an existing salt concentration gradient.