An official website of the United States government
G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (β 2 AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (G s ). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to β 2 AR by G s protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from β 2 AR and the conformational interconversions of G s between closed and open conformations in the NE(+)–β 2 AR–G s ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter G s α subunit (G s α) conformational transitions. Our simulations showed that the interdomain movement and the stacking of G s α α1 and α5 helices are significant for increasing the distance between the G s α and β 2 AR, which may indicate a partial dissociation of G s α The distance increase commences when G s α is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of β 2 AR interacting with G s α, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.
more »
« less
To Probe Full and Partial Activation of Human Peroxisome Proliferator-Activated Receptors by Pan-Agonist Chiglitazar Using Molecular Dynamics Simulations
Chiglitazar is a promising new-generation insulin sensitizer with low reverse effects for the treatment of type II diabetes mellitus (T2DM) and has shown activity as a nonselective pan-agonist to the human peroxisome proliferator-activated receptors (PPARs) (i.e., full activation of PPAR γ and a partial activation of PPAR α and PPAR β / δ ). Yet, it has no high-resolution complex structure with PPARs and its detailed interactions and activation mechanism remain unclear. In this study, we docked chiglitazar into three experimentally resolved crystal structures of hPPAR subtypes, PPAR α , PPAR β / δ , and PPAR γ , followed by 3 μ s molecular dynamics simulations for each system. Our MM-GBSA binding energy calculation revealed that chiglitazar most favorably bound to hPPAR γ (-144.6 kcal/mol), followed by hPPAR α (-138.0 kcal/mol) and hPPAR β (-135.9 kcal/mol), and the order is consistent with the experimental data. Through the decomposition of the MM-GBSA binding energy by residue and the use of two-dimensional interaction diagrams, key residues involved in the binding of chiglitazar were identified and characterized for each complex system. Additionally, our detailed dynamics analyses support that the conformation and dynamics of helix 12 play a critical role in determining the activities of the different types of ligands (e.g., full agonist vs. partial agonist). Rather than being bent fully in the direction of the agonist versus antagonist conformation, a partial agonist can adopt a more linear conformation and have a lower degree of flexibility. Our finding may aid in further development of this new generation of medication.
more »
« less
- Award ID(s):
- 1904797
- PAR ID:
- 10147966
- Date Published:
- Journal Name:
- PPAR Research
- Volume:
- 2020
- ISSN:
- 1687-4757
- Page Range / eLocation ID:
- 1 to 24
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Structural characterization of the complex [B(β-pinane) 3 ] (1) reveals non-covalent H⋯H contacts that are consistent with the generation of London dispersion energies involving the β-pinane ligand frameworks. The homolytic fragmentations of 1 , and camphane and sabinane analogues ([B(camphane) 3 ] (2) and [B(sabinane) 3 ] (3)) were studied computationally. Isodesmic exchange results showed that London dispersion interactions are highly dependent on the terpene's stereochemistry, with the β-pinane framework providing the greatest dispersion free energy (Δ G = −7.9 kcal mol −1 ) with Grimme's dispersion correction (D3BJ) employed. PMe 3 was used to coordinate to [B(β-pinane) 3 ], giving the complex [Me 3 P–B(β-pinane) 3 ] ( 4 ), which displayed a dynamic coordination equilibrium in solution. The association process was found to be slightly endergonic at 302 K (Δ G = +0.29 kcal mol −1 ).more » « less
-
Particulate methane monooxygenase (pMMO), a membrane-bound enzyme having three subunits (α, β, and γ) and copper-containing centers, is found in most of the methanotrophs that selectively catalyze the oxidation of methane into methanol. Active sites in the pMMO of Methylosinus trichosporium OB3b were determined by docking the modeled structure with ethylbenzene, toluene, 1,3-dibutadiene, and trichloroethylene. The docking energy between the modeled pMMO structure and ethylbenzene, toluene, 1,3-dibutadiene, and trichloroethylene was −5.2, −5.7, −4.2, and −3.8 kcal/mol, respectively, suggesting the existence of more than one active site within the monomeric subunits due to the presence of multiple binding sites within the pMMO monomer. The evaluation of tunnels and cavities of the active sites and the docking results showed that each active site is specific to the radius of the substrate. To increase the catalysis rates of methane in the pMMO of M. trichosporium OB3b, selected amino acid residues interacting at the binding site of ethylbenzene, toluene, 1,3-dibutadiene, and trichloroethylene were mutated. Based on screening the strain energy, docking energy, and physiochemical properties, five mutants were downselected, B:Leu31Ser, B:Phe96Gly, B:Phe92Thr, B:Trp106Ala, and B:Tyr110Phe, which showed the docking energy of −6.3, −6.7, −6.3, −6.5, and −6.5 kcal/mol, respectively, as compared to the wild type (−5.2 kcal/mol) with ethylbenzene. These results suggest that these five mutants would likely increase methane oxidation rates compared to wild-type pMMO.more » « less
-
How far can we push the limits in removing stereoelectronic protection from an unstable intermediate? We address this question by exploring the interplay between the primary and secondary stereoelectronic effects in the Baeyer–Villiger (BV) rearrangement by experimental and computational studies of γ-OR-substituted γ-peroxylactones, the previously elusive non-strained Criegee intermediates (CI). These new cyclic peroxides were synthesized by the peroxidation of γ-ketoesters followed by in situ cyclization using a BF 3 ·Et 2 O/H 2 O 2 system. Although the primary effect (alignment of the migrating C–R m bond with the breaking O–O bond) is active in the 6-membered ring, weakening of the secondary effect (donation from the OR lone pair to the breaking C–R m bond) provides sufficient kinetic stabilization to allow the formation and isolation of stable γ-hydroperoxy-γ-peroxylactones with a methyl-substituent in the C6-position. Furthermore, supplementary protection is also provided by reactant stabilization originating from two new stereoelectronic factors, both identified and quantified for the first time in the present work. First, an unexpected boat preference in the γ-hydroperoxy-γ-peroxylactones weakens the primary stereoelectronic effects and introduces a ∼2 kcal mol −1 Curtin–Hammett penalty for reacquiring the more reactive chair conformation. Second, activation of the secondary stereoelectronic effect in the TS comes with a ∼2–3 kcal mol −1 penalty for giving up the exo-anomeric stabilization in the 6-membered Criegee intermediate. Together, the three new stereoelectronic factors (inverse α-effect, misalignment of reacting bonds in the boat conformation, and the exo-anomeric effect) illustrate the richness of stereoelectronic patterns in peroxide chemistry and provide experimentally significant kinetic stabilization to this new class of bisperoxides. Furthermore, mild reduction of γ-hydroperoxy-γ-peroxylactone with Ph 3 P produced an isolable γ-hydroxy-γ-peroxylactone, the first example of a structurally unencumbered CI where neither the primary nor the secondary stereoelectronic effect are impeded. Although this compound is relatively unstable, it does not undergo the BV reaction and instead follows a new mode of reactivity for the CI – a ring-opening process.more » « less
-
This work probed the thermal “switchability” from ethylene coordination/insertion to controlled radical polymerization of methyl acrylate (MA) for Brookhart-type α-diimine PdII catalysts. The investigation focused on the extremely bulky 2,6-bis(3,5-dimethylphenyl)-4-methylphenyl (Xyl4Ph) α-diimine N-substituents to probe reversible PdII–C bond activation in the MA-quenched Pd-capped PE intermediate and reversible trapping during radical MA polymerization. The substituent steric effect on the relative stability of various [PE–MA–PdII(ArN═CMeCMe═NAr)]+ chain-end structures and on the bond dissociation-free energy (BDFE) for the homolytic PdII–C bond cleavage has been assessed by DFT calculations at the full quantum mechanics (QM) and QM/molecular mechanics (QM/MM) methods. The structures comprise ester-chelated forms with the Pd atom bonded to the α, β, and γ C atoms as a result of 2,1 MA insertion into the PE–Pd bond and of subsequent chain walking, as well as related monodentate (ring-opened) forms resulting from the addition of MA or acetonitrile. The opened Cα-bonded form is electronically favored for smaller N-substituents, including 2,6-diisopropylphenyl (Dipp), particularly when MeCN is added, but the open Cγ-bonded form is preferred for the extremely bulky system with Ar = Xyl4Ph. The Pdα–C bond is the weakest one to cleave, with the BDFE decreasing as the Ar steric bulk is increased (31.8, 25.8, and 12.6 kcal mol–1 for Ph, Dipp, and Xyl4Ph, respectively). However, experimental investigations on the [PE–MA–PdII(ArN═CMeCMe═NAr)]+ (Ar = Xyl4Ph) macroinitiator do not show any evidence of radical formation under thermal activation conditions, while photolytic activation produces both TEMPO-trapped (TEMPO = 2,2,6,6-tetramethylpiperidinyloxy) and unsaturated MA-containing PE chains. The DFT investigation has highlighted a low-energy pathway for termination of the PE–MA• radicals by disproportionation, promoted by β-H elimination/dissociation and H-atom abstraction from the PdII–H intermediate by a second radical. This phenomenon appears to be the main reason for the failure of this PdII system to control the radical polymerization of MA by the OMRP (OMRP = organometallic-mediated radical polymerization) mechanism.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1155/2020/5314187
