Advantage of the ether linkage between the positive charge and the cholesteryl skeleton in cholesterol-based amphiphiles as vectors for gene delivery

Twelve novel cationic cholesterol derivatives with different linkage types between the cationic headgroup and the cholesteryl backbone have been developed. These have been tested for their efficacies as gene transfer agents as mixtures with dioleoyl phosphatidylethanolamine (DOPE). A pronounced improvement in transfection efficiency was observed when the cationic center was linked to the steroid backbone using an ether type bond. Among these, cholest-5-en-3b-oxyethane-N, N,N-trimethylammonium bromide (2a) and cholest-5-en-3b-oxyethane-N, N-dimethyl-N-2-hydroxyethylammonium bromide (3d) showed transfection efficiencies considerably greater than commercially available reagents such as Lipofectin or Lipofectamine. To achieve transfection, 3d did not require DOPE. Increasing hydration at the headgroup level for both ester- and ether-linked amphiphiles resulted in progressive loss of transfection efficiency. Transfection efficiency was also greatly reduced when a 'disorder'-inducing chain like an oleyl (cis-9-octadecenyl) segment was added to these cholesteryl amphiphiles. Importantly, the transfection ability of 2a with DOPE in the presence of serum was significantly greater than for a commercially available reagent, Lipofectamine. This suggests that these novel cholesterol-based amphiphiles might prove promising in applications involving liposome-mediated gene transfection. This investigation demonstrates the importance of structural features at the molecular level for the design of cholesterol-based gene delivery reagents that would aid the development of newer, more efficient formulations based on this class of molecules.

Vibrational phase relaxation of O–H stretch in bulk water: Role of large amplitude angular jumps and negative cross-correlations among the forces on the O–H bond

The ultrafast vibrational phase relaxation of O–H stretch in bulk water is investigated in molecular dynamics simulations. The dephasing time (T2) of the O–H stretch in bulk water calculated from the frequency fluctuation time correlation function (Cω(t)) is in the range of 70–80 femtosecond (fs), which is comparable to the characteristic timescale obtained from the vibrational echo peak shift measurements using infrared photon echo [W.P. de Boeij, M.S. Pshenichnikov, D.A. Wiersma, Ann. Rev. Phys. Chem. 49 (1998) 99]. The ultrafast decay of Cω(t) is found to be responsible for the ultrashort T2 in bulk water. Careful analysis reveals the following two interesting reasons for the ultrafast decay of Cω(t). (A) The large amplitude angular jumps of water molecules (within 30–40 fs time duration) provide a large scale contribution to the mean square vibrational frequency fluctuation and gives rise to the rapid spectral diffusion on 100 fs time scale. (B) The projected force, due to all the atoms of the solvent molecules on the oxygen (FO(t)) and hydrogen (FH(t)) atom of the O–H bond exhibit a large negative cross-correlation (NCC). We further find that this NCC is partly responsible for a weak, non-Arrhenius temperature dependence of the dephasing rate.

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.

Hydrogen-bond dynamics near a micellar surface: Origin of the universal slow relaxation at complex aqueous interfaces

The dynamics of hydrogen bonds among water molecules themselves and with the polar head groups (PHG) at a micellar surface have been investigated by long molecular dynamics simulations. The lifetime of the hydrogen bond between a PHG and a water molecule is found to be much longer than that between any two water molecules, and is likely to be a general feature of hydrophilic surfaces of organized assemblies. Analyses of individual water trajectories suggest that water molecules can remain bound to the micellar surface for more than 100 ps. The activation energy for such a transition from the bound to a free state for the water molecules is estimated to be about 3.5 kcal/mol.

Increased glycosidic bond stabilities in 4-C-hydroxymethyl linked disaccharides

Three new hydroxymethyl-linked non-natural disaccharide analogues, containing an additional methylene group in between the glycosidic linkage, were synthesized by utilizing 4-C-hydroxymethyl-alpha-D-glucopyranoside as the glycosyl donor. A kinetic study was undertaken to assess the hydrolytic stabilities of these new disaccharide analogues toward acid-catalyzed hydrolysis, at 60 degrees C and 70 degrees C. The studies showed that the disaccharide analogues were stable, by an order of magnitude, than naturally-occurring disaccharides, such as, cellobiose, lactose, and maltose. The first order rate constants were lower than that of methyl glycosides and the trend of hydrolysis rate constants followed that of naturally-occurring disaccharides. alpha-Anomer showed faster hydrolysis than the beta-anomer and the presence of axial hydroxyl group also led to faster hydrolysis among the disaccharide analogues. Energy minimized structures, derived through molecular modeling, showed that dihedral angles around the glycosidic bond in disaccharide analogues were nearly similar to that of naturally-occurring disaccharides. (C) 2011 Elsevier Ltd. All rights reserved.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Natural Convection in a Vertical Channel with Arrays of Flush-Mounted Heaters on Opposite Conductive Walls

Three-dimensional numerical study of natural convection in a vertical channel with flush-mounted discrete heaters on opposite conductive substrate walls is carried out in the present work. Detailed flow and heat transfer characteristics are presented for various Grashof numbers. The heat transfer effects on one wall by the presence of heaters on its opposite wall is examined. It is found that heat transfer rates on one wall are increased by the presence of heaters on its opposite wall. The thermal boundary layers on the opposite walls complement each other for enhanced heat transfer. The effects of spacing between the heated walls, spacings between heaters and substrate conductivity on flow and heat transfer are examined. Existence of optimum spacings between the heated walls for maximum heat transfer and mass flow are observed. It is found that the heat transfer and fluid flow do not follow the same optimum spacings. Mass flow rate reaches maximum value at a wall spacing greater than the spacing for maximum heat transfer. This is because the interaction of thermal boundary layers on individual walls ceases at a lower spacing before the velocity boundary layers separate each other. It is found that increased spacings between heaters reduce individual heater temperatures provided the heaters close to exit on both substrates avail sufficient substrate potions on the exit side. Insufficient substrate portions between the exit heaters and the exit cause abnormal local temperature rise in the exit heaters which are the hottest ones among all the heaters. Optimal heater spacings exist for minimum hottest heater temperature rise. Correlations are presented for dimensionless mass flow rate, temperature maximum, and average Nusselt number.

Tautomeric Preference and Conformation Locking in Fenobam, Thiofenobam, and Their Analogues: The Decisive Role of Hydrogen Bond Hierarchy

The crystal and molecular structures of the potential antidepressant drug fenobam and its derivatives are examined in terms of the preferred form among the two possible tautomeric structures. In this study, chemical derivatization has been utilized as a means to ``experimentally simulate'' the tautomeric preference and conformational variability in fenobam. Eight new derivatives of fenobam have been synthesized, and structural features have been characterized by single-crystal X-ray diffraction and NMR spectroscopy. The specific tautomeric preference found in all of these compounds and their known crystal forms have been construed in terms of the stabilizing intramolecular N-H center dot center dot center dot O and N-H center dot center dot center dot S hydrogen bonding. The hierarchy of intramolecular hydrogen bonds evidenced as the preference of the C-H center dot center dot center dot O hydrogen bond over C-H center dot center dot center dot N and that of the C-H center dot center dot center dot N hydrogen bond over C-H center dot center dot center dot S explains the two distinct conformations adopted by fenobam and thiofenobam derivatives. The relative energy values of different molecular conformations have been calculated and compared.

Use of polydispersity index as control parameter to study melting/freezing of Lennard-Jones system: Comparison among predictions of bifurcation theory with Lindemann criterion, inherent structure analysis and Hansen-Verlet rule

Using polydispersity index as an additional order parameter we investigate freezing/melting transition of Lennard-Jones polydisperse systems (with Gaussian polydispersity in size), especially to gain insight into the origin of the terminal polydispersity. The average inherent structure (IS) energy and root mean square displacement (RMSD) of the solid before melting both exhibit quite similar polydispersity dependence including a discontinuity at solid-liquid transition point. Lindemann ratio, obtained from RMSD, is found to be dependent on temperature. At a given number density, there exists a value of polydispersity index (delta (P)) above which no crystalline solid is stable. This transition value of polydispersity(termed as transition polydispersity, delta (P) ) is found to depend strongly on temperature, a feature missed in hard sphere model systems. Additionally, for a particular temperature when number density is increased, delta (P) shifts to higher values. This temperature and number density dependent value of delta (P) saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (delta (TP)). This value (delta (TP) similar to 0.11) is in excellent agreement with the experimental value of 0.12, but differs from hard sphere transition where this limiting value is only 0.048. Terminal polydispersity (delta (TP)) thus has a quasiuniversal character. Interestingly, the bifurcation diagram obtained from non-linear integral equation theories of freezing seems to provide an explanation of the existence of unique terminal polydispersity in polydisperse systems. Global bond orientational order parameter is calculated to obtain further insights into mechanism for melting.

The role of lysine-41 in RNase A catalysis - A quantum chemical study on the active site ligand complex

Several mechanisms have been proposed to explain the action of enzymes at the atomic level. Among them, the recent proposals involving short hydrogen bonds as a step in catalysis by Gerlt and Gassman [1] and proton transfer through low barrier hydrogen bonds (LBHBs) [2, 3] have attracted attention. There are several limitations to experimentally testing such hypotheses, Recent developments in computational methods facilitate the study of active site-ligand complexes to high levels of accuracy, Our previous studies, which involved the docking of the dinucleotide substrate UpA to the active site of RNase A [4, 5], enabled us to obtain a realistic model of the ligand-bound active site of RNase A. From these studies, based on empirical potential functions, we were able to obtain the molecular dynamics averaged coordinates of RNase A, bound to the ligand UpA. A quantum mechanical study is required to investigate the catalytic process which involves the cleavage and formation of covalent bonds. In the present study, we have investigated the strengths of some of the hydrogen bonds between the active site residues of RNase A and UpA at the ab initio quantum chemical level using the molecular dynamics averaged coordinates as the starting point. The 49 atom system and other model systems were optimized at the 3-21G level and the energies of the optimized systems were obtained at the 6-31G* level. The results clearly indicate the strengthening of hydrogen bonds between neutral residues due to the presence of charged species at appropriate positions. Such a strengthening manifests itself in the form of short hydrogen bonds and a low barrier for proton transfer. In the present study, the proton transfer between the 2'-OH of ribose (from the substrate) and the imidazole group from the H12 of RNase A is influenced by K41, which plays a crucial role in strengthening the neutral hydrogen bond, reducing the barrier for proton transfer.

Facile C-O bond scission in alcohols on Zn surfaces

Interaction of methanol, ethanol, and 2-propanol with polycrystalline as well as (0001) surfaces of Zn has been investigated by photoelectron spectroscopy and vibrational energy loss spectroscopy. All the alcohols show evidence for the condensed species along with the chemisorbed species at 80 K. With increase in temperature to similar to 120 K, the condensed species desorbs, leaving the chemisorbed species which decomposes to give the alkoxy species. The alkoxy species is produced increasingly at lower temperatures as we go from methanol to 2-propanol, the 2-propoxy species occurring even at 80 K. The alkoxy species undergo C-O bond scission giving rise to a hydrocarbon species and oxygen. The C-O bond cleavage occurs at a relatively low temperature of similar to 150 K. The effect of preadsorbed oxygen is to stabilize the methoxy species and prevent C-O bond scission. On the other hand, coadsorption of oxygen with methanol favors the formation of the methoxy species and gives rise to hydrocarbon species arising from the C-O bond scission even at 80 K.

Durability of rammed earth walls exposed for 20 years to natural weathering

This paper presents a study on the durability of different types of stabilised and unstabilised rammed earth walls. These rammed earth walls were constructed and exposed for 20 years to natural weathering, in a wet continental climate. None of these walls have shown complete collapse to date. A method to measure the rammed earth walls erosion by stereo-photogrammetry has been developed. The result shows that the mean erosion depth of the studied walls is about 2 mm (0.5% wall thickness) in the case of rammed earth wall stabilised with 5% by dry weight of hydraulic lime and about 6.4 mm (1.6% wall thickness) in the case of unstabilised rammed earth walls. The stabilisation enables to not use any plaster to protect the walls. In the case of the unstabilised rammed earth walls, an extrapolated lifetime longer than 60 years can be assessed. This shows a potential for the use of unstabilised rammed earth in the similar climatic conditions with this study. The method of stereo-photogrammetry used to measure the erosion of rammed earth walls on site may also help to calibrate and develop more pertinent laboratory test to assess the durability of rammed earth wall.

Modeling of a pressure modulated desalination system using bond graph methodology

A desalination system is a complex multi energy domain system comprising power/energy flow across several domains such as electrical, thermal, and hydraulic. The dynamic modeling of a desalination system that comprehensively addresses all these multi energy domains is not adequately addressed in the literature. This paper proposes to address the issue of modeling the various energy domains for the case of a single stage flash evaporation desalination system. This paper presents a detailed bond graph modeling of a desalination unit with seamless integration of the power flow across electrical, thermal, and hydraulic domains. The paper further proposes a performance index function that leads to the tracking of the optimal chamber pressure giving the optimal flow rate for a given unit of energy expended. The model has been validated in steady state conditions by simulation and experimentation.

Bond graph toolbox for handling complex variable

Bond graph is an apt modelling tool for any system working across multiple energy domains. Power electronics system modelling is usually the study of the interplay of energy in the domains of electrical, mechanical, magnetic and thermal. The usefulness of bond graph modelling in power electronic field has been realised by researchers. Consequently in the last couple of decades, there has been a steadily increasing effort in developing simulation tools for bond graph modelling that are specially suited for power electronic study. For modelling rotating magnetic fields in electromagnetic machine models, a support for vector variables is essential. Unfortunately, all bond graph simulation tools presently provide support only for scalar variables. We propose an approach to provide complex variable and vector support to bond graph such that it will enable modelling of polyphase electromagnetic and spatial vector systems. We also introduced a rotary gyrator element and use it along with the switched junction for developing the complex/vector variable's toolbox. This approach is implemented by developing a complex S-function tool box in Simulink inside a MATLAB environment This choice has been made so as to synthesise the speed of S-function, the user friendliness of Simulink and the popularity of MATLAB.

Chemical bond approach to determining conductivity band gaps in amorphous chalcogenides and pnictides

The minimum energy required for the formation of conjugate pair of charged defects is found to be approximately equal to the experimental activation energy for d.c. conductivity in a number of amorphous chalcoganides and pnictides. This observation implies that the defect pair formation energy represents an intrinsic gap for transport in amorphous chalcogenides.