Intermolecular Interactions in G Protein-Coupled Receptor Allosteric Sites at the Membrane Interface from Molecular Dynamics Simulations and Quantum Chemical Calculations.

Allosteric modulators are called promising candidates in G protein-coupled receptor (GPCR) drug development by displaying subtype selectivity and more specific receptor modulation. Among the allosteric sites known to date, cavities at the receptor-lipid interface represent an uncharacteristic binding location that raises many questions about the ligand interactions and stability, the binding site structure, and how all of these are affected by lipid molecules. In this work, we analyze interactions in the allosteric sites of the PAR2, C5aR1, and GCGR receptors in three lipid compositions using molecular dynamics simulations. In addition, we performed quantum chemical calculations involving the symmetry-adapted perturbation theory (SAPT) and the natural population analysis to quantify the strength of intermolecular interactions. We show that besides classical hydrogen bonds, weak polar interactions such as O-HC, O-Br, and long-range electrostatics with the backbone amides contribute to the stability of allosteric modulators at the receptor-lipid interface. The allosteric cavities are detectable in various membrane compositions. The availability of polar atoms for interactions in such cavities can be assessed by water molecules from simulations. Although ligand-lipid interactions are weak, lipid tails play a role in ligand binding pose stability and the size of allosteric cavities. We discuss physicochemical aspects of ligand binding at the receptor-lipid interface and suggest a compound library enriched by weak donor groups for ligand search in such sites.

What Makes a Good Pore Former: A Study of Synthetic Melittin Derivatives.

Pore formation by membrane-active peptides, naturally encountered in innate immunity and infection, could have important medical and technological applications. Recently, the well-studied lytic peptide melittin has formed the basis for the development of combinatorial libraries from which potent pore-forming peptides have been derived, optimized to work under different conditions. We investigate three such peptides, macrolittin70, which is most active at neutral pH; pHD15, which is active only at low pH; and MelP5_Δ6, which was rationally designed to be active at low pH but formed only small pores. There are three, six, and six acidic residues in macrolittin70, pHD15, and MelP5_Δ6, respectively. We perform multi-microsecond simulations in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) of hexamers of these peptides starting from transmembrane orientations at neutral pH (all residues at standard protonation), low pH (acidic residues and His protonated), and highly acidic environments in which C-termini are also protonated. Previous simulations of the parent peptides melittin and MelP5 in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) are repeated in POPC. We find that the most potent pore-forming peptides exhibit strong interpeptide interactions, including salt bridges, H-bonds, and polar interactions. Protonation of the C-terminus promotes helicity and pore size. The proximity of the peptides allows fewer lipid headgroups to line the pores than in previous simulations, making the pores intermediate between barrel stave and toroidal. Based on these structures and geometrical arguments, we attempt to rationalize the factors that under different conditions can increase or decrease pore stability and propose mutations that could be tested experimentally.

Effects of Peptide Charge, Orientation, and Concentration on Melittin Transmembrane Pores.

Melittin is a short cationic peptide that exerts cytolytic effects on bacterial and eukaryotic cells. Experiments suggest that in zwitterionic membranes, melittin forms transmembrane toroidal pores supported by four to eight peptides. A recently constructed melittin variant with a reduced cationic charge, MelP5, is active at 10-fold lower concentrations. In previous work, we performed molecular dynamics simulations on the microsecond timescale to examine the supramolecular pore structure of a melittin tetramer in zwitterionic and partially anionic membranes. We now extend that study to include the effects of peptide charge, initial orientation, and number of monomers on the pore formation and stabilization processes. Our results show that parallel transmembrane orientations of melittin and MelP5 are more consistent with experimental data. Whereas a MelP5 parallel hexamer forms a large stable pore during the 5-µs simulation time, a melittin hexamer and an octamer are not fully stable, with several monomers dissociating during the simulation time. Interaction-energy analysis shows that this difference in behavior between melittin and MelP5 is not due to stronger electrostatic repulsion between neighboring melittin peptides but to peptide-lipid interactions that disfavor the isolated MelP5 transmembrane monomer. The ability of melittin monomers to diffuse freely in the 1,2-dimyristoyl-SN-glycero-3-phosphocholine membrane leads to dynamic pores with varying molecularity.

Pore Structure and Synergy in Antimicrobial Peptides of the Magainin Family.

Magainin 2 and PGLa are among the best-studied cationic antimicrobial peptides. They bind preferentially to negatively charged membranes and apparently cause their disruption by the formation of transmembrane pores, whose detailed structure is still unclear. Here we report the results of 5-9 µs all-atom molecular dynamics simulations starting from tetrameric transmembrane helical bundles of these two peptides, as well as their stoichiometric mixture, and the analog MG-H2 in DMPC or 3:1 DMPC/DMPG membranes. The simulations produce pore structures that appear converged, although some effect of the starting peptide arrangement (parallel vs. antiparallel) is still observed on this timescale. The peptides remain mostly helical and adopt tilted orientations. The calculated tilt angles for PGLa are in excellent agreement with recent solid state NMR experiments. The antiparallel dimer structure in the magainin 2 simulations resembles previously determined NMR and crystal structures. More transmembrane orientations and a larger and more ordered pore are seen in the 1:1 heterotetramer with an antiparallel helix arrangement. Insights into the mechanism of synergy between these two peptides are obtained via implicit solvent modeling of homo- and heterodimers and analysis of interactions in the atomistic simulations. This analysis suggests stronger pairwise interactions in the heterodimer than in the two homodimers.

Structural and functional analogies and differences between histidine decarboxylase and aromatic l-amino acid decarboxylase molecular networks: Biomedical implications.

Human histidine decarboxylase (HDC) and dopa decarboxilase (DDC) are highly homologous enzymes responsible for the synthesis of biogenic amines (BA) like histamine, and serotonin and dopamine, respectively. The enzymes share many structural and functional analogies, while their product metabolisms also follow similar patterns that are confluent in some metabolic steps. They are involved in common physiological functions, such as neurotransmission, gastrointestinal track function, immunity, cell growth and cell differentiation. As a consequence, metabolic elements of both BA subfamilies are also co-participants in a long list of human diseases. This review summarizes the analogies and differences in their origin (HDC and DDC) as well as their common pathophysiological scenarios. The major gaps of information are also underlined, as they delay the possibility of holistic approaches that would help personalized medicine and pharmacological initiatives for prevalent and rare diseases.

The structure of a melittin-stabilized pore.

Melittin has been reported to form toroidal pores under certain conditions, but the atomic-resolution structure of these pores is unknown. A 9-µs all-atom molecular-dynamics simulation starting from a closely packed transmembrane melittin tetramer in DMPC shows formation of a toroidal pore after 1 µs. The pore remains stable with a roughly constant radius for the rest of the simulation. Surprisingly, one or two melittin monomers frequently transition between transmembrane and surface states. All four peptides are largely helical. A simulation in a DMPC/DMPG membrane did not lead to a stable pore, consistent with the experimentally observed lower activity of melittin on anionic membranes. The picture that emerges from this work is rather close to the classical toroidal pore, but more dynamic with respect to the configuration of the peptides.

An Insertion Within SIRPß1 Shows a Dual Effect Over Alzheimer's Disease Cognitive Decline Altering the Microglial Response.

Background: Microglial dysfunction plays a causative role in Alzheimer's disease (AD) pathogenesis. Here we focus on a germline insertion/deletion variant mapping SIRPß1, a surface receptor that triggers amyloid-ß(Aß) phagocytosis via TYROBP. Objective: To analyze the impact of this copy-number variant in SIRPß1 expression and how it affects AD molecular etiology. Methods: Copy-number variant proxy rs2209313 was evaluated in GERALD and GR@ACE longitudinal series. Hippocampal specimens of genotyped AD patients were also examined. SIRPß1 isoform-specific phagocytosis assays were performed in HEK393T cells. Results: The insertion alters the SIRPß1 protein isoform landscape compromising its ability to bind oligomeric Aß and its affinity for TYROBP. SIRPß1 Dup/Dup patients with mild cognitive impairment show an increased cerebrospinal fluid t-Tau/Aß ratio (p = 0.018) and a higher risk to develop AD (OR = 1.678, p = 0.018). MRIs showed that Dup/Dup patients exhibited a worse initial response to AD. At the moment of diagnosis, all patients showed equivalent Mini-Mental State Examination scores. However, AD patients with the duplication had less hippocampal degeneration (p < 0.001) and fewer white matter hyperintensities. In contrast, longitudinal studies indicate that patients bearing the duplication allele show a slower cognitive decline (p = 0.013). Transcriptional analysis also shows that the SIRPß1 duplication allele correlates with higher TREM2 expression and an increased microglial activation. Conclusions: The SIRPß1 internal duplication has opposite effects over MCI-to-Dementia conversion risk and AD progression, affecting microglial response to Aß. Given the pharmacological approaches focused on the TREM2-TYROBP axis, we believe that SIRPß1 structural variant might be considered as a potential modulator of this causative pathway.

Bicyclic derivatives of L-idonojirimycin as pharmacological chaperones for neuronopathic forms of Gaucher disease.

New human ß-glucocerebrosidase (GCase) ligands with rigid 1,6-anhydro-ß-L-idonojirimycin cores have been designed with the aid of molecular modeling. Efficient pharmacological chaperones for the L444P (trafficking-incompetent) mutant GCase enzyme associated with type 2 and 3 Gaucher disease (GD) were identified.

Histamine: an undercover agent in multiple rare diseases?

Histamine is a biogenic amine performing pleiotropic effects in humans, involving tasks within the immune and neuroendocrine systems, neurotransmission, gastric secretion, cell life and death, and development. It is the product of the histidine decarboxylase activity, and its effects are mainly mediated through four different G-protein coupled receptors. Thus, histamine-related effects are the results of highly interconnected and tissue-specific signalling networks. Consequently, alterations in histamine-related factors could be an important part in the cause of multiple rare/orphan diseases. Bearing this hypothesis in mind, more than 25 rare diseases related to histamine physiopathology have been identified using a computationally assisted text mining approach. These newly integrated data will provide insight to elucidate the molecular causes of these rare diseases. The data can also help in devising new intervention strategies for personalized medicine for multiple rare diseases.

Structural perspective on the direct inhibition mechanism of EGCG on mammalian histidine decarboxylase and DOPA decarboxylase.

Histidine decarboxylase (HDC) and l-aromatic amino acid decarboxylase (DDC) are homologous enzymes that are responsible for the synthesis of important neuroactive amines related to inflammatory, neurodegenerative, and neoplastic diseases. Epigallocatechin-3-gallate (EGCG), the most abundant catechin in green tea, has been shown to target histamine-producing cells and to promote anti-inflammatory, antitumor, and antiangiogenic effects. Previous experimental work has demonstrated that EGCG has a direct inhibitory effect on both HDC and DDC. In this study, we investigated the binding modes of EGCG to HDC and DDC as a first step for designing new polyphenol-based HDC/DDC-specific inhibitors.

Histamine, Metabolic Remodelling and Angiogenesis: A Systems Level Approach.

Histamine is a highly pleiotropic biogenic amine involved in key physiological processes including neurotransmission, immune response, nutrition, and cell growth and differentiation. Its effects, sometimes contradictory, are mediated by at least four different G-protein coupled receptors, which expression and signalling pathways are tissue-specific. Histamine metabolism conforms a very complex network that connect many metabolic processes important for homeostasis, including nitrogen and energy metabolism. This review brings together and analyses the current information on the relationships of the "histamine system" with other important metabolic modules in human physiology, aiming to bridge current information gaps. In this regard, the molecular characterization of the role of histamine in the modulation of angiogenesis-mediated processes, such as cancer, makes a promising research field for future biomedical advances.

AMMO-Prot: amine system project 3D-model finder.

BACKGROUND: Amines are biogenic amino acid derivatives, which play pleiotropic and very important yet complex roles in animal physiology. For many other relevant biomolecules, biochemical and molecular data are being accumulated, which need to be integrated in order to be effective in the advance of biological knowledge in the field. For this purpose, a multidisciplinary group has started an ontology-based system named the Amine System Project (ASP) for which amine-related information is the validation bench. RESULTS: In this paper, we describe the Ontology-Based Mediator developed in the Amine System Project (http://asp.uma.es) using the infrastructure of Semantic Directories, and how this system has been used to solve a case related to amine metabolism-related protein structures. CONCLUSIONS: This infrastructure is used to publish and manage not only ontologies and their relationships, but also metadata relating to the resources committed with the ontologies. The system developed is available at http://asp.uma.es/WebMediator.

Transmembrane Pore Structures of ß-Hairpin Antimicrobial Peptides by All-Atom Simulations.

Protegrin-1 is an 18-residue ß-hairpin antimicrobial peptide (AMP) that has been suggested to form transmembrane ß-barrels in biological membranes. However, alternative structures have also been proposed. Here, we performed multimicrosecond, all-atom molecular dynamics simulations of various protegrin-1 oligomers on the membrane surface and in transmembrane topologies. The membrane surface simulations indicated that protegrin dimers are stable, whereas trimers and tetramers break down. Tetrameric arcs remained stably inserted in lipid membranes, but the pore water was displaced by lipid molecules. Unsheared protegrin ß-barrels opened into ß-sheets that surrounded stable aqueous pores, whereas tilted barrels with sheared hydrogen bonding patterns were stable in most topologies. A third type of observed pore consisted of multiple small oligomers surrounding a small, partially lipidic pore. We also considered the ß-hairpin AMP tachyplesin, which showed less tendency to oligomerize than protegrin: the octameric bundle resulted in small pores surrounded by six peptides as monomers and dimers, with some peptides returning to the membrane surface. The results imply that multiple configurations of protegrin oligomers may produce aqueous pores and illustrate the relationship between topology and putative steps in protegrin-1's pore formation. However, the long-term stability of these structures needs to be assessed further.

The noni anthraquinone damnacanthal is a multi-kinase inhibitor with potent anti-angiogenic effects.

The natural bioactive compound damnacanthal inhibits several tyrosine kinases. Herein, we show that -in fact- damancanthal is a multi kinase inhibitor. A docking and molecular dynamics simulation approach allows getting further insight on the inhibitory effect of damnacanthal on three different kinases: vascular endothelial growth factor receptor-2, c-Met and focal adhesion kinase. Several of the kinases targeted and inhibited by damnacanthal are involved in angiogenesis. Ex vivo and in vivo experiments clearly demonstrate that, indeed, damnacanthal is a very potent inhibitor of angiogenesis. A number of in vitro assays contribute to determine the specific effects of damnacanthal on each of the steps of the angiogenic process, including inhibition of tubulogenesis, endothelial cell proliferation, survival, migration and production of extracellular matrix remodeling enzyme. Taken altogether, these results suggest that damancanthal could have potential interest for the treatment of cancer and other angiogenesis-dependent diseases.

Candidate gene study of TRAIL and TRAIL receptors: association with response to interferon beta therapy in multiple sclerosis patients.

TRAIL and TRAIL Receptor genes have been implicated in Multiple Sclerosis pathology as well as in the response to IFN beta therapy. The objective of our study was to evaluate the association of these genes in relation to the age at disease onset (AAO) and to the clinical response upon IFN beta treatment in Spanish MS patients. We carried out a candidate gene study of TRAIL, TRAILR-1, TRAILR-2, TRAILR-3 and TRAILR-4 genes. A total of 54 SNPs were analysed in 509 MS patients under IFN beta treatment, and an additional cohort of 226 MS patients was used to validate the results. Associations of rs1047275 in TRAILR-2 and rs7011559 in TRAILR-4 genes with AAO under an additive model did not withstand Bonferroni correction. In contrast, patients with the TRAILR-1 rs20576-CC genotype showed a better clinical response to IFN beta therapy compared with patients carrying the A-allele (recessive model: p = 8.88×10(-4), pc = 0.048, OR = 0.30). This SNP resulted in a non synonymous substitution of Glutamic acid to Alanine in position 228 (E228A), a change previously associated with susceptibility to different cancer types and risk of metastases, suggesting a lack of functionality of TRAILR-1. In order to unravel how this amino acid change in TRAILR-1 would affect to death signal, we performed a molecular modelling with both alleles. Neither TRAIL binding sites in the receptor nor the expression levels of TRAILR-1 in peripheral blood mononuclear cell subsets (monocytes, CD4+ and CD8+ T cells) were modified, suggesting that this SNP may be altering the death signal by some other mechanism. These findings show a role for TRAILR-1 gene variations in the clinical outcome of IFN beta therapy that might have relevance as a biomarker to predict the response to IFN beta in MS.