Andrew Toyi Banyikwa1, Alan Goos1, David J Kiemle2, Michael A C Foulkes1, Mark S Braiman1. 1. Chemistry Department, Center for Science & Technology, Syracuse University, Room 1-014, Syracuse, New York 13244, United States. 2. Analytical and Technical Services, SUNY College of Environmental Science and Forestry, 123 Jahn Laboratory, Syracuse, New York 13210, United States.
Abstract
H-bonds between neutral tyrosine and arginine in nonpolar environments are modeled by small-molecule phenol/guanidine complexes. From the temperature and concentration dependence of UV spectra, a value of ΔH° = -74 ± 4 kJ mol-1 is deduced for the formation of H-bonded p-cresol/dodecylguanidine in hexane. ΔE = -71 kJ mol-1 is computed with density functional theory (in vacuo). In dimethyl sulfoxide or crystals, (p-phenolyl)alkylguanidines form head-to-tail homodimers with two strong H-bonding interactions, as evidenced by UV, IR, and NMR spectral shifts, strong IR continuum absorbance bands, and short O···N distances in X-ray crystal structures. Phenol/alkylguanidine H-bonded complexes consist of polarizable rapidly interconverting tautomers, with the proton shift from phenol to guanidine increasing with increase in the polarity of the aprotic solvent. As measured by NMR, both groups in these strongly H-bonded neutral complexes can simultaneously appear to be predominantly protonated. These systems serve as models for the hypothetical hydrogen-Bonded Uncharged (aRginine + tYrosine), or "BU(RY)", motifs in membrane proteins.
H-bonds between neutral pan class="Chemical">tyrosinen> and pan class="Chemical">arginine in nonpolar environments are modeled by small-molecule pan class="Chemical">phenol/guanidine complexes. From the temperature and concentration dependence of UV spectra, a value of ΔH° = -74 ± 4 kJ mol-1 is deduced for the formation of H-bonded p-cresol/dodecylguanidine in hexane. ΔE = -71 kJ mol-1 is computed with density functional theory (in vacuo). In dimethyl sulfoxide or crystals, (p-phenolyl)alkylguanidines form head-to-tail homodimers with two strong H-bonding interactions, as evidenced by UV, IR, and NMR spectral shifts, strong IR continuum absorbance bands, and short O···N distances in X-ray crystal structures. Phenol/alkylguanidine H-bonded complexes consist of polarizable rapidly interconverting tautomers, with the proton shift from phenol to guanidine increasing with increase in the polarity of the aprotic solvent. As measured by NMR, both groups in these strongly H-bonded neutral complexes can simultaneously appear to be predominantly protonated. These systems serve as models for the hypothetical hydrogen-Bonded Uncharged (aRginine + tYrosine), or "BU(RY)", motifs in membrane proteins.
There has been little
experimental modeling of H-bonding interactions
between the n class="Chemical">pan class="Chemical">guanidine groupn> of pan class="Chemical">arginine and the pan class="Chemical">phenolic group of
tyrosine, especially with the guanidine in a deprotonated state. Geometries
suitable for H-bonding between tyrosine and arginine were for a long
time quite uncommon in observed crystal structures of proteins, making
such H-bonds a relatively uninteresting field of study. (See section
titled “Possible New Motif for Membrane
Protein Structures”, for more details on the relative
infrequency of such groupings in the crystallographic database.)
However, this has changed recently with the publication of a number
of crystal structures of active intermediates of memn class="Chemical">pan class="Chemical">brane pn>roteins,
in which conserved pan class="Chemical">arg–pan class="Chemical">tyr pairs approach each other within
a hydrophobic region of the protein interior, in a geometry that would
permit H-bonding between them. Examples were seen first in the M state
of bacteriorhodopsin (bR), a microbial rhodopsin,[1] and subsequently in bovine opsin, a G-protein-coupled receptor
(GPCR).[2] (For an explanation of the distinction
between these the two major superfamilies of heptahelical transmembrane
(7-TM) proteins, see a recent review of microbial rhodopsins.[3]) In examples from both these superfamilies,[1,2] the conformational change that forms the active intermediate involves
significant displacement of a highly conserved arginine in the third
transmembrane helix, away from an asp or glusalt-bridge partner and
into a less-polar protein environment that includes a conserved tyr.
Evidence for arg–tyr interaction can be seen in greater
detail in more recent crystal structures of metarhodopsin II (Meta
II), with and without the C-terminal binding portion of the α
subunit of Gt, GαCT (see Figure ). These structures show that one of the
key interactions in Meta II involves the protrusion of GαCT
deep into the protein interior, allowing a direct H-bonded interaction
of the backbone of GαCT with the highly conserved arginine (R135
in rhodopsin).[4]
Figure 1
H-bonding interactions
of the conserved arginine–tyrosine
grouping in activated rhodopsin (Meta II), without (left) and with
(right) the peptide backbone of the bound C-terminal portion of transducin
(GαCT, with a backbone shown in green in the right panel). Coordinates
were obtained from data published with 3.0 Å (3XPO) and 2.85
Å (3PQR) resolutions, respectively.[4] In Meta II without GαCT (at left), one of terminal nitrogens
of Arg135 from TM3 is only 2.7 Å from the Tyr223 oxygen atom.
The other two nitrogens are not within the H-bond distance to donor
or acceptor group detected in the crystal structure. Upon binding
of Meta II with GαCT (right panel), the latter protein’s
Val347 backbone C=O group serves as an H-bond acceptor for
the other terminal nitrogen of Arg135, as well as for the intrachain
nitrogen, with N–O distances of 3.0 and 3.1 Å, respectively.
(Only the former distance is indicated by a dashed line in this figure.)
Furthermore, the C-terminal carboxylate of GαCT moves in to
become the closest anion, ∼9 Å away from Arg135. Formation
of these ionic and H-bonding interactions between GαCT and Arg135
apparently weakens the H-bonding interaction of the arginine with
tyr223, as evidenced by an increased N–O distance (3.3 Å
in the right panel).
H-bonding interactions
of the conserved pan class="Chemical">argininen>–pan class="Chemical">tyrosine
grouping in activated pan class="Gene">rhodopsin (Meta II), without (left) and with
(right) the peptide backbone of the bound C-terminal portion of transducin
(GαCT, with a backbone shown in green in the right panel). Coordinates
were obtained from data published with 3.0 Å (3XPO) and 2.85
Å (3PQR) resolutions, respectively.[4] In Meta II without GαCT (at left), one of terminal nitrogens
of Arg135 from TM3 is only 2.7 Å from the Tyr223oxygen atom.
The other two nitrogens are not within the H-bond distance to donor
or acceptor group detected in the crystal structure. Upon binding
of Meta II with GαCT (right panel), the latter protein’s
Val347 backbone C=O group serves as an H-bond acceptor for
the other terminal nitrogen of Arg135, as well as for the intrachain
nitrogen, with N–O distances of 3.0 and 3.1 Å, respectively.
(Only the former distance is indicated by a dashed line in this figure.)
Furthermore, the C-terminal carboxylate of GαCT moves in to
become the closest anion, ∼9 Å away from Arg135. Formation
of these ionic and H-bonding interactions between GαCT and Arg135
apparently weakens the H-bonding interaction of the arginine with
tyr223, as evidenced by an increased N–O distance (3.3 Å
in the right panel).
Each time a Gt binds to a Mpan class="Gene">etan> II in this fashion,
the
Gt nucleotide-binding site opens upn>, permitting the spontaneous
release of GDP, and replacement with GTP from solution. This cyclical,
catalytic G-protein activation process clearly depends on the specificity,
tightness, and reversibility of the binding of Gt to Meta
II. Many aspects of the binding interaction are conserved in other
GPCR systems.
Both pan class="Chemical">argininen> and pan class="Chemical">tyrosine side chains are generally
thought to
be protonated (cationic and neutral, respectively) in most biological
environments, including the highly conserved pan class="Chemical">arginine–tyrosine
grouping within Meta II (Figure ), as well as other activated GPCRs.[6−9] However, interpretation of these
side chains’ protonation and H-bonding state in these structures
is complicated by their unusually aprotic environment. In the published
crystal structure of Meta II in the absence of GαCT (3XPO, Figure ), numerous structural
water molecules are observable in the protein interior, but the closest
is at a distance of ∼8.5 Å from arg135.[4,5] In
fact, the environment within 5 Å of the arg135–tyr223
H-bond consists almost exclusively of hydrocarbon side chains. Such
an environment is very different from those that have previously been
used to model the arginine side chain spectroscopically.[10,11] Deprotonated arginine side chains, and H-bonded arginine–tyrosine
dyads, should be spectroscopically modeled in aprotic environments
similar to these proteins’ interiors. The current study,
and the accompanying paper,[12] provides
some key steps toward achieving this.
Scheme summarizes
three generalized proton configurations corresn class="Chemical">ponding to H-bonded
arginine–tyrosine dyads that would most likely be observable
in nonpolar environments. The first (A) involves neutral tyrosine
acting as an H-bond acceptor for protonated arginine,
which additionally participates in a stabilizing Coulombic interaction
with a nearby counterion (X–). The second (B) involves
a zwitterion comprising deprotonated tyrosine and protonated arginine,
in which the tyrosineoxygen again acts as an H-bond acceptor. The
third (C) is fully neutral, comprising neutral tyrosine serving as
an H-bond donor, as well as neutral (deprotonated)
arginine taking on the new role of H-bond acceptor (right).
Scheme 1
Three Possible Protomeric Configurations for Arg−Tyr
Dyads
in Nonpolar Environments
Scheme B,C are
the proton-limiting structures of an overall-neutral
H-bonded n class="Chemical">pair. That is, they differ only by a very small movement
of a proton within an H-bond. Therefore, if either is present, the
other is expected also to be present, in rapid equilipan class="Chemical">brium.
We modeled structures B/C computationally; and A–C exn class="Chemical">perimentally
by synthesizing model compounds that exhibited these three structures
in properly selected aprotic solvents, based on spectroscopic evidence.
Three common aprotic solvents, and one protic solvent used as a control,
were chosen to span a wide range of dielectric constants (ε):
hexane (ε = 1.8), CCl4 (ε = 2.2), dimethyl
sulfoxide (DMSO, ε = 47), and methanol (ε = 33). This
span includes the smaller ranges generally assumed for protein and
membrane interiors. These solvents were also selected for several
other important properties, including ready availability with needed
isotope variations for NMR; maintaining chemical stability of our
samples; and allowing sufficient solubility of them to obtain UV,
IR, and/or NMR spectra with suitable signal/noise ratios. Water itself
is not generally a suitable solvent for spectroscopy of any of the
deprotonated guanidine compounds because hydroxide-catalyzed hydrolysis
occurs within minutes at high pH. We were also able to crystallize
several of the model compounds in the forms corresponding to both
A and B.
Using our model systems, we demonstrate that the net-neutral
pair
(structures B and C in Scheme ) can indeed be formed with well-defined stoichiometry in
the three an class="Chemical">protic solvents examined, and structure B can be observed
in the crystalline state. The H-bond between guanidine and phenol
in such overall-neutral systems is not only unusually strong but also
highly polarizable. In the least-polar environments, both protomers
participate nearly equally.
These results allow more accurate
spectroscon class="Chemical">pic modeling of membrane-buried
arg–tyr dyads that might undergo deprotonation, that is, with
overall loss of a proton shared between arginine and a nearby tyrosine.
We term this an H-Bonded Uncharged (aRginine–tYrosine), or
“BU(RY)” dyad. We also present a specific hypothesis
for the role such a BU(RY) dyad might play in G-protein activation
by GPCRs.
Results
Computational Studies of the Phenol–Guanidine
H-Bond
The complex betweenn class="Chemical">pan class="Chemical">p-cresol and
pn>an class="Chemical">monoalkylguanidines
was modeled by using density functional theory (DFT). Rather than
dodecylamine itself, the computational modeling used the shorter alkyl
chain of ethylguanidine. The energy-optimized structure of the 1:1
complex of p-cresol with ethylguanidine in vacuo
is shown in Figure A (see Supporting Information for details).
The computed binding energy is −71 kJ mol–1, which corresponds to the energy of the optimized geometry (Figure A), relative to the
sum of the energies of the same two molecules optimized individually.
Even when a larger basis set was used (e.g., 6-311G** as compared
to 6-31G**), the computed binding energy remained unchanged to within
∼1%.
Figure 2
DFT-computed H-bonding geometries. (A) p-Cresol/ethylguanidine
complex. (B) Internally H-bonded p-phenol-dodecylguanidine.
See Supporting Information for computational
details and energy-optimized Cartesian coordinates.
DFT-computed H-bonding geometries. (A) n class="Chemical">pan class="Chemical">p-Cresol/pn>an class="Chemical">ethylguanidine
complex. (B) Internally H-bonded p-phenol-dodecylguanidine.
See Supporting Information for computational
details and energy-optimized Cartesian coordinates.
In the optimized geometry (Figure A), both H-bonds are significantly nonlinear.
For examn class="Chemical">ple,
there is a 158° O–H–N angle for the primary H-bond,
which has a 2.60 Å distance from the phenolic oxygen to the iminonitrogen. There is also a more distorted secondary (2.90 Å) H-bond,
from an aminenitrogen back to the oxygen. Formation of the H-bond(s)
resulted in significant lengthening of the O–H covalent bond
length, to 1.03 Å (from 0.98 Å in isolated p-cresol).
A similar H-bonding geometry was obtained when the
pan class="Chemical">ethylguanidinen>
and pan class="Chemical">cresol were linked together with a dodecyl chain, as shown in Figure B. One interesting
difference is the involvement of the intrachain NH group in the H-bond,
instead of only the two terminal pan class="Chemical">nitrogens. Despite the overall similar
bond lengths and angles for the H-bonding groups as in Figure A, the binding energy in Figure B is only −39
kJ mol–1. In this case, the reference state was
the optimized local energy minimum, with an extended
dodecyl chain in an all-s-trans conformation. Presumably, the H-bond
energy in Figure B,
which by itself would not be expected to be greatly different than
that in Figure A,
was significantly canceled by unfavorable strain introduced in the
alkane chain.
UV Spectroscopic Detection of Phenol–Guanidine
H-Bonding
in Hexane
The formation of pan class="Chemical">p-cresoln>/pan class="Chemical">alkylguanidine
complexes in nonpolar solvents was confirmed experimentally by UV
spectral measurements. The ∼280 nm absorption band of the pan class="Chemical">phenol
group is strongly affected by the formation of H-bonded complexes,
as has been shown previously using amines as the H-bond acceptor groups.[13,14] Similar effects are clearly seen in the UV spectra shown in Figure , which demonstrates
the interaction between p-cresol and dodecylguanidine
in hexane. From these data, we determined a value for the enthalpy
of complex formation in hexane (ΔH° =
−74 ± 4 kJ mol–1) that closely matched
the computed in vacuo ΔE = −71 kJ mol–1. (The closeness of this match is fortuitous; see Discussion.) Details of this experimental determination
are as follows.
Figure 3
UV absorption spectra demonstrating the reversible complexation
of p-cresol with equimolar dodecylguanidine in hexane.
Spectra were obtained as a function of total concentration of p-cresol (A, in which all seven spectra were measured at
58 °C; the inset shows a log–log plot of A308 vs total concentration); or as a function of temperature
(B, in which all eight spectra were measured with 0.50 mM total concentrations
of p-cresol and dodecylguanidine; the inset shows
a van’t Hoff plot). Temperatures in (B) are indicated by trace
colors: blue, 58 °C; green, 60 °C; orange, 62 °C; red,
64 °C. Two nearly superimposed traces of each color are plotted,
measured during up-then-down temperature steps. Spectral traces in
black (overlapping curves near 0 absorbance) in (A) and (B) represent
the residuals after subtracting the best-fitted linear combination
of the spectra of free p-cresol and complexed p-cresol/dodecylguanidine (C). The contribution of free p-cresol was computed as the measured extinction-coefficient
spectrum of pure p-cresol in hexane (C, black trace),
multiplied by its measured total concentration and the best-fit mole
fraction for it. These 11 mole fractions of free p-cresol, 1 for each of the 11 measured spectra, were calculated as
1 – Xi, where the values of Xi were 11 of the 512 adjustable parameters during
the least-squares fit. These Xi were the
mole fractions of the p-cresol present as a 1:1 complex
with dodecylguanidine. All 501 values of the extinction-coefficient
spectrum of this 1:1 complex (C, red trace) were also optimized during
the fit. See Supporting Information for
further details.
UV absorption sn class="Chemical">pectra demonstrating the reversible complexation
of p-cresol with equimolar dodecylguanidine in hexane.
Spectra were obtained as a function of total concentration of p-cresol (A, in which all seven spectra were measured at
58 °C; the inset shows a log–log plot of A308 vs total concentration); or as a function of temperature
(B, in which all eight spectra were measured with 0.50 mM total concentrations
of p-cresol and dodecylguanidine; the inset shows
a van’t Hoff plot). Temperatures in (B) are indicated by trace
colors: blue, 58 °C; green, 60 °C; orange, 62 °C; red,
64 °C. Two nearly superimposed traces of each color are plotted,
measured during up-then-down temperature steps. Spectral traces in
black (overlapping curves near 0 absorbance) in (A) and (B) represent
the residuals after subtracting the best-fitted linear combination
of the spectra of freep-cresol and complexed p-cresol/dodecylguanidine (C). The contribution of freep-cresol was computed as the measured extinction-coefficient
spectrum of pure p-cresol in hexane (C, black trace),
multiplied by its measured total concentration and the best-fit mole
fraction for it. These 11 mole fractions of freep-cresol, 1 for each of the 11 measured spectra, were calculated as
1 – Xi, where the values of Xi were 11 of the 512 adjustable parameters during
the least-squares fit. These Xi were the
mole fractions of the p-cresol present as a 1:1 complex
with dodecylguanidine. All 501 values of the extinction-coefficient
spectrum of this 1:1 complex (C, red trace) were also optimized during
the fit. See Supporting Information for
further details.
First, absorption sn class="Chemical">pectra
of the unassociated species were determined.
When dissolved by itself in hexane, p-cresol gives
a strong UV absorption with a maximum at 278 nm as well as vibronic
structure (Figure C, black trace). This spectrum shows no sign of concentration dependence
below ∼10 mM. Similarly, dilute dodecylguanidinefree base
in hexane shows no measurable extinction coefficient (<1 M –1 cm–1) anywhere in this 250–350
nm range. However, at dodecylguanidine concentrations above ∼5
mM, significant light scattering occurs (data not shown), indicating
self-association as reverse micelles. Unassociated dodecylguanidinium
cation was presumed also not to have any significant absorption in
hexane, based on the low absorbance (<0.001) measured for ∼10
mM solutions of dodecylguanidinium bromide in other solvents (DMSO,
CHCl3; data not shown). However, dodecylguanidine–HBr
was not sufficiently soluble in hexane to measure its extinction-coefficient
spectrum directly.
Combining pan class="Chemical">ethylguanidinen> and pan class="Chemical">p-cresol in pan class="Chemical">hexane
clearly results in an H-bonded complex. This is evident from the concentration
and temperature dependence of the p-cresol absorption
band in the presence of an equimolar amount of dodecylguanidine (Figure ). Exact stoichiometric
equivalence is easily obtained by mixing together ∼1 mM hexane
solutions of the two components, and then collecting the resulting
crystalline precipitate (“dodecylguanidinium p-cresolate”).
This precin class="Chemical">pitated salt has a 1:1 fixed
stoichiometry regardless
of relative initial concentrations, as shown by relative 1H NMR peak areas after the salt is collected and redissolved in DMSO.
The solid dodecylguanidinium p-cresolate salt has
very low solubility in hexane at room temperature, but it can easily
be redissolved to a concentration of ∼10–3 M in hexane at temperatures >55 °C. UV absorbance spectra
were
measured as function of further dilution with hexane (Figure A), or further temperature
elevation (Figure B). These spectra clearly show a broad absorption band centered near
310 nm. This species likely has a structure related to phenolate anion,
because pure sodium p-cresolate also gives a broad
absorption band near 310 nm in DMSO[15] (see
also Supporting Information).
It
is a priori unlikely that such a low-dielectric solvent as n class="Chemical">pan class="Chemical">hexane
(ε = 1.8) would contain indepn>endently solvated pan class="Chemical">dodecylguanidinium
and pan class="Chemical">p-cresolate ions. This hypothesis is supported
experimentally by the observed dilution-dependent increases in the A280/A310 spectral
absorbance ratio in Figure A. This is not what is expected, if a simple Brønsted–Lowry
acid–base equilibrium prevails, that is, ϕ-OH + :N-guan
⇄ ϕ-O– + HN+-guan. For such
an equilibrium, the ratio of [ϕ-OH]/[ϕ-O–] should be unaffected by dilution with additional solvent. (As a
control, such Brønsted–Lowry behavior was shown to occur
when methanol was used as a solvent for the same equilibrium; see
below.)
Instead, we determined that concentration- and temperature-den class="Chemical">pendent
spectral changes in hexane could be better modeled by the equilibrium,
ϕ-OH + :N-guan ⇄ ϕ-OH:N-guan, which describes the
formation of a Lewis acid–base complex. Note that what is written
as a unitary complexation product on the right-hand side probably
involves multiple interconverting protomers, such as the ones shown
in Scheme B,C.
The simplest demonstration of this is a log–log n class="Chemical">plot of pan class="Chemical">A308 versus dissolved concentration of pan class="Chemical">dodecylguanidinium p-cresolate (Figure A inset). This shows a limiting slope of ∼2 at low
concentrations, a clear sign that the 310 nm absorption band results
from the formation of a binary complex. A higher-order complex would
be expected to show a higher-power concentration dependence.
It was possible to determine more n class="Chemical">precise values of equilibrium
concentrations by fitting a wider range of wavelengths simultaneously.
At a total concentration near 1 mM and tempn>eratures near 60 °C,
the 1:1 mixture of pn>an class="Chemical">p-cresol and dodecylguanidine
shows a broad spectral shape with a variable plateau from 298 to 310
nm (Figure A,B). The
relative size of this plateau decreases either upon overall dilution
with hexane (Figure A) or upon heating (Figure B). The reversibility of the complexation reaction is demonstrated
by the nearly perfect superposition of spectra in Figure B obtained during heating (from
58 to 64 °C) and then cooling (from 64 to 58 °C), at intervals
of 2 °C There was evidence for ∼1% overall decrease in p-cresol content over the course of the entire heating–cooling
cycle, as evidenced by a uniform decrease in the entire absorption
spectrum in the final spectrum at 58 °C, as compared with the
initial spectrum at the same temperature. This continual slow loss
of sample was corrected (to first order) during the subsequent data
analysis by averaging the identical-temperature spectra pairwise from
the heating–cooling cycle.
By subtracting out differently
weighted amounts of the pure n class="Chemical">pan class="Chemical">p-cresol spn>ectrum in
pan class="Chemical">hexane (Figure C,
black trace), it required only a simple
least-squares fitting procedure to compute the shape of the UV spectrum
of the single additional component (Figure C, red trace) that could best account for
all the measured spectra in Figure A,B. The fitting procedure utilized singular value
decomposition (pan class="Disease">SVD) to minimize the sum-of-squares in the residual
absorbance spectra (thin black traces in Figure A,B) that could not be accounted for by linear
combinations of just p-cresol and one other species.
This pan class="Disease">SVDn>-based fitting procedure gave a unique spectral shape,
but not its vertical scale. That is, based solely on the SVD procedure,
the fitted spectral bandshape (Figure C, red trace) could be rescaled vertically by any arbitrary
factor. Choosing the scale for the red trace in Figure C required one additional assumption, namely,
that the eight different concentrations of the 1:1 complex, deduced
from the properly scaled extinction-coefficient spectrum, should give
the same equilibrium constant for 1:1 complex formation, Kcomplex, for all eight dilutions measured at 58 °C
(seven measurements in Figure A, plus one in Figure B). In fact, the scale was optimized to minimize the standard
deviation among eight resultant values of ΔG°, rather than among the values of Kcomplex itself (see Table 2-S in Supporting Information
for details).
The optimization gave a value of ΔG°
= −23.5 ± 0.8 kJ mol–1 at 58 °C,
corresn class="Chemical">ponding to Kcomplex = 5100 ±
1500 M –1. The reciprocal of this Kcomplex represents the dissociation constant in pan class="Chemical">hexane
at 58 °C, Kd = 0.195 ± 0.06
mM. The moderate standard deviation is generally consistent with the
initial hypothesis that 1:1 compn>lexation is the main reaction occurring
over the 12-fold concentration range used in these experiments, although
clearly this is not the whole story.
The best-fit extinction-coefficient
spectrum (Figure C,
red trace) can be most easily
intern class="Chemical">preted as involving two species in rapid internal proton-transfer
equilibrium, corresponding to structures B and C in Scheme . These structures give rise,
respectively, to “guanidinium/phenolate” absorbance
near 310 nm and “guanidine/phenol” absorbance near 280
nm in Figure C. However,
even the latter subcomponent shows a significantly different vibronic
substructure than the isolated p-cresol in hexane
(Figure C, black trace).
The scaled spectrum (Figure C, red trace) was also used to determine the temn class="Chemical">perature dependence
of the concentrations in Figure B, and thus of the equilibrium constant. A linear least-squares
fit of the resulting van’t Hoff plot (Figure B, inset) gave a value of ΔH° = −74 ± 4 kJ mol–1 for the formation of the 1:1 complex. (The stated uncertainty here
is simply R times the standard error of the slope,
obtained from the linear regression of ln K vs 1/T. Concentration errors from the SVD fit contribute to errors
in the four individual values for Kcomplex, but these end up being insignificant compared with 4 kJ mol–1, as can be judged by smaller size of the error in
ΔG°.) Combining this value of ΔH° with the value of ΔG°
at 58 °C, we obtain the value of ΔS°
= −152 ± 12 J mol–1 K–1 at 58 °C. We can further extrapolate these results to physiological
temperature (37 °C), obtaining a value of ΔG° = −26.7 ± 0.8 kJ mol–1 and a
corresponding Kd = 0.03 ± 0.01 mM
for the formation of a heterodimer complex in hexane.
UV Measurements
in More Polar Solvents
pan class="Chemical">DMSOn> offers
distinct adpan class="Disease">vantages for studying the interactions between pan class="Chemical">phenol and
guanidine groups. The most important is significantly greater solubility,
for example, for the covalently linked free bases (p-phenolyl)-dodecylguanidine and (p-phenolyl)ethylguanidine
and their HBr salts, as well as for potassium p-cresolate.
None of these compounds dissolved sufficiently in alkanes to permit
measurement of an interpretable UV spectrum. Furthermore, dodecylguanidine
itself could be dissolved to a much higher concentration in DMSO (∼50
mM), without producing strongly scattering micelles as was consistently
observed above ∼2 mM in hexane. Similar high concentrations
could also be obtained in the protic solvent MeOH, which served as
an important control.
UV spectral measurements inn class="Chemical">pan class="Chemical">DMSO (Figure A–D) also
pn>rovide evidence for pan class="Chemical">Lewis acid–base complexation rather than
pan class="Chemical">Brønsted–Lowry H+ transfer. The 1:1 dodecylguanidinium p-cresolate salt, obtained as above, when redissolved at
0.1 and 1 mM in DMSO (Figure B, black and red traces) gave extinction-coefficient spectra
superimposable on those of pure p-cresol (Figure A, black and red
traces). Thus, in DMSO, the presence of <1 mM concentrations of
the very strong base, dodecylguanidine, was insufficient to deprotonate
a similar concentration of p-cresol, a weak acid,
contrary to expectation if Brønsted–Lowry acid–base
equilibrium predominated.
Figure 4
UV spectral detection of phenol–guanidine
interactions,
measured using DMSO (A–D) or MeOH (E–H) as solvent.
Each column represents spectra obtained on the solute molecule shown
above (or two solutes, in the case of the second column). Spectra
are presented as extinction-coefficient data, to correct for varying
path lengths and concentrations used during measurements, and are
color-coded according to solute concentration(s). Every panel (A–H)
includes one spectrum of the indicated phenol-containing compound
at 0.1 mM (black), and one at 1.0 mM (red); in (B) and (F), those
red and black traces also include 0.1 or 1.0 mM dodecylguanidine free
base, respectively. (B) and (F) each also include three additional
traces in green, blue, and purple, taken respectively with 5, 10,
and 50 mM total concentration of dodecylguanidine free base; all of
these still have only 1.0 mM p-cresol. All spectra
were measured at 25 °C with 1 cm pathlength sample and reference
cuvettes on a Shimadzu UV-265 spectrometer. Absolute vertical scales
of the spectra are accurate only to ±20% due to weighing and
dilution imprecision and were adjusted within this range to facilitate
comparisons.
UV spectral detection of n class="Chemical">pan class="Chemical">phenol–pn>an class="Chemical">guanidine
interactions,
measured using DMSO (A–D) or MeOH (E–H) as solvent.
Each column represents spectra obtained on the solute molecule shown
above (or two solutes, in the case of the second column). Spectra
are presented as extinction-coefficient data, to correct for varying
path lengths and concentrations used during measurements, and are
color-coded according to solute concentration(s). Every panel (A–H)
includes one spectrum of the indicated phenol-containing compound
at 0.1 mM (black), and one at 1.0 mM (red); in (B) and (F), those
red and black traces also include 0.1 or 1.0 mM dodecylguanidinefree
base, respectively. (B) and (F) each also include three additional
traces in green, blue, and purple, taken respectively with 5, 10,
and 50 mM total concentration of dodecylguanidinefree base; all of
these still have only 1.0 mM p-cresol. All spectra
were measured at 25 °C with 1 cm pathlength sample and reference
cuvettes on a Shimadzu UV-265 spectrometer. Absolute vertical scales
of the spectra are accurate only to ±20% due to weighing and
dilution imprecision and were adjusted within this range to facilitate
comparisons.
We did not obtain measurements
of the 1:1 pan class="Chemical">saltn> in pan class="Chemical">DMSO at concentrations
above 1 mM, due to too high pan class="Chemical">p-cresol absorption for
the 1 cm pathlength cuvette we were using. However, when additional
crystalline dodecylguanidinefree base was added to the 1 mM salt,
a new absorption maximum could be seen at ∼318 nm (Figure B). This increased
in intensity as the total dodecylguanidine concentration was raised
successively to 5, 10, and 50 mM (Figure B, green, blue, and violet traces, respectively).
This 318 nm peak is consistent with Lewis acid–base complex
formation, analogous to that seen in hexane (Figure A). However, the appearance of a new band
at 272 nm at the highest (50 mM) dodecylguanidine concentration (Figure B) makes clear that,
compared with hexane, it is less reasonable in DMSO to analyze the
equilibrium based on the assumption of only two species (isolated p-cresol and its 1:1 complex with dodecylguanidine). From
the spectra in Figure B, we can only estimate an apparent Kd for the p-cresol/dodecylguanidine complex of 10–50
mM in DMSO at 25 °C. This is the dodecylguanidine concentration
for which 50% of the maximum possible 318 nm cresolate absorbance
is achieved. This apparent Kd is at least
∼100-fold greater than in hexane (extrapolated to a similar
temperature), indicating a weaker H-bonding interaction in DMSO.
Complexation inn class="Chemical">pan class="Chemical">DMSO clearly occurs at a lower concentration for
covalently linked pn>an class="Chemical">(p-phenolyl)alkylguanidines (Figure C,D) than is required
for p-cresol/dodecylguanidine (Figure B). That is, the spectra of (p-phenolyl)ethylguanidine and (p-phenolyl)dodecylguanidine
in DMSO (Figure C,D)
show clear evidence of Lewis acid–base complexation, as evidenced
by the appearance of a 318 nm band even at concentrations of ∼1
mM but not at 0.1 mM. This complexation is probably intermolecular
rather than intramolecular, given the absence of the 318 nm peak at
the lower (0.1 mM) concentration, and its similar size at 1 mM, for
both alkyl chain lengths. The short ethyl linker clearly cannot permit
the formation of an intramolecular H-bond. The dodecyl linker also
does not appear to favor intramolecular H-bonding in DMSO solvent,
despite computations (Figure B), indicating that a conformational rearrangement to permit
intramolecular H-bond formation is energetically favored in vacuo.
Empirically, the extended chain form appears to allow much more favorable
external interactions in condensed phases, that is, both in solid
state (as shown in crystal structures below) and for either dimers
or monomers solvated in DMSO, as suggested by Figure C,D. There is a barely perceptible 318 nm
absorbance above the baseline for 0.1 M (p-phenolyl)dodecylguanidine
in DMSO (Figure D,
black trace), which is not present for the ethyl-linked compound (Figure C, black trace).
This could be a sign of a small amount of intramolecular H-bonding
in the former. However, we cannot rule out that this could simply
be due to a minor difference in light scattering.
pan class="Chemical">Freen>-energy
and enthalpy changes for (p-phenolyl)ethylguanidine
dimerization in DMSO were determined by using concentration- and temperature-dependent
UV spectral anapan class="Chemical">lyses, as described above for the p-cresol/dodecylguanidine system in hexane. Due to considerably higher
solubility, the spectra in DMSO could be measured over a wider temperature
range (19–45 °C), as well as a wider concentration range
(0.1–10 mM; Figure C and Supporting Information).
The thermodynamic values obtained for (p-phenolyl)ethylguanidine
dimerization in DMSO are ΔG°298 = −8.4 ± 0.6 kJ mol–1 and ΔH° = −22 ± 3 kJ mol–1. Because there are two H-bonding interactions per dimer, the enthalpy
change for forming each phenol/guanidine H-bond is only −11
± 2 kJ mol–1, that is, barely one-sixth of
the value in hexane.
In contrast to the complexation behavior
inn class="Chemical">pan class="Chemical">hexane and pn>an class="Chemical">DMSO, the
same molecules in MeOH simply undergo Brønsted–Lowry acid–base
equilibria. This is seen in Figure F, where sufficiently high relative concentrations
of the base (dodecylguanidine) fully deprotonate the acid (p-cresol), as evidenced by the complete replacement of the
280 nm absorption band by the 296 nm band characteristic of cresolate
in MeOH. That is, the purple trace in Figure F, corresponding to 1 mM p-cresol + 50 mM dodecylguanidine, matches almost perfectly the spectrum
of 1 mM potassium p-cresolate in MeOH (data not shown).
There is a fairly clean isosbestic point in Figure F, indicating that only two species (solvated p-cresol and solvated p-cresolate) contribute
to the spectra in this panel.
Likewise, the lack of concentration
dependence in the sn class="Chemical">pectra in Figure G,H is indicative
of the presence of buffering concentrations of conjugate acid and
base forms of both groups within a single amphiprotic solute. The
concentration-independent acid–base equilibrium of these (p-phenolyl)alkylguandine compounds is analogous to other
amphiprotic zwitterions, such as amino acids. When present as the
only solute in water or other protic solvents, such zwitterions buffer
the pH at the midpoint between the two pKa values and show almost no concentration dependence for the extent
of H+ transfer.
The presence of both n class="Chemical">pan class="Chemical">phenol and pn>an class="Chemical">phenolate
forms in Figure F–H
indicates that in
methanol, the pKa value of the phenol
group in this solvent must be quite similar to the pKa of the guanidinium. In the limit of phenol and guanidinium
groups sitting on separate molecules in MeOH, ΔpKa ≈1.6. This is seen in Figure F, where 50% conversion of p-cresol to p-cresolate (corresponding roughly to
the green trace), requires approximately a 20-fold excess of total
guanidine (10 mM) over total p-cresol (0.5 mM). Thus,
surprisingly, in MeOH, which is often thought of as only a little
less polar than water, the relative acidity of these two groups is
already swapped, with guanidinium clearly more acidic than phenol.
However, as the two groups are held n class="Chemical">progressively closer, for example,
by a dodecyl or an ethyl linker, the phenol becomes progressively
more acidic, eventually becoming about as acidic as guanidine, as
seen by the nearly 50% contribution of phenolate to the spectrum (compare
296 nm peaks for black traces in Figure F–H). That is, the zwitterionic form
is favored by progressively closer distances between the charges,
a clear sign of electrostatic perturbation of the pKa values. This is analogous to the well-known phenomenon
in amino acids, where the pKa values of
amino and carboxylic acid groups are pushed farther apart, favoring
the zwitterionic form, as a result of these groups’ proximity
and mutual electrostatic stabilization.
IR Spectroscopy
Direct IR evidence for 1:1 complexation
was not as strong for IR as for UV sn class="Chemical">pectra, because measurement of
the concentration dependence of the IR solution spectra was not possible.
In nonpolar solvents (hexane, CCl4), the maximum room-temperature
solubility of (p-phenolyl)alkylguanidine compounds
was too low (<0.5 mM) to measure spectra accurately. In DMSO, the
background solvent absorbances were too high to permit precise quantitative
measurements even at saturating concentrations (∼150 mM). Nevertheless,
IR spectra are consistent with the presence of H-bonded complexes
of phenol with guanidinefree bases. Such spectra show clearly that
the H-bond between these groups is polarizable, that is, the equilibrium
position of the proton shifts moving increasingly toward the guanidine
group as the polarity of the environment increases from CCl4 to DMSO to MeOH to tightly packed crystals.
Nonpolar Solvents
Due to its low solubility, we failed
to obtain adequate-quality IR spectra of n class="Chemical">pan class="Chemical">p-cresol-dodecylguanidine
1:1 in pn>an class="Chemical">hexane, which would be directly analogous to the UV spectra
in Figure . The closest
solvent we could use is CCl4 (Figure ). This spectrum supports formation of a
binary H-bonded complex, as indicated by characteristic bands due
to both p-cresol (1514 and 1255 cm–1) and p-cresolate (1498 and 1272 cm–1).[16,18] The characteristic C–N stretch region
of the guanidino group, between 1500 and 1700 cm–1, also shows features characteristic of both deprotonated guanidine
(1629 and 1555 cm–1) and guanindinium ion (1665
cm–1). A key feature observed is the strong broad
absorption band between 3500 and 2500 cm–1, a “continuum”
absorbance characteristic of polarizable H-bonded systems.
Figure 5
IR absorbance
spectrum of p-cresol/dodecylguanidine
1:1 dissolved in CCl4. The strong continuum absorbance
from 3500 to 2500 cm–1 is a clear evidence of the
formation of a strongly H-bonded complex, rather than a simple Brønsted–Lowry
proton transfer equilibrium. The exact position and intensity of the
most characteristic peak of the deprotonated guanidine group (1555
cm–1) are difficult to quantify due to interference
by a strong solvent absorption band at this frequency. Spectral region
near 2350 cm–1 is blanked due to interference by
CO2 absorption.
IR absorbance
spectrum of n class="Chemical">pan class="Chemical">p-cresol/pn>an class="Chemical">dodecylguanidine
1:1 dissolved in CCl4. The strong continuum absorbance
from 3500 to 2500 cm–1 is a clear evidence of the
formation of a strongly H-bonded complex, rather than a simple Brønsted–Lowry
proton transfer equilibrium. The exact position and intensity of the
most characteristic peak of the deprotonated guanidine group (1555
cm–1) are difficult to quantify due to interference
by a strong solvent absorption band at this frequency. Spectral region
near 2350 cm–1 is blanked due to interference by
CO2 absorption.
DMSO and MeOH Solutions
It was possible to obtain sn class="Chemical">pectra
of (p-phenolyl)dodecylguanidine in DMSO (Figure A) and in MeOH (Figure C), but only at ∼100-fold
higher concentrations than for the corresponding UV spectra (Figure D,H). The strongest
characteristic IR bands of the protonated forms of both guanidinium
(∼1670 cm–1 shoulder) and phenol (1514 cm–1) are present in DMSO (Figure A). The additional bands at 1642 and (especially)
1548 cm–1 are characteristic of deprotonated guanidine
in either DMSO[12,17] or chloroform.[19] Meanwhile, characterisitic phenolate bands are observed
at 1599 and 1494 cm–1 (Figure A), that is, shifted somewhat relative to
this compound in MeOH (Figure B, 1602 and 1501 cm–1). The latter wavenumber
values are identical with those seen for high-pH aqueous phenolate[18] and tyrosine.[16]
Figure 6
IR spectra
of (p-phenolyl)dodecylguanidine, as
the “free base” (green; A, B) or as the HBr salt (purple;
C, D), with the former two showing clear evidence of H-bonding and/or
proton transfer from the phenol. Each sample was recrystallized from
methanol, dried, and then redissolved to a concentration of ∼100
mM in DMSO (A, C) or methanol (B, D). The absorbance scale on the y axis is approximate; the measured spectra had somewhat
different concentrations and path lengths and were rescaled for optimum
visual comparison.
IR spectra
of n class="Chemical">pan class="Chemical">(p-phenolyl)dodecylguanidine, as
the “pn>an class="Chemical">free base” (green; A, B) or as the HBr salt (purple;
C, D), with the former two showing clear evidence of H-bonding and/or
proton transfer from the phenol. Each sample was recrystallized from
methanol, dried, and then redissolved to a concentration of ∼100
mM in DMSO (A, C) or methanol (B, D). The absorbance scale on the y axis is approximate; the measured spectra had somewhat
different concentrations and path lengths and were rescaled for optimum
visual comparison.
Thus, these IR spectra
confirm that inn class="Chemical">pan class="Chemical">DMSO, this pn>an class="Chemical">free base compound
exists in an equilibrium between a guanidinium/phenolate zwitterion
(Scheme B) and a fully
neutral guanidine/phenol form (Scheme C). The relative sizes of the bands suggest that the
zwitterion accounts for some 50–75% of the total concentration,
more than seen in the UV absorption spectra in Figure D. This can be explained by the presence
of a higher proportion of the H-bonded complex, as a result of a ∼100-fold
higher overall concentration in the IR samples.
An even greater
contribution from the zwitterionic form is observed
in pan class="Chemical">MeOHn> (Figure B)
(e.g., lpan class="Chemical">arger characteristic pan class="Chemical">guanidinium peak near 1670 cm–1 in Figure B compared
with an unlabeled shoulder in Figure A). This vibration is characteristic of guanidinium;
it is clearly seen in the spectra of the corresponding HBr salt in
either DMSO or MeOH (Figure C,D). Likewise, in MeOH, phenolate peaks are stronger (1602
and 1501 cm–1 in Figure B). Similar features were observed in the
spectrum of (p-phenolyl)ethylguanidinefree base
dissolved in methanol (data not shown). These spectra also indicate
that the zwitterion contribution to the spectrum is increased significantly
by the use of >100-fold higher concentration for the IR spectrum
(Figure B), as compared
with
the <1 mM concentrations used to measure UV spectra (Figure H). This might contradict somewhat
the hypothesis of simple Brønsted–Lowry acid–base
equilibrium. Alternatively, it could simply result from favorable
long-range electrostatic (dipole–dipole) interactions coming
into play at the highest concentrations.
Very different IR results
are obtained with the pan class="Chemical">saltn> (p-phenolyl) dodecylguanidiniumpan class="Chemical">bromide, which corresponds to the structure
in Scheme A, with
X– = Br–. This HBr salt has a
much higher solubility in most solvents compared with the free-base
compound, especially in DMSO. This makes it easier to measure accurately
over a wider spectral range (Figure ). In both DMSO (Figure C; and Figure , purple trace) and MeOH (Figure D), this salt contains a solvated fully protonated
guanidinium, with no evidence of interaction with the phenolic group.
The monoalkylguanidinium C–N stretch bands at 1668 and 1651
cm–1 in DMSO (Figure C; Figure , purple trace) are largely unchanged from their appearance
in dodecylguanidinium bromide in the same solvent (Figure , red trace).[5,17]
Figure 7
IR
evidence that the phenol and guanidinium groups in (p-phenolyl)dodecylguanidinium bromide are independently
solvated by DMSO, that is, there is no significant interaction between
them. The IR spectrum of this salt (purple) resembles a simple sum
of absorbances due to separate samples of dodecylguanidinium bromide
(red) and p-cresol (blue). All spectra were obtained
in DMSO.
IR
evidence that the pan class="Chemical">phenoln> and pan class="Chemical">guanidinium groups in pan class="Chemical">(p-phenolyl)dodecylguanidinium bromide are independently
solvated by DMSO, that is, there is no significant interaction between
them. The IR spectrum of this salt (purple) resembles a simple sum
of absorbances due to separate samples of dodecylguanidinium bromide
(red) and p-cresol (blue). All spectra were obtained
in DMSO.
Likewise, characteristic pan class="Chemical">phenoln>
peaks are observed near 1615, 1594,
1515, 1266, and 1236 cm–1 in (p-phenolyl)dodecylguanidinium bromide (Figure C,D; Figure , purple trace), virtually unchanged from p-cresol (Figure ,
blue trace). In fact, these characteristic pan class="Chemical">phenol peaks are generally
very insensitive to the environment. This is shown in greater detail
in Figure . The C–O
stretch vibration is consistently observed as a strong IR absorbance
band for p-alkylated phenols in a fairly narrow range
between 1245 and 1275 cm–1. The lowest frequencies
in this range have been reported for poly-tyrosine in dried films
or aqueous solvents.[18] In Figures and 8A–C, which are limited to crystalline samples and solution
measurements in methanol and DMSO, the variation is within a narrower
range of 1265–1268 cm–1.
Figure 8
Phenolic C–O stretch
region (1250–1310 cm–1) is less sensitive
to solvent environment when the phenolic group
is protonated (A–C) than when it is unprotonated (D–G).
(A) p-Cresol in DMSO; (B) (p-phenolyl)dodecylguanidinium
bromide in DMSO, continuation of spectrum in Figure C; (C) (p-phenolyl)ethylguanidinium
bromide in MeOH, continuation of Figure D; (D) (p-phenolyl)ethylguanidine
free base in methanol, continuation of Figure B; (E) (p-phenolyl)dodecylguanidine
free base in DMSO, continuation of Figure A; (F) p-cresol and dodecylguanidine
free base (1:1) in CCl4, expanded region from Figure ; (G) potassium p-cresolate in DMSO.
pan class="Chemical">Phenoln>ic C–O stretch
region (1250–1310 cm–1) is less sensitive
to solvent environment when the pan class="Chemical">phenolic group
is protonated (A–C) than when it is unprotonated (D–G).
(A) pan class="Chemical">p-Cresol in DMSO; (B) (p-phenolyl)dodecylguanidiniumbromide in DMSO, continuation of spectrum in Figure C; (C) (p-phenolyl)ethylguanidiniumbromide in MeOH, continuation of Figure D; (D) (p-phenolyl)ethylguanidinefree base in methanol, continuation of Figure B; (E) (p-phenolyl)dodecylguanidinefree base in DMSO, continuation of Figure A; (F) p-cresol and dodecylguanidinefree base (1:1) in CCl4, expanded region from Figure ; (G) potassium p-cresolate in DMSO.
In contrast, deprotonation of the n class="Chemical">pan class="Chemical">phenolic groupn> leads to
a much
lpan class="Chemical">arger range of C–O stretch frequencies, as seen in Figure D–G. The lowest
observed value, 1264 cm–1, is for pan class="Chemical">(p-phenolyl)ethylguanidine free base in MeOH (Figure D). This value overlaps with the range seen
for protonated phenol (Figure A–C), yet it was measured under conditions where the
phenolate form predominates, as shown by the C=C ring modes
(1602, 1501 cm–1 in Figure B). At the same time, the structurally similar
compound, (p-phenolyl)dodecylguanidine in DMSO, appears
to show a 40 cm–1 higher frequency for its phenolate
C–O stretch (1309 cm–1, Figure E).
These big differences
in C–O stretch frequency for pan class="Chemical">phenolaten>
in various environments are likely attributable to a strong dependence
of the C–O bond order on solvation and H-bond strength. These
also produce changes in the degree of delocalization of negative charge
into the ring, weakly altering other characteristic vibrational frequencies.
For example, the phenolate C=C ring stretching vipan class="Chemical">brations are
lower in frequency for the H-bonded complex of (p-phenolyl)dodecylguanidine in DMSO (1599 and 1494 cm–1, Figure A) than
in MeOH (1602 and 1501 cm–1, Figure B). The corresponding frequencies for tyrosinate
in water at a high pH[16,18] are intermediate, that is, 1600
and 1498 cm–1, suggesting that the H-bonding interactions
of phenolate with water are stronger than with MeOH, but that neither
type of solvent produces an H-bonding perturbation on phenol that
is as strong overall as a single H-bond to guanidine in an aprotic
solvent such as DMSO. The protonated phenol group again serves as
a useful control; its solvent interactions are weaker than those of
phenolate, that is, only small solvent-dependent effects are observed
for the phenol ring vibrations (near 1615, 1595, and 1515 cm–1 in Figures A,C,D
and 7).
Solid-State IR Spectra
An important question is, what
do Fourier transform infrared (FTIR) sn class="Chemical">pectra indicate for the protomeric
form(s) present in the crystalline state, specifically in crystals
that are suitable for determining high-resolution X-ray structures?
In such crystals, the polar ends of neighboring dimers are packed
close to each other, and also within the H-bonding distance of structural
MeOH molecules (see below). Thus, the environment around the guanidine–phenol
H-bond is quite polar.
Indeed, in the IR spectrum of these crystals
(green trace in Figure ), only the n class="Chemical">pan class="Chemical">phenolate pn>eak at 1498 cm–1 is observed;
the peak at 1516 cm–1, characteristic of protonated
(neutral) phenol, is completely missing. Additionally, in the crystalline
forms of both the free base (Figure , green trace) and the corresponding HBr salt (Figure , purple trace),
the strongest C–N stretch bands are very similar in frequency,
near 1655 cm–1. This is a lower wavenumber than
is typically seen for pan class="Chemical">monoalkylguanidinium salts in most solvents
(cf. Figures and 7). Nonetheless, it agrees with the measurements
in the crystalline state of the highest-frequency C–N stretch
vibration for dodecylguanidinium bromide, anisolyl dodecylguanidiniumbromide, phenolyl ethylguanidine free base, and phenolyl ethylguanidiniumbromide, which are consistently observed between 1659 and 1669 cm–1 (Figure 6-S). In all these
crystals, therefore, the free base exists mostly as the zwitterionic
form, with a guanidinium group (Scheme B). In contrast, the highest C–N stretch frequencies
for the crystalline states of authentic monoalkylguanidinefree bases,
such as dodecylguanidine or anisolyldodecylguanidine, are 1637 and
1636 cm–1, respectively (see spectra C and E in Figure 6-S).
Figure 9
Broad
IR continuum absorbance 3500–2500 cm–1 in
FTIR spectrum of a crystalline sample of p-phenol-dodecylguanidine
free base (green) is the evidence that it is more strongly H-bonded
than the corresponding crystalline p-phenol-dodecylguanidinium
bromide salt (purple). Both samples were measured in KBr disks. Blanked
regions had errors due to background correction, for example, CO2 absorbance near 2350 cm–1.
Nevertheless, the IR spectrum
of the n class="Chemical">pan class="Chemical">p-phenolyl-dodecylguanidine
pn>an class="Chemical">free base (Figure , green trace) may be consistent with the presence of small amounts
of the strongly H-bonded but fully neutral structure (Scheme C). The new peaks, which are
observed only in the free base crystals (Figure , green trace), and not the HBr salt (Figure , purple trace),
include a characteristic deprotonated guanidine frequency (1554 cm–1), as well as a broad continuum absorption band from
2500 to 3500 cm–1. Both of these features can also
be discerned in the crystalline samples of either (p-phenolyl)ethylguanidine, or (less prominently) dodecylguanidinium p-cresolate salt (see Supporting Information). Because of the absence of characteristic vibrational bands of
MeOH, at most, a small portion of the continuum bands (above ∼3200
cm–1) can be attributed to O–H stretch vibrations
from the cocrystallized MeOH in the phenolyl alkylguanidinefree bases
(see crystal structures below). Instead, these broad absorption bands
are likely indicative of proton polarizability,[20] which is associated with a strong H-bonding.[21,22]
pan class="Chemical">Brn>oad
IR continuum absorbance 3500–2500 cm–1 in
FTIR spectrum of a crystalline sampn>le of p-phenol-dodecylguanidinefree base (green) is the evidence that it is more strongly H-bonded
than the corresponding crystalline p-phenol-pan class="Chemical">dodecylguanidinium
bromide salt (purple). Both samples were measured in KBr disks. Blanked
regions had errors due to background correction, for example, CO2 absorbance near 2350 cm–1.
1H and 13C NMR Measurements
of the Phenol–Guanidine
Interaction
Clear supn class="Chemical">port for the formation of an H-bonding
interaction between phenol and guanidine groups was also observed
in the solution-state 1H and 13C NMR spectra. Figure shows that at
a low concentration of (p-phenolyl)ethylguanidine,
the chemical shifts of both the guanidine group and the phenol group
are much the same as when measured individually in dodecylguanidine
or p-cresol, respectively. However, as the concentration
is increased above ∼30 mM, there is a significant change, especially
in the value for the guanidine protons.
Figure 10
Concentration dependence
of the 1H chemical shifts of
the phenolic ring (gray solid circle, orange solid circle) and of
the guanidine group (blue solid circle), in (p-phenolyl)ethylguanidine.
Average values are shown for peaks split by j–j coupling. Likewise, for the monoalkylguanidinium bromide salts,
there were two to three partially resolved resonances for the guanidino
group protons in the range of 6.9–7.5; these were averaged
with weightings of 1:2:2 to give the plotted value of 7.25 ppm (blue
line next to the right axis). Well below ∼30 mM, the chemical
shifts in (p-phenolyl)ethylguanidine are nearly identical
with those seen for isolated phenol and guanidine groups, indicated
by lines next to the left vertical axis. The changing values over
the 20–150 mM concentration range are consistent with the formation
of an H-bonded dimer complex. That is, at the highest concentrations,
the averaged chemical shift of the guanidine protons approaches that
of dodecylguanidinium cation, whereas the phenolic ring protons move
somewhat toward the values for phenolate. The latter are represented
by the chemical shifts of Br– or K+ salt
solutions in DMSO, respectively, as indicated next to the right axis.
Concentration dependence
of the n class="Chemical">pan class="Chemical">1H chemical shifts of
the pn>an class="Chemical">phenolic ring (gray solid circle, orange solid circle) and of
the guanidine group (blue solid circle), in (p-phenolyl)ethylguanidine.
Average values are shown for peaks split by j–j coupling. Likewise, for the monoalkylguanidinium bromide salts,
there were two to three partially resolved resonances for the guanidino
group protons in the range of 6.9–7.5; these were averaged
with weightings of 1:2:2 to give the plotted value of 7.25 ppm (blue
line next to the right axis). Well below ∼30 mM, the chemical
shifts in (p-phenolyl)ethylguanidine are nearly identical
with those seen for isolated phenol and guanidine groups, indicated
by lines next to the left vertical axis. The changing values over
the 20–150 mM concentration range are consistent with the formation
of an H-bonded dimer complex. That is, at the highest concentrations,
the averaged chemical shift of the guanidine protons approaches that
of dodecylguanidinium cation, whereas the phenolic ring protons move
somewhat toward the values for phenolate. The latter are represented
by the chemical shifts of Br– or K+ salt
solutions in DMSO, respectively, as indicated next to the right axis.
Figure shows
additional 13C NMR results in DMSO solution, which also
demonstrate an H-bonding interaction in (p-phenolyl)ethylguanidine
dimers at elevated concentrations. The chemical shift of the phenol
ring carbons is clearly dependent on concentration. The dependence
is strongest for the ipsocarbon, which shifts from 157.0 ppm at 3
mM (Figure A) to
158.7 ppm at 200 mM (Figure C), and for the para carbon (128.1–126.3 ppm; same
spectra).
Figure 11
13C NMR spectra (112–172 ppm portion) of the
indicated samples, calibrated to the chemical shift value for the
solvent (d6-DMSO, 39.5 ppm). (A)–(G)
were measured with a 600 MHz Bruker NMR and (H)–(K), with a
300 MHz instrument. The chemical shift values labeled in blue are
assigned to the guanidine carbon; in black, to the six phenolic ring
carbons. The two largest phenolic peaks always correspond to ortho
and meta positions, each representing 2 equiv carbons. Horizontally
expanded and horizontally shifted insets for (D), (F), and (G) aid
in visualizing spectral details near 157 ppm. Spectrum (F) was measured
at 2× finer resolution than the others, the better to resolve
the nearly overlapping phenolic ipso and guanidine resonances at 157.2
ppm.
pan class="Chemical">13Cn> NMR spectra (112–172 ppm portion) of the
indicated sampn>les, calibrated to the chemical shift value for the
solvent (d6-DMSO, 39.5 ppm). (A)–(G)
were measured with a 600 MHz Bruker NMR and (H)–(K), with a
300 MHz instrument. The chemical shift values labeled in blue are
assigned to the guanidinecarbon; in black, to the six phenolic ring
pan class="Chemical">carbons. The two largest phenolic peaks always correspond to ortho
and meta positions, each representing 2 equiv carbons. Horizontally
expanded and horizontally shifted insets for (D), (F), and (G) aid
in visualizing spectral details near 157 ppm. Spectrum (F) was measured
at 2× finer resolution than the others, the better to resolve
the nearly overlapping phenolic ipso and guanidine resonances at 157.2
ppm.
Pron class="Chemical">perly assigning the phenolic
ipsocarbons within their overlapping
ranges with the guanidinecarbon required careful measurements of
a number of control spectra (Figure D–K). To minimize perturbations on the guanidine
group from H2O impurities, which could lead to formation
of guanidinium hydroxide rather than the free base, these were measured
on samples that had been recrystallized from acetonitrile, protected
from moisture during storage, weighed quickly, then redissolved at
100–200 mM by adding each directly to a freshly opened glass
ampule of d6-DMSO. Comparisons between
monoalkylguanidinefree bases lacking other H-bond donor groups (Figure F,J), and their
corresponding HBr salts (Figure G,K) show firstly that protonation of guanidine produces
only a ∼0.5 ppm change in its 13C chemical shift.
These control measurements made clear that regardless of the protonation
state, the 13C chemical shifts for monoalkylguanidine compounds
in DMSO can safely be assigned within the narrow range of 156.7–157.3
ppm.
The lowest concentration of pan class="Chemical">(p-phenolyl)ethylguanidinen>
(Figure A) showed
closely spaced 13C resonances at 157.0 and 157.2 ppm. These
could be definitively assigned based on a comparison with the similar
compound (p-phenolyl)dodecylguanidine (Figure D). The latter
was measured as a 50/50 mix of the natural-abundance compound with
an otherwise-identical compound that had been prepared with both terminal
pan class="Chemical">nitrogens labeled with >97% 15N. Such labeling produced
a triplet pattern due to j–j coupling between 13C and two attached 15N
nuclei, which could easily be seen for the peak at 157.2 (Figure D inset). The fourth
component peak at 157.2 ppm, which is the largest and nearest to the
center of the triplet pattern, is due to the ∼50% of sample
with 14N on all three of the guanidinenitrogens. The coupling
pattern in this sample allows the 157.2 peak in Figure D to be assigned to the guanidine
group, and by analogy also in Figure A. The minor differences in phenolic ipso-carbon chemical
shifts between Figure A,D might be due to the difference in the H-binding geometries, as
suggested by different X-ray crystal structures (see below).
As with the UV spectra in Figure D and the IR sn class="Chemical">pectra in Figure A above, the concentration-dependent NMR
results in both Figures and 11 are consistent with the formation
of an H-bonded Lewis acid–base dimer of (p-phenolyl)ethylguanidine, rather than the formation of solvated zwitterion
monomers with complete H+ transfer. In particular, even
at the highest concentrations, the 1H and 13C chemical shift values for the phenol group never come close to
those seen for potassium p-cresolate. The concentration
range over which the chemical shifts are seen to change is approximately
consistent with the dissociation constant (34 mM) computed from the
concentration dependence of the UV spectra of (p-phenolyl)ethylguanidine
(see ΔG° value given above; and raw data
in Supporting Information). Thus, the NMR
results confirm that in DMSO, phenol and guanidine groups readily
form strong H-bonded Lewis acid–base complexes, with partial
H+ transfer, that is, a hybrid of the structures shown
in Scheme B,C.
At the highest concentration for pan class="Chemical">(p-phenolyl)ethylguanidinen>
in Figure (0.05
M), the average chemical shift of the pan class="Chemical">guanidine protons (5.9) had
risen just over 50% of the way from the value for pan class="Chemical">dodecylguanidine
free base (4.3) to the value for protonated dodecylguanidinium bromide
(7.3). Overall, there is less sensitivity of the guanidine 13C NMR signal to its protonation state. Nevertheless, for (p-phenolyl)ethylguanidine at 0.05 M (Figure C), it is also shifted to 157.0, that is,
about halfway from the ∼157.3 ppm value for deprotonated guanidine
(Figure F,J) to
the ∼156.7 ppm value seen for protonated guanidinium salts
(Figure E,G,K).
This change is consistent with the predicted dimer fraction of 55%
(of the total mass) at 0.05 M (p-phenolyl)ethylguanidine,
based on Kd = 34 mM. Together, these chemical
shifts suggest that the guanidinium group in the dimer complex behaves
as if it is fully protonated, that is, the proton transfer from phenol
to guanidine is nearly complete in the dimer in DMSO.
However,
at the same pan class="Chemical">(p-phenolyl)ethylguanidinen>
concentration, the pan class="Chemical">phenolic ring protons moved only ∼20% of
the way from their values for pan class="Chemical">p-cresol to p-cresolate (Figure ). The phenolic ipsocarbon 13C shift in Figure C, at 158.7 ppm,
is also only 25% of the way from its value in p-cresol
(155.2 ppm, Figure G) to the value in p-cresolate (Figure H). Likewise, the 13C shifts for the other ring carbons in the (p-phenolyl)ethylguanidine
dimer (Figure C)
are very similar to those of p-cresol (Figure I), and quite different
from those of the potassium salt of p-cresolate (Figure H). Based purely
on an assumption of time-averaging of chemical shifts, this would
suggest that the phenol group is less than 50% deprotonated in the
dimer complex, a conclusion that would be somewhat at odds with the
results from optical spectroscopy (both UV and IR).
In summary,
the NMR chemical shifts in Figure A make the pan class="Chemical">guanidinen> group look about 50%
protonated, whereas the phenolic group looks nearly ∼75% protonated.
That is, by NMR both groups participating in the H-bond look predominantly
protonated, despite the lack of HBr or other possible external proton
donor. This is not due to contamination of the NMR solvent (pan class="Chemical">DMSO-d6) with water (or another acid), which could
indeed cause both groups of the free base to appear protonated because
of the formation of guanidinium hydroxide (or another salt). Contrary
to what would be expected in such hypothetical cases, simple dilution
with the same DMSO-d6 causes a breakup
of the H-bonded complex and results in the guanidine group becoming less protonated, that is, both its 1H and 13C signals move closer to the values seen in guanidinefree
bases without strong H-bond donor groups.
15N NMR
We were unable to observe pan class="Chemical">15Nn> chemical shifts of the terminal
pan class="Chemical">guanidine pan class="Chemical">nitrogens for the phenolyl-dodecylguanidinefree base in DMSO, despite the use of 50% 15N isotopic
enrichment. This was true over a range of temperatures from 25 to
50 °C, even after sufficient signal-averaging that we could observe
the 15N resonance from the natural-isotope-abundance intrachain
nitrogen at 63 ppm (see Figure 11-S). However,
we were able to obtain high-quality solid-state cross-polarized magic-angle-spinning
(CP-MAS) signals from crystalline samples of this compound, as well
as a number of control samples (Figure ).
Figure 12
15N solid-state CP-MAS NMR spectrum
of the free bases
of (p-anisolyl)-dodecylguanidine-[15N2] (A); (p-phenolyl)dodecylguanidine-[15N2] (B); (p-phenolyl)ethylguanidine-[15N2] (C); and the corresponding HBr salts (D-F,
respectively). The x axis scale represents chemical
shift, in ppm, relative to solid ammonium-[15N] sulfate.
15N solid-state CP-MAS NMR spectrum
of the free bases
of (p-anisolyl)-dodecylguanidine-[15N2] (A); (p-phenolyl)dodecylguanidine-[15N2] (B); (p-phenolyl)ethylguanidine-[15N2] (C); and the corresponding HBr salts (D-F,
respectively). The x axis scale represents chemical
shift, in ppm, relative to solid ammonium-[15N] sulfate.The first control spectrum is
noteworthy. Crystalline n class="Chemical">pan class="Chemical">(p-anisolyl)-dodecylguanidine-[15N2] pn>an class="Chemical">free base (Figure A) shows widely separated resonances for
imine (110 ppm) and amine
(52 and 46 ppm) nitrogens. A similar separation was also seen for
the dodecylguanidine-[15N2]free base.[12] However, the dodecylguanidine spectrum showed
only a single amine resonance at 49 ppm; addition of the anisolyl
group (Figure A)
splits this into two peaks. Only the two terminal nitrogens in the
(p-anisolyl)-dodecylguanidine were 15N-labeled,
so it is not immediately obvious why these crystals should show three
distinct resonances of similar intensity (Figure A). The likeliest explanation is the possible
presence of crystalline polymorphs. Another explanation is that in
the crystal unit cell, there could be two nonequivalent molecules
with different environments for their amino nitrogen, as was seen
for a number of crystalline argininesalts.[25] A full explanation of the amine resonance splitting in Figure A is not possible
in the absence of a crystal structure, which we have not yet been
able to obtain for (p-anisolyl)-dodecylguanidine.
Nevertheless, it is clear that the anisole group, which lacks an H-bond
donor, produces only a minor perturbation on the 15N spectrum
of the deprotonated dodecylguanidine.[12]
A very different situation prevails when a n class="Chemical">pan class="Chemical">phenol groupn> is
present
instead of anisole. The free bases of crystalline (p-phenolyl)dodecylguanidine-[15N2] (Figure B) and pan class="Chemical">(p-phenolyl)-ethylguanidine-[15N2]
(Figure C) each
show only two distinct resonances, at ∼50 and ∼45 ppm.
The 5 ppm chemical shift splittings likely correspond to the distinct
H-bonding patterns for the two terminal nitrogens, as seen in the
corresponding crystal structures (see below). The distinct 15N chemical shift ranges confirm the distinct chemical identities
of the free bases and their HBr salts.
Among the three pan class="Chemical">HBr saltsn>
(Figure D–F),
there are varying degrees of separation
for the two terminal −NH2 resonances. However, none
show any sign of an imino-type nitrogen with a strongly downfield-shifted pan class="Chemical">15N resonance near 110 ppm. That is, the guanidine groups in
these compounds behave almost fully protonated, as their 15N NMR spectra are generally similar to those of the corresponding
HBr salts (Figure E,F). In (p-phenolyl)ethylguanidinium-[15N2] bromide (Figure F), the 6 ppm separation likely corresponds to the
different H-bonding environment of the two −NH2 groups,
which is clearly visible in the crystal structure (see below). Somewhat
surprisingly, in the corresponding dodecyl-linked compounds (Figure E), the resonances
are superimposed (Figure E). This indicates that H-bonding environments of the two
terminal nitrogens are very similar. In the corresponding (p-anisolyl)dodecylguanidine-[15N2]
compound (Figure D), the asymmetry in the environment of the two terminal nitrogens
appears to be highest among the three HBr salts. Without the crystal
structures for either of the dodecyl-linked HBr salts (Figure D,F), we cannot provide any
detailed rationalization for these results.
All of the pan class="Chemical">HBr saltsn>
(Figure D–F)
give resonances near 54 ppm, that is, shifted
a bit downfield from the amino resonances of the corresponding free
bases (Figure A–C).
However, the chemical shift of all the −NH2 resonances
in Figure B–F,
as well as the splittings between resonances within any one molecule,
is well within the range of those seen previously seen among various
argininepan class="Chemical">salts.[25] We thus conclude that
these five spectra represent fully protonated guanidinium groups.
Only Figure A, as
well as the similar result obtained for dodecylguanidine,[12] represent good models for the 15N
spectrum of a fully deprotonated monoalkylguanidine.
X-ray Crystallography
High-resolution crystal structures
were obtained from pan class="Chemical">(p-phenolyl)-dodecylguanidinen>
pan class="Chemical">free base (Figure A), pan class="Chemical">(p-phenolyl)ethylguanidine free base (Figure B) and (p-phenolyl)guanidinium bromide (Figure C). These three compounds all crystallize
from MeOH as head-to-tail dimers in the unit cell. Both free bases
(Figure A,B) actually
crystallize as guanidinium + phenolate zwitterions. In these, but
not the HBr salt (Figure C), multiple solvent (MeOH) molecules per unit cell are clearly
included, clustered about the phenol–guanidine interaction
site. As a result of the presence of protic solvent, these zwitterion
crystals have significantly weaker individual H-bonds between guanidinium
and phenolate than are present in for dimer complexes formed in aprotic
solvents (see Discussion). Thus, the X-ray
crystal geometries are merely suggestive of the solution complexes.
Figure 13
Crystal
structures of compounds containing both phenol and guanidine
groups, showing H-bond network and distances. (A) (p-phenolyl)dodecylguanidine free base. Formula: C21 H41N3O3. Space group: P – 1. Cell length (Å): a = 6.7626, b = 7.9144, c = 21.892. Cell angles: α
= 94.407°, β = 98.175°, γ = 106.478°. R-factor (%): 4.27. (B) (p-Phenolyl)ethylguanidine
free base. Formula: C10 H17N3O2. Space group: P21/c. Cell length (Å): a = 10.1822 (5), b = 10.0597 (5), c = 11.9084 (6). Cell
angles: α = 90.00°, β = 109.8880° (10), γ
= 90.00°. R-factor (%): 3.03. (C) (p-Phenolyl)ethylguanidinium bromide. Formula: C9H14BrN3O2. Space group: P21/c. Cell length (Å): a = 11.1931 (7), b = 10.0659 (6), c = 9.8327 (6). Cell angles: α = 90.00, β = 101.508(2),
γ = 90.00. R-factor (%): 1.95.
Crystal
structures of compounds containing both n class="Chemical">pan class="Chemical">phenol and pn>an class="Chemical">guanidine
groups, showing H-bond network and distances. (A) (p-phenolyl)dodecylguanidinefree base. Formula: C21 H41N3O3. Space group: P – 1. Cell length (Å): a = 6.7626, b = 7.9144, c = 21.892. Cell angles: α
= 94.407°, β = 98.175°, γ = 106.478°. R-factor (%): 4.27. (B) (p-Phenolyl)ethylguanidinefree base. Formula: C10 H17N3O2. Space group: P21/c. Cell length (Å): a = 10.1822 (5), b = 10.0597 (5), c = 11.9084 (6). Cell
angles: α = 90.00°, β = 109.8880° (10), γ
= 90.00°. R-factor (%): 3.03. (C) (p-Phenolyl)ethylguanidinium bromide. Formula: C9H14BrN3O2. Space group: P21/c. Cell length (Å): a = 11.1931 (7), b = 10.0659 (6), c = 9.8327 (6). Cell angles: α = 90.00, β = 101.508(2),
γ = 90.00. R-factor (%): 1.95.
These structures are even imperfect ren class="Chemical">presentations
of the free-base
sampn>les that we measured using the IR spn>ectra in Kpan class="Chemical">Br pellets and solid-state
NMR, which did not include pan class="Chemical">MeOH. That is, when crystals prepared for
X-ray measurements were directly redissolved in DMSO, spectral signals
from the cocrystallized MeOH could usually be seen (see, e.g., the 13C NMR signal from MeOH in Figure 10-S-E). However, even when the X-ray crystals were thoroughly ground with
KBr to make pellets for IR measurements, no IR absorption bands from
MeOH were ever detected. It appears the MeOH was lost during KBr grinding.
Because stoichiometric MeOH could not be included in the IR samples,
acetonitrile was used as the recrystallization solvent for those samples;
and this choice was extended to solid-state NMR samples as well.
Nevertheless, these crystal structures (Figure ) supn class="Chemical">port the conclusion from solid-state 15N NMR (see above) that protons in the crystalline states
of the pn>an class="Chemical">(p-phenolyl)alkylguanidine free base compounds
are essentially fully transferred from the phenol to form guanidinium.
In all of them, the three C–N bonds of the guanidinium group
are of equal length, lying in a narrow range of 1.33–1.34 Å.
However, the free base crystal structures (Figure A,B) still show clear signs of stronger
H-bonding than is present in the HBr salt (Figure C).
In pan class="Chemical">(p-phenolyl)dodecylguanidinen>
(Figure A), the
two terminal pan class="Chemical">guanidinium
pan class="Chemical">nitrogens form H-bonds with three different O atoms. The shortest,
with an N–O distance of 2.75 Å, is with the phenolate
of the partner molecule in the head-to-tail dimer. The others are
to 2 distinct cocrystallized MeOH molecules. One MeOH accepts H-bonds
from both terminal nitrogens, of very similar lengths (2.87 and 2.95
Å). The other MeOH accepts an even stronger H-bond (2.83 Å)
from only one of the two terminal nitrogens. The two MeOH molecules
also serve as H-bond donors to the phenolateoxygen, with O–O
distances of 2.36 and 2.62 Å. Somewhat surprisingly, the intrachain
nitrogen does not participate in any H-bonding in this crystal.
By contrast, in pan class="Chemical">(p-phenolyl)ethylguanidinen> pan class="Chemical">free
base (Figure B),
all of the three pan class="Chemical">nitrogens participate in H-bonding. There are four
H-bond acceptor groups for each guanidinium. Three of these are distinct
MeOHoxygens, accepting H-bonds from the terminal nitrogens, with
the O–N bond distances in the narrow range of 2.85–2.90
Å. The fourth H-bond acceptor, phenolate, shows an H-bonding
pattern that is strikingly different from that in the dodecyl-linked
compound (Figure A), in that it accepts moderately strong H-bonds of nearly equal
length (2.87 and 2.88 Å) from two nitrogens of the same guanidinium
group, one being the intrachain nitrogen. The phenolate in the ethyl-linked
compound (Figure B) simultaneously accepts even stronger H-bonds from two other groups,
as measured by their short lengths. One is MeOH (O–O distance
of 1.69 Å); the other is a guanidinenitrogen from a neighboring
dimer (O–N distance of 2.85 Å).
One desired control
compound, the n class="Chemical">pan class="Chemical">HBr salt of the dodecyl-linked
compn>ound in Figure A, has not yet been amenable to crystallography. However, the HBrsalt of the corresponding pan class="Chemical">ethyl-linked compound (Figure C) again clearly demonstrates
head-to-tail dimer formation. In this case, Br– counterions
are present, and no cocrystallized MeOH molecules. This results in
somewhat different H-bonding patterns. As with the corresponding free
base (Figure B),
the intrachain nitrogen and one of the terminal nitrogens still form
equal-length H-bonds with the phenolic oxygen of the dimer partner
molecule. However, in the HBr salt (Figure C), the N–O bond distances are both
larger than any N–O distances in the free bases. This indicates
weaker H-bonding between guanidinium and protonated phenol, as compared
to guanidinium + phenolate. Furthermore, as a result of complete proton
transfer to the phenol group, the phenolic C–O bond length
is also clearly 0.2–0.4 Å longer in the HBr salt than
in the free bases.
Discussion
Phenol and Guanidine Free-Base
Groups Form an Unusually Strong
and Polarizable H-Bond in Nonpolar Environments
Computations
in vacuo (Figure )
and UV sn class="Chemical">pectroscopic experiments in hexane (Figure ) demonstrate that phenol and guanidine groups
form an overall-neutral H-bonded complex in nonpolar environments,
with a binding ΔE or ΔH of −70 kJ mol–1. The closeness of the agreement
between the computational and experimental numbers is fortuitous,
because the in vacuo computation did not take into account changes
in solvent interfacial interactions that could easily contribute up
to ±10–20 kJ mol–1 to the binding energy.
Furthermore, Figure presents only the single H-bonded structure with the very lowest
energy. Several other local minima (e.g., involving H-bonding with
a deprotonated intrachain nitrogen) were only a few kJ mol–1 higher in energy and would likely be represented in the physiological-temperature
ensemble measured in the experimental work. Such alternative structures
could also easily play a role inside proteins.
The experimental
UV sn class="Chemical">pectrum of the H-bonded complex in hexane (Figure C) was compn>uted by a least-squares fitting
procedure, based on the assumpn>tion of a simple compn>lexation equilibrium
that could be shifted by dilution or heating. This UV spectrum shows
distinct features attributable to both phenol (∼280 nm) and
phenolate (∼310 nm) structures, indicating a rapid protomeric
equilipan class="Chemical">brium within the complex corresponding to the structures in Scheme B,C. There is no
certain way to extract the individual extinction coefficients and
individual concentrations of individual protomers. However, in hexane,
they appear to provide roughly equal mole fractions on the basis of
an assumption of roughly equal integrated extinction coefficients
for the phenol (280 nm) and phenolate (310 nm) bands.
Even in
a more polar an class="Chemical">protic solvent (DMSO), H-bonding is detectable
in UV spectroscopn>ic measurements, which show concentration-dependent
formation of (p-phenolyl)alkylguanidine dimer complexes
(Figure ; additional
data in Supporting Information). The experimental
values of ΔG°298 = −4
kJ mol–1 and ΔH° = −11
kJ mol–1 per phenol–pan class="Chemical">guanidine H-bonding pair
are considerably smaller in DMSO than in hexane. This is attributable
to the much stronger solvating power of DMSO than of hexane. As a
result of solvation, formation of the dimer is expected to disrupt
numerous energetically favorable interactions between DMSO molecules
and each (p-phenolyl)ethylguanidine monomer. These
include interactions in which DMSOoxygens act as H-bond acceptors
from phenol, as well as other dipole–dipole interactions.
The UV spectrum of n class="Chemical">pan class="Chemical">(p-phenolyl)ethylguanidine
dimer in pn>an class="Chemical">DMSO was also estimated (see Supporting Information). In this case, the phenolate-like absorbance band,
with a maximum at ∼308 nm, was clearly predominant, and only
much smaller features near 280 nm might be attributable to a contribution
from phenol. Thus, shifting from hexane to a more polar solvent such
as DMSO is sufficient to favor the zwitterionic protomer within the
H-bonded dimer. Even taken by themselves, these UV spectral results
demonstrate that the H-bond between phenol and guanidine is indeed
highly polarizable. That is, the favored protomer can easily be changed
by alterations in solvent environment. With widely differing dielectric
constants, DMSO and hexane produce big changes in proton polarization.
The existence of pan class="Chemical">phenoln>–pan class="Chemical">guanidine complexes exhibiting rapid
protomeric equilipan class="Chemical">bria and large proton polarizability is also supported
by IR spectroscopy in aprotic solvents. In both CCl4 (Figure ) and DMSO (Figure ), the phenol/guanidine
H-bonded complexes simultaneously exhibit characteristic peaks of
protonated and deprotonated phenol groups, as well as protonated and
deprotonated monoalkylguanidine groups. In DMSO (Figure ), the IR bands representative
of the zwitterionic limiting structure (Scheme B) are more favored. However, in the less-polar
CCl4 (Figure ), there are nearly equal contributions of both forms (Scheme B,C).
Mpan class="Chemical">ethyln>ating the
pan class="Chemical">phenol, or protonating pan class="Chemical">guanidine with a strong
acid such as HBr, eliminates the possibility of this type of strong
H-bonding interaction. Consistent with this, in crystals or in DMSO,
the IR and NMR signals from the guanidine group of (p-anisolyl)dodecylguanidinefree base are nearly the same as those
from dodecylguanidinefree base. The same similarities apply to the
corresponding HBr salts, (p-anisolyl)dodecylguanidinium-Br
and dodecylguanidinium bromide. In the latter case, the similarities
also extend to (p-phenolyl)dodecylguanidinium bromide.
This shows that a protonated (neutral) phenol group interacts only
weakly with guanidinium. Spectral features of phenol are also reciprocally
unaffected by the presence of guanidinium (Figure ; see also Supporting Information). These control experiments show that the O of
anisole or neutral phenol serves as only a very weak
H-bond acceptor from guanidinium (and also, as shown
in the case of anisole, from guanidinefree base).
Phenolate Spectral
Properties Are Strongly Perturbed When It
Accepts an H-Bond from Guanidinium, but This Is Only Seen in Aprotic
Environments
Strong perturbation of the pan class="Chemical">phenolate group by
H-bonding to a pan class="Chemical">guanidinium counterion in aprotic solvents is detectable
with UV, IR, and NMR spectroscopy.
UV
The pan class="Chemical">phenolaten>
component of the H-bonded compn>lex
(Scheme B) appears
to have λmax near 310 nm in hexane (Figure ) and double maxima at 293
and 310 nm in DMSO (see Figure 4-S). In
contrast, cresolate salts in protic solvents show a single broad peak
with λmax near ∼296 nm. This value is seen
for MeOH in Figure F–H, and is nearly the same as that reported for pan class="Chemical">tyrosinate
in water at high pH (λmax ≈ 294 nm).[18] That is, when combined in protic solvents, phenol
and guanidine simply form solvated phenolate and guanidinium ions
in the familiar Brønsted–Lowry type of acid–base
equilibrium. Even at concentrations exceeding ∼50 mM (Figure ), there is no evidence
of significant formation of H-bonded complexes between these ions
in protic solvents.
The neutral (pan class="Chemical">phenoln>) species within the
very same H-bonded compn>lexes serves as a control. It shows considerably
less solvent-depn>endent variability for λmax, which
appears at 280–283 nm for both protic and aprotic solvents
(Figures and 4). Only the vipan class="Chemical">bronic bandshape changes, with considerable
sharpening as the solvent polarity is reduced.
IR
A clear IR demonstration of H-bonding perturbation
on the n class="Chemical">pan class="Chemical">phenolate is seen for the strongest C=C stretch mode.
For pn>an class="Chemical">(p-phenolyl)dodecylguanidine in DMSO, at a concentration
expected to produce mainly the dimer complex, it is at 1494 cm–1 (Figure A), compared with 1501 cm–1 in MeOH (Figure B). The latter frequency
is also observed in tyrosine solutions at high pH.[18] This same vibration is at an intermediate frequency of
1498 cm–1 both in CCl4 (Figure ) and in a solid state measurement
(Figure ).
In
all these spectra, the n class="Chemical">pan class="Chemical">phenolate pn>eak at 1494–1501 cm–1 consistently has a lpan class="Chemical">arger half-width than the pan class="Chemical">phenol peak at ∼1515
cm–1, likely due to inhomogeneous broadening. This
is consistent with a greater degree of sensitivity of the phenolate
vibrational frequency to the solvent environment as compared to phenol
itself.
An even greater variability of pan class="Chemical">phenolaten> vipan class="Chemical">brational
frequencies
is seen for the C–O stretch. This vipan class="Chemical">bration ranges from 1264
to 1309 cm–1 depending on the solvent and counterion
environments (Figure D–G). The lowest value, for (p-phenolyl)ethylguanidine
in MeOH, is similar to the value of 1270 cm–1 reported
for aqueous tyrosine and tyrosine-containing peptides at a high pH.[18] The highest value of 1309 cm–1, for (p-phenolyl)ethylguanidine in DMSO, is probably
related to the unusual H-bonding environment in the dimer complex.
Once again, the corresponding C–O vibration in phenol itself
is more constant (Figure A–C). It varies in intensity, sometimes resulting in
the inability to resolve it from a nearby peak at 1235 cm–1, but appears to fall in a narrow (∼3-cm–1) range around 1246 cm–1. The latter is its value
in MeOH (Figure C)
and in water-solvated tyrosine.[18] In summary,
in the IR as in the UV region, phenolate spectral bands are much more
sensitive to solvent environment than the corresponding features of
phenol.
NMR
The pan class="Chemical">phenoln>/pan class="Chemical">phenolate ring pan class="Chemical">1H resonances
in DMSO show a consistent dependence on the presence and concentration
of guanidine/guanidinium H-bonding partners, but they always remain
far from the values seen in potassium p-cresolate.
This can be seen most clearly in Figure . The asymptotic values at low (<5 mM)
concentration of (p-phenolyl)dodecylguanidine match
those of p-cresol. The asymptotic values at high
(>100 mM) concentration shift barely a quarter of the way to those
seen for potassium p-cresolate in DMSO. Similarly,
the spectra of dodecylguanidinium-phenolate salts in DMSO show ring 1H resonances (see Figures 9-S-B,C in Supporting Information) that are much closer to those of p-cresol than of potassium p-cresolate.
This is seen with a dodecylguanidine/phenol ratio of either 1:1 or
1.6:1, either of which would be expected to fully deprotonate 0.06
M p-cresol in DMSO if Brønsted–Lowry
acid–base equilibrium was present. That is, the 1H NMR spectra of phenolate when strongly H-bonded to a guanidinium
counterion largely resemble those of p-cresol itself,
and are very different than with a weakly H-bonded K+ counterion.
Similar conclusions apply to the 13C resonance of the phenolate
ring in aprotic solvents. With a guanidinium counterion, these phenolatesalts remain closer to values seen for p-cresol,
than to the corresponding potassium salt (see Figure ).
Proton Polarizability
Based on UV, IR, and NMR spectroscon class="Chemical">pic
measurements, the proton in the phenol–guanidine H-bond rapidly
shifts position between the two bases (guanidine and phenolate), and
is highly polarizable. That is, the fraction of time it spends on
each of the two groups is strongly dependent on the solvent environment,
with the zwitterionic form being favored in more polar solvents. Additionally,
IR spectra of the H-bonded complexes in aprotic solvents (Figure ), as well as in
crystalline state (Figures and 6-S-G) display a strong broad
absorption band between 3500 and 2500 cm–1. Such
a continuum absorption is characteristic of
a polarizable H-bond. It is notably absent from the control spectra
of the HBr salts of the same compounds in crystalline state (also
shown in Figures and 6-S-G in Supporting Information). This confirms
that the phenol–guanidinium H-bonding interaction is much weaker
than for phenolate–guanidinium.
The different degrees
of polarization in various solvent environments can be correlated
with solvent-den class="Chemical">pendent pKa values of the
phenol and guanidine groups. In aprotic solvents less polar than DMSO,
the acidity of guanidinium is expected to increase, whereas that of
phenol is expected to decrease, effectively bringing the pKa values of guanidine and phenol closer together.
This promotes the formation of an even stronger H-bond, with an even
greater proton polarizability, because H-bonding is generally strongest
when the pKa of the donor is closest to
that of (the conjugate-acid form of) the acceptor.[22] That is, the proton affinities of the groups competing
for the H+ must be large to get the strongest H-bonds.
This is precisely the case for guanidine and phenolate in our H-bonded
complexes, most particularly in nonpolar solvents such as hexane and
CCl4. A much weaker interaction is expected when guanidinium
serves as the H-bond donor and phenol serves as the H-bond acceptor
because of the much lower proton affinity of (neutral) phenol compared
with phenolate.
No single parameter such as dielectric constant
is a n class="Chemical">perfect predictor
of solvent-dependent pKa, nor of whether
phenol and guanidine are likely to participate in the H-bond interactions.
Indeed, when considering such H-bonding in crystals and other solid
states, including the interiors of any particular protein or lipid
bilayer, only wide ranges are available for the dielectric constant.
Furthermore, trying to predict the interactions of a highly directional
and localized H-bond with its environment by using a continuum model
is intrinsically unreliable. Therefore, we cannot claim to be modeling
the arg–tyr system inside a protein or membrane with great
precision by sampling their side chains’ behavior in only four
homogeneous solvents. Observing our model compounds in a wider range
of solvent environments will likely reveal interesting additional
aspects of their interactions.
Additionally, the concept of
n class="Chemical">pKa in
nonaqueous environments is not rigorous and depends strongly on the
particular counterions present. This lack of rigor explains some experimentally
contradictory behavior. In water, guanidinium is unquestionably a
weaker Brønsted–Lowry acid than phenol. However, we have
shown quite directly (Figure F) that in MeOH, guanidinium is the stronger acid by about
1.3 pKa units. This contradicts the much
lower published pKa value for p-cresol than for guanidinium in DMSO (18.5[23] and 28.5,[24] respectively). These
published values suggest that guanidinefree base should be capable
of completely deprotonating phenol in either MeOH or DMSO. Besides
showing an opposite Brønsted–Lowry behavior in MeOH (Figure F), our spectra measured
in DMSO (e.g., Figure B–D) indicate that Brønsted–Lowry acid–base
equilibria do not occur with this pair of solutes. Instead, Lewis
acid–base complexation predominates. Within such a complex,
it is not easy to assign pKa values to
the individual components.
Moderately Strong H-Bonding between Guanidine
and Phenol Groups
Is Retained in the Crystalline State
So far, we have only
been able to obtain crystal structures of the zwitterion (Scheme B) in the presence
of cocrystallized n class="Chemical">pan class="Chemical">MeOH. Even in the absence of pn>an class="Chemical">MeOH, the IR spectra
of KBr pellets (Figure ) support the conclusion that the zwitterion predominates in the
solid state. These IR spectra nevertheless demonstrate that even this
zwitterion form (predominantly Scheme B) retains a more polarizable H-bond than the phenol
+ guanidinium salt form (Scheme A), as evidenced by the stronger continuum absorption
band (compare Figure green and purple spectra). Thus, the phenolate + guanidinium H-bond
remains much stronger than the phenol + guanidinium H-bond. This conclusion
also finds support in the crystal structures, which show guanidinium-phenolate
H-bonds that are shorter (2.75–2.9 Å, Figure A,B), than for guanidinium
+ phenol (2.96–3.00 Å; Figure C). The value of 2.75 Å is indicative
of quite a strong H-bonding. Indeed, it is close to the range (<2.65
Å) for very strong N–H···O bonds that have
been interpreted as having partial covalent character.[26]
Interpretation of 15N NMR Measurements
The
solid-state 15N chemical shifts of the terminal nitrogens
of (p-phenolyl)alkylguanidines (Figure B,C) are all within the range
of 44–51 ppm. These are similar to those of the corresponding
HBr salts (Figure E,F) and are also well within the 30–70 ppm range previously
observed in 15N NMR chemical shifts of arginine and other
guanidinium salts in solutions,[17] as well
as in the solid state.[25] These results
are consistent with the conclusion from X-ray crystallography and
IR spectra that the zwitterion, that is, the guanidinium-phenolate
protomer, predominates in the crystals of (p-phenolyl)alkylguanidinefree bases.In contrast, actual stoichiometric deprotonation
of the n class="Chemical">pan class="Chemical">guanidine groupn> produces changes far outside the range for
guanidinium. Thus, the ∼60 ppm separation between amino and
pan class="Chemical">imino 15N chemical shifts for crystalline (p-anisolyl)dodecylguanidinefree base (Figure A) falls on a continuum with values previously
published for alkylguanidinefree bases in various solvent environments.
The least H-bonded of these are the penta-alkylated guanidinefree
bases, which in aprotic solvents completely lack H-bond donors. In
such solvents, the chemical shift splitting between the imino and
amino nitrogens is very large, typically near 150 ppm.[27] These splittings can be observed directly in
the solution-state measurements in CHCl3 because there
are no rapid H+ transfers that can make the terminal amino
and imino nitrogens chemically equivalent. The same prior work showed
that when even a single N–H bond is introduced, that is, in
tetra-alkylated guanidinefree bases in CHCl3, the 15N chemical shift splitting drops to 126 ppm. As we have now
shown, with four N–H bonds per guanidine, specifically in monoalkylated
guanidinefree bases such as dodecylguanidine, the splitting drops
further to 60 ppm. However, due to rapid tautomerization in a solution,
this chemical shift splitting can only be observed with the solid-state
measurements (see ref (12) and Figure A).
We conclude that lpan class="Chemical">argn>e variations in pan class="Chemical">15N chemical shift
splittings in deprotonated pan class="Chemical">guanidines are produced by the presence
or absence of H-bond donation to the imino nitrogen, mainly from neighboring
amino nitrogens within the same alkylguanidine molecule. The presence
of such H-bonding renders the chemical environment of the imino nitrogen
more “amino-like”. Additionally, in the presence of
a sufficiently strong proton donor(s) external to the guanidine group,
the imino group becomes completely or nearly completely protonated,
that is, nearly indistinguishable from the nearby amino group, leading
to small splittings such as in Figure B,C.
The availability of internal
and external H-bond pan class="Species">donorn>(s) to the
pan class="Chemical">imino nitrogen thus can account for most of the enormous variation
in pan class="Chemical">15N chemical-shift splittings seen in guanidine-free
bases. The weaker the overall H+ donation to the imino
group, the greater the splitting. If this hypothesis is true, then
it should be possible to systematically adjust the H+-donating
strength of the environment to obtain just about any value of splitting
for the η-nitrogens of arginine, from 0 to 60 ppm. NMR spectroscopic
modeling of possible arginine environments in biology cannot likely
be considered complete until a good range of intermediate values have
been observed experimentally.
In the pan class="Chemical">(p-phenolyl)alkylguanidinen>
model system
we have developed, nearly compn>lete H+ transfer prevails,
and the splittings of only 5–6 ppm are obtained (Figure B,C). However,
making the environment around the neutral phenol–guanidine
H-bond less polar, so as to emphasize Scheme C more, should result in less H+ transfer to the imino group, and thus to greater splitting. The
simplest way to accomplish this, while still allowing X-ray crystallography
and solid-state 15N measurements, might be to crystallize
our (p-phenolyl)alkylguanidine model compounds without
including MeOH. Alternatively, it should be possible to resolve a
large splitting of pan class="Chemical">15N chemical shifts of the two terminal
nitrogens of the H-bonded dimer complex in DMSO, either in a liquid
state or as a glassy solid. Our attempts at these measurements have
not been successful to date (see Supporting Information for 15N NMR measurements in DMSO at 25 and 50 °C).
Possible New Motif for Membrane Protein Structures
The 70
kJ mol–1 binding enthalpy of the (n class="Chemical">pan class="Chemical">phenol
+ pn>an class="Chemical">guanidine) complex is quite large—nearly double that of (phenol
+ alkylamine) complexes, which have been measured by direct calorimetry,[28] as well as by UV measurements similar to those
shown in Figure above.[13,14]
H-bonded dyads (e.g., pan class="Chemical">argn>–pan class="Chemical">asp pairs with closest contacts
commonly below 3.0 Å) have long been an identifiable motif in
crystal structures, as shown in a statistical anapan class="Chemical">lysis of amino acid
geometries in such structures.[29] In this
publication, no common H-bonded tyr–arg patterns were clearly
discernable; however, no membrane proteins were included among the
62 proteins analyzed, as it was limited to structures determined prior
to 1988. The database used for a more recently updated online version
of this catalog[30] now includes ∼40×
as many structures. This probably includes a number of membrane proteins;
however, the database is not organized according to protein type or
localization. This recent database shows recurring geometric motifs
that appear to be suitable for H-bonding not only for dyads of tyr–lys
(“Cluster 6”) but also for tyr–arg (“Cluster
1” and “Cluster 5”). The H-bonded tyr–arg
grouping in Figure appears possibly to fit in the latter Cluster 5.
Still, few
pan class="Chemical">tyrn>–pan class="Chemical">arg environments have been clearly identified
in proteins that are sufficiently aprotic and nonpolar, for example,
as in Mpan class="Gene">eta II (Figure ) that they might, based on our modeling, support deprotonation of
arg; certainly not as readily as deprotonation of lys. Nevertheless,
we hypothesize, based on our modeling, that the overall neutral form
of the (arg + tyr) dyad may play a significant role in a small number
of protein structures, specifically in very nonpolar aprotic environments.
We designate this hypothetical grouping a Bonded Uncharged (aRginine
+ tYrosine) or BU(RY) motif, signifying any strongly H-bonded arg–tyr
dyad that is overall neutral, that is, in zwitterionic (Scheme B), fully neutral (Scheme C), or intermediate
protomeric states.
Our experimental modeling clearly shows this
dyad is not likely
allowed in n class="Chemical">protic environments. Thus, a functional BU(RY) grouping
resembling our model compounds will probably occur only in biological
environments deficient in H-bonding groups, predominantly membranes,
where no proton accepn>tor for the tyrosine-OH is available, besides
a deprotonated arginine; and no proton donor is available for that
deprotonated arginine, besides the tyrosine-OH. A specific example
of a sufficiently nonpolar aprotic environment is provided in the
region surrounding the conserved R–Y grouping in Meta II (Figure ).[4] Conserved adjacent RY groupings in protein primary sequences
are quite common in small-peptide hormones, such as members of the
neuropeptide Y family. In fact, the powerful endogenous brain analgesic
kyotorphin is simply Y–R dipeptide. The fact that kyotorphin
and other RY-rich hormones tend to act at or within membranes suggests
that the formation of a BU(RY) dyads within the nonpolar interior
of a membrane bilayer should be considered as a possible physiological
structure for such peptides.
Other candidates for BU(RY) motifs
are in the activated intermediate
states of heptahelical memn class="Chemical">pan class="Chemical">brane pn>roteins. X-ray crystallography shows
a conserved arg135–tyr223 dyad in an aprotic environment in
Mpan class="Gene">eta II. When the G-protein transducin reversibly dissociates from
this activated receptor (Figure A), there is no clear evidence of any counteranion
for arg135 remaining behind or entering the binding pocket to substitute
for transducin’s C-terminal carboxylate. In fact, the only
protic or oxygen-containing group within 5 Å of the arg135guanidine
group is tyr227. Specifically, no structural water molecules are observable
within this distance. Instead, nearby hydrocarbon side chains almost
completely fill the space surrounding arg135 and tyr223. Their H-bonding
distance of 2.7 Å matches well that seen in our best crystallographic
model for a BU(RY) grouping, (p-phenolyl)dodecylguanidine
dimer (2.75 Å, Figure ).
Solid-state NMR measurements have previously beenn class="Chemical">performed, which
detected a clear perturbation of tyr227 during Meta II formation,
but seemed to exclude the possibility that either it or arg135 underwent
deprotonation.[10] However, it may not be
safe to rely on expectations based on prior model compounds,[10,25,27,31,32] which did not include any with a guanidine–phenol
pairing. In fact, (p-phenolyl)dodecylguanidine and
(p-phenolyl)ethylguanidine in aprotic environments
deviate significantly in their spectroscopic behavior from all of
the earlier model compounds. Specifically, based only on earlier model
compounds, deprotonation of tyr227 would be expected to result in
a ∼15 ppm downfield shift for the 13C resonance
of its ipsocarbon.[10] However, our new
model compounds for such a BU(RY) dyad show 13C and 1H chemical shift values for the phenolic ring that are substantially
closer to those of neutral p-cresol than potassium p-cresolate, while still showing the 13C and 1H chemical shifts closer to protonated guanidinium than to
guanidine (see Figures and 11). That is, when strongly H-bonded,
both guanidine and phenol can appear spectroscopically to be mostly
protonated, even though overall they are neutral, that is, there is
a full stoichiometric deprotonation of the combined grouping.
Rather than NMR or IR, the best spectroscon class="Chemical">pic method to detect
a BU(RY) grouping in a membrane protein may be UV absorpn>tion, using
the spn>ectral marker band for the pan class="Chemical">phenolate component near 310–315
nm (Figures and 4), or possibly UV resonance Raman. Unfortunately,
the strong UV–visible absorption bands of pan class="Chemical">retinylidene chromophores
of rhodopsin and Meta II overlap this 310 nm region. As a result,
Meta II-minus-rhodopsin UV–visible difference spectra provide
no clear evidence, or lack of evidence, for a 310 nm difference band
attributable to a BU(RY) dyad.[33] However,
substitution of a retinal analogue such as retinal2 could
shift the absorption spectrum of both rhodopsin and Meta II chromophores
to longer wavelengths, reducing spectral interference in the 310 nm
region. Alternatively, many other GPCRs and their agonists have no
competing absorptions in this region, so UV difference spectroscopy
should provide a rather simple and unambiguous test of whether a BU(RY)
dyad can form upon addition of agonist to receptor.
Materials and
Methods
UV–Visible Spectra
Measurements were performed
on a Shimadzu UV-265 dual-beam spectrometer, using 0.5 or 1 cm pathlength
fused pan class="Chemical">silica cuvettes as indicated.
IR Spectra
Spectra
were obtained using a Nicolet Magna
IR 860 sn class="Chemical">pectrometer. For solution spectra, a refillable variable pathlength
cell with 15 μm spacers and 2 mm thick pan class="Gene">CaF2 windows
was used. For crystalline sampn>les, the sampn>les were prepn>ared as Kpn>an class="Chemical">Br
pellets.
UV/Visible/IR Spectral Manipulations and Fitting
Spectral
manin class="Chemical">pulations (e.g., water vapor subtraction, smoothing, pn>eak fitting)
were performed using GRAMS32/AI software, excepn>t for SVD, which was
performed by transferring the data as ASCII files and then manipulating
them as matrices in Matlab. The SVD routine in Matlab was used to
compute the best-fit spectrum of the 1:1 complex formed between p-cresol and pan class="Chemical">n-dodecylguanidine in hexane,
at varying dilutions of an initial total concentration of ∼700
μM; as well as the best-fit concentration of this species. The
fitting was done in a least-squares fashion by fitting a set of 11
measured spectra of the sample at different concentrations and temperatures
to weighted sums of two components: the separately measured spectrum
of pure p-cresol in hexane and a single additional
spectral component, which was allowed to vary. Further details are
provided in Results and Supporting Information.
Computational Modeling
All computations were done by
using GAUSSIAN09 imn class="Chemical">plemented with the GAUSSVIEW 4.1 visualization
interface. DFT energy computations used the B3LYP functional with
6-31G++ or 6-311G++ basis set.
Solution-State 1H and 13C NMR Spectra
Except where another solvent
was indicated, all of the n class="Chemical">pan class="Chemical">1H and pn>an class="Chemical">13C NMR spectra
were recorded in DMSO-d6 using a 400 MHz
Bruker NMR spectrometer. Chemical shifts
are presented in ppm and are referenced to the residual DMSO peak
at 2.50 ppm for 1H and to the natural-abundance 13C-DMSO peak at 39.5 ppm peak.
15N Solid-State
NMR Spectra
All solid-state pan class="Chemical">15Nn> NMR cross-polarization
magic-angle spinning (CP-MAS) measurements
were performed at the Analytical and Technical services at the SUNY
ESF using a 400 MHz Bruker spectrometer. The spinning rate was 4000
Hz, with a delay between pulse 5 s, 400 ppm sweep width (12 165.450
Hz), and 1H decoupling field of 47 200 Hz. Number
of data points was 484, acquisition time was 0.02 s; processed with
a line pan class="Chemical">broadening of 40 Hz reference to ammonium-15N2 sulfate standard at 0 ppm.
Sources of Chemicals
pan class="Chemical">p-Cresol, as
well as all of the reagents and solvents used in the synthesis and
anapan class="Chemical">lysis of other compounds were obtained from Sigma-Aldrich and used
directly unless otherwise indicated. The solvents were dried by the
addition of molecular sieves. Kpan class="Chemical">Br for solid-state IR samples was dried
at 110 °C over several days before use.
Syntheses of Model Compounds
Dpan class="Gene">etan>iled procedures for
synthesis of p-phenol-dodecylguanidinium bromide
and dodecylguanidinium bromide are described elsewhere.[12] Analogous procedures were followed to prepare pan class="Chemical">p-phenol-ethylguanidinium bromide, using commercially available p-phenol-ethylamine, that is, tyramine, as the starting
point and converting the amine to a guanidine by reaction with S-methyl thiourea hydroiodide, just as was done to p-phenol-dodecylamine in the previously described synthesis
of p-phenol-dodecylguanidinium bromide.[12] The procedure for incorporating enriched (∼50
atom %) 15N at the guanidine terminal nitrogens, in this
case to aid in obtaining solid-state 15N NMR spectra, was
also described in the same publication.
Deprotonation of Alkylguanidine
Compounds
A general
method for deprotonating n class="Chemical">pan class="Chemical">alkylguanidinium bromides to form pn>an class="Chemical">alkylguanidine
free bases is described in the accompanying paper.[12] This procedure reproducibly gave crystalline samples, with
elemental analyses that match the predicted values for the corresponding
free bases. In brief, a dry crystalline alkylguanidinium bromide (or p-phenol alkylguanidinium bromide) was first recrystallized
from acetonitrile and then redissolved in dry methanol (MeOH). Potassium tert-butoxide (1 mol equiv), also dissolved in dry methanol
just before use, was added into the solution dropwise, causing the
solution to become milky due to the precipitation of KBr. After the
precipitation had stopped, the solution was heated briefly to boiling,
then centrifuged while still hot (2 min, 3000g).
The clear supernatant was transferred to a clean test tube and then
concentrated ∼20-fold under a gentle stream of dry N2. The evaporative cooling aided in promoting a high yield of crystals.
These were harvested while still cold by pipetting off the bulk of
residual tert-butanol and MeOH and then drying under
N2.
Sample Crystallization and X-ray Crystallography
Crystals
of the free bases of (p-phenolyl)-dodecylguanidine
and (p-phenolyl)ethylguanidine, as well as their
HBr salts, were obtained by slowly evaporating and cooling concentrated
solutions of them in dry MeOH, under a gentle stream of room-temperature
N2. Making a sufficiently concentrated solution of (p-phenolyl)dodecylguanidinefree base in dry MeOH required
gentle heating, whereas the other compounds dissolved readily at room
temperature. The resulting ∼0.1 mm size crystals were stored
for up to ∼1 week in a minimal volume of MeOH at 4 °C
before mounting in the diffractometer. Structures were recorded at
90 K.Crystals of (p-phenolyl)-class="Chemical">n>an class="Chemical">dodecylguanidinium
bromide, (p-anisolyl)-dodecylguanidinefree base
and its HBr salt, as well as dodecylguanidinefree base and its HBrsalt could also be made by the approach described above, but were
not of sufficient quality to produce good diffraction.