pK Shifts Research Articles

Forced Gating Motions by a Substituted Titratable Side Chain at the Bundle Crossing of a Potassium Channel

Numerous inwardly rectifying potassium (Kir) channels possess an aromatic residue in the helix bundle crossing region, forming the narrowest pore constriction in crystal structures. However, the role of the Kir channel bundle crossing as a functional gate remains uncertain. We report a unique phenotype of Kir6.2 channels mutated to encode glutamate at this position (F168E). Despite a prediction of four glutamates in close proximity, Kir6.2(F168E) channels are predominantly closed at physiological pH, whereas alkalization causes rapid and reversible channel activation. These findings suggest that F168E glutamates are uncharged at physiological pH but become deprotonated at alkaline pH, forcing channel opening due to mutual repulsion of nearby negatively charged side chains. The potassium channel pore scaffold likely brings these glutamates close together, causing a significant pK(a) shift relative to the free side chain (as seen in the KcsA selectivity filter). Alkalization also shifts the apparent ATP sensitivity of the channel, indicating that forced motion of the bundle crossing is coupled to the ATP-binding site and may resemble conformational changes involved in wild-type Kir6.2 gating. The study demonstrates a novel mechanism for engineering extrinsic control of channel gating by pH and shows that conformational changes in the bundle crossing region are involved in ligand-dependent gating of Kir channels.

Influence of pK a Shifts on the Calculated Dipole Moments of Proteins

The protein dipole moment is a low-resolution parameter that characterizes the second-order charge organization of a biomolecule. Theoretical approaches to calculate protein dipole moments rely on pK(a) values, which are either computed individually for each ionizable residue or obtained from model compounds. The influence of pK(a) shifts are evaluated first by comparing calculated and measured dipole moments of β-lactoglobulin. Second, calculations are made on a dataset of 66 proteins from the Protein Data Bank, and average differences are determined between dipole moments calculated with model pK(a)s, pK(a)s derived using a Poisson-Boltzmann approach, and empirically-calculated pK(a)s. Dipole moment predictions that neglect pK(a) shifts are consistently larger than predictions in which they are included. The importance of pK(a) shifts are observed to vary with protein size, internal permittivity, and solution pH.

Nonlocal continuum electrostatic theory predicts surprisingly small energetic penalties for charge burial in proteins

We study the energetics of burying charges, ion pairs, and ionizable groups in a simple protein model using nonlocal continuum electrostatics. Our primary finding is that the nonlocal response leads to markedly reduced solvent screening, comparable to the use of application-specific protein dielectric constants. Employing the same parameters as used in other nonlocal studies, we find that for a sphere of radius 13.4 Å containing a single +1e charge, the nonlocal solvation free energy varies less than 18 kcal/mol as the charge moves from the surface to the center, whereas the difference in the local Poisson model is ∼35 kcal/mol. Because an ion pair (salt bridge) generates a comparatively more rapidly varying Coulomb potential, energetics for salt bridges are even more significantly reduced in the nonlocal model. By varying the central parameter in nonlocal theory, which is an effective length scale associated with correlations between solvent molecules, nonlocal-model energetics can be varied from the standard local results to essentially zero; however, the existence of the reduction in charge-burial penalties is quite robust to variations in the protein dielectric constant and the correlation length. Finally, as a simple exploratory test of the implications of nonlocal response, we calculate glutamate pK(a) shifts and find that using standard protein parameters (ε(protein) = 2-4), nonlocal results match local-model predictions with much higher dielectric constants. Nonlocality may, therefore, be one factor in resolving discrepancies between measured protein dielectric constants and the model parameters often used to match titration experiments. Nonlocal models may hold significant promise to deepen our understanding of macromolecular electrostatics without substantially increasing computational complexity.

Proton storage site in bacteriorhodopsin: new insights from quantum mechanics/molecular mechanics simulations of microscopic pK(a) and infrared spectra.

Identifying the group that acts as the proton storage/loading site is a challenging but important problem for understanding the mechanism of proton pumping in biomolecular proton pumps, such as bacteriorhodopsin (bR) and cytochrome c oxidase. Recent experimental studies of bR propelled the idea that the proton storage/release group (PRG) in bR is not an amino acid but a water cluster embedded in the protein. We argue that this idea is at odds with our knowledge of protein electrostatics, since invoking the water cluster as the PRG would require the protein to raise the pK(a) of a hydronium by almost 11 pK(a) units, which is difficult considering known cases of pK(a) shifts in proteins. Our recent quantum mechanics/molecular mechanics (QM/MM) simulations suggested an alternative "intermolecular proton bond" model in which the stored proton is shared between two conserved Glu residues (194 and 204). Here we show that this model leads to microscopic pK(a) values consistent with available experimental data and the functional requirement of a PRG. Extensive QM/MM simulations also show that, independent of a number of technical issues, such as the influence of QM region size, starting X-ray structure, and nuclear quantum effects, the "intermolecular proton bond" model is qualitatively consistent with available spectroscopic data. Potential of mean force calculations show explicitly that the stored proton strongly prefers the pair of Glu residues over the water cluster. The results and analyses help highlight the importance of considering protein electrostatics and provide arguments for why the "intermolecular proton bond" model is likely applicable to the PRG in biomolecular proton pumps in general.

Exploring conformational changes coupled to ionization states using a hybrid Rosetta‐MCCE protocol

A hybrid protocol combining Rosetta fullatom refinement and Multi-Conformation Continuum Electrostatics (MCCE) to estimate pK(a) is applied to the blind prediction of 94 mutated residues in Staphylococcal nuclease (SNase), as part of the pK(a)-cooperative benchmark test. The standard MCCE method is limited to sidechain conformational changes. The Rosetta refinement protocol is used to add the backbone conformational changes in pK(a) calculations. The non-electrostatic energy component from Rosetta and the electrostatic energy from MCCE are combined to weight the calculated ionization states. Of 63 measured pK(a)s, the root mean squared deviation (RMSD) between the calculated pK(a)s and the measured values is 4.3, showing an improvement compared to the RMSD of 6.6 in the standard MCCE calculations, using a low protein dielectric constant of 4. The breakdown of pK(a) shift from the solution values (ΔpK(a)) shows that the desolvation energy contributes the most in the standard MCCE calculations. Lowering desolvation penalties and optimizing electrostatic interactions with the Rosetta/MCCE protocol reduces the ΔpK(a) to favor the charged states. Analysis also showed that the Rosetta/MCCE protocol samples conformations with pK(a)s close to the solution values. The question remains whether the correct conformational changes coupled to the ionization changes are found here. Nevertheless, a challenge emerges to accurately estimate the reorganization energy, which is not directly measured from the electrostatic environment of the site of interest. Possible improvements to the protocol are also discussed.

Toward accurate prediction of pKa values for internal protein residues: The importance of conformational relaxation and desolvation energy

Proton uptake or release controls many important biological processes, such as energy transduction, virus replication, and catalysis. Accurate pK(a) prediction informs about proton pathways, thereby revealing detailed acid-base mechanisms. Physics-based methods in the framework of molecular dynamics simulations not only offer pK(a) predictions but also inform about the physical origins of pK(a) shifts and provide details of ionization-induced conformational relaxation and large-scale transitions. One such method is the recently developed continuous constant pH molecular dynamics (CPHMD) method, which has been shown to be an accurate and robust pK(a) prediction tool for naturally occurring titratable residues. To further examine the accuracy and limitations of CPHMD, we blindly predicted the pK(a) values for 87 titratable residues introduced in various hydrophobic regions of staphylococcal nuclease and variants. The predictions gave a root-mean-square deviation of 1.69 pK units from experiment, and there were only two pK(a)'s with errors greater than 3.5 pK units. Analysis of the conformational fluctuation of titrating side-chains in the context of the errors of calculated pK(a) values indicate that explicit treatment of conformational flexibility and the associated dielectric relaxation gives CPHMD a distinct advantage. Analysis of the sources of errors suggests that more accurate pK(a) predictions can be obtained for the most deeply buried residues by improving the accuracy in calculating desolvation energies. Furthermore, it is found that the generalized Born implicit-solvent model underlying the current CPHMD implementation slightly distorts the local conformational environment such that the inclusion of an explicit-solvent representation may offer improvement of accuracy.

Proton Transfer between Tryptophan and Ionic Liquid Solvents Studied with Molecular Dynamics Simulations

The reaction free energies and associated pK(a) values for proton transfer from positively charged tryptophan (HTrp(+)) to the two pure ionic liquids (ILs) BMIM-PF6 and BMIM-BF4 are derived from molecular simulations. IL solvation effects are examined with molecular dynamics simulations together with an empirical force field in which the average charge distribution in the actual IL is taken into account. A combination of molecular mechanical and quantum mechanical description (QM/MM) is used to examine the protonation of the anion constituents of the ILs. A dissociation of the protonated anions is observed into hydrogen fluoride and BF3 or PF5. Finally, pK(a) values of 16.5 and 21.5 in BMIM-BF4 and BMIM-PF6, respectively, are found for proton transfer from HTrp(+) to PF6(-) and BF4(-) anions, which indicates that a deprotonation of HTrp(+) is highly unfavorable compared to aqueous solutions. An examination of the contributions to the reaction free energies demonstrates that a deprotonation of tryptophan is impeded because two ions need to be annihilated for the reaction to occur: HTrp(+) and an anion. While the solvation effects induced by the two ILs are similar, the low proton acceptance of PF6(-) anions leads to the larger pK(a) value in BMIM-PF6. Also, estimates suggest that IL-induced pK(a) shifts are comparably small in proton transfer reactions where the total number of ions remains unchanged. For the first time, pK(a) values of acids were determined computationally in ILs. The obtained results elucidate the role of solvation effects on proton transfer between amino acids and ILs and improve our understanding of the observed pH memory of proteins that are solvated in ILs.

Conformational relaxation and water penetration coupled to ionization of internal groups in proteins.

Molecular dynamics simulations were used to examine the effects of ionization of internal groups on the structures of eighteen variants of staphylococcal nuclease (SNase) with internal Lys, Asp, or Glu. In most cases the RMSD values of internal ionizable side chains were larger when the ionizable moieties were charged than when they were neutral. Calculations of solvent-accessible surface area showed that the internal ionizable side chains were buried in the protein interior when they were neutral and moved toward crevices and toward the protein-water interface when they were charged. The only exceptions are Lys-36, Lys-62, and Lys-103, which remained buried even after charging. With the exception of Lys-38, the number of internal water molecules surrounding the ionizable group increased upon charging: the average number of water oxygen atoms within the first hydration shell increased by 1.7 for Lys residues, by 5.2 for Asp residues, and by 3.2 for Glu residues. The polarity of the microenvironment of the ionizable group also increased when the groups were charged: the average number of polar atoms of any kind within the first hydration shell increased by 2.7 for Lys residues, by 4.8 for Asp residues, and by 4.0 for Glu residues. An unexpected correlation was observed between the absolute value of the shifts in pK(a) values measured experimentally, and several parameters of structural relaxation: the net difference in the polarity of the microenvironment of the charged and neutral forms of the ionizable groups, the net difference in hydration of the charged and neutral forms of the ionizable groups, and the difference in RMSD values of the charged and neutral forms of the ionizable groups. The effects of ionization of internal groups on the conformation of the backbone were noticeable but mostly small and localized to the area immediately next to the internal ionizable moiety. Some variants did exhibit local unfolding.

Equilibrium Unfolding Studies of Monellin: The Double-Chain Variant Appears To Be More Stable Than the Single-Chain Variant

To improve our understanding of the contributions of different stabilizing interactions to protein stability, including that of residual structure in the unfolded state, the small sweet protein monellin has been studied in both its two variant forms, the naturally occurring double-chain variant (dcMN) and the artificially created single-chain variant (scMN). Equilibrium guanidine hydrochloride-induced unfolding studies at pH 7 show that the standard free energy of unfolding, ΔG°(U), of dcMN to unfolded chains A and B and its dependence on guanidine hydrochloride (GdnHCl) concentration are both independent of protein concentration, while the midpoint of unfolding has an exponential dependence on protein concentration. Hence, the unfolding of dcMN like that of scMN can be described as two-state unfolding. The free energy of dissociation, ΔG°(d), of the two free chains, A and B, from dcMN, as measured by equilibrium binding studies, is significantly lower than ΔG°(U), apparently because of the presence of residual structure in free chain B. The value of ΔG°(U), at the standard concentration of 1 M, is found to be ∼5.5 kcal mol(-1) higher for dcMN than for scMN in the range from pH 4 to 9, over which unfolding appears to be two-state. Hence, dcMN appears to be more stable than scMN. It seems that unfolded scMN is stabilized by residual structure that is absent in unfolded dcMN and/or that native scMN is destabilized by strain that is relieved in native dcMN. The value of ΔG°(U) for both protein variants decreases with an increase in pH from 4 to 9, apparently because of the thermodynamic coupling of unfolding to the protonation of a buried carboxylate side chain whose pK(a) shifts from 4.5 in the unfolded state to 9 in the native state. Finally, it is shown that although the thermodynamic stabilities of dcMN and scMN are very different, their kinetic stabilities with respect to unfolding in GdnHCl are very similar.

Method for Estimating the Internal Permittivity of Proteins Using Dielectric Spectroscopy

Protein charge organization is dependent on the low-permittivity region in the hydrophobic core of the molecule. We suggest a novel approach to estimate the dielectric constant of this region by comparing measured and simulated first- and second-order charge moments. Here, the dipole moment is measured as a function of pH using dielectric spectroscopy. The results are compared to dipole moments based on Poisson-Boltzmann estimates of pK(a) shifts calculated from structures in the Protein Data Bank. Structures are additionally refined using CHARMM molecular dynamics simulations. The best estimate for the internal permittivity is found by minimizing the root-mean-square residual between measured and predicted charge moments. Using the protein β-lactoglobulin, a core dielectric constant in the range of 6-7 is estimated.

Influence of End-Capping on the Self-Assembly of Model Amyloid Peptide Fragments

The influence of charge and aromatic stacking interactions on the self-assembly of a series of four model amyloid peptides has been examined. The four model peptides are based on the KLVFF motif from the amyloid β peptide, Aβ(16-20) extended at the N terminus with two β-alanine residues. We have studied NH(2)-βAβAKLVFF-COOH (FF), NH(2)-βAβAKLVF-COOH (F), CH(3)CONH-βAβAKLVFF-CONH(2) (CapF), and CH(3)CONH-βAβAKLVFF-CONH(2) (CapFF). The former two are uncapped (net charge +2) and differ by one hydrophobic phenylalanine residue; the latter two are the analogous capped peptides (net charge +1). The self-assembly characteristics of these peptides are remarkably different and strongly dependent on concentration. NMR shows a shift from carboxylate to carboxylic acid forms upon increasing concentration. Saturation transfer measurements of solvent molecules indicate selective involvement of phenylalanine residues in driving the self-assembly process of CapFF due presumably to the effect of aromatic stacking interactions. FTIR spectroscopy reveals β-sheet features for the two peptides containing two phenylalanine residues but not the single phenylalanine residue, pointing again to the driving force for self-assembly. Circular dichroism (CD) in dilute solution reveals the polyproline II conformation, except for F which is disordered. We discuss the relationship of this observation to the significant pH shift observed for this peptide when compared the calculated value. Atomic force microscopy and cryogenic-TEM reveals the formation of twisted fibrils for CapFF, as previously also observed for FF. The influence of salt on the self-assembly of the model β-sheet forming capped peptide CapFF was investigated by FTIR. Cryo-TEM reveals that the extent of twisting decreases with increased salt concentration, leading to the formation of flat ribbon structures. These results highlight the important role of aggregation-induced pK(a) shifts in the self-assembly of model β-sheet peptides.

Supramolecular assembly of hoechst-33258 with cucurbit[7]uril macrocycle

Molecular assemblies of potential guest molecules through non-covalent host-guest interactions have found immense use in many applied areas. In this study supramolecular interaction of a biologically important dye Hoechst-33258 (H33258) has been investigated in aqueous solutions at different pHs, in the presence of a macrocyclic host, namely, cucurbit[7]uril (CB7). The pH dependent emission behaviour of H33258 is inherently connected with its protolytic equilibria which allow the dye in different geometrical conformations. This pH dependent structural orientation is greatly affected by the complexation with CB7. The significant structural changes in the monocationic H33258 brought out by CB7 at pH 7 have been documented in the fluorescence emission and lifetime data, which are comparatively less affected in case of the dicationic form, which is prominent in dye solutions at pH 4.5. The strong ion-dipole interactions provided by the carbonyl portals of the CB7 host adequately stabilize the CB7-H33258 complex, both in 1:1 and 2:1 stoichiometries at both the pH conditions. The Job's plot method, fluorescence anisotropy, NMR measurements and geometry optimization calculations confirm the stoichiometric arrangement and are found to be tunable with the addition of metal ions. The non-covalently stabilized assembly brings out large enhancement in the fluorescence emission due to the unique structural orientation attained by H33258, which reduces the non-radiative relaxation pathways. Comparison of the spectral data of the dye at different pH conditions in the absence and presence of CB7 proposes a large upward pK(a) shift due to CB7 encapsulation, thus providing a handy tool to modulate the photophysical characteristics of the guest molecules.

NMR determination of pKa values in α‐synuclein

The intrinsically unfolded protein α-synuclein has an N-terminal domain with seven imperfect KTKEGV sequence repeats and a C-terminal domain with a large proportion of acidic residues. We characterized pK(a) values for all 26 sites in the protein that ionize below pH 7 using 2D (1) H-(15) N HSQC and 3D C(CO)NH NMR experiments. The N-terminal domain shows systematically lowered pK(a) values, suggesting weak electrostatic interactions between acidic and basic residues in the KTKEGV repeats. By contrast, the C-terminal domain shows elevated pK(a) values due to electrostatic repulsion between like charges. The effects are smaller but persist at physiological salt concentrations. For α-synuclein in the membrane-like environment of sodium dodecylsulfate (SDS) micelles, we characterized the pK(a) of His50, a residue of particular interest since it is flanked within one turn of the α-helix structure by the Parkinson's disease-linked mutants E46K and A53T. The pK(a) of His50 is raised by 1.4 pH units in the micelle-bound state. Titrations of His50 in the micelle-bound states of the E46K and A53T mutants show that the pK(a) shift is primarily due to interactions between the histidine and the sulfate groups of SDS, with electrostatic interactions between His50 and Glu46 playing a much smaller role. Our results indicate that the pK(a) values of uncomplexed α-synuclein differ significantly from random coil model peptides even though the protein is intrinsically unfolded. Due to the long-range nature of electrostatic interactions, charged residues in the α-synuclein sequence may help nucleate the folding of the protein into an α-helical structure and confer protection from misfolding.

Switchable Cucurbituril–Bipyridine Beacons

4-Aminobipyridine derivatives form strong inclusion complexes with cucurbit[6]uril, exhibiting remarkably large enhancements in fluorescence intensity and quantum yields. The remarkable complexation-induced pK(a) shift (DeltapK(a)=3.3) highlights the strong charge-dipole interaction upon binding. The reversible binding phenomenon can be used for the design of switchable beacons that can be incorporated into cascades of binding networks. This concept is demonstrated herein by three different applications: 1) a switchable fluorescent beacon for chemical sensing of transition metals and other ligands; 2) direct measurement of binding constants between cucurbit[6]uril and various nonfluorescent guest molecules; and 3) quantitative monitoring of biocatalytic reactions and determination of their kinetic parameters. The latter application is illustrated by the hydrolysis of an amide catalyzed by penicillin G acylase and by the elimination reaction of a beta-cabamoyloxy ketone catalyzed by aldolase antibody 38C2.

Determination of p Ka values of benzimidazole derivatives from mobility obtained by capillary electrophoresis

Dissociation constants of benzimidazole derivatives have been determined using capillary zone electrophoresis (CZE). Since CZE is a separation method, high purity and known concentration for the samples is not necessary because only mobilities are measured. The precision of pK(a) measurements of seven compounds is useful to observe pK(a) shifts induced by chemical variations. Some of them were compared to potentiometry and spectroscopy experiments. Good correlated pK(a) values are observed between the three analytical techniques.

Investigation of the pH‐dependent complex formation between β‐cyclodextrin and dipeptide enantiomers by capillary electrophoresis and calorimetry

The effect of pH on complex formation between beta-CD and the enantiomers of the dipeptides Ala-Phe, Ala-Tyr and Asp-PheOMe was investigated at 298.15 K by CE and calorimetry. Beta-CD displayed a higher enantioselectivity toward the protonated peptides compared to their zwitterionic forms. While stronger binding of the DD-enantiomers than the LL-stereoisomers were found by calorimetry regardless of the ionization state of the peptides, essentially equal complexation constants of the enantiomers were determined by CE for the zwitterionic species of the peptides. The reversal of the enantiomer migration order observed in CE was attributed primarily to a stereoselective complexation-induced pK(a) shift. In calorimetry, complexation of the protonated DD-enantiomers by beta-CD was accompanied by higher enthalpy and entropy changes resulting in more stable complexes compared to the LL-peptides. The enthalpy and entropy of complexation was affected by pH and peptide structure.

NMR Investigations of the Rieske Protein from Thermus thermophilus Support a Coupled Proton and Electron Transfer Mechanism

The Rieske protein component of the cytochrome bc complex contains a [2Fe−2S] cluster ligated by two cysteines and two histidines. We report here the pKa values of each of the imidazole rings of the two ligating histidines (His134 and His154) in the oxidized and reduced states of the Rieske protein from Thermus thermophilus (TtRp) as determined by NMR spectroscopy. Knowledge of these pKa values is of critical interest because of their pertinence to the mechanism of electron and proton transfer in the bifurcated Q-cycle. Although we earlier had observed the pH dependence of a 15N NMR signal from each of the two ligand histidines in oxidized TtRp (Lin, I. J.; Chen, Y.; Fee, J. A.; Song, J.; Westler, W. M.; Markley, J. L.J. Am. Chem. Soc.2006, 132, 10672−10673), the strong paramagnetism of the [2Fe−2S] cluster prevented the assignment of these signals by conventional methods. Our approach here was to take advantage of the unique histidine−leucine (His134−Leu135) sequence and to use residue-selective labeling to establish a key sequence-specific assignment, which was then extended. Analysis of the pH dependence of assigned 13C′, 13Cα, and 15Nε2 signals from the two histidine cluster ligands led to unambiguous assignment of the pKa values of oxidized and reduced TtRp. The results showed that the pKa of His134 changes from 9.1 in oxidized to ∼12.3 in reduced TtRp, whereas the pKa of His154 changes from 7.4 in oxidized to ∼12.6 in reduced TtRp. This establishes His154, which is close to the quinone when the Rieske protein is in the cytochrome b site, as the residue experiencing the remarkable redox-dependent pKa shift. Secondary structural analysis of oxidized and reduced TtRp based upon our extensive chemical shift assignments rules out a large conformational change between the oxidized and reduced states. Therefore, TtRp likely translocates between the cytochrome b and cytochrome c sites by passive diffusion. Our results are most consistent with a mechanism involving the coupled transfer of an electron and transfer of the proton across the hydrogen bond between the hydroquinone and His154 at the cytochrome b site.

Direct Observation of the pH-Dependent Equilibrium between Metarhodopsins I and II and the pH-Independent Interaction of Metarhodopsin II with Transducin C-Terminal Peptide

Bovine rhodopsin contains 11-cis-retinal as a light-absorbing chromophore that binds to a lysine residue of the apoprotein opsin via a protonated Schiff base linkage. Light isomerizes 11-cis-retinal into the all-trans form, which eventually leads to the formation of an enzymatically active state, metarhodopsin II (MII). It is widely believed that MII forms a pH-dependent equilibrium with metarhodopsin I (MI), but direct evidence for this equilibrium has not been reported. Here, we confirmed this equilibrium by direct observation of the mutual conversions of MI and MII upon changing the pH of the MI/MII mixture. We also observed a reversible binding of the synthetic peptide constituting the C-terminal 11 amino acids of the transducin alpha-subunit to MII, which resulted in change of the amounts of MI and MII in the equilibrium. Interestingly, addition of the peptide did not induce a simple pK(a) shift but rather induced an increase of the MII fraction at high pH. These results indicate that in addition to the MII that is formed from MI in a pH-dependent manner there also exists another MII, which is in equilibrium with MI in a pH-independent manner and can bind to the peptide. Therefore, there is no need for proton uptake by the protein moiety of opsin for the binding to the peptide.

Contrasting guest binding interaction of cucurbit[7-8]urils with neutral red dye: controlled exchange of multiple guests

Interactions among macrocyclic hosts and dyes/drugs have been explored extensively for their direct usage in controlled uptake and release of large number of potential drug molecules. In this paper we report the non-covalent interaction of cucurbit[8]uril macrocycle (CB8) with a biologically important dye, neutral red, by absorption and fluorescence spectroscopy. A comparative analysis with the complexation behaviour of the dye with CB7, the lower homologue of CB8, indicates contrasting guest binding behaviour with significant changes in the photophysical characteristics of the dye. While CB7 interaction leads to a 1 ratio 1 stoichiometry resulting in approximately 6 fold enhancement in the fluorescence emission of the dye, CB8 displays signatures for a 1 ratio 2 host-guest stoichiometry with drastic reduction in the fluorescence emission. Apart from the evaluation of approximately 2 unit shift in the protolytic equilibrium on complexation (pK(a) shift), the measurements with tryptophan established a selective guest exchange to favour a co-localized dimer inside the CB8 cavity. In a protein medium (BSA), the 1 ratio 2 complex was converted to a 1 ratio 1 ratio 1 CB8-NRH(+)-BSA complex. The finding that NRH(+) can be transferred from CB8 to BSA, even though the binding constant for NRH(+)-CB8 is much higher than NRH(+)-BSA, is projected for a controlled slow release of NRH(+) towards BSA. Since the release and activity of drugs can be controlled by regulating the protolytic equilibrium, the macromolecular encapsulation and release of NRH(+) demonstrated here provide information relevant to host-guest based drug delivery systems and its applications.

Conformational Sampling of Influenza Fusion Peptide in Membrane Bilayers as a Function of Termini and Protonation States

Influenza fusion peptide is critical for mediating the fusion of viral and host cell membranes during viral entry. The interaction of monomeric influenza fusion peptide with membranes is studied with replica exchange molecular dynamics simulations using a new implicit membrane model to effectively reach microsecond to millisecond time scales. The conformational sampling of the fusion peptide was studied as a function of different N- and C-termini, including an experimental construct with an additional C-terminal tag, as well as a function of protonation of acidic residues. It is found that the influenza fusion peptide mostly adopts helical structures with a pronounced kink at residues 11-13 with both N-terminal and C-terminal helices oriented mostly parallel to the membrane surface. A charged C-terminus and the presence of a charge C-terminal tag significantly alters the conformational sampling of the fusion peptide and results in more diverse conformational ensembles that include obliquely inserted N-terminal peptide structures. Protonation of acidic residues also affects the conformational sampling, however, based on pK(a) shift estimates the overall effect of pH = 5 on the conformational sampling of the influenza fusion peptide appears to be only minor.