The role of charge interactions in muscarinic agonist binding, and receptor-response coupling

Site-directed mutagenesis has been used to evaluate the roles of the key aspartate and arginine residues in transmembrane domain three of the muscarinic receptors. The results suggest that the formation of an ionic bond between the Asp carboxylate group and the onium headgroup is essential to anchor acetylcholine in its active, bound conformation in both binary agonist-receptor and ternary agonist-receptor-G-protein complexes, but that secondary, non-productive binding modes, promoted by non-polar forces, may contribute to binary complex formation by other ligands. The positive charge of the arginyl side-chain is central to the recognition, and subsequent activation of G-proteins by the agonist-M1 mAChR complex.     

Structure and dynamics of primary hydration shell of phosphatidylcholine bilayers at subzero temperatures.

Deuterium NMR relaxation and intensity measurements of the 2H-labeled H2O/dimyristoyl phosphatidylcholine bilayer were performed to understand the molecular origin of the freezing event of phospholipid headgroup and the structure and dynamics of unfrozen water molecules in the interbilayer space at subzero temperatures. The results suggest that about one to two water molecules associated with the phosphate group freeze during the freezing event of phospholipid headgroups, whereas about five to six waters near the trimethylammonium group behave as a water cluster and remain unfrozen at temperatures as low as -70 degrees C. In addition, temperature-dependent T1 and T2 relaxation times suggest that dynamic coupling occurs not only between the phosphate group and its bound water, but also between the methyl group and the adjacent water molecules. Based on these observations, the primary hydration shell of phosphatidylcholine headgroup at subzero temperatures is suggested to consist of two distinct regions: a clathrate-like water cluster, most likely a water pentamer, near the hydrophobic methyl group, and hydration water molecules associated with the phosphate group.     

Cholesterol effects on the phosphatidylcholine bilayer polar region: a molecular simulation study

A molecular dynamics (MD) simulation of a fully hydrated, liquid-crystalline dimyristoylphosphatidylcholine (DMPC)-Chol bilayer membrane containing approximately 22 mol% Chol was carried out for 4.3 ns. The bilayer reached thermal equilibrium after 2.3 ns of MD simulation. A 2.0-ns trajectory generated during 2.3-4.3 ns of MD simulation was used for analyses to determine the effects of Chol on the membrane/water interfacial region. In this region, 70% of Chol molecules are linked to DMPC molecules via short-distance interactions, where the Chol hydroxyl group (OH-Chol) is 1) charge paired to methyl groups of the DMPC choline moiety ( approximately 34%), via the hydroxyl oxygen atom (Och); 2) water bridged to carbonyl ( approximately 19%) and nonester phosphate ( approximately 14%) oxygen atoms, via both Och and the hydroxyl hydrogen atom (Hch); and 3) directly hydrogen (H) bonded to carbonyl ( approximately 11%) and nonester phosphate ( approximately 5%) oxygen atoms, via Hch ( approximately 17% of DMPC-Chol links are multiple). DMPC's gamma-chain carbonyl oxygen atom is involved in 44% of water bridges and 51% of direct H bonds formed between DMPC and Chol. On average, a Chol molecule forms 0.9 links with DMPC molecules, while a DMPC molecule forms 2.2 and 0.3 links with DMPC and Chol molecules, respectively. OH-Chol makes hydrogen bonds with 1.1 water molecules, preferentially via Hch. The average number of water molecules H bonded to the DMPC headgroup is increased by 7% in the presence of Chol. These results indicate that inclusion of Chol decreases interlipid links and increases hydration in the polar region of the membrane.     

Why do water molecules around small hydrophobic solutes form stronger hydrogen bonds than in the bulk?

Molecular solutes are known to have a strong effect on the structural and dynamical properties of the surrounding water. In our recent study (PNAS, 114, 322 (2017)) we have identified the presence of strengthened water hydrogen bonds near hydrophobic solutes by using both IR spectroscopy and ab-initio molecular dynamics simulations. The water molecules involved in the enhanced hydrogen bonding have been shown to display extensive structural ordering and restricted mobility. We observed that an individual pair of water molecules can make stronger hydrogen bond to each other if it is not surrounded by intercalating water molecules. Here we present compelling simulation results which unravel a simple mechanistic picture of the emergence of the hydrogen bond (HB) strengthening around solvated methane. We show explicitly that actual absence of water molecules within the excluded volume due to the hydrophobic molecule assures smaller residual torque on neighboring water molecules enabling the formation of stronger HBs between them.     

Dielectric study of the cooperative order–disorder transition in aqueous solutions of schizophyllan,a triple-helical polysaccharide

Schizophyllan exists in aqueous solution as a triple helix, which is intact at room temperature. Its aqueous solution forms some ordered structure at low temperatures but undergoes a sharp transition to a disordered structure as the temperature is raised. The transition temperature T_c is about 7 and 18°C for H₂O and D₂O solutions, respectively. This transition was followed by time-domain reflectometry to investigate dynamic aspects of the transition. In addition to a major peak around 10 GHz, the dielectric dispersion curve of a 20 wt % schizophyllan in D₂O exhibited a small peak around 100 MHz below T_c and around 10 MHz above T_c. The major peak is due to bulk water, whereas the 100 MH_z peak is assigned to “bound” or “structured” water, and that around 10 MHz to side-chain glucose residues. However, unlike usual bound water reported for biopolymer solutions, this “structured” water disappears abruptly when the temperature becomes close to T_c without accompanying a conformational transition of the main chain. The above assignment is consistent with the structure of the ordered phase derived from previous static data that it consists of side-chain glucose residues along with nearby water molecules surrounding the helix core that are interacting with each other loosely through hydrogen bonds, and spreads radially only a layer of one or two water molecules but a long distance along the helix axis. © 1995 John Wiley & Sons, Inc.     

On the properties of aqueous amide solutions through classical molecular dynamics simulations

The structure and dynamics of infinitely diluted aqueous amide solutions is studied for 13 compounds in the NVT ensemble using classical molecular dynamics simulations. The aim of this work is to provide valuable insights into the effect of amides on liquid water properties in order to understand the amides role in the kinetic inhibition of clathrate hydrate formation in natural gas mixtures. The OPLS-AA forcefield is used to describe the amides, with parameters obtained through fitting of computed B3LYP/6-311++g^{* *} data when not available in the literature, and the SPC-E model is applied for water molecules. Structural properties of the solutions are analyzed via calculated radial distribution functions and dynamic properties are studied with the computed mean square displacements and velocity autocorrelation functions. Most of the studied compounds show a remarkable structuring effect on the surrounding water with strong interactions resulting from hydrogen bonding between solute and solvent molecules. Hydrophobic and hydrophilic synergistic effects influence the amide–water interaction and the properties of the water solvation shells around amides.     

Ab initio molecular dynamics study of the hydration of the formohydroxamate anion

We apply ab initio molecular dynamics (AIMD) to study the hydration structures and electronic properties of the formohydroxamate anion in liquid water. We consider the cis- nitrogen-deprotonated, cis- oxygen-deprotonated, and trans- oxygen-deprotonated formohydroxamate tautomers. They form an average of 6.3, 6.9, and 6.0 hydrogen bonds with water molecules, respectively. The predicted pair correlation functions and time dependence of the hydration numbers suggest that water is highly structured around the nominally negatively charged oxime oxygen in O-deprotonated tautomers but significantly less so around the nitrogen atom in the N-deprotonated species. Wannier function analysis suggests that, in the O-deprotonated anions, the negative charge is concentrated on the oxime oxygen, while in the N-deprotonated case, it is partially delocalized between the nitrogen and the adjoining oxime oxygen atom.     

Bilayer interfacial properties modulate the binding of amphipathic peptides

The free energy of transfer (DeltaG degrees ) from water to lipid bilayers was measured for two amphipathic peptides, the presequence of the mitochondrial peptide rhodanese (MPR) and melittin. Experiments were designed to determine the effects on peptide partitioning of the addition of lipids that produce structural modifications to the bilayer/water interface. In particular, the addition of cholesterol or the cholesterol analog 6-ketocholestanol increases the bilayer area compressibility modulus, indicating that these molecules modify lipid-lipid interactions in the plane of the bilayer. The addition of 6-ketocholestanol or lipids with attached polyethylene glycol chains (PEG-lipids) modify the effective thickness of the interfacial region; 6-ketocholestanol increases the width of hydrophilic headgroup region in the direction of the acyl chains whereas the protruding PEG chains of PEG-lipids increase the structural width of the headgroup region into the surrounding aqueous phase. The incorporation of PEG-lipids with PEG molecular weights of 2000 or 5000 had no appreciable effect on peptide partitioning that could not be accounted for by the presence of surface charge. However, for both MPR and melittin DeltaG degrees decreased linearly with increasing bilayer compressibility modulus, demonstrating the importance of bilayer mechanical properties in the binding of amphipathic peptides.     

Molecular mechanics and dynamics simulations of various dispersant models on the water surface (001)

Five dispersant-molecule models of succinimide, acrylate, imide, phenylsulfonic and salicyl were used to study their interactions with the water surface (001). The interaction energy, molecular configuration, charge distribution and radial distribution function (RDF) curve for each of the dispersant molecules were analyzed from the molecular mechanics (MM) and molecular dynamics (MD) simulation results. It can be seen that the system energies, mostly electrostatic and hydrogen bond energies, were reduced significantly when the dispersant molecules interacted with the water surface. The hydrophilic group of a dispersant molecule can attach itself to the water surface firmly and reach for a stable energy-minimized configuration, which is helpful to the dispersants' dispersancy. The influence exerted by the hydrophobic group of the dispersant molecule, which was the substituted hydrocarbon chain of n-octadecanyl in this paper, is discussed in comparison with the naked polar headgroup. Steric configuration, charge distribution and substitute hydrocarbon chain of the dispersant molecule influenced the interaction between dispersants and polar water surface.     

What makes the bioactive lipids phosphatidic acid and lysophosphatidic acid so special?

Phosphatidic acid and lysophosphatidic acid are minor but important anionic bioactive lipids involved in a number of key cellular processes, yet these molecules have a simple phosphate headgroup. To find out what is so special about these lipids, we determined the ionization behavior of phosphatidic acid (PA) and lysophosphatidic acid (LPA) in extended (flat) mixed lipid bilayers using magic angle spinning 31P NMR. Our data show two surprising results. First, despite identical phosphomonoester headgroups, LPA carries more negative charge than PA when present in a phosphatidylcholine bilayer. Dehydroxy-LPA [1-oleoyl-3-(phosphoryl)propanediol] behaves in a manner identical to that of PA, indicating that the difference in negative charge between LPA and PA is caused by the hydroxyl on the glycerol backbone of LPA and its interaction with the phosphomonoester headgroup. Second, deprotonation of phosphatidic acid and lysophosphatidic acid was found to be strongly stimulated by the inclusion of phosphatidylethanolamine in the bilayer, indicating that lipid headgroup charge depends on local lipid composition and will vary between the different subcellular locations of (L)PA. Our findings can be understood in terms of a hydrogen bond formed within the phosphomonoester headgroup of (L)PA and its destabilization by competing intra- or intermolecular hydrogen bonds. We propose that this hydrogen bonding property of (L)PA is involved in the various cellular functions of these lipids.     

17O NMR investigation of the hydration of L-alanine and L-proline in water/Me2SO mixtures

The hydration state of L-Alanine and L-Proline has been assessed via 17O NMR. At neutral and basic pH, two water molecules are hydrogen bonded at the carboxylate group, one to each oxygen, whereas a third water molecule is hydrogen bonded to the protonated COOH group at acidic pH, via the hydroxyl hydrogen. The possible formation of dimers and/or higher complexes in DMSO is indicated not only from the chemical shift but also from the linewidth of the amino acids.     

Computer simulation of hydration networks around amino acids

The networks of solvent hydrogen bonds around polar and apolar amino acids have been studied by computer simulation techniques using a non-pair additive model for the water molecules interactions. Analysis of the simulated aqueous solutions has shown the presence of water molecules which (a) form a bridge around individual polar solute atoms (self-bridging loops) and (b) form chains between different polar solute atoms (polar bridging chains). Some of these networks associated with polar solute atoms from pentagons but 4, 6 and 7 sided polygons are also seen. The water molecule close to apolar solute atoms (<4.0 Å) also form irregular networks with polygons of 4, 5, 6 and 7 sides. These networks are compared with those found experimentally in ice, clathrates and crystal hydrates of macromolecules.     

Interaction among the twelve-residue segment of antifreeze protein type I,or its mutants,water and a hexagonal ice crystal

We have carried out a molecular dynamics analysis on a mixture of supercooled water, a hexagonal ice crystal and segments of winter flounder antifreeze protein. The segment consists of nine alanine residues, two threonine residues and one asparagine residue. Mutant segments, in which the threonine residues are replaced with valine residues, or serine residues, are also used. It is found that the threonine residue near the asparagine residue of the original segment is located in the vicinity of the prism face of the ice crystal. This is due to the hydrogen bond between the hydrophilic sites of these residues and water molecules, and the hydrogen bond between these water molecules and the water molecules on the ice surface. The valine and serine residues in the mutant segments do not approach the prism face of the ice crystal compared with the threonine residue near the asparagine residue. The motion of five segments, closely located side by side, is not remarkable. This is because of the gathering of water molecules caused by hydrophobic hydration, not only around alanine residues but also around the methyl sites of threonine residues.     

A Possible Evolutionary Origin for the Mn<Subscript>4</Subscript> Cluster in Photosystem II: From Manganese Superoxide Dismutase to Oxygen Evolving Complex

The recently published X-ray absorption fine structure of photosystem II provides a more detailed architecture of the oxygen-evolving complex (OEC) and the surrounding amino acids. In this paper, a comparison between manganese superoxide dismutase, dinuclear manganese catalase enzymes and the oxygen evolving complex in photosystem II is reported. The author suggests that the development of oxygenic photosynthesis occurred in steps, the first of which involved only one manganese ion (Mn(II)) that oxidized two water molecules to hydrogen peroxide and then oxygen.     

Molecular dynamics investigation of the structure of a fully hydrated gel-phase dipalmitoylphosphatidylcholine bilayer.

We report the results of a constant pressure and temperature molecular dynamics simulation of a gel-phase dipalmitoylphosphatidylcholine bilayer with nw = 11.8 water molecules/lipid at 19 degrees C. The results of the simulation were compared in detail with a variety of x-ray and neutron diffraction data. The average positions of specific carbon atoms along the bilayer normal and the interlamellar spacing and electron density profile were in very good agreement with neutron and x-ray diffraction results. The area per lipid and the details of the in-plane hydrocarbon chain structure were in excellent agreement with wide-angle x-ray diffraction results. The only significant deviation is that the chains met in a pleated arrangement at the bilayer center, although they should be parallel. Novel discoveries made in the present work include the observation of a bimodal headgroup orientational distribution. Furthermore, we found that there are a significant number of gauche conformations near the ends of the hydrocarbon chains and, in addition to verifying a previous suggestion that there is partial rotational ordering in the hydrocarbon chains, that the two chains in a given molecule are inequivalent with respect to rotations. Finally, we have investigated the lipid/water interface and found that the water penetrates beneath the headgroups, but not as far as the carbonyl groups, that the phosphates are strongly hydrated almost exclusively at the nonesterified oxygen atoms, and that the hydration of the ammonium groups is more diffuse, with some water molecules concentrated in the grooves between the methyl groups.     

Molecular simulation study of structural and dynamic properties of mixed DPPC/DPPE bilayers

Molecular dynamics simulations have been used to study structural and dynamic properties of fully hydrated mixed 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE) bilayers at 0, 25, 50, 75, and 100 mol % DPPE. Simulations were performed for 50 ns at 350 K and 1 bar for the liquid-crystalline state of the mixtures. Results show that the average area per headgroup reduces from 0.65 +/- 0.01 nm(2) in pure DPPC to 0.52 +/- 0.01 nm(2) in pure DPPE systems. The lipid tails become more ordered with increasing DPPE concentration, resulting in a slight increase in membrane thickness (3.43 +/- 0.01 nm in pure DPPC to 4.00 +/- 0.01 nm in pure DPPE). The calculated area per headgroup and order parameter for pure DPPE deviates significantly from available experimental measurements, suggesting that the force field employed requires further refinement. In-depth analysis of the hydrogen-bond distribution in DPPE molecules shows that the amine groups strongly interact with the phosphate and carbonyl groups through inter/intramolecular hydrogen bonds. This yields a bilayer structure with DPPE headgroups preferentially located near the lipid phosphate and ester oxygens. It is observed that increasing DPPE concentrations causes competitive hydrogen bonding between the amine groups (hydrogen-donor) and the phosphate/carbonyl groups or water (hydrogen-acceptor). Due to the increasing number of hydrogen-donors from DPPE molecules with increasing concentration, DPPE becomes more hydrated. Trajectory analysis shows that DPPE molecules in the lipid mixtures move laterally and randomly around the membrane surface and the movement becomes more localized with increasing DPPE concentrations. For the conditions and simulation time considered, no aggregation or phase separation was observed between DPPC and DPPE.     

A non-invasive method for measuring preimplantation embryo physiology

The physiology of the early embryo may be indicative of embryo vitality and therefore methods for non-invasively monitoring physiological parameters from embryos could improve preimplantation diagnoses. The self-referencing electrophysiological technique is capable of non-invasive measurement of the physiology of individual cells by monitoring the movement of ions and molecules between the cell and the surrounding media. Here we use this technique to monitor gradients of calcium, potassium, oxygen and hydrogen peroxide around individual mouse preimplantation embryos. The calcium-sensitive electrode in self-referencing mode identified a region of elevated calcium concentration (approximately 0.25 pmol) surrounding each embryo. The calcium gradient surrounding embryos was relatively steep, such that the region of elevated calcium extended into the medium only 4 microns from the embryo. By contrast, using an oxygen-sensitive electrode an extensive gradient of reduced dissolved oxygen concentration was measured surrounding the embryo and extended tens of micrometres into the medium. A gradient of neither potassium nor hydrogen peroxide was observed around unperturbed embryos. We also demonstrate that monitoring the physiology of embryos using the self-referencing technique does not compromise their subsequent development. Blastocyts studied with the self-referencing technique implanted and developed to term at the same frequency as did unexamined, control embryos. Therefore, the self-referencing electrode provides a valuable non-invasive technique for studying the physiology and pathophysiology of individual embryos without hindering their subsequent development.     

Oxygen-oxygen distances in protein-bound crystallographic water suggest the presence of protonated clusters

BackgroundThe availability of high-resolution X-ray structures has shown that proteins contain numerous water molecules, but their role is still not fully understood. Protonated and deprotonated water species are often involved in biochemical reactions. However protons are exceedingly difficult to detect directly because they are electron-poor species.MethodsThe oxygen‑oxygen distance of the crystallographic water molecules was analyzed in a large high-resolution data set. Moreover, a detailed analysis was carried out on the protein-bound water in the available structures of carbonic anhydrase II and cytochrome c oxidase, chosen as protein models in which protonated and deprotonated water species play a significant role.ResultsThe analysis shows an excess of water-water distances below the expected value for hydrogen bond. In the cavities and on the surface of the considered model proteins, clusters of water molecules are found, whose structure suggests the presence of chemical species deriving from self-ionization of water.ConclusionsThe presence of a small maximum below the hydrogen bond threshold in the oxygen‑oxygen distance distribution of crystallographic water molecules, along with the location of many of these water clusters, suggest the presence of Zundel-like structures in, or near, the proteins. Particularly significant is the presence of such structures in protein regions which have been identified as proton antennae or channels.General significanceThis work shows the possibilities, still unexplored, offered by this type of analysis in detecting in structures obtained by X-ray diffraction the presence of aqueous protons or hydroxide ions, which are chemical species as important as elusive.     

Hydration of cis and trans N-methylformamide as revealed by the use of 17O-NMR,molecular mechanics,and ab initio calculations

The solvation of cis and trans N-methylformamide (NMF) by water was investigated using a combination of ¹⁷O-nmr spectroscopy, classical molecular mechanics [MM2(77) and MM2(87)] force field, and ab initio 4-31G* gradient optimization calculations. In dilute aqueous solution, the ¹⁷O-nmr spectra of NMF indicate strong shielding by 66.9 and 66.1 ppm for the cis and trans amide oxygens, respectively, compared to those values obtained in dilute toluene solution. This demonstrates that both isomers are equally solvated by molecules of water, which are further hydrogen bonded to molecules of water of the bulk solvent. Molecular mechanics simulations were carried out for cis and trans NMF in a cluster of water molecules. Radial distribution functions show structural contacts by several water molecules at the amide CO and NH group, which are significantly more pronounced with MM2 (87) calculations. Ab initio 4-31G* gradient optimization calculations on the supermolecule trans NMF-(H₂O)₃ indicates the presence of more than two hydrogen-bond contacts at the carbonyl oxygen. This is in agreement with MM2 calculations and provides further evidence for multiple acceptor properties of the amide oxygen and an out of the amide plane arrangement of the bound molecules of water. Comparison of the integration data to the first radial distribution function (rdf) minima shows that the local solvation of the CO and NH groups is very similar for both cis and trans isomers. The intermolecular geometric parameters of the supermolecule trans NMA–(H₂O)₃ and the first rdf maxima resulting from MM2 (87) and MM2 (77) calculations are compared with distribution of water molecules around the CO and NH groups of peptides and proteins resulting from x-ray and neutron diffraction experiments. The rdfs involving the methyl group of NMF demonstrate the nonrandom distribution of solvent sites with first maxima in reasonable agreement with distribution of water molecules around the apolar side chain of amino acid residues in proteins. © 1995 John Wiley & Sons, Inc.