Search

Atypical Chemokine Receptor 3 (ACKR3) belongs to the G protein-coupled receptor family but it does not signal through G proteins. The structural properties that govern the functional selectivity and the conformational dynamics of ACKR3 activation are poorly understood. Here, we combined hydrogen/deuterium exchange mass spectrometry, site-directed mutagenesis, and molecular dynamics simulations to examine the binding mode and mechanism of action of ACKR3 ligands of different efficacies. Our results show that activation or inhibition of ACKR3 is governed by intracellular conformational changes of helix 6, intracellular loop 2, and helix 7, while the DRY motif becomes protected during both processes. Moreover, we identified the binding sites and the allosteric modulation of ACKR3 upon β-arrestin 1 binding. In summary, this study highlights the structure-function relationship of small ligands, the binding mode of β-arrestin 1, the activation dynamics, and the atypical dynamic features in ACKR3 that may contribute to its inability to activate G proteins. [ABSTRACT FROM AUTHOR]

Cryptococcus neoformans is a fungal pathogen responsible for cryptococcosis and cryptococcal meningitis. The C. neoformans' capsular polysaccharide and its shed exopolysaccharide function both as key virulence factors and to protect the fungal cell from phagocytosis. Currently, a glycoconjugate of these polysaccharides is being explored as a vaccine to protect against C. neoformans infection. In this study, NOE and J-coupling values from NMR experiments were consistent with a converged structure of the synthetic decasaccharide, GXM10-Ac3, calculated from MD simulations. GXM10-Ac3 was designed as an extension of glucuronoxylomannan (GXM) polysaccharide motif (M2) which is common in the clinically predominant serotype A strains and is recognized by protective forms of GXM-specific monoclonal antibodies. The M2 motif is a hexasaccharide with a three-residue α-mannan backbone, modified by β-(1→2)-xyloses (Xyl) on the first two mannoses (Man) and a β-(1→2)-glucuronic acid (GlcA) on the third Man. Combined NMR and MD analyses reveal that GXM10-Ac3 adopts an extended structure, with Xyl/GlcA branches alternating sides along the α-mannan backbone. O-acetyl esters also alternate sides and are grouped in pairs. MD analysis of a twelve M2-repeating unit polymer supports the notion that the GXM10-Ac3 structure is uniformly represented throughout the polysaccharide. This derived GXM model displays high flexibility while maintaining a structural identity, yielding insights to further explore intermolecular interactions between polysaccharides, interactions with anti-GXM mAbs, and the cryptococcal polysaccharide architecture. [ABSTRACT FROM AUTHOR]

Prostaglandin D2 (PGD2) signals through the G protein-coupled receptor (GPCR) CRTH2 to mediate various inflammatory responses. CRTH2 is the only member of the prostanoid receptor family that is phylogenetically distant from others, implying a nonconserved mechanism of lipid action on CRTH2. Here, we report a crystal structure of human CRTH2 bound to a PGD2 derivative, 15R-methyl-PGD2 (15mPGD2), by serial femtosecond crystallography. The structure revealed a "polar group in"-binding mode of 15mPGD2 contrasting the "polar group out"-binding mode of PGE2 in its receptor EP3. Structural comparison analysis suggested that these two lipid-binding modes, associated with distinct charge distributions of ligand-binding pockets, may apply to other lipid GPCRs. Molecular dynamics simulations together with mutagenesis studies also identified charged residues at the ligand entry port that function to capture lipid ligands of CRTH2 from the lipid bilayer. Together, our studies suggest critical roles of charge environment in lipid recognition by GPCRs.

Conformational dynamics play essential roles in RNA function. However, detailed structural characterization of excited states of RNA remains challenging. Here, we apply high hydrostatic pressure (HP) to populate excited conformational states of tRNALys3, and structurally characterize them using a combination of HP 2D-NMR, HP-SAXS (HP-small-angle X-ray scattering), and computational modeling. HP-NMR revealed that pressure disrupts the interactions of the imino protons of the uridine and guanosine U-A and G-C base pairs of tRNALys3. HP-SAXS profiles showed a change in shape, but no change in overall extension of the transfer RNA (tRNA) at HP. Configurations extracted from computational ensemble modeling of HP-SAXS profiles were consistent with the NMR results, exhibiting significant disruptions to the acceptor stem, the anticodon stem, and the D-stem regions at HP. We propose that initiation of reverse transcription of HIV RNA could make use of one or more of these excited states. [ABSTRACT FROM AUTHOR]

Despite the negative charge of the DNA backbone, acidic residues (Asp/Glu) commonly participate in the base readout, with a strong preference for cytosine. In fact, in the solved DNA/protein structures, cytosine is recognized almost exclusively by Asp/Glu through a direct hydrogen bond, while at the same time, adenine, regardless of its amino group, shows no propensity for Asp/Glu. Here, we analyzed the contribution of Asp/Glu to sequence-specific DNA binding using classical and ab initio simulations of selected transcription factors and found that it is governed by a fine balance between the repulsion from backbone phosphates and attractive interactions with cytosine. Specifically, Asp/Glu lower the affinity for noncytosine sites and thus act as negative selectors preventing off-target binding. At cytosine-containing sites, the favorable contribution does not merely rely on the formation of a single H-bond but usually requires the presence of positive potential generated by multiple cytosines, consistently with the observed excess of cytosine in the target sites. Finally, we show that the preference of Asp/Glu for cytosine over adenine is a result of the repulsion from the adenine imidazole ring and a tendency of purine--purine dinucleotides to adopt the BII conformation. [ABSTRACT FROM AUTHOR]

DNA recombination is a ubiquitous process that ensures genetic diversity. Contrary to textbook pictures, DNA recombination, as well as generic DNA translocations, occurs in a confined and highly entangled environment. Inspired by this observation, here, we investigate a solution of semiflexible polymer rings undergoing generic cutting and reconnection operations under spherical confinement. Our setup may be realized using engineered DNA in the presence of recombinase proteins or by considering micelle-like components able to form living (or reversibly breakable) polymer rings. We find that in such systems, there is a topological gelation transition, which can be triggered by increasing either the stiffness or the concentration of the rings. Flexible or dilute polymers break into an ensemble of short, unlinked, and segregated rings, whereas sufficiently stiff or dense polymers self-assemble into a network of long, linked, and mixed loops, many of which are knotted. We predict that the two phases should behave qualitatively differently in elution experiments monitoring the escape dynamics from a permeabilized container. Besides shedding some light on the biophysics and topology of genomes undergoing DNA reconnection in vivo, our findings could be leveraged in vitro to design polymeric complex fluids--e.g., DNA-based complex fluids or living polymer networks--with desired topologies. [ABSTRACT FROM AUTHOR]

The use of lipid nanoparticles (LNPs) for therapeutic RNA delivery has gained significant interest, particularly highlighted by recent milestones such as the approval of Onpattro and two mRNA-based SARS-CoV-2 vaccines. However, despite substantial advancements in this field, our understanding of the structure and internal organization of RNA-LNPs -and their relationship to efficacy, both in vitro and in vivo- remains limited. In this study, we present a coarse-grained molecular dynamics (MD) approach that allows for the simulations of full-size LNPs. By analyzing MD-derived structural characteristics in conjunction with cellular experiments, we investigate the effect of critical parameters, such as pH and composition, on LNP structure and potency. Additionally, we examine the mobility and chemical environment within LNPs at a molecular level. Our findings highlight the significant impact that LNP composition and internal molecular mobility can have on key stages of LNP-based intracellular RNA delivery., Competing Interests: Competing interests statement:All authors are current employees of Bayer AG or Nuvisan ICB GmbH.

Dynamic motions of enzymes occurring on a broad range of timescales play a pivotal role in all steps of the reaction pathway, including substrate binding, catalysis, and product release. However, it is unknown whether structural information related to conformational flexibility can be exploited for the directed evolution of enzymes with higher catalytic activity. Here, we show that mutagenesis of residues exclusively located at flexible regions distal to the active site of Homo sapiens kynureninase (HsKYNase) resulted in the isolation of a variant (BF-HsKYNase) in which the rate of the chemical step toward kynurenine was increased by 45-fold. Mechanistic pre--steady-state kinetic analysis of the wild type and the evolved enzyme shed light on the underlying effects of distal mutations (>10 Å from the active site) on the rate-limiting step of the catalytic cycle. Hydrogen-deuterium exchange coupled to mass spectrometry and molecular dynamics simulations revealed that the amino acid substitutions in BF-HsKYNase allosterically affect the flexibility of the pyridoxal-50 -phosphate (PLP) binding pocket, thereby impacting the rate of chemistry, presumably by altering the conformational ensemble and sampling states more favorable to the catalyzed reaction. [ABSTRACT FROM AUTHOR]

ASCT2 (SLC1A5) is a sodium-dependent neutral amino acid transporter that controls amino acid homeostasis in peripheral tissues. In cancer, ASCT2 is up-regulated where it modulates intracellular glutamine levels, fueling cell proliferation. Nutrient deprivation via ASCT2 inhibition provides a potential strategy for cancer therapy. Here, we rationally designed stereospecific inhibitors exploiting specific subpockets in the substrate binding site using computational modeling and cryo-electron microscopy (cryo-EM). The final structures combined with molecular dynamics simulations reveal multiple pharmacologically relevant conformations in the ASCT2 binding site as well as a previously unknown mechanism of stereospecific inhibition. Furthermore, this integrated analysis guided the design of a series of unique ASCT2 inhibitors. Our results provide a framework for future development of cancer therapeutics targeting nutrient transport via ASCT2, as well as demonstrate the utility of combining computational modeling and cryo-EM for solute carrier ligand discovery. [ABSTRACT FROM AUTHOR]

Defining the molecular details and consequences of the association of water-soluble proteins with membranes is fundamental to understanding protein-lipid interactions and membrane functioning. Phospholipase A2 (PLA2) enzymes, which catalyze the hydrolysis of phospholipid substrates that compose the membrane bilayers, provide the ideal system for studying protein-lipid interactions. Our study focuses on understanding the catalytic cycle of two different human PLA2s: the cytosolic Group IVA cPLA2 and calcium-independent Group VIA iPLA2. Computer-aided techniques guided by deuterium exchange mass spectrometry data, were used to create structural complexes of each enzyme with a single phospholipid substrate molecule, whereas the substrate extraction process was studied using steered molecular dynamics simulations. Molecular dynamic simulations of the enzyme-substrate-membrane systems revealed important information about the mechanisms by which these enzymes associate with the membrane and then extract and bind their phospholipid substrate. Our data support the hypothesis that the membrane acts as an allosteric ligand that binds at the allosteric site of the enzyme's interfacial surface, shifting its conformation from a closed (inactive) state in water to an open (active) state at the membrane interface.

Prostaglandin D2 (PGD2) signals through the G protein-coupled receptor (GPCR) CRTH2 to mediate various inflammatory responses. CRTH2 is the only member of the prostanoid receptor family that is phylogenetically distant from others, implying a nonconserved mechanism of lipid action on CRTH2. Here, we report a crystal structure of human CRTH2 bound to a PGD2 derivative, 15R-methyl-PGD2 (15mPGD2), by serial femtosecond crystallography. The structure revealed a "polar group in"-binding mode of 15mPGD2 contrasting the "polar group out"-binding mode of PGE2 in its receptor EP3. Structural comparison analysis suggested that these two lipid-binding modes, associated with distinct charge distributions of ligand-binding pockets, may apply to other lipid GPCRs. Molecular dynamics simulations together with mutagenesis studies also identified charged residues at the ligand entry port that function to capture lipid ligands of CRTH2 from the lipid bilayer. Together, our studies suggest critical roles of charge environment in lipid recognition by GPCRs. [ABSTRACT FROM AUTHOR]

Allostery can be manifested as a combination of repression and activation in multidomain proteins allowing for fine tuning of regulatory mechanisms. Here we have used single molecule fluorescence resonance energy transfer (smFRET) and molecular dynamics simulations to study the mechanism of allostery underlying negative cooperativity between the two agonists glutamate and glycine in the NMDA receptor. These data show that binding of one agonist leads to conformational flexibility and an increase in conformational spread at the second agonist site. Mutational and cross-linking studies show that the dimer-dimer interface at the agonist-binding domain mediates the allostery underlying the negative cooperativity. smFRET on the transmembrane segments shows that they are tightly coupled in the unliganded and single agonist-bound form and only upon binding both agonists the transmembrane domain explores looser packing which would facilitate activation. [ABSTRACT FROM AUTHOR]

Small-angle X-ray scattering (SAXS) measurements reveal a striking difference in intermolecular interactions between two short highly charged peptides--deca-arginine (R10) and deca-lysine (K10). Comparison of SAXS curves at high and low salt concentration shows that R10 self-associates, while interactions between K10 chains are purely repulsive. The self-association of R10 is stronger at lower ionic strengths, indicating that the attraction between R10 molecules has an important electrostatic component. SAXS data are complemented by NMR measurements and potentials of mean force between the peptides, calculated by means of umbrella-sampling molecular dynamics (MD) simulations. All-atom MD simulations elucidate the origin of the R10-R10 attraction by providing structural information on the dimeric state. The last two C-terminal residues of R10 constitute an adhesive patch formed by stacking of the side chains of two arginine residues and by salt bridges formed between the like-charge ion pair and the C-terminal carboxyl groups. A statistical analysis of the Protein Data Bank reveals that this mode of interaction is a common feature in proteins. [ABSTRACT FROM AUTHOR]

We report experimental and computational studies investigating the effects of three osmolytes, trimethylamine N-oxide (TMAO), betaine, and glycine, on the hydrophobic collapse of an elastin-like polypeptide (ELP). All three osmolytes stabilize collapsed conformations of the ELP and reduce the lower critical solution temperature (LSCT) linearly with osmolyte concentration. As expected from conventional preferential solvation arguments, betaine and glycine both increase the surface tension at the air-water interface. TMAO, however, reduces the surface tension. Atomically detailed molecular dynamics (MD) simulations suggest that TMAO also slightly accumulates at the polymer-water interface, whereas glycine and betaine are strongly depleted. To investigate alternative mechanisms for osmolyte effects, we performed FTIR experiments that characterized the impact of each cosolvent on the bulk water structure. These experiments showed that TMAO red-shifts the OH stretch of the IR spectrum via a mechanism that was very sensitive to the protonation state of the NO moiety. Glycine also caused a red shift in the OH stretch region, whereas betaine minimally impacted this region. Thus, the effects of osmolytes on the OH spectrum appear uncorrelated with their effects upon hydrophobic collapse. Similarly, MD simulations suggested that TMAO disrupts the water structure to the least extent, whereas glycine exerts the greatest influence on the water structure. These results suggest that TMAO stabilizes collapsed conformations via amechanismthat is distinct from glycine and betaine. In particular, we propose that TMAO stabilizes proteins by acting as a surfactant for the heterogeneous surfaces of folded proteins. [ABSTRACT FROM AUTHOR]

Conformational dynamics play essential roles in RNA function. However, detailed structural characterization of excited states of RNA remains challenging. Here, we apply high hydrostatic pressure (HP) to populate excited conformational states of tRNA Lys3 , and structurally characterize them using a combination of HP 2D-NMR, HP-SAXS (HP-small-angle X-ray scattering), and computational modeling. HP-NMR revealed that pressure disrupts the interactions of the imino protons of the uridine and guanosine U-A and G-C base pairs of tRNA Lys3 . HP-SAXS profiles showed a change in shape, but no change in overall extension of the transfer RNA (tRNA) at HP. Configurations extracted from computational ensemble modeling of HP-SAXS profiles were consistent with the NMR results, exhibiting significant disruptions to the acceptor stem, the anticodon stem, and the D-stem regions at HP. We propose that initiation of reverse transcription of HIV RNA could make use of one or more of these excited states.

The simple process of a liquid wetting a solid surface is controlled by a plethora of factors--surface texture, liquid droplet size and shape, energetics of both liquid and solid surfaces, as well as their interface. Studying these events at the nanoscale provides insights into the molecular basis of wetting. Nanotube wetting studies are particularly challenging due to their unique shape and small size. Nonetheless, the success of nanotubes, particularly inorganic ones, as fillers in composite materials makes it essential to understand how common liquids wet them. Here, we present a comprehensive wetting study of individual tungsten disulfide nanotubes by water. We reveal the nature of interaction at the inert outer wall and show that remarkably high wetting forces are attained on small, open-ended nanotubes due to capillary aspiration into the hollow core. This study provides a theoretical and experimental paradigm for this intricate problem. [ABSTRACT FROM AUTHOR]

Using all-atom molecular dynamics (MD) simulations at constant water chemical potential in combination with basic theoretical arguments, we study hydration-induced interactions between two overall charge-neutral yet polar planar surfaces with different wetting properties. Whether the water film between the two surfaces becomes unstable below a threshold separation and cavitation gives rise to long-range attraction, depends on the sum of the two individual surface contact angles. Consequently, cavitation-induced attraction also occurs for a mildly hydrophilic surface interacting with a very hydrophobic surface. If both surfaces are very hydrophilic, hydration repulsion dominates at small separations and direct attractive force contribution can--if strong enough--give rise to wet adhesion in this case. In between the regimes of cavitation-induced attraction and hydration repulsion we find a narrow range of contact angle combinations where the surfaces adhere at contact in the absence of cavitation. This dry adhesion regime is driven by direct surface-surface interactions. We derive simple laws for the cavitation transition as well as for the transition between hydration repulsion and dry adhesion, which favorably compare with simulation results in a generic adhesion state diagram as a function of the two surface contact angles. [ABSTRACT FROM AUTHOR]

Defining the molecular details and consequences of the association of water-soluble proteins with membranes is fundamental to understanding protein-lipid interactions and membrane functioning. Phospholipase A2 (PLA2) enzymes, which catalyze the hydrolysis of phospholipid substrates that compose the membrane bilayers, provide the ideal system for studying protein-lipid interactions. Our study focuses on understanding the catalytic cycle of two different human PLA2s: the cytosolic Group IVA cPLA2 and calcium-independent Group VIA iPLA2. Computer-aided techniques guided by deuterium exchange mass spectrometry data, were used to create structural complexes of each enzyme with a single phospholipid substrate molecule, whereas the substrate extraction process was studied using steered molecular dynamics simulations. Molecular dynamic simulations of the enzyme-substrate-membrane systems revealed important information about the mechanisms by which these enzymes associate with the membrane and then extract and bind their phospholipid substrate. Our data support the hypothesis that the membrane acts as an allosteric ligand that binds at the allosteric site of the enzyme's interfacial surface, shifting its conformation from a closed (inactive) state in water to an open (active) state at the membrane interface. [ABSTRACT FROM AUTHOR]

Halogenated inhaled general anesthetic agents modulate voltagegated ion channels, but the underlying molecular mechanisms are not understood. Many general anesthetic agents regulate voltagegated Na+ (NaV) channels, including the commonly used drug sevoflurane. Here, we investigated the putative binding sites and molecular mechanisms of sevoflurane action on the bacterial NaV channel NaChBac by using a combination of molecular dynamics simulation, electrophysiology, and kinetic analysis. Structural modeling revealed multiple sevoflurane interaction sites possibly associated with NaChBac modulation. Electrophysiologically, sevoflurane favors activation and inactivation at low concentrations (0.2 mM), and additionally accelerates current decay at high concentrations (2 mM). Explaining these observations, kinetic modeling suggests concurrent destabilization of closed states and low-affinity open channel block. We propose that the multiple effects of sevoflurane on NaChBac result from simultaneous interactions at multiple sites with distinct affinities. This multiple-site, multiple-mode hypothesis offers a framework to study the structural basis of general anesthetic action. [ABSTRACT FROM AUTHOR]

The M2 protein of influenza virus A is a proton-selective ion channel activated by pH. Structure determination by solid-state and solution NMR and X-ray crystallography has contributed significantly to our understanding, but channel activation may involve conformations not captured by these studies. Indeed, solid-state NMR data demonstrate that the M2 protein possesses significant conformational heterogeneity. Here, we report molecular dynamics (MD) simulations of the M2 transmembrane domain (TMD) in the absence and presence of the antiviral drug amantadine. The ensembles of MD conformations for both apo and bound forms reproduced the NMR data well. The TMD helix was found to kink around Gly-34, where water molecules penetrated deeply into the backbone. The amantadine-bound form exhibited a single peak ≈10° in the distribution of helix-kink angle, but the apo form exhibited 2 peaks, ≈0° and 40°. Conformations of the apo form with small and large kink angles had narrow and wide pores, respectively, around the primary gate formed by His-37 and Trp-41. We propose a structural model for channel activation, in which the small-kink conformations dominate before proton uptake by His-37 from the exterior, and proton uptake makes the large-kink conformations more favorable, thereby priming His-37 for proton release to the interior. [ABSTRACT FROM AUTHOR]

Protein α-helices are ubiquitous secondary structural elements, seldom considered to be stable without tertiary contacts. However, amino acid sequences in proteins that are based on alternating repeats of four glutamic acid (E) residues and four positively charged residues, a combination of arginine (R) and lysine (K), have been shown to form stable α-helices in a few proteins, in the absence of tertiary interactions. Here, we find that this ER/K motif is more prevalent than previously reported, being represented in proteins of diverse function from archaea to humans. By using molecular dynamics (MD) simulations, we characterize a dynamic pattern of side-chain interactions that extends along the backbone of ER/K α-helices. A simplified model predicts that side-chain interactions alone contribute substantial bending rigidity (0.5 pN/ nm) to ER/K α-helices. Results of small-angle x-ray scattering (SAXS) and single-molecule optical-trap analyses are consistent with the high bending rigidity predicted by our model. Thus, the ER/K α-helix is an isolated secondary structural element that can efficiently span long distances in proteins, making it a promising tool in designing synthetic proteins. We propose that the significant rigidity of the ER/K α-helix can help regulate protein function, as a force transducer between protein subdomains. [ABSTRACT FROM AUTHOR]

The short-time orientational relaxation of water is studied by ultrafast infrared pump-probe spectroscopy of the hydroxyl stretching mode (OD of dilute HOD in H2O). The anisotropy decay displays a sharp drop at very short times caused by inertial orientational motion, followed by a much slower decay that fully randomizes the orientation. Investigation of temperatures from 1°C to 65°C shows that the amplitude of the inertial component (extent of inertial angular displacement) depends strongly on the stretching frequency of the OD oscillator at higher temperatures, although the slow component is frequency-independent. The inertial component becomes frequency-independent at low temperatures. At high temperatures there is a correlation between the amplitude of the inertial decay and the strength of the O-D—O hydrogen bond, but at low temperatures the correlation disappears, showing that a single hydrogen bond (OD—O) is no longer a significant determinant of the inertial angular motion. It is suggested that the loss of correlation at lower temperatures is caused by the increased importance of collective effects of the extended hydrogen bonding network. By using a new harmonic cone model, the experimentally measured amplitudes of the inertial decays yield estimates of the characteristic frequencies of the intermolecular angular potential for various strengths of hydrogen bonds. The frequencies are in the range of ≈400 cm-1. A comparison with recent molecular dynamics simulations employing the simple point charge-extended water model at room temperature shows that the simulations qualitatively reflect the correlation between the inertial decay and the OD stretching frequency. [ABSTRACT FROM AUTHOR]

Prostaglandin D 2 (PGD 2 ) signals through the G protein-coupled receptor (GPCR) CRTH2 to mediate various inflammatory responses. CRTH2 is the only member of the prostanoid receptor family that is phylogenetically distant from others, implying a nonconserved mechanism of lipid action on CRTH2. Here, we report a crystal structure of human CRTH2 bound to a PGD 2 derivative, 15R-methyl-PGD 2 (15mPGD 2 ), by serial femtosecond crystallography. The structure revealed a "polar group in"-binding mode of 15mPGD 2 contrasting the "polar group out"-binding mode of PGE 2 in its receptor EP3. Structural comparison analysis suggested that these two lipid-binding modes, associated with distinct charge distributions of ligand-binding pockets, may apply to other lipid GPCRs. Molecular dynamics simulations together with mutagenesis studies also identified charged residues at the ligand entry port that function to capture lipid ligands of CRTH2 from the lipid bilayer. Together, our studies suggest critical roles of charge environment in lipid recognition by GPCRs., Competing Interests: The authors declare no competing interest.

We report experimental and computational studies investigating the effects of three osmolytes, trimethylamine N -oxide (TMAO), betaine, and glycine, on the hydrophobic collapse of an elastin-like polypeptide (ELP). All three osmolytes stabilize collapsed conformations of the ELP and reduce the lower critical solution temperature (LSCT) linearly with osmolyte concentration. As expected from conventional preferential solvation arguments, betaine and glycine both increase the surface tension at the air-water interface. TMAO, however, reduces the surface tension. Atomically detailed molecular dynamics (MD) simulations suggest that TMAO also slightly accumulates at the polymer-water interface, whereas glycine and betaine are strongly depleted. To investigate alternative mechanisms for osmolyte effects, we performed FTIR experiments that characterized the impact of each cosolvent on the bulk water structure. These experiments showed that TMAO red-shifts the OH stretch of the IR spectrum via a mechanism that was very sensitive to the protonation state of the NO moiety. Glycine also caused a red shift in the OH stretch region, whereas betaine minimally impacted this region. Thus, the effects of osmolytes on the OH spectrum appear uncorrelated with their effects upon hydrophobic collapse. Similarly, MD simulations suggested that TMAO disrupts the water structure to the least extent, whereas glycine exerts the greatest influence on the water structure. These results suggest that TMAO stabilizes collapsed conformations via a mechanism that is distinct from glycine and betaine. In particular, we propose that TMAO stabilizes proteins by acting as a surfactant for the heterogeneous surfaces of folded proteins.

Cyclic nucleotide-gated (CNG) ion channels, despite a significant homology with the highly selective K(+) channels, do not discriminate among monovalent alkali cations and are permeable also to several organic cations. We combined electrophysiology, molecular dynamics (MD) simulations, and X-ray crystallography to demonstrate that the pore of CNG channels is highly flexible. When a CNG mimic is crystallized in the presence of a variety of monovalent cations, including Na(+), Cs(+), and dimethylammonium (DMA(+)), the side chain of Glu66 in the selectivity filter shows multiple conformations and the diameter of the pore changes significantly. MD simulations indicate that Glu66 and the prolines in the outer vestibule undergo large fluctuations, which are modulated by the ionic species and the voltage. This flexibility underlies the coupling between gating and permeation and the poor ionic selectivity of CNG channels.

Halogenated inhaled general anesthetic agents modulate voltage-gated ion channels, but the underlying molecular mechanisms are not understood. Many general anesthetic agents regulate voltage-gated Na(+) (NaV) channels, including the commonly used drug sevoflurane. Here, we investigated the putative binding sites and molecular mechanisms of sevoflurane action on the bacterial NaV channel NaChBac by using a combination of molecular dynamics simulation, electrophysiology, and kinetic analysis. Structural modeling revealed multiple sevoflurane interaction sites possibly associated with NaChBac modulation. Electrophysiologically, sevoflurane favors activation and inactivation at low concentrations (0.2 mM), and additionally accelerates current decay at high concentrations (2 mM). Explaining these observations, kinetic modeling suggests concurrent destabilization of closed states and low-affinity open channel block. We propose that the multiple effects of sevoflurane on NaChBac result from simultaneous interactions at multiple sites with distinct affinities. This multiple-site, multiple-mode hypothesis offers a framework to study the structural basis of general anesthetic action.