The structure and backbone dynamics of an extended second transmembrane segment (TM2e) of the human neuronal glycine receptor alpha(1) subunit in sodium dodecyl sulfate micelles were studied by (1)H and (15)N solution-state NMR. The 28-amino acid segment contained the consensus TM2 domain plus part of the linker between the second and third transmembrane domains. The presence of a well-structured helical region of at least 13 amino acids long and an unstructured region near the linker was evident from the proton chemical shifts and the pattern of midrange nuclear Overhauser effects (NOE). (15)N relaxation rate constants, R(1) and R(2), and (15)N-[(1)H] NOE indicated restricted internal motions in the helical region with NOE values between 0.6 and 0.8. The squared order parameter (S(2)), the effective correlation time for fast internal motions (tau(e)), and the global rotational correlation time (tau(m)) were calculated for all TM2e backbone N-H bonds using the model-free approach. The S(2) values ranged about 0.75-0.86, and the tau(e) values were below 100 ps for most of the residues in the helical region. The tau(m) value, calculated from the dynamics of the helical region, was 5.1 ns. The S(2) values decreased to 0.1, and the tau(e) values sharply increased up to 1.2 ns at the linker near the C-terminus, indicating that the motion of this region is unrestricted. The results suggest a relatively high degree of motional freedom of TM2e in micelles and different propensities of the N- and C-terminal moieties of the transmembrane domain to assume stable helical structures.

The potential of mean forces (PMF) governing Na+ permeation through gramicidin A (gA) channels with explicit water and membrane was characterized using steered molecular dynamics (SMD) simulations. Constant-force SMD with a steering force parallel to the channel axis revealed at least seven energy wells in each monomer of the channel dimer. Except at the channel dimer interface, each energy well is associated with at least three and at most four backbone carbonyl oxygens and two water oxygens in a pseudo-hexahedral or pseudo-octahedral coordination with the Na+ ion. Repeated constant-velocity SMD by dragging a Na+ ion from each energy well in opposite directions parallel to the channel axis allowed the computation of the PMF across the gA channel, revealing a global minimum corresponding to Na+ binding sites near the entrance of gA at +/-9.3 A from the geometric center of the channel. The effect of volatile anesthetics on the PMF was also analyzed in the presence of halothane molecules. Although the accuracy of the current PMF calculation from SMD simulations is not yet sufficient to quantify the PMF difference with and without anesthetics, the comparison of the overall PMF profiles nevertheless confirms that the anesthetics cause insignificant changes to the structural makeup of the free energy wells along the channel and the overall permeation barrier. On average, the PMF appears less rugged in the outer part of the channel in the presence of anesthetics, consistent with our earlier finding that halothane interaction with anchoring residues makes the gA channel more dynamic. A causal relationship was observed between the reorientation of the coordinating backbone carbonyl oxygen and Na+ transit from one energy well to another, suggesting the possibility that even minute changes in the conformation of pore-lining residues due to dynamic motion could be sufficient to trigger the ion permeation. Because some of the carbonyl oxygens contribute to Na+ coordination in two adjacent energy wells, our SMD results reveal that the atomic picture of ion "hopping" through a gA channel actually involves a Na+ ion being carried in a relay by the coordinating oxygens from one energy well to the next. Steered molecular dynamics complements other computational approaches as an attractive means for the atomistic interpretation of experimental permeation studies.

Ab initio and empirical methods were combined to optimize the partial atomic charges and Lennard-Jones parameters for two halogenated compounds, halothane (CF3CHClBr, a potent volatile anesthetic) and hexafluoroethane (CF3CF3, a nonanesthetic). Charge optimization was achieved using empirical calculations by systematically adjusting the charge assignments to fit minimum interaction energies and geometries between a TIP3 water molecule and the halogenated compounds to the corresponding values from the ab initio calculations, which were carried out at the HF/6-311+G(2d,p) and HF/6-31G(d) levels for halothane and hexafluoroethane, respectively. To optimize the Lennard-Jones parameters, the initial estimates were obtained from scaling the values from the ab initio minimum interaction energies and geometries between neon and the halogenated compounds calculated at the MP3/6-311++G(3d,3p) level. The Lennard-Jones parameters were further refined by fitting the empirical interaction energies to the corresponding ab initio values. The refined parameters were finalized by reproducing experimental values of the heats of vaporization and densities for liquid halothane and hexafluoroethane, using molecular dynamics simulations. The calculated heats of vaporization and liquid densities using the optimized parameters are in excellent agreement with the experimental values. The results indicate that the combination of ab initio and empirical approaches works well for obtaining the nonbonded parameters of molecules with heavy halogen atoms, such as Cl and Br. The refined nonbonded parameters are readily applicable in molecular dynamics simulations involving these halogenated compounds.

A 61-residue polypeptide resembling the second and third transmembrane domains (TM23) of the alpha-1 subunit of human glycine receptor and its truncated form, both with the wild-type loop linking the two TM domains (the "23" loop), were studied using high-resolution NMR. Well-defined domain structures can be identified for the TM2, 23 loop, and TM3 regions. Contrary to the popular model of a long and straight alpha-helical structure for the pore-lining TM2 domain for the Cys-loop receptor family, the last three residues of the TM2 domain and the first eight residues of the 23 loop (S16-S26) seem to be intrinsically nonhelical and highly flexible even in trifluoroethanol, a solvent known to promote and stabilize alpha-helical structures. The six remaining residues of the 23 loop and most of the TM3 domain exhibit helical structures with a kinked pi-helix (or a pi-turn) from W34 to C38 and a kink angle of 159 +/- 3 degrees . The tertiary fold of TM3 relative to TM2 is defined by several unambiguously identified long-range NOE cross-peaks within the loop region and between TM2 and TM3 domains. The 20 lowest-energy structures show a left-handed tilt of TM3 relative to TM2 with a tilting angle of 44 +/- 2 degrees between TM2 (V1-Q14) and TM3 (L39-E48) helix axes. This left-handed TM2-TM3 arrangement ensures a neatly packed right-handed quaternary structure of five subunits to form an ion-conducting pore. This is the first time that two TM domains of the glycine receptor linked by the important 23 loop have ever been analyzed at atomistic resolution. Many structural characteristics of the receptor can be inferred from the structural and dynamical features identified in this study.