Electrostatic force microscopy and potentiometry of realistic nanostructured systems

We investigate the dependency of electrostatic interaction forces on applied potentials in electrostatic force microscopy (EFM) as well as in related local potentiometry techniques such as Kelvin probe microscopy (KPM). The approximated expression of electrostatic interaction between two conductors, usually employed in EFM and KPM, may loose its validity when probe-sample distance is not very small, as often realized when realistic nanostructured systems with complex topography are investigated. In such conditions, electrostatic interaction does not depend solely on the potential difference between probe and sample, but instead it may depend on the bias applied to each conductor. For instance, electrostatic force can change from repulsive to attractive for certain ranges of applied potentials and probe-sample distances, and this fact cannot be accounted for by approximated models. We propose a general capacitance model, even applicable to more than two conductors, considering values of potentials applied to each of the conductors to determine the resulting forces and force gradients, being able to account for the above phenomenon as well as to describe interactions at larger distances. Results from numerical simulations and experiments on metal stripe electrodes and semiconductor nanowires supporting such scenario in typical regimes of EFM investigations are presented, evidencing the importance of a more rigorous modeling for EFM data interpretation. Furthermore, physical meaning of Kelvin potential as used in KPM applications can also be clarified by means of the reported formalism. © 2009 American Institute of Physics.

The application of electrostatic dry powder deposition technology to coat drug-eluting stents

Purpose: A novel methodology has been introduced to effectively coat intravascular stents with sirolimus-loaded polymeric microparticles. Methods: Dry powders of the microparticulate formulation, consisting of non-erodible polymers, were produced by a supercritical, aerosol, solvent extraction system (ASES). A design of experiment (DOE) approach was conducted on the independent variables, such as organic/CO2 phase volume ratio, polymer weight and stirring-rate, while regression analysis was utilized to interpret the influence of all operational parameters on the dependent variable of particle size. The dry powders, so formed, entered an electric field created by corona charging and were sprayed on the earthed metal stent. Furthermore, the thermal stability of sirolimus was investigated to define the optimum conditions for fusion to the metal surfaces. Results: The electrostatic dry powder deposition technology (EDPDT) was used on the metal strut followed by fusion to produce uniform, reproducible and accurate coatings. The coated stents exhibited sustained release profiles over 25 days, similar to commercial products. EDPDT-coated stents displayed significant reduced platelet adhesion. Conclusions: EDPDT appeared to be a robust accurate and reproducible technology to coat eluting stents.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Reactivity characteristics of cytochrome c, cytochrome oxidase and the electrostatic cytochrome c - cytochrome oxidase complex /

The mechanism whereby cytochrome £ oxidase catalyses elec-. tron transfer from cytochrome £ to oxygen remains an unsolved problem. Polarographic and spectrophotometric activity measurements of purified, particulate and soluble forms of beef heart mitochondrial cytochrome c oxidase presented in this thesis confirm the following characteristics of the steady-state kinetics with respect to cytochrome £: (1) oxidation of ferrocytochrome c is first order under all conditions. -(2) The relationship between sustrate concentration and velocity is of the Michaelis-Menten type over a limited range of substrate. concentrations at high ionic strength. (3) ~he reaction rate is independent from oxygen concentration until very low levels of oxygen. (4) "Biphasic" kinetic plots of enzyme activity as a function of substrate concentration are found when the range of cytochrome c concentrations is extended; the biphasicity ~ is more apparent in low ionic strength buffer. These results imply two binding sites for cytochrome £ on the oxidase; one of high affinity and one of low affinity with Km values of 1.0 pM and 3.0 pM, respectively, under low ionic strength conditions. (5) Inhibition of the enzymic rate by azide is non-c~mpetitive with respect to cytochrome £ under all conditions indicating an internal electron transfer step, and not binding or dissociation of £ from the enzyme is rate limiting. The "tight" binding of cytochrome '£ to cytochrome c oxidase is confirmed in column chromatographic experiments. The complex has a cytochrome £:oxidase ratio of 1.0 and is dissociated in media of high ionic strength. Stopped-flow spectrophotometric studies of the reduction of equimolar mixtures and complexes of cytochrome c and the oxidase were initiated in an attempt to assess the functional relevance of such a complex. Two alternative routes -for reduction of the oxidase, under conditions where the predominant species is the £ - aa3 complex, are postulated; (i) electron transfer via tightly bound cytochrome £, (ii) electron transfer via a small population of free cytochrome c interacting at the "loose" binding site implied from kinetic studies. It is impossible to conclude, based on the results obtained, which path is responsible for the reduction of cytochrome a. The rate of reduction by various reductants of free cytochrome £ in high and low ionic strength and of cytochrome £ electrostatically bound to cytochrome oxidase was investigated. Ascorbate, a negatively charged reagent, reduces free cytochrome £ with a rate constant dependent on ionic strength, whereas neutral reagents TMPD and DAD were relatively unaffected by ionic strength in their reduction of cytochrome c. The zwitterion cysteine behaved similarly to uncharged reductants DAD and TI~PD in exhibiting only a marginal response to ionic strength. Ascorbate reduces bound cytochrome £ only slowly, but DAD and TMPD reduce bound cytochrome £ rapidly. Reduction of cytochrome £ by DAD and TMPD in the £ - aa3 complex was enhanced lO-fold over DAD reduction of free £ and 4-fold over TMPD reduction of free c. Thus, the importance of ionic strength on the reactivity of cytochrome £ was observed with the general conclusion being that on the cytochrome £ molecule areas for anion (ie. phosphate) binding, ascorbate reduction and complexation to the oxidase overlap. The increased reducibility for bound cytochrome £ by reductants DAD and TMPD supports a suggested conformational change of electrostatically bound c compare.d to free .£. In addition, analysis of electron distribution between cytochromes £ and a in the complex suggest that the midpotential of cytochrome ~ changes with the redox state of the oxidase. Such evidence supports models of the oxidase which suggest interactions within the enzyme (or c - enzyme complex) result in altered midpoint potentials of the redox centers.

Statistical potentials for evolutionary studies

Les séquences protéiques naturelles sont le résultat net de l’interaction entre les mécanismes de mutation, de sélection naturelle et de dérive stochastique au cours des temps évolutifs. Les modèles probabilistes d’évolution moléculaire qui tiennent compte de ces différents facteurs ont été substantiellement améliorés au cours des dernières années. En particulier, ont été proposés des modèles incorporant explicitement la structure des protéines et les interdépendances entre sites, ainsi que les outils statistiques pour évaluer la performance de ces modèles. Toutefois, en dépit des avancées significatives dans cette direction, seules des représentations très simplifiées de la structure protéique ont été utilisées jusqu’à présent. Dans ce contexte, le sujet général de cette thèse est la modélisation de la structure tridimensionnelle des protéines, en tenant compte des limitations pratiques imposées par l’utilisation de méthodes phylogénétiques très gourmandes en temps de calcul. Dans un premier temps, une méthode statistique générale est présentée, visant à optimiser les paramètres d’un potentiel statistique (qui est une pseudo-énergie mesurant la compatibilité séquence-structure). La forme fonctionnelle du potentiel est par la suite raffinée, en augmentant le niveau de détails dans la description structurale sans alourdir les coûts computationnels. Plusieurs éléments structuraux sont explorés : interactions entre pairs de résidus, accessibilité au solvant, conformation de la chaîne principale et flexibilité. Les potentiels sont ensuite inclus dans un modèle d’évolution et leur performance est évaluée en termes d’ajustement statistique à des données réelles, et contrastée avec des modèles d’évolution standards. Finalement, le nouveau modèle structurellement contraint ainsi obtenu est utilisé pour mieux comprendre les relations entre niveau d’expression des gènes et sélection et conservation de leur séquence protéique.

Effect of headgroup size, charge, and solvent structure on polymer−micelle interactions, studied by molecular dynamics simulations

We performed atomistic molecular dynamics simulations of anionic and cationic micelles in the presence of poly(ethylene oxide) (PEO) to understand why nonionic water-soluble polymers such as PEO interact strongly with anionic micelles but only weakly with cationic micelles. Our micelles include sodium n-dodecyl sulfate (SDS), n-dodecyl trimethylammonium chloride (DTAC), n-dodecyl ammonium chloride (DAC), and micelles in which we artificially reverse the sign of partial charges in SDS and DTAC. We observe that the polymer interacts hydrophobically with anionic SDS but only weakly with cationic DTAC and DAC, in agreement with experiment. However, the polymer also interacts with the artificial anionic DTAC but fails to interact hydrophobically with the artificial cationic SDS, illustrating that large headgroup size does not explain the weak polymer interaction with cationic micelles. In addition, we observe through simulation that this preference for interaction with anionic micelles still exists in a dipolar "dumbbell" solvent, indicating that water structure and hydrogen bonding alone cannot explain this preferential interaction. Our simulations suggest that direct electrostatic interactions between the micelle and polymer explain the preference for interaction with anionic micelles, even though the polymer overall carries no net charge. This is possible given the asymmetric distribution of negative charges on smaller atoms and positive charges oil larger units in the polymer chain.

Deciphering the role of the electrostatic interactions in the alpha-tropomyosin head-to-tail complex

Skeletal alpha-tropomyosin (Tm) is a dimeric coiled-coil protein that forms linear assemblies under low ionic strength conditions in vitro through head-to-tail interactions. A previously published NMR structure of the Tin head-to-tail complex revealed that it is formed by the insertion of the N-terminal coiled-coil of one molecule into a cleft formed by the separation of the helices at the C-terminus of a second molecule. To evaluate the contribution of charged residues to complex stability, we employed single and double-mutant Tm fragments in which specific charged residues were changed to alanine in head-to-tail binding assays, and the effects of the mutations were analyzed by thermodynamic double-mutant cycles and protein-protein docking. The results show that residues K5, K7, and D280 are essential to the stability of the complex. Though D2, K6, D275, and H276 are exposed to the solvent and do not participate in intermolecular contacts in the NMR structure, they may contribute to head-to-tail complex stability by modulating the stability of the helices at the Tm termini.

Solvent assisted decomposition of the tetrahedral intermediate of the transesterification reaction to biodiesel production. A density functional study

B3LYP/6-31 + G(d) calculations were employed to investigate the mechanism of the transesterification reaction between a model monoglyceride and the methoxide and ethoxide anions. The gas-phase results reveal that both reactions have essentially the same activation energy (5.9 kcal mol(-1)) for decomposition of the key tetrahedral intermediate. Solvent effects were included by means of both microsolvation and the polarizable continuum solvation model CPCM. Both solvent approaches reduce the activation energy, however, only the microsolvation model is able to introduce some differentiation between methanol and ethanol, yielding a lower activation energy for decomposition of the tetrahedral intermediate in the reaction with methanol (1.1 kcal mol(-1)) than for the corresponding reaction with ethanol (2.8 kcal mol(-1)), in line with experimental evidences. Analysis of the individual energy components within the CPCM approach reveals that electrostatic interactions are the main contribution to stabilization of the transition state. (C) 2009 Elsevier Ltd. All rights reserved.

Properties and one-step synthesis of (2-acetylpyridine)tetraammineruthenium(II), [Ru-II(2-acpy)(NH3)(4)](2+) and tetraammine(2-benzoylpyridine)ruthenium(II), [Ru-II(NH3)(4)(2-bzpy)](2+) Redox potentials, UV-Vis and NMR spectra

A more direct and efficient route to the syntheses of [Ru(NH3)(4)(X-Y)](BF4)(2), where X-Y can be 2-acetylpyridine (2-acpy) or 2-benzoylpyridine (2-bzpy), based on the reactions of [RuCl(NH3)(5)]Cl-2 with these ortho-substituted azines is described. The [Ru(2-acpy)(NH3)(4)](BF4)(2) and [Ru(NH3)(5)(2-bzpy)](BF4)(2) complexes have a molar conductance of 328 and 292 Ohm(-1) cm(2) mol(-1), respectively, corresponding to a 1:2 species in solution. These complexes showed two intense absorption bands around 620-650 and 380 nm, the energies of which are solvent dependent, decreasing with the increase of the Gutman's donor number of the solvent, and were assigned as metal-to-ligand charge transfer (MLCT). The complexes have oxidation potentials (Ru-II/III) of +0.380 V vs. Ag/AgCl (2-acpy) and +0.400 V vs. Ag/AgCl (2-bzpy), and reduction potentials (X-Y0/-) of -1.10 V vs. Ag/AgCl (2-acpy) and -0.950 V vs. Ag/AgCl (2-bzpy) on CF3COOH/NaCF3COO at pH=3.0, scan rate 100 mV s(-1), [Ru]=1.0x10(-3) mol l(-1). Both processes show a coupled chemical reaction. Upon oxidation of the metal center, the MLCT absorption bands are bleached and restored upon subsequent reduction. In order to confirm the structure of the complexes a detailed LH NMR investigation was performed in d(6)-acetone. Further confirmation of the structure was obtained by recording the N-15 NMR spectrum of [Ru(NH3)(4)(2-bzpy)](2+) in d(6)-DMSO using the INEPT pulse sequence improving the sensitivity of N-15 by polarization transfer from the protons to the N-15. The Nuclear Overhauser Effect (NOE) experiments were made qualitatively for [Ru(NH3)(4)(2-acpy)](2+), and showed that H-6 of the pyridine is close to a NH3 proton, which should then be in a cis position, and, hence, confirming that acpy is acting as a bidentate ligand. (C) 1999 Elsevier B.V. Ltd. All rights reserved.

Using an effective surface charge to explain surface potentials of Langmuir monolayers from dialkyldimethylammonium halides with the Gouy-Chapman theory

The use of an effective surface charge density has allowed the Gouy-Chapman (CC) theory to explain surface potential isotherms of Langmuir monolayers of dioctadecyldimethylammonium bromide (DODAB). The effective surface charge density of DODAB monolayer increases with the electronegativity of the counterions in the subphase. The pressure-area isotherms indicate a very condensed monolayer for DODAB spread on an I--containing subphase, which exhibits the lowest surface charge density, whereas the monolayer on a F-containing subphase is extremely expanded owing to the high surface charge density or electrostatic repulsion between headgroups. (C) 2001 Published by Elsevier B.V. B.V.

Mapping deformed hyperbolic potentials into nondeformed ones

In this work we introduce a mapping between the so-called deformed hyperbolic potentials, which are presenting a continuous interest in the last few years, and the corresponding nondeformed ones. As a consequence, we conclude that these deformed potentials do not pertain to a new class of exactly solvable potentials, but to the same one of the corresponding nondeformed ones. Notwithstanding, we can reinterpret this type of deformation as a kind of symmetry of the nondeformed potentials. © 2005 Elsevier B.V. All rights reserved.

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions?

The collapse of linear polyelectrolyte chains in a poor solvent: When does a collapsing polyelectrolyte collect its counter ions? The collapse of polyions in a poor solvent is a complex system and is an active research subject in the theoretical polyelectrolyte community. The complexity is due to the subtle interplay between hydrophobic effects, electrostatic interactions, entropy elasticity, intrinsic excluded volume as well as specific counter-ion and co-ion properties. Long range Coulomb forces can obscure single molecule properties. The here presented approach is to use just a small amount of screening salt in combination with a very high sample dilution in order to screen intermolecular interaction whereas keeping intramolecular interaction as much as possible (polyelectrolyte concentration cp ≤ 12 mg/L, salt concentration; Cs = 10^-5 mol/L). This is so far not described in literature. During collapse, the polyion is subject to a drastic change in size along with strong reduction of free counterions in solution. Therefore light scattering was utilized to obtain the size of the polyion whereas a conductivity setup was developed to monitor the proceeding of counterion collection by the polyion. Partially quaternized PVP’s below and above the Manning limit were investigated and compared to the collapse of their uncharged precursor. The collapses were induced by an isorefractive solvent/non-solvent mixture consisting of 1-propanol and 2-pentanone, with nearly constant dielectric constant. The solvent quality for the uncharged polyion could be quantified which, for the first time, allowed the experimental investigation of the effect of electrostatic interaction prior and during polyion collapse. Given that the Manning parameter M for QPVP4.3 is as low as lB / c = 0.6 (lB the Bjerrum length and c the mean contour distance between two charges), no counterion binding should occur. However the Walden product reduces with first addition of non solvent and accelerates when the structural collapse sets in. Since the dielectric constant of the solvent remains virtually constant during the chain collapse, the counterion binding is entirely caused by the reduction in the polyion chain dimension. The collapse is shifted to lower wns with higher degrees of quaternization as the samples QPVP20 and QPVP35 show (M = 2.8 respectively 4.9). The combination of light scattering and conductivity measurement revealed for the first time that polyion chains already collect their counter ions well above the theta-dimension when the dimensions start to shrink. Due to only small amounts of screening salt, strong electrostatic interactions bias dynamic as well as static light scattering measurements. An extended Zimm formula was derived to account for this interaction and to obtain the real chain dimensions. The effective degree of dissociation g could be obtained semi quantitatively using this extrapolated static in combination with conductivity measurements. One can conclude the expansion factor a and the effective degree of ionization of the polyion to be mutually dependent. In the good solvent regime g of QPVP4.3, QPVP20 and QPVP35 exhibited a decreasing value in the order 1 > g4.3 > g20 > g35. The low values of g for QPVP20 and QPVP35 are assumed to be responsible for the prior collapse of the higher quaternized samples. Collapse theory predicts dipole-dipole attraction to increase accordingly and even predicts a collapse in the good solvent regime. This could be exactly observed for the QPVP35 sample. The experimental results were compared to a newly developed theory of uniform spherical collapse induced by concomitant counterion binding developed by M. Muthukumar and A. Kundagrami. The theory agrees qualitatively with the location of the phase boundary as well as the trend of an increasing expansion with an increase of the degree of quaternization. However experimental determined g for the samples QPVP4.3, QPVP20 and QPVP35 decreases linearly with the degree of quaternization whereas this theory predicts an almost constant value.

The theory of electrostatic probes in strong magnetic fields : Thesis submitted to the Faculty of The Graduate School of the University of Colorado for the Degree Doctor of Philosophy Department of Aerospace Engineering Science

A theory is developed of an electrostatic probe in a fully-ionized plasma in the presence of a strong magnetic field. The ratio of electron Larmor radius to probe transverse dimension is assumed to be small. Poisson's equation, together with kinetic equations for ions and electrons are considered. An asymptotic perturbation method of multiple scales is used by considering the characteristic lengths appearing in the problem. The leading behavior of the solution is found. The results obtained appear to apply to weaker fields also, agreeing with the solutions known in the limit of no magnetic field. The range of potentials for wich results are presented is limited. The basic effects produced by the field are a depletion of the plasma near the probe and a non-monotonic potential surrounding the probe. The ion saturation current is not changed but changes appear in both the floating potential Vf and the slope of the current-voltage diagram at Vf. The transition region extends beyond the space potential Vs,at wich point the current is largely reduced. The diagram does not have an exponential form in this region as commonly assumed. There exists saturation in electron collection. The extent to which the plasma is disturbed is determined. A cylindrical probe has no solution because of a logarithmic singularity at infinity. Extensions of the theory are considered.

Electrostatic steering and ionic tethering in enzyme–ligand binding: Insights from simulations

To bind at an enzyme’s active site, a ligand must diffuse or be transported to the enzyme’s surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and β-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as “ionic tethering.” We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme’s surroundings even when the substrate is nonpolar.