Seungho Choe

The Mechanism of Sodium and Substrate Release from the Binding Pocket of Vsglt

Membrane co-transport proteins that use a five-helix inverted repeat motif have recently emerged as one of the largest structural classes of secondary active transporters. However, despite many structural advances there is no clear evidence of how ion and substrate transport are coupled. Here we report a comprehensive study of the sodium/galactose transporter from Vibrio parahaemolyticus (vSGLT), consisting of molecular dynamics simulations, biochemical characterization and a new crystal structure of the inward-open conformation at a resolution of 2.7 Å. Our data show that sodium exit causes a reorientation of transmembrane helix 1 that opens an inner gate required for substrate exit, and also triggers minor rigid-body movements in two sets of transmembrane helical bundles. This cascade of events, initiated by sodium release, ensures proper timing of ion and substrate release. Once set in motion, these molecular changes weaken substrate binding to the transporter and allow galactose readily to enter the intracellular space. Additionally, we identify an allosteric pathway between the sodium-binding sites, the unwound portion of transmembrane helix 1 and the substrate-binding site that is essential in the coupling of co-transport.

A continuum method for determining membrane protein insertion energies and the problem of charged residues

Continuum electrostatic approaches have been extremely successful at describing the charged nature of soluble proteins and how they interact with binding partners. However, it is unclear whether continuum methods can be used to quantitatively understand the energetics of membrane protein insertion and stability. Recent translation experiments suggest that the energy required to insert charged peptides into membranes is much smaller than predicted by present continuum theories. Atomistic simulations have pointed to bilayer inhomogeneity and membrane deformation around buried charged groups as two critical features that are neglected in simpler models. Here, we develop a fully continuum method that circumvents both of these shortcomings by using elasticity theory to determine the shape of the deformed membrane and then subsequently uses this shape to carry out continuum electrostatics calculations. Our method does an excellent job of quantitatively matching results from detailed molecular dynamics simulations at a tiny fraction of the computational cost. We expect that this method will be ideal for studying large membrane protein complexes.

Water Permeation Through the Sodium-Dependent Galactose Cotransporter VSGLT

It is well accepted that cotransporters facilitate water movement by two independent mechanisms: osmotic flow through a water channel in the protein and flow driven by ion/substrate cotransport. However, the molecular mechanism of transport-linked water flow is controversial. Some researchers believe that it occurs via cotransport, in which water is pumped along with the transported cargo, while others believe that flow is osmotic in response to an increase in intracellular osmolarity. In this letter, we report the results of a 200-ns molecular dynamics simulation of the sodium-dependent galactose cotransporter vSGLT. Our simulation shows that a significant number of water molecules cross the protein through the sugar-binding site in the presence as well as the absence of galactose, and 70-80 water molecules accompany galactose as it moves from the binding site into the intracellular space. During this event, the majority of water molecules in the pathway are unable to diffuse around the galactose, resulting in water in the inner half of the transporter being pushed into the intracellular space and replaced by extracellular water. Thus, our simulation supports the notion that cotransporters act as both passive water channels and active water pumps with the transported substrate acting as a piston to rectify the motion of water.

Structural Determinants of Water Permeation through the Sodium-Galactose Transporter vSGLT

Sodium-glucose transporters (SGLTs) facilitate the movement of water across the cell membrane, playing a central role in cellular homeostasis. Here, we present a detailed analysis of the mechanism of water permeation through the inward-facing state of vSGLT based on nearly 10 μs of molecular dynamics simulations. These simulations reveal the transient formation of a continuous water channel through the transporter that permits water to permeate the protein. Trajectories in which spontaneous release of galactose is observed, as well as those in which galactose remains in the binding site, show that the permeation rate, although modulated by substrate occupancy, is not tightly coupled to substrate release. Using a, to our knowledge, novel channel-detection algorithm, we identify the key residues that control water flow through the transporter and show that solvent gating is regulated by side-chain motions in a small number of residues on the extracellular face. A sequence alignment reveals the presence of two insertion sites in mammalian SGLTs that flank these outer-gate residues. We hypothesize that the absence of these sites in vSGLT may account for the high water permeability values for vSGLT determined via simulation compared to the lower experimental estimates for mammalian SGLT1.

Stochastic steps in secondary active sugar transport

Significance The potential energy stored in ion gradients across cell membranes drives nutrients in and out of cells by cotransport proteins, e.g., uphill glucose accumulation in cells by sodium cotransporters. Insight into the mechanism of cotransport has been obtained from high-resolution atomic structures of the transporters, but further progress requires dynamic information about ion and substrate movements through the proteins. We have used multiple long molecular-dynamic simulations and electrophysiological assays to explore the dynamics of the transport cycle. Ligands bound to sodium-dependent glucose transporters are released to the cytoplasm stochastically, whereas release to the external solution is ordered with sugar first. The order of events is intimately tied to how the protein converts the energy stored in an ion gradient into a sugar gradient. Secondary active transporters, such as those that adopt the leucine-transporter fold, are found in all domains of life, and they have the unique capability of harnessing the energy stored in ion gradients to accumulate small molecules essential for life as well as expel toxic and harmful compounds. How these proteins couple ion binding and transport to the concomitant flow of substrates is a fundamental structural and biophysical question that is beginning to be answered at the atomistic level with the advent of high-resolution structures of transporters in different structural states. Nonetheless, the dynamic character of the transporters, such as ion/substrate binding order and how binding triggers conformational change, is not revealed from static structures, yet it is critical to understanding their function. Here, we report a series of molecular simulations carried out on the sugar transporter vSGLT that lend insight into how substrate and ions are released from the inward-facing state of the transporter. Our simulations reveal that the order of release is stochastic. Functional experiments were designed to test this prediction on the human homolog, hSGLT1, and we also found that cytoplasmic release is not ordered, but we confirmed that substrate and ion binding from the extracellular space is ordered. Our findings unify conflicting published results concerning cytoplasmic release of ions and substrate and hint at the possibility that other transporters in the superfamily may lack coordination between ions and substrate in the inward-facing state.

Conformational dynamics of the inner pore helix of voltage-gated potassium channels

Voltage-gated potassium (Kv) channels control the electrical excitability of neurons and muscles. Despite this key role, how these channels open and close or gate is not fully understood. Gating is usually attributed to the bending and straightening of pore-lining helices at glycine and proline residues. In this work we focused on the role of proline in the Pro-Val-Pro (PVP) motif of the inner S6 helix in the Kv1.2 channel. We started by developing a simple hinged-rod model to fully explore the configurational space of bent helices and we related these configurations to the degree of pore opening. We then carried out fully atomistic simulations of the S6 helices and compared these simulations to the hinged-rod model. Both methods suggest that Kv1 channels are not tightly closed when the inner helices are straight, unlike what is seen in the non-PVP containing channels KcsA and KirBac. These results invite the possibility that the S6 helices may be kinked when Kv1 channels are closed. Our simulations indicate that the wild-type helix adopts multiple spatially distinct configurations, which is consistent with its role in adopting a closed state and an open state. The two most dominant configurational basins correspond to a 6 A movement of the helix tail accompanied by the PVP region undergoing a local alpha-helix to 3(10)-helix transition. We explored how single point mutations affect the propensity of the S6 helix to adopt particular configurations. Interestingly, mutating the first proline, P405 (P473 in Shaker), to alanine completely removed the bistable nature of the S6 helix possibly explaining why this mutation compromises the channel. Next, we considered four other mutations in the area known to affect channel gating and we saw similarly dramatic changes to the helixs dynamics and range of motion. Our results suggest a possible mechanism of helix pore closure and they suggest differences in the closed state of glycine-only channels, like KcsA, and PVP containing channels.

CMB spectral mu-distortion of multiple inflation scenario

Abstract In multiple inflation scenario having two inflations with an intermediate matter-dominated phase, the power spectrum is estimated to be enhanced on scales smaller than the horizon size at the beginning of the second inflation, k > k b . We require k b > 10 Mpc − 1 to make sure that the enhanced power spectrum is consistent with large scale observation of cosmic microwave background (CMB). We consider the CMB spectral distortions generated by the dissipation of acoustic waves to constrain the power spectrum. The μ-distortion value can be 102 times larger than the expectation of the standard ΛCDM model ( μ Λ CDM ≃ 2 × 10 − 8 ) for k b ≲ 10 3 Mpc − 1 , while the y-distortion is hardly affected by the enhancement of the power spectrum.