The second intracellular loop of the calcitonin gene-related peptide receptor provides molecular determinants for signal transduction and cell surface expression

The calcitonin gene-related peptide (CGRP) receptor is a heterodimer of a family B G-protein-coupled receptor, calcitonin receptor-like receptor (CLR), and the accessory protein receptor activity modifying protein 1. It couples to Gs, but it is not known which intracellular loops mediate this. We have identified the boundaries of this loop based on the relative position and length of the juxtamembrane transmembrane regions 3 and 4. The loop has been analyzed by systematic mutagenesis of all residues to alanine, measuring cAMP accumulation, CGRP affinity, and receptor expression. Unlike rhodopsin, ICL2 of the CGRP receptor plays a part in the conformational switch after agonist interaction. His-216 and Lys-227 were essential for a functional CGRP-induced cAMP response. The effect of (H216A)CLR is due to a disruption to the cell surface transport or surface stability of the mutant receptor. In contrast, (K227A)CLR had wild-type expression and agonist affinity, suggesting a direct disruption to the downstream signal transduction mechanism of the CGRP receptor. Modeling suggests that the loop undergoes a significant shift in position during receptor activation, exposing a potential G-protein binding pocket. Lys-227 changes position to point into the pocket, potentially allowing it to interact with bound G-proteins. His-216 occupies a position similar to that of Tyr-136 in bovine rhodopsin, part of the DRY motif of the latter receptor. This is the first comprehensive analysis of an entire intracellular loop within the calcitonin family of G-protein-coupled receptor. These data help to define the structural and functional characteristics of the CGRP-receptor and of family B G-protein-coupled receptors in general. © 2006 by The American Society for Biochemistry and Molecular Biology, Inc.

HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.

Some aspects of the interaction of humic substances with metal ions

Humic substances are the major organic constituents of soils and sediments. They are heterogeneous, polyfunctional, polydisperse, macromolecular and have no accurately known chemical structure. Their interactions with radionuclides are particularly important since they provide leaching mechanisms from disposal sites. The central theme to this research is the interaction of heavy metal actinide analogues with humic materials. Studies described focus on selected aspects of the characteristics and properties of humic substances. Some novel approaches to experiments and data analysis are pursued. Several humic substances are studied; all but one are humic acids, and those used most extensively were obtained commercially. Some routine characterisation techniques are applied to samples in the first instance. Humic substances are coloured, but their ultra-violet and visible absorption spectra are featureless. Yet, they fluoresce over a wide range of wavelengths. Enhanced fluorescence in the presence of luminescent europium(III) ions is explained by energy transfer from irradiated humic acid to the metal ion in a photophysical model. Nuclear magnetic resonance spectroscopy is applied to the study of humic acids and their complexes with heavy metals. Proton and carbon-13 NMR provides some structural and functionality information; Paramagnetic lanthanide ions affect these spectra. Some heavy metals are studied as NMR nuclei, but measurements are restricted by their sensitivity. A humic acid is fractionated yielding a broad molecular weight distribution. Electrophoretic mobilities and particle radii determined by Laser Doppler Electrophoretic Light Scattering are sensitive to the conditions of the supporting media, and the concentration and particle size distribution of humic substances. In potentiometric titrations of humate dispersions, the organic matter responds slowly and the mineral acid addition is buffered. Proton concentration data is modelled and a mechanism is proposed involving two key stages, both resulting in proton release after some conformational changes.

Interaction of peptides, peptidomimetics and other drug candidates with the DI/tripeptide transport system in intestinal epithelial cells, using the in vitro caco-2 cell culture system

An uptake system was developed using Caco-2 cell monolayers and the dipeptide, glycyl-[3H]L-proline, as a probe compound. Glycyl-[3H]L-proline uptake was via the di-/tripeptide transport system (DTS) and, exhibited concentration-, pH- and temperature-dependency. Dipeptides inhibited uptake of the probe, and the design of the system allowed competitors to be ranked against one another with respect to affinity for the transporter. The structural features required to ensure or increase interaction with the DTS were defined by studying the effect of a series of glycyl-L-proline and angiotensin-converting enzyme (ACE)-inhibitor (SQ-29852) analogues on the uptake of the probe. The SQ-29852 structure was divided into six domains (A-F) and competitors were grouped into series depending on structural variations within specific regions. Domain A was found to prefer a hydrophobic function, such as a phenyl group, and was intolerant to positive charges and H+ -acceptors and donors. SQ-29852 analogues were more tolerant of substitutions in the C domain, compared to glycyl-L-proline analogues, suggesting that interactions along the length of the SQ-29852 molecule may override the effects of substitutions in the C domain. SQ-29852 analogues showed a preference for a positive function, such as an amine group in this region, but dipeptide structures favoured an uncharged substitution. Lipophilic substituents in domain D increased affinity of SQ-29852 analogues with the DTS. A similar effect was observed for ACE-NEP inhibitor analogues. Domain E, corresponding to the carboxyl group was found to be tolerant of esterification for SQ-29852 analogues but not for dipeptides. Structural features which may increase interaction for one series of compounds, may not have the same effect for another series, indicating that the presence of multiple recognition sites on a molecule may override the deleterious effect of anyone change. Modifying current, poorly absorbed peptidomimetic structures to fit the proposed hypothetical model may improve oral bioavailability by increasing affinity for the DTS. The stereochemical preference of the transporter was explored using four series of compounds (SQ-29852, lysylproline, alanylproline and alanylalanine enantiomers). The L, L stereochemistry was the preferred conformation for all four series, agreeing with previous studies. However, D, D enantiomers were shown in some cases to be substrates for the DTS, although exhibiting a lower affinity than their L, L counterparts. All the ACE-inhibitors and β-lactam antibiotics investigated, produced a degree of inhibition of the probe, and thus show some affinity for the DTS. This contrasts with previous reports that found several ACE inhibitors to be absorbed via a passive process, thus suggesting that compounds are capable of binding to the transporter site and inhibiting the probe without being translocated into the cell. This was also shown to be the case for oligodeoxynucleotide conjugated to a lipophilic group (vitamin E), and highlights the possibility that other orally administered drug candidates may exert non-specific effects on the DTS and possibly have a nutritional impact. Molecular modelling of selected ACE-NEP inhibitors revealed that the three carbonyl functions can be oriented in a similar direction, and this conformation was found to exist in a local energy-minimised state, indicating that the carbonyls may possibly be involved in hydrogen-bond formation with the binding site of the DTS.

Interaction of metal salts with polyethers

The dielectric properties of pure low to medium molecular weight poly(ethylene glycol) and poly(propylene glycol) and a variety of their salt complexes have been studied through the measurement of the dielectric permittivity and dielectric loss over a range of frequency and temperature. The major proportion of this study has been concerned with the examination of the nature of the interaction between mercuric chloride and poly(propylene glycol) (PPG). Other salt-poly-ether combinations have also been considered such as cobalt chloride-PPG cadmium chloride-PPG zinc chloride-PPG and ferric chloride-PEG (polyethylene glycol). Some of this work was also supported by chemical shift and spin-lattice Nuclear Magnetic Resonance (N.M.R.) spectroscopy. The dielectric permittivity data were analysed using the Onsager relation to calculate the mean dipole moment per dipolar unit. This approach was employed in the discussion of various models proposed for the structure of salt-polyether complexes. The effect of mercuric chloride on the statistical conformations of poly(propylene-glycol) was studied in a quantitative manner using the relationships of marchal-Benoit. The dielectric relaxation activation energy and mean energy difference between gauche and trans conformations of poly(propylene glycol) in the presence of mercuric chloride, both showed a distinct minimum when the concentration of mercuric chloride was close to 5 mole %. Opposite behaviour was observed for the Cole-Cole parameter. It was concluded that the majority of the dielectric data could be rationalised in terms of a 5-membered cyclic complex formed between mercuric chloride and PPG in which the complexed segment of the polyether-(OMeCH2CH2O)- adopted either gauche or cis conformations.

Studies of molecular interactions using physical techniques

The investigations described in this thesis concern the molecular interactions between polar solute molecules and various aromatic compounds in solution. Three different physical methods were employed. Nuclear magnetic resonance (n.m.r.) spectroscopy was used to determine the nature and strength of the interactions and the geometry of the transient complexes formed. Cryoscopic studies were used to provide information on the stoichiometry of the complexes. Dielectric constant studies were conducted in an attempt to confirm and supplement the spectroscopic investigations. The systems studied were those between nitromethane, chloroform, acetonitrile (solutes) and various methyl substituted benzenes. In the n.m.r. work the dependence of the solute chemical shift upon the compositions of the solutions was determined. From this the equilibrium quotients (K) for the formation of each complex and the shift induced in the solute proton by the aromatic in the complex were evaluated. The thermodynamic parameters for the interactions were obtained from the determination of K at several temperatures. The stoichiometries of the complexes obtained from cryoscopic studies were found to agree with those deduced from spectroscopic investigations. For most systems it is suggested that only one type of complex, of 1:1 stiochiometry, predominates except that for the acetonitrile-benzene system a 1:2 complex is formed. Two sets of dielectric studies were conducted, the first to show that the nature of the interaction is dipole-induced dipole and the second to calculate K. The equilibrium quotients obtained from spectroscopic and dielectric studies are compared. Time-averaged geometries of the complexes are proposed. The orientation of solute, with respect to the aromatic for the 1:1 complexes, appears to be the one in which the solute lies symmetrically about the aromatic six-fold axis whereas for the 1:2 complex, a sandwich structure is proposed. It is suggested that the complexes are formed through a dipole-induced dipole interaction and steric factors play some part in the complex formation.

Molecular regulation of iron uptake in pseudomonas aeruginosa

The microbial demand for iron is often met by the elaboration of siderophores into the surrounding medium and expression of cognate outer membrane receptors for the ferric siderophore complexes. Conditions of iron limitation, such as those encountered in vivo, cause Pseudomonas aeruginosa to express two high-affinity iron-uptake systems based on pyoverdin and pyochelin. These systems will operate both in the organism's natural habitat, soil and water, where the solubility of iron at neutral pH is extremely low, and in the human host where the availability of free iron is too low to sustain bacterial growth due to the iron-binding glycoproteins transferrin and lactoferrin. Cross-feeding and radiolabelled iron uptake experiments demonstrated that pyoverdin biosynthesis and uptake were highly heterogeneous amongst P.aeruginosa strains, that growth either in the presence of pyoverdin or pyochelin resulted in induction of specific IROMPs, and that induction of iron uptake is siderophore-specific. The P.aeruginosa Tn5 mutant PH1 is deficient in ferripyoverdin uptake and resistant to pyocin Sa, suggesting that the site of interaction of pyocin Sa is a ferripyoverdin receptor. Additional Tn5 mutants appeared to exploit different strategies to achieve pyocin Sa-resistance, involving modifications in expression of pyoverdin-mediated iron uptake, indicating that complex regulatory systems exist to enable these organisms to compete effectively for iron. Modulation of expression of IROMPs prompted a study of the mechanism of uptake of a semi-synthetic C(7) α-formamido substituted cephalosporin BRL 41897A. Sensitivity to this agent correlated with expression of the 75 kDa ferri-pyochelin receptor and demonstrated the potential of high-affinity iron uptake systems for targeting of novel antibiotics. Studies with ferri-pyoverdin uptake-deficient mutant PH1 indicated that expression of outer membrane protein G (OprG), which is usually expressed under iron-rich conditions and repressed under iron-deficient conditions, was perturbed. Attempts were made to clone the oprG gene using a degenerate probe based on the N-terminal amino acid sequence. A strongly hybridising HindIll restriction fragment was cloned and sequenced, but failed to reveal an open reading frame correspondmg to OprG. However, there appears to be good evidence that a part of the gene codmg for the hydrophilic membrane-associated ATP-binding component of a hitherto uncharacterised periplasmic- binding-protein-dependent transport system has been isolated. The full organisation and sequence of the operon, and substrate for this putative transport system, are yet: to be elucidated,

Molecular architecture influences on material properties of pharmaceutical compounds

Salt formation has extensively been studied as a strategy to improve drug solubility but it has not been explored as a strategy to improve mechanical properties. A better understanding of which factors of the solid state can have an influence in the mechanical properties of pharmaceutical powders can help to optimise and reduce cost of tablet manufacturing. The aim of this study was to form different series of amine salts of flurbiprofen, gemfibrozil and diclofenac and to establish predictive relationships between architectural characteristics and physicochemical and mechanical properties of the salts. For this purpose, three different carboxylic acid drugs were selected: flurbiprofen, gemfibrozil and diclofenac, similar in size but varying in flexibility and shape and three different series of counterions were also chosen: one with increasing bulk and no hydroxyl groups to limit the hydrogen bonding potential; a second one with increasing number of hydroxyl groups and finally a third series, related to the latter in number of hydroxyl groups but with different molecular shape and flexibility. Physico-chemical characterization was performed (DSC, TGA, solubility, intrinsic dissolution rate, particle size, true density) and mechanical properties measured using a compaction replicator. Strained molecular conformations produce weaker compacts as they have higher energy than preferred conformations that usually lie close to energy minimums and oppose plastic deformation. It was observed that slip planes, which correspond to regions of weakest interaction between the planes, were associated with improved plasticity and stronger compacts. Apart from hydrogen bonds, profuse van der Waals forces can result in ineffective slip planes. Salts displaying two-dimensional densely hydrogen bonded layers produced stronger compacts than salts showing one-dimensional networks of non-bonded columns, probably by reducing the attachment energy between layers. When hydrogen bonds are created intramolecularly, it is possible that the mechanical properties are compromised as they do not contribute so much to create twodimensional densely bonded layers and they can force molecules into strained conformations. Some types of hydrogen bonding network may be associated with improved mechanical properties, such as type II, or R (10) 3 4 using graph-set notation, versus type III, or R (12) 4 8 , columns. This work clearly demonstrates the potential of investigating crystal structure-mechanical property relationship in pharmaceutical materials.

T-cell epitope prediction and immune complex simulation using molecular dynamics:state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

Evidence that interaction between conserved residues in transmembrane helices 2, 3, and 7 are crucial for human VPAC1 receptor activation

The VPAC(1) receptor belongs to family B of G protein-coupled receptors (GPCR-B) and is activated upon binding of the vasoactive intestinal peptide (VIP). Despite the recent determination of the structure of the N terminus of several members of this receptor family, little is known about the structure of the transmembrane (TM) region and about the molecular mechanisms leading to activation. In the present study, we designed a new structural model of the TM domain and combined it with experimental mutagenesis experiments to investigate the interaction network that governs ligand binding and receptor activation. Our results suggest that this network involves the cluster of residues Arg(188) in TM2, Gln(380) in TM7, and Asn(229) in TM3. This cluster is expected to be altered upon VIP binding, because Arg(188) has been shown previously to interact with Asp(3) of VIP. Several point mutations at positions 188, 229, and 380 were experimentally characterized and were shown to severely affect VIP binding and/or VIP-mediated cAMP production. Double mutants built from reciprocal residue exchanges exhibit strong cooperative or anticooperative effects, thereby indicating the spatial proximity of residues Arg(188), Gln(380), and Asn(229). Because these residues are highly conserved in the GPCR-B family, they can moreover be expected to have a general role in mediating function.

Toward an atomistic understanding of the immune synapse:large-scale molecular dynamics simulation of a membrane-embedded TCR-pMHC-CD4 complex

T-cell activation requires interaction of T-cell receptors (TCR) with peptide epitopes bound by major histocompatibility complex (MHC) proteins. This interaction occurs at a special cell-cell junction known as the immune or immunological synapse. Fluorescence microscopy has shown that the interplay among one agonist peptide-MHC (pMHC), one TCR and one CD4 provides the minimum complexity needed to trigger transient calcium signalling. We describe a computational approach to the study of the immune synapse. Using molecular dynamics simulation, we report here on a study of the smallest viable model, a TCR-pMHC-CD4 complex in a membrane environment. The computed structural and thermodynamic properties are in fair agreement with experiment. A number of biomolecules participate in the formation of the immunological synapse. Multi-scale molecular dynamics simulations may be the best opportunity we have to reach a full understanding of this remarkable supra-macromolecular event at a cell-cell junction.

Molecular dynamics study of solvent transport in nanoscale osmosis

An ideal of osmotic equilibrium between an ideal solution and pure solvent separated by a semi-permeable membrane is studied numerically using the method of molecular dynamics. The osmotic flow is observed as the inflow of the solvent across the membrane from the dilute to the concentrated side. The validity of van't Hoff's law for osmotic pressure is confirmed over a wide range of concentrations. It is found that the law is established by a balance between non-uniform partial pressures of solute and solvent. Furthermore, the present model permits an understanding of the mechanism of the osmotic flow in the relaxation process as the liquids evolve from the initial state to the equilibrium state. We focus in particular on the interaction between solute and solvent. ©2008 The Physical Society of Japan.

Mapping interaction sites within the N-terminus of the calcitonin gene-related peptide receptor; the role of residues 23-60 of the calcitonin receptor-like receptor

The calcitonin receptor-like receptor (CLR) acts as a receptor for the calcitonin gene-related peptide (CGRP) but in order to recognize CGRP, it must form a complex with an accessory protein, receptor activity modifying protein 1 (RAMP1). Identifying the protein/protein and protein/ligand interfaces in this unusual complex would aid drug design. The role of the extreme N-terminus of CLR (Glu23-Ala60) was examined by an alanine scan and the results were interpreted with the help of a molecular model. The potency of CGRP at stimulating cAMP production was reduced at Leu41Ala, Gln45Ala, Cys48Ala and Tyr49Ala; furthermore, CGRP-induced receptor internalization at all of these receptors was also impaired. Ile32Ala, Gly35Ala and Thr37Ala all increased CGRP potency. CGRP specific binding was abolished at Leu41Ala, Ala44Leu, Cys48Ala and Tyr49Ala. There was significant impairment of cell surface expression of Gln45Ala, Cys48Ala and Tyr49Ala. Cys48 takes part in a highly conserved disulfide bond and is probably needed for correct folding of CLR. The model suggests that Gln45 and Tyr49 mediate their effects by interacting with RAMP1 whereas Leu41 and Ala44 are likely to be involved in binding CGRP. Ile32, Gly35 and Thr37 form a separate cluster of residues which modulate CGRP binding. The results from this study may be applicable to other family B GPCRs which can associate with RAMPs.

Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

Encapsulated membrane proteins:a simplified system for molecular simulation

Over the past 50 years there has been considerable progress in our understanding of biomolecular interactions at an atomic level. This in turn has allowed molecular simulation methods employing full atomistic modeling at ever larger scales to develop. However, some challenging areas still remain where there is either a lack of atomic resolution structures or where the simulation system is inherently complex. An area where both challenges are present is that of membranes containing membrane proteins. In this review we analyse a new practical approach to membrane protein study that offers a potential new route to high resolution structures and the possibility to simplify simulations. These new approaches collectively recognise that preservation of the interaction between the membrane protein and the lipid bilayer is often essential to maintain structure and function. The new methods preserve these interactions by producing nano-scale disc shaped particles that include bilayer and the chosen protein. Currently two approaches lead in this area: the MSP system that relies on peptides to stabilise the discs, and SMALPs where an amphipathic styrene maleic acid copolymer is used. Both methods greatly enable protein production and hence have the potential to accelerate atomic resolution structure determination as well as providing a simplified format for simulations of membrane protein dynamics.