Ionizable Amino Acid Side Chains Research Articles

Computational study of peptide plane stacking with polar and ionizable amino acid side chains.

Parallel and T-shaped stacking interactions of the peptide plane with polar and ionizable amino acid side chains (including aspartic/glutamic acid, asparagine/glutamine, and arginine) are investigated using the quantum mechanical MP2 and CCSD computational methods. It is found that the electrostatic interaction plays an essential role in determining the optimal stacking configurations for all investigated stacking models. For certain complexes, the dispersion interaction also contributes considerably to stacking. In the gas phase, the stacking interaction of the charged system is stronger than that of the neutral system, and T-shaped stacking is generally more preferred than parallel stacking, with the stacking energy in the range of -4 to -18 kcal/mol. The solvation effect overall weakens stacking, especially for the charged system and the T-shaped stacking configurations. In water, the interaction energies of different stacking models are comparable.

Interactions between Ionizable Amino Acid Side Chains at a Lipid Bilayer–Water Interface

Potentials of mean force (PMF) between ionizable amino acid side chains (Arg, Lys, His, Glu) in the headgroup area of a palmitoyl oleoyl phosphatidylcholine lipid bilayer were obtained from all-atom molecular dynamics simulations and the adaptive biasing force method. Simulations in bulk water were also performed for comparison. Side chains were constrained in collinear, stacking, and orthogonal (T-shaped) orientations. The most structured and attractive PMFs were observed for hydrogen-bonded side chains. Contact minima occurred at a distance of 2.6-3.1 Å between selected atoms or centers of mass with the most attractive interaction (-9.6 kcal/mol) observed between Arg(+) and Glu(-). Hydrogen bonds play a significant role in stabilizing these interactions. Interactions between like charged side chains can also be very attractive if the charges are screened by surrounding molecules or groups (e.g., the PMF value at the contact minimum for Arg(+)···Arg(+) is -7.6 kcal/mol). Like charged side chains can have contact minima as close as 3.6 Å. The PMFs depend strongly on the relative orientation of the side chains. In agreement with experimental studies and other simulations, we found the stacking arrangement of like charged side chains to be the most favorable orientation. Interaction energies and Lennard-Jones energies between side chains, headgroups, and water molecules were analyzed in order to rationalize the observed PMFs and their dependence on orientation. In general, the results cannot be explained by simple dielectric arguments.

A possible CO2 conducting and concentrating mechanism in plant stomata SLAC1 channel.

Background The plant SLAC1 is a slow anion channel in the membrane of stomatal guard cells, which controls the turgor pressure in the aperture-defining guard cells, thereby regulating the exchange of water vapour and photosynthetic gases in response to environmental signals such as drought, high levels of carbon dioxide, and bacterial invasion. Recent study demonstrated that bicarbonate is a small-molecule activator of SLAC1. Higher CO2 and HCO3 – concentration activates S-type anion channel currents in wild-type Arabidopsis guard cells. Based on the SLAC1 structure a theoretical model is derived to illustrate the activation of bicarbonate to SLAC1 channel. Meanwhile a possible CO2 conducting and concentrating mechanism of the SLAC1 is proposed.MethodologyThe homology structure of Arabidopsis thaliana SLAC1 (AtSLAC1) provides the structural basis for study of the conducting and concentrating mechanism of carbon dioxide in SLAC1 channels. The pKa values of ionizable amino acid side chains in AtSLAC1 are calculated using software PROPKA3.0, and the concentration of CO2 and anion HCO3 – are computed based on the chemical equilibrium theory.ConclusionsThe AtSLAC1 is modeled as a five-region channel with different pH values. The top and bottom layers of channel are the alkaline residue-dominated regions, and in the middle of channel there is the acidic region surrounding acidic residues His332. The CO2 concentration is enhanced around 104 times by the pH difference between these regions, and CO2 is stored in the hydrophobic region, which is a CO2 pool. The pH driven CO2 conduction from outside to inside balances the back electromotive force and maintain the influx of anions (e.g. Cl– and NO3 –) from inside to outside. SLAC1 may be a pathway providing CO2 for photosynthesis in the guard cells.

An aqueous H+ permeation pathway in the voltage-gated proton channel Hv1.

Hv1 voltage-gated proton channels mediate rapid and selective transmembrane H(+) flux and are gated by both voltage and pH gradients. Selective H(+) transfer in membrane proteins is commonly achieved by Grotthuss proton 'hopping' in chains of ionizable amino acid side chains and intraprotein water molecules. To identify whether ionizable residues are required for proton permeation in Hv1, we neutralized candidate residues and measured expressed voltage-gated H(+) currents. Unexpectedly, charge neutralization was insufficient to abrogate either the Hv1 conductance or coupling of pH gradient and voltage-dependent activation. Molecular dynamics simulations revealed water molecules in the central crevice of Hv1 model structures but not in homologous voltage-sensor domain (VSD) structures. Our results indicate that Hv1 most likely forms an internal water wire for selective proton transfer and that interactions between water molecules and S4 arginines may underlie coupling between voltage- and pH-gradient sensing.

Potential of Mean Force Between Ionizable Amino Acid Side Chains in Lipid Bilayer

Potentials of mean force (PMF) between ionizable amino acid side chains (Arg, Lys, His, Glu/Asp) in different protonation states in palmitoyl oleoyl phosphatidylcholine lipid bilayer were obtained from all-atom explicit solvent molecular dynamics simulations and the adaptive biasing force approach available with NAMD.1,2 Side chains (SC) were constrained in different orientations: collinear, stacked and T-shaped and placed into the bilayer interface. The most structured PMFs were observed for unlike-charged ions or pairs with neutral SCs in collinear orientation. Contact pairs (CP) occurred at a distance of 2.6-3.1 A with the strongest interaction of −9.6 kcal/mol between Arg+ and Glu- ions. Like-charged SCs in this orientation displayed less stable contact minima at greater distances or solvent separated minima. All pairs in stacking approach showed similar, well-structured PMF profiles with CPs at ∼3.8 A. The strongest interaction between like-charged pairs was observed for stacked arginines. Like-charged pairs constrained in T-shaped geometry mostly displayed slightly stable solvent separated minima. A relationship between water and phosphate coordination numbers, contact pair minima and free energy barriers was found. There is also dependence of PMF shapes on H-bonding between amino acids. Generally, interactions between ionizable SCs are more attractive and the PMFs are more structured in a lipid bilayer than in water.31 Darve, E.; Pohorille, A. J. Chem. Phys. 2001, 115, 9169.2 Darve, E.; Wilson, M.; Pohorille, A. Mol. Simul.2002, 28, 113.3 Masunov, A.; Lazaridis, T.J. Am. Chem. Soc.2003125, 1722.

Assessing Atomistic and Coarse-Grained Force Fields for Protein−Lipid Interactions: the Formidable Challenge of an Ionizable Side Chain in a Membrane

Ionizable amino acid side chains play important roles in membrane protein structure and function, including the activation of voltage-gated ion channels, where it has been previously suggested that charged side chains may move through the hydrocarbon core of the membrane. However, all-atom molecular dynamics simulations have demonstrated large free energy barriers for such lipid-exposed motions. These simulations have also revealed that the membrane will deform due to the presence of a charged side chain, leading to a complex solvation microenvironment for which empirical force fields were not specifically parametrized. We have tested the ability of the all-atom CHARMM, Drude polarizable CHARMM, and a recent implementation of a coarse-grained force field to measure the thermodynamics of arginine-membrane interactions as a function of protonation state. We have employed model systems to attempt to match experimental bulk partitioning and quantum mechanical interactions within the membrane and found that free energy profiles from nonpolarizable and polarizable CHARMM simulations are accurate to within 1-2 kcal/mol. In contrast, the coarse-grained simulations failed to reproduce the same membrane deformations, exhibit interactions that are an order of magnitude too small, and thus, have incorrect free energy profiles. These results illustrate the need for careful parametrization of coarse-grained force fields and demonstrate the utility of atomistic molecular dynamics for providing quantitative thermodynamic and mechanistic analysis of protein-lipid interactions.

A neutron crystallographic analysis of a cubic porcine insulin at pD 6.6

The p K a values of ionizable amino acid side chains are tabulated in standard textbooks. However, whether a certain amino acid side chain in a protein is charged or not cannot be estimated from standard pH values measured from protein solutions. Protonation and deprotonation of various ionizable amino acid residues were observed by a neutron diffraction experiment and discussed on the basis of the charged states estimated by the p K a values of the amino acid residues. The neutron diffraction study has been carried out at 2.7 Å resolution on cubic porcine insulin at pD 6.6 using the BIX-4 single crystal diffractometer at the JRR-3 reactor of the Japan Atomic Energy Agency. For the present work, a large single crystal of 2.7 mm 3 (=2.0 × 1.7 × 0.8 mm) was obtained by dialysis. The structure refinement was carried out using the program CNS. The resulting R cryst is 21.6% and the R free is 29.1% at a resolution of 2.7 Å. In the case of His B5, both N π and N τ of an imidazole ring are protonated at pD 6.6, but at pD 9 only N π is protonated. In contrast, for His B10, both N π and N τ are protonated at pD 6.6 as well as at pD 9. The ionization states of several amino acids in porcine insulin have been obtained at pD 6.6 and they are compared with those at pD 9 obtained by neutron diffraction as well as those at pH 6.50 and 6.98 obtained by X-ray diffraction. In this manuscript, the difference between these forms will be discussed.

Proton Movement and Photointermediate Kinetics in Rhodopsin Mutants

The role of ionizable amino acid side chains in the bovine rhodopsin activation mechanism was studied in mutants E134Q, E134R/R135E, H211F, and E122Q. All mutants exhibited bathorhodopsin stability on the 30 ns to 1 micros time scale similar to that of the wild type. Lumirhodopsin decay was also similar to that of the wild type except for the H211F mutant where early decay (20 micros) to a second form of lumirhodopsin was seen, followed by formation of an extremely long-lived Meta I(480) product (34 ms), an intermediate which forms to a much reduced extent, if at all, in dodecyl maltoside suspensions of wild-type rhodopsin. A smaller amount of a similar long-lived Meta I(480) product was seen after photolysis of E122Q, but E134Q and E134R/R135Q displayed kinetics much more similar to those of the wild type under these conditions (i.e., no Meta I(480) product). These results support the idea that specific interaction of His211 and Glu122 plays a significant role in deprotonation of the retinylidene Schiff base and receptor activation. Proton uptake measurements using bromcresol purple showed that E122Q was qualitatively similar to wild-type rhodopsin, with at least one proton being released during lumirhodopsin decay per Meta I(380) intermediate formed, followed by uptake of at least two protons per rhodopsin bleached on a time scale of tens of milliseconds. Different results were obtained for H211F, E134Q, and E134R/R135E, which all released approximately two protons per rhodopsin bleached. These results show that several ionizable groups besides the Schiff base imine are affected by the structural changes involved in rhodopsin activation. At least two proton uptake groups and probably at least one proton release group in addition to the Schiff base are present in rhodopsin.

Stabilization of Internal Charges in a Protein: Water Penetration or Conformational Change?

The ionizable amino acid side chains of proteins are usually located at the surface. However, in some proteins an ionizable group is embedded in an apolar internal region. Such buried ionizable groups destabilize the protein and may trigger conformational changes in response to pH variations. Because of the prohibitive energetic cost of transferring a charged group from water to an apolar medium, other stabilizing factors must be invoked, such as ionization-induced water penetration or structural changes. To examine the role of water penetration, we have measured the 17O and 2H magnetic relaxation dispersions (MRD) for the V66E and V66K mutants of staphylococcal nuclease, where glutamic acid and lysine residues are buried in predominantly apolar environments. At neutral pH, where these residues are uncharged, we find no evidence of buried water molecules near the mutation site. This contrasts with a previous cryogenic crystal structure of the V66E mutant, but is consistent with the room-temperature crystal structure reported here. MRD measurements at different pH values show that ionization of Glu-66 or Lys-66 is not accompanied by penetration of long-lived water molecules. On the other hand, the MRD data are consistent with a local conformational change in response to ionization of the internal residues.

Electrostatic influence on ion transport through the alphaHL channel.

The current-voltage relationship of a single Staphylococcus aureus alpha-hemolysin (alphaHL) channel is nonlinear, rectifying, and depends on the bulk pH and the ionic strength. The data are described qualitatively by a simple one-dimensional Nernst-Planck analysis in which the fixed charges inside and near the pore's entrances affect the transport of ions through the channel. The distances of these fixed charges from one of the channel's entrances are obtained from the channel's crystal structure. The model demonstrates that rectification of monovalent ion flow through the alphaHL channel can be related to the asymmetry in the location of the ionizable amino acid side chains.

Potentials of mean force between ionizable amino acid side chains in water.

Potentials of mean force (PMF) between all possible ionizable amino acid side chain pairs in various protonation states were calculated using explicit solvent molecular dynamics simulations with umbrella sampling and the weighted histogram analysis method. The side chains were constrained in various orientations inside a spherical cluster of 200 water molecules. Beglov and Roux's Spherical Solvent Boundary Potential was used to account for the solvent outside this sphere. This approach was first validated by calculating PMFs between monatomic ions (K(+), Na(+), Cl(-)) and comparing them to results from the literature and results obtained using Ewald summation. The strongest interaction (-4.5 kcal/mol) was found for the coaxial Arg(+).Glu(-) pair. Many like-charged side chains display a remarkable lack of repulsion, and occasionally a weak attraction. The PMFs are compared to effective energy curves obtained with common implicit solvation models, namely Generalized Born (GB), EEF1, and uniform dielectric of 80. Overall, the EEF1 curves are too attractive, whereas the GB curves in most cases match the minima of the PMF curves quite well. The uniform dielectric model, despite some fortuitous successes, is grossly inadequate.

Patterns in ionizable side chain interactions in protein structures.

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In alpha-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of alpha-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 angstrom. The sch:mch interactions involving ionizable side chains that belong either to beta-strands or to the central part of alpha-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.

What are enzyme structures telling us?

Most globular proteins are waxy inside and soapy outside. Their compact structures are stabilised by the hydrophobic effect which is mainly entropic, and by hydrogen bonds and dispersion forces which are mainly enthalpic. Structurally homologous proteins tend to share a common set of internal sites from which polar residues are excluded, even if they share little sequence homology. Electrostatic effects are dominant in enzyme catalysis. The active sites of enzymes are generally buried in clefts or cavities where dipoles tend to be oriented so as to optimise the pKas of ionizable amino acid side chains for catalysis; substrates are clamped close by to maximise electrostatic interactions. The activities of many enzymes have long been known to be controlled by allosteric effects or by induced fit. Some serine proteinase inhibitors have evolved yet another control mechanism: this is a spring-loaded safety catch that makes them revert to their latent, stable, inactive form unless the catch is kept in the 'loaded' position by another molecule.