Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

In addition to the biologically active monomer of the protein insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The insulin monomer contains two distinct surfaces, namely, the dimer forming surface (DFS) and the hexamer forming surface (HFS), that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behavior of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. In particular, we find that more structured water molecules, with higher residence times (similar to 300-500 ps), are present near HFS than those near DFS. A significant slowing down is observed in the decay of associated rotational auto time correlation functions of O-H bond vector of water in the vicinity of HFS. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly an outward H-atoms directed, clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inward H-atoms directed orientation owing to the flat curvature of hydrophobic surface and the interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces that reveal the significant differences between the two hydration layers.

Charge density distribution and electrostatic interactions of ethionamide: an inhibitor of the enoyl acyl carrier protein reductase (inhA) enzyme of Mycobacterium tuberculosis

An experimental charge density analysis of an anti-TB drug ethionamide was carried out from high resolution X-ray diffraction at 100 K to understand its charge density distribution and electrostatic properties. The experimental results were validated from periodic theoretical charge density calculations performed using CRYSTAL09 at the B3LYP/6-31G** level of theory. The electron density rho(bcp)(r) and the Laplacian of electron density del(2)(rho bcp)(r) of the molecule calculated from both the methods display the charge density distribution of the ethionamide molecule in the crystal field. The electrostatic potential map shows a large electropositive region around the pyridine ring and a large electronegative region at the vicinity of the thiol atom. The calculated experimental dipole moment is 10.6D, which is higher than the value calculated from theory (8.2D). The topological properties of C-H center dot center dot center dot S, N-H center dot center dot center dot N and N-H center dot center dot center dot S hydrogen bonds were calculated, revealing their strength. The charge density analysis of the ethionamide molecule determined from both the experiment and theory gives the topological and electrostatic properties of the molecule, which allows to precisely understand the nature of intra and intermolecular interactions.

A simple computational model for predicting .pi.-facial selectivity in reductions of sterically unbiased ketones. Relative importance of electrostatic and orbital interactions

Various factore controlling the preferred facial selectivity in the reductions of a number of sterically unbiased ketones have been evaluated using a semiempirical MO procedure. MNDO optimized geometries do not reveal any significant ground-state distortions which can be correlated with the observed face selectivities. Electrostatic effecta due to an approaching reagent were modeled by placing a test negative charge at a fixed distance from the carbonyl carbon on each of the two faces. A second series of calculations was carried out using the hydride ion as a test nucleophile. The latter calculations effectively include orbital interactions involving the u and u* orbitals of the newly formed bond in the reaction. The computed energy differences with the charge model are generally much larger compared to those with the hydride ion. However, both models lead to predictions which are qualitatively consistent with the experimentally determined facial preferences for most of the systems. Thus, electrostatic interactions between the nucleophile and the substrate seem to effectively determine the face selectivities in these molecules. However, there are a few exceptions in which orbital interactions are found to contribute significantly and occasionally reverse the preference dictated by electrostatic effecta. The remarkable succew of the hydride model calculations, in spite of retaining the unperturbed geometries of the substrates, points to the unimportance of torsional effeds and orbital distortions associated with the pyramidalized carbonyl unit in the transition state in most of the substrates considered. Additional experimental results are reported which provide useful calibration for the present computational approach.

Interactions of linear dicationic molecules with lipid A: structural requisites for optimal binding affinity

The structural determinants of the binding affinity of linear dicationic molecules toward lipid A have been examined with respect to the distance between the terminal cationic functions, the basicity, and the type of cationic moieties using a series of spermidine derivatives and pentamidine analogs by fluorescence spectroscopic methods, The presence of two terminal cationic groups corresponds to enhanced affinity, A distinct sigmoidal relationship between the intercationic distance and affinity was observed with a sharp increase at 11 Angstrom, levelling off at about 13 Angstrom. The basicity (pK) and nature of the cationic functions are poor correlates of binding potency, since molecules bearing primary amino, imidazolino, or guanido termini are equipotent, The interaction of pentamidine, a bisamidine drug, with lipid A, characterized in considerable detail employing the putative intermolecular excimerization of the drug, suggests a stoichiometry of 1:1 in the resultant complex, The binding is driven almost exclusively by electrostatic forces, and is dependent on the ionization states of both lipid A and the drug, Under conditions when lipid A is highly disaggregated, pentamidine binds specifically to bis-phosphoryl- but not to monophosphoryl-lipid A indicating that both phosphate groups of lipid A are necessary for electrostatic interactions by the terminal amidininium groups of the drug, Based on these data, a structural model is proposed for the pentamidine-lipid A complex, which may be of value in designing endotoxin antagonists from first principles.

Evaluation of solute-solvent interactions from solvent blue-shifts of n → π* transitions of C=O, C=S, NO2 and N=N groups: hydrogen bond energies of various donor-acceptor systems

Effects of non-polar, polar and proton-donating solvents on the n → π* transitions of C=O, C=S, NO2 and N=N groups have been investigated. The shifts of the absorption maxima in non-polar and polar solvents have been related to the electrostatic interactions between solute and solvent molecules, by employing the theory of McRAE. In solvents which can donate protons the solvent shifts are mainly determined by solute-solvent hydrogen bonding. Isobestic points have been found in the n → π* bonds of ethylenetrithio-carbonate in heptane-alcohol and heptane-chloroform solvent systems, indicating the existence of equilibria between the hydrogen bonded and the free species of the solute. Among the different proton-donating solvents studied water produces the largest blue-shifts. The blue-shifts in alcohols decrease in the order 2,2,2-trifluoroethanol, methanol, ethanol, isopropanol and t-butanol, the blue-shift in trifluoroethanol being nearly equal to that in water. This trend is exactly opposite to that for the self-association of alcohols. It is suggested that electron-withdrawing groups not merely decrease the extent of self-association of alcohols, but also increase the ability to donate hydrogen bonds. The approximate hydrogen-bond energies for several donor-acceptor systems have been estimated. In a series of aliphatio ketones and nitro compounds studied, the blue-shifts and consequently the hydrogen bond energies decrease with the decrease in the electron-withdrawing power of the alkyl groups. It is felt that electron-withdrawing groups render the chromophores better proton acceptors, and the alcohols better donors. A linear relationship between n → π* transition frequency and the infrared frequency of ethylenetrithiocarbonate has been found. It is concluded that stabilization of the electronic ground states of solute molecules by electrostatic and/or hydrogen-bond interactions determines the solvent shifts.

Ions in water: Role of attractive interactions in size dependent diffusivity maximum

A molecular dynamics study of model ions in water is reported. The van der Waals diameter of both the cations and anions is varied. We have carried out two sets of simulations-with and without dispersion interaction-between the ion and water. Self-diffusivity of the ions exhibits an anomalous maximum as a function of the van der Waals diameter for both these sets. This existence of a maximum in self-diffusivity when there is no dispersion interaction between the ion and the water is attributed to the attractive term from electrostatic interactions. Detailed analysis of this effect shows that the solvent shell is more strongly defined in the presence of dispersion interactions. A smaller ion exhibits biexponential decay while a single exponential decay is seen for the ion with maximum diffusivity in the self-part of the intermediate scattering function. The solvent structure around the ion appears to determine much of the dynamics of the ion. Interesting trends are seen in the activation energies and these can be understood in terms of the levitation effect. (C) 2010 American Institute of Physics. doi:10.1063/1.3481656]

Interactions of tris buffer with nucleotides: the crystal structure of tris(hydroxymethyl)methylammonium adenosine 5'-diphosphate dihydrate

The crystal and molecular structures of the Tris salt of adenosine 5'-diphosphate were determined from X-ray diffraction data. The crystals are monoclinic, space P21, and Z = 2 with a=9.198 (2) A, b=6.894 (1) A, c=18.440 (4) A, and beta = 92.55 (2) degrees. Intensity data were collected on an automated diffractometer. The structure was solved by the heavy-atom technique and refined by least squares to R = 0.047. The ADP molecule adopts a folded conformation. The conformation about the glycosidic bond is anti. The conformation of the ribose ring is close to a perfect C(2')-endo-C-(3')-exo puckering. The conformation about C(4')-C(5') is gauche-gauche, similar to other nucleotide structures. The pyrophosphate chain displays a nearly eclipsed geometry when viewed down the P-P vector, unlike the staggered conformation observed in crystal structures of other pyrophosphates. The less favorable eclipsed conformation probably results from the observed association of Tris molecules with the polar diphosphate chain through electrostatic interactions and hydrogen bonds. Such interactions may play an important role in Tris-buffered aqueous solutions of nucleotides and metal ions.

Unusual Salt-Induced Color Modulation through Aggregation-Induced Emission Switching of a Bis-cationic Phenylenedivinylene-Based pi Hydrogelator

The synthesis, hydrogelation, and aggregation-induced emission switching of the phenylenedivinylene bis-N-octyl pyridinium salt is described. Hydrogelation occurs as a consequence of pi-stacking, van der Waals, and electrostatic interactions that lead to a high gel melting temperature and significant mechanical properties at a very low weight percentage of the gelator. A morphology transition from fiber-to-coil-to-tube was observed depending on the concentration of the gelator. Variation in the added salt type, salt concentrations, or temperature profoundly influenced the order of aggregation of the gelator molecules in aqueous solution. Formation of a novel chromophore assembly in this way leads to an aggregation-induced switch of the emission colors. The emission color switches from sky blue to white to orange depending upon the extent of aggregation through mere addition of external inorganic salts. Remarkably, the salt effect on the assembly of such cationic phenylenedivinylenes in water follow the behavior predicted from the well-known Hofmeister effects. Mechanistic insights for these aggregation processes were obtained through the counterion exchange studies. The aggregation-induced emission switching that leads to a room-temperature white-light emission from a single chromophore in a single solvent (water) is highly promising for optoelectronic applications.

Copper(II) Complexes of l-Arginine as Netropsin Mimics Showing DNA Cleavage Activity in Red Light

Copper(II) complexes [Cu(L-arg)(2)](NO3)(2) (1) and [Cu(L-arg)(B)Cl]Cl (2-5), where B is a heterocyclic base, namely, 2,2'-bipyridine (bpy, 2), 1,10-phenanthroline (phen, 3), dipyrido[3,2-d:2',3'-f]quinoxaline (dpq, 4), and dipyrido[3,2-a:2',3'-c)phenazine (dppz, 5), are prepared and their DNA binding and photoinduced DNA cleavage activity studied. Ternary complex 3, structurally characterized using X-ray crystallography, shows a square-pyramidal (4 + 1) coordination geometry in which the N,O-donor L-arginine and N,N-donor 1,10-phenanthroline form the basal plane with one chloride at the elongated axial site. The complex has a pendant cationic guanidinium moiety. The one-electron paramagnetic complexes display a metal-centered d-d band in the range of 590-690 nm in aqueous DMF They show quasireversible cyclic voltammetric response due to the Cu(II)/Cu(I) couple in the range of -0.1 to -0.3 V versus a saturated calomel electrode in a DMF-Tris HCl buffer (pH 7.2). The DNA binding propensity of the complexes is studied using various techniques. Copper(II) bis-arginate 1 mimics the minor groove binder netropsin by showing preferential binding to the AT-rich sequence of double-strand (ds) DNA. DNA binding study using calf thymus DNA gives an order: 5 (L-arg-dppz) >= 1 (biS-L-arg) > 4 (L-arg-dpq) > 3 (L-arg-phen) >> 2 (L-arg-bpy). Molecular docking calculations reveal that the complexes bind through extensive hydrogen bonding and electrostatic interactions with ds-DNA. The complexes cleave supercoiled pUC19 DNA in the presence of 3-mercaptopropionic acid as a reducing agent forming hydroxyl ((OH)-O-center dot) radicals. The complexes show oxidative photoinduced DNA cleavage activity in UV-A light of 365 nm and red light of 647.1 nm (Ar-Kr mixed-gas-ion laser) in a metal-assisted photoexcitation process forming singlet oxygen (O-1(2)) species in a type-II pathway. All of the complexes, barring complex 2, show efficient DNA photocleavage activity. Complexes 4 and 5 exhibit significant double-strand breaks of DNA in red light of 647.1 nm due to the presence of two photosensitizers, namely, L-arginine and dpq or dppz in the molecules.

Structure and Rheology of Cationic Molecular Hydrogels of Quinuclidine Grafted Bile Salts. Influence of the Ionic Strength and Counter-Ion type

Quinuclidine grafted cationic bile salts are forming salted hydrogels. An extensive investigation of the effect of the electrolyte and counterions on the gelation has been envisaged. The special interest of the quinuclidine grafted bile salt is due to its broader experimental range of gelation to study the effect of electrolyte. Rheological features of the hydrogels are typical of enthalpic networks exhibiting a scaling law of the elastic shear modulus with the concentration (scaling exponent 2.2) modeling cellular solids in which the bending modulus is the dominant parameter. The addition of monovalent salt (NaCl) favors the formation of gels in a first range (0.00117 g cm-3 (0.02 M) < TNaCl < 0.04675 g cm-3 (0.8 M)). At larger salt concentrations, the gels become more heterogeneous with nodal zones in the micron scale. Small-angle neutron scattering experiments have been used to characterize the rigid fibers ( ≈ 68 Å) and the nodal zones. Stress sweep and creeprecovery measurements are used to relate the lack of linear viscoelastic domain to a mechanism of disentanglement of the fibers from their associations into fagots. The electrostatic interactions can be screened by addition of salt to induce a progressive evolution toward flocculation. SEM, UV absorbance, and SAXS study of the Bragg peak at large Q-values complete the investigation.

Electrochemical Reduction of Oriented Graphene Oxide Films: An in Situ Raman Spectroelectrochemical Study

Graphene oxide (GO) is assembled on a gold substrate by a layer-by-layer technique using a self-assembled cystamine monolayer. The negatively charged GO platelets are attached to the positively charged cystamine monolayer through electrostatic interactions. Subsequently, it is shown that the GO can be reduced electrochemically using applied DC bias by scanning the potential from 0 to -1 V vs a saturated calomel electrode in an aqueous electrolyte. The GO and reduced graphene oxide (RGO) are characterized by Raman spectroscopy and atomic force microscopy (AFM). A clear shift of the G band from 1610 cm-1 of GO to 1585 cm-1 of RGO is observed. The electrochemical reduction is followed in situ by micro Raman spectroscopy by carrying out Raman spectroscopic studies during the application of DC bias. The GO and RGO films have been characterized by conductive AFM that shows an increase in the current flow by at least 3 orders of magnitude after reduction. The electrochemical method of reducing GO may open up another way of controlling the reduction of GO and the extent of reduction to obtain highly conducting graphene on electrode materials.

Synthetic peptides as chemoattractants for bull spermatozoa structure activity correlations

The ability of various synthetic peptide analogs of. Formyl-Met-Leu-Phe to induce chemotaxis in bull sperm is compared using an inverted capillary assay. The formyl group is essential for chemotactic activity and corresponding t-butyloxycarbonyl tripeptides are inactive. Sequence analogs, Formyl-Met-Phe-Leu, Formyl-Leu-Met-Phe and Formyl-Leu-Phe-Met are active. Replacement of Met and Leu by Pro does not diminish activity. Formyl-Met-Leu-Phe-NH2 is active suggesting that electrostatic interactions involving the carboxyl group may be unimportant in receptor interactions. The studies establish the importance of an amino terminal formyl group and a sequence of at least three hydrophobic residues, for inducing sperm chemotaxis.

Elastic Properties of Boron Nitride Nanotubes and Their Comparison with Carbon Nanotubes

Boron Nitride Nanotubes (BNNTs) have alternating boron and nitrogen atoms in graphite like network and are strongly polar in nature due to a large charge on boron and nitrogen atoms. Hence electrostatic interactions are expected to play an important role in determining the elastic properties of BNNTs. In the absence of specific partial atomic charge information for boron and nitrogen, we have studied the elastic properties BNNTs varying the partial atomic charges on boron and nitrogen. We have computed Young modulus (Y) and Shear modulus (G) of BNNT as a function of the tube radius and number of walls using molecular mechanics calculation. Our calculation shows that Young modulus of BNNTs increases with increase in magnitude of the partial atomic charge on B and N and can be larger than the Young modulus of CNTs of same radius. This is in contrast to the earlier finding that CNTs has the largest tensile strength (PRL, 80, 4502, 1998). Shear modulus, on the other hand depends weakly on the magnitude of partial atomic charge and is less than the shear modulus of the CNT. The values obtained for Young modulus and Shear modulus are in excellent agreement with the available experimental results.

Reversible electrochemistry of cytochrome c on recompressed, binderless exfoliated graphite electrodes

The direct electrochemistry of cytochrome c (cyt-c) has been investigated on exfoliated graphite (EG) electrodes. The as-polished and roughened (using SiC emery sheet) EG surfaces are inactive for the direct electron transfer. However, when the EG electrode was sonicated before the experiment, a pair of redox waves were obtained for freely diffusing cyt-c in the solution phase. The formal potential was found to be 0.01 V (vs. SCE) in 0.1 M phosphate buffer at a pH of 7.1. The electrochemical response for the adsorbed cyt-c on sonicated EG electrodes, which is shown to have carbonyl functional groups on its surface, shows nearly reversible voltammograms in the same electrolyte. However, the formal potential in the adsorbed state is more negative than that observed for the solution phase cyt-c. A structure based on an open heme conformation proposed by Hildebrandt and Stockburger is probably present on the EG surface. It is suggested that the electrochemistry at the EG electrode is essentially governed by favourable electrostatic interactions.

Effect of coordinated ions on structure and flexibility of parallel G-quadruplexes: A molecular dynamics study

ingle tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. There quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Angstrom from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Nai counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. in the absence of any coordinated ion. due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.