SAMPL7: Host-guest binding prediction by molecular dynamics and quantum mechanics.

Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) challenges provide routes to compare chemical quantities determined using computational chemistry approaches to experimental measurements that are shared after the competition. For this effort, several computational methods have been used to calculate the binding energies of Octa Acid (OA) and exo-Octa Acid (exoOA) host-guest systems for SAMPL7. The initial poses for molecular dynamics (MD) were generated by molecular docking. Binding free energy calculations were performed using molecular mechanics combined with Poisson-Boltzmann or generalized Born surface area solvation (MMPBSA/MMGBSA) approaches. The factors that affect the utility of the MMPBSA/MMGBSA approaches including solvation, partial charge, and solute entropy models were also analyzed. In addition to MD calculations, quantum mechanics (QM) calculations were performed using several different density functional theory (DFT) approaches. From SAMPL6 results, B3PW91-D3 was found to overestimate binding energies though it was effective for geometry optimizations, so it was considered for the DFT geometry optimizations in the current study, with single-point energy calculations carried out with B2PLYP-D3 with double-, triple-, and quadruple-ζ level basis sets. Accounting for dispersion effects, and solvation models was deemed essential for the predictions. MMGBSA and MMPBSA correlated better to experiment when used in conjunction with an empirical/linear correction.

A novel series of cysteine-dependent, allosteric inverse agonists of the nuclear receptor RORγt.

Inhibition of the nuclear receptor Retinoic Acid Receptor-Related Orphan Receptor γt (RORγt) is a promising strategy for the treatment of autoimmune diseases. In this paper, we describe a series of allosteric, cysteine-dependent, inverse agonists of RORγt. Site-directed mutagenesis and molecular dynamics simulations are supportive of a mechanism of action through specific binding to Cys476 on alpha helix 11 of the ligand binding domain (LBD). Representative compounds in the series selectively inhibit RORγt, potently suppress interleukin-17A (IL-17A) production by human CD4+ T cells, and inhibit T helper 17 (Th17) differentiation from human naïve CD4+ T cells. The advanced compound 13 is orally bioavailable and active at a dose of 3 mg/kg in a murine collagen-induced model of rheumatoid arthritis. Collectively, these data are supportive of the development of compound 13 in autoimmune diseases.

Binding of Per- and Polyfluoroalkyl Substances to the Human Pregnane X Receptor.

Per- and polyfluoroalkyl substances (PFASs) are a class of environmentally persistent industrial compounds that disrupt various metabolic pathways. Among the protein receptors to which PFASs bind, the human pregnane X receptor (hPXR) is found to be a host for a variety of long- and short-chain PFASs that lead to its overactivation. Overactivation of hPXR is linked to potential endocrine disruption, oxidative stress, hepatic steatosis, and adverse drug interactions. In this study, molecular dynamics (MD) is used to study the binding between hPXR and a number of PFAS compounds, including alternatives whose activity on hPXR has not been experimentally tested. This is the first-time MD is used to study the interactions between PFASs and hPXR, showing how relative binding free energies of PFASs relate to hPXR agonism. Binding free energy calculations, hydrogen bond analysis, per-residue decomposition calculations, and alanine scanning studies are done to provide further insight. Activities on hPXR for several short-chain and alternative PFAS compounds to long-chain PFASs that have yet to be reported will also be considered. These short-chain and alternative species include perfluorobutane sulfonic acid (PFBS), Gen-X (trade name for 2,3,3,3-tetrafluoro-2-heptafluoropropoxy propanoic acid), ADONA (trade name for 4,8-dioxa-3H-perfluorononanoic acid), and 6:2 fluorotelomer carboxylic acid (6:2 FTCA). The study shows key aspects of PFAS recognition on the hPXR, the link between PFAS binding to hPXR and the hPXR activity change observed upon the PFAS exposure, and the potential effects of alternative PFASs on hPXR activity.

SAMPL6 host-guest challenge: binding free energies via a multistep approach.

In this effort in the SAMPL6 host-guest binding challenge, a combination of molecular dynamics and quantum mechanical methods were used to blindly predict the host-guest binding free energies of a series of cucurbit[8]uril (CB8), octa-acid (OA), and tetramethyl octa-acid (TEMOA) hosts bound to various guest molecules in aqueous solution. Poses for host-guest systems were generated via molecular dynamics (MD) simulations and clustering analyses. The binding free energies for the structures obtained via cluster analyses of MD trajectories were calculated using the MMPBSA method and density functional theory (DFT) with the inclusion of Grimme's dispersion correction, an implicit solvation model to model the aqueous solution, and the resolution-of-the-identity (RI) approximation (MMPBSA, RI-B3PW91-D3, and RI-B3PW91, respectively). Among these three methods tested, the results for OA and TEMOA systems showed MMPBSA and RI-B3PW91-D3 methods can be used to qualitatively rank binding energies of small molecules with an overbinding by 7 and 37 kcal/mol respectively, and RI-B3PW91 gave the poorest quality results, indicating the importance of dispersion correction for the binding free energy calculations. Due to the complexity of the CB8 systems, all of the methods tested show poor correlation with the experimental results. Other quantum mechanical approaches used for the calculation of binding free energies included DFT without the RI approximation, utilizing truncated basis sets to reduce the computational cost (memory, disk space, CPU time), and a corrected dielectric constant to account for ionic strength within the implicit solvation model.

Binding of Per- and Polyfluoro-alkyl Substances to Peroxisome Proliferator-Activated Receptor Gamma.

Peroxisome proliferator receptor gamma (PPARγ), a type II nuclear receptor, fundamental in the regulation of genes, glucose metabolism, and insulin sensitization has been shown to be impacted by per- and poly-fluoroalkyl substances (PFASs). To consider the influence of PFASs upon PPARγ, the molecular interactions of 27 PFASs have been investigated. Two binding sites have been identified on the PPARγ homodimer structure: the dimer pocket and the ligand binding pocket, the former has never been studied prior. Molecular dynamics calculations were performed to gain insights about PFASs-PPARγ binding and the role of acidic and basic residues. The electrostatic interactions for acidic and basic residues far from the binding site were probed, together with their effect on PPARγ recognition. Short-range electrostatic and van der Waals interactions with nearby residues and their influence on binding energies were investigated. As the negative effects of perfluorooctane sulfonate acid were previously shown to be alleviated by one of its natural ligands, l-carnitine, here, the utility of l-carnitine as a possible inhibitor for other PFASs has been considered. A comparison of the binding patterns of l-carnitine and PFASs provides insights toward mitigation strategies for PFASs.

Chemical synthesis of human syndecan-4 glycopeptide bearing O-, N-sulfation and multiple aspartic acids for probing impacts of the glycan chain and the core peptide on biological functions.

Proteoglycans are a family of complex glycoproteins with glycosaminoglycan chains such as heparan sulfate (HS) attached to the core protein backbone. Due to the high structural heterogeneity of HS in nature, it is challenging to decipher the respective roles of the HS chain and the core protein on proteoglycan functions. While the sulfation patterns of HS dictate many activities, the core protein can potentially impact HS functions. In order to decipher this, homogeneous proteoglycan glycopeptides are needed. Herein, we report the first successful synthesis of proteoglycan glycopeptides bearing multiple aspartic acids in the core peptide and O- and N-sulfations in the glycan chain, as exemplified by the syndecan-4 glycopeptides. To overcome the high acid sensitivities of sulfates and base sensitivities of the glycopeptide during synthesis, a new synthetic approach has been developed to produce a sulfated glycan chain on a peptide sequence prone to the formation of aspartimide side products. The availability of the structurally well-defined synthetic glycopeptide enabled the investigation of their biological functions including cytokine, growth factor binding and heparanase inhibition. Interestingly, the glycopeptide exhibited context dependent enhancement or decrease of biological activities compared to the peptide or the glycan alone. The results presented herein suggest that besides varying the sulfation patterns of HS, linking the HS chain to core proteins as in proteoglycans may be an additional approach to modulate biological functions of HS in nature.

Chemoenzymatic synthesis of glycopeptides bearing rare N-glycan sequences with or without bisecting GlcNAc.

N-Linked glycopeptides have highly diverse structures in nature. Herein, we describe the first synthesis of rare multi-antennary N-glycan bearing glycan chains on 6-OH of both α1,6- and α1,3-linked mannose arms. To expedite divergent generation of N-glycan structures, four orthogonal protective groups were installed at the branching points on the core tetrasaccharide, which could be removed individually without affecting one another. In addition, the synthetic route is flexible, allowing a bisecting glucosamine moiety to be introduced at a late stage of the synthesis, further expanding the diversity of sequences that could be achieved. The bisecting glucosamine unit significantly reduced the glycosylation yields of adjacent mannoses, which was attributed to steric hindrance imposed by the glucosamine based on molecular modelling analysis. The N-glycans were then transformed to oxazoline donors and ligated with a glycopeptide acceptor from haptoglobin promoted by the wild type Arthrobacter endo-ß-N-acetylglucosaminidase (Endo-A). Endo-A exhibited interesting substrate preferences depending on donor sizes, which was rationalized through molecular dynamics studies. This is the first time that a glycopeptide bearing a bisecting N-acetyl glucosamine (GlcNAc), the rare N-glycan branch, and two LewisX trisaccharide antennae was synthesized, enabling access to this class of complex glycopeptide structures.