Interference between the hydrogen bonds to the two rings of nicotine

The biochemical transport and binding of nicotine depends on the hydrogen bonding between water and binding site residues to the pyridine ring and the protonated pyrrolidinium ring. To test the independence of these two moderately separated hydrogen-bonding sites, we have calculated the structures of clusters of protonated nicotine with water and a bicarbonate anion, benzene, indole, or a second water molecule. Unprotonated nicotine-water clusters have also been studied for contrast. The potential energy surfaces are first explored with an intermolecular anisotropic atom-atom model potential. Full geometry optimizations are then carried out using density functional theory to include nonadditive terms in the interaction energies. The presence of the charge on the pyrrolidine nitrogen removes the conventional hydrogen-bonding site on the pyridine ring. The hydrogen-bond ability of this site is nearly recovered when the protonated pyrrolidinium ring is bound to a bicarbonate anion, whereas its interaction with benzene shows a much smaller effect. Indole appears to partially restore the hydrogen-bond ability of the pyridine nitrogen, although indole and benzene both pi-bond to the pyrrolidinium ring. A second hydrogen-bonding water produces a significant conformational distortion of the nicotine. This demonstrates the limitations of the conventional qualitative predictions of hydrogen bonding based on the independence of molecular fragments. It also provides benchmarks for the development of atomistic modeling of biochemical systems.     

Strong solute-solute dispersive interactions in a protein-ligand complex

The contributions of solute-solute dispersion interactions to binding thermodynamics have generally been thought to be small, due to the surmised equality between solute-solvent dispersion interactions prior to the interaction versus solute-solute dispersion interactions following the interaction. The thermodynamics of binding of primary alcohols to the major urinary protein (MUP-I) indicate that this general assumption is not justified. The enthalpy of binding becomes more favorable with increasing chain length, whereas the entropy of binding becomes less favorable, both parameters showing a linear dependence. Despite the hydrophobicity of the interacting species, these data show that binding is not dominated by the classical hydrophobic effect, but can be attributed to favorable ligand-protein dispersion interactions.     

Automated linkage of proteins and payloads producing monodisperse conjugates

Controlled protein functionalization holds great promise for a wide variety of applications. However, despite intensive research, the stoichiometry of the functionalization reaction remains difficult to control due to the inherent stochasticity of the conjugation process. Classical approaches that exploit peculiar structural features of specific protein substrates, or introduce reactive handles via mutagenesis, are by essence limited in scope or require substantial protein reengineering. We herein present equimolar native chemical tagging (ENACT), which precisely controls the stoichiometry of inherently random conjugation reactions by combining iterative low-conversion chemical modification, process automation, and bioorthogonal trans-tagging. We discuss the broad applicability of this conjugation process to a variety of protein substrates and payloads.

Controlled protein functionalization holds great promise for a wide variety of applications.

Applications of protein conjugates are limitless, including imaging, diagnostics, drug delivery, and sensing.^1–4 In many of these applications, it is crucial that the conjugates are homogeneous.⁵ The site-selectivity of the conjugation process and the number of functional labels per biomolecule, known as the degree of conjugation (DoC), are crucial parameters that define the composition of the obtained products and are often the limiting factors to achieving adequate performance of the conjugates. For instance, immuno-PCR, an extremely sensitive detection technique, requires rigorous control of the average number of oligonucleotide labels per biomolecule (its DoC) in order to achieve high sensitivity.⁶ In optical imaging, the performance of many super-resolution microscopy techniques is directly defined by the DoC of fluorescent tags.⁷ For therapeutics, an even more striking example is provided by antibody–drug conjugates, which are prescribed for the treatment of an increasing range of cancer indications.⁸ A growing body of evidence from clinical trials indicates that bioconjugation parameters, DoC and DoC distribution, directly influence the therapeutic index of these targeted agents and hence must be tightly controlled.⁹Standard bioconjugation techniques, which rely on nucleophile–electrophile reactions, result in a broad distribution of different DoC species (Fig. 1a), which have different biophysical parameters, and consequently different functional properties.¹⁰Open in a separate windowFig. 1Schematic representation of the types of protein conjugates.To address this key issue and achieve better DoC selectivity, a number of site-specific conjugation approaches have been developed (Fig. 1b). These techniques rely on protein engineering for the introduction of specific motifs (e.g., free cysteines,¹¹ selenocysteines,¹² non-natural amino acids,^13,14 peptide tags recognized by specific enzymes^15,16) with distinct reactivity compared to the reactivity of the amino acids present in the native protein. These motifs are used to simultaneously control the DoC (via chemo-selective reactions) and the site of payload attachment. Both parameters are known to influence the biological and biophysical parameters of the conjugates,¹¹ but so far there has been no way of evaluating their impact separately.The influence of DoC is more straightforward, with a lower DoC allowing the minimization of the influence of payload conjugation on the properties of the protein substrate. The lowest DoC that can be achieved for an individual conjugate is 1 (corresponding to one payload attached per biomolecule). It is noteworthy that DoC 1 is often difficult to achieve through site-specific conjugation techniques due to the symmetry of many protein substrates (e.g., antibodies). Site selection is a more intricate process, which usually relies on a systematic screening of conjugation sites for some specific criteria, such as stability or reactivity.¹⁷Herein, we introduce a method of accessing an entirely new class of protein conjugates with multiple conjugation sites but strictly homogenous DoCs (Fig. 1c). To achieve this, we combined (a) iterative low conversion chemical modification, (b) process automation, and (c) bioorthogonal trans-tagging in one workflow.The method has been exemplified for protein substrates, but it is applicable to virtually any native bio-macromolecule and payload. Importantly, this method allows for the first time the disentangling of the effects of homogeneous DoC and site-specificity on conjugate properties, which is especially intriguing in the light of recent publications revealing the complexity of the interplay between payload conjugation sites and DoC for in vivo efficacy of therapeutic bioconjugates.¹⁸ Finally, it is noteworthy that this method can be readily combined with an emerging class of site-selective bioconjugation reagents to produce site-specific DoC 1 conjugates, thus further expanding their potential for biotechnology applications.¹⁹     

Dipole-bound anions of highly polar molecules: ethylene carbonate and vinylene carbonate

Results of experimental and theoretical studies of dipole-bound negative ions of the highly polar molecules ethylene carbonate (EC, C3H4O3, mu=5.35 D) and vinylene carbonate (VC, C3H2O3, mu=4.55 D) are presented. These negative ions are prepared in Rydberg electron transfer (RET) reactions in which rubidium (Rb) atoms, excited to ns or nd Rydberg states, collide with EC or VC molecules to produce EC- or VC- ions. In both cases ions are produced only when the Rb atoms are excited to states described by a relatively narrow range of effective principal quantum numbers, n*; the greatest yields of EC- and VC- are obtained for n*(max)=9.0+/-0.5 and 11.6+/-0.5, respectively. Charge transfer from low-lying Rydberg states of Rb is characteristic of a large excess electron binding energy (Eb) of the neutral parent; employing the previously derived empirical relationship Eb=23/n*(max)(2.8) eV, the electron binding energies are estimated to be 49+/-8 meV for EC and 24+/-3 meV for VC. Electron photodetachment studies of EC- show that the excess electron is bound by 49+/-5 meV, in excellent agreement with the RET results, lending credibility to the empirical relationship between Eb and n*(max). Vertical electron affinities for EC and VC are computed employing aug-cc-pVDZ atom-centered basis sets supplemented with a (5s5p) set of diffuse Gaussian primitives to support the dipole-bound electron; at the CCSD(T) level of theory the computed electron affinities are 40.9 and 20.1 meV for EC and VC, respectively.     

The matching of electrostatic extrema: A useful method in drug design? A study of phosphodiesterase III inhibitors

Summary Ligands which bind to a specific protein binding site are often expected to have a similar electrostatic environment which complements that of the binding site. One method of assessing molecular electrostatic similarity is to examine the possible overlay of the maxima and minima in the electrostatic potential outside the molecules and thereby match the regions where strong electrostatic interactions, including hydrogen bonds, with the residues of the binding site may be possible. This approach is validated with accurate calculations of the electrostatic potential, derived from a distributed multipole analysis of an ab initio charge density of the molecule, so that the effects of lone pair and -electron density are correctly included. We have applied this method to the phosphodiesterase (PDE) III substrate adenosine-3,5-cyclic monophosphate (cAMP) and a range of nonspecific and specific PDE III inhibitors. Despite the structural variation between cAMP and the inhibitors, it is possible to match three or four extrema to produce relative orientations in which the inhibitors are sufficiently sterically and electrostatically similar to the natural substrate to account for their affinity for PDE III. This matching of extrema is more apparent using the accurate electrostatic models than it was when this approach was first applied, using semiempirical point charge models. These results reinforce the hypothesis of electrostatic similarity and give weight to the technique of extrema matching as a useful tool in drug design.     

A Product Formed from Glycylglycine in the Presence of Vanadate and Hydrogen Peroxide: The (Glycylde-N-hydroglycinato-kappa(3)N(2),N(N)(),O(1))oxoperoxovanadate(V) Anion

A crystalline glycylglycine complex of monoperoxovanadate has been obtained and its X-ray structure determined. The coordination is pentagonal bipyramidal with the peroxo group and a tridentate glycylglycine occupying the equatorial positions. The axial positions of the anion are occupied by the oxo ligand and by one oxygen of the peroxo group of the adjacent anion. The latter interaction establishes the seventh bond and produces a dimeric structure in the crystalline material. NMR studies of its dissolution in water combined with previously reported results from equilibrium measurements show that the dimer dissociates in water to the monomeric precursor. It is proposed that this monomer corresponds to the complex responsible for the inhibition of the vanadium-catalyzed decomposition of hydrogen peroxide by glycylglycine. Crystal structure of [NEt(4)][VO(O(2))(GlyGly)].1.58H(2)O: monoclinic, space group P2(1); Z = 4; a = 10.618(2) ?; b = 14.803(2) ?; c = 11.809(2) ?; beta = 101.37(2) degrees; V = 1819.7 ?(3); T = 198 K; R(F)() = 0.029 for 2664 data (I(o) >/= 2.5sigma(I(o))) and 431 variables.     

Mononuclear nitrogen/sulfur-ligated zinc methoxide and hydroxide complexes: investigating ligand effects on the hydrolytic stability of zinc alkoxide species

The synthesis and properties of mononuclear zinc methoxide ([(ebnpa)Zn-OCH3]ClO4) (1) and hydroxide ([(ebnpa)Zn-OH]ClO4) (2) complexes of a new mixed nitrogen/sulfur ligand (ebnpa = N-2-(ethylthio)ethyl-N,N-bis(6-neopentylamino-2-pyridylmethyl)amine) are reported. The structures of 1 and 2 were determined by X-ray diffraction. Each possesses a single zinc-coordinated anion (methoxide or hydroxide) and exhibits an overall trigonal bipyramidal geometry. Structural and spectroscopic studies indicate the presence of two hydrogen-bonding interactions involving the oxygen atom of the zinc-bound anion in each complex. Treatment of [(ebnpa)Zn-OH]ClO4 with CH3OH results in the formation of an equilibrium mixture of 1 and 2. 1H NMR spectroscopic methods were used to examine the equilibrium as a function of temperature, yielding KMe (304 K) = 0.30(8), DeltaHMe = -0.9(1) kcal/mol, and DeltaSMe = -5(1) eu. The negative enthalpy indicates that spontaneous zinc alkoxide formation from a hydroxide precursor occurs in this system at low temperature. Using the experimentally determined DeltaHMe value, we found the homolytic Zn-O bond dissociation energy (BDE) in the Zn-OCH3 unit to be approximately -14 kcal/mol relative to the Zn-O BDE in the Zn-OH unit.     

Enzyme electrokinetics: electrochemical studies of the anaerobic interconversions between active and inactive states of Allochromatium vinosum [NiFe]-hydrogenase

The cycling between active and inactive states of the catalytic center of [NiFe]-hydrogenase from Allochromatium vinosum has been investigated by dynamic electrochemical techniques. Adsorbed on a rotating disk pyrolytic graphite "edge" electrode, the enzyme is highly electroactive: this allows precise manipulations of the complex redox chemistry and facilitates quantitative measurements of the interconversions between active catalytic states and the inactive oxidized form Ni(r) (also called Ni-B or "ready") as functions of pH, H(2) partial pressure, temperature, and electrode potential. Cyclic voltammograms for catalytic H(2) oxidation (current is directly related to turnover rate) are highly asymmetric (except at pH > 8 and high temperature) due to inactivation being much slower than activation. Controlled potential-step experiments show that the rate of oxidative inactivation increases at high pH but is independent of potential, whereas the rate of reductive activation increases as the potential becomes more negative. Indeed, at 45 degrees C, activation takes just a few seconds at -288 mV. The cyclic asymmetry arises because interconversion is a two-stage reaction, as expected if the reduced inactive Ni(r)-S state is an intermediate. The rate of inactivation depends on a chemical process (rearrangement and uptake of a ligand) that is independent of potential, but sensitive to pH, while activation is driven by an electron-transfer process, Ni(III) to Ni(II), that responds directly to the driving force. The potentials at which fast activation occurs under different conditions have been analyzed to yield the potential-pH dependence and the corresponding entropies and enthalpies. The reduced (active) enzyme shows a pK of 7.6; thus, when a one-electron process is assumed, reductive activation at pH < 7 involves a net uptake of one proton (or release of one hydroxide), whereas, at pH > 8, there is no net exchange of protons with solvent. Activation is favored by a large positive entropy, consistent with the release of a ligand and/or relaxation of the structure around the active site.     

Determination of sugar structures in solution from residual dipolar coupling constants: methodology and application to methyl beta-D-xylopyranoside

We have developed methodology for the determination of solution structures of small molecules from residual dipolar coupling constants measured in dilute liquid crystals. The power of the new technique is demonstrated by the determination of the structure of methyl beta-d-xylopyranoside (I) in solution. An oriented sample of I was prepared using a mixture of C(12)E(5) and hexanol in D(2)O. Thirty residual dipolar coupling constants, ranging from -6.44 to 4.99 Hz, were measured using intensity-based J-modulated NMR techniques. These include 15 D(HH), 4 (1)D(CH), and 11 (n)D(CH) coupling constants. The accuracy of the dipolar coupling constants is estimated to be < +/- 0.02 Hz. New constant-time HMBC NMR experiments were developed for the measurement of (n)D(CH) coupling constants, the use of which was crucial for the successful structure determination of I, as they allowed us to increase the number of fitted parameters. The structure of I was refined using a model in which the directly bonded interatom distances were fixed at their ab initio values, while 16 geometrical and 5 order parameters were optimized. These included 2 CCC and 6 CCH angles, and 2 CCCC and 6 CCCH dihedral angles. Vibrationally averaged dipolar coupling constants were used during the refinement. The refined solution structure of I is very similar to that obtained by ab initio calculations, with 11 bond and dihedral angles differing by 0.8 degrees or less and the remaining 5 parameters differing by up to 3.3 degrees . Comparison with the neutron diffraction structure showed larger differences attributable to crystal packing effects. Reducing the degree of order by using dilute liquid crystalline media in combination with precise measurement of small residual dipolar coupling constants, as shown here, is a way of overcoming the limitation of strongly orienting liquid crystals associated with the complexity of (1)H NMR spectra for molecules with more than 12 protons.     

Low frequency Raman spectroscopic study on lamellar polymer crystals

The recently observed very low-frequencyRaman bands in polyethylene open up the possibility of determining the straight-chain stem lengths within the crystal as opposed to the full periodicity constituted by crystal core plus amorphous surface layer given by lowangle X-ray diffraction. Past works (6–10) have established the correlation with lamellar thickness but leave doubts about the possibility of making the distinction just indicated. The first objective of the present work was to examine this point on crystals where the chains are inclined to different extents with respect to the lamellar surface. The samples include single-crystals grown from solution at different temperatures and oriented bulk polyethylenes. It is concluded that the low-frequencyRaman peaks correspond to the inclined chain length as opposed to the layer thickness, hence providing a measure of a chain as opposed to a platelet property. TheRaman based length figures are only slightly below those of the appropriately inclined straight traverses across the full lamellae, hence on the basis of existing understanding non-trans chain sequences are confined to a surface region which must be much narrower than required by the usually envisaged amorphous crystalline two-phase structure. Similar studies were carried out also on single crystals of the linear polyester poly(decamethylene sebacate). These demonstrate in the first place that polymers other than polyethylene can also be amenable to suchRaman studies and secondly that in this particular instance, changes in the relative thicknesses of crystal core and surface layer could be registered. The two studies in conjunction affirm the potential of the low-frequencyRaman technique for the study of polymer crystal morphology.