G protein inactive and active forms investigated by simulation methods

Maize and Common Bean Seed Exudates Mediate Part of Nonhost Resistance to Phytophthora sojae Prior to Infection

Phytophthora sojae does not infect nonhost maize (Zea mays) but infects nonhost common bean (Phaseolus vulgaris) under inoculation. Soybean seed exudates participate in mediating host resistance to P. sojae before infection. This study aims to elucidate the role of seed exudates in mediating the nonhost resistance of maize and common bean to P. sojae before infection. The behaviors of P. sojae zoospores in response to the seed exudates were determined using an assay chamber and a concave slide. The proteomes of P. sojae zoospores in response to the seed exudates were analyzed with the tandem mass tag method. The key proteins were quantitatively verified by parallel reaction monitoring. Maize seed exudates exerted a repellent effect on zoospores of P. sojae. This result explains why zoospores sense repelling signaling molecules in maize seed exudates that weaken and strongly inhibit chemotaxis signals in the phosphatidylinositol signaling pathway and arachidonic acid metabolism pathway. Common bean seed exudates did not exhibit any attraction to the zoospores because the guanine nucleotide-binding protein signaling pathway, which is responsible for transmitting chemotactic signals, had no significant change. The proteins protecting the cell membrane structure were significantly downregulated, and the early apoptosis signal glutathione was enhanced in zoospores responding to common bean seed exudates, which resulted in dissolution of the cysts. Maize and common bean seed exudates mediate part of the nonhost resistance to P. sojae via different mechanisms before infection. The immunity of maize to P. sojae is caused by the repellent effect of maize seed exudates on zoospores. Common bean seed exudates participate in mediating nonhost resistance by dissolving the cysts.

Genetics of Adaptation of the Ascomycetous Fungus Podospora anserina to Submerged Cultivation

Podospora anserina is a model ascomycetous fungus which shows pronounced phenotypic senescence when grown on solid medium but possesses unlimited lifespan under submerged cultivation. In order to study the genetic aspects of adaptation of P. anserina to submerged cultivation, we initiated a long-term evolution experiment. In the course of the first 4 years of the experiment, 125 single-nucleotide substitutions and 23 short indels were fixed in eight independently evolving populations. Six proteins that affect fungal growth and development evolved in more than one population; in particular, in the G-protein alpha subunit FadA, new alleles fixed in seven out of eight experimental populations, and these fixations affected just four amino acid sites, which is an unprecedented level of parallelism in experimental evolution. Parallel evolution at the level of genes and pathways, an excess of nonsense and missense substitutions, and an elevated conservation of proteins and their sites where the changes occurred suggest that many of the observed fixations were adaptive and driven by positive selection.

Invited review: Activation of G proteins by GTP and the mechanism of Gα-catalyzed GTP hydrolysis

This review addresses the regulatory consequences of the binding of GTP to the alpha subunits (Gα) of heterotrimeric G proteins, the reaction mechanism of GTP hydrolysis catalyzed by Gα and the means by which GTPase activating proteins (GAPs) stimulate the GTPase activity of Gα. The high energy of GTP binding is used to restrain and stabilize the conformation of the Gα switch segments, particularly switch II, to afford stable complementary to the surfaces of Gα effectors, while excluding interaction with Gβγ, the regulatory binding partner of GDP-bound Gα. Upon GTP hydrolysis, the energy of these conformational restraints is dissipated and the two switch segments, particularly switch II, become flexible and are able to adopt a conformation suitable for tight binding to Gβγ. Catalytic site pre-organization presents a significant activation energy barrier to Gα GTPase activity. The glutamine residue near the N-terminus of switch II (Glncat ) must adopt a conformation in which it orients and stabilizes the γ phosphate and the water nucleophile for an in-line attack. The transition state is probably loose with dissociative character; phosphoryl transfer may be concerted. The catalytic arginine in switch I (Argcat ), together with amide hydrogen bonds from the phosphate binding loop, stabilize charge at the β-γ bridge oxygen of the leaving group. GAPs that harbor "regulator of protein signaling" (RGS) domains, or structurally unrelated domains within G protein effectors that function as GAPs, accelerate catalysis by stabilizing the pre-transition state for Gα-catalyzed GTP hydrolysis, primarily by restraining Argcat and Glncat to their catalytic conformations. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 449-462, 2016.

Revealing Atomic-Level Mechanisms of Protein Allostery with Molecular Dynamics Simulations

Molecular dynamics (MD) simulations have become a powerful and popular method for the study of protein allostery, the widespread phenomenon in which a stimulus at one site on a protein influences the properties of another site on the protein. By capturing the motions of a protein's constituent atoms, simulations can enable the discovery of allosteric binding sites and the determination of the mechanistic basis for allostery. These results can provide a foundation for applications including rational drug design and protein engineering. Here, we provide an introduction to the investigation of protein allostery using molecular dynamics simulation. We emphasize the importance of designing simulations that include appropriate perturbations to the molecular system, such as the addition or removal of ligands or the application of mechanical force. We also demonstrate how the bidirectional nature of allostery-the fact that the two sites involved influence one another in a symmetrical manner-can facilitate such investigations. Through a series of case studies, we illustrate how these concepts have been used to reveal the structural basis for allostery in several proteins and protein complexes of biological and pharmaceutical interest.

SIGNAL TRANSDUCTION. Structural basis for nucleotide exchange in heterotrimeric G proteins

G protein-coupled receptors (GPCRs) relay diverse extracellular signals into cells by catalyzing nucleotide release from heterotrimeric G proteins, but the mechanism underlying this quintessential molecular signaling event has remained unclear. Here we use atomic-level simulations to elucidate the nucleotide-release mechanism. We find that the G protein α subunit Ras and helical domains-previously observed to separate widely upon receptor binding to expose the nucleotide-binding site-separate spontaneously and frequently even in the absence of a receptor. Domain separation is necessary but not sufficient for rapid nucleotide release. Rather, receptors catalyze nucleotide release by favoring an internal structural rearrangement of the Ras domain that weakens its nucleotide affinity. We use double electron-electron resonance spectroscopy and protein engineering to confirm predictions of our computationally determined mechanism.

G-protein αq participates in the steroid hormone 20-hydroxyecdysone nongenomic signal transduction

The nuclear receptor-mediated genomic pathways of the animal steroid hormones are well known. However, the cell membrane receptor-mediated nongenomic pathways of the animal steroid hormones are little understood. In this study, we report the participation of a G-protein alpha q (Gαq)(1) subunit in the 20E nongenomic pathway in the cell membrane and regulating gene expression during molting and metamorphosis in a lepidopteran insect, Helicoverpa armigera. 20E-induced phosphorylation of Gαq was detected using two-dimensional electrophoresis techniques. Knockdown of Gαq by injecting double-stranded RNA suppressed the development of larvae, delayed metamorphosis, and inhibited 20E-induced gene expression. Gαq was distributed throughout the cell, and migrated toward the plasma membrane upon 20E induction. Gαq was necessary in the 20E-induced intracellular Ca(2+) release and extracellular Ca(2+) influx. The protein kinase C (PKC) inhibitor could repress 20E-induced phosphorylation of cyclin-dependent kinase 10 (CDK10) and transcription factor ultraspiracle (USP1). PKC inhibitor could repress the Gαq phosphorylation and membrane trafficking. These results suggest that Gαq participates in 20E signaling in the cell membrane at the pre-genomic stage by modulating the increase of the intracellular Ca(2+) and phosphorylation of CDK10 and USP1 in 20E transcription complex to regulate gene transcription.

Domain-opening and dynamic coupling in the α-subunit of heterotrimeric G proteins

Heterotrimeric G proteins are conformational switches that turn on intracellular signaling cascades in response to the activation of G-protein-coupled receptors. Receptor activation by extracellular stimuli promotes a cycle of GTP binding and hydrolysis on the G protein α-subunit (Gα). Important conformational transitions occurring during this cycle have been characterized from extensive crystallographic studies of Gα. However, the link between the observed conformations and the mechanisms involved in G-protein activation and effector interaction remain unclear. Here we describe a comprehensive principal component analysis of available Gα crystallographic structures supplemented with extensive unbiased conventional and accelerated molecular dynamics simulations that together characterize the response of Gα to GTP binding and hydrolysis. Our studies reveal details of activating conformational changes as well as the intrinsic flexibility of the α-helical domain that includes a large-scale 60° domain opening under nucleotide-free conditions. This result is consistent with the recently reported open crystal structure of Gs, the stimulatory G protein for adenylyl cyclase, in complex with the α2 adrenergic receptor. Sets of unique interactions potentially important for the conformational transition are also identified. Moreover simulations reveal nucleotide-dependent dynamical couplings of distal regions and residues potentially important for the allosteric link between functional sites.

Conformational fluctuations of UreG, an intrinsically disordered enzyme

UreG proteins are small GTP binding (G) proteins that catalyze the hydrolysis of GTP necessary for the maturation of urease, a virulence factor in bacterial pathogenesis. UreG proteins are the first documented cases of intrinsically disordered enzymes. The comprehension of the dynamics of folding-unfolding events occurring in this protein could shed light on the enzymatic mechanism of UreG. Here, we used the recently developed replica exchange with solute tempering (REST2) computational methodology to explore the conformational space of UreG from Helicobacter pylori (HpUreG) and to identify its structural fluctuations. The same simulation and analysis protocol has been applied to HypB from Methanocaldococcus jannaschii (MjHypB), which is closely related to UreG in both sequence and function, even though it is not intrinsically disordered. A comparison of the two systems reveals that both HpUreG and MjHypB feature a substantial rigidity of the protein regions involved in catalysis, justifying its residual catalytic activity. On the other hand, HpUreG tends to unfold more than MjHypB in portions involved in protein-protein interactions with metallochaperones necessary for the formation of multiprotein complexes known to be involved in urease activation.

Pharmacology of nonsteroidal antiinflammatory drugs and opioids

Chronic pain affects up to 50 million Americans every day. Traditional treatment has included acetaminophen, nonsteroidal antiinflammatory drugs (NSAIDs), or opioids. The combination of NSAIDs and opioids can provide effective treatment for up to 90% of patients with chronic pain, but the NSAIDs have the potential for significant, even life-threatening side effects. Additionally, the nonselective cyclooxygenase inhibitors with 16,000 deaths per year in the United States might not be any safer. The opioids are great for short-term pain, but may need to be adjusted or changed frequently due to the development of tolerance. Understanding of the mechanism of opioids and NSAIDs has improved greatly over the past decade, but is still incomplete.

The use of G-protein coupled receptor models in lead optimization

With the emerging new crystal structures of G-protein coupled receptors (GPCRs), the number of reported in silico receptor models vastly increases every year. The use of these models in lead optimization (LO) is investigated here. Although there are many studies where GPCR models are used to identify new chemotypes by virtual screening, the classical application in LO is rarely reported. The reason for this may be that the quality of a model, which is appropriate for atomistic modeling, must be very high, and the biology of GPCR ligand-dependent signaling is still not fully understood. However, the few reported studies show that GPCR models can be used efficiently in LO for various problems, such as affinity optimization or tuning of physicochemical parameters.

A concerted mechanism for opening the GDP binding pocket and release of the nucleotide in hetero-trimeric G-proteins

G-protein hetero-trimers play a fundamental role in cell function. Their dynamic behavior at the atomic level remains to be understood. We have studied the Gi hetero-trimer through a combination of molecular dynamics simulations and normal mode analyses. We showed that these big proteins could undergo large-amplitude conformational changes, without any energy penalty and with an intrinsic dynamics centered on their GDP binding pocket. Among the computed collective motions, one of the modes (mode 17) was particularly able to significantly open both the base and the phosphate sides of the GDP binding pocket. This mode describing mainly a motion between the Ras-like and the helical domains of G(α) was in close agreement with some available X-ray data and with many other biochemical/biophysical observations including the kink of helix α5. The use of a new protocol, which allows extraction of the GDP ligand along the computed normal modes, supported that the exit of GDP was largely coupled to an opening motion along mode 17. We propose for the first time a "concerted mechanism" model in which the opening of the GDP pocket and the kink of the α5 helix occur concomitantly and favor GDP release from G(αβγ) complexes. This model is discussed in the context of the G-protein-coupled receptor/G-protein interaction close to the cell membrane.

Conformational flexibility and binding interactions of the G protein βγ heterodimer

Previous NMR experiments on unbound G protein βγ heterodimer suggested that particular residues in the binding interface are mobile on the nanosecond timescale. In this work we performed nanosecond-timescale molecular dynamics simulations to investigate conformational changes and dynamics of Gβγ in the presence of several binding partners: a high-affinity peptide (SIGK), phosducin, and the GDP-bound α subunit. In these simulations, the high mobility of GβW99 was reduced by SIGK, and it appeared that a tyrosine might stabilize GβW99 by hydrophobic or aromatic stacking interactions in addition to hydrogen bonds. Simulations of the phosducin-Gβγ complex showed that the mobility of GβW99 was restricted, consistent with inferences from NMR. However, large-scale conformational changes of Gβγ due to binding, which were hypothesized in the NMR study, were not observed in the simulations, most likely due to their short (nanosecond) duration. A pocket consisting of hydrophobic amino acids on Gα appears to restrict GβW99 mobility in the crystal structure of the Gαβγ? heterotrimer. The simulation trajectories are consistent with this idea. However, local conformational changes of residues GβW63, GβW211, GβW297, GβW332, and GβW339 were detected during the MD simulations. As expected, the magnitude of atomic fluctuations observed in simulations was greater for α than for the βγ subunits, suggesting that α has greater flexibility. These observations support the notion that to maintain the high mobility of GβW99 observed by solution NMR requires that the Gβ-α interface must open up on time scale longer than can be observed in nanosecond scale simulations.

Using entropy of drug and protein graphs to predict FDA drug-target network: theoretic-experimental study of MAO inhibitors and hemoglobin peptides from Fasciola hepatica

There are many drugs described with very different affinity to a large number of receptors. In this work, we selected Drug-Target pairs (DTPs/nDTPs) of drugs with high affinity/non-affinity for different targets like proteins. Quantitative Structure-Activity Relationships (QSAR) models become a very useful tool in this context to substantially reduce time and resources consuming experiments. Unfortunately, most QSAR models predict activity against only one protein. To solve this problem, we developed here a multi-target QSAR (mt-QSAR) classifier using the MARCH-INSIDE technique to calculate structural parameters of drug and target plus one Artificial Neuronal Network (ANN) to seek the model. The best ANN model found is a Multi-Layer Perceptron (MLP) with profile MLP 32:32-15-1:1. This MLP classifies correctly 623 out of 678 DTPs (Sensitivity = 91.89%) and 2995 out of 3234 nDTPs (Specificity = 92.61%), corresponding to training Accuracy = 92.48%. The validation of the model was carried out by means of external predicting series. The model classifies correctly 313 out of 338 DTPs (Sensitivity = 92.60%) and 1411 out of 1534 nDTP (Specificity = 91.98%) in validation series, corresponding to total Accuracy = 92.09% for validation series (Predictability). This model favorably compares with other LDA and ANN models developed in this work and Machine Learning classifiers published before to address the same problem in different aspects. These mt-QSARs offer also a good opportunity to construct drug-protein Complex Networks (CNs) that can be used to explore large and complex drug-protein receptors databases. Finally, we illustrated two practical uses of this model with two different experiments. In experiment 1, we report prediction, synthesis, characterization, and MAO-A and MAO-B pharmacological assay of 10 rasagiline derivatives promising for anti-Parkinson drug design. In experiment 2, we report sampling, parasite culture, SEC and 1DE sample preparation, MALDI-TOF MS and MS/MS analysis, MASCOT search, MM/MD 3D structure modeling, and QSAR prediction for different peptides of hemoglobin found in the proteome of the human parasite Fasciola hepatica; which is promising for anti-parasite drug targets discovery.

GoLoco motif proteins binding to Galpha(i1): insights from molecular simulations

Molecular dynamics simulations, computational alanine scanning and sequence analysis were used to investigate the structural properties of the Gα_i1/GoLoco peptide complex. Using these methodologies, binding of the GoLoco motif peptide to the Gα_i1 subunit was found to restrict the relative movement of the helical and catalytic domains in the Gα_i1 subunit, which is in agreement with a proposed mechanism of GDP dissociation inhibition by GoLoco motif proteins. In addition, the results provide further insights into the role of the “Switch IV” region located within the helical domain of Gα, the conformation of which might be important for interactions with various Gα partners.