High-Throughput Prediction of the Impact of Genetic Variability on Drug Sensitivity and Resistance Patterns for Clinically Relevant Epidermal Growth Factor Receptor Mutations from Atomistic Simulations.

Mutations in the kinase domain of the epidermal growth factor receptor (EGFR) can be drivers of cancer and also trigger drug resistance in patients receiving chemotherapy treatment based on kinase inhibitors. A priori knowledge of the impact of EGFR variants on drug sensitivity would help to optimize chemotherapy and design new drugs that are effective against resistant variants before they emerge in clinical trials. To this end, we explored a variety of in silico methods, from sequence-based to "state-of-the-art" atomistic simulations. We did not find any sequence signal that can provide clues on when a drug-related mutation appears or the impact of such mutations on drug activity. Low-level simulation methods provide limited qualitative information on regions where mutations are likely to cause alterations in drug activity, and they can predict around 70% of the impact of mutations on drug efficiency. High-level simulations based on nonequilibrium alchemical free energy calculations show predictive power. The integration of these "state-of-the-art" methods into a workflow implementing an interface for parallel distribution of the calculations allows its automatic and high-throughput use, even for researchers with moderate experience in molecular simulations.

Genome-wide pathway analysis identifies VEGF pathway association with oral ulceration in systemic lupus erythematosus.

BACKGROUND: Systemic lupus erythematosus (SLE) is a genetically complex rheumatic disease characterized by heterogeneous clinical manifestations of unknown etiology. Recent studies have suggested the existence of a genetic basis for SLE heterogeneity. The objective of the present study was to identify new genetic variation associated with the clinically relevant phenotypes in SLE. METHODS: A two-stage pathway-based approach was used to identify the genetic variation associated with the main clinical phenotypes in SLE. In the discovery stage, 482 SLE patients were genotyped using Illumina Human Quad610 microarrays. Association between 798 reference genetic pathways from the Molecular Signatures Database and 11 SLE phenotypes was tested using the set-based method implemented in PLINK software. Pathways significantly associated after multiple test correction were subsequently tested for replication in an independent cohort of 425 SLE patients. Using an in silico approach, we analyzed the functional effects of common SLE therapies on the replicated genetic pathways. The association of known SLE risk variants with the development of the clinical phenotypes was also analyzed. RESULTS: In the discovery stage, we found a significant association between the vascular endothelial growth factor (VEGF) pathway and oral ulceration (P value for false discovery rate (P FDR) < 0.05), and between the negative regulation signaling pathway of retinoic acid inducible gene-I/melanoma differentiation associated gene 5 and the production of antinuclear antibodies (P FDR < 0.05). In the replication stage, we validated the association between the VEGF pathway and oral ulceration. Therapies commonly used to treat mucocutaneous phenotypes in SLE were found to strongly influence VEGF pathway gene expression (P = 4.60e-4 to 5.38e-14). Analysis of known SLE risk loci identified a strong association between PTPN22 and the risk of hematologic disorder and with the development of antinuclear antibodies. CONCLUSIONS: The present study has identified VEGF genetic pathway association with the risk of oral ulceration in SLE. New therapies targeting the VEGF pathway could be more effective in reducing the severity of this phenotype. These findings represent a first step towards the understanding of the genetic basis of phenotype heterogeneity in SLE.

A deletion at ADAMTS9-MAGI1 locus is associated with psoriatic arthritis risk.

OBJECTIVE: Copy number variants (CNVs) have been associated with the risk to develop multiple autoimmune diseases. Our objective was to identify CNVs associated with the risk to develop psoriatic arthritis (PsA) using a genome-wide analysis approach. METHODS: A total of 835 patients with PsA and 1498 healthy controls were genotyped for CNVs using the Illumina HumanHap610 BeadChip genotyping platform. Genomic CNVs were characterised using CNstream analysis software and analysed for association using the χ(2) test. The most significant genomic CNV associations with PsA risk were independently tested in a validation sample of 1133 patients with PsA and 1831 healthy controls. In order to test for the specificity of the variants with PsA aetiology, we also analysed the association to a cohort of 822 patients with purely cutaneous psoriasis (PsC). RESULTS: A total of 165 common CNVs were identified in the genome-wide analysis. We found a highly significant association of an intergenic deletion between ADAMTS9 and MAGI1 genes on chromosome 3p14.1 (p=0.00014). Using the independent patient and control cohort, we validated the association between ADAMTS9-MAGI1 deletion and PsA risk (p=0.032). Using next-generation sequencing, we characterised the 26 kb associated deletion. Finally, analysing the PsC cohort we found a lower frequency of the deletion compared with the PsA cohort (p=0.0088) and a similar frequency to that of healthy controls (p>0.3). CONCLUSIONS: The present genome-wide scan for CNVs associated with PsA risk has identified a new deletion associated with disease risk and which is also differential from PsC risk.

Identification of risk loci for Crohn's disease phenotypes using a genome-wide association study.

BACKGROUND & AIMS: Crohn's disease is a highly heterogeneous inflammatory bowel disease comprising multiple clinical phenotypes. Genome-wide association studies (GWASs) have associated a large number of loci with disease risk but have not associated any specific genetic variants with clinical phenotypes. We performed a GWAS of clinical phenotypes in Crohn's disease. METHODS: We genotyped 576,818 single-nucleotide polymorphisms in a well-characterized cohort of 1090 Crohn's disease patients of European ancestry. We assessed their association with 17 phenotypes of Crohn's disease (based on disease location, disease behavior, disease course, age at onset, and extraintestinal manifestations). A total of 57 markers with strong associations to Crohn's disease phenotypes (P < 2 × 10(-4)) were subsequently analyzed in an independent replication cohort of 1296 patients of European ancestry. RESULTS: We replicated the association of 4 loci with different Crohn's disease phenotypes. Variants in MAGI1, CLCA2, 2q24.1, and LY75 loci were associated with a complicated stricturing disease course (Pcombined = 2.01 × 10(-8)), disease location (Pcombined = 1.3 × 10(-6)), mild disease course (Pcombined = 5.94 × 10(-7)), and erythema nodosum (Pcombined = 2.27 × 10(-6)), respectively. CONCLUSIONS: In a GWAS, we associated 4 loci with clinical phenotypes of Crohn's disease. These findings indicate a genetic basis for the clinical heterogeneity observed for this inflammatory bowel disease.

Effects of cadmium and mercury on the upper part of skeletal muscle glycolysis in mice.

The effects of pre-incubation with mercury (Hg(2+)) and cadmium (Cd(2+)) on the activities of individual glycolytic enzymes, on the flux and on internal metabolite concentrations of the upper part of glycolysis were investigated in mouse muscle extracts. In the range of metal concentrations analysed we found that only hexokinase and phosphofructokinase, the enzymes that shared the control of the flux, were inhibited by Hg(2+) and Cd(2+). The concentrations of the internal metabolites glucose-6-phosphate and fructose-6-phosphate did not change significantly when Hg(2+) and Cd(2+) were added. A mathematical model was constructed to explore the mechanisms of inhibition of Hg(2+) and Cd(2+) on hexokinase and phosphofructokinase. Equations derived from detailed mechanistic models for each inhibition were fitted to the experimental data. In a concentration-dependent manner these equations describe the observed inhibition of enzyme activity. Under the conditions analysed, the integral model showed that the simultaneous inhibition of hexokinase and phosphofructokinase explains the observation that the concentrations of glucose-6-phosphate and fructose-6-phosphate did not change as the heavy metals decreased the glycolytic flux.

A genome-wide association study on a southern European population identifies a new Crohn's disease susceptibility locus at RBX1-EP300.

OBJECTIVE: Genome-wide association studies (GWAS) have identified multiple risk loci for Crohn's disease (CD). However, the cumulative risk exerted by these loci is low, and the likelihood that additional, as-yet undiscovered loci contribute to the risk of CD is very high. We performed a GWAS on a southern European population to identify new CD risk loci. DESIGN: We genotyped 620 901 genome markers on 1341 CD patients and 1518 controls from Spain. The top association signals representing new candidate risk loci were subsequently analysed in an independent replication cohort of 1365 CD patients and 1396 controls. RESULTS: We identified a genome-wide significant association on chromosome 22q13.2 in the intergenic region between the RBX1 and EP300 genes (single nucleotide polymorphism rs4820425, OR 1.27, 95% CI 1.17 to 1.38, p=3.42E-8). We also found suggestive evidence for the association of the IFNGR2 (21q22.11), FOXP2 (7q31), MACROD2 (20p12.1) and AIF1 (6p21.3) loci with CD risk. CONCLUSIONS: In this GWAS performed on a southern European cohort, we have identified a new risk locus for CD between RBX1 and EP300. This study demonstrates that using populations of different ancestry is a useful strategy to identify new risk loci for CD.

Dissection of the recognition properties of p38 MAP kinase. Determination of the binding mode of a new pyridinyl-heterocycle inhibitor family.

The main recognition characteristics of the ATP binding site of p38 mitogen activated protein kinase alpha (p38alpha MAPK) have been explored by a combination of modeling and bioinformatics techniques, making special emphasis in the characteristics of the site that justifies binding specificity with respect to other MAP kinases. Particularly, we have analyzed the binding mode of a new family of p38 MAPK inhibitors based on the pyridinyl-heterocycle core. This family of compounds has a marked pseudosymmetry and can exist in different tautomeric forms, which makes the determination of the binding mode especially challenging. A combination of homology modeling, quantum mechanics, classical docking, and molecular dynamics calculations allowed us to determine the main characteristics defining the binding mode of this new scaffold in the ATP binding site of p38alpha. A set of free energy calculations allowed us to verify the binding mode proposed, giving an overall excellent agreement with the experimental values. Finally, the binding mode of this new family of compounds was compared to that of other members of the pyridinyl and pyrimidinyl heterocycle class.

Linear response theory: an alternative to PB and GB methods for the analysis of molecular dynamics trajectories?

We explore the use of classical Linear Response Theory (LRT) as an alternative strategy to the use of Molecular Mechanics/Poisson-Boltzmann strategies to compute the solvation free energy of macromolecules from molecular dynamics simulations using an explicit representation of solvent. The method reproduces well the free energy of solvation of standard amino acid side chains, small peptides, and proteins. The use of a fully discrete representation of solvent avoids the possible problems of continuum models to represent the solvation of systems containing tightly bound water molecules.