Arg Side Chain Research Articles

The Role of Active-Site Residues Phe98, His239, and Arg243 in DNA Binding and in the Catalysis of Human Uracil-DNA Glycosylase SMUG1.

Human SMUG1 (hSMUG1) hydrolyzes the N-glycosidic bond of uracil and some uracil lesions formed in the course of epigenetic regulation. Despite the functional importance of hSMUG1 in the DNA repair pathway, the damage recognition mechanism has been elusive to date. In the present study, our objective was to build a model structure of the enzyme–DNA complex of wild-type hSMUG1 and several hSMUG1 mutants containing substitution F98W, H239A, or R243A. Enzymatic activity of these mutant enzymes was examined by polyacrylamide gel electrophoresis analysis of the reaction product formation and pre-steady-state analysis of DNA conformational changes during enzyme–DNA complex formation. It was shown that substitutions F98W and H239A disrupt specific contacts generated by the respective wild-type residues, namely stacking with a flipped out Ura base in the damaged base-binding pocket or electrostatic interactions with DNA in cases of Phe98 and His239, respectively. A loss of the Arg side chain in the case of R243A reduced the rate of DNA bending and increased the enzyme turnover rate, indicating facilitation of the product release step.

Adsorption of RGD Tripeptide on Au (111) Surface

To date, many studies have been performed on interactions between amino acids and inorganic surfaces. However, the interactions are still in debate. Au is one of important biometals with inertness and subsequent biocompatibility in a living body. Because amino acids tend to adsorb on Au surface, it is required to investigate interactions between amino acids and Au surface for safe application of Au. In the present work, adsorption of RGD (arginine-glycine-aspartic acid, Arg-Gly-Asp) sequence on Au (111) surface is investigated by first principles calculations. When two oxygen atoms of carboxyl group in Asp were located very close to the atop sites of the Au surface, RGD was the most stable on the Au surface, like the adsorption of a single molecule of Asp on the Au surface. The adsorption energy of RGD on the Au surface was lower than that of Asp on the Au surface. This is because the negative hyperpolarization of Asp side chain and the positive hyperpolarization of Arg side chain are enhanced by the RGD sequence.

Energetic switch of the proline effect in collision‐induced dissociation of singly and doubly protonated peptide Ala‐Ala‐Arg‐Pro‐Ala‐Ala

Suppression of the selective cleavage at N-terminal of proline is observed in the peptide cleavage by proteolytic enzyme trypsin and in the fragment ion mass spectra of peptides containing Arg-Pro sequence. An insight into the fragmentation mechanism of the influence of arginine residue on the proline effect can help in prediction of mass spectra and in protein structure analysis. In this work, collision-induced dissociation spectra of singly and doubly charged peptide AARPAA were studied by ESI MS/MS and theoretical calculation methods. The proline effect was evaluated by comparing the experimental ratio of fragments originated from cleavage of different amide bonds. The results revealed that the backbone amide bond cleavage was selected by the energy barrier height of the fragmentation pathway although the strong proton affinity of the Arg side chain affected the stereostructure of the peptide and the dissociation mechanism. The thermodynamic stability of the fragment ions played a secondary role in the abundance ratio of fragments generated via different pathways. Fragmentation studies of protonated peptide AACitPAA supported the energy-dependent hypothesis. The results provide an explanation to the long-term arguments between the steric conflict and the proton mobility mechanisms of proline effect.

Recent Updates on Protein N-Phosphoramidate Hydrolases.

In contrast to well-recognized protein phosphorylation on the side-chain oxygen of Ser, Thr, or Tyr residues, analogous phosphoramidation of the nitrogen of His, Lys, and Arg side chains remains much less investigated, mainly due to the instability of post-translational modifications and technical difficulties involved in their analysis. For example, reports on the enzyme activities responsible for the formation and hydrolysis of these phosphoramidates date back to as early as the 1950s, but some of these enzymes have only recently been identified and functionally characterized; this has been aided by the development of novel research tools. In this review, we summarize current knowledge of the enzymes that hydrolyze protein N-phosphoramidates, in terms of their structure, activities, and biological functions, as well as the chemical tools used to investigate them.

NMR Methods for Characterizing the Basic Side Chains of Proteins: Electrostatic Interactions, Hydrogen Bonds, and Conformational Dynamics.

NMR spectroscopy is a powerful tool for studying protein dynamics. Conventionally, NMR studies on protein dynamics have probed motions of protein backbone NH, side-chain aromatic, and CH3 groups. Recently, there has been remarkable progress in NMR methodologies that can characterize motions of cationic groups in protein side chains. These NMR methods allow investigations of the dynamics of positively charged lysine (Lys) and arginine (Arg) side chains and their hydrogen bonds as well as their electrostatic interactions important for protein function. Here, describing various practical aspects, we provide an overview of the NMR methods for dynamics studies of Lys and Arg side chains. Some example data on protein-DNA complexes are shown. We will also explain how molecular dynamics (MD) simulations can facilitate the interpretation of the NMR data on these basic side chains. Studies combining NMR and MD have revealed the highly dynamic nature of short-range electrostatic interactions via ion pairs, especially those involving Lys side chains.

Phosphate-Catalyzed Succinimide Formation from an NGR-Containing Cyclic Peptide: A Novel Mechanism for Deammoniation of the Tetrahedral Intermediate.

Spontaneous deamidation in the Asn-Gly-Arg (NGR) motif that yields an isoAsp-Gly-Arg (isoDGR) sequence has recently attracted considerable attention because of the possibility of application to dual tumor targeting. It is well known that Asn deamidation reactions in peptide chains occur via the five-membered ring succinimide intermediate. Recently, we computationally showed by the B3LYP density functional theory method, that inorganic phosphate and the Arg side chain can catalyze the NGR deamidation using a cyclic peptide, c[CH2CO–NGRC]–NH2. In this previous study, the tetrahedral intermediate of the succinimide formation was assumed to be readily protonated at the nitrogen originating from the Asn side chain by the solvent water before the release of an NH3 molecule. In the present study, we found a new mechanism for the decomposition of the tetrahedral intermediate that does not require the protonation by an external proton source. The computational method is the same as in the previous study. In the new mechanism, the release of an NH3 molecule occurs after a proton exchange between the peptide and the phosphate and conformational changes. The rate-determining step of the overall reaction course is the previously reported first step, i.e., the cyclization to form the tetrahedral intermediate.

The Role of the Arginine in the Conserved N-Terminal Domain RLFDQxFG Motif of Human Small Heat Shock Proteins HspB1, HspB4, HspB5, HspB6, and HspB8.

Although the N-terminal domain of vertebrate small heat shock proteins (sHsp) is poorly conserved, it contains a core motif preserved in many members of the sHsp family. The role of this RLFDQxFG motif remains elusive. We analyzed the specific role of the first arginine residue of this conserved octet sequence in five human sHsps (HspB1, HspB4, HspB5, HspB6, and HspB8). Substitution of this arginine with an alanine induced changes in thermal stability and/or intrinsic fluorescence of the related HspB1 and HspB8, but yielded only modest changes in the same biophysical properties of HspB4, HspB5, and HspB6 which together belong to another clade of vertebrate sHsps. Removal of the positively charged Arg side chain resulted in destabilization of the large oligomers of HspB1 and formation of smaller size oligomers of HspB5. The mutation induced only minor changes in the structure of HspB4 and HspB6. In contrast, the mutation in HspB8 was accompanied by shifting the equilibrium from dimers towards the formation of larger oligomers. We conclude that the RLFDQxFG motif plays distinct roles in the structure of several sHsp orthologs. This role correlates with the evolutionary relationship of the respective sHsps, but ultimately, it reflects the sequence context of this motif.

Enhancing Oral Bioavailability of Cyclic RGD Hexa-peptides by the Lipophilic Prodrug Charge Masking Approach: Redirection of Peptide Intestinal Permeability from a Paracellular to Transcellular Pathway.

Hydrophilic peptides constitute most of the active peptides. They mostly permeate via tight junctions (paracellular pathway) in the intestine. This permeability mechanism restricts the magnitude of their oral absorption and bioavailability. We hypothesized that concealing the hydrophilic residues of the peptide using the lipophilic prodrug charge masking approach (LPCM) can improve the bioavailability of hydrophilic peptides. To test this hypothesis, a cyclic N-methylated hexapeptide containing Arg-Gly-Asp (RGD) and its prodrug derivatives, masking the Arg and Asp charged side chains, were synthesized. The library was evaluated for intestinal permeability in vitro using the Caco-2 model. Further investigation of metabolic stability ex vivo models in rat plasma, brush border membrane vesicles (BBMVs), and isolated CYP3A4 microsomes and pharmacokinetic studies was performed on a selected peptide and its prodrug (peptide 12). The parent drug analogues were found to have a low permeability rate in vitro, corresponding to atenolol, a marker for paracellular permeability. Moreover, palmitoyl carnitine increased the Papp of peptide 12 by 4-fold, indicating paracellular permeability. The Papp of the prodrug derivatives was much higher than that of their parent peptides. For instance, the Papp of the prodrug 12P was 20-fold higher than the Papp of peptide 12 in the apical to basolateral (AB) direction. Whereas the permeability in the opposite direction (BA of the Caco-2 model) was significantly faster than the Papp AB, indicating the involvement of an efflux system. These results were corroborated when verapamil, a P-gp inhibitor, was added to the Caco-2 model and increased the Papp AB of prodrug 12P by 3-fold. The prodrug 12P was stable in the BBMVs environment, yet degraded quickly (less than 5 min) in the plasma into the parent peptide 12. Pharmacokinetic studies in rats showed an increase in the bioavailability of peptide 12 > 70-fold (from 0.58 ± 0.11% to 43.8 ± 14.9%) after applying the LPCM method to peptide 12 and converting it to the prodrug 12P. To conclude, the LPCM approach converted the absorption mechanism of the polar peptides from a paracellular to transcellular pathway that tremendously affects their oral bioavailability. The LPCM method provides a solution for the poor bioavailability of RGD cyclohexapeptides and paves the way for other active hydrophilic and charged peptides with poor oral bioavailability.

Impact of two-bond 15N-15N scalar couplings on 15N transverse relaxation measurements for arginine side chains of proteins.

NMR relaxation of arginine (Arg) 15Nε nuclei is useful for studying side-chain dynamics of proteins. In this work, we studied the impact of two geminal 15N-15N scalar couplings on measurements of transverse relaxation rates (R 2 ) for Arg side-chain 15Nε nuclei. For 12 Arg side chains of the DNA-binding domain of the Antp protein, we measured the geminal 15N-15N couplings ( 2 J NN ) of the 15Nε nuclei and found that the magnitudes of the 2 J NN coupling constants were virtually uniform with an average of 1.2Hz. Our simulations, assuming ideal 180° rotations for all 15N nuclei, suggested that the two 2 J NN couplings of this magnitude could in principle cause significant modulation in signal intensities during the Carr-Purcell-Meiboom-Gill (CPMG) scheme for Arg 15Nε R 2 measurements. However, our experimental data show that the expected modulation via two 2 J NN couplings vanishes during the 15N CPMG scheme. This quenching of J modulation can be explained by the mechanism described in Dittmer and Bodenhausen (Chemphyschem 7:831-836, 2006). This effect allows for accurate measurements of R 2 relaxation rates for Arg side-chain 15Nε nuclei despite the presence of two 2 J NN couplings. Although the so-called recoupling conditions may cause overestimate of R 2 rates for very mobile Arg side chains, such conditions can readily be avoided through appropriate experimental settings.

Zn2+ Modulates Hv1 Proton Channel Gating via Conformational Coupling to an Intracellular Coulombic Network

Gating of the intrinsic ‘aqueous’ H+ conductance (GAQ) in the Hv1 proton channel is sensitive to extracellular Zn2+ and other divalent cations in the micromolar range. Previous studies showed that the extracellularly-accessible residues H140 and H193 are required for Zn2+ to positively shift the GAQ-V relation. However, the molecular mechanism by which extracellular Zn2+ occupancy is coupled to voltage-dependent GAQ gating remains known. We used site-directed mutagenesis and electrophysiology to identify Zn2+-interacting residues and used the experimental data to build new models of the Hv1 voltage sensor (VS) domain in Zn2+-liganded (Hv1F·Zn2+) and unliganded (Hv1F + Zn2+) conformations. Although imidazole nitrogen atoms (H140-Nδ1 and H193-Nδ1) are close enough to simultaneously coordinate a single Zn2+ ion, we did not observe a stable Zn2+ coordination sphere during all-atom molecular dynamics (MD) simulations using the CHARMM36 force field. To simulate electronic interactions between Zn2+ and H140-Nδ1 and H193-Nδ1, we imposed harmonic constraints (5 kcal/mol/Å) and equilibrated the Hv1F·Zn2+ model during a 50 ns MD simulation. 20 ns after releasing constraints, ZN13 spontaneously dissociates from H140-Nδ1 and H193-Nδ1 but remains associated with the carboxyl group of E119. The VS also undergoes backbone reorganization involving the S1, S3 and S4 helices after ZN13 dissociation that suggest Zn2+-dependent conformational changes. Concomitant changes in the organization of an electrostatic network involving intracellular acidic residues (E153 and D174) and S4 ‘gating charge’ Arg side chains may reveal how voltage-dependent gating is altered by Zn2+. In support of this hypothesis, we find that E153 neutralization also alters Zn2+ potency. In summary, extracellular Zn2+ occupancy is coupled to reorganization of an intracellular Coulombic network that is required for voltage-dependent gating in Hv1.

Modeling of a C-end rule peptide adsorbed onto gold nanoparticles.

The RPAR peptide, a prototype C-end Rule (CendR) sequence that binds to neuropilin-1 (NRP-1), has potential therapeutic uses as internalization trigger in anticancer nanodevices. Recently, the functionalization of gold nanoparticles with CendR peptides has been proved to be a successful strategy to target the NRP-1 receptor in prostate cancer cells. In this work, we investigate the influence of two gold surface facets, (100) and (111), on the conformational preferences of RPAR using molecular dynamics simulations. Both clustering and conformational analyses revealed that the peptide backbone becomes very rigid upon adsorption onto gold, which is a very fast and favored process, the only flexibility being attributed to the side chains of the two Arg residues. Thus, the different components of RPAR tend to adopt an elongated shape, which is characterized by the pseudo-extended conformation of both the backbone and the Arg side chains. This conformation is very different from the already known bioactive conformation, indicating that RPAR is drastically affected by the substrate. Interestingly, the preferred conformations of the peptide adsorbed onto gold facets are not stabilized by salt bridges and/or specific intramolecular hydrogen bonds, which represent an important difference with respect to the conformations found in other environments (e.g. the peptide in solution and interacting with NRP-1 receptor). However, the conformational changes induced by the substrate are not detrimental for the use of gold nanoparticles as appropriate vehicles for the transport and targeted delivery of the RPAR. Thus, once their high affinity for the NRP-1 receptor induces the targeted delivery of the elongated peptide molecules from the gold nanoparticles, the lack of intramolecular interactions facilitates their evolution towards the bioactive conformation, increasing the therapeutic efficacy of the peptide.

Internal Motions of Basic Side Chains of the Antennapedia Homeodomain in the Free and DNA-Bound States.

Basic side chains play crucial roles in protein-DNA interactions. In this study, using NMR spectroscopy, we investigated the dynamics of Arg and Lys side chains of the fruit fly Antennapedia homeodomain in the free state and in the complex with target DNA. We measured 15N relaxation for Arg and Lys side chains at two magnetic fields, from which generalized order parameters for the cationic groups were determined. Mobility of the R5 side chain, which makes hydrogen bonds with a thymine base in the DNA minor groove, was greatly dampened. Several Lys and Arg side chains that form intermolecular ion pairs with DNA phosphates were found to retain high mobility with the order parameter being <0.6 in the DNA-bound state. Interestingly, some of the interfacial cationic groups in the complex were more mobile than in the free protein. The retained or enhanced mobility of the Arg and Lys side chains in the complex should mitigate the overall loss of conformational entropy in the protein-DNA association and allow dynamic molecular recognition.

NMR Insights into the Structure-Function Relationships in the Binding of Melanocortin Analogues to the MC1R Receptor.

Linear and cyclic analogues of the α-melanocyte stimulating hormone (α-MSH) targeting the human melanocortin receptor 1 (MC1R) are of pharmacological interest for detecting and treating melanoma. The central sequence of α-MSH (His-Phe-Arg-Trp) has been identified as being essential for receptor binding. To deepen current knowledge on the molecular basis for α-MSH bioactivity, we aimed to understand the effect of cycle size on receptor binding. To that end, we synthesised two macrocyclic isomeric α-MSH analogues, c[NH-NO2-C6H3-CO-His-DPhe-Arg-Trp-Lys]-Lys-NH2 (CycN-K6) and c[NH-NO2-C6H3-CO-His-DPhe-Arg-Trp-Lys-Lys]-NH2 (CycN-K7). Their affinities to MC1R receptor were determined by competitive binding assays, and their structures were analysed by 1H and 13C NMR. These results were compared to those of the previously reported analogue c[S-NO2-C6H3-CO-His-DPhe-Arg-Trp-Cys]-Lys-NH2 (CycS-C6). The MC1R binding affinity of the 22-membered macrocyclic peptide CycN-K6 (IC50 = 155 ± 16 nM) is higher than that found for the 25-membered macrocyclic analogue CycN-K7 (IC50 = 495 ± 101 nM), which, in turn, is higher than that observed for the 19-membered cyclic analogue CycS-C6 (IC50 = 1770 ± 480 nM). NMR structural study indicated that macrocycle size leads to changes in the relative dispositions of the side chains, particularly in the packing of the Arg side chain relative to the aromatic rings. In contrast to the other analogues, the 22-membered cycle’s side chains are favorably positioned for receptor interaction.

Chemical and biological properties of a supramolecular complex of tuftsin and cucurbit[7]uril

Cucurbit[7]uril (CB7) is an uncharged and water-soluble macrocyclic host. CB7 binds to doubly protonated tuftsin, which is the tetrapeptide Thr-Lys-Pro-Arg, with moderate affinity (Ka=2.1×103M−1). In this study, the host–guest complexation was confirmed by fluorescence titration. This affinity would allow for easy release of the peptide under physiological conditions. According to density functional theory calculations, the structural binding motif involves hydrogen bonding. The most energetically stable form had the Arg side chain inside the CB7 cavity. The effects of the tuftsin–CB7 complex on the proliferation and cytokine activity of immune cells were studied. The complex had broader spectrum immunomodulation than free peptides, and caused statistically significant (p<0,05) changes in cytokine production (tumor necrosis factor-α, interleukin-2, interferon-γ, and interleukin-10) by mononuclear cells. By contrast, the free peptide only activated tumor necrosis factor-α production.

Unique Contributions of an Arginine Side Chain to Ligand Recognition in a Glutamate-gated Chloride Channel

Glutamate recognition by neurotransmitter receptors often relies on Arg residues in the binding site, leading to the assumption that charge-charge interactions underlie ligand recognition. However, assessing the precise chemical contribution of Arg side chains to protein function and pharmacology has proven to be exceedingly difficult in such large and complex proteins. Using the in vivo nonsense suppression approach, we report the first successful incorporation of the isosteric, titratable Arg analog, canavanine, into a neurotransmitter receptor in a living cell, utilizing a glutamate-gated chloride channel from the nematode Haemonchus contortus. Our data unveil a surprisingly small contribution of charge at a conserved arginine side chain previously suggested to form a salt bridge with the ligand, glutamate. Instead, our data show that Arg contributes crucially to ligand sensitivity via a hydrogen bond network, where Arg interacts both with agonist and with a conserved Thr side chain within the receptor. Together, the data provide a new explanation for the reliance of neurotransmitter receptors on Arg side chains and highlight the exceptional capacity of unnatural amino acid incorporation for increasing our understanding of ligand recognition.

Succinimide Formation from an NGR-Containing Cyclic Peptide: Computational Evidence for Catalytic Roles of Phosphate Buffer and the Arginine Side Chain.

The Asn-Gly-Arg (NGR) motif and its deamidation product isoAsp-Gly-Arg (isoDGR) have recently attracted considerable attention as tumor-targeting ligands. Because an NGR-containing peptide and the corresponding isoDGR-containing peptide target different receptors, the spontaneous NGR deamidation can be used in dual targeting strategies. It is well known that the Asn deamidation proceeds via a succinimide derivative. In the present study, we computationally investigated the mechanism of succinimide formation from a cyclic peptide, c[CH2CO-NGRC]-NH2, which has recently been shown to undergo rapid deamidation in a phosphate buffer. An H2PO4− ion was explicitly included in the calculations. We employed the density functional theory using the B3LYP functional. While geometry optimizations were performed in the gas phase, hydration Gibbs energies were calculated by the SM8 (solvation model 8) continuum model. We have found a pathway leading to the five-membered ring tetrahedral intermediate in which both the H2PO4− ion and the Arg side chain act as catalyst. This intermediate, once protonated at the NH2 group on the five-membered ring, was shown to easily undergo NH3 elimination leading to the succinimide formation. This study is the first to propose a possible catalytic role for the Arg side chain in the NGR deamidation.

Structure and Function of FS50, a salivary protein from the flea Xenopsylla cheopis that blocks the sodium channel NaV1.5.

Naturally occurring toxins have been invaluable tools for the study of structural and functional relationships of voltage-gated sodium channels (VGSC). Few studies have been made of potential channel-modulating substances from blood-feeding arthropods. He we describe the characterization FS50, a salivary protein from the flea, Xenopsylla cheopis, that exhibits an inhibitory activity against the NaV1.5 channel with an IC50 of 1.58 μM. The pore-blocking mechanism of this toxin is evident from the kinetics of activation and inactivation suggesting that FS50 does not interfere with the voltage sensor of NaV1.5. FS50 exhibits high specificity for NaV1.5, since 10 μM FS50 had no discernable effect on voltage-gated Na+, K+ and Ca2+ channels in rat dorsal root ganglia or VGSC forms individually expressed in HEK 293T cells. Furthermore, intravenous injection of FS50 into rats and monkeys elicited recovery from arrhythmia induced by BaCl2, as would be expected from a blockade of NaV1.5. The crystal structure of FS50 revealed a βαββ domain similar to that of scorpion β toxin and a small N-terminal βαβ domain. Site-directed mutagenesis experiments have implicated a basic surface including the side chains of Arg 6, His 11 and Lys 32 as potentially important in the FS50 NaV1.5 interaction.

Bypass of DNA-Protein Cross-links Conjugated to the 7-Deazaguanine Position of DNA by Translesion Synthesis Polymerases

DNA-protein cross-links (DPCs) are bulky DNA lesions that form both endogenously and following exposure to bis-electrophiles such as common antitumor agents. The structural and biological consequences of DPCs have not been fully elucidated due to the complexity of these adducts. The most common site of DPC formation in DNA following treatment with bis-electrophiles such as nitrogen mustards and cisplatin is the N7 position of guanine, but the resulting conjugates are hydrolytically labile and thus are not suitable for structural and biological studies. In this report, hydrolytically stable structural mimics of N7-guanine-conjugated DPCs were generated by reductive amination reactions between the Lys and Arg side chains of proteins/peptides and aldehyde groups linked to 7-deazaguanine residues in DNA. These model DPCs were subjected to in vitro replication in the presence of human translesion synthesis DNA polymerases. DPCs containing full-length proteins (11-28 kDa) or a 23-mer peptide blocked human polymerases η and κ. DPC conjugates to a 10-mer peptide were bypassed with nucleotide insertion efficiency 50-100-fold lower than for native G. Both human polymerase (hPol) κ and hPol η inserted the correct base (C) opposite the 10-mer peptide cross-link, although small amounts of T were added by hPol η. Molecular dynamics simulation of an hPol κ ternary complex containing a template-primer DNA with dCTP opposite the 10-mer peptide DPC revealed that this bulky lesion can be accommodated in the polymerase active site by aligning with the major groove of the adducted DNA within the ternary complex of polymerase and dCTP.

Ionization Properties of Histidine Residues in the Lipid Bilayer Membrane Environment

We address the critically important ionization properties of histidine side chains of membrane proteins, when exposed directly to lipid acyl chains within lipid bilayer membranes. The problem is important for addressing general principles that may underlie membrane protein function. To this end, we have employed a favorable host peptide framework provided by GWALP23 (acetyl-GGALW(5)LALALALALALALW(19)LAGA-amide). We inserted His residues into position 12 or 14 of GWALP23 (replacing either Leu(12) or Leu(14)) and incorporated specific [(2)H]Ala labels within the helical core sequence. Solid-state (2)H NMR spectra report the folding and orientation of the core sequence, revealing marked differences in the histidine-containing transmembrane helix behavior between acidic and neutral pH conditions. At neutral pH, the GWALP23-H12 and GWALP23-H14 helices exhibit well defined tilted transmembrane orientations in dioleoylphosphatidylcholine (DOPC)and dilauroylphosphatidylcholine (DLPC) bilayer membranes. Under acidic conditions, when His(12) is protonated and charged, the GWALP23-H12 helix exhibits a major population that moves to the DOPC bilayer surface and a minor population that occupies multiple transmembrane states. The response to protonation of His(14) is an increase in helix tilt, but GWALP23-H14 remains in a transmembrane orientation. The results suggest pKa values of less than 3 for His(12) and about 3-5 for His(14) in DOPC membranes. In the thinner DLPC bilayers, with increased water access, the helices are less responsive to changes in pH. The combined results enable us to compare the ionization properties of lipid-exposed His, Lys, and Arg side chains in lipid bilayer membranes.

Lysine and the Na+/K+ Selectivity in Mammalian Voltage-Gated Sodium Channels.

Voltage-gated sodium (Nav) channels are critical in the generation and transmission of neuronal signals in mammals. The crystal structures of several prokaryotic Nav channels determined in recent years inspire the mechanistic studies on their selection upon the permeable cations (especially between Na+ and K+ ions), a property that is proposed to be mainly determined by residues in the selectivity filter. However, the mechanism of cation selection in mammalian Nav channels lacks direct explanation at atomic level due to the difference in amino acid sequences between mammalian and prokaryotic Nav homologues, especially at the constriction site where the DEKA motif has been identified to determine the Na+/K+ selectivity in mammalian Nav channels but is completely absent in the prokaryotic counterparts. Among the DEKA residues, Lys is of the most importance since its mutation to Arg abolishes the Na+/K+ selectivity. In this work, we modeled the pore domain of mammalian Nav channels by mutating the four residues at the constriction site of a prokaryotic Nav channel (NavRh) to DEKA, and then mechanistically investigated the contribution of Lys in cation selection using molecular dynamics simulations. The DERA mutant was generated as a comparison to understand the loss of ion selectivity caused by the K-to-R mutation. Simulations and free energy calculations on the mutants indicate that Lys facilitates Na+/K+ selection by electrostatically repelling the cation to a highly Na+-selective location sandwiched by the carboxylate groups of Asp and Glu at the constriction site. In contrast, the electrostatic repulsion is substantially weakened when Lys is mutated to Arg, because of two intrinsic properties of the Arg side chain: the planar geometric design and the sparse charge distribution of the guanidine group.