Evaluating the performance of MM/PBSA for binding affinity prediction using class A GPCR crystal structures.

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.

Development and Validation of Decision Forest Model for Estrogen Receptor Binding Prediction of Chemicals Using Large Data Sets.

Some chemicals in the environment possess the potential to interact with the endocrine system in the human body. Multiple receptors are involved in the endocrine system; estrogen receptor α (ERα) plays very important roles in endocrine activity and is the most studied receptor. Understanding and predicting estrogenic activity of chemicals facilitates the evaluation of their endocrine activity. Hence, we have developed a decision forest classification model to predict chemical binding to ERα using a large training data set of 3308 chemicals obtained from the U.S. Food and Drug Administration's Estrogenic Activity Database. We tested the model using cross validations and external data sets of 1641 chemicals obtained from the U.S. Environmental Protection Agency's ToxCast project. The model showed good performance in both internal (92% accuracy) and external validations (∼ 70-89% relative balanced accuracies), where the latter involved the validations of the model across different ER pathway-related assays in ToxCast. The important features that contribute to the prediction ability of the model were identified through informative descriptor analysis and were related to current knowledge of ER binding. Prediction confidence analysis revealed that the model had both high prediction confidence and accuracy for most predicted chemicals. The results demonstrated that the model constructed based on the large training data set is more accurate and robust for predicting ER binding of chemicals than the published models that have been developed using much smaller data sets. The model could be useful for the evaluation of ERα-mediated endocrine activity potential of environmental chemicals.

Molecular dynamics simulations of the adenosine A2a receptor in POPC and POPE lipid bilayers: effects of membrane on protein behavior.

Analysis of 300 ns (ns) molecular dynamics (MD) simulations of an adenosine A2a receptor (A2a AR) model, conducted in triplicate, in 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) bilayers reveals significantly different protein dynamical behavior. Principal component analysis (PCA) shows that the dissimilarities stem from interhelical rather than intrahelical motions. The difference in the hydrophobic thicknesses of these simulated lipid bilayers is potentially a significant reason for the observed difference in results. The distinct lipid headgroups might also lead to different molecular interactions and hence different protein loop motions. Overall, the A2a AR shows higher mobility and flexibility in POPC as compared to POPE.

Flavonoids with M1 muscarinic acetylcholine receptor binding activity.

Muscarinic acetylcholine receptor-active compounds have potential for the treatment of Alzheimer's disease. In this study, a series of natural and synthetic flavones and flavonols was assayed in vitro for their ability to inhibit radioligand binding at human cloned M1 muscarinic receptors. Several compounds were found to possess competitive binding affinity (Ki=40-110 µM), comparable to that of acetylcholine (Ki=59 µM). Despite the fact that these compounds lack a positively-charged ammonium group under physiological conditions, molecular modelling studies suggested that they bind to the orthosteric site of the receptor, mainly through non-polar interactions.

Conversion of a non-selective adenosine receptor antagonist into A3-selective high affinity fluorescent probes using peptide-based linkers.

Advances in fluorescence-based imaging technologies have helped propel the study of real-time biological readouts and analysis across many different areas. In particular the use of fluorescent ligands as chemical tools to study proteins such as G protein-coupled receptors (GPCRs) has received ongoing interest. Methods to improve the efficient chemical synthesis of fluorescent ligands remain of paramount importance to ensure this area of bioanalysis continues to advance. Here we report conversion of the non-selective GPCR adenosine receptor antagonist Xanthine Amine Congener into higher affinity and more receptor subtype-selective fluorescent antagonists. This was achieved through insertion and optimisation of a dipeptide linker between the adenosine receptor pharmacophore and the fluorophore. Fluorescent probe 27 containing BODIPY 630/650 (pK(D) = 9.12 ± 0.05 [hA3AR]), and BODIPY FL-containing 28 (pK(D) = 7.96 ± 0.09 [hA3AR]) demonstrated clear, displaceable membrane binding using fluorescent confocal microscopy. From in silico analysis of the docked ligand-receptor complexes of 27, we suggest regions of molecular interaction that could account for the observed selectivity of these peptide-linker based fluorescent conjugates. This general approach of converting a non-selective ligand to a selective biological tool could be applied to other ligands of interest.

Molecular dynamics simulations of the adenosine A2a receptor: structural stability, sampling, and convergence.

Molecular dynamics (MD) simulations of membrane-embedded G-protein coupled receptors (GPCRs) have rapidly gained popularity among the molecular simulation community in recent years, a trend which has an obvious link to the tremendous pharmaceutical importance of this group of receptors and the increasing availability of crystal structures. In view of the widespread use of this technique, it is of fundamental importance to ensure the reliability and robustness of the methodologies so they yield valid results and enable sufficiently accurate predictions to be made. In this work, 200 ns simulations of the A2a adenosine receptor (A2a AR) have been produced and evaluated in the light of these requirements. The conformational dynamics of the target protein, as obtained from replicate simulations in both the presence and absence of an inverse agonist ligand (ZM241385), have been investigated and compared using principal component analysis (PCA). Results show that, on this time scale, convergence of the replicates is not readily evident and dependent on the types of the protein motions considered. Thus rates of inter- as opposed to intrahelical relaxation and sampling can be different. When studied individually, we find that helices III and IV have noticeably greater stability than helices I, II, V, VI, and VII in the apo form. The addition of the inverse agonist ligand greatly improves the stability of all helices.

In vitro functional evaluation of isolaureline, dicentrine and glaucine enantiomers at 5-HT2 and α1 receptors.

Compounds with activity at serotonin (5-hydroxytryptamine) 5-HT2 and α1 adrenergic receptors have potential for the treatment of central nervous system disorders, drug addiction or overdose. Isolaureline, dicentrine and glaucine enantiomers were synthesized, and their in vitro functional activities at human 5-HT2 and adrenergic α1 receptor subtypes were evaluated. The enantiomers of isolaureline and dicentrine acted as antagonists at 5-HT2 and α1 receptors with (R)-isolaureline showing the greatest potency (pKb  = 8.14 at the 5-HT2C receptor). Both (R)- and (S)-glaucine also antagonized α1 receptors, but they behaved very differently to the other compounds at 5-HT2 receptors: (S)-glaucine acted as a partial agonist at all three 5-HT2 receptor subtypes, whereas (R)-glaucine appeared to act as a positive allosteric modulator at the 5-HT2A receptor.

Assessing GPCR homology models constructed from templates of various transmembrane sequence identities: Binding mode prediction and docking enrichment.

GPCR crystal structures have become more readily accessible in recent years. However, homology models of GPCRs continue to play an important role as many GPCR structures remain unsolved. The new crystal structures now available provide not only additional templates for homology modelling but also the opportunity to assess the performance of homology models against their respective crystal structures and gain insight into the performance of such models. In this study we have constructed homology models from templates of various transmembrane sequence identities for eight GPCR targets to better understand the relationship between transmembrane sequence identity and model quality. Model quality was assessed relative to the crystal structure in terms of structural accuracy as well as performance in two typical structure-based drug design applications: ligand binding pose prediction and docking enrichment in virtual screening. Crystal structures significantly outperformed homology models in both assessments. Accurate ligand binding pose prediction was possible but difficult to achieve using homology models, even with the use of induced fit docking. In virtual screening using homology models still conferred significant enrichment compared to random selection, with a clear benefit also observed in using models optimized through induced fit docking. Our results indicate that while homology models that are reasonably accurate structurally can be constructed, without significant refinement homology models will be outperformed by crystal structures in ligand binding pose prediction and docking enrichment regardless of the template used, primarily due to the extremely high level of structural accuracy needed for such applications.

Correction: Synthesis and evaluation of nuciferine and roemerine enantiomers as 5-HT2 and α1 receptor antagonists.

[This corrects the article DOI: 10.1039/C7MD00629B.].

Synthesis and evaluation of nuciferine and roemerine enantiomers as 5-HT2 and α1 receptor antagonists.

In this study, the (S)-enantiomers of the aporphine alkaloids, nuciferine and roemerine, were prepared via a synthetic route involving catalytic asymmetric hydrogenation and both stereoisomers were evaluated in vitro for functional activity at human 5-HT2 and adrenergic α1 receptor subtypes using a transforming growth factor-α shedding assay. Both enantiomers of each of the compounds were found to act as antagonists at 5-HT2 and α1 receptors. (R)-roemerine was the most potent compound at 5-HT2A and 5-HT2C receptors (pKb = 7.8-7.9) with good selectivity compared to (S)-roemerine at these two receptors and compared to its activity at 5-HT2B, α1A, α1B and α1D receptors.

Inhibition of cobalamin-dependent methionine synthase by substituted benzo-fused heterocycles.

The cobalamin-dependent cytosolic enzyme, methionine synthase (EC.2.1.1.13), catalyzes the remethylation of homocysteine to methionine using 5-methyltetrahydrofolate as the methyl donor. The products of this remethylation--methionine and tetrahydrofolate--participate in the active methionine and folate pathways. Impaired methionine synthase activity has been implicated in the pathogenesis of anaemias, cancer and neurological disorders. Although the need for potent and specific inhibitors of methionine synthase has been recognized, there is a lack of such agents. In this study, we designed, synthesized and evaluated the inhibitory activity of a series of substituted benzimidazoles and small benzothiadiazoles. Kinetic analysis revealed that the benzimidazoles act as competitive inhibitors of the rat liver methionine synthase, whilst the most active benzothiadiazole (IC(50) = 80 microm) exhibited characteristics of uncompetitive inhibition. A model of the methyltetrahydrofolate-binding site of the rat liver methionine synthase was constructed; docking experiments were designed to elucidate, in greater detail, the binding mode and reveal structural requirements for the design of inhibitors of methionine synthase. Our results indicate that the potency of the tested compounds is related to a planar region of the inhibitor that can be positioned in the centre of the active site, the presence of a nitro functional group and two or three probable hydrogen-bonding interactions.

Modelling the restoration of wild-type dynamic behaviour in DeltaF508-CFTR NBD1 by 8-cyclopentyl-1,3-dipropylxanthine.

Cystic fibrosis (CF) is the most frequently occurring severe, genetic disease in western populations with an incidence as high as 1 in 2500. The principal biochemical defect in CF is a mutation in a membrane transport protein, namely the cystic fibrosis transmembrane conductance regulator (CFTR), which is responsible for the conductance of chloride ions across cell membranes. In 70% of cases a single mutation in CFTR, namely the deletion of amino acid 508 (called DeltaF508) is sufficient to cause severe disease. This mutation manifests as a failure of the protein to be effectively targeted to the membrane. Recently, it has been shown that small molecule drug therapy can restore the membrane-targeting of DeltaF508-CFTR, where the mutant channel functions adequately. We have created models of the first nucleotide-binding domain (NBD1) region (which houses the proposed binding site of these restorative drugs) of the wild-type and mutant forms of human CFTR. We have simulated the dynamical behaviour of these proteins in the presence of drugs that restore trafficking of the protein. Our results indicate that there are particular modes of dynamic motion that are distinguishable between wild-type and mutant CFTR. These regions of motion are localized in the regions of the DeltaF508 mutation and the drug-binding regions. The simulations of drug binding indicate that wild-type dynamic motions are restored in these regions. We conclude therefore that these drugs are able to alter the dynamic properties of DeltaF508-CFTR such that the drug-bound mutant protein more closely resembles the wild-type protein dynamic behaviour, and hence we hypothesize that it is this that allows for correct targeting to the membrane.

Crystallization and preliminary X-ray characterization of the Bacillus amyloliquefaciens YwrO enzyme.

CB1954 is an anticancer prodrug that is currently in clinical trials coupled with the Escherichia coli flavoenzyme nitroreductase (NTR) for use in directed-enzyme prodrug therapy (DEPT). The NTR enzyme is responsible for the conversion of the prodrug into a cytotoxic agent. The bifunctional alkylating agent produced by this bioactivation process leads to DNA damage and death of cancer cells. Recently, a novel flavoenzyme from Bacillus amyloliquefaciens, YwrO (Bam YwrO), was reported to be able to reduce CB1954 from its noncytotoxic form into its active form. The crystallization and preliminary X-ray diffraction analysis of two crystal forms of Bam YwrO are reported. The first crystal form is orthorhombic, with space group P22(1)2(1), and diffracts X-rays to 2.18 A resolution. The second crystal form is tetragonal, with space group P4(1), and diffracts X-rays to 3.4 A. Determination of the Bam YwrO crystal structure will provide an understanding of the molecular recognition between this enzyme and the anticancer prodrug CB1954.

A dual-application poly (dl-lactic-co-glycolic) acid (PLGA)-chitosan composite scaffold for potential use in bone tissue engineering.

The development of patient-friendly alternatives to bone-graft procedures is the driving force for new frontiers in bone tissue engineering. Poly (dl-lactic-co-glycolic acid) (PLGA) and chitosan are well-studied and easy-to-process polymers from which scaffolds can be fabricated. In this study, a novel dual-application scaffold system was formulated from porous PLGA and protein-loaded PLGA/chitosan microspheres. Physicochemical and in vitro protein release attributes were established. The therapeutic relevance, cytocompatibility with primary human mesenchymal stem cells (hMSCs) and osteogenic properties were tested. There was a significant reduction in burst release from the composite PLGA/chitosan microspheres compared with PLGA alone. Scaffolds sintered from porous microspheres at 37 °C were significantly stronger than the PLGA control, with compressive strengths of 0.846 ± 0.272 MPa and 0.406 ± 0.265 MPa, respectively (p < 0.05). The formulation also sintered at 37 °C following injection through a needle, demonstrating its injectable potential. The scaffolds demonstrated cytocompatibility, with increased cell numbers observed over an 8-day study period. Von Kossa and immunostaining of the hMSC-scaffolds confirmed their osteogenic potential with the ability to sinter at 37 °C in situ.

Binding of the anticancer prodrug CB1954 to the activating enzyme NQO2 revealed by the crystal structure of their complex.

CB1954 is an attractive prodrug for directed-enzyme prodrug therapy (DEPT) and a conventional prodrug against tumors in which the enzyme NQO2 is highly expressed. We have determined the crystal structure of the NQO2-CB1954 complex to 2.0 A resolution. The binding of the prodrug is governed by hydrophobic forces, while two key electrostatic contacts determine the specific orientation of the ligand. The structure also reveals an unfavorable interaction, therefore suggesting possible avenues for DEPT-tailored engineering studies.

Physicomechanical properties of sintered scaffolds formed from porous and protein-loaded poly(DL-lactic-co-glycolic acid) microspheres for potential use in bone tissue engineering.

An injectable poly(DL-lactic-co-glycolic acid) (PLGA) system comprising both porous and protein-loaded microspheres capable of forming porous scaffolds at body temperature was developed for tissue regeneration purposes. Porous and non-porous (lysozyme loaded) PLGA microspheres were formulated to represent 'low molecular weight' 22-34 kDa, 'intermediate molecular weight' (IMW) 53 kDa and 'high molecular weight' 84-109 kDa PLGA microspheres. The respective average size of the microspheres was directly related to the polymer molecular weight. An initial burst release of lysozyme was observed from both microspheres and scaffolds on day 1. In the case of the lysozyme-loaded microspheres, this burst release was inversely related to the polymer molecular weight. Similarly, scaffolds loaded with 1 mg lysozyme/g of scaffold exhibited an inverse release relationship with polymer molecular weight. The burst release was highest amongst IMW scaffolds loaded with 2 and 3 mg/g. Sustained lysozyme release was observed after day 1 over 50 days (microspheres) and 30 days (scaffolds). The compressive strengths of the scaffolds were found to be inversely proportional to PLGA molecular weight at each lysozyme loading. Surface analysis indicated that some of the loaded lysozyme was distributed on the surfaces of the microspheres and thus responsible for the burst release observed. Overall the data demonstrates the potential of the scaffolds for use in tissue regeneration.

Toward activated homology models of the human M1 muscarinic acetylcholine receptor.

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.

The influence of substituted phenols on the sol:gel transition of hydroxypropyl methylcellulose (HPMC) aqueous solutions.

The influence of the physicochemical parameters of substituted aromatic molecules on the phase transition from sol to gel of hydroxypropyl methylcellulose (HPMC) has been investigated using a homologous series of substituted phenols. Using a turbimetric methodology, concentration dependent suppression of phase transition temperature of HPMC was observed for phenol and its derivatives, including methyl-, nitro- and chloro-substituted molecules. Although no strong direct relationship between single molecular physicochemical properties of the phenolic compounds (such as pKa, LogP and other molecular descriptors) and ΔCPT was found for the compounds tested, a successful prediction of behaviour could be obtained by using a combination of parameters. This suggested that the interaction mechanism between HPMC and the substituted aromatic moiety is a complex summation of the different molecular physicochemical properties. Identification of these potentially deleterious chemical moieties may be of value in a pharmaceutical context when considering preformulation of drug structures containing them. An incompatibility between drug and polymer may be indicative of deleterious effects resulting from formulation with hydrophilic matrix dosage forms containing cellulose ethers such as HPMC.

Structure-based identification of aporphines with selective 5-HT(2A) receptor-binding activity.

Selective blockade of the serotonin 5-HT(2A) receptor is a useful therapeutic approach for a number of disorders, including schizophrenia, insomnia and ischaemic heart disease. A series of aporphines were docked into a homology model of the rat 5-HT(2A) receptor using AutoDock. Selected compounds with high in silico binding affinities were screened in vitro using radioligand-binding assays against rat serotonin (5-HT(1A) and 5-HT(2A)) and dopamine (D1 and D2) receptors. (R)-Roemerine and (±)-nuciferine were found to have high affinity for the 5-HT(2A) receptor (K(i) = 62 and 139 nM, respectively), with (R)-roemerine showing 20- to 400-fold selectivity for the 5-HT(2A) receptor over the 5-HT(1A), D1 and D2 receptors. Investigation into the ligand-receptor interactions suggested that the selectivity of (R)-roemerine is due to it having stronger H-bonding and dipole-dipole interactions with several of the key residues in the 5-HT(2A) receptor-binding site.

Homology modeling of the human 5-HT1A, 5-HT 2A, D1, and D2 receptors: model refinement with molecular dynamics simulations and docking evaluation.

5-HT(1A) serotonin and D1 dopamine receptor agonists have been postulated to be able to improve negative and cognitive impairment symptoms of schizophrenia, while partial agonists and antagonists of the D2 and 5-HT(2A) receptors have been reported to be effective in reducing positive symptoms. There is therefore a need for well-defined homology models for the design of more selective antipsychotic agents, since no three-dimensional (3D) crystal structures of these receptors are currently available. In this study, homology models were built based on the high-resolution crystal structure of the ß(2)-adrenergic receptor (2RH1) and further refined via molecular dynamics simulations in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer system with the GROMOS96 53A6 united atom force field. Docking evaluations with representative agonists and antagonists using AutoDock 4.2 revealed binding modes in agreement with experimentally determined site-directed mutagenesis data and significant correlations between the computed and experimental pK (i) values. The models are also able to distinguish between antipsychotic agents with different selectivities and binding affinities for the four receptors, as well as to differentiate active compounds from decoys. Hence, these human 5-HT(1A), 5-HT(2A), D1 and D2 receptor homology models are capable of predicting the activities of novel ligands, and can be used as 3D templates for antipsychotic drug design and discovery.