Nuclear magnetic resonance studies of molecular interactions and structures

This thesis is concerned with the investigation, by nuclear magnetic resonance spectroscopy, of the molecular interactions occurring in mixtures of benzene and cyclohexane to which either chloroform or deutero-chloroform has been added. The effect of the added polar molecule on the liquid structure has been studied using spin-lattice relaxation time, 1H chemical shift, and nuclear Overhauser effect measurements. The main purpose of the work has been to validate a model for molecular interaction involving local ordering of benzene around chloroform. A chemical method for removing dissolved oxygen from samples has been developed to encompass a number of types of sample, including quantitative mixtures, and its supremacy over conventional deoxygenation technique is shown. A set of spectrometer conditions, the use of which produces the minimal variation in peak height in the steady state, is presented. To separate the general diluting effects of deutero-chloroform from its effects due to the production of local order a series of mixtures involving carbon tetrachloride, instead of deutero-chloroform, have been used as non-interacting references. The effect of molecular interaction is shown to be explainable using a solvation model, whilst an approach involving 1:1 complex formation is shown not to account for the observations. It is calculated that each solvation shell, based on deutero-chloroform, contains about twelve molecules of benzene or cyclohexane. The equations produced to account for the T1 variations have been adapted to account for the 1H chemical shift variations in the same system. The shift measurements are shown to substantiate the solvent cage model with a cage capacity of twelve molecules around each chloroform molecule. Nuclear Overhauser effect data have been analysed quantitatively in a manner consistent with the solvation model. The results show that discrete shells only exist when the mole fraction of deutero-chloroform is below about 0.08.

Studies of the effects of molecular encounters on N.M.R. spin-lattice relaxation times

This thesis is concerned with investigations of the effects of molecular encounters on nuclear magnetic resonance spin-lattice relaxation times, with particular reference to mesitylene in mixtures with cyclohexane and TMS. The purpose of the work was to establish the best theoretical description of T1 and assess whether a recently identified mechanism (buffeting), that influences n.m.r. chemical shifts, governs Tl also. A set of experimental conditions are presented that allow reliable measurements of Tl and the N. O. E. for 1H and 13C using both C. W. and F.T. n.m.r. spectroscopy. Literature data for benzene, cyclohexane and chlorobenzene diluted by CC14 and CS2 are used to show that the Hill theory affords the best estimation of their correlation times but appears to be mass dependent. Evaluation of the T1 of the mesitylene protons indicates that a combined Hill-Bloembergen-Purcell-Pound model gives an accurate estimation of T1; subsequently this was shown to be due to cancellation of errors in the calculated intra and intemolecular components. Three experimental methods for the separation of the intra and intermolecular relaxation times are described. The relaxation times of the 13C proton satellite of neat bezene, 1,4 dioxane and mesitylene were measured. Theoretical analyses of the data allow the calculation of Tl intra. Studies of intermolecular NOE's were found to afford a general method of separating observed T1's into their intra and intermolecular components. The aryl 1H and corresponding 13C T1 values and the NOE for the ring carbon of mesitylene in CC14 and C6H12-TMS have been used in combination to determine T1intra and T1inter. The Hill and B.P.P. models are shown to predict similarly inaccurate values for T1linter. A buffeting contribution to T1inter is proposed which when applied to the BPP model and to the Gutowsky-Woessner expression for T1inter gives an inaccuracy of 12% and 6% respectively with respect to theexperimentally based T1inter.

Ionic and molecular diffusion in cementitious materials

The work described in this thesis is an attempt to provide improved understanding of the effects of several factors affecting diffusion in hydrated cement pastes and to aid the prediction of ionic diffusion processes in cement-based materials. Effect of pore structure on diffusion was examined by means of comparative diffusion studies of quaternary ammonium ions with different ionic radii. Diffusivities of these ions in hydrated pastes of ordinary portland cement with or without addition of fly ash were determined by a quasi-steady state technique. The restriction of the pore geometry on diffusion was evaluated from the change of diffusivity in response to the change of ionic radius. The pastes were prepared at three water-cement ratios, 0.35, 0.50 and 0.65. Attempts were made to study the effect of surface charge or the electrochemical double layer at the pore/solution interface on ionic diffusion. An approach was to evaluate the zeta potentials of hydrated cement pastes through streaming potential measurements. Another approach was the comparative studies of the diffusion kinetics of chloride and dissolved oxygen in hydrated pastes of ordinary portland cement with addition of 0 and 20% fly ash. An electrochemical technique for the determination of oxygen diffusivity was also developed. Non-steady state diffusion of sodium potassium, chloride and hydroxyl ions in hydrated ordinary portland cement paste of water-cement ratio 0.5 was studied with the aid of computer-modelling. The kinetics of both diffusion and ionic binding were considered for the characterization of the concentration profiles by Fick's first and second laws. The effect of the electrostatic interactions between ions on the overall diffusion rates was also considered. A general model concerning the prediction of ionic diffusion processes in cement-based materials has been proposed.

Toward an atomistic understanding of the immune synapse:large-scale molecular dynamics simulation of a membrane-embedded TCR-pMHC-CD4 complex

T-cell activation requires interaction of T-cell receptors (TCR) with peptide epitopes bound by major histocompatibility complex (MHC) proteins. This interaction occurs at a special cell-cell junction known as the immune or immunological synapse. Fluorescence microscopy has shown that the interplay among one agonist peptide-MHC (pMHC), one TCR and one CD4 provides the minimum complexity needed to trigger transient calcium signalling. We describe a computational approach to the study of the immune synapse. Using molecular dynamics simulation, we report here on a study of the smallest viable model, a TCR-pMHC-CD4 complex in a membrane environment. The computed structural and thermodynamic properties are in fair agreement with experiment. A number of biomolecules participate in the formation of the immunological synapse. Multi-scale molecular dynamics simulations may be the best opportunity we have to reach a full understanding of this remarkable supra-macromolecular event at a cell-cell junction.

Reduction of circulating soluble Flt-1 alleviates preeclampsia-like symptoms in a mouse model

Preeclampsia (PE) is characterized by widespread endothelial damage with hypertension, proteinuria, glomeruloendotheliosis and elevated soluble Flt-1 (sFlt-1), a natural occurring antagonist of vascular endothelial growth factor (VEGF). Cancer patients receiving anti-VEGF therapy exhibit similar symptoms. We suggested that a decrease in circulating sFlt-1 would alleviate the symptoms associated with PE. Adenoviral (Adv) overexpression of sFlt-1 induced proteinuria, caused glomerular damage and increase in blood pressure in female Balb/c mice. Circulating level of sFlt-1 above 50 ng/ml plasma induced severe vascular damage and glomerular endotheliosis. Albumin concentration in urine was elevated up to 30-fold, compared to control AdvGFP-treated animals. The threshold of kidney damage was in the range of 20-30 ng/ml sFlt-1 in plasma (8-15 ng/ml in urine). Co-administration of AdvsFlt-1 with AdvVEGF to neutralize circulating sFlt-1 resulted in more than a 70% reduction in free sFlt-1 in plasma, more than 80% reduction in urine and rescued the damaging effect of sFlt-1 on the kidneys. This demonstrates that below a critical threshold sFlt-1 fails to elicit damage to the fenestrated endothelium and that co-expression of VEGF is able to rescue effects mediated by sFlt-1 overexpression.

Molecular dynamics study of solvent transport in nanoscale osmosis

An ideal of osmotic equilibrium between an ideal solution and pure solvent separated by a semi-permeable membrane is studied numerically using the method of molecular dynamics. The osmotic flow is observed as the inflow of the solvent across the membrane from the dilute to the concentrated side. The validity of van't Hoff's law for osmotic pressure is confirmed over a wide range of concentrations. It is found that the law is established by a balance between non-uniform partial pressures of solute and solvent. Furthermore, the present model permits an understanding of the mechanism of the osmotic flow in the relaxation process as the liquids evolve from the initial state to the equilibrium state. We focus in particular on the interaction between solute and solvent. ©2008 The Physical Society of Japan.

Uncatalysed and potassium-catalysed pyrolysis of the cell-wall constituents of biomass and their model compounds

Cell-wall components (cellulose, hemicellulose (oat spelt xylan), lignin (Organosolv)), and model compounds (levoglucosan (an intermediate product of cellulose decomposition) and chlorogenic acid (structurally similar to lignin polymer units)) have been investigated to probe in detail the influence of potassium on their pyrolysis behaviours as well as their uncatalysed decomposition reaction. Cellulose and lignin were pretreated to remove salts and metals by hydrochloric acid, and this dematerialized sample was impregnated with 1% of potassium as potassium acetate. Levoglucosan, xylan and chlorogenic acid were mixed with CHCOOK to introduce 1% K. Characterisation was performed using thermogravimetric analysis (TGA) and differential thermal analysis (DTA). In addition to the TGA pyrolysis, pyrolysis-gas chromatography-mass spectrometry (PY-GC-MS) analysis was introduced to examine reaction products. Potassium-catalysed pyrolysis has a huge influence on the char formation stage and increases the char yields considerably (from 7.7% for raw cellulose to 27.7% for potassium impregnated cellulose; from 5.7% for raw levoglucosan to 20.8% for levoglucosan with CHCOOK added). Major changes in the pyrolytic decomposition pathways were observed for cellulose, levoglucosan and chlorogenic acid. The results for cellulose and levoglucosan are consistent with a base catalysed route in the presence of the potassium salt which promotes complete decomposition of glucosidic units by a heterolytic mechanism and favours its direct depolymerization and fragmentation to low molecular weight components (e.g. acetic acid, formic acid, glyoxal, hydroxyacetaldehyde and acetol). Base catalysed polymerization reactions increase the char yield. Potassium-catalysed lignin pyrolysis is very significant: the temperature of maximum conversion in pyrolysis shifts to lower temperature by 70 K and catalysed polymerization reactions increase the char yield from 37% to 51%. A similar trend is observed for the model compound, chlorogenic acid. The addition of potassium does not produce a dramatic change in the tar product distribution, although its addition to chlorogenic acid promoted the generation of cyclohexane and phenol derivatives. Postulated thermal decomposition schemes for chlorogenic acid are presented. © 2008 Elsevier B.V. All rights reserved.

Tissue-specific selection of reference genes is required for expression studies in the mouse model of maternal protein undernutrition

Suboptimal maternal nutrition during gestation results in the establishment of long-term phenotypic changes and an increased disease risk in the offspring. To elucidate how such environmental sensitivity results in physiological outcomes, the molecular characterisation of these offspring has become the focus of many studies. However, the likely modification of key cellular processes such as metabolism in response to maternal undernutrition raises the question of whether the genes typically used as reference constants in gene expression studies are suitable controls. Using a mouse model of maternal protein undernutrition, we have investigated the stability of seven commonly used reference genes (18s, Hprt1, Pgk1, Ppib, Sdha, Tbp and Tuba1) in a variety of offspring tissues including liver, kidney, heart, retro-peritoneal and inter-scapular fat, extra-embryonic placenta and yolk sac, as well as in the preimplantation blastocyst and blastocyst-derived embryonic stem cells. We find that although the selected reference genes are all highly stable within this system, they show tissue, treatment and sex-specific variation. Furthermore, software-based selection approaches rank reference genes differently and do not always identify genes which differ between conditions. Therefore, we recommend that reference gene selection for gene expression studies should be thoroughly validated for each tissue of interest. © 2011 Elsevier Inc.

A novel concentration dependent amino acid ion pair strategy to mediate drug permeation using indomethacin as a model insoluble drug

Assessment of oral drug bioavailability is an important parameter for new chemical entities (NCEs) in drug development cycle. After evaluating the pharmacological response of these new molecules, the following critical stage is to investigate their in vitro permeability. Despite the great success achieved by prodrugs, covalent linking the drug molecule with a hydrophobic moiety might result in a new entity that might be toxic or ineffective. Therefore, an alternative that would improve the drug uptake without affecting the efficacy of the drug molecule would be advantageous. The aim of the current study is to investigate the effect of ion-pairing on the permeability profile of a model drug: indomethacin (IND) to understand the mechanism behind the permeability improvement across Caco-2 monolayers. Arginine and lysine formed ion-pairs with IND at various molar ratios 1:1, 1:2, 1:4 and 1:8 as reflected by the double reciprocal graphs. The partitioning capacities of the IND were evaluated using octanol/water partitioning studies and the apparent permeabilities (P app) were measured across Caco-2 monolayers for the different formulations. Partitioning studies reflected the high hydrophobicity of IND (Log P = 3) which dropped upon increasing the concentrations of arginine/lysine in the ion pairs. Nevertheless, the prepared ion pairs improved IND permeability especially after 60 min of the start of the experiment. Coupling partitioning and permeability results suggest a decrease in the passive transcellular uptake due to the drop in IND portioning capacities and a possible involvement of active carriers. Future work will investigate which transport gene might be involved in the absorption of the ion paired formulations using molecular biology technologies. © 2014 Elsevier B.V. All rights reserved.

Adsorption of dimethyl ether (DME) on zeolite molecular sieves

In recent years there has been growing interest in the use of dimethyl ether (DME) as an alternative fuel. In this study, the adsorption of DME on molecular sieves 4Å (Mol4A) and 5Å (Mol5A) has been experimentally investigated using the volumetric adsorption method. Data on the adsorption isotherms, heats of adsorption, and adsorption kinetic have been obtained and used to draw conclusions and compare the performance of the two adsorbents. Within the conditions considered, the adsorption capacity of Mol5A was found to be around eight times higher than the capacity of Mol4A. Low temperature adsorption and thermal pre-treatment of the adsorbents in vacuum were observed to be favourable for increased adsorption capacity. The adsorption isotherms for both adsorbent were fitted to the Freundlich model and the corresponding model parameters are proposed. The adsorption kinetic analysis suggest that the DME adsorption on Mol5A is controlled by intracrystalline diffusion resistance, while on Mol4A it is mainly controlled by surface layering resistance with the diffusion only taking place at the start of adsorption and for a very limited short time. The heats of adsorption were calculated by a calorimetric method based on direct temperature measurements inside the adsorption cell. Isosteric heats, calculated by the thermodynamic approach (Clasius-Clapeyron equation), have consistently shown lower values. The maximum heat of adsorption was found to be 25.9kJmol-1 and 20.1kJmol-1 on Mol4A and Mol5A, respectively; thus indicating a physisorption type of interactions. © 2014 Elsevier B.V.

Evidence-based nanoscopic and molecular framework for excipient functionality in compressed orally disintegrating tablets

The work investigates the adhesive/cohesive molecular and physical interactions together with nanoscopic features of commonly used orally disintegrating tablet (ODT) excipients microcrystalline cellulose (MCC) and D-mannitol. This helps to elucidate the underlying physico-chemical and mechanical mechanisms responsible for powder densification and optimum product functionality. Atomic force microscopy (AFM) contact mode analysis was performed to measure nano-adhesion forces and surface energies between excipient-drug particles (6-10 different particles per each pair). Moreover, surface topography images (100 nm2-10 μm2) and roughness data were acquired from AFM tapping mode. AFM data were related to ODT macro/microscopic properties obtained from SEM, FTIR, XRD, thermal analysis using DSC and TGA, disintegration testing, Heckel and tabletability profiles. The study results showed a good association between the adhesive molecular and physical forces of paired particles and the resultant densification mechanisms responsible for mechanical strength of tablets. MCC micro roughness was 3 times that of D-mannitol which explains the high hardness of MCC ODTs due to mechanical interlocking. Hydrogen bonding between MCC particles could not be established from both AFM and FTIR solid state investigation. On the contrary, D-mannitol produced fragile ODTs due to fragmentation of surface crystallites during compression attained from its weak crystal structure. Furthermore, AFM analysis has shown the presence of extensive micro fibril structures inhabiting nano pores which further supports the use of MCC as a disintegrant. Overall, excipients (and model drugs) showed mechanistic behaviour on the nano/micro scale that could be related to the functionality of materials on the macro scale. © 2014 Al-khattawi et al.

Bridging large and small scales of water models using hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework

This thesis presents a two-dimensional water model investigation and development of a multiscale method for the modelling of large systems, such as virus in water or peptide immersed in the solvent. We have implemented a two-dimensional ‘Mercedes Benz’ (MB) or BN2D water model using Molecular Dynamics. We have studied its dynamical and structural properties dependence on the model’s parameters. For the first time we derived formulas to calculate thermodynamic properties of the MB model in the microcanonical (NVE) ensemble. We also derived equations of motion in the isothermal–isobaric (NPT) ensemble. We have analysed the rotational degree of freedom of the model in both ensembles. We have developed and implemented a self-consistent multiscale method, which is able to communicate micro- and macro- scales. This multiscale method assumes, that matter consists of the two phases. One phase is related to micro- and the other to macroscale. We simulate the macro scale using Landau Lifshitz-Fluctuating Hydrodynamics, while we describe the microscale using Molecular Dynamics. We have demonstrated that the communication between the disparate scales is possible without introduction of fictitious interface or approximations which reduce the accuracy of the information exchange between the scales. We have investigated control parameters, which were introduced to control the contribution of each phases to the matter behaviour. We have shown, that microscales inherit dynamical properties of the macroscales and vice versa, depending on the concentration of each phase. We have shown, that Radial Distribution Function is not altered and velocity autocorrelation functions are gradually transformed, from Molecular Dynamics to Fluctuating Hydrodynamics description, when phase balance is changed. In this work we test our multiscale method for the liquid argon, BN2D and SPC/E water models. For the SPC/E water model we investigate microscale fluctuations which are computed using advanced mapping technique of the small scales to the large scales, which was developed by Voulgarakisand et. al.

Modeling the peptide-T cell receptor interaction by the comparative molecular similarity indices analysis-soft independent modeling of class analogy technique

A set of 38 epitopes and 183 non-epitopes, which bind to alleles of the HLA-A3 supertype, was subjected to a combination of comparative molecular similarity indices analysis (CoMSIA) and soft independent modeling of class analogy (SIMCA). During the process of T cell recognition, T cell receptors (TCR) interact with the central section of the bound nonamer peptide; thus only positions 4−8 were considered in the study. The derived model distinguished 82% of the epitopes and 73% of the non-epitopes after cross-validation in five groups. The overall preference from the model is for polar amino acids with high electron density and the ability to form hydrogen bonds. These so-called “aggressive” amino acids are flanked by small-sized residues, which enable such residues to protrude from the binding cleft and take an active role in TCR-mediated T cell recognition. Combinations of “aggressive” and “passive” amino acids in the middle part of epitopes constitute a putative TCR binding motif

Non-Markov state model of peptide dynamics

A hidden Markov state model has been applied to classical molecular dynamics simulated small peptide in explicit water. The methodology allows increasing the time resolution of the model and describe the dynamics with the precision of 0.3 ps (comparing to 6 ps for the standard methodology). It also permits the investigation of the mechanisms of transitions between the conformational states of the peptide. The detailed description of one of such transitions for the studied molecule is presented. © 2012 Elsevier B.V. All rights reserved.

A hybrid molecular dynamics/fluctuating hydrodynamics method for modelling liquids at multiple scales in space and time

A new 3D implementation of a hybrid model based on the analogy with two-phase hydrodynamics has been developed for the simulation of liquids at microscale. The idea of the method is to smoothly combine the atomistic description in the molecular dynamics zone with the Landau-Lifshitz fluctuating hydrodynamics representation in the rest of the system in the framework of macroscopic conservation laws through the use of a single "zoom-in" user-defined function s that has the meaning of a partial concentration in the two-phase analogy model. In comparison with our previous works, the implementation has been extended to full 3D simulations for a range of atomistic models in GROMACS from argon to water in equilibrium conditions with a constant or a spatially variable function s. Preliminary results of simulating the diffusion of a small peptide in water are also reported.