Partial agonism, neutral antagonism, and inverse agonism at the human wild-type and constitutively active cholecystokinin-2 receptors

Cholecystokinin receptor-2 (CCK2R) is a G protein receptor that regulates a number of physiological functions. Activation of CCK2R and/or expression of a constitutively active CCK2R variant may contribute to human diseases, including digestive cancers. Search for antagonists of the CCK2R has been an important challenge during the last few years, leading to discovery of a set of chemically distinct compounds. However, several early-discovered antagonists turned out to be partial agonists. In this context, we carried out pharmacological characterization of six CCK2R antagonists using COS-7 cells expressing the human CCK2R or a CCK2R mutant having a robust constitutive activity on inositol phosphates production, and we investigated the molecular mechanisms which, at a CCK2R binding site, account for these features. Results indicated that three compounds, 3R(+)-N-(2,3-dihydro-1-methyl-2-oxo-5-phenyl-1H-1,4-benzodiazepin-3- yl)-N'-(3-methylphenyl)urea (L365,260), 4-{[2-[[3-(lH-indol-3-yl)-2- methyl-1-oxo-2-[[[1.7.7-trimethyl-bicyclo[2.2.1]hept-2-yl)-oxy]carbonyl]amino] propyl]amino]-1-phenylethyl]amino-4-oxo-[lS-la.2[S*(S*)]4a]} -butanoate N-methyl-D-glucamine (PD135, 158), and (R)-1-naphthalenepropanoic acid, b-[2-[[2-(8-azaspiro-[4.5]dec-8-ylcarbonyl)-4,6-dimethylphenyl]amino]-2- oxoethyl] (CR2945), were partial agonists; one molecule, 1-[(R)-2,3-dihydro-1- (2,3-dihydro-1-(2-methylphenacyl)-2-oxo-5-phenyl-1H-1,4-benzodiazepin-3-yl] -3-(3-methylphenyl)urea (YM022), was a neutral antagonist; and two compounds, N-(+)-[1-(adamant-1-ylmethyl)-2,4-dioxo-5-phenyl2,3,4,5-tetrahydro-1H-1, 5-benzodiazepin-3-yl]-N'-phenylurea (GV150,013X) and ([(N-[methoxy-3 phenyl] N-[N-methyl N-phenyl carbamoylmethyl], carbomoylmethyl)-3 ureido]-3-phenyl)2-propionic acid (RPR101,048), were inverse agonists. Furthermore, target- and pharmacophore-based docking of ligands followed by molecular dynamic simulation experiments resulted in consistent motion of aromatic residues belonging to a network presumably important for activation, thus providing the first structural explanations for the different pharmacological profiles of tested compounds. This study confirms that several referenced so-called antagonists are in fact partial agonists, and because of this undesired activity, we suggest that newly generated molecules should be preferred to efficiently block CCK2R-related physiological effects. Furthermore, data on the structural basis for the different pharmacological features of CCK2R ligands will serve to further clarify CCK2R mechanism of activation. Copyright © 2006 The American Society for Pharmacology and Experimental Therapeutics.

Solid and liquid charge-transfer complex formation between 1-methylnaphthalene and 1-alkyl-cyanopyridinium bis{(trifluoromethyl)sulfonyl}imide ionic liquids

Liquid charge-transfer (CT) complexes were observed to form on contacting electron-rich aromatics with electron withdrawing group appended 1-alkyl-4-cyanopyridinium ionic liquids (ILs). Cooling below the melting point of the ionic liquid resulted in crystallisation of ionic liquid from the complex for 2-cyano and 3-cyano pyridinium isomers and in the formation of a 1 : 1 IL : aromatic crystalline CT-complex with the 4-cyanopyridinium isomer. The liquid structure of a 1 : 1 mixture of 1-methyl-4-cyanopyridinium bis{(trifluoromethyl)sulfonyl} imide with 1-methylnaphthalene has been probed by neutron diffraction experiments and molecular dynamics simulations. A high degree of correlation between the experimental data and the simulations was found with a significant displacement of the anions from around the cation by the aromatic species and the resulting structure having pi-pi stacks between the cations and the aromatic.

Liquid Structure and Dynamics of Aqueous Isopropanol over gamma-Alumina

The liquid structures of thin films of aqueous solutions of 0, 7, 19, 50, and 100 mol % isopropanol above O/Al-terminated gamma-alumina surfaces have been investigated by means of classical molecular dynamics simulations. The structuring effect of the oxide oil the liquid mixtures is strong and heavily dependent on the local structure of the oxide. Two distinct re-ions are found oil the oxide Surface characterized by the degree of coordination of Al atoms. Above octahedral Al atoms, water and isopropanol molecules adsorb via the oxygen atoms to maximize the electrostatic interaction, whereas above tetrahedral Al sites the solvent molecules adsorb via hydrogen atoms with the oxygen atoms away front the surface. More mobility is found in the second layer compared with the first; however, its structure is still influenced significantly by the orientation of molecules in the first adsorbed layer. Qualitatively, the displacement of water from the surface by the adsorption of isopropanol occurs with 2.6 Water molecules lost for every alcohol molecule present based on the effective surface areas of the two species calculated from the pure simulations.

A crystallographic and modelling study of a human telomeric RNA (TERRA) quadruplex

DNA telomeric repeats in mammalian cells are transcribed to guanine-rich RNA sequences, which adopt parallel-stranded G-quadruplexes with a propeller-like fold. The successful crystallization and structure analysis of a bimolecular human telomeric RNA G-quadruplex, folded into the same crystalline environment as an equivalent DNA oligonucleotide sequence, is reported here. The structural basis of the increased stability of RNA telomeric quadruplexes over DNA ones and their preference for parallel topologies is described here. Our findings suggest that the 2'-OH hydroxyl groups in the RNA quadruplex play a significant role in redefining hydration structure in the grooves and the hydrogen bonding networks. The preference for specific nucleotides to populate the C3'-endo sugar pucker domain is accommodated by alterations in the phosphate backbone, which leads to greater stability through enhanced hydrogen bonding networks. Molecular dynamics simulations on the DNA and RNA quadruplexes are consistent with these findings. The computations, based on the native crystal structure, provide an explanation for RNA G-quadruplex ligand binding selectivity for a group of naphthalene diimide ligands as compared to the DNA G-quadruplex.

Computational Verification of Two Universal Relations for Simple Ionic Liquids. Kinetic Properties of a Model 2:1 Molten Salt

Two semianalytical relations [Nature, 1996, 381, 137 and Phys. Rev. Lett. 2001, 87, 245901] predicting dynamical coefficients of simple liquids on the basis of structural properties have been tested by extensive molecular dynamics simulations for an idealized 2:1 model molten salt. In agreement with previous simulation studies, our results support the validity of the relation expressing the self-diffusion coefficient as a Function of the radial distribution functions for all thermodynamic conditions such that the system is in the ionic (ie., fully dissociated) liquid state. Deviations are apparent for high-density samples in the amorphous state and in the low-density, low-temperature range, when ions condense into AB(2) molecules. A similar relation predicting the ionic conductivity is only partially validated by our data. The simulation results, covering 210 distinct thermodynamic states, represent an extended database to tune and validate semianalytical theories of dynamical properties and provide a baseline for the interpretation of properties of more complex systems such as the room-temperature ionic liquids.

Effect of hydrophobic nanopatches within an ionic surface on the structure of liquids

The structures of liquid water and isopropanol have been studied as a function of the size of a hydrophobic patch present in a model hydrophilic surface via molecular dynamics simulations. A significant anisotropy extending into the first few solvent layers is found over the patch which suggests implications for many real-world systems in which nanoscale heterogeneity is found.

Nonadiabatic forces in ion-solid interactions:The initial stages of radiation damage

The Born-Oppenheimer approximation is the keystone for molecular dynamics simulations of radiation damage processes; however, actual materials response involves nonadiabatic energy exchange between nuclei and electrons. In this work, time dependent density functional theory is used to calculate the electronic excitations produced by energetic protons in Al. We study the influence of these electronic excitations on the interatomic forces and find that they differ substantially from the adiabatic case, revealing a nontrivial connection between electronic and nuclear stopping that is absent in the adiabatic case. These results unveil new effects in the early stages of radiation damage cascades.

On thin ice: surface order and disorder during pre-melting

The effect of temperature on the structure of the ice Ih (0001) surface is considered through a series of molecular dynamics simulations on an ice slab. At relatively low temperatures (200K) a small fraction of surface self-interstitials (i.e. admolecules) appear that are formed exclusively from molecules leaving the outermost bilayer. At higher temperatures (ca. 250 K), vacancies start to appear in the inner part of the outermost bilayer exposing the underlying bilayer and providing sites with a high concentration of the dangling hydrogen bonds. Around 250-260 K aggregates of molecules formed on top of the outermost bilayer from self-interstitials become more mobile and have diffusivities approaching that of liquid water. At similar to 270-280 K the inner bilayer of one surface noticeably destructures and it appears that at above 285 K both surfaces are melting. The observed disparity in the onset of melting between the two sides of the slab is rationalised by considering the relationship between surface energy and the spatial distribution of protons at the surface; thermodynamic stability is conferred on the surface by maximising separations between dangling protons at the crystal exterior. Local hotspots associated with a high dangling proton density are suggested to be susceptible to pre-melting and may be more efficient at trapping species at the external surface than regions with low concentrations of protons thus potentially helping ice particles to catalyse reactions. A preliminary conclusion of this work is that only about 10-20 K below the melting temperature of the particular water potential employed is major disruption of the crystalline lattice noted which could be interpreted as being "liquid", the thickness of this film being about a nanometre.

The Origin of Phragmoplast Asymmetry

The phragmoplast coordinates cytokinesis in plants [1]. It directs vesicles to the midzone, the site where they coalesce to form the new cell plate. Failure in phragmoplast function results in aborted or incomplete cytokinesis leading to embryo lethality, morphological defects, or multinucleate cells [2, 3]. The asymmetry of vesicular traffic is regulated by microtubules [1, 4, 5, 6], and the current model suggests that this asymmetry is established and maintained through treadmilling of parallel microtubules. However, we have analyzed the behavior of microtubules in the phragmoplast using live-cell imaging coupled with mathematical modeling and dynamic simulations and report that microtubules initiate randomly in the phragmoplast and that the majority exhibit dynamic instability with higher turnover rates nearer to the midzone. The directional transport of vesicles is possible because the majority of the microtubules polymerize toward the midzone. Here, we propose the first inclusive model where microtubule dynamics and phragmoplast asymmetry are consistent with the localization and activity of proteins known to regulate microtubule assembly and disassembly.

The Role of Arg228 in the phosphorylation of galactokinase: The mechanism of GHMP kinases by QM/MM studies

GHMP kinases are a group of structurally-related small molecule kinases. They have been found in all kingdoms of life and are mostly responsible for catalysing the ATP-dependent phosphorylation of intermediary metabolites. Although the GHMP kinases are of clinical, pharmaceutical and biotechnological importance, the mechanism of GHMP-kinases is controversial. A catalytic base mechanism was suggested for mevalonate kinase that has a structural feature of the ?-phosphate of ATP close to an aspartate residue; however, for one GHMP member, homoserine kinase, where the residue acting as general base is absent, a direct phosphorylation mechanism was suggested. Furthermore, it has been proposed by some authors that all the GHMP kinases function via the direct phosphorylation mechanism. This controversy in mechanism has limited our ability to exploit these enzymes as drug targets and in biotechnology. Here the phosphorylation reaction mechanism of the human galactokinase, a member of GHMP kinase was investigated using molecular dynamics simulations and density functional theory-based QM/MM calculations (B3LYP-D/AMBER99). The reaction coordinates were localized by potential energy scan using adiabatic mapping method. Our results indicate that a highly conserved Glu174 captures Arg105 to the proximity of the a-phosphate of ATP forming a H-bond network, therefore the mobility of ATP in the large oxyanion hole is restricted. Arg228 functions to stabilize the negative charge developed at the ß,?-bridging oxygen of the ATP during bond cleavage. The reaction occurs via direct phosphorylation mechanism and the Asp186 in proximity of ATP does not directly participate in the reaction pathway. Since Arg228 is not conserved among GHMP kinases, reagents which form interactions with Arg228, and therefore can interrupt its function in phosphorylation may be developed into potential selective inhibitors for galactokinase.

Comparison of dynamics of wildtype and V94M human UDP-galactose 4-epimerase-A computational perspective on severe epimerase-deficiency galactosemia

UDP-galactose 4'-epimerase (GALE) catalyzes the interconversion of UDP-galactose and UDP-glucose, an important step in galactose catabolism. Type III galactosemia, an inherited metabolic disease, is associated with mutations in human GALE. The V94M mutation has been associated with a very severe form of type III galactosemia. While a variety of structural and biochemical studies have been reported that elucidate differences between the wildtype and this mutant form of human GALE, little is known about the dynamics of the protein and how mutations influence structure and function. We performed molecular dynamics simulations on the wildtype and V94M enzyme in different states of substrate and cofactor binding. In the mutant, the average distance between the substrate and both a key catalytic residue (Tyr157) and the enzyme-bound NAD(+) cofactor and the active site dynamics are altered making substrate binding slightly less stable. However, overall stability or dynamics of the protein is not altered. This is consistent with experimental findings that the impact is largely on the turnover number (kcat), with less substantial effects on Km. Active site fluctuations were found to be correlated in enzyme with substrate bound to just one of the subunits in the homodimer suggesting inter-subunit communication. Greater active site loop mobility in human GALE compared to the equivalent loop in Escherichia coli GALE explains why the former can catalyze the interconversion of UDP-N-acetylgalactosamine and UDP-N-acetylglucosamine while the bacterial enzyme cannot. This work illuminates molecular mechanisms of disease and may inform the design of small molecule therapies for type III galactosemia.

Impact of Heat Pump Load on Distribution Networks

Heat pumps can provide domestic heating at a cost that is competitive with oil heating in particular. If the electricity supply contains a significant amount of renewable generation, a move from fossil fuel heating to heat pumps can reduce greenhouse gas emissions. The inherent thermal storage of heat pump installations can also provide the electricity supplier with valuable flexibility. The increase in heat pump installations in the UK and Europe in the last few years poses a challenge for low-voltage networks, due to the use of induction motors to drive the pump compressors. The induction motor load tends to depress voltage, especially on starting. The paper includes experimental results, dynamic load modelling, comparison of experimental results and simulation results for various levels of heat pump deployment. The simulations are based on a generic test network designed to capture the main characteristics of UK distribution system practice. The simulations employ DIgSlILENT to facilitate dynamic simulations that focus on starting current, voltage variations, active power, reactive power and switching transients.

PLUMED 2: New feathers for an old bird

Enhancing sampling and analyzing simulations are central issues in molecular simulation. Recently, we introduced PLUMED, an open-source plug-in that provides some of the most popular molecular dynamics (MD) codes with implementations of a variety of different enhanced sampling algorithms and collective variables (CVs). The rapid changes in this field, in particular new directions in enhanced sampling and dimensionality reduction together with new hardware, require a code that is more flexible and more efficient. We therefore present PLUMED 2 here a,complete rewrite of the code in an object-oriented programming language (C++). This new version introduces greater flexibility and greater modularity, which both extends its core capabilities and makes it far easier to add new methods and CVs. It also has a simpler interface with the MD engines and provides a single software library containing both tools and core facilities. Ultimately, the new code better serves the ever-growing community of users and contributors in coping with the new challenges arising in the field.

Program summary

Program title: PLUMED 2

Catalogue identifier: AEEE_v2_0

Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEEE_v2_0.html

Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland

Licensing provisions: Yes

No. of lines in distributed program, including test data, etc.: 700646

No. of bytes in distributed program, including test data, etc.: 6618136

Distribution format: tar.gz

Programming language: ANSI-C++.

Computer: Any computer capable of running an executable produced by a C++ compiler.

Operating system: Linux operating system, Unix OSs.

Has the code been vectorized or parallelized?: Yes, parallelized using MPI.

RAM: Depends on the number of atoms, the method chosen and the collective variables used.

Classification: 3, 7.7, 23. Catalogue identifier of previous version: AEEE_v1_0.

Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 1961.

External routines: GNU libmatheval, Lapack, Bias, MPI. (C) 2013 Elsevier B.V. All rights reserved.

Probing the Unfolded Configurations of a β-Hairpin Using Sketch-Map

This work examines the conformational ensemble involved in β-hairpin folding by means of advanced molecular dynamics simulations and dimensionality reduction. A fully atomistic description of the protein and the surrounding solvent molecules is used, and this complex energy landscape is sampled by means of parallel tempering metadynamics simulations. The ensemble of configurations explored is analyzed using the recently proposed sketch-map algorithm. Further simulations allow us to probe how mutations affect the structures adopted by this protein. We find that many of the configurations adopted by a mutant are the same as those adopted by the wild-type protein. Furthermore, certain mutations destabilize secondary-structure-containing configurations by preventing the formation of hydrogen bonds or by promoting the formation of new intramolecular contacts. Our analysis demonstrates that machine-learning techniques can be used to study the energy landscapes of complex molecules and that the visualizations that are generated in this way provide a natural basis for examining how the stabilities of particular configurations of the molecule are affected by factors such as temperature or structural mutations.

General, PDB-based collective variables for protein folding

New, automated forms of data-analysis are required in order to understand the high-dimensional trajectories that are obtained from molecular dynamics simulations on proteins. Dimensionality reduction algorithms are particularly appealing in this regard as they allow one to construct unbiased, low-dimensional representations of the trajectory using only the information encoded in the trajectory. The downside of this approach is that different sets of coordinates are required for each different chemical systems under study precisely because the coordinates are constructed using information from the trajectory. In this paper we show how one can resolve this problem by using the sketch-map algorithm that we recently proposed to construct a low-dimensional representation of the structures contained in the protein data bank (PDB). We show that the resulting coordinates are as useful for analysing trajectory data as coordinates constructed using landmark configurations taken from the trajectory and that these coordinates can thus be used for understanding protein folding across a range of systems.