Computer simulation of ordering and dynamics in liquid crystals in the bulk and close to the surface

The aim of this PhD thesis is to investigate the orientational and dynamical properties of liquid crystalline systems, at molecular level and using atomistic computer simulations, to reach a better understanding of material behavior from a microscopic point view. In perspective this should allow to clarify the relation between the micro and macroscopic properties with the objective of predicting or confirming experimental results on these systems. In this context, we developed four different lines of work in the thesis. The first one concerns the orientational order and alignment mechanism of rigid solutes of small dimensions dissolved in a nematic phase formed by the 4-pentyl,4 cyanobiphenyl (5CB) nematic liquid crystal. The orientational distribution of solutes have been obtained with Molecular Dynamics Simulation (MD) and have been compared with experimental data reported in literature. we have also verified the agreement between order parameters and dipolar coupling values measured in NMR experiments. The MD determined effective orientational potentials have been compared with the predictions of Maier­Saupe and Surface tensor models. The second line concerns the development of a correct parametrization able to reproduce the phase transition properties of a prototype of the oligothiophene semiconductor family: sexithiophene (T6). T6 forms two crystalline polymorphs largely studied, and possesses liquid crystalline phases still not well characterized, From simulations we detected a phase transition from crystal to liquid crystal at about 580 K, in agreement with available experiments, and in particular we found two LC phases, smectic and nematic. The crystal­smectic transition is associated to a relevant density variation and to strong conformational changes of T6, namely the molecules in the liquid crystal phase easily assume a bent shape, deviating from the planar structure typical of the crystal. The third line explores a new approach for calculating the viscosity in a nematic through a virtual exper- iment resembling the classical falling sphere experiment. The falling sphere is replaced by an hydrogenated silicon nanoparticle of spherical shape suspended in 5CB, and gravity effects are replaced by a constant force applied to the nanoparticle in a selected direction. Once the nanoparticle reaches a constant velocity, the viscosity of the medium can be evaluated using Stokes' law. With this method we successfully reproduced experimental viscosities and viscosity anisotropy for the solvent 5CB. The last line deals with the study of order induction on nematic molecules by an hydrogenated silicon surface. Gaining predicting power for the anchoring behavior of liquid crystals at surfaces will be a very desirable capability, as many properties related to devices depend on molecular organization close to surfaces. Here we studied, by means of atomistic MD simulations, the flat interface between an hydrogenated (001) silicon surface in contact with a sample of 5CB molecules. We found a planar anchoring of the first layers of 5CB where surface interactions are dominating with respect to the mesogen intermolecular interactions. We also analyzed the interface 5CB­vacuum, finding a homeotropic orientation of the nematic at this interface.

Intrinsic uncoupling in the ATP synthase of Escherichia coli. Studies on WT and ε-truncated mutants

The H+/ATP ratio in the catalysis of ATP synthase has generally been considered a fixed parameter. However, Melandri and coworkers have recently shown that, in the ATP synthase of the photosynthetic bacterium Rb.capsulatus, this ratio can significantly decrease during ATP hydrolysis when the concentration of either ADP or Pi is maintained at a low level (Turina et al., 2004). The present work has dealt with the ATP synthase of E.coli, looking for evidence of this phenomenon of intrinsic uncoupling in this organism as well. First of all, we have shown that the DCCD-sensitive ATP hydrolysis activity of E.coli internal membranes was strongly inhibited by ADP and Pi, with a half-maximal effect in the submicromolar range for ADP and at 140 µM for Pi. In contrast to this monotonic inhibition, however, the proton pumping activity of the enzyme, as estimated under the same conditions by the fluorescence quenching of the ΔpH-sensitive probe ACMA, showed a clearly biphasic progression, both for Pi, increasing from 0 up to approximately 200 µM, and for ADP, increasing from 0 up to a few µM. We have interpreted these results as indicating that the occupancy of ADP and Pi binding sites shifts the enzyme from a partially uncoupled state to a fully coupled state, and we expect that the ADP- and Pi-modulated intrinsic uncoupling is likely to be a general feature of prokaryotic ATP synthases. Moreover, the biphasicity of the proton pumping data suggested that two Pi binding sites are involved. In order to verify whether the same behaviour could be observed in the isolated enzyme, we have purified the ATP synthase of E.coli and reconstituted it into liposomes. Similarly as observed in the internal membrane preparation, in the isolated and reconstituted enzyme it was possible to observe inhibition of the hydrolytic activity by ADP and Pi (with half-maximal effects at few µM for ADP and at 400 µM for Pi) with a concomitant stimulation of proton pumping. Both the inhibition of ATP hydrolysis and the stimulation of proton pumping as a function of Pi were lost upon ADP removal by an ADP trap. These data have made it possible to conclude that the results obtained in E.coli internal membranes are not due to the artefactual interference of enzymatic activities other than the ones of the ATP synthase. In addition, data obtained with liposomes have allowed a calibration of the ACMA signal by ΔpH transitions of known extent, leading to a quantitative evaluation of the proton pumping data. Finally, we have focused our efforts on searching for a possible structural candidate involved in the phenomenon of intrinsic uncoupling. The ε-subunit of the ATP-synthase is known as an endogenous inhibitor of the hydrolysis activity of the complex and appears to undergo drastic conformational changes between a non-inhibitory form (down-state) and an inhibitory form (up-state)(Rodgers & Wilce, 2000; Gibbons et al., 2000). In addition, the results of Cipriano & Dunn (2006) indicated that the C-terminal domain of this subunit played an important role in the coupling mechanism of the pump, and those of Capaldi et al. (2001), Suzuki et al. (2003) were consistent with the down-state showing a higher hydrolysis-to-synthesis ratio than the up-state. Therefore, we decided to search for modulation of pumping efficiency in a C-terminally truncated ε mutant. A low copy number expression vector has been built, carrying an extra copy of uncC, with the aim of generating an ε-overexpressing E.coli strain in which normal levels of assembly of the mutated ATP-synthase complex would be promoted. We have then compared the ATP hydrolysis and the proton pumping activity in membranes prepared from these ε-overexpressing E.coli strains, which carried either the WT ε subunit or the ε88-stop truncated form. Both strains yielded well energized membranes. Noticeably, they showed a marked difference in the inhibition of hydrolysis by Pi, this effect being largely lost in the truncated mutant. However, pre-incubation of the mutated enzyme with ADP at low nanomolar concentrations (apparent Kd = 0.7nM) restored the hydrolysis inhibition, together with the modulation of intrinsic uncoupling by Pi, indicating that, contrary to wild-type, during membrane preparation the truncated mutant had lost the ADP bound at this high-affinity site, evidently due to a lower affinity (and/or higher release) for ADP of the mutant relative to wild type. Therefore, one of the effects of the C-terminal domain of ε appears to be to modulate the affinity of at least one of the binding sites for ADP. The lack of this domain does not appear so much to influence the modulability of coupling efficiency, but instead the extent of this modulation. At higher preincubated ADP concentrations (apparent Kd = 117nM), the only observed effects were inhibition of both hydrolysis and synthesis, providing a direct proof that two ADP-binding sites on the enzyme are involved in the inhibition of hydrolysis, of which only the one at higher affinity also modulates the coupling efficiency.

Contribution of Non-Covalent Interactions and Electronic Effects on the Conformational Landscape and Tautomeric Equilibria of Molecules and Molecular Complexes: Structural and Dynamical Data from Rotational Spectroscopy

This thesis concerns the study of complex conformational surfaces and tautomeric equilibria of molecules and molecular complexes by quantum chemical methods and rotational spectroscopy techniques. In particular, the focus of this research is on the effects of substitution and noncovalent interactions in determining the energies and geometries of different conformers, tautomers or molecular complexes. The Free-Jet Absorption Millimeter Wave spectroscopy and the Pulsed-Jet Fourier Transform Microwave spectroscopy have been applied to perform these studies and the obtained results showcase the suitability of these techniques for the study of conformational surfaces and intermolecular interactions. The series of investigations of selected medium-size molecules and complexes have shown how different instrumental setups can be used to obtain a variety of results on molecular properties. The systems studied, include molecules of biological interest such as anethole and molecules of astrophysical interest such as N-methylaminoethanol. Moreover halogenation effects have been investigated on halogen substituted tautomeric systems (5-chlorohydroxypyridine and 6-chlorohydroxypyridine), where it has shown that the position of the inserted halogen atom affects the prototropic equilibrium. As for fluorination effects, interesting results have been achieved investigating some small complexes where a molecule of water is used as a probe to reveal the changes on the electrostatic potential of different fluorinated compounds: 2-fluoropyridine, 3-fluoropyridine and penta-fluoropyridine. While in the case of the molecular complex between water and 2-fluoropyridine and 3-fluoropyridine the geometry of the complex with one water molecule is analogous to that of pyridine with the water molecule linked to the pyridine nitrogen, the case of pentafluoropyridine reveals the effect of perfluorination and the water oxygen points towards the positive center of the pyridine ring. Additional molecular adducts with a molecule of water have been analyzed (benzylamine-water and acrylic acid-water) in order to reveal the stabilizing driving forces that characterize these complexes.