Accurate theoretical prediction of vibrational frequencies in an inhomogeneous dynamic environment: a case study of a glutamate molecule in water solution and in a protein-bound form

Ion Permeation Mechanism in Epithelial Calcium Channel TRVP6

Calcium is the most abundant metal in the human body that plays vital roles as a cellular electrolyte as well as the smallest and most frequently used signaling molecule. Calcium uptake in epithelial tissues is mediated by tetrameric calcium-selective transient receptor potential (TRP) channels TRPV6 that are implicated in a variety of human diseases, including numerous forms of cancer. We used TRPV6 crystal structures as templates for molecular dynamics simulations to identify ion binding sites and to study the permeation mechanism of calcium and other ions through TRPV6 channels. We found that at low Ca²⁺ concentrations, a single calcium ion binds at the selectivity filter narrow constriction formed by aspartates D541 and allows Na⁺ permeation. In the presence of ions, no water binds to or crosses the pore constriction. At high Ca²⁺ concentrations, calcium permeates the pore according to the knock-off mechanism that includes formation of a short-lived transition state with three calcium ions bound near D541. For Ba²⁺, the transition state lives longer and the knock-off permeation occurs slower. Gd³⁺ binds at D541 tightly, blocks the channel and prevents Na⁺ from permeating the pore. Our results provide structural foundations for understanding permeation and block in tetrameric calcium-selective ion channels.

Role of the Ion Channel Extracellular Collar in AMPA Receptor Gating

AMPA subtype ionotropic glutamate receptors mediate fast excitatory neurotransmission and are implicated in numerous neurological diseases. Ionic currents through AMPA receptor channels can be allosterically regulated via different sites on the receptor protein. We used site-directed mutagenesis and patch-clamp recordings to probe the ion channel extracellular collar, the binding region for noncompetitive allosteric inhibitors. We found position and substitution-dependent effects for introduced mutations at this region on AMPA receptor gating. The results of mutagenesis suggested that the transmembrane domains M1, M3 and M4, which contribute to the ion channel extracellular collar, undergo significant relative displacement during gating. We used molecular dynamics simulations to predict an AMPA receptor open state structure and rationalize the results of mutagenesis. We conclude that the ion channel extracellular collar plays a distinct role in gating and represents a hub for powerful allosteric modulation of AMPA receptor function that can be used for developing novel therapeutics.

An intrinsically disordered photosystem II subunit, PsbO, provides a structural template and a sensor of the hydrogen-bonding network in photosynthetic water oxidation

Photosystem II (PSII) is a membrane-bound enzyme that utilizes solar energy to catalyze the photooxidation of water. Molecular oxygen is evolved after four sequential light-driven oxidation reactions at the Mn4CaO5 oxygen-evolving complex, producing five sequentially oxidized states, Sn. PSII is composed of 17 membrane-spanning subunits and three extrinsic subunits, PsbP, PsbQ, and PsbO. PsbO is intrinsically disordered and plays a role in facilitation of the water oxidizing cycle. Native PsbO can be removed and substituted with recombinant PsbO, thereby restoring steady-state activity. In this report, we used reaction-induced Fourier transform infrared spectroscopy to obtain information concerning the role of PsbP, PsbQ, and PsbO during the S state cycle. Light-minus-dark difference spectra were acquired, monitoring structural changes associated with each accessible flash-induced S state transition in a highly purified plant PSII preparation (Triton X-100, octylthioglucoside). A comparison of S2 minus S1 spectra revealed that removal of PsbP and PsbQ had no significant effect on the data, whereas amide frequency and intensity changes were associated with PsbO removal. These data suggest that PsbO acts as an organizational template for the PSII reaction center. To identify any coupled conformational changes arising directly from PsbO, global (13)C-PsbO isotope editing was employed. The reaction-induced Fourier transform infrared spectra of accessible S states provide evidence that PsbO spectral contributions are temperature (263 and 277 K) and S state dependent. These experiments show that PsbO undergoes catalytically relevant structural dynamics, which are coupled over long distance to hydrogen-bonding changes at the Mn4CaO5 cluster.

A density functional and quantum Monte Carlo study of glutamic acid in vacuo and in a dielectric continuum medium

We present density functional theory (DFT) and quantum Monte Carlo (QMC) calculations of the glutamic acid and glutamate ion in vacuo and in various dielectric continuum media within the polarizable continuum model (PCM). In DFT, we employ the integral equation formalism variant of PCM while, in QMC, we use a PCM scheme we have developed to include both surface and volume polarization. We investigate the gas-phase protonation thermochemistry of the glutamic acid using a large set of structural conformations, and find that QMC is in excellent agreement with the best available theoretical and experimental results. For the solvated glutamic acid and glutamate ion, we perform DFT calculations for dielectric constants, ε, between 4 and 78. We find that the glutamate ion in the zwitterionic form is more stable than the non-zwitterionic form over the whole range of dielectric constants, while the glutamic acid is more stable in its non-zwitterionic form at ε = 4. The dielectric constant at which the two glutamic acid species have the same energy depends on the cavity size and lies between 5 and 12.5. We validate these results with QMC for the two limiting values of the dielectric constant, and find qualitative agreement with DFT even though the solvent polarization is less pronounced at the QMC level.

Car-Parrinello simulation of hydrogen bond dynamics in sodium hydrogen bissulfate

We studied proton dynamics of a short hydrogen bond of the crystalline sodium hydrogen bissulfate, a hydrogen-bonded ferroelectric system. Our approach was based on the established Car-Parrinello molecular dynamics (CPMD) methodology, followed by an a posteriori quantization of the OH stretching motion. The latter approach is based on snapshot structures taken from CPMD trajectory, calculation of proton potentials, and solving of the vibrational Schrodinger equation for each of the snapshot potentials. The so obtained contour of the OH stretching band has the center of gravity at about 1540 cm(-1) and a half width of about 700 cm(-1), which is in qualitative agreement with the experimental infrared spectrum. The corresponding values for the deuterated form are 1092 and 600 cm(-1), respectively. The hydrogen probability densities obtained by solving the vibrational Schrodinger equation allow for the evaluation of potential of mean force along the proton transfer coordinate. We demonstrate that for the present system the free energy profile is of the single-well type and features a broad and shallow minimum near the center of the hydrogen bond, allowing for frequent and barrierless proton (or deuteron) jumps. All the calculated time-averaged geometric parameters were in reasonable agreement with the experimental neutron diffraction data. As the present methodology for quantization of proton motion is applicable to a variety of hydrogen-bonded systems, it is promising for potential use in computational enzymology.

Computational Modeling of the Nitrile Stretching Vibration of 5-Cyanoindole in Water

The bandwidth of the nitrile (C≡N) stretching vibration of 5-cyanotryptophan shows a significant broadening upon hydration. Thus, it has been proposed to be a useful infrared probe of the local hydration environment of proteins. However, the molecular mechanism underlying this hydration-induced spectral broadening is not known, making interpretation of the experimental results difficult. Herein, we investigate how interactions of water with various sites of 5-cyanoindole, the sidechain of 5-cyanotryptophan, affect its C≡N stretching vibration via a combined electronic structure/molecular dynamics approach. It is found that, besides those interactions with the nitrile group, interactions of water with the indole ring also play a significant role in mediating the C≡N stretching frequency. Thus, this study provides a molecular basis for understanding how hydration affects the C≡N stretching band of 5-cyanotryptophan. In addition, an empirical model, which includes interactions of water with both the nitrile and indole groups, is developed for predicting the C≡N stretching vibrational band via molecular dynamics simulations.

Energetics of the cleft closing transition and the role of electrostatic interactions in conformational rearrangements of the glutamate receptor ligand binding domain

The ionotropic glutamate receptors are localized in the pre- and postsynaptic membrane of neurons in the brain. Activation by the principal excitatory neurotransmitter glutamate allows the ligand binding domain to change conformation, communicating opening of the channel for ion conduction. The free energy of the GluR2 S1S2 ligand binding domain (S1S2) closure transition was computed using a combination of thermodynamic integration and umbrella sampling modeling methods. A path that involves lowering the charge on E705 was chosen to clarify the role of this binding site residue. A continuum electrostatics approach in S1S2 is used to show E705, located in the ligand binding cleft, stabilizes the closed conformation of S1S2 via direct interactions with other protein residues, not through the ligand. In the closed conformation, in the absence of a ligand, S1S2 is somewhat more closed than what has been reported in X-ray structures. A semiopen conformation has been identified which is characterized by disruption of a single cross-cleft interaction and differs only slightly in energy from the fully closed S1S2. The fully open S1S2 conformation exhibits a wide energy well and shares structural similarity with the apo S1S2 crystal structure. Hybrid continuum electrostatics/MD calculations along the chosen closure transition pathway reveal solvation energies, and electrostatic interaction energies between two lobes of the protein increase the relative energetic difference between the open and closed conformational states. By analyzing the role of several cross-cleft contacts as well as other binding site residues, we demonstrate how S1S2 interactions facilitate formation of the closed conformation of the GluR2 ligand binding domain.

Approaches for the calculation of vibrational frequencies in liquids: comparison to benchmarks for azide/water clusters

Ultrafast vibrational spectroscopy experiments, together with molecular-level theoretical interpretation, can provide important information about the structure and dynamics of complex condensed phase systems, including liquids. The theoretical challenge is to calculate the instantaneous vibrational frequencies of a molecule in contact with a molecular environment, accurately and quickly, and to this end a number of different methods have been developed. In this paper we critically analyze these different methods by comparing their results to accurate benchmark calculations on azide/water clusters. We also propose an optimized quantum mechanics/molecular mechanics method, which for this problem is superior to the other methods.

Protein flexibility using constraints from molecular dynamics simulations

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.

On the Binding Determinants of the Glutamate Agonist with the Glutamate Receptor Ligand Binding Domain

Ionotropic glutamate receptors (GluRs) are ligand-gated membrane channel proteins found in the central neural system that mediate a fast excitatory response of neurons. In this paper, we report theoretical analysis of the ligand-protein interactions in the binding pocket of the S1S2 (ligand binding) domain of the GluR2 receptor in the closed conformation. By utilizing several theoretical methods ranging from continuum electrostatics to all-atom molecular dynamics simulations and quantum chemical calculations, we were able to characterize in detail glutamate agonist binding to the wild-type and E705D mutant proteins. A theoretical model of the protein-ligand interactions is validated via direct comparison of theoretical and Fourier transform infrared spectroscopy (FTIR) measured frequency shifts of the ligand's carboxylate group vibrations [Jayaraman et al. (2000) Biochemistry 39, 8693-8697; Cheng et al. (2002) Biochemistry 41, 1602-1608]. A detailed picture of the interactions in the binding site is inferred by analyzing contributions to vibrational frequencies produced by protein residues forming the ligand-binding pocket. The role of mobility and hydrogen-bonding network of water in the ligand-binding pocket and the contribution of protein residues exposed in the binding pocket to the binding and selectivity of the ligand are discussed. It is demonstrated that the molecular surface of the protein in the ligand-free state has mainly positive electrostatic potential attractive to the negatively charged ligand, and the potential produced by the protein in the ligand-binding pocket in the closed state is complementary to the distribution of the electrostatic potential produced by the ligand itself. Such charge complementarity ensures specificity to the unique charge distribution of the ligand.