Ordering of apolar and polar solutes in nematic solvents

Prediction of anisotropic NMR data without knowledge of alignment medium structure by surface decomposition

Prediction of anisotropic NMR data directly from solute-medium interaction is of significant theoretical and practical interest, particularly for structure elucidation, configurational analysis and conformational studies of complex organic molecules and natural products. Current prediction methods require an explicit structural model of the alignment medium: a requirement either impossible or impractical on a scale necessary for small organic molecules. Here we formulate a comprehensive mathematical framework for a parametrization protocol that deconvolutes an arbitrary surface of the medium into several simple local landscapes that are distributed over the medium's surface by specific orientational order parameters. The shapes and order parameters of these local landscapes are determined via fitting that maximizes the congruence between experimentally determined anisotropic NMR measurables and their predicted counterparts, thus avoiding the need for an a priori knowledge of the global medium morphology. This method achieves substantial improvements in the accuracy of predicted anisotropic NMR values compared to current methods, as demonstrated herein with sixteen natural products. Furthermore, because this formalism extracts structural commonalities of the medium by combining anisotropic NMR data from different compounds, its robustness and accuracy are expected to improve as more experimental data become available for further re-optimization of fitting parameters.

Effect of the Functional Groups of Racemic Rodlike Schiff Base Mesogens on the Stabilization of Blue Phase in Binary Mixture Systems

Four series of rodlike racemic Schiff base mesogens possessing different alkyl chains and two types of linkages, ester and alkynyl linkages, were synthesized and applied to induce cubic blue phases (BPs) in simple binary mixture systems. The mesophases of these Schiff base mesogens were confirmed by variable-temperature X-ray diffraction and the characteristic texture from polarized optical microscopy (POM). In general, when chiral additive S-(+)-2-octyl 4-(4-hexyloxybenzoyloxy)benzoate (S811; 20-40 wt %) is added into the rodlike racemic salicylaldimine-based mesogens, the cubic BPs could be observed and its temperature range is larger than 20 K. The widest temperature range of the cubic BP (35 K) can be observed in the blending mixture composed of rodlike racemic salicylaldimine-based mesogen OH-TI_n possessing alkynyl linkage and 35-40 wt % S811. However, Schiff base mesogens possessing alkynyl linkage show a direct isotropic to chiral nematic transition when equal amount of chiral dopant is added. Notably, the termination temperature of BPs is very close to room temperature (ca. 35 °C) after 40.0 wt % S811 is added into the salicylaldimine-based mesogens possessing terminal alkyl chains and ester linkage. Interestingly, wide BPs (>30 K) can also be induced by adding chiral additive 1,4:3,6-dianhydro-2,5-bis[4-(n-hexyl-1-oxy)benzoic acid]sorbitol (ISO(6OBA)₂) with a high helical twisting power into the racemic Schiff base mesogen possessing ester linkage. Cubic BPI and BPII can be confirmed by reflectance spectra and POM. The results of reflectance spectra indicate that the binary mixture composed of salicylaldimine-based mesogens and S811 easily exhibits a supercooling effect and induces BPI. However, only BPII can be observed in all binary mixtures containing Schiff base mesogens. On the basis of our experimental results and molecular modeling, we suppose that the values of biaxiality, polarizability, and the dipole moment of molecular geometry are the main factors that affect BP stabilization.

Smectic order parameters via liquid crystal NMR spectroscopy: Application to a partial bilayer smectic A phase

Solute molecules were dissolved in the liquid crystal 4-cyano-4'-n-octyloxybiphenyl (8OCB), known to form a partial bilayer smectic-A phase. Through measurement of solutes' and solvent's orientational order parameters via nuclear magnetic resonance spectroscopy, and their analysis via a statistical thermodynamic density functional theory, values of the solvent's positional order parameters and solutes' positional-orientational distribution functions were obtained. Near to the transition to the nematic phase, the main positional order parameter of the smectic liquid crystal turned out to be comprised in the interval 0.4-0.6, though the quality of the fittings assuming the phase as nematic all across the temperature range investigated was only slightly worse. This may be ascribed to the looseness of the partial bilayer smectic structure. Solutes were found to preferentially lie in those regions where liquid crystal molecule terminal chains are located.

Synthesis and evaluation of 1,5-disubstituted tetrazoles as rigid analogues of combretastatin A-4 with potent antiproliferative and antitumor activity

Tubulin, the major structural component of microtubules, is a target for the development of anticancer agents. Two series of 1,5-diaryl substituted 1,2,3,4-tetrazoles were concisely synthesized, using a palladium-catalyzed cross-coupling reaction, and identified as potent antiproliferative agents and novel tubulin polymerization inhibitors that act at the colchicine site. SAR analysis indicated that compounds with a 4-ethoxyphenyl group at the N-1 or C-5 position of the 1,2,3,4-tetrazole ring exhibited maximal activity. Several of these compounds also had potent activity in inhibiting the growth of multidrug resistant cells overexpressing P-glycoprotein. Active compounds induced apoptosis through the mitochondrial pathway with activation of caspase-9 and caspase-3. Furthermore, compound 4l significantly reduced in vivo the growth of the HT-29 xenograft in a nude mouse model, suggesting that 4l is a promising new antimitotic agent with clinical potential.

Domain cooperativity in multidomain proteins: what can we learn from molecular alignment in anisotropic media?

Many proteins have modular design with multiple globular domains connected via flexible linkers. As a simple model of such system, we study a tandem construct consisting of two identical SH3 domains and a variable-length Gly/Ser linker. When the linker is short, this construct represents a dumbbell-shaped molecule with limited amount of domain-domain mobility. Due to its elongated shape, this molecule efficiently aligns in steric alignment media. As the length of the linker increases, the two domains become effectively uncoupled and begin to behave as independent entities. Consequently, their degree of alignment drops, approaching that found in the (near-spherical) isolated SH3 domains. To model the dependence of alignment parameters on the length of the interdomain linker, we have generated in silico a series of conformational ensembles representing SH3 tandems with different linker length. These ensembles were subsequently used as input for alignment prediction software PALES. The predicted alignment tensors were compared with the results of experimental measurements using a series of tandem-SH3 samples in PEG/hexanol alignment media. This comparison broadly confirmed the expected trends. At the same time, it has been found that the isolated SH3 domain aligns much stronger than expected. This finding can be attributed to complex morphology of the PEG/hexanol media and/or to weak site-specific interactions between the protein and the media. In the latter case, there are strong indications that electrostatic interactions may play a role. The fact that PEG/hexanol does not behave as a simple steric media should serve as a caution for studies that use PALES as a quantitative prediction tool (especially for disordered proteins). Further progress in this area depends on our ability to accurately model the anisotropic media and its site-specific interactions with protein molecules. Once this ability is improved, it should be possible to use the alignment parameters as a measure of domain-domain cooperativity, thus identifying the situations where two domains transiently interact with each other or become coupled through a partially structured linker.

Orientational mechanisms in liquid crystalline systems. 1. A reaction field analytical description of the interaction between the electric quadrupole moment of a probe solute and the electric field gradient of a nematic solvent

In the present paper, the fundamental problem of calculating the electric field gradient (EFG) experienced by a highly idealized solute, represented by a general point quadrupole immersed in an anisotropic uniaxial medium, has been tackled. Following a generalized reaction field approach (based upon the original ideas and the "mean-field philosophy" due to Kirkwood and Onsager) in the linear response approximation, a closed analytical expression of the EFG has been derived (to the best of our knowledge, for the first time). The obtained expression is particularly simple and elegant, also thanks to the oversimplifying approximation that the virtual cavity containing the solute is assumed to be perfectly spherical. This compact and manageable formula, obtained by a rigorous mathematical derivation (unlike other mean-field phenomenological models previously suggested in literature) can be useful to investigate and better understand a likely orientational mechanism, partly responsible for the ordering of small solutes dissolved in nematic mesophases, based on the interaction between the electric quadrupole of the solute and the electric field gradient of the anisotropic uniaxial medium (in the next paper of this issue, the formulation obtained in this work is widely tested on a variety of uniaxial and biaxial solutes dissolved in different nematic solvents).

Orientational mechanisms in liquid crystalline systems. 2. The contribution to solute ordering from the reaction field interaction between the solute electric quadrupole moment and the solvent electric field gradient

In the previous paper of this issue, [Celebre, G.; Ionescu, A. J. Phys. Chem. B doi: 10.1021/jp907310g], following a generalized reaction field approach in the linear response approximation, we were successful in obtaining an analytical compact expression for the mean-field anisotropic orientational potential U(Q-EFG) theoretically experienced by a highly idealized nonionic and apolar solute, considered as a point quadrupole immersed in a uniaxial polarizable continuum medium (model of a nematic solvent comprised of dipolar mesogenic molecules). The term U(Q-EFG) describes the electrostatic interaction between the electric quadrupole of the solute and the electric field gradient induced at the solute by the surrounding medium polarized by the distribution of electric charges representing the quadrupolar solute itself. In the present paper, the obtained potential has been considered as an additional orientational interaction contributing to the solute ordering, besides the well-recognized and very effective "short-range" (size-and-shape-dictated) mechanisms. Since in our theory the solvent is characterized by its dielectric tensor, the model has been widely tested by taking as references the experimental order parameters of several uniaxial and biaxial different small rigid probe molecules (H(2), N(2), acetylene, allene, propyne, benzene, hexafluorobenzene, 1,4-difluorobenzene, and norbornadiene) dissolved in the nematic solvents ZLI1132 (Deltaepsilon >> 0) and EBBA (Deltaepsilon < 0); moreover, the order parameters of the same solutes in the so-called nematic "magic mixture" (45 wt % EBBA + 55 wt % ZLI1132), where the short-range orientational effects are commonly believed to be very dominant, have been conventionally assumed as reference of the absence of electrostatic orientational effects. The experimental order parameters of the treated solutes, obtained in the past by liquid crystal NMR and available from literature, have been then compared with those theoretically predicted by our theoretical approach in order to obtain useful hints about two basic points, (a) the real physical nature of the interactions (other than the "size-and-shape") involved in the orientational mechanisms and (b) the conceptual effectiveness of the suggested mean-field approach in describing this kind of phenomena. Successes and failures of the approach in the predictions are discussed at length, along with their possible reasons and implications.

Solvent smectic order parameters from solute nematic order parameters

In liquid crystals, while the second and fourth rank orientational order parameters characterizing a nematic phase can be experimentally determined via several techniques, there is no straightforward experiment rendering the positional order parameters characterizing a smectic A phase. This work illustrates a novel method to estimate the positional order parameters of a smectogenic liquid crystal solvent from knowledge of the orientational order parameters of a number of solutes dissolved therein. The latter order parameters can be experimentally determined via liquid crystal NMR spectroscopy. These data can be then analyzed with a statistical-thermodynamic density functional theory, whose basic ingredient is a model for solute-solvent intermolecular interactions. Its parametrization and the subsequent fitting procedure eventually permit one to obtain the positional order parameters of the solvent besides the positional-orientational distribution function of the solutes. The method is applied to the smectogen 4,4(')-di-n-heptyl-azoxybenzene, in which the solutes 1,4-dichlorobenzene and naphthalene have been dissolved. With the help of this exploratory practical example, pros and cons of the method are pointed out and further developments prospected.

Protein dynamics from NMR: the slowly relaxing local structure analysis compared with model-free analysis

(15)N-(1)H spin relaxation is a powerful method for deriving information on protein dynamics. The traditional method of data analysis is model-free (MF), where the global and local N-H motions are independent and the local geometry is simplified. The common MF analysis consists of fitting single-field data. The results are typically field-dependent, and multifield data cannot be fit with standard fitting schemes. Cases where known functional dynamics has not been detected by MF were identified by us and others. Recently we applied to spin relaxation in proteins the slowly relaxing local structure (SRLS) approach, which accounts rigorously for mode mixing and general features of local geometry. SRLS was shown to yield MF in appropriate asymptotic limits. We found that the experimental spectral density corresponds quite well to the SRLS spectral density. The MF formulas are often used outside of their validity ranges, allowing small data sets to be force-fitted with good statistics but inaccurate best-fit parameters. This paper focuses on the mechanism of force-fitting and its implications. It is shown that MF analysis force-fits the experimental data because mode mixing, the rhombic symmetry of the local ordering and general features of local geometry are not accounted for. Combined multifield multitemperature data analyzed with the MF approach may lead to the detection of incorrect phenomena, and conformational entropy derived from MF order parameters may be highly inaccurate. On the other hand, fitting to more appropriate models can yield consistent physically insightful information. This requires that the complexity of the theoretical spectral densities matches the integrity of the experimental data. As shown herein, the SRLS spectral densities comply with this requirement.

Assessment of solvent effects: do weak alignment media affect the structure of the solute?

Alignment media used for measuring residual dipolar couplings, such as solutions of filamentous phages, phospholipid mixtures, polyacrylamide gels and various lyotropic liquid crystalline systems were investigated with respect to solvent effects on molecular structure. Structural parameters of the small rigid model compound 13C-acetonitrile were calculated from dipolar couplings and variations from expectation values were used for assessment of solvent effects. Only minor solvent effects were observed for most of the media employed and the measured structural data are in good agreement with microwave data and theoretical predictions.

Expected and Unexpected Behavior of the Orientational Order and Dynamics Induced by Azobenzene Solutes in a Nematic

We have explored the changes in the phase stability, orientational order, and dynamics of the nematic 4-cyano-4'-n-pentylbiphenyl (5CB) doped with either the trans or the cis form of different p-azobenzene derivatives using the ESR spin-probe technique. In particular, we have studied the effects induced by each of the seven nonmesogenic 4-R-phenylazobenzenes (R = H, F, Br, CH3, CF3, On-Bu, Ot-Bu) at 1% and 7% mole fraction on the order parameter <P2> and on the shift of the nematic-isotropic transition temperature (TNI), as reported by a nitroxide spin probe, and we have tried to relate them to the solute shape and charge distribution. In all the cases the presence of the azo-derivative causes a depression of T(NI), more pronounced for the cis isomers. The dependence of <P2> on the reduced temperature T* = T/T(NI) remains the same as that of pure 5CB in all trans-doped samples at 1% and 7% and decreases only slightly in the cis at 1%. However, we observe different and in some cases large variations (up to 25%) in <P2> for the cis at 7%, showing solute effects that go beyond the shift in T(NI). Surprisingly enough, even at the highest concentration, the probe dynamics appears to be essentially independent of the nature, the configuration, and the concentration of the different solutes and very similar to that observed in the pure 5CB.

Methyl Dynamics in Proteins from NMR Slowly Relaxing Local Structure Spin Relaxation Analysis: A New Perspective

NMR spin relaxation of (2)H nuclei in (13)CH(2)D groups is a powerful method for studying side-chain motion in proteins. The analysis is typically carried out with the original model-free (MF) approach adapted to methyl dynamics. The latter is described in terms of axial local motions around, and of, the methyl averaging axis, mutually decoupled and independent of the global motion of the protein. Methyl motion is characterized primarily by the axial squared order parameter, S(axis)2, associated with fluctuations of the methyl averaging axis. This view is shown to be oversimplified by applying to typical experimental data the slowly relaxing local structure (SRLS) approach of Polimeno and Freed (Adv. Chem. Phys. 1993, 83, 89) which can be considered the generalization of the MF approach. Neglecting mode coupling and the asymmetry of the local ordering and treating approximately features of local geometry imply inaccurate values of S(axis)2, hence of the residual configurational entropy derived from it. S(axis)2, interpreted as amplitude of motion, was found to range from near disorder to almost complete order. Contrary to this picture, we find with the SRLS approach a moderate distribution in the magnitude of asymmetric local ordering and significant variation in its symmetry. The latter important property can be associated implicitly with the contribution of side-chain rotamer jumps. This is consistent with experimental residual dipolar coupling studies and theoretical work based on molecular dynamics simulations and molecular mechanics considerations. Configurational entropy is obtained in the SRLS approach directly from experimentally determined asymmetric potentials. Inconsistency between order parameters from 2H relaxation and from eta(HC-HH) cross-correlation and increase in order parameters with increasing temperature were observed with the MF approach. These discrepancies are reconciled, and physically tenable temperature dependence is obtained with the SRLS approach.

Solvation Dynamics by Computer Simulation: Coumarin C153 in 1,4-Dioxane

Computer simulation results are presented for an atomistic pair potential model of 1,4-dioxane which takes into account molecular flexibility. The model has been conceived to be applied to the study of solvatochromism and solvation dynamics in the presence of the polar probe coumarin C153. Computer simulations on the pure liquid have produced thermodynamical, structural, and dynamical data in good agreement with available experimental measures. This constitutes a valuable test of the 1,4-dioxane all-atom model employed. The study of solute-solvent interactions for C153 in 1,4-dioxane has been motivated by the aim of casting light, through simulations, on the interesting experimental findings according to which such a solvent behaves as a "polar" solvent with respect to dynamic solvation properties. Molecular dynamics is particularly suitable to model the process and provides an interpretation of the so-called "dioxane anomaly". An investigation of the structure of the solvation shell and of the dynamics of solvation is presented and discussed. In particular, the satisfactory accordance between simulated and experimental solvation response implies that the simulations give a reliable description of both solute and solvent at a molecular level and reinforces the idea that the explicit inclusion of discrete solvent molecules is needed for a realistic treatment of solvation phenomena in which the local structure of the liquid plays a key role.

Molecular Ordering and Structure of Quasi-spherical Solutes by Liquid Crystal NMR and Monte Carlo Simulations: The Case of Norbornadiene

Norbornadiene (a C2v symmetry bicyclic rigid hydrocarbon) dissolved in three different nematic mesophases has been studied by liquid crystal NMR, to contribute to a better understanding of the influence of solvents on the solute's ordering and structure. The main results achieved by this work can be summarized as follows: (i) the order parameters obtained by the analysis of the 1H NMR spectra (at different temperatures) were successfully reproduced by a recently proposed model of solute/liquid crystal interactions, by using Monte Carlo numerical simulations; (ii) the theoretical (B3LYP/6-31++G**) "equilibrium" geometry of norbornadiene, vibrationally corrected by using the force field calculated at the same level, is compatible (within, at most, a 5% error) with experimental LXNMR data. This leads to the conclusion that the structure is not significantly distorted by the tested solvents.

A Combined Theoretical and Experimental Approach to Determining Order Parameters of Solutes in Liquid Crystals from 13C NMR Data

The ordering properties of an anisotropic liquid crystal can be studied by recording 13C NMR spectra at different temperatures for a number of rigid solutes. The traditional difficulty in analyzing 13C data comes from the scarcity of experimental information about the carbon shielding tensors and from their limited transferability among different solutes. We show that these obstacles can be overcome by computing high-level ab initio shielding tensors, also including the solvent effects by the polarizable continuum model. The reliability of this combined approach is carefully verified, and the order parameters of several solutes are obtained by reanalyzing previously published spectra. The quality of the results is shown to be comparable to deuterium NMR without the need of isotopic substitution.

Prediction of charge-induced molecular alignment of biomolecules dissolved in dilute liquid-crystalline phases

Alignment of macromolecules in nearly neutral aqueous lyotropic liquid-crystalline media such as bicelles, commonly used in macromolecular NMR studies, can be predicted accurately by a steric obstruction model (Zweckstetter and Bax, 2000). A simple extension of this model is described that results in improved predictions for both the alignment orientation and magnitude of protein and DNA solutes in charged nematic media, such as the widely used medium of filamentous phage Pf1. The extended model approximates the electrostatic interaction between a solute and an ordered phage particle as that between the solute's surface charges and the electric field of the phage. The model is evaluated for four different proteins and a DNA oligomer. Results indicate that alignment in charged nematic media is a function not only of the solute's shape, but also of its electric multipole moments of net charge, dipole, and quadrupole. The relative importance of these terms varies greatly from one macromolecule to another, and evaluation of the experimental data indicates that these terms scale differently with ionic strength. For several of the proteins, the calculated alignment is sensitive to the precise position of the charged groups on the protein surface. This suggests that NMR alignment measurements can potentially be used to probe protein electrostatics. Inclusion of electrostatic interactions in addition to steric effects makes the extended model applicable to all liquid crystals used in biological NMR to date.