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
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
- NSF-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
-
-
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
-
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
-
Proteins involved in immune checkpoint pathways, such as CTLA4, PD1, and PD-L1, have become important targets for cancer immunotherapy; however, development of small molecule drugs targeting these pathways has proven difficult due to the nature of their protein–protein interfaces. Here, using a hierarchy of computational techniques, we design a cyclic peptide that binds CTLA4 and follow this with experimental verification of binding and biological activity, using bio-layer interferometry, cell culture, and a mouse tumor model. Beginning from a template excised from the X-ray structure of the CTLA4:B7-2 complex, we generate several peptide sequences using flexible docking and modeling steps. These peptides are cyclized head-to-tail to improve structural and proteolytic stability and screened using molecular dynamics simulation and MM-GBSA calculation. The standard binding free energies for shortlisted peptides are then calculated in explicit-solvent simulation using a rigorous multistep technique. The most promising peptide, cyc(EIDTVLTPTGWVAKRYS), yields the standard free energy −6.6 ± 3.5 kcal mol^−1, which corresponds to a dissociation constant of ∼15 μmol L^−1. The binding affinity of this peptide for CTLA4 is measured experimentally (31 ± 4 μmol L^−1) using bio-layer interferometry. Treatment with this peptide inhibited tumor growth in a co-culture of Lewis lung carcinoma (LLC) cells and antigen primed T cells, as well as in mice with an orthotropic Lewis lung carcinoma allograft model.more » « less
-
more » « less
Abstract Accurate estimation of solvation free energy (SFE) lays the foundation for accurate prediction of binding free energy. The Poisson‐Boltzmann (PB) or generalized Born (GB) combined with surface area (SA) continuum solvation method (PBSA and GBSA) have been widely used in SFE calculations because they can achieve good balance between accuracy and efficiency. However, the accuracy of these methods can be affected by several factors such as the charge models, polar and nonpolar SFE calculation methods and the atom radii used in the calculation. In this work, the performance of the ABCG2 (AM1‐BCC‐GAFF2) charge model as well as other two charge models, that is, RESP (Restrained Electrostatic Potential) and AM1‐BCC (Austin Model 1‐bond charge corrections), on the SFE prediction of 544 small molecules in water by PBSA/GBSA was evaluated. In order to improve the performance of the PBSA prediction based on the ABCG2 charge, we further explored the influence of atom radii on the prediction accuracy and yielded a set of atom radius parameters for more accurate SFE prediction using PBSA based on the ABCG2/GAFF2 by reproducing the thermodynamic integration (TI) calculation results. The PB radius parameters of carbon, oxygen, sulfur, phosphorus, chloride, bromide and iodine, were adjusted. New atom types,
on oi hn1 hn2 hn3
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1155/2020/5314187
