Detailed mechanism for AmtB conducting NH4+/NH3: molecular dynamics simulations

The mechanism by which the ammonium transporter, AmtB, conducts NH4+/NH3 into the cytoplasm was investigated using conventional molecular dynamics (MD) simulations. These simulations revealed that the neutral molecule, NH3, passes automatically through the channel upon its arrival at the Am2 site and that the function of the channel as a one-way valve for passage of NH3 is not determined by the cytoplasmic exit gate but, rather, by the periplasmic entrance gate of the channel. The NH3, produced by deprotonation of NH4+ at the entrance gate, is spontaneously conveyed to the central region of the channel via a hydrogen-bond network comprising His-168, His-318, Tyr-32, and the NH3 molecule. Finally, the NH3 molecule exits the channel through the exit gate formed by Phe-31, Ile-266, Val-314, and His-318. In addition, Ser-263 is shown to play a critical role in the final stages, acting as a pivoting arm to shunt the NH3 molecule from the cytoplasmic exit gate of the channel out into the cytoplasm. This finding is further supported by another simulation which shows that NH3 fails to be translocated through the channel formed by the Ser-263-Ala mutation. Thus, this study casts new insights on the mechanism of AmtB-mediated passage of NH3 across cellular membranes.     

Site-directed mutagenesis in haemoglobin. Functional role of tyrosine-42(C7) alpha at the alpha 1-beta 2 interface

To clarify the functional role of Tyr-42(C7) alpha, which forms a hydrogen bond with Asp-99(G1) beta at the alpha 1-beta 2 interface of human deoxyhaemoglobin, we engineered two artificial mutant haemoglobins (Hb), in which Tyr-42 alpha was replaced by Phe (Hb Phe-42 alpha) or His (Hb His-42 alpha), and investigated their oxygen binding properties together with structural consequences of the mutations by using various spectroscopic probes. Like most of the natural Asp-99 beta mutants, Hb Phe-42 alpha showed a markedly increased oxygen affinity, a reduced Bohr effect and diminished co-operativity. Structural probes such as ultraviolet-region derivative and oxy-minus-deoxy difference spectra, resonance Raman scattering and proton nuclear magnetic resonance spectra indicate that, in Hb Phe-42 alpha, the deoxy T quaternary structure is highly destabilized and the strain imposed on the Fe-N epsilon (proximal His) bond is released, stabilizing the oxy tertiary structure. In contrast with Hb Phe-42 alpha, Hb His-42 alpha showed an intermediately impaired function and only moderate destabilization of the T-state, which can be explained by the formation of a new, weak hydrogen bond between His-42 alpha and Asp-99 beta. Spectroscopic data were consistent with this assumption. The present study proves that the hydrogen bond between Tyr-42 alpha and Asp-99 beta plays a key role in stabilizing the deoxy T structure and consequently in co-operative oxygen binding.     

Role of the carboxyl terminal region of H(+)-ATPase (F0F1) a subunit from Escherichia coli.

The effects of amino acid substitutions in the carboxyl terminal region of the H(+)-ATPase a subunit (271 amino acid residues) of Escherichia coli were studied using a defined expression system for uncB genes coded by recombinant plasmids. The a subunits with the mutations, Tyr-263----end, Trp-231----end, Glu-219----Gln, and Arg-210----Lys (or Gln) were fully defective in ATP-dependent proton translocation, and those with Gln-252----Glu (or Leu), His-245----Glu, Pro-230----Leu, and Glu-219----His were partially defective. On the other hand, the phenotypes of the Glu-269----end, Ser-265----Ala (or end), and Tyr-263----Phe mutants were essentially similar to that of the wild-type. These results suggested that seven amino acid residues between Ser-265 and the carboxyl terminus were not required for the functional proton pathway but that all the other residues except Arg-210, Glu-219, and His-245 were required for maintaining the correct conformation of the proton pathway. The results were consistent with a report that Arg-210 is directly involved in proton translocation.     

The aromatic residues of bovine pancreatic ribonuclease studied by 1H nuclear magnetic resonance.

1. The aromatic proton resonances in the 360-MHz 1H nuclear magnetic resonance (NMR) spectrum of bovine pancreatic ribonuclease were divided into histidine, tyrosine and phenylalanine resonances by means of pH titrations and double resonance experiments. 2. Photochemically induced dynamic nuclear polarization spectra showed that one histidine (His-119) and two tyrosines are accessibly to photo-excited flavin. This permitted the identification of the C-4 proton resonance of His-119. 3. The resonances of the ring protons of Tyr-25, Tyr-76 and Tyr-115 and the C-4 proton of His-12 were identified by comparison with subtilisin-modified and nitrated ribonucleases. Other resonances were assigned tentatively to Tyr-73, Tyr-92 and Phe-46. 4. On addition of active-site inhibitors, all phenylalanine resonances broadened or disappeared. The resonance that was most affected was assigned tentatively to Phe-120. 5. Four of the six tyrosines of bovine RNase, identified as Tyr-76, Tyr-115 and, tentatively, Tyr-73 and Tyr-92, are titratable above pH 9. The rings of Tyr-73 and Tyr-115 are rapidly rotating or flipping by 180 degrees about their C beta--C gamma bond and are accessible to flavin in photochemically induced dynamic nuclear polarization experiments. Tyr-25 is involved in a pH-dependent conformational transition, together with Asp-14 and His-48. A scheme for this transition is proposed. 6. Binding of active-site inhibitors to bovine RNase only influences the active site and its immediate surroundings. These conformational changes are probably not connected with the pH-dependent transition in the region of Asp-14, Tyr-25 and His-48. 7. In NMR spectra of RNase A at elevated temperatures, no local unfolding below the temperature of the thermal denaturation was observed. NMR spectra of thermally unfolded RNase A indicated that the deviations from a random coil are small and might be caused by interactions between neighbouring residues.     

Rationalization of poor solubility of TGF-β3 using MD simulation

Despite sequence and structural similarity, TGF-β3 has low solubility among other isoforms of TGF-β. We used nanosecond of molecular dynamic simulations (MD) with explicit solvent, alone and in presence of urea, to investigate the intermediates resulting from the unfolding process of TGF-β3 and TGF-β1. MD simulations of the full-length proteins show a very early loss of α-helix in TGF-β3 compared to the one in the TGF-β1. MD simulation of a small fragment consisting of H3 α-helix of TGF-β3 shows conversion of this segment to β-sheet. Relative instability of H3 α-helix in TGF-β3 and its propensity to form β-sheet may explain the poor solubility of TGF-β3 compared to TGF-β1. The other reasons for poor solubility of TGF-β3 may be the hydrophobic patches on its surface and low charge over the entire range of pH.     

Desensitization mechanism in prokaryotic ligand-gated ion channel

Crystal structures of Gloeobacter violaceus ligand-gated ion channel (GLIC), a proton-gated prokaryotic homologue of pentameric ligand-gated ion channel (LGIC) from G. violaceus, have provided high-resolution models of the channel architecture and its role in selective ion conduction and drug binding. However, it is still unclear which functional states of the LGIC gating scheme these crystal structures represent. Much of this uncertainty arises from a lack of thorough understanding of the functional properties of these prokaryotic channels. To elucidate the molecular events that constitute gating, we have carried out an extensive characterization of GLIC function and dynamics in reconstituted proteoliposomes by patch clamp measurements and EPR spectroscopy. We find that GLIC channels show rapid activation upon jumps to acidic pH followed by a time-dependent loss of conductance because of desensitization. GLIC desensitization is strongly coupled to activation and is modulated by voltage, permeant ions, pore-blocking drugs, and membrane cholesterol. Many of these properties are parallel to functions observed in members of eukaryotic LGIC. Conformational changes in loop C, measured by site-directed spin labeling and EPR spectroscopy, reveal immobilization during desensitization analogous to changes in LGIC and acetylcholine binding protein. Together, our studies suggest conservation of mechanistic aspects of desensitization among LGICs of prokaryotic and eukaryotic origin.     

Conformational changes associated with L16P and T118M mutations in the membrane-embedded PMP22 protein,consequential in CMT-1A

Peripheral myelin protein 22 (PMP22) resides in the plasma membrane and is required for myelin formation in the peripheral nervous system. Excess PMP22 mutants accumulate in the endoplasmic reticulum (ER) resulting in the inherited neuropathies of Charcot–Marie–Tooth disease. However, there was no evidence of the structure of PMP22 or how mutations affect its folding. Therefore, in this study, we combined bioinformatics and homology modeling approaches to obtain three-dimensional native and mutated PMP22 models and its anchoring to a POPC membrane, submitted to .5-μs MD simulations, to determine how the L16P and T118M mutations affect the conformational behavior of PMP22. In addition, we investigated the ability of the native and mutated species to accumulate in the ER, via interaction with RER1, by combining protein–protein docking and MD simulations, taking the conformations that were most representative of the native and mutated PMP22 systems and RER1 conformations. Principal component analysis over MD simulations revealed that L16P and T118M mutations resulted in increased structural instability compared to the native form, which is consistent with previous experimental findings of increased structural fluctuations along a loop connecting transmembrane α-helix1 and α-helix2. Docking and MD simulations coupled with the MMGBSA approach allowed the identification that the binding interface for the PMP22-RER1 complex takes place through transmembrane α-helix1 and α-helix2, with higher effective binding free energy values between the mutated PMP22 systems and RER1 than for the native PMP22, mainly through van der Waals interactions.     

Kinetic steps for α-helix formation

The kinetics of α-helix formation in polyalanine and polyglycine eicosamers (20-mers) were examined using torsional-coordinate molecular dynamics (MD). Of one hundred fifty-five MD experiments on extended (Ala)₂₀ carried out for 0.5 ns each, 129 (83%) formed a persistent α-helix. In contrast, the extended state of (Gly)₂₀ only formed a right-handed α-helix in two of the 20 MD experiments (10%), and these helices were not as long or as persistent as those of polyalanine. These simulations show helix formation to be a competition between the rates of (a) forming local hydrogen bonds (i.e. hydrogen bonds between any residue i and its i + 2, i + 3, i + 4, or i + 5th neighbor) and (b) forming nonlocal hydrogen bonds (HBs) between residues widely separated in sequence. Local HBs grow rapidly into an α-helix; but nonlocal HBs usually retard helix formation by “trapping” the polymer in irregular, “balled-up” structures. Most trajectories formed some nonlocal HBs, sometimes as many as eight. But, for (Ala)₂₀, most of these eventually rearranged to form local HBs that lead to α-helices. A simple kinetic model describes the rate of converting nonlocal HBs into α-helices. Torsional-coordinate MD speeds folding by eliminating bond and angle degrees of freedom and reducing dynamical friction. Thus, the observed 210 ps half-life for helix formation is likely to be a lower bound on the real rate. However, we believe the sequential steps observed here mirror those of real systems. Proteins 33:343–357, 1998. © 1998 Wiley-Liss, Inc.     

Deuterium isotope effects during carbon-hydrogen bond cleavage by trimethylamine dehydrogenase. Implications for mechanism and vibrationally assisted hydrogen tunneling in wild-type and mutant enzymes

His-172 and Tyr-169 are components of a triad in the active site of trimethylamine dehydrogenase (TMADH) comprising Asp-267, His-172, and Tyr-169. Stopped-flow kinetic studies with trimethylamine as substrate have indicated that mutation of His-172 to Gln reduces the limiting rate constant for flavin reduction approximately 10-fold (Basran, J., Sutcliffe, M. J., Hille, R., and Scrutton, N. S. (1999) Biochem. J. 341, 307-314). A kinetic isotope effect (KIE = k(H)/k(D)) accompanies flavin reduction by H172Q TMADH, the magnitude of which varies significantly with solution pH. With trimethylamine, flavin reduction by H172Q TMADH is controlled by a single macroscopic ionization (pK(a) = 6.8 +/- 0.1). This ionization is perturbed (pK(a) = 7.4 +/- 0.1) in reactions with perdeuterated trimethylamine and is responsible for the apparent variation in the KIE with solution pH. At pH 9.5, where the functional group controlling flavin reduction is fully ionized, the KIE is independent of temperature in the range 277-297 K, consistent with vibrationally assisted hydrogen tunneling during breakage of the substrate C-H bond. Y169F TMADH is approximately 4-fold more compromised than H172Q TMADH for hydrogen transfer, which occurs non-classically. Studies with Y169F TMADH suggest partial thermal excitation of substrate prior to hydrogen tunneling by a vibrationally assisted mechanism. Our studies illustrate the varied effects of compromising mutations on tunneling regimes in enzyme molecules.     

Structure of the M2 transmembrane segment of GLIC, a prokaryotic Cys loop receptor homologue from Gloeobacter violaceus, probed by substituted cysteine accessibility

GLIC is a homopentameric proton-gated, prokaryotic homologue of the Cys loop receptor family of neurotransmitter-gated ion channels. Recently, crystal structures of GLIC hypothesized to represent an open channel state were published. To explore the channel structure in functional GLIC channels, we tested the ability of p-chloromercuribenzenesulfonate to react with 30 individual cysteine substitution mutants in and flanking the M2 channel-lining segment in the closed state (pH 7.5) and in a submaximally activated state (pH 5.0). Nine mutants did not tolerate cysteine substitution and were not functional. From positions 10' to 27', p-chloromercuribenzenesulfonate significantly modified the currents at pH 7.5 and 5.0 in all mutants except H234C (11'), I235C (12'), V241C (18'), T243C (20'), L245C (22'), and Y250C (27'), which were not functional, except for 12'. Currents for P246C (23') and K247C (24') were only significantly altered at pH 5.0. The reaction rates were all >1000 m(-1) s(-1). The reactive residues were more accessible in the activated than in the resting state. We infer that M2 is tightly associated with the adjacent transmembrane helices at the intracellular end but is more loosely packed from 10' to the extracellular end than the x-ray structures suggest. We infer that the charge selectivity filter is in the cytoplasmic half of the channel. We also show that below pH 5.0, GLIC desensitizes on a time scale of minutes and infer that the crystal structures may represent a desensitized state.     

Roles of the histidine and tryptophan side chains in the M2 proton channel from influenza A virus

The M2 protein form influenza A virus forms a tetrameric ion channel, which enables proton passage across biological membranes when the N-terminal side is acidified. Among the amino acid residues in the transmembrane domain of the M2 protein, His37 and Trp41 are essential for the pH-regulated proton conductance. Current knowledge about the structures and interactions of His37 and Trp41 suggests a model for the M2 ion channel, in which the channel is closed by a network of His37 hydrogen bonds at neutral pH and is opened by a His37-Trp41 cation-pi interaction at acidic pH.     

Mechanistic characterization of N-formimino-L-glutamate iminohydrolase from Pseudomonas aeruginosa

N-Formimino-l-glutamate iminohydrolase (HutF) from Pseudomonas aeruginosa catalyzes the deimination of N-formimino-l-glutamate in the histidine degradation pathway. An amino acid sequence alignment between HutF and members of the amidohydrolase superfamily containing mononuclear metal centers indicated that residues Glu-235, His-269, and Asp-320 are involved in substrate binding and activation of the nucleophilic water molecule. The purified enzyme contained up to one equivalent of zinc. The metal was removed by dialysis against the metal chelator dipicolinate with the complete loss of catalytic activity. Enzymatic activity was restored by incubation of the apoprotein with Zn2+, Cd2+, Ni2+, or Cu2+. The mutation of Glu-235, His-269, or Asp-320 resulted in the diminution of catalytic activity by two to six orders of magnitude. Bell-shaped profiles were observed for kcat and kcat/Km as a function of pH. The pKa of the group that must be unprotonated for catalytic activity was consistent with the ionization of His-269. This residue is proposed to function as a general base in the abstraction of a proton from the metal-bound water molecule. In the proposed catalytic mechanism, the reaction is initiated by the abstraction of a proton from the metal-bound water molecule by the side chain imidazole of His-269 to generate a tetrahedral intermediate of the substrate. The collapse of the tetrahedral intermediate commences with the abstraction of a second proton via the side chain carboxylate of Asp-320. The C-N bond of the substrate is subsequently cleaved with proton transfer from His-269 to form ammonia and the N-formyl product. The postulated role of the invariant Glu-235 is to ion pair with the positively charged formimino group of the substrate.     

Mechanism of proton gating of a urea channel

The size and complexity of many pH-gated channels have frustrated the development of specific structural models. The small acid-activated six-membrane segment urea channel of Helicobacter hepaticus (HhUreI), homologous to the essential UreI of the gastric pathogen Helicobacter pylori, enables identification of all the periplasmic sites of proton gating by site-directed mutagenesis. Exposure to external acidity enhances [(14)C]urea uptake by Xenopus oocytes expressing HhUreI, with half-maximal activity (pH(0.5)) at pH 6.8. A downward shift of pH(0.5) in single site mutants identified four of six protonatable periplasmic residues (His-50 at the boundary of the second transmembrane segment TM2, Glu-56 in the first periplasmic loop, Asp-59 at the boundary of TM3, and His-170 at the boundary of TM6) that affect proton gating. Asp-59 was the only site at which a protonatable residue appeared to be essential for pH gating. Mutation of Glu-110 or Glu-114 in PL2 did not affect the pH(0.5) of gating. A chimera, where the entire periplasmic domain of HhUreI was fused to the membrane domain of Streptococcus salivarius UreI (SsUreI), retained the pH-independent properties of SsUreI. Hence, proton gating of HhUreI likely depends upon the formation of hydrogen bonds by periplasmic residues that in turn produce conformational changes of the transmembrane domain. Further studies on HhUreI may facilitate understanding of other physiologically important pH-responsive channels.     

Stabilization of the imidazole ring of His-195 at the active site of chloramphenicol acetyltransferase

The imidazole of His-195 plays an essential role in the proposed general base mechanism of chloramphenicol acetyltransferase (CAT). The structure of the binary complex of CATIII and chloramphenicol suggests that two unusual interactions might determine the conformation of the side chain of His-195: (i) an intraresidue hydrogen bond between its main chain carbonyl and the protonated N delta 1 of the imidazole ring and (ii) face-to-face van der Waals contact between the His-195 imidazole group and the aromatic side chain of Tyr-25. Tyr-25 also makes a hydrogen bond, via its phenolic hydroxyl, to the carbonyl oxygen of the substrate chloramphenicol. Replacement of Tyr-25 of CATIII by phenylalanine results in a modest increase in the Km for chloramphenicol (from 11.6 to 14.6 microM) and a 2-fold fall in kcat (599 to 258 s-1), indicative of a free energy contribution to transition state binding of 0.6 kcal mol-1 for the hydrogen bond between Tyr-25 and chloramphenicol. In contrast, substitution of Tyr-25 by alanine yields an enzyme that is dramatically impaired in its ability to bind chloramphenicol (Km = 173 microM). As kcat for Ala-25 CAT is also reduced (130 s-1), the loss of the aryl group results in a 69-fold decrease in kcat/Km, corresponding to a free energy contribution to binding and catalysis of 2.5 kcal mol-1. In addition to the loss of the hydrogen bond between Tyr-25 and chloramphenicol, the loss of substrate affinity in Ala-25 CAT may be a direct consequence of reduced hydrophobicity of the chloramphenicol-binding site and/or the loss of critical constraints on the precise conformation of the catalytic imidazole. However, as with wild type CAT, inactivation of Ala-25 CAT by the affinity reagent 3-(bromoacetyl) chloramphenicol is accompanied by modification solely at N epsilon 2 of His-195. Hence, the results demonstrate that tautomeric stabilization of the imidazole ring persists in the absence of van der Waals interactions with the side chain of Tyr-25, probably as a consequence of hydrogen bonding between the protonated N delta 1 and the carbonyl oxygen of His-195.     

Water as a structural element in a channel: gating in the Kcsa channel, and implications for voltage-gated ion channels

Water is becoming understood as a structural element in proteins. Here we are concerned with one particular type of protein, ion channels. The S. Lividans KcsA K(+) channel, the X-ray structure of which is known, is gated by protons (i.e, by a drop in pH). Ab initio calculations suggest that an H(5)O(2) group, partially charged, connects the E118 residues in the gating region, when the four residues have a -2 net charge, but that the hydrogen bonding is not strong enough to do this when the charge becomes -1. The H(5)O(2) group would block the channel, in the -2 state, and prevent motion of the four transmembrane (TM) segments of the protein, by binding them. With the weaker bond in the -1 state, the TM segments would be able to separate (as they have been found to do experimentally, opening the channel. Voltage gated channels have four additional TM segments for each of the four domains of the channel protein. These appear to allow motion of protons; in fact there is evidence that the initial step in gating must be the transfer of a proton. We have earlier shown that the transfer of a single proton between two methylamines under the influence of a field is possible, as proton tunneling. Subsequent steps are hypothesized to result from four proton transfer cascades of about three protons each, triggered by the initial proton transfer. We suggest that the extra 4 TM segments of the voltage gated channel act as a voltage to proton-current transducer. Water, held by hydrogen bonds, is also suggested as the source of the accessibility data found with MTS reagents, based largely on simulations, our earlier Monte Carlo simulations as well as molecular dynamics studies reported by others. These waters may also play a structural role in the protein.     

Epidermal growth factor from the mouse. Structural characterization by proton nuclear magnetic resonance and nuclear overhauser experiments at 500 MHz

Mouse epidermal growth factor (mEGF), a protein hormone effector molecule that regulates cellular development and division, has been investigated by using proton nuclear magnetic resonance techniques at 500 MHz. Well-resolved downfield aromatic and alpha-CH proton resonances (5.0-8.0 ppm) and an upfield ring current shifted isoleucine delta-methyl resonance (approximately 0.5 ppm) have been examined by using principally nuclear Overhauser methods. The data are analyzed in terms of model building based on the predictive Chou-Fasman secondary structure algorithm applied to mEGF [Holladay, L. A., Savage, C. R., Cohen, S., & Puett, D. (1976) Biochemistry 15, 2624-2633], which suggests the existence of some beta-structure and little or no alpha-helicity. Proximity relationships derived from nuclear Overhauser data among Tyr-3, -10, and -13, His-22, and Ile-23 allow refinement of some aspects of the predicted secondary structure and render additional information on how the protein backbone in mEGF is folded (i.e., tertiary structure). Nuclear Overhauser effects (NOEs) from irradiation of several alpha-CH proton resonances give evidence for tiered beta-sheet structure in mEGF. Such proximity relationships derived from NOE data place stringent limitations on possible models for the molecule. pH titration data demonstrate a His-22 pKa of 7.1, indicating either a salt bridge or hydrogen-bond formation between His-22 and another residue. The His-22 pKa is also reflected in the chemical shift changes of several other resonances as a function of pH. Nuclear Overhauser methods, used to differentiate direct (protonation) and indirect (conformation) effects on the chemical shift changes in the spectra of mEGF by varying the pH, yield evidence for a pH-induced conformational transition in the protein hormone associated with the breaking of the His-22 salt bridge or hydrogen bond.     

Independent Contribution of Extracellular Proton Binding Sites to ASIC1a Activation

Acid-sensing ion channels (ASICs) are a group of trimeric cation permeable channels gated by extracellular protons that are mainly expressed in the nervous system. Despite the structural information available for ASIC1, there is limited understanding of the molecular mechanism that allows these channels to sense and respond to drops in extracellular pH. In this report, we employed the substituted cysteine accessibility method and site-directed mutagenesis to examine the mechanism of activation of ASIC1a by extracellular protons. We found that the modification of E238C and D345C channels by MTSET reduced proton apparent affinity for activation. Furthermore, the introduction of positively charged residues at position 345 rendered shifted biphasic proton activation curves. Likewise, channels bearing mutations at positions 79 and 416 in the palm domain of the channel showed reduced proton apparent affinity and biphasic proton activation curves. Of significance, the effect of the mutations at positions 79 and 345 on channel activation was additive. E79K-D345K required a change to a pH lower than 2 for maximal activation. In summary, this study provides direct evidence for the presence of two distinct proton coordination sites in the extracellular region of ASIC1a, which jointly facilitate pore opening in response to extracellular acidification.     

Anesthetic Binding in a Pentameric Ligand-Gated Ion Channel: GLIC

Cys-loop receptors are molecular targets of general anesthetics, but the knowledge of anesthetic binding to these proteins remains limited. Here we investigate anesthetic binding to the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), a structural homolog of cys-loop receptors, using an experimental and computational hybrid approach. Tryptophan fluorescence quenching experiments showed halothane and thiopental binding at three tryptophan-associated sites in the extracellular (EC) domain, transmembrane (TM) domain, and EC-TM interface of GLIC. An additional binding site at the EC-TM interface was predicted by docking analysis and validated by quenching experiments on the N200W GLIC mutant. The binding affinities (K_D) of 2.3 ± 0.1 mM and 0.10 ± 0.01 mM were derived from the fluorescence quenching data of halothane and thiopental, respectively. Docking these anesthetics to the original GLIC crystal structure and the structures relaxed by molecular dynamics simulations revealed intrasubunit sites for most halothane binding and intersubunit sites for thiopental binding. Tryptophans were within reach of both intra- and intersubunit binding sites. Multiple molecular dynamics simulations on GLIC in the presence of halothane at different sites suggested that anesthetic binding at the EC-TM interface disrupted the critical interactions for channel gating, altered motion of the TM23 linker, and destabilized the open-channel conformation that can lead to inhibition of GLIC channel current. The study has not only provided insights into anesthetic binding in GLIC, but also demonstrated a successful fusion of experiments and computations for understanding anesthetic actions in complex proteins.     

Distinct histidine residues control the acid-induced activation and inhibition of the cloned K(ATP) channel

The modulation of K(ATP) channels during acidosis has an impact on vascular tone, myocardial rhythmicity, insulin secretion, and neuronal excitability. Our previous studies have shown that the cloned Kir6.2 is activated with mild acidification but inhibited with high acidity. The activation relies on His-175, whereas the molecular basis for the inhibition remains unclear. To elucidate whether the His-175 is indeed the protonation site and what other structures are responsible for the pH-induced inhibition, we performed these studies. Our data showed that the His-175 is the only proton sensor whose protonation is required for the channel activation by acidic pH. In contrast, the channel inhibition at extremely low pH depended on several other histidine residues including His-186, His-193, and His-216. Thus, proton has both stimulatory and inhibitory effects on the Kir6.2 channels, which attribute to two sets of histidine residues in the C terminus.     

Leukocyte integrin αLβ2 transmembrane association dynamics revealed by coarse-grained molecular dynamics simulations

Integrins are transmembrane (TM) proteins that mediate bidirectional mechanical signaling between the extracellular matrix and the cellular cytoskeletal network. Each integrin molecule consists of non-covalently associated α- and β-subunits, with each subunit consisting of a large ectodomain, a single-pass TM helix, and a short cytoplasmic tail. Previously we found evidence for a polar interaction (hydrogen bond) in the outer membrane clasp (OMC) of the leukocyte integrin αLβ2 TMs that is absent in the platelet integrin αIIβ3 OMC. Here, we compare the self-assembly dynamics of αLβ2 and αIIβ3 TM helices in a model membrane using coarse-grained molecular dynamics simulations. We found that although αIIβ3 TM helices associate more easily, packing is suboptimal. In contrast, αLβ2 TM helices achieve close-to-optimal packing. This suggests that αLβ2 TM packing is more specific, possibly due to the interhelix hydrogen bond. Theoretical association free energy profiles show a deeper minimum at a smaller helix-helix separation for αLβ2 compared with αIIβ3. The αIIβ3 profile is also more rugged with energetic barriers whereas that of αLβ2 is almost without barriers. Disruption of the interhelix hydrogen bond in αLβ2 via the β2T686G mutation results in poorer association and a similar profile as αIIβ3. The OMC polar interaction in αLβ2 thus plays a significant role in the packing of the TM helices.