Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation.

The computational prediction of relative binding free energies is a crucial goal for drug discovery, and G protein-coupled receptors (GPCRs) are arguably the most important drug target class. However, they present increased complexity to model compared to soluble globular proteins. Despite breakthroughs, experimental X-ray crystal and cryo-EM structures are challenging to attain, meaning computational models of the receptor and ligand binding mode are sometimes necessary. This leads to uncertainty in understanding ligand-protein binding induced changes such as, water positioning and displacement, side chain positioning, hydrogen bond networks, and the overall structure of the hydration shell around the ligand and protein. In other words, the very elements that define structure activity relationships (SARs) and are crucial for accurate binding free energy calculations are typically more uncertain for GPCRs. In this work we use free energy perturbation (FEP) to predict the relative binding free energies for ligands of two different GPCRs. We pinpoint the key aspects for success such as the important role of key water molecules, amino acid ionization states, and the benefit of equilibration with specific ligands. Initial calculations following typical FEP setup and execution protocols delivered no correlation with experiment, but we show how results are improved in a logical and systematic way. This approach gave, in the best cases, a coefficient of determination (R2) compared with experiment in the range of 0.6-0.9 and mean unsigned errors compared to experiment of 0.6-0.7 kcal/mol. We anticipate that our findings will be applicable to other difficult-to-model protein ligand data sets and be of wide interest for the community to continue improving FE binding energy predictions.

Glucagon-Like Peptide-1 and Its Class B G Protein-Coupled Receptors: A Long March to Therapeutic Successes.

The glucagon-like peptide (GLP)-1 receptor (GLP-1R) is a class B G protein-coupled receptor (GPCR) that mediates the action of GLP-1, a peptide hormone secreted from three major tissues in humans, enteroendocrine L cells in the distal intestine, α cells in the pancreas, and the central nervous system, which exerts important actions useful in the management of type 2 diabetes mellitus and obesity, including glucose homeostasis and regulation of gastric motility and food intake. Peptidic analogs of GLP-1 have been successfully developed with enhanced bioavailability and pharmacological activity. Physiologic and biochemical studies with truncated, chimeric, and mutated peptides and GLP-1R variants, together with ligand-bound crystal structures of the extracellular domain and the first three-dimensional structures of the 7-helical transmembrane domain of class B GPCRs, have provided the basis for a two-domain-binding mechanism of GLP-1 with its cognate receptor. Although efforts in discovering therapeutically viable nonpeptidic GLP-1R agonists have been hampered, small-molecule modulators offer complementary chemical tools to peptide analogs to investigate ligand-directed biased cellular signaling of GLP-1R. The integrated pharmacological and structural information of different GLP-1 analogs and homologous receptors give new insights into the molecular determinants of GLP-1R ligand selectivity and functional activity, thereby providing novel opportunities in the design and development of more efficacious agents to treat metabolic disorders.

Improved inhibition of the histone acetyltransferase PCAF by an anacardic acid derivative.

Several lines of evidence indicate that histone acetyltransferases (HATs) are novel drug targets for treatment of diseases like, for example, cancer and inflammation. The natural product anacardic acid is a starting point for development of small molecule inhibitors of the histone acetyltransferase (HAT) p300/CBP associated factor (PCAF). In order to optimize the inhibitory potency, a binding model for PCAF inhibition by anacardic acid was proposed and new anacardic acid derivatives were designed. Ten new derivatives were synthesized using a novel synthetic route. One compound showed a twofold improved inhibitory potency for the PCAF HAT activity and a twofold improved inhibition of histone acetylation in HEP G2 cells.

Conformationally constrained opioid ligands: the Dmt-Aba and Dmt-Aia versus Dmt-Tic scaffold.

Replacement of the constrained phenylalanine analogue 1,2,3,4-tetrahydroisoquinoline-3-carboxylic acid (Tic) in the opioid Dmt-Tic-Gly-NH-Bn scaffold by the 4-amino-1,2,4,5-tetrahydro-indolo[2,3-c]azepin-3-one (Aia) and 4-amino-1,2,4,5-tetrahydro-2-benzazepin-3-one (Aba) scaffolds has led to the discovery of novel potent mu-selective agonists (Structures 5 and 12) as well as potent and selective delta-opioid receptor antagonists (Structures 9 and 15). Both stereochemistry and N-terminal N,N-dimethylation proved to be crucial factors for opioid receptor selectivity and functional bioactivity in the investigated small peptidomimetic templates. In addition to the in vitro pharmacological evaluation, automated docking models of Dmt-Tic and Dmt-Aba analogues were constructed in order to rationalize the observed structure-activity data.