Backbone Atoms Research Articles

Polylutidines: Multifunctional Surfaces through Vapor-Based Polymerization of Substituted Pyridinophanes.

We report a new class of functionalized polylutidine polymers that are prepared by chemical vapor deposition polymerization of substituted [2](1,4)benzeno[2](2,5)pyridinophanes. To prepare sufficient amounts of monomer for CVD polymerization, a new synthesis route for ethynylpyridinophane has been developed in three steps with an overall yield of 59 %. Subsequent CVD polymerization yielded well-defined films of poly(2,5-lutidinylene-co-p-xylylene) and poly(4-ethynyl-2,5-lutidinylene-co-p-xylylene). All polymers were characterized by infrared reflection-absorption spectroscopy, ellipsometry, contact angle studies, and X-ray photoelectron spectroscopy. Moreover, ζ-potential measurements revealed that polylutidine films have higher isoelectric points than the corresponding poly-xylylene surfaces owing to the nitrogen atoms in the polymer backbone. The availability of reactive alkyne groups on the surface of poly(4-ethynyl-2,5-lutidinylene-co-p-xylylene) coatings was confirmed by spatially controlled surface modification by means of Huisgen 1,3-dipolar cycloaddition. Compared to the more hydrophobic poly-p-xylylyenes, the presence of the heteroatom in the polymer backbone of polylutidine polymers resulted in surfaces that supported an increased adhesion of primary human umbilical vein endothelial cells (HUVECs). Vapor-based polylutidine coatings are a new class of polymers that feature increased hydrophilicity and increased cell adhesion without limiting the flexibility in selecting appropriate functional side groups.

High-Performance Multifunctional Fuel Additives Derived from Renewable Fatty Acids and Phosphazene

The increasing need for high-quality diesel fuel can be fulfilled by using environmental friendly, high-performance additives. We have synthesized novel fatty acid-derived multifunctional additives coupled with phosphazene moiety having the backbone of alternative nitrogen and phosphorous atoms which provided stability, while the additives possessed affinity for the metal surface to form a surface-complex film leading to the hindered metal contact and providing friction-reducing and antiwear properties. In addition, the applicability of the synthesized additives was investigated by testing the anti-oxidant and anticorrosion properties in polyethylene glycol (PEG200) which was taken as the reference base fluid. A rotating bomb oxidation test was used to evaluate the anti-oxidant characteristics, and a four-ball test for the tribological properties. Among the synthesized additives, phosphazene ester of erucic acid exhibited the best antiwear and antifriction properties; however, phosphazene ester of oleic acid was found to be the best anti-oxidant additive.

On the relationship between residue structural environment and sequence conservation in proteins.

Residues that are crucial to protein function or structure are usually evolutionarily conserved. To identify the important residues in protein, sequence conservation is estimated, and current methods rely upon the unbiased collection of homologous sequences. Surprisingly, our previous studies have shown that the sequence conservation is closely correlated with the weighted contact number (WCN), a measure of packing density for residue's structural environment, calculated only based on the Cα positions of a protein structure. Moreover, studies have shown that sequence conservation is correlated with environment-related structural properties calculated based on different protein substructures, such as a protein's all atoms, backbone atoms, side-chain atoms, or side-chain centroid. To know whether the Cα atomic positions are adequate to show the relationship between residue environment and sequence conservation or not, here we compared Cα atoms with other substructures in their contributions to the sequence conservation. Our results show that Cα positions are substantially equivalent to the other substructures in calculations of various measures of residue environment. As a result, the overlapping contributions between Cα atoms and the other substructures are high, yielding similar structure-conservation relationship. Take the WCN as an example, the average overlapping contribution to sequence conservation is 87% between Cα and all-atom substructures. These results indicate that only Cα atoms of a protein structure could reflect sequence conservation at the residue level.

A computational study to identify the key residues of peroxisome proliferator-activated receptor gamma in the interactions with its antagonists

Peroxisome proliferator-activated receptors (PPARs) compose a family of nuclear receptors, PPARα, PPARβ, and PPARγ, which mediate the effects of lipidic ligands at the transcriptional level. Among these, the PPARγ has been known to regulate adipocyte differentiation, fatty acid storage and glucose metabolism, and is a target of antidiabetic drugs. In this work, the interactions between PPARγ and its six known antagonists were investigated using computational methods such as molecular docking, molecular dynamics (MD) simulations, and the hybrid quantum mechanics/molecular mechanics (QM/MM). The binding energies evaluated by molecular docking varied between −22.59 and −35.15 kJ mol− 1. In addition, MD simulations were performed to investigate the binding modes and PPARγ conformational changes upon binding of antagonists. Analysis of the root-mean-square fluctuations (RMSF) of backbone atoms shows that H3 of PPARγ has a higher mobility in the absence of antagonists and moderate conformational changes were observed. The interaction energies between antagonists and each PPARγ residue involved in the interactions were studied by QM/MM calculations. These calculations reveal that antagonists with different structures show different interaction energies with the same residue of PPARγ. Therefore, it can be concluded that the key residues vary depending on the structure of the ligand, which binds to PPARγ.

1H, 13C and 15N NMR chemical shift assignments of A. thaliana RCD1 RST.

The A. thaliana RCD1 (radical-induced cell death1) protein is a cellular signaling hub protein which interacts with numerous plant transcription factors from different families. It consists of three conserved domains and intervening unstructured regions, the C-terminal RST domain being responsible for the interactions with the transcription factors. It has been shown that many partner proteins interact with RCD1 RST via their intrinsically disordered regions, and that the domain is able to house partners with divergent folds. We aim to structurally characterize the RCD1 RST domain and its complexes [complex with DREB2A]. Here we report the 1H, 15N and 13C chemical shift assignments of the backbone and sidechain atoms for RCD1 (468-589) containing the RST (510-567) domain.

Lid opening and conformational stability of T1 Lipase is mediated by increasing chain length polar solvents.

The dynamics and conformational landscape of proteins in organic solvents are events of potential interest in nonaqueous process catalysis. Conformational changes, folding transitions, and stability often correspond to structural rearrangements that alter contacts between solvent molecules and amino acid residues. However, in nonaqueous enzymology, organic solvents limit stability and further application of proteins. In the present study, molecular dynamics (MD) of a thermostable Geobacillus zalihae T1 lipase was performed in different chain length polar organic solvents (methanol, ethanol, propanol, butanol, and pentanol) and water mixture systems to a concentration of 50%. On the basis of the MD results, the structural deviations of the backbone atoms elucidated the dynamic effects of water/organic solvent mixtures on the equilibrium state of the protein simulations in decreasing solvent polarity. The results show that the solvent mixture gives rise to deviations in enzyme structure from the native one simulated in water. The drop in the flexibility in H2O, MtOH, EtOH and PrOH simulation mixtures shows that greater motions of residues were influenced in BtOH and PtOH simulation mixtures. Comparing the root mean square fluctuations value with the accessible solvent area (SASA) for every residue showed an almost correspondingly high SASA value of residues to high flexibility and low SASA value to low flexibility. The study further revealed that the organic solvents influenced the formation of more hydrogen bonds in MtOH, EtOH and PrOH and thus, it is assumed that increased intraprotein hydrogen bonding is ultimately correlated to the stability of the protein. However, the solvent accessibility analysis showed that in all solvent systems, hydrophobic residues were exposed and polar residues tended to be buried away from the solvent. Distance variation of the tetrahedral intermediate packing of the active pocket was not conserved in organic solvent systems, which could lead to weaknesses in the catalytic H-bond network and most likely a drop in catalytic activity. The conformational variation of the lid domain caused by the solvent molecules influenced its gradual opening. Formation of additional hydrogen bonds and hydrophobic interactions indicates that the contribution of the cooperative network of interactions could retain the stability of the protein in some solvent systems. Time-correlated atomic motions were used to characterize the correlations between the motions of the atoms from atomic coordinates. The resulting cross-correlation map revealed that the organic solvent mixtures performed functional, concerted, correlated motions in regions of residues of the lid domain to other residues. These observations suggest that varying lengths of polar organic solvents play a significant role in introducing dynamic conformational diversity in proteins in a decreasing order of polarity.

Exploring the RNA-bound and RNA-free human Argonaute-2 by molecular dynamics simulation method.

Argonaute 2 (Ago2) protein is the major vehicle of microRNAs (miRNAs)-guided gene repression and silencing processes. Although the crystal structure of human Ago2 (hAgo2) has recently been disclosed, the information of dynamically structural character of protein-RNA recognition is still lacking. Molecular dynamics simulations were used to systematically explore hAgo2 in the presence and absence of RNA duplex. Stable direct and water-mediated hydrogen bonds were observed between guide RNA backbone atoms and hAgo2, especially for nucleotides 2-7. In addition, water-mediated hydrogen bonds are indicated to be critical in the specific recognition between hAgo2 and the conserved adenine in position 1 of target RNA. The core domains (N, PAZ, MID, and PIWI) possess rigid body movements during the simulations. The motions of N-PAZ and PIWI-MID are negatively correlated with or without RNA binding and PAZ domain is identified as the most mobile domain in both systems. The reorientation of PAZ domain not only influences the binding of helix-7 and RNA duplex, the initial pairing process, but also the shape of N-PAZ cleft, where the supplemental base pairing occurs. It is speculated that PAZ domain could be a key regulator in hAgo2-mediated miRNA-induced gene regulation.

Deciphering the Molecular Effects of Mutations on ATRX Cause ATRX Syndrome: A Molecular Dynamics Study.

α-thalassemia mental retardation X-linked (ATRX) syndrome is caused by the dysfunction of ATRFfigX protein. The present study explored the structural consequences influenced by two observed mutations V194I and C220R on ADD domain of ATRX protein by applying all atom molecular dynamics (MD) simulation. MD result showed that both the mutants exhibited wide variations in their backbone dynamics, as a result, mutant V210I showed complete distortion on α3 and the mutant C220R displayed a biased disruption on α2-3. The interference in the local folding of α-helices in both the mutants resulted by the loss of hydrogen bonds mediated by the backbone atoms. Principle component analysis (PCA) elucidated that both the mutants endured a diverse conformational dynamics, consequently adopted thermodynamically different conformational state. Besides, binding residues in both the mutants showed more structural disorder, thereby unable to recognize the hallmark modification, K9me3 (tri-methylated lysine at position 9) of histone H3 peptide and it was not conducive for the wild type ADD domain like functionality. Altogether, our findings provide knowledge to understand the structural and functional relationship of disease-associated mutations, V194I and C220R on ADD domain as well as gain further insights into the molecular pathogenesis of ATRX syndrome. J. Cell. Biochem. 118: 3318-3327, 2017. © 2017 Wiley Periodicals, Inc.

Design and Development of Highly Potent HIV-1 Protease Inhibitors with a Crown-Like Oxotricyclic Core as the P2-Ligand To Combat Multidrug-Resistant HIV Variants.

Design, synthesis, and evaluation of a new class of exceptionally potent HIV-1 protease inhibitors are reported. Inhibitor 5 displayed superior antiviral activity and drug-resistance profiles. In fact, this inhibitor showed several orders of magnitude improved antiviral activity over the FDA approved drug darunavir. This inhibitor incorporates an unprecedented 6-5-5 ring-fused crown-like tetrahydropyranofuran as the P2 ligand and an aminobenzothiazole as the P2' ligand with the (R)-hydroxyethylsulfonamide isostere. The crown-like P2 ligand for this inhibitor has been synthesized efficiently in an optically active form using a chiral Diels-Alder catalyst providing a key intermediate in high enantiomeric purity. Two high resolution X-ray structures of inhibitor-bound HIV-1 protease revealed extensive interactions with the backbone atoms of HIV-1 protease and provided molecular insight into the binding properties of these new inhibitors.

In-depth study of DNA binding of Cys2His2 finger domains in testis zinc-finger protein.

Previously, we identified that both fingers 1 and 2 in the three Cys2His2 zinc-finger domains (TZD) of testis zinc-finger protein specifically bind to its cognate DNA; however, finger 3 is non-sequence–specific. To gain insights into the interaction mechanism, here we further investigated the DNA-binding characteristics of TZD bound to non-specific DNAs and its finger segments bound to cognate DNA. TZD in non-specific DNA binding showed smaller chemical shift perturbations, as expected. However, the direction of shift perturbation, change of DNA imino-proton NMR signal, and dynamics on the 15N backbone atom significantly differed between specific and non-specific binding. Using these unique characteristics, we confirmed that the three single-finger segments (TZD1, TZD2 and TZD3) and the two-finger segment (TZD23) non-specifically bind to the cognate DNA. In comparison, the other two-finger segment (TZD12) binding to the cognate DNA features simultaneous non-specific and semi-specific binding, both slowly exchanged in terms of NMR timescale. The process of TZD binding to the cognate DNA is likely stepwise: initially TZD non-specifically binds to DNA, then fingers 1 and 2 insert cooperatively into the major groove of DNA by semi-specific binding, and finally finger 3 non-specifically binds to DNA, which promotes the specific binding on fingers 1 and 2 and stabilizes the formation of a specific TZD–DNA complex.

The hydrogen peroxide hypersensitivity of OxyR2 in Vibrio vulnificus depends on conformational constraints

Most Gram-negative bacteria respond to excessive levels of H2O2 using the peroxide-sensing transcriptional regulator OxyR, which can induce the expression of antioxidant genes to restore normality. Vibrio vulnificus has two distinct OxyRs (OxyR1 and OxyR2), which are sensitive to different levels of H2O2 and induce expression of two different peroxidases, Prx1 and Prx2. Although OxyR1 has both high sequence similarity and H2O2 sensitivity comparable with that of other OxyR proteins, OxyR2 exhibits limited sequence similarity and is more sensitive to H2O2 To investigate the basis for this difference, we determined crystal structures and carried out biochemical analyses of OxyR2. The determined structure of OxyR2 revealed a flipped conformation of the peptide bond before Glu-204, a position occupied by glycine in other OxyR proteins. Activity assays showed that the sensitivity to H2O2 was reduced to the level of other OxyR proteins by the E204G mutation. We solved the structure of the OxyR2-E204G mutant with the same packing environment. The structure of the mutant revealed a dual conformation of the peptide bond before Gly-204, indicating the structural flexibility of the region. This structural duality extended to the backbone atoms of Gly-204 and the imidazole ring of His-205, which interact with H2O2 and invariant water molecules near the peroxidatic cysteine, respectively. Structural comparison suggests that Glu-204 in OxyR2 provides rigidity to the region that is important in H2O2 sensing, compared with the E204G structure or other OxyR proteins. Our findings provide a structural basis for the higher sensitivity of OxyR2 to H2O2 and also suggest a molecular mechanism for bacterial regulation of expression of antioxidant genes at divergent concentrations of cellular H2O2.

Solution structure and dynamics of the outer membrane cytochrome OmcF from Geobacter sulfurreducens

Gene knock-out studies on Geobacter sulfurreducens cells showed that the outer membrane-associated monoheme cytochrome OmcF is involved in respiratory pathways leading to the extracellular reduction of Fe(III) and U(VI). In addition, microarray analysis of an OmcF-deficient mutant revealed that many of the genes with decreased transcript level were those whose expression is up-regulated in cells grown with a graphite electrode as electron acceptor, suggesting that OmcF also regulates the electron transfer to electrode surfaces and the concomitant electricity production by G. sulfurreducens in microbial fuel cells. 15N,13C–labeled OmcF was produced and NMR spectroscopy was used to determine the solution structure of the protein in the fully reduced state and the pH-dependent conformational changes. In addition, 15N relaxation NMR experiments were used to characterize the overall and internal backbone dynamics of OmcF. The structure obtained is well-defined, with an average pairwise root mean square deviation of 0.37Å for the backbone atoms and 0.98Å for all heavy atoms. For the first time a solution structure and the protein motions were determined for an outer membrane cytochrome from G. sulfurreducens, which constitutes an important step to understand the extracellular electron transfer mechanism in Geobacter cells.

Backbone assignment of perdeuterated proteins by solid-state NMR using proton detection and ultrafast magic-angle spinning.

Solid-state NMR (ssNMR) is a technique that allows the study of protein structure and dynamics at atomic detail. In contrast to X-ray crystallography and cryo-electron microscopy, proteins can be studied under physiological conditions-for example, in a lipid bilayer and at room temperature (0-35 °C). However, ssNMR requires considerable amounts (milligram quantities) of isotopically labeled samples. In recent years, 1H-detection of perdeuterated protein samples has been proposed as a method of alleviating the sensitivity issue. Such methods are, however, substantially more demanding to the spectroscopist, as compared with traditional 13C-detected approaches. As a guide, this protocol describes a procedure for the chemical shift assignment of the backbone atoms of proteins in the solid state by 1H-detected ssNMR. It requires a perdeuterated, uniformly 13C- and 15N-labeled protein sample with subsequent proton back-exchange to the labile sites. The sample needs to be spun at a minimum of 40 kHz in the NMR spectrometer. With a minimal set of five 3D NMR spectra, the protein backbone and some of the side-chain atoms can be completely assigned. These spectra correlate resonances within one amino acid residue and between neighboring residues; taken together, these correlations allow for complete chemical shift assignment via a 'backbone walk'. This results in a backbone chemical shift table, which is the basis for further analysis of the protein structure and/or dynamics by ssNMR. Depending on the spectral quality and complexity of the protein, data acquisition and analysis are possible within 2 months.

Size Matters: Influence of Gold‐to‐Ligand Ratio and Sulfur‐Sulfur Distance of Linear Thioether Heptamers on the Size of Gold Nanoparticles

A systematic investigation of two parameters steering the size of linear octadentate heptamer‐coated gold nanoparticles (AuNPs) is presented, being i) the chemical structure (sulfur‐sulfur distance) of the coating thioether heptamer ligand and ii) the ratio of ligand to tetrachloroauric acid (HAuCl4) reduced during the formation of the AuNPs. For this purpose, a novel terphenyl‐based thioether heptamer (Ter) is synthesized via an end‐capping oligomerization strategy, comprising an increased distance between neighboring sulfur atoms in the ligand backbone compared to the meta‐xylene‐ (Xyl) and tetraphenylmethane‐ (TPM) based heptamers. While for both investigated parameters a clear trend to various‐sized NPs is shown, a stronger influence in the resulting sizes is observed by alteration of ligand to gold‐ratio. Remarkable processability‐ and long‐term stability‐features were observed for AuNPs stabilized by the bulky tetraphenylmethane‐based heptamer (TPM).

Residual Structure in the Denatured State of a Three-Helix Bundle Protein

The denatured state of a protein is canonically described as an ensemble of non-interacting random coil conformations. By studying a helical bundle protein under denaturing conditions, we have found that residual structure persists in its denatured state. Our work uses the upstream Ubiquitin-Associated domain, UBA(1), of the HHR23A protein as a model system to study residual structure in the denatured state ensemble. Multi-dimensional nuclear magnetic resonance (NMR) experiments 1H-15N HSQC, HNCO, and HNCA were used to find backbone atom secondary chemical shifts which correspond to protein secondary structure with amino acid resolution. By comparing chemical shifts in moderate- to highly denaturing conditions of 4M, 5M, and 6M guanidine HCl to a reference state of 7M GdnHCl, UBA(1) shows residual helical content among all three helices. In contrast, the helical propensity is predicted by the program Agadir to be higher in helix 2 and negligible elsewhere, with the assumption that helices do not interact. The presence of relatively uniform helical content among all three helices may indicate stabilizing tertiary interactions despite highly denaturing conditions. To further probe tertiary interactions, single-residue mutants of UBA(1) will be used to determine if residues in the hydrophobic core or turn regions are contributing to stability of residual structure in the denatured state.

Error assessment in molecular dynamics trajectories using computed NMR chemical shifts

Accurate chemical shifts for the atoms in molecular mechanics (MD) trajectories can be obtained from quantum mechanical (QM) calculations that depend solely on the coordinates of the atoms in the localized regions surrounding atoms of interest. If these coordinates are correct and the sample size is adequate, the ensemble average of these chemical shifts should be equal to the chemical shifts obtained from NMR spectroscopy. If this is not the case, the coordinates must be incorrect. We have utilized this fact to quantify the errors associated with the backbone atoms in MD simulations of proteins. A library of regional conformers containing 169,499 members was constructed from 6 model proteins. The chemical shifts associated with the backbone atoms in each of these conformers was obtained from QM calculations using density functional theory at the B3LYP level with a 6-311+G(2d,p) basis set. Chemical shifts were assigned to each backbone atom in each MD simulation frame using a template matching approach. The ensemble average of these chemical shifts was compared to chemical shifts from NMR spectroscopy. A large systematic error was identified that affected the 1H atoms of the peptide bonds involved in hydrogen bonding with water molecules or peptide backbone atoms. This error was highly sensitive to changes in electrostatic parameters. Smaller errors affecting the 13Ca and 15N atoms were also detected. We believe these errors could be useful as metrics for comparing the force-fields and parameter sets used in MD simulation because they are directly tied to errors in atomic coordinates.

Backbone resonance assignments of the F-actin binding domain of mouse αN-catenin.

α-Catenin is a filamentous actin (F-actin) binding protein that links the classical cadherin-catenin complex to the actin cytoskeleton at adherens junctions (AJs). Its C-terminal F-actin binding domain is required for regulating the dynamic interaction between AJs and the actin cytoskeleton during tissue development. Thus, obtaining the molecular details of this interaction is a crucial step towards understanding how α-catenin plays critical roles in biological processes, such as morphogenesis, cell polarity, wound healing and tissue maintenance. Here we report the backbone atom (1HN, 15N, 13Cα, 13Cβ and 13C') resonance assignments of the C-terminal F-actin binding domain of αN-catenin.

Examination of Trifluoroethanol Interactions with Trp-Cage through MD Simulations and Intermolecular Nuclear Overhauser Effects.

Molecular dynamics simulations of the protein model Trp-cage in 42% 2,2,2-trifluoroethanol (TFE)-water at 318 K have been carried out with the goal of exploring solvent fluorine-peptide hydrogen nuclear spin cross-relaxation. TFE5 and TFE6 models of TFE developed in previous work from this laboratory were used with the TIP5PE model of water. System densities and component translational diffusion coefficients were well predicted by the simulations, as were many of the cross-relaxation parameters (ΣHF) for which experimental values are available. However, the calculated ΣHF values were too small for some hydrogens of the Trp6 indole ring and for amino acid hydrogens near this residue in the native structure. Simulations carried out with unfolded versions of Trp-cage suggest that underestimates of ΣHF are the result of insufficient sampling of conformational motions that expose these hydrogens to interactions with solvent molecules during simulations of the native peptide. Consideration of the relative amounts of TFE and water surrounding the Trp-cage structure indicates that the composition of the solvent mixture at distances beyond ∼1.5 nm from the surface of the peptide is close to the composition of the bulk solvent, but, as observed by others, TFE tends to accumulate preferentially near the surface of the peptide. Both TFE and water molecules make contacts with the peptide surface; water molecules predominate in contacts with the peptide backbone atoms, and TFE molecules generally preferentially interact with side chains. Translational diffusion of solvent molecules appears to slow down near the surface of the peptide.

Conformation of Pyrene-Labeled Amylose in DMSO Characterized with the Fluorescence Blob Model

Amylose and poly(methyl acrylate) were randomly labeled with pyrene to yield a series of Py-Amylose and Py-PMA constructs, and their ability to form excimer in DMSO was characterized quantitatively by steady-state and time-resolved fluorescence. First, the ratio of the fluorescence intensity of the excimer over that of the monomer, namely the IE/IM ratio, was obtained from the fluorescence spectra. Second, the product ⟨kblob × Nblob⟩ was obtained from the fluorescence blob model (FBM) analysis of the fluorescence decays. Both IE/IM and ⟨kblob × Nblob⟩ yielded similar values when expressed in terms of moles of pyrene per backbone atom for Py-Amylose and Py-PMA. Since IE/IM and ⟨kblob × Nblob⟩ reflect the efficiency of pyrene excimer formation, the similar behavior observed for both parameters obtained for rigid amylose and flexible PMA could only be rationalized by postulating that amylose adopted a compact helical conformation in DMSO. To confirm whether this was possible, molecular mechanics optimization...

Identifying the elusive link between amino acid sequence and charge selectivity in pentameric ligand-gated ion channels

Among neurotransmitter-gated ion channels, the superfamily of pentameric ligand-gated ion channels (pLGICs) is unique in that its members display opposite permeant-ion charge selectivities despite sharing the same structural fold. Although much effort has been devoted to the identification of the mechanism underlying the cation-versus-anion selectivity of these channels, a careful analysis of past work reveals that discrepancies exist, that different explanations for the same phenomenon have often been put forth, and that no consensus view has yet been reached. To elucidate the molecular basis of charge selectivity for the superfamily as a whole, we performed extensive mutagenesis and electrophysiological recordings on six different cation-selective and anion-selective homologs from vertebrate, invertebrate, and bacterial origin. We present compelling evidence for the critical involvement of ionized side chains-whether pore-facing or buried-rather than backbone atoms and propose a mechanism whereby not only their charge sign but also their conformation determines charge selectivity. Insertions, deletions, and residue-to-residue mutations involving nonionizable residues in the intracellular end of the pore seem to affect charge selectivity by changing the rotamer preferences of the ionized side chains in the first turn of the M2 α-helices. We also found that, upon neutralization of the charged residues in the first turn of M2, the control of charge selectivity is handed over to the many other ionized side chains that decorate the pore. This explains the long-standing puzzle as to why the neutralization of the intracellular-mouth glutamates affects charge selectivity to markedly different extents in different cation-selective pLGICs.