Contribution Of Side Chains Research Articles

Identification of hydrogen bonding network for proton transfer at the quinol oxidation site of Rhodobacter capsulatus cytochrome bc1

Cytochrome bc1 catalyzes electron transfer from quinol (QH2) to cytochrome c in reactions coupled to proton translocation across the energy-conserving membrane. Energetic efficiency of the catalytic cycle is secured by a two-electron, two-proton bifurcation reaction leading to oxidation of QH2 and reduction of the Rieske cluster and heme bL. The proton paths associated with this reaction remain elusive. Here we used site-directed mutagenesis and quantum mechanical (QM) calculations to analyze the contribution of protonable side chains located at the heme bL side of the quinol oxidation site (Qo) in Rhodobacter capsulatus cytochrome bc1. We observe that the proton path is effectively switched off when H276 and E295 are simultaneously mutated to the non-protonable residues in the H276F/E295V double mutant. The two single mutants, H276F or E295V, are less efficient, but still transfer protons at functionally-relevant rates. Natural selection exposed two single mutations, N279S and M154T, that restored the functional proton transfers in H276F/E295V. QM calculations indicated that H276F/E295V traps the side chain of Y147 in a position distant from QH2, while either N279S or M154T induce local changes releasing Y147 from that position. This shortens the distance between the protonable groups of Y147 and D278 and/or increases mobility of the Y147 side chain, which makes Y147 efficient in transferring protons from QH2 toward D278 in H276F/E295V. Overall, our study identified an extended hydrogen bonding network, build up by E295, H276, D278 and Y147, involved in efficient removal of the proton from QH2 at the heme bL side of Qo.

Biophysical Interpretation of Evolutionary Consequences on the SARS-CoV2 Main Protease through Molecular Dynamics Simulations and Network Topology Analysis.

In this study, we present a combined analysis procedure between atomistic molecular dynamics (MD) simulations and network topology to obtain more understanding on the evolutionary consequences on protein stability and substrate binding of the main protease enzyme of SARS-CoV2. Communicability matrices of the protein residue networks (PRNs) were extracted from MD trajectories of both Mpro enzymes in complex with the nsp8/9 peptide substrate to compare the local communicability within both proteases that would affect the enzyme function, along with biophysical details on global protein conformation, flexibility, and contribution of amino acid side chains to both intramolecular and intermolecular interactions. The analysis displayed the significance of the mutated residue 46 with the highest communicability gain to the binding pocket closure. Interestingly, the mutated residue 134 with the highest communicability loss corresponded to a local structural disruption of the adjacent peptide loop. The enhanced flexibility of the disrupted loop connecting to the catalytic residue Cys145 introduced an extra binding mode that brought the substrate in proximity and could facilitate the reaction. This understanding might provide further help in the drug development strategy against SARS-CoV2 and prove the capability of the combined techniques of MD simulations and network topology analysis as a "reverse" protein engineering tool.

Evaluation of interactions between the hepatitis C virus NS3/4A and sulfonamidobenzamide based molecules using molecular docking, molecular dynamics simulations and binding free energy calculations.

The Hepatitis C Virus (HCV) NS3/4A is an attractive target for the treatment of Hepatitis C infection. Herein, we present an investigation of HCV NS3/4A inhibitors based on a sulfonamidobenzamide scaffold. Inhibitor interactions with HCV NS3/4A were explored by molecular docking, molecular dynamics simulations, and MM/PBSA binding free energy calculations. All of the inhibitors adopt similar molecular docking poses in the catalytic site of the protease that are stabilized by hydrogen bond interactions with G137 and the catalytic S139, which are known to be important for potency and binding stability. The quantitative assessments of binding free energies from MM/PBSA correlate well with the experimental results, with a high coefficient of determination, R2 of 0.92. Binding free energy decomposition analyses elucidate the different contributions of Q41, F43, H57, R109, K136, G137, S138, S139, A156, M485, and Q526 in binding different inhibitors. The importance of these sidechain contributions was further confirmed by computational alanine scanning mutagenesis. In addition, the sidechains of K136 and S139 show crucial but distinct contributions to inhibitor binding with HCV NS3/4A. The structural basis of the potency has been elucidated, demonstrating the importance of the R155 sidechain conformation. This extensive exploration of binding energies and interactions between these compounds and HCV NS3/4A at the atomic level should benefit future antiviral drug design.

Y98 Mutation Leads to the Loss of RsfS Anti-Association Activity in Staphylococcus aureus.

Ribosomal silencing factor S (RsfS) is a conserved protein that plays a role in the mechanisms of ribosome shutdown and cell survival during starvation. Recent studies demonstrated the involvement of RsfS in the biogenesis of the large ribosomal subunit. RsfS binds to the uL14 ribosomal protein on the large ribosomal subunit and prevents its association with the small subunit. Here, we estimated the contribution of RsfS amino acid side chains at the interface between RsfS and uL14 to RsfS anti-association function in Staphylococcus aureus through in vitro experiments: centrifugation in sucrose gradient profiles and an S. aureus cell-free system assay. The detected critical Y98 amino acid on the RsfS surface might become a new potential target for pharmacological drug development and treatment of S. aureus infections.

Gaining Insights on the Interactions of a Class of Decorated (2-([2,2'-Bipyridin]-6-yl)phenyl)platinum Compounds with c-Myc Oncogene Promoter G-Quadruplex and Other DNA Structures.

Organometallic molecules offer some of the most promising scaffolds for interaction with G-quadruplex nucleic acids. We report the efficient synthesis of a family of organoplatinum(II) complexes, featuring a 2-([2,2'-bipyridin]-6-yl)phenyl tridentate (N∧ N∧ C) ligand, that incorporates peripheral side-chains aiming at enhancing and diversifying its interaction capabilities. These include a di-isopropyl carbamoyl amide, a morpholine ethylenamide, two enantiomeric proline imides and an oxazole. The binding affinities of the Pt-complexes were evaluated via UV-vis and fluorescence titrations, against 5 topologically-distinct DNA structures, including c-myc G-quadruplex, two telomeric (22AG) G-quadruplexes, a duplex (ds26) and a single-stranded (polyT) DNA. All compounds exhibited binding selectivity in favour of c-myc, with association constants (Ka ) in the range of 2-5×105 M-1 , lower affinity for both folds of 22AG and for ds26 and negligible affinity for polyT. Remarkable emission enhancements (up to 200-fold) upon addition of excess DNA were demonstrated by a subset of the compounds with c-myc, providing a basis for optical selectivity, since optical response to all other tested DNAs was low. A c-myc DNA-melting experiment showed significant stabilizing abilities for all compounds, with the most potent binder, the morpholine-Pt-complex, exhibiting a ΔTm >30 °C, at 1 : 5 DNA-to-ligand molar ratio. The same study implied contributions of the diverse side-chains to helix stabilization. To gain direct evidence of the nature of the interactions, mixtures of c-myc with the four most promising compounds were studied via UV Resonance Raman (UVRR) spectroscopy, which revealed end-stacking binding mode, combined with interactions of side-chains with loop nucleobase residues. Docking simulations were conducted to provide insights into the binding modes for the same four Pt-compounds, suggesting that the binding preference for two alternative orientations of the c-myc G-quadruplex thymine 'cap' ('open' vs. 'closed'), as well as the relative contributions to affinity from end-stacking and H-bonding, are highly dependent on the nature of the interacting Pt-complex side-chain.

Free-energy landscapes and insertion pathways for peptides in membrane environment.

Free-energy landscapes for short peptides-specifically for variants of the pH low insertion peptide (pHLIP)-in the heterogeneous environment of a lipid bilayer or cell membrane are constructed, taking into account a set of dominant interactions and the conformational preferences of the peptide backbone. Our methodology interprets broken internal H-bonds along the backbone of a polypeptide as statistically interacting quasiparticles, activated from the helix reference state. The favored conformation depends on the local environment (ranging from polar to nonpolar), specifically on the availability of external H-bonds (with H_{2}O molecules or lipid headgroups) to replace internal H-bonds. The dominant side-chain contribution is accounted for by residue-specific transfer free energies between polar and nonpolar environments. The free-energy landscape is sensitive to the level of pH in the aqueous environment surrounding the membrane. For high pH, we identify pathways of descending free energy that suggest a coexistence of membrane-adsorbed peptides with peptides in solution. A drop in pH raises the degree of protonation of negatively charged residues and thus increases the hydrophobicity of peptide segments near the C terminus. For low pH, we identify insertion pathways between the membrane-adsorbed state and a stable trans-membrane state with the C terminus having crossed the membrane.

Phosphatidylserine and Phosphatidylethanolamine Asymmetry Have a Negligible Effect on the Global Structure, Dynamics, and Interactions of the KRAS Lipid Anchor.

The intrinsically disordered C-terminus of the prominent oncogenic protein KRAS-4B (KRAS) selectively interacts and clusters with phosphatidylserine (PS) lipids in the plasma membrane (PM). This 11-residue segment, called tK, contains a polybasic domain (PBD) of six contiguous lysine residues and a farnesylated cysteine. Previous molecular dynamics (MD) simulation studies of tK in phosphatidylcholine (PC)/PS bilayers have suggested that backbone conformational dynamics modulate tK-PS interactions. These simulations have been conducted in symmetric membranes whereas the PM is compositionally asymmetric, with the inner leaflet, where KRAS is localized, being enriched with PS and phosphatidylethanolamine (PE) lipids. To examine if bilayer asymmetry affects tK conformational dynamics and interaction with lipids, we conducted two 10 μs long MD simulations of tK bound to a PC/PS and a PC/PS/PE bilayer in which the PS and PE lipids are distributed in one leaflet. We found that, first, these compositional asymmetries caused differences in acyl chain dynamics between leaflets, but the equilibrium structural and dynamic properties of the two asymmetric bilayers are similar; second, in both systems tK is highly dynamic and samples at least two distinct conformational states; third, PS-tK hydrogen-bonding interactions vary with peptide backbone conformations, and lysine side chains in the PBD predominantly interact with the serine oxygens of PS. These results are in good agreement with previous observations of tK in symmetric membranes. The effects of POPS asymmetry or the presence of POPE on tK are limited to modulating the relative contribution of individual side chains to interactions with lipids and redistributing conformational substates. Additional observations include the larger flexibility of tK in the current simulations, which we attribute to the longer duration of the simulations and the use of the CHARMM36m force field, which more accurately models intrinsically disordered peptides such as tK.

A Comparison between the Lower Critical Solution Temperature Behavior of Polymers and Biomacromolecules

All-atom molecular dynamics (MD) simulations are employed to compare the lower critical solution temperature (LCST) behaviors of poly(N-isopropylacrylamide) (PNIPAM) and elastin-like polypeptides (ELPs) with the canonical Val-Pro-Gly-Val-Gly ((VPGVG)n) sequence over a range of temperatures from 280 K to 380 K. Our simulations suggest that the structure of proximal water dictates the conformation of both the (VPGVG)n ELPs and PNIPAM chains. Specifically, the LCST transition in ELPs can be attributed to a combination of thermal disruption of the network of the proximal water near both hydrophilic and hydrophobic groups in the backbone and side-chain of (VPGVG)n, resulting in a reduction in solvent accessible surface area (SASA). This is accompanied with an increase in the secondary structure above its LCST. In the case of PNIPAM, the LCST transition is a result of a combination of a reduction in the hydrophobic SASA primarily due to the contributions of isopropyl side-chain and less to the backbone and the formation of intra-chain hydrogen bonds between the amide groups on the side-chain above its LCST.

Functional mapping of the 14-3-3 hub protein as a guide to design 14-3-3 molecular glues.

Molecular glues represent an evolution in drug discovery, however, targeted stabilization of protein complexes remains challenging, owing to a paucity of drug design rules. The functional mapping of hotspots has been critical to protein-protein interaction (PPI) inhibitor research, however, the orthogonal approach to stabilize PPIs has not exploited this information. Utilizing the hub protein 14-3-3 as a case study we demonstrate that functional mapping of hotspots provides a triage map for 14-3-3 molecular glue development. Truncation and mutation studies allowed deconvoluting the energetic contributions of sidechain and backbone interactions of a 14-3-3-binding non-natural peptide. Three central 14-3-3 hotspots were identified and their thermodynamic characteristics profiled. In addition to the phospho-binding pocket; (i) Asn226, (ii) Lys122 and (iii) the hydrophobic patch formed by Leu218, Ile219 and Leu222 were critical for protein complex formation. Exploiting this hotspot information allowed a peptide-based molecular glue that elicits high cooperativity (α = 36) and selectively stabilizes the 14-3-3/ChREBP PPI to be uniquely developed.

Physicochemical study on molecular interactions of the active pharmaceutical ingredient ionic liquid domiphen salicylate with amino acids at different temperatures

The molecular interactions between the active pharmaceutical ingredient ionic liquid domiphen salicylate ([Dom][Sal]) and some amino acids (glycine, L-alanine, L-valine, L-leucine) in aqueous solution were studied. In order to accomplish this objective, the volume, conductivity and UV–vis spectroscopy of ternary systems were measured at different temperatures. The infinite dilution apparent molar volume, V2,φ0, transfer volume, ΔtV0, hydration number, nH and infinite dilution apparent molar expansibility, Eφ0 of amino acids in aqueous [Dom][Sal] solution were derived from the measured density data. Side chain contributions to V2,φ0 have been estimated. Critical micelle concentration, cmc, and thermodynamic properties of micellization (ΔGm0, ΔHm0,ΔSm0) for [Dom][Sal] in aqueous amino acids solutions were calculated using electrical conductivity data. The binding constant Kb of the amino acid with [Dom][Sal] was estimated based on the maximum absorbance data of the UV–vis spectroscopy. The impacts of temperature, concentration and the length of aliphatic chain of amino acids on these parameters were discussed. The mentioned parameters was employed to unveil the solute–solvent interactions in ternary systems.

Site-Specific Backbone and Side-Chain Contributions to Thermodynamic Stabilizing Forces of the WW Domain.

The native structure of a protein is stabilized by a number of interactions such as main-chain hydrogen bonds and side-chain hydrophobic contacts. However, it has been challenging to determine how these interactions contribute to protein stability at single amino acid resolution. Here, we quantified site-specific thermodynamic stability at the molecular level to extend our understanding of the stabilizing forces in protein folding. We derived the free energy components of individual amino acid residues separately for the folding of the human Pin WW domain based on simulated structures. A further decomposition of the thermodynamic properties into contributions from backbone and side-chain groups enabled us to identify the critical residues in the secondary structure and hydrophobic core formation, without introducing physical modifications to the system as in site-directed mutagenesis methods. By relating the structural and thermodynamic changes upon folding for each residue, we find that the simultaneous formation of the backbone hydrogen bonds and side-chain contacts cooperatively stabilizes the folded structure. The identification of stabilizing interactions in a folding protein at atomic resolution will provide molecular insights into understanding the origin of the protein structure and into engineering a more stable protein.

Amino acid side chain contribution to protein FTIR spectra: impact on secondary structure evaluation

Prediction of protein secondary structure from FTIR spectra usually relies on the absorbance in the amide I–amide II region of the spectrum. It assumes that the absorbance in this spectral region, i.e., roughly 1700–1500 cm−1 is solely arising from amide contributions. Yet, it is accepted that, on the average, about 20% of the absorbance is due to amino acid side chains. The present paper evaluates the contribution of amino acid side chains in this spectral region and the potential to improve secondary structure prediction after correcting for their contribution. We show that the β-sheet content prediction is improved upon subtraction of amino acid side chain contributions in the amide I–amide II spectral range. Improvement is relatively important, for instance, the error of prediction of β-sheet content decreases from 5.42 to 4.97% when evaluated by ascending stepwise regression. Other methods tested such as partial least square regression and support vector machine have also improved accuracy for β-sheet content evaluation. The other structures such as α-helix do not significantly benefit from side chain contribution subtraction, in some cases prediction is even degraded. We show that co-linearity between secondary structure content and amino acid composition is not a main limitation for improving secondary structure prediction. We also show that, even though based on different criteria, secondary structures defined by DSSP and XTLSSTR both arrive at the same conclusion: only the β-sheet structure clearly benefits from side chain subtraction. It must be concluded that side chain contribution subtraction benefit for the evaluation of other secondary structure contents is limited by the very rough description of side chain absorbance which does not take into account the variations related to their environment. The study was performed on a large protein set. To deal with the large number of proteins present, we worked on protein microarrays deposited on BaF2 slides and FTIR spectra were acquired with an imaging system.

Permut - Spatially Resolved Hydration Entropies from Atomistic Simulations

The thermodynamics of the first few surface water layers profoundly affects the overall free energy change in protein folding and unfolding, bilayer self-assembly, and ligand binding, and many other biomolecular processes. Whereas the enthalpy contribution can, in principle, be directly calculated from atomistic molecular dynamics simulations, the solvation shell entropy calculation suffers particularly from the poor sampling of the high dimensional configuration space. To furthermore understand, e.g., the effects of water molecules inside a binding pocket or the contribution of individual side chains to the hydration entropy, spatial resolution is required. To address these problems, we developed Per|Mut, a new method to calculate spatially resolved hydration entropies. In Per|Mut, we first employ a permutation reduction, which exploits the permutation symmetry of the identical water molecules to alleviate the sampling problem by the Gibbs factor N! without changing the physics. Next, we use a third-order mutual information expansion to obtain the hydration shell entropy. The spatial resolution and the decomposition into one- two- and three-body correlations of the translational and rotational entropies, as well as of the translation-rotation correlation, allows for an easy physical interpretation and characterization of the resolved entropy changes. In test-simulations, Per|Mut yielded accurate entropies for an argon gas test system and solvated alkanes. Applied to the solvation statistical mechanics of hydrated octanol and the protein Crambin, Per|Mut revealed the local effects of individual chemical groups or side chains on the solvation entropy.

Contributions of primary coordination ligands and importance of outer sphere interactions in UFsc, a de novo designed protein with high affinity for metal ions

Metalloproteins constitute nearly half of all proteins and catalyze some of the most complex chemical reactions. Recently, we reported a design of 4G-UFsc (Uno Ferro single chain), a single chain four-helical bundle with extraordinarily high (30 pM) affinity for zinc. We evaluated the contribution of different side chains to binding of Co(II), Ni(II), Zn(II) and Mn(II) using systematic mutagenesis of the amino acids that constitute the primary metal coordination and outer spheres. The binding affinity of proteins for metals was then measured using isothermal titration calorimetry. Our results show that both primary metal coordination environment and side chains in the outer sphere of UFsc are highly sensitive to even slight changes and can be adapted to binding different 3d metals, including hard-to-tightly bind metal ions such as Mn(II). The studies on the origins of tight metal binding will guide future metalloprotein design efforts.

Dueling Backbones: Comparing Peptoid and Peptide Analogues of a Mussel Adhesive Protein

Ensembles of amino acid side chains often dominate the interfacial interactions of intrinsically disordered proteins; however, backbone contributions are far from negligible. Using a combination of nanoscale force measurements and molecular dynamics simulations, we demonstrated with analogous mussel-mimetic adhesive peptides and peptoids both 34 residues long that highly divergent adhesive/cohesive outcomes can be achieved on mica surfaces by altering backbone chemistry only. The Phe, Tyr, and Dopa containing peptoid variants used in this study deposited as dehydrated and incompressible films that facilitated analysis of peptoid side chain contributions to adhesion and cohesion. For example, whereas Phe and Dopa peptoids exhibited similar cohesion, Dopa peptoids were ∼3 times more adhesive than Phe peptoids on mica. Compared with the peptides, Phe peptoid achieved only ∼20% of Phe containing peptide adhesion, but the Dopa peptoids were >2-fold more adhesive than the Dopa peptides. Cation−π interactions accounted for some but not all of the cohesive interactions. Our results were corroborated by molecular dynamics simulations and highlight the importance of backbone chemistry and the potential of peptoids or peptoid/peptide hybrids as wet adhesives and primers.

ADD Force Field for Sugars and Polyols: Predicting the Additivity of Protein-Osmolyte Interaction.

The protein–osmolyte interaction has been shown experimentally to follow an additive construct, where the individual osmolyte–backbone and osmolyte–side-chain interactions contribute to the overall conformational stability of proteins. Here, we computationally reconstruct this additive relation using molecular dynamics simulations, focusing on sugars and polyols, including sucrose and sorbitol, as model osmolytes. A new set of parameters (ADD) is developed for this purpose, using the individual Kirkwood–Buff integrals for sugar–backbone and sugar–side-chain interactions as target experimental data. We show that the ADD parameters can reproduce the additivity of protein–sugar interactions and correctly predict sucrose and sorbitol self-association and their interaction with water. The accurate description of the separate osmolyte–backbone and osmolyte–side-chain contributions also automatically translates into a good prediction of preferential exclusion from the surface of ribonuclease A and α-chymotrypsinogen A. The description of sugar polarity is improved compared to previous force fields, resulting in closer agreement with the experimental data and better compatibility with charged groups, such as the guanidinium moiety. The ADD parameters are developed in combination with the CHARMM36m force field for proteins, but good compatibility is also observed with the AMBER 99SB-ILDN and the OPLS-AA force fields. Overall, exploiting the additivity of protein–osmolyte interactions is a promising approach for the development of new force fields.

Relationship between activity and stability: Design and characterization of stable variants of human frataxin

The relationships between conformational dynamics, stability and protein function are not obvious. Frataxin (FXN) is an essential protein that forms part of a supercomplex dedicated to the iron-sulfur (Fe–S) cluster assembly within the mitochondrial matrix. In humans, the loss of FXN expression or a decrease in its functionality results in Friedreich's Ataxia, a cardio-neurodegenerative disease. Recently, the way in which FXN interacts with the rest of the subunits of the supercomplex was uncovered. This opens a window to explore relationships between structural dynamics and function. In this study, we prepared a set of FXN variants spanning a broad range of conformational stabilities. Variants S160I, S160M and A204R were more stable than the wild-type and showed similar biological activity. Additionally, we prepared SILCAR, a variant that combines S160I, L203C and A204R mutations. SILCAR was 2.4 kcal mol−1 more stable and equally active. Some of the variants were significantly more resistant to proteolysis than the wild-type FXN. SILCAR showed the highest resistance, suggesting a more rigid structure. It was corroborated by means of molecular dynamics simulations. Relaxation dispersion NMR experiments comparing SILCAR and wild-type variants suggested similar internal motions in the microsecond to millisecond timescale. Instead, variant S157I showed higher denaturation resistance but a significant lower function, similarly to that observed for the FRDA variant N146K. We concluded that the contribution of particular side chains to the conformational stability of FXN might be highly subordinated to their impact on both the protein function and the stability of the functional supercomplex.

Volumetric, acoustic and infrared spectroscopic study of amino acids in aqueous solutions of pyrrolidinium based ionic liquid, 1–butyl–1–methyl pyrrolidinium bromide

The effect of ionic liquid, 1–butyl–1–methyl pyrrolidinium bromide [C4C1Pyrr]Br on the thermodynamic properties of glycine, L–alanine and L–valine in aqueous solution was studied at the temperature range of 288.15–308.15 K and experimental pressure of p = 0.1 MPa by density, ultrasonic speed and FTIR spectroscopic techniques. Apparent molar volumes, standard partial molar volumes, standard partial molar volumes of transfer and hydration number are evaluated using density data. Side chain contributions and temperature dependence of partial molar volume have been calculated. Apparent molar isentropic compression, standard partial molar isentropic compression, standard partial molar isentropic compression of transfer and hydration number have been calculated using ultrasonic speed data. Fourier transform infrared (FTIR) spectral studies have also have been carried out for the present mixture. These volumetric, acoustic and spectroscopic studies can help to understand the mixing effects and other complex biological processes between amino acids and ionic liquid aqueous solution.

SESCA: Predicting Circular Dichroism Spectra from Protein Molecular Structures.

Circular dichroism (CD) spectroscopy is a highly sensitive but low-resolution technique to study the structure of proteins. Combined with molecular modeling or other complementary techniques, CD spectroscopy can provide essential information at higher resolution. To this end, we introduce a new computational method to calculate the electronic circular dichroism spectra of proteins from a structural model or ensemble using the average secondary structure composition and a precalculated set of basis spectra. The method is designed for model validation to estimate the error of a given protein structural model based on the measured CD spectrum. We compared the predictive power of our method to that of existing algorithms, namely, DichroCalc and PDB2CD, and found that it predicts CD spectra more accurately. Our results indicate that the derived basis sets are robust to both experimental errors in the reference spectra and the choice of the secondary structure classification algorithm. For over 80% of the globular reference proteins, our basis sets accurately predict the experimental spectrum solely from their secondary structure composition. For the remaining 20%, correcting for intensity normalization considerably improves the prediction power. Additionally, we show that the predictions for short peptides and an example complex of intrinsically disordered proteins strongly benefit from accounting for side-chain contributions and structural flexibility.

Replacing voltage sensor arginines with citrulline provides mechanistic insight into charge versus shape.

Voltage-dependent activation of voltage-gated cation channels results from the outward movement of arginine-bearing helices within proteinaceous voltage sensors. The voltage-sensing residues in potassium channels have been extensively characterized, but current functional approaches do not allow a distinction between the electrostatic and steric contributions of the arginine side chain. Here we use chemical misacylation and in vivo nonsense suppression to encode citrulline, a neutral and nearly isosteric analogue of arginine, into the voltage sensor of the Shaker potassium channel. We functionally characterize the engineered channels and compare them with those bearing conventional mutations at the same positions. We observe effects on both voltage sensitivity and gating kinetics, enabling dissection of the roles of residue structure versus positive charge in channel function. In some positions, substitution with citrulline causes mild effects on channel activation compared with natural mutations. In contrast, substitution of the fourth S4 arginine with citrulline causes substantial changes in the conductance-voltage relationship and the kinetics of the channel, which suggests that a positive charge is required at this position for efficient voltage sensor deactivation and channel closure. The encoding of citrulline is expected to enable enhanced precision for the study of arginine residues located in crowded transmembrane environments in other membrane proteins. In addition, the method may facilitate the study of citrullination in vivo.