The recent expansion of GPCR crystal structures provides the opportunity to assess the performance of structure-based drug design methods for the GPCR superfamily. Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA)-based methods are commonly used for binding affinity prediction, as they provide an intermediate compromise of speed and accuracy between the empirical scoring functions used in docking and more robust free energy perturbation methods. In this study, we systematically assessed the performance of MM/PBSA in predicting experimental binding free energies using twenty Class A GPCR crystal structures and 934 known ligands. Correlations between predicted and experimental binding free energies varied significantly between individual targets, ranging from r = - 0.334 in the inactive-state CB1 cannabinoid receptor to r = 0.781 in the active-state CB1 cannabinoid receptor, while average correlation across all twenty targets was relatively poor (r = 0.183). MM/PBSA provided better predictions of binding free energies compared to docking scores in eight out of the twenty GPCR targets while performing worse for four targets. MM/PBSA binding affinity predictions calculated using a single, energy minimized structure provided comparable predictions to sampling from molecular dynamics simulations and may be more efficient when computational cost becomes restrictive. Additionally, we observed that restricting MM/PBSA calculations to ligands with a high degree of structural similarity to the crystal structure ligands improved performance in several cases. In conclusion, while MM/PBSA remains a valuable tool for GPCR structure-based drug design, its performance in predicting the binding free energies of GPCR ligands remains highly system-specific as demonstrated in a subset of twenty Class A GPCRs, and validation of MM/PBSA-based methods for each individual case is recommended before prospective use.
Structure-based virtual screening offers a good opportunity for the discovery of selective M1 muscarinic acetylcholine receptor (mAChR) agonists for the treatment of Alzheimer's disease. However, no 3-D structure of an M1 mAChR is yet available and the homology models that have been previously reported are only able to identify antagonists in virtual screening experiments. In this study, we generated a homology model of the human M1 mAChR, based on the crystal structure of an M3 mAChR as the template. This initial model was modified, using the agonist-bound crystal structure of a β2-adrenergic receptor as a guide, to give two possible activated structures. The T192 side chain was adjusted in both structures and one of the structures also had the whole of transmembrane (TM) 5 rotated and tilted toward the inner channel of the transmembrane region. The binding sites of all three structures were then refined by induced-fit docking (IFD) with acetylcholine. Virtual screening experiments showed that all three refined models could efficiently differentiate agonists from decoy molecules, with the TM5-modified models also giving good agonist/antagonist selectivity. The whole range of agonists and antagonists was observed to bind within the orthosteric site of the structure obtained by IFD refinement alone, implying that it has inactive state character. In contrast, the two TM5-modified structures were unable to accommodate the antagonists, supporting the proposition that they possess activated state character.
Human sweet taste receptor (hSTR) recognizes a wide array of sweeteners, resulting in sweet taste perception. Maltitol and lactitol have been extensively used in place of sucrose due to their capability to prevent dental caries. Herein, several molecular modeling approaches were applied to investigate the structural and energetic properties of these two polyols/hSTR complexes. Triplicate 500 ns molecular dynamics (MD) simulations and molecular mechanics/generalized Born surface area (MM/GBSA)-based free energy calculations revealed that the TAS1R2 monomer is the preferential binding site for maltitol and lactitol rather than the TAS1R3 region. Several polar residues (D142, S144, Y215, D278, E302, R383, and especially N143) were involved in polyols binding through electrostatic attractions and H-bond formations. The molecular complexation process not only induced the stable form of ligands but also stimulated the conformational adaptation of the TAS1R2 monomer to become a close-packed structure through an induced-fit mechanism. Notably, the binding affinity of the maltitol/TAS1R2 complex (ΔGbind of -17.93 ± 1.49 kcal/mol) was significantly higher than that of the lactitol/TAS1R2 system (-8.53 ± 1.78 kcal/mol), in line with the experimental relative sweetness. These findings provide an in-depth understanding of the differences in the sweetness response between maltitol and lactitol, which could be helpful to design novel polyol derivatives with higher sweet taste perception.
Non-rod non-cone photopigments in the eyes and the brain can directly mediate non-visual functions of light in non-mammals. This was supported by our recent findings on vertebrate ancient long (VAL)-opsin photopigments encoded by the val-opsinA (valopa) and val-opsinB (valopb) genes in zebrafish. However, the physiological functions of valop isoforms remain unknown. Here, we generated valop-mutant zebrafish using CRISPR/Cas genome editing, and examined the phenotypes of loss-of-function mutants. F0 mosaic mutations and germline transmission were confirmed via targeted insertions and/or deletions in the valopa or valopb gene in F1 mutants. Based on in silico analysis, frameshift mutations converted VAL-opsin proteins to non-functional truncated forms with pre-mature stop codons. Most F1 eggs or embryos from F0 female valopa/b mutants showed either no or only partial chorion elevation, and the eggs or embryos died within 26 hour-post-fertilization. However, most F1 embryos from F0 male valopa mutant developed but hatched late compared to wild-type embryos, which hatched at 4 day-post-fertilization. Late-hatched F1 offspring included wild-type and mutants, indicating the parental effects of valop knockout. This study shows valop gene knockout affects chorion formation and embryonic hatching in the zebrafish.
Phytochemical investigation on the bark of Goniothalamus tapis Miq. and G. uvaroides King has resulted in the isolation of eight styryl-lactones, (-)-cryptomeridiol, liriodenine, 3-methyl-1H-benz[f]indole-4,9-dione, (-)-stigmasterol and dimethyl terephthalate. The structures of the compounds were elucidated by spectroscopic techniques. The compounds were evaluated for their effect on platelet-activating factor (PAF) receptor binding on rabbit platelets using (3) H-PAF as a ligand. Among the compounds tested, (-)-cryptomeridiol, (+)-goniothalamin and (+)-isoaltholactone exhibited a significant and concentration-dependent inhibitory effect on PAF receptor binding, with inhibitory concentration (IC)(50) values of 17.5, 19.7 and 46.5 µm, respectively. The inhibitory effects of the first two compounds were comparable to that obtained from the positive control, cedrol. The results indicated that these compounds were strong PAF receptor binding inhibitors.