Molecular dynamics simulations of transducin: interdomain and front to back communication in activation and nucleotide exchange

Structure and dynamics of GPCR signaling complexes

G-protein-coupled receptors (GPCRs) relay numerous extracellular signals by triggering intracellular signaling through coupling with G proteins and arrestins. Recent breakthroughs in the structural determination of GPCRs and GPCR-transducer complexes represent important steps toward deciphering GPCR signal transduction at a molecular level. A full understanding of the molecular basis of GPCR-mediated signaling requires elucidation of the dynamics of receptors and their transducer complexes as well as their energy landscapes and conformational transition rates. Here, we summarize current insights into the structural plasticity of GPCR-G-protein and GPCR-arrestin complexes that underlies the regulation of the receptor's intracellular signaling profile.

Suppressing allostery in epitope mapping experiments using millisecond hydrogen / deuterium exchange mass spectrometry

Localization of the interface between the candidate antibody and its antigen target, commonly known as epitope mapping, is a critical component of the development of therapeutic monoclonal antibodies. With the recent availability of commercial automated systems, hydrogen / deuterium eXchange (HDX) is rapidly becoming the tool for mapping epitopes preferred by researchers in both industry and academia. However, this approach has a significant drawback in that it can be confounded by ‘allosteric’ structural and dynamic changes that result from the interaction, but occur far from the point(s) of contact. Here, we introduce a ‘kinetic’ millisecond HDX workflow that suppresses allosteric effects in epitope mapping experiments. The approach employs a previously introduced microfluidic apparatus that enables millisecond HDX labeling times with on-chip pepsin digestion and electrospray ionization. The ‘kinetic’ workflow also differs from conventional HDX-based epitope mapping in that the antibody is introduced to the antigen at the onset of HDX labeling. Using myoglobin / anti-myoglobin as a model system, we demonstrate that at short ‘kinetic’ workflow labeling times (i.e., 200 ms), the HDX signal is already fully developed at the ‘true’ epitope, but is still largely below the significance threshold at allosteric sites. Identification of the ‘true’ epitope is supported by computational docking predictions and allostery modeling using the rigidity transmission allostery algorithm.

Single-molecule analysis of ligand efficacy in β2AR-G-protein activation

G protein-coupled receptor (GPCR)-mediated signal transduction is central to human physiology and disease intervention, yet the molecular mechanisms responsible for ligand-dependent signaling responses remain poorly understood. In Class A GPCRs, receptor activation and G protein coupling entail outward movements of transmembrane segment 6 (TM6). Using single-molecule Fluorescence Resonance Energy Transfer (smFRET) imaging, we examine TM6 motions in the β₂ adrenergic receptor (β₂AR) upon exposure to orthosteric ligands with different efficacies, in the absence and presence of the G_s heterotrimer. We show that partial and full agonists affect TM6 motions in a manner that differentially regulates the rate at which GDP-bound β₂AR-G_s complexes are formed and the efficiency of nucleotide exchange leading to G_s activation. These data also reveal transient nucleotide-bound β₂AR-G_s species distinct from known structures and single-molecule perspectives on the allosteric link between ligand and nucleotide binding pockets that shed new light on the G protein activation mechanism.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

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.

A Conserved Hydrophobic Core in Gαi1 Regulates G Protein Activation and Release from Activated Receptor

G protein-coupled receptor-mediated heterotrimeric G protein activation is a major mode of signal transduction in the cell. Previously, we and other groups reported that the α5 helix of Gαi1, especially the hydrophobic interactions in this region, plays a key role during nucleotide release and G protein activation. To further investigate the effect of this hydrophobic core, we disrupted it in Gαi1 by inserting 4 alanine amino acids into the α5 helix between residues Gln(333) and Phe(334) (Ins4A). This extends the length of the α5 helix without disturbing the β6-α5 loop interactions. This mutant has high basal nucleotide exchange activity yet no receptor-mediated activation of nucleotide exchange. By using structural approaches, we show that this mutant loses critical hydrophobic interactions, leading to significant rearrangements of side chain residues His(57), Phe(189), Phe(191), and Phe(336); it also disturbs the rotation of the α5 helix and the π-π interaction between His(57) and Phe(189) In addition, the insertion mutant abolishes G protein release from the activated receptor after nucleotide binding. Our biochemical and computational data indicate that the interactions between α5, α1, and β2-β3 are not only vital for GDP release during G protein activation, but they are also necessary for proper GTP binding (or GDP rebinding). Thus, our studies suggest that this hydrophobic interface is critical for accurate rearrangement of the α5 helix for G protein release from the receptor after GTP binding.

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.

A conserved phenylalanine as a relay between the α5 helix and the GDP binding region of heterotrimeric Gi protein α subunit

G protein activation by G protein-coupled receptors is one of the critical steps for many cellular signal transduction pathways. Previously, we and other groups reported that the α5 helix in the G protein α subunit plays a major role during this activation process. However, the precise signaling pathway between the α5 helix and the guanosine diphosphate (GDP) binding pocket remains elusive. Here, using structural, biochemical, and computational techniques, we probed different residues around the α5 helix for their role in signaling. Our data showed that perturbing the Phe-336 residue disturbs hydrophobic interactions with the β2-β3 strands and α1 helix, leading to high basal nucleotide exchange. However, mutations in β strands β5 and β6 do not perturb G protein activation. We have highlighted critical residues that leverage Phe-336 as a relay. Conformational changes are transmitted starting from Phe-336 via β2-β3/α1 to Switch I and the phosphate binding loop, decreasing the stability of the GDP binding pocket and triggering nucleotide release. When the α1 and α5 helices were cross-linked, inhibiting the receptor-mediated displacement of the C-terminal α5 helix, mutation of Phe-336 still leads to high basal exchange rates. This suggests that unlike receptor-mediated activation, helix 5 rotation and translocation are not necessary for GDP release from the α subunit. Rather, destabilization of the backdoor region of the Gα subunit is sufficient for triggering the activation process.

A transient interaction between the phosphate binding loop and switch I contributes to the allosteric network between receptor and nucleotide in Gαi1

Receptor-mediated activation of the Gα subunit of heterotrimeric G proteins requires allosteric communication between the receptor binding site and the guanine nucleotide binding site, which are separated by >30 Å. Structural changes in the allosteric network connecting these sites are predicted to be transient in the wild-type Gα subunit, making studies of these connections challenging. In the current work, site-directed mutants that alter the energy barriers between the activation states are used as tools to better understand the transient features of allosteric signaling in the Gα subunit. The observed differences in relative receptor affinity for intact Gαi1 subunits versus C-terminal Gαi1 peptides harboring the K345L mutation are consistent with this mutation modulating the allosteric network in the protein subunit. Measurement of nucleotide exchange rates, affinity for metarhodopsin II, and thermostability suggest that the K345L Gαi1 variant has reduced stability in both the GDP-bound and nucleotide-free states as compared with wild type but similar stability in the GTPγS-bound state. High resolution x-ray crystal structures reveal conformational changes accompanying the destabilization of the GDP-bound state. Of these, the conformation for Switch I was stabilized by an ionic interaction with the phosphate binding loop. Further site-directed mutagenesis suggests that this interaction between Switch I and the phosphate binding loop is important for receptor-mediated nucleotide exchange in the wild-type Gαi1 subunit.

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.

Energetic analysis of the rhodopsin-G-protein complex links the α5 helix to GDP release

We present a model of interaction of Gi protein with the activated receptor (R*) rhodopsin, which pinpoints energetic contributions to activation and reconciles the β2 adrenergic receptor-Gs crystal structure with new and previously published experimental data. In silico analysis demonstrated energetic changes when the Gα C-terminal helix (α5) interacts with the R* cytoplasmic pocket, thus leading to displacement of the helical domain and GDP release. The model features a less dramatic domain opening compared with the crystal structure. The α5 helix undergoes a 63° rotation, accompanied by a 5.7-Å translation, that reorganizes interfaces between α5 and α1 helices and between α5 and β6-α5. Changes in the β6-α5 loop displace αG. All of these movements lead to opening of the GDP-binding pocket. The model creates a roadmap for experimental studies of receptor-mediated G-protein activation.

Conformational flexibility and structural dynamics in GPCR-mediated G protein activation: a perspective

Structure and dynamics of G proteins and their cognate receptors, both alone and in complex, are becoming increasingly accessible to experimental techniques. Understanding the conformational changes and timelines that govern these changes can lead to new insights into the processes of ligand binding and associated G protein activation. Experimental systems may involve the use of, or otherwise stabilize, non-native environments. This can complicate our understanding of structural and dynamic features of processes such as the ionic lock, tryptophan toggle, and G protein flexibility. While elements in the receptor's transmembrane helices and the C-terminal α5 helix of Gα undergo well-defined structural changes, regions subject to conformational flexibility may be important in fine-tuning the interactions between activated receptors and G proteins. The pairing of computational and experimental approaches will continue to provide powerful tools to probe the conformation and dynamics of receptor-mediated G protein activation.

Differences in intradomain and interdomain motion confer distinct activation properties to structurally similar Gα proteins

Proteins with similar crystal structures can have dissimilar rates of substrate binding and catalysis. Here we used molecular dynamics simulations and biochemical analysis to determine the role of intradomain and interdomain motions in conferring distinct activation rates to two Gα proteins, Gα(i1) and GPA1. Despite high structural similarity, GPA1 can activate itself without a receptor, whereas Gα(i1) cannot. We found that motions in these proteins vary greatly in type and frequency. Whereas motion is greatest in the Ras domain of Gα(i1), it is greatest in helices αA and αB from the helical domain of GPA1. Using protein chimeras, we show that helix αA from GPA1 is sufficient to confer rapid activation to Gα(i1). Gα(i1) has less intradomain motion than GPA1 and instead displays interdomain displacement resembling that observed in a receptor-heterotrimer crystal complex. Thus, structurally similar proteins can have distinct atomic motions that confer distinct activation mechanisms.

A constitutively active Gα subunit provides insights into the mechanism of G protein activation

The activation of Gα subunits of heterotrimeric G proteins by G protein-coupled receptors (GPCRs) is a critical event underlying a variety of biological responses. Understanding how G proteins are activated will require structural and biochemical analyses of GPCRs complexed to their G protein partners, together with structure-function studies of Gα mutants that shed light on the different steps in the activation pathway. Previously, we reported that the substitution of a glycine for a proline at position 56 within the linker region connecting the helical and GTP-binding domains of a Gα chimera, designated αT*, yields a more readily exchangeable state for guanine nucleotides. Here we show that GDP-GTP exchange on αT*(G56P), in the presence of the light-activated GPCR, rhodopsin (R*), is less sensitive to the β1γ1 subunit complex than to wild-type αT*. We determined the X-ray crystal structure for the αT*(G56P) mutant and found that the G56P substitution leads to concerted changes that are transmitted to the conformationally sensitive switch regions, the α4-β6 loop, and the β6 strand. The α4-β6 loop has been proposed to be a GPCR contact site that signals to the TCAT motif and weakens the binding of the guanine ring of GDP, whereas the switch regions are the contact sites for the β1γ1 complex. Collectively, these biochemical and structural data lead us to suggest that αT*(G56P) may be adopting a conformation that is normally induced within Gα subunits by the combined actions of a GPCR and a Gβγ subunit complex during the G protein activation event.

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.

Structural basis for the specific inhibition of heterotrimeric Gq protein by a small molecule

Heterotrimeric GTP-binding proteins (G proteins) transmit extracellular stimuli perceived by G protein-coupled receptors (GPCRs) to intracellular signaling cascades. Hundreds of GPCRs exist in humans and are the targets of a large percentage of the pharmaceutical drugs used today. Because G proteins are regulated by GPCRs, small molecules that directly modulate G proteins have the potential to become therapeutic agents. However, strategies to develop modulators have been hampered by a lack of structural knowledge of targeting sites for specific modulator binding. Here we present the mechanism of action of the cyclic depsipeptide YM-254890, which is a recently discovered Gq-selective inhibitor. YM-254890 specifically inhibits the GDP/GTP exchange reaction of alpha subunit of Gq protein (Galphaq) by inhibiting the GDP release from Galphaq. X-ray crystal structure analysis of the Galphaqbetagamma-YM-254890 complex shows that YM-254890 binds the hydrophobic cleft between two interdomain linkers connecting the GTPase and helical domains of the Galphaq. The binding stabilizes an inactive GDP-bound form through direct interactions with switch I and impairs the linker flexibility. Our studies provide a novel targeting site for the development of small molecules that selectively inhibit each Galpha subunit and an insight into the molecular mechanism of G protein activation.

Deciphering the deformation modes associated with function retention and specialization in members of the Ras superfamily

The evolutionary and physical deformability patterns of members of the Ras GTPase superfamily were investigated by Principal Component and Elastic Network-Normal Mode analyses. The study helped to decipher the dynamics information encrypted into the conserved core and to separate the trans-family intrinsic flexibility associated with a common function from the protein motions related to functional specialization of selected families or family members. The conserved core is dynamically divided into two lobes. The deformation modes, which allow the Ras GTPases to accomplish their switching function, are conserved along evolution and are localized in lobe 1 portions close to the nucleotide. These modes lead to functional specialization when associated with evolution-driven deformations of protein portions essentially located in lobe 2, distal from the nucleotide, and involved in peculiar interactions with membrane, guanine nucleotide exchange factors, or effectors. Overall, a complete picture of the functional and evolutionary dynamics of the Ras superfamily emerges.

G protein inactive and active forms investigated by simulation methods

Molecular dynamics and computational alanine scanning techniques have been used to investigate G proteins in their inactive state (the Galpha(i1)beta(1)gamma(2) heterotrimer) as well as in their empty and monomeric active states (Galpha(i1) subunit). We find that: (i) the residue Q204 of Galpha(i1) plays a key role for binding Gbeta(1)gamma(2) and is classified among the most relevant in the interaction with a key cellular partner, the so-called regulator of G protein signaling protein. The mutation of this residue to L, which is observed in a variety of diseases, provides still fair stability to the inactive state because of the formation of van der Waals interactions. (ii) The empty state turns out to adopt some structural features of the active one, including a previously unrecognized rearrangement of a key residue (K46). (iii) The so-called Switch IV region increases its mobility on passing from the empty to the active state, and, even more, to the inactive state. Such change in mobility could be important for its several structural and functional roles. (iv) A large scale motion of the helical domain in the inactive state might be important for GDP release upon activation by GPCR, consistently with experimental data.

Structural evidence for a sequential release mechanism for activation of heterotrimeric G proteins

Heptahelical G-protein (heterotrimeric guanine nucleotide-binding protein)-coupled receptors couple to heterotrimeric G proteins to relay extracellular signals to intracellular signaling networks, but the molecular mechanism underlying guanosine 5'-diphosphate (GDP) release by the G protein alpha-subunit is not well understood. Amino acid substitutions in the conserved alpha5 helix of G(i), which extends from the C-terminal region to the nucleotide-binding pocket, cause dramatic increases in basal (receptor-independent) GDP release rates. For example, mutant Galpha(i1)-T329A shows an 18-fold increase in basal GDP release rate and, when expressed in culture, it causes a significant decrease in forskolin-stimulated cAMP accumulation. The crystal structure of Galpha(i1)-T329A.GDP shows substantial conformational rearrangement of the switch I region and additional striking alterations of side chains lining the catalytic pocket that disrupt the Mg(+2) coordination sphere and dislodge bound Mg(+2). We propose a "sequential release" mechanism whereby a transient conformational change in the alpha5 helix alters switch I to induce GDP release. Interestingly, this mechanistic model for heterotrimeric G protein activation is similar to that suggested for the activation of the plant small G protein Rop4 by RopGEF8.

Combining an Elastic Network With a Coarse-Grained Molecular Force Field: Structure, Dynamics, and Intermolecular Recognition

Structure-based and physics-based coarse-grained molecular force fields have become attractive approaches to gain mechanistic insight into the function of large biomolecular assemblies. Here, we study how both approaches can be combined into a single representation, that we term ELNEDIN. In this representation an elastic network is used as a structural scaffold to describe and maintain the overall shape of a protein and a physics-based coarse-grained model (MARTINI-2.1) is used to describe both inter- and intramolecular interactions in the system. The results show that when used in molecular dynamics simulations ELNEDIN models can be built so that the resulting structural and dynamical properties of a protein, including its collective motions, are comparable to those obtained using atomistic protein models. We then evaluate the behavior of such models in (1) long, microsecond time-scale, simulations, (2) the modeling of very large macromolecular assemblies, a viral capsid, and (3) the study of a protein-protein association process, the reassembly of the ROP homodimer. The results for this series of tests indicate that ELNEDIN models allow microsecond time-scale molecular dynamics simulations to be carried out readily, that large biological entities such as the viral capsid of the cowpea mosaic virus can be stably modeled as assemblies of independent ELNEDIN models, and that ELNEDIN models show significant promise for modeling protein-protein association processes.

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.

Computational molecular biology approaches to ligand-target interactions

Binding of small molecules to their targets triggers complex pathways. Computational approaches are keys for predictions of the molecular events involved in such cascades. Here we review current efforts at characterizing the molecular determinants in the largest membrane-bound receptor family, the G-protein-coupled receptors (GPCRs). We focus on odorant receptors, which constitute more than half GPCRs. The work presented in this review uncovers structural and energetic aspects of components of the cellular cascade. Finally, a computational approach in the context of radioactive boron-based antitumoral therapies is briefly described.

A novel 9-bp insertion detected in steroid 21-hydroxylase gene (CYP21A2): prediction of its structural and functional implications by computational methods

BACKGROUND

Steroid 21-hydroxylase deficiency is the most common cause of congenital adrenal hyperplasia (CAH). Detection of underlying mutations in CYP21A2 gene encoding steroid 21-hydroxylase enzyme is helpful both for confirmation of diagnosis and management of CAH patients. Here we report a novel 9-bp insertion in CYP21A2 gene and its structural and functional consequences on P450c21 protein by molecular modeling and molecular dynamics simulations methods.

METHODS

A 30-day-old child was referred to our laboratory for molecular diagnosis of CAH. Sequencing of the entire CYP21A2 gene revealed a novel insertion (duplication) of 9-bp in exon 2 of one allele and a well-known mutation I172N in exon 4 of other allele. Molecular modeling and simulation studies were carried out to understand the plausible structural and functional implications caused by the novel mutation.

RESULTS

Insertion of the nine bases in exon 2 resulted in addition of three valine residues at codon 71 of the P450c21 protein. Molecular dynamics simulations revealed that the mutant exhibits a faster unfolding kinetics and an overall destabilization of the structure due to the triple valine insertion was also observed.

CONCLUSION

The novel 9-bp insertion in exon 2 of CYP21A2 genesignificantly lowers the structural stability of P450c21 thereby leading to the probable loss of its function.

How do receptors activate G proteins?

Heterotrimeric G proteins couple the activation of heptahelical receptors at the cell surface to the intracellular signaling cascades that mediate the physiological responses to extracellular stimuli. G proteins are molecular switches that are activated by receptor-catalyzed GTP for GDP exchange on the G protein alpha subunit, which is the rate-limiting step in the activation of all downstream signaling. Despite the important biological role of the receptor-G protein interaction, relatively little is known about the structure of the complex and how it leads to nucleotide exchange. This chapter will describe what is known about receptor and G protein structure and outline a strategy for assembling the current data into improved models for the receptor-G protein complex that will hopefully answer the question as to how receptors flip the G protein switch.

Exploring allosteric coupling in the alpha-subunit of heterotrimeric G proteins using evolutionary and ensemble-based approaches

BACKGROUND

Allosteric coupling, which can be defined as propagation of a perturbation at one region of the protein molecule (such as ligand binding) to distant sites in the same molecule, constitutes the most general mechanism of regulation of protein function. However, unlike molecular details of ligand binding, structural elements involved in allosteric effects are difficult to diagnose. Here, we identified allosteric linkages in the alpha-subunits of heterotrimeric G proteins, which were evolved to transmit membrane receptor signals by allosteric mechanisms, by using two different approaches that utilize fundamentally different and independent information.

RESULTS

We analyzed: 1) correlated mutations in the family of G protein alpha-subunits, and 2) cooperativity of the native state ensemble of the Galphai1 or transducin. The combination of these approaches not only recovered already-known details such as the switch regions that change conformation upon nucleotide exchange, and those regions that are involved in receptor, effector or Gbetagamma interactions (indicating that the predictions of the analyses can be viewed with a measure of confidence), but also predicted new sites that are potentially involved in allosteric communication in the Galpha protein. A summary of the new sites found in the present analysis, which were not apparent in crystallographic data, is given along with known functional and structural information. Implications of the results are discussed.

CONCLUSION

A set of residues and/or structural elements that are potentially involved in allosteric communication in Galpha is presented. This information can be used as a guide to structural, spectroscopic, mutational, and theoretical studies on the allosteric network in Galpha proteins, which will provide a better understanding of G protein-mediated signal transduction.

Mechanisms of inter- and intramolecular communication in GPCRs and G proteins

This study represents the first attempt to couple, by computational experiments, the mechanisms of intramolecular and intermolecular communication concerning a guanidine nucleotide exchange factor (GEF), the thromboxane A2 receptor (TXA2R), and the cognate G protein (Gq) in its heterotrimeric GDP-bound state. Two-way pathways mediate the communication between the receptor-G protein interface and both the agonist binding site of the receptor and the nucleotide binding site of the G protein. The increase in solvent accessibility in the neighborhoods of the highly conserved E/DRY receptor motif, in response to agonist binding, is instrumental in favoring the penetration of the C-terminus of Gqalpha in between the cytosolic ends of H3, H5, and H6. The arginine of the E/DRY motif is predicted to be an important mediator of the intramolecular and intermolecular communication involving the TXA2R. The receptor-G protein interface is predicted to involve multiple regions from the receptor and the G protein alpha-subunit. However, receptor contacts with the C-terminus of the alpha5-helix seem to be the major players in the receptor-catalyzed motion of the alpha-helical domain with respect to the Ras-like domain and in the formation of a nucleotide exit route in between the alphaF-helix and beta6/alpha5 loop of Gqalpha. The inferences from this study are of wide interest, as they are expected to apply to the whole rhodopsin family, given also the considerable G protein promiscuity.

Heterotrimeric G protein activation by G-protein-coupled receptors

Heterotrimeric G proteins have a crucial role as molecular switches in signal transduction pathways mediated by G-protein-coupled receptors. Extracellular stimuli activate these receptors, which then catalyse GTP-GDP exchange on the G protein alpha-subunit. The complex series of interactions and conformational changes that connect agonist binding to G protein activation raise various interesting questions about the structure, biomechanics, kinetics and specificity of signal transduction across the plasma membrane.

Probing the activation-promoted structural rearrangements in preassembled receptor-G protein complexes

Activation of heterotrimeric G proteins by their cognate seven transmembrane domain receptors is believed to involve conformational changes propagated from the receptor to the G proteins. However, the nature of these changes remains unknown. We monitored the conformational rearrangements at the interfaces between receptors and G proteins and between G protein subunits by measuring bioluminescence resonance energy transfer between probes inserted at multiple sites in receptor-G protein complexes. Using the data obtained for the alpha(2A)AR-G alpha(i1) beta1gamma2 complex and the available crystal structures of G alpha(i1) beta1gamma2, we propose a model wherein agonist binding induces conformational reorganization of a preexisting receptor-G protein complex, leading the G alpha-G betagamma interface to open but not dissociate. This conformational change may represent the movement required to allow nucleotide exit from the G alpha subunit, thus reflecting the initial activation event.

Mechanism of the receptor-catalyzed activation of heterotrimeric G proteins

Heptahelical receptors activate intracellular signaling pathways by catalyzing GTP for GDP exchange on the heterotrimeric G protein alpha subunit (G alpha). Despite the crucial role of this process in cell signaling, little is known about the mechanism of G protein activation. Here we explore the structural basis for receptor-mediated GDP release using electron paramagnetic resonance spectroscopy. Binding to the activated receptor (R*) causes an apparent rigid-body movement of the alpha5 helix of G alpha that would perturb GDP binding at the beta6-alpha5 loop. This movement was not observed when a flexible loop was inserted between the alpha5 helix and the R*-binding C terminus, which uncouples R* binding from nucleotide exchange, suggesting that this movement is necessary for GDP release. These data provide the first direct observation of R*-mediated conformational changes in G proteins and define the structural basis for GDP release from G alpha.

Molecular modeling of the effects of mutant alleles on chalcone synthase protein structure

Chalcone synthase (CHS) catalyzes the first committed step in flavonoid biosynthesis, a major pathway of plant secondary metabolism. An allelic series for the Arabidopsis CHS locus, tt4, was previously characterized at the gene, protein, and end-product levels. In an effort to deduce the molecular basis for the observed phenotypes, homology models were generated for five of the tt4 proteins based on the crystal structure of CHS2 from Medicago. Molecular dynamics simulations provided insights into how even those substitutions that are not in close spatial proximity to key functional residues may still alter the architecture and dynamic movement of the enzyme, with dramatic effects on enzyme function. Simulations carried out at different temperatures pointed to optimized positioning of key residues in the active site or dimerization domain, rather than enhancement of overall structure, as underlying the higher activity of two temperature-sensitive variants at lower temperatures. Extending this type of analysis to account for protein-protein interactions may offer additional insights into the mechanisms by which single amino-acid substitutions can affect diverse aspects of protein function.

Open interface and large quaternary structure movements in 3D domain swapped proteins: insights from molecular dynamics simulations of the C-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the beta-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.