Pan-cancer functional analysis of somatic mutations in G protein-coupled receptors.

G Protein-coupled receptors (GPCRs) are the most frequently exploited drug target family, moreover they are often found mutated in cancer. Here we used a dataset of mutations found in patient samples derived from the Genomic Data Commons and compared it to the natural human variance as exemplified by data from the 1000 genomes project. We explored cancer-related mutation patterns in all GPCR classes combined and individually. While the location of the mutations across the protein domains did not differ significantly in the two datasets, a mutation enrichment in cancer patients was observed among class-specific conserved motifs in GPCRs such as the Class A "DRY" motif. A Two-Entropy Analysis confirmed the correlation between residue conservation and cancer-related mutation frequency. We subsequently created a ranking of high scoring GPCRs, using a multi-objective approach (Pareto Front Ranking). Our approach was confirmed by re-discovery of established cancer targets such as the LPA and mGlu receptor families, but also discovered novel GPCRs which had not been linked to cancer before such as the P2Y Receptor 10 (P2RY10). Overall, this study presents a list of GPCRs that are amenable to experimental follow up to elucidate their role in cancer.

Open PHACTS computational protocols for in silico target validation of cellular phenotypic screens: knowing the knowns.

Phenotypic screening is in a renaissance phase and is expected by many academic and industry leaders to accelerate the discovery of new drugs for new biology. Given that phenotypic screening is per definition target agnostic, the emphasis of in silico and in vitro follow-up work is on the exploration of possible molecular mechanisms and efficacy targets underlying the biological processes interrogated by the phenotypic screening experiments. Herein, we present six exemplar computational protocols for the interpretation of cellular phenotypic screens based on the integration of compound, target, pathway, and disease data established by the IMI Open PHACTS project. The protocols annotate phenotypic hit lists and allow follow-up experiments and mechanistic conclusions. The annotations included are from ChEMBL, ChEBI, GO, WikiPathways and DisGeNET. Also provided are protocols which select from the IUPHAR/BPS Guide to PHARMACOLOGY interaction file selective compounds to probe potential targets and a correlation robot which systematically aims to identify an overlap of active compounds in both the phenotypic as well as any kinase assay. The protocols are applied to a phenotypic pre-lamin A/C splicing assay selected from the ChEMBL database to illustrate the process. The computational protocols make use of the Open PHACTS API and data and are built within the Pipeline Pilot and KNIME workflow tools.

Improving the accuracy of protein pKa calculations: conformational averaging versus the average structure.

Several methods for including the conformational flexibility of proteins in the calculation of titration curves are compared. The methods use the linearized Poisson-Boltzmann equation to calculate the electrostatic free energies of solvation and are applied to bovine pancreatic trypsin inhibitor (BPTI) and hen egg-white lysozyme (HEWL). An ensemble of conformations is generated by a molecular dynamics simulation of the proteins with explicit solvent. The average titration curve of the ensemble is calculated in three different ways: an average structure is used for the pKa calculation; the electrostatic interaction free energies are averaged and used for the pKa calculation; and the titration curve for each structure is calculated and the curves are averaged. The three averaging methods give very similar results and improve the pKa values to approximately the same degree. This suggests, in contrast to implications from other work, that the observed improvement of pKa values in the present studies is due not to averaging over an ensemble of structures, but rather to the generation of a single properly averaged structure for the pKa calculation.

The role of polar interactions in the molecular recognition of CD40L with its receptor CD40.

CD40 Ligand (CD40L) is transiently expressed on the surface of T-cells and binds to CD40, which is expressed on the surface of B-cells. This binding event leads to the differentiation, proliferation, and isotype switching of the B-cells. The physiological importance of CD40L has been demonstrated by the fact that expression of defective CD40L protein causes an immunodeficiency state characterized by high IgM and low IgG serum levels, indicating faulty T-cell dependent B-cell activation. To understand the structural basis for CD40L/CD40 association, we have used a combination of molecular modeling, mutagenesis, and X-ray crystallography. The structure of the extracellular region of CD40L was determined by protein crystallography, while the CD40 receptor was built using homology modeling based upon a novel alignment of the TNF receptor superfamily, and using the X-ray structure of the TNF receptor as a template. The model shows that the interface of the complex is composed of charged residues, with CD40L presenting basic side chains (K143, R203, R207), and CD40 presenting acidic side chains (D84, E114, E117). These residues were studied experimentally through site-directed mutagenesis, and also theoretically using electrostatic calculations with the program Delphi. The mutagenesis data explored the role of the charged residues in both CD40L and CD40 by switching to Ala (K143A, R203A, R207A of CD40L, and E74A, D84A, E114A, E117A of CD40), charge reversal (K143E, R203E, R207E of CD40L, and D84R, E114R, E117R of CD40), mutation to a polar residue (K143N, R207N, R207Q of CD40L, and D84N, E117N of CD40), and for the basic side chains in CD40L, isosteric substitution to a hydrophobic side chain (R203M, R207M). All the charge-reversal mutants and the majority of the Met and Ala substitutions led to loss of binding, suggesting that charged interactions stabilize the complex. This was supported by the Delphi calculations which confirmed that the CD40/CD40L residue pairs E74-R203, D84-R207, and E117-R207 had a net stabilizing effect on the complex. However, the substitution of hydrophilic side chains at several of the positions was tolerated, which suggests that although charged interactions stabilize the complex, charge per se is not crucial at all positions. Finally, we compared the electrostatic surface of TNF/TNFR with CD40L/CD40 and have identified a set of polar interactions surrounded by a wall of hydrophobic residues that appear to be similar but inverted between the two complexes.

Titration calculations of foot-and-mouth disease virus capsids and their stabilities as a function of pH.

Foot-and-mouth disease virus (FMDV), a non-enveloped picornavirus, is sensitive to acidic conditions. At pH values below 7 the icosahedral virus capsid, formed from 60 copies of a protomer containing four polypeptides (VP1 to 4), dissociates into 12 pentamers, releasing the viral RNA. Evidence suggests that this acid lability may assist FMDV cell entry via an endosomal pathway. Calculations of titration curves and pH-stability profiles are presented for three different strains of FMDV, O1BFS, A10(61) and A22 Iraq, and compared with experimental data for complete virions and empty capsids (which lack RNA). The finite difference Poisson-Boltzmann method was used for the calculation of electrostatic free energies with the solvent treated as a dielectric continuum. The inter-pentamer interface in the virus is formed by two protomers related by 2-fold icosahedral symmetry. As a simple model for inter-pentamer interactions, a dimer and two separate protomers were compared. The association free energy was computed by integrating the difference between the titration curves of the two species. The calculations reproduced the observed decrease in capsid stability at acidic pH but not the difference in pH sensitivities of the two type A viruses. It is shown that only residues within 15 A of the interface play a significant role in determining acid lability. For the experimentally studied pH range (5 to 7.6), histidine residues were found to dominate the pH-dependence of the stability. Two histidine residues in VP3, H142 and H145, are shown to have the greatest effect by virtue of their interactions with many polar residues across the inter-pentamer interface; the interaction of H142 with an alpha-helix in the opposite pentamer contributes only a small proportion of the destabilization energy.

PDB-based protein loop prediction: parameters for selection and methods for optimization.

An approach to loop prediction that starts with a database search is presented and analyzed. To obtain meaningful statistics, 130 loops from 21 proteins were studied. The correlation between the internal conformation of the loop and the conformation of the neighboring stem residues was examined. Distances between C(alpha) and C(beta) of the immediate neighbor residues at each end select template loops as well as more complex (e.g. three residues on either side) matching criteria. To have a high probability that the best possible loop candidate in the database is included in the set, relatively large cutoffs for matching the interatomic distances of the stem residues have to be used in the template loop selection procedure; for loops of length 5, this results in an average of 1000 loops and for loops of length 9, the number is about 1500. The required number increases only slowly with loop length, in contrast to the exponential time increase involved in direct searches of the conformational space. The best loops among the large number of candidates can be determined by ranking them with the standard CHARMM non-bonded energy function (without electrostatics) applied to the backbone and C(beta) atoms. The same representation (backbone plus C(beta)) can be used to optimize the loop orientations relative to the rest of the protein by constrained energy minimization. Target loops that have many non-bonded contacts with the protein yield better results so that analysis of the non-bonded contacts of the selected template loops is useful in determining the expected accuracy of a prediction. The method for loop selection and optimization predicted eight (out of 18) loops of up to nine residues to an RMSD better than 1.07 A relative to the crystal structure; for 17 of the 18 loops, one of the three lowest energy template loops had an RMSD of less than 1.79 A. The prediction of antibody loops from a database search is more effective than that for non-antibody loops. Provided that they belong to one of the canonical classes, very similar antibody loops are certain to exist in the database. Superposition of the stem residues for antibody loops also results in a better orientation than with arbitrary target loops because the neighboring residues tend to have a more similar beta-strand structure. Two H3 loops (for which no canonical structures have been proposed) were predicted with reasonable accuracy (RMSD of 0.49 A and 1.07 A) even though no corresponding antibody loops were in the database.

Evaluation of comparative protein modeling by MODELLER.

We evaluate 3D models of human nucleoside diphosphate kinase, mouse cellular retinoic acid binding protein I, and human eosinophil neurotoxin that were calculated by MODELLER, a program for comparative protein modeling by satisfaction of spatial restraints. The models have good stereochemistry and are at least as similar to the crystallographic structures as the closest template structures. The largest errors occur in the regions that were not aligned correctly or where the template structures are not similar to the correct structure. These regions correspond predominantly to exposed loops, insertions of any length, and non-conserved side chains. When a template structure with more than 40% sequence identity to the target protein is available, the model is likely to have about 90% of the mainchain atoms modeled with an rms deviation from the X-ray structure of approximately 1 A, in large part because the templates are likely to be that similar to the X-ray structure of the target. This rms deviation is comparable to the overall differences between refined NMR and X-ray crystallography structures of the same protein.

An immunohistochemical study of endothelial cell heterogeneity in the rat: observations in "en face" Häutchen preparations.

"En face" sheets of endothelium, in which cellular spatial relationships were maintained, were prepared from proximal pulmonary and femoral arteries and aortae of the Wistar rat, in order to visualise patterns of heterogeneity in populations of endothelial cells. These preparations, termed Häutchen, were immunolabelled with antibodies to angiotensin II, endothelin and Factor-VIII-related antigen, and visualised by an avidin/biotin peroxidase complex. Clusters of cells, which accounted for approximately 30% of the total endothelial cell population, and which were positively immunostained for angiotensin II, were found perpendicular to the longitudinal axis of the aorta and femoral artery. Cells in the pulmonary artery were immunonegative for angiotensin II. The majority of cells in all three vessels were immunopositive for endothelin; groups of intensely stained cells were present in both the femoral artery and aorta, but not in the pulmonary artery. Immunoreactivity to Factor-VIII-related antigen was heterogeneous, with intensely stained amorphous patches of endothelial cells present in the femoral artery and aorta. Häutchen preparations present an opportunity for the investigation of endothelial cell heterogeneity, both within and between vessels; this may provide a basis for the interpretation of the heterogeneity of endothelium-dependent responses in vessels of differing origin.

A model for the antagonist binding site on the adenosine A1 receptor, based on steric, electrostatic, and hydrophobic properties.

With the aid of molecular modeling, both adenosine and adenosine A1 receptor antagonists belonging to various chemical classes were compared with respect to their minimum-energy conformations and molecular electrostatic potentials, as computed by the semiempirical molecular orbital program MOPAC. Distinct steric and electrostatic similarities between adenosine and the prototypic adenosine antagonist theophylline are evident when both compounds are superimposed, with theophylline in a so-called flipped orientation. Similar patterns were found for all other A1 antagonists investigated in this study. A model for the antagonist binding site on the adenosine A1 receptor, based on steric, electrostatic, and hydrophobic properties contributing to potency, is proposed.

A study of model energetics and conformational properties of polynucleotide triplexes.

The formation of triple-stranded nucleic acid helices is studied by molecular mechanics and molecular dynamics calculations. Using standard TAT and CGG homopolymers, single, triple, and quintuple molecular replacements are made. Some of these replacements are expected to form Hoogsteen bonds and some are not. While the electrostatic and total energetic differences for base triplet mismatches were dependent on the electrostatic model chosen, clear trends in the local geometric distortions were apparent. Relationships between these model-built strand geometries and chemical probe experiments are discussed.

Molecular modeling of a putative antagonist binding site on helix III of the beta-adrenoceptor.

In recent biochemical studies it was demonstrated that residue Asp113 of the beta-adrenoceptor (beta-AR) is an indispensable amino acid for the binding of beta-AR antagonists. Earlier fluorescence studies showed that a tryptophan-rich region of the beta-AR is involved in the binding of propranolol, the prototype beta-AR antagonist. Bearing these two biochemical findings in mind, we explored the beta-AR part containing Asp113, for an energetically favorable antagonist binding site. This was done by performing molecular docking studies with the antagonist propranolol and a specific beta-AR peptide which included, besides Asp113, two possibly relevant tryptophan residues. In the docking calculations, the propranolol molecule was allowed to vary all its internal torsional angles. The receptor peptide was kept in an alpha-helix conformation, while side chains relevant to ligand binding were flexible to enable optimal adaptations to the ligand's binding conformation. By means of force-field calculations the total energy was minimized, consisting of the intramolecular energies of both ligand and receptor peptide, and the intermolecular energy. We found an antagonist binding site, consisting of amino acids Asp113 and Trp109, which enabled energetically favorable interactions with the receptor-binding groups of propranolol. According to these results, binding involves three main interaction points: (i) a reinforced ionic bond; (ii) a hydrogen bond; and (iii) a hydrophobic/charge transfer interaction. The deduced binding site shows a difference in affinity between the levo- and dextrorotatory isomers of propranolol caused by a difference in ability to form a hydrogen bond, which is in conformity with the experimentally observed stereoselectivity. Moreover, it also provides an explanation for the beta 1-selectivity of p-phenyl substituted phenoxypropanolamines like betaxolol. The p-phenyl substituent of betaxolol was shown to be sterically hindered upon binding to the beta 2-AR peptide, whereas this hindrance is very likely to be much less with the beta 1-AR peptide. Finally, the proposed antagonist binding site is discussed in the light of some recent biochemical findings and theories.

A molecular graphics study exploring a putative ligand binding site of the beta-adrenoceptor.

The recent elucidation of the primary structure of the cell membrane-bound beta-adrenoceptor has prompted us to explore putative ligand binding sites on this physiologically important receptor. By minimizing the energies of the 'prototype' ligand propranolol, (part of) the receptor and the proposed ligand-receptor complex with the aid of force field and quantum chemical calculations, we identified amino acid residue Trp313 as a highly probable candidate for interaction with the aromatic moiety of propranolol. The charge distribution on the indole nucleus of another beta-blocker, pindolol, with higher affinity for the beta-adrenoceptor, enables an even stronger interaction with the tryptophan residue. The carboxylic amino acid residue Glu306, located near the extracellular space of the cell membrane, interacts favorably with the positively charged nitrogen atom in the aliphatic side chain of the ligands. Finally, this putative model is discussed in the light of recent findings in mutagenesis studies, and compared to other ideas with respect to ligand-receptor interactions.