In silico characterization of ligand-receptor interactions for U-47700, N,N-didesmethyl-U-47700, U-50488 at mu- and kappa-opioid receptors.

The increasing misuse of novel synthetic opioids (NSOs) represents a serious public health concern. In this regard, U-47700 (trans-3,4-dichloro-N-[2-(dimethylamino)cyclohexyl]-N-methylbenzamide) and related "U-compounds" emerged on recreational drug markets as synthetic substitutes for illicit heroin and constituents of counterfeit pain medications. While the pharmacology of U-compounds has been investigated using in vitro and in vivo methods, there is still a lack of understanding about the details of ligand-receptor interactions at the molecular level. To this end, we have developed a molecular modeling protocol based on docking and molecular dynamics simulations to assess the nature of ligand-receptor interactions for U-47700, N,N-didesmethyl U-47700, and U-50488 at the mu-opioid receptor (MOR) and kappa-opioid receptor (KOR). The evaluation of ligand-receptor and ligand-receptor-membrane interaction energies enabled the identification of subtle conformational shifts in the receptors induced by ligand binding. Interestingly, the removal of two key methyl groups from U-47700, to form N,N-didesmethyl U-47700, caused a loss of hydrogen bond contact with tryptophan (Trp)229, which may underlie the lower interaction energy and reduced MOR affinity for the compound. Taken together, our results are consistent with the reported biological findings for U-compounds and provide a molecular basis for the MOR selectivity of U-47700 and KOR selectivity of U-50488.

The impact of missense mutation in PIGA associated to paroxysmal nocturnal hemoglobinuria and multiple congenital anomalies-hypotonia-seizures syndrome 2: A computational study.

Paroxysmal nocturnal hemoglobinuria (PNH) is an acquired clonal blood disorder that manifests with hemolytic anemia, thrombosis, and peripheral blood cytopenias. The disease is caused by the deficiency of two glycosylphosphatidylinositols (GPI)-anchored proteins (CD55 and CD59) in the hemopoietic stem cells. The deficiency of GPI-anchored proteins has been associated with the somatic mutations in phosphatidylinositol glycan class A (PIGA). However, the mutations that do not cause PNH is associated with the multiple congenital anomalies-hypotonia-seizures syndrome 2 (MCAHS2). To best of our knowledge, no computational study has been performed to explore at an atomistic level the impact of PIGA missense mutations on the structure and dynamics of the protein. Therefore, we focused our study to provide molecular insights into the changes in protein structural dynamics upon mutation. In the initial step, screening for the most pathogenic mutations from the pool of publicly available mutations was performed. Further, to get a better understanding, pathogenic mutations were mapped to the modeled structure and the resulting protein was subjected to 100 ns molecular dynamics simulation. The residues close to C- and N-terminal regions of the protein were found to exhibit greater flexibility upon mutation. Our study suggests that four mutations are highly effective in altering the structural conformation and stability of the PIGA protein. Among them, mutant G48D was found to alter protein's structural dynamics to the greatest extent, both on a local and a global scale.

Novel 2-pheynlbenzofuran derivatives as selective butyrylcholinesterase inhibitors for Alzheimer's disease.

Alzheimer's disease (AD) is a neurodegenerative disorder representing the leading cause of dementia and is affecting nearly 44 million people worldwide. AD is characterized by a progressive decline in acetylcholine levels in the cholinergic systems, which results in severe memory loss and cognitive impairments. Expression levels and activity of butyrylcholinesterase (BChE) enzyme has been noted to increase significantly in the late stages of AD, thus making it a viable drug target. A series of hydroxylated 2-phenylbenzofurans compounds were designed, synthesized and their inhibitory activities toward acetylcholinesterase (AChE) and BChE enzymes were evaluated. Two compounds (15 and 17) displayed higher inhibitory activity towards BChE with IC50 values of 6.23 µM and 3.57 µM, and a good antioxidant activity with EC50 values 14.9 µM and 16.7 µM, respectively. The same compounds further exhibited selective inhibitory activity against BChE over AChE. Computational studies were used to compare protein-binding pockets and evaluate the interaction fingerprints of the compound. Molecular simulations showed a conserved protein residue interaction network between the compounds, resulting in similar interaction energy values. Thus, combination of biochemical and computational approaches could represent rational guidelines for further structural modification of these hydroxy-benzofuran derivatives as future drugs for treatment of AD.

2-Phenylbenzofuran derivatives as butyrylcholinesterase inhibitors: Synthesis, biological activity and molecular modeling.

A series of 2-phenylbenzofurans compounds was designed, synthesized and evaluated as cholinesterase inhibitors. The biological assay experiments showed that most of the compounds displayed a clearly selective inhibition for butyrylcholinesterase (BChE), while a weak or no effect towards acetylcholinesterase (AChE) was detected. Among these benzofuran derivatives, compound 16 exhibited the highest BChE inhibition with an IC50 value of 30.3 µM. This compound was found to be a mixed-type inhibitor as determined by kinetic analysis. Moreover, molecular dynamics simulations revealed that compound 16 binds to both the catalytic anionic site (CAS) and peripheral anionic site (PAS) of BChE and it displayed the best interaction energy value, in agreement with our experimental data.

Metabolomics analysis and modeling suggest a lysophosphocholines-PAF receptor interaction in fibromyalgia.

Fibromyalgia Syndrome (FMS) is a chronic disease characterized by widespread pain, and difficult to diagnose and treat. We analyzed the plasma metabolic profile of patients with FMS by using a metabolomics approach combining Liquid Chromatography-Quadrupole-Time Of Flight/Mass Spectrometry (LC-Q-TOF/MS) with multivariate statistical analysis, aiming to discriminate patients and controls. LC-Q-TOF/MS analysis of plasma (FMS patients: n = 22 and controls: n = 21) identified many lipid compounds, mainly lysophosphocholines (lysoPCs), phosphocholines and ceramides. Multivariate statistical analysis was performed to identify the discriminating metabolites. A protein docking and molecular dynamic (MD) study was then performed, using the most discriminating lysoPCs, to validate the binding to Platelet Activating Factor (1-alkyl-2-acetyl-sn-glycero-3-phosphocholine, PAF) Receptor (PAFr). Discriminating metabolites between FMS patients and controls were identified as 1-tetradecanoyl-sn-glycero-3-phosphocholine [PC(14:0/0:0)] and 1-hexadecanoyl-sn-glycero-3-phosphocholine [PC(16:0/0:0)]. MD and docking indicate that the ligands investigated have similar potentialities to activate the PAFr receptor. The application of a metabolomic approach discriminated FMS patients from controls, with an over-representation of PC(14:0/0:0) and PC(16:0/0:0) compounds in the metabolic profiles. These results and the modeling of metabolite-PAFr interaction, allowed us to hypothesize that lipids oxidative fragmentation might generate lysoPCs in abundance, that in turn will act as PAF-like bioactivators. Overall results suggest disease biomarkers and potential therapeutical targets for FMS.

Antigenic peptide molecular recognition by the DRB1-DQB1 haplotype modulates multiple sclerosis susceptibility.

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system that has a notably high incidence in Sardinia. Our study focuses on two HLA class II haplotypes associated with the disease in Sardinia, the rare predisposing DRB1*15:01-DQB1*06:02 and the widespread protective DRB1*16:01-DQB1*05:02. This framework enabled the highlighting of HLA binding pocket specificity and peptide recognition mechanisms by employing molecular dynamics simulations of the whole DRB1-DQB1 haplotype interacting with MBP- and EBV-derived peptides. We analyzed peptide-protein interaction networks and temporal evolution of the original complexes and after key amino acid mutations. The mutation G86V of the protective DRB1 allele exerted its effect mainly in the presence of the EBV viral peptide, with local and long range outcomes. However, the V38A mutation of the protective DQB1 showed a long range effect only in the case of the MBP myelin peptide. Our findings also demonstrate a DRB1/DQB1 complementary molecular recognition of peptides. This mechanism could provide a robust synergistic action and a differential role of DRB1 and DQB1 in tissues and in the time-steps towards autoimmunity. In addition, we demonstrate that negatively charged residues in pockets 4 and 9 play a role in MS susceptibility. Our findings are supported by recent experiments using a closely related MS animal model. Overall, our analysis confirms the role of the DRB1-DQB1 haplotype in conferring disease predisposition and could provide a valuable aid in designing optimal therapeutic peptides for MS therapy.

Identification of calcium binding sites on calsequestrin 1 and their implications for polymerization.

Biophysical studies have shown that each molecule of calsequestrin 1 (CASQ1) can bind about 70-80 Ca(2+) ions. However, the nature of Ca(2+)-binding sites has not yet been fully characterized. In this study, we employed in silico approaches to identify the Ca(2+) binding sites and to understand the molecular basis of CASQ1-Ca(2+) recognition. We built the protein model by extracting the atomic coordinates for the back-to-back dimeric unit from the recently solved hexameric CASQ1 structure (PDB id: ) and adding the missing C-terminal residues (aa350-364). Using this model we performed extensive 30 ns molecular dynamics simulations over a wide range of Ca(2+) concentrations ([Ca(2+)]). Our results show that the Ca(2+)-binding sites on CASQ1 differ both in affinity and geometry. The high affinity Ca(2+)-binding sites share a similar geometry and interestingly, the majority of them were found to be induced by increased [Ca(2+)]. We also found that the system shows maximal Ca(2+)-binding to the CAS (consecutive aspartate stretch at the C-terminus) before the rest of the CASQ1 surface becomes saturated. Simulated data show that the CASQ1 back-to-back stacking is progressively stabilized by the emergence of an increasing number of hydrophobic interactions with increasing [Ca(2+)]. Further, this study shows that the CAS domain assumes a compact structure with an increase in Ca(2+) binding, which suggests that the CAS domain might function as a Ca(2+)-sensor that may be a novel structural motif to sense metal. We propose the term "Dn-motif" for the CAS domain.

Structural and dynamical insights on HLA-DR2 complexes that confer susceptibility to multiple sclerosis in Sardinia: a molecular dynamics simulation study.

Sardinia is a major Island in the Mediterranean with a high incidence of multiple sclerosis, a chronic autoimmune inflammatory disease of the central nervous system. Disease susceptibility in Sardinian population has been associated with five alleles of major histocompatibility complex (MHC) class II DRB1 gene. We performed 120 ns of molecular dynamics simulation on one predisposing and one protective alleles, unbound and in complex with the two relevant peptides: Myelin Basic Protein and Epstein Barr Virus derived peptide. In particular we focused on the MHC peptide binding groove dynamics. The predisposing allele was found to form a stable complex with both the peptides, while the protective allele displayed stability only when bound with myelin peptide. The local flexibility of the MHC was probed dividing the binding groove into four compartments covering the well known peptide anchoring pockets. The predisposing allele in the first half cleft exhibits a narrower and more rigid groove conformation in the presence of myelin peptide. The protective allele shows a similar behavior, while in the second half cleft it displays a narrower and more flexible groove conformation in the presence of viral peptide. We further characterized these dynamical differences by evaluating H-bonds, hydrophobic and stacking interaction networks, finding striking similarities with super-type patterns emerging in other autoimmune diseases. The protective allele shows a defined preferential binding to myelin peptide, as confirmed by binding free energy calculations. All together, we believe the presented molecular analysis could help to design experimental assays, supports the molecular mimicry hypothesis and suggests that propensity to multiple sclerosis in Sardinia could be partly linked to distinct peptide-MHC interaction and binding characteristics of the antigen presentation mechanism.

Interaction between HLA-DRB1-DQB1 haplotypes in Sardinian multiple sclerosis population.

We performed a case-control study in 2,555 multiple sclerosis (MS) Sardinian patients and 1,365 healthy ethnically matched controls, analyzing the interactions between HLA-DRB1-DQB1 haplotypes and defining a rank of genotypes conferring a variable degree of risk to the disease. Four haplotypes were found to confer susceptibility (*13:03-*03:01 OR = 3.3, Pc 5.1 × 10(-5), *04:05-*03:01 OR = 2.1, Pc 9.7 × 10(-8), *15:01-*06:02 OR = 2.0, Pc = 9.1 × 10(-3), *03:01-*02:01 OR = 1.7 Pc = 7.9 × 10(-22)) and protection (*11, OR = 0.8, Pc = 2.7 × 10(-2), *16:01-*05:02 OR = 0.6, Pc = 4.8 × 10(-16), *14:01-4-*05:031 = OR = 0.5, Pc = 9.8 × 10(-4) and *15:02-*06:01 OR = 0.4, Pc = 5.1 × 10(-4)). The relative predispositional effect method confirms all the positively associated haplotypes and showed that also *08 and *04 haplotypes confers susceptibility, while the *11 was excluded as protective haplotype. Genotypic ORs highlighted two typologies of interaction between haplotypes: i) a neutral interaction, in which the global risk is coherent with the sum of the single haplotype risks; ii) a negative interaction, in which the genotypic OR observed is lower than the sum of the OR of the two haplotypes. The phylogenic tree of the MS-associated DRB1 alleles found in Sardinian patients revealed a cluster represented by *14:01, *04:05, *13∶03, *08:01 and *03:01 alleles. Sequence alignment analysis showed that amino acids near pocket P4 and pocket P9 differentiated protective from predisposing alleles under investigation. Furthermore, molecular dynamics simulation performed on alleles revealed that position 70 is crucial in binding of MBP 85-99 peptide. All together, these data suggest that propensity to MS observed in Sardinian population carried by the various HLA-DRB1-DQB1 molecules can be due to functional peculiarity in the antigen presentation mechanisms.

HLA-DRB1-DQB1 haplotypes confer susceptibility and resistance to multiple sclerosis in Sardinia.

INTRODUCTION: Genetic predisposition to multiple sclerosis (MS) in Sardinia (Italy) has been associated with five DRB1*-DQB1* haplotypes of the human leukocyte antigen (HLA). Given the complexity of these associations, an in-depth re-analysis was performed with the specific aims of confirming the haplotype associations; establishing the independence of the associated haplotypes; and assessing patients' genotypic risk of developing MS. METHODS AND RESULTS: A transmission disequilibrium test (TDT) of the DRB1*-DQB1* haplotypes in 943 trio families, confirmed a higher than expected transmission rate (over-transmission) of the *13:03-*03:01 (OR = 2.9, P = 7.6×10(-3)), *04:05-*03:01 (OR = 2.4, P = 4.4×10(-6)) and *03:01-*02:01 (OR = 2.1, P = 1.0×10(-15)) haplotype. In contrast, the *16:01-*05:02 (OR = 0.5, P = 5.4×10(-11)) and the *15:02-*06:01 (OR = 0.3, P = 1.5×10(-3)) haplotypes exhibited a lower than expected transmission rate (under-transmission). The independence of the transmission of each positively and negatively associated haplotype was confirmed relative to all positively associated haplotypes, and to the negatively associated *16:01-*05:02 haplotype. In patients, carriage of two predisposing haplotypes, or of protective haplotypes, respectively increased or decreased the patient's risk of developing MS. The risk of MS followed a multiplicative model of genotypes, which was, in order of decreasing ORs: *04:05-*0301/*03:01-*02:01 (OR = 4.5); *03:01-*02:01/*03:01-*02:01 (OR = 4.1); and the *16:01-*05:02/*16:01-*0502 (OR = 0.2) genotypes. Analysis of DRB1 and DQB1 protein chain residues showed that the Val/Gly residue at position 86 of the DRB1 chain was the only difference between the protective *16:01- *15:02 alleles and the predisposing *15:01 one. Similarly, the Ala/Val residue at position 38 of the DQB1 chain differentiated the positively associated *06:02 allele and the negatively associated *05:02, *06:01 alleles. CONCLUSIONS: These findings show that the association of specific, independent DRB1*-DQB1* haplotypes confers susceptibility or resistance to MS in the MS-prone Sardinian population. The data also supports a functional role for specific residues of the DRB1 and DQB1 proteins in predisposing patients to MS.

Inferring gene networks: dream or nightmare?

We describe several algorithms with winning performance in the Dialogue for Reverse Engineering Assessments and Methods (DREAM2) Reverse Engineering Competition 2007. After the gold standards for the challenges were released, the performance of the algorithms could be thoroughly evaluated under different parameters or alternative ways of solving systems of equations. For the analysis of Challenge 4, the "In-silico" challenges, we employed methods to explicitly deal with perturbation data and time-series data. We show that original methods used to produce winning submissions could easily be altered to substantially improve performance. For Challenge 5, the genome-scale Escherichia coli network, we evaluated a variety of measures of association. These data are troublesome, and no good solutions could be produced, either by us or by any other teams. Our best results were obtained when analyzing subdatasets instead of considering the dataset as a whole.

Protein networking: insights into global functional organization of proteomes.

The formulation of network models from global protein studies is essential to understand the functioning of organisms. Network models of the proteome enable the application of Complex Network Analysis, a quantitative framework to investigate large complex networks using techniques from graph theory, statistical physics, dynamical systems and other fields. This approach has provided many insights into the functional organization of the proteome so far and will likely continue to do so. Currently, several network concepts have emerged in the field of proteomics. It is important to highlight the differences between these concepts, since different representations allow different insights into functional organization. One such concept is the protein interaction network, which contains proteins as nodes and undirected edges representing the occurrence of binding in large-scale protein-protein interaction studies. A second concept is the protein-signaling network, in which the nodes correspond to levels of post-translationally modified forms of proteins and directed edges to causal effects through post-translational modification, such as phosphorylation. Several other network concepts were introduced for proteomics. Although all formulated as networks, the concepts represent widely different physical systems. Therefore caution should be taken when applying relevant topological analysis. We review recent literature formulating and analyzing such networks.