Fundamental symmetry aspects of chirality

Physical systems which exhibit distinguishable enantiomers under space inversion are not necessarily chiral. A new definition of chirality is proposed that enables true and false chirality to be distinguished. Although spatial enantiomorphism is sufficient to guarantee chirality in a stationary object, enantiomorphous systems are not necessarily chiral when motion is involved. Only a truly chiral influence can induce absolute asymmetric synthesis in a reaction mixture at thermodynamic equilibrium, but false chirality might suffice if equilibrium is not attained. Parity violation lifts only the degeneracy of enantiomers of truly chiral systems, the true enantiomers (i.e. strictly degenerate) being interconverted by space inversion together with charge conjugation. The time-independence of optical activity arising from parity violation is contrasted with the time-dependence of that arising from spontaneous parity breaking.     

Antropologic factors in prediction of ovary cancer

The aim of the study was to determine a combination of anthropometric variables that would enable better differentiation between benign and malignant ovarian masses. Prospective study has been performed in a two year period in which 208 women with ovarian lesions were analyzed and correlated with histopathologic surgical findings. We examined the relation between self-reported anthropometric and other variables (height, weight, body mass index--BMI, parity, marital status, education, age, rural versus urban residence, menopausal status) and incidence of ovarian cancer. Age, parity, marital status and menopausal status individually showed statistical significance.     

Calculation of the stability of mini-duplexes of DNA in an electrolyte solution using a molecular mechanics method in approximating a statistical model of the environment

Approximation of the statistical model of environment (SME) to estimate the energy of the macromolecule in electrolyte solution has been developed and used for calculating the conformational energy of nucleic acids by means of the molecular mechanics method. Calculation of base pairs opening delta Hcalop enthalpies and enthalpies of DNA miniduplexes dissociation delta Hcaldis were performed for 10 types of diduplexes. The approximation SME enable to perform calculations of the absolute base-dependent values of delta Hcalop which coincide in the range of 1 kcal/mol with the experimental base-dependent values of the helix-coil transition enthalpies. Values of dissociation enthalpies delta Hcalop greater than delta Hcaldis for all miniduplexes, the difference of delta Hcalop--delta Hcaldis determine the base-dependent energy of the helix-coil boundary. The values of activation barriers for the strands dissociation delta H d not equal to congruent to 8 kcal/mol and association delta H not equal to as congruent to 4 kcal/mol were obtained for GG/CC and AA/TT duplexes. It is concluded that the approximation SME enables to increase substantially the accuracy of the calculations of the macromolecule conformational rearrangement enthalpies in the electrolytes solution.     

Parity violation and the evolution of biomolecular homochirality

Parity violation at the level of terrestrial biopolymers, as seen in proteins, DNAs, and RNAs, and parity violation at the level of nuclear processes, as evident in longitudinally polarized beta-particles and parity-violating energy differences (PVEDs), are discussed and their fundamental importances are emphasized. Attempts to find a causal connection between the unique homochirality of biopolymers and parity violation at the nuclear level, and speculations that the former is a consequence of the latter, are reviewed. Consideration of all lines of evidence leads to the conclusion that there is no substantiation for such a causal connection, and that the two levels of parity violation are entirely independent of each other.     

The origin of chirality in protein amino acids

We discuss the origin of the chirality of protein amino acids from the point of view of a phase transition from a racemic mixture into an optically pure state. We assume that Bose–Einstein condensation may act as an amplification mechanism. The original theory is due to Salam. We suggest a new role for the phase transition. Following Quack we distinguish parity violation of two kinds (de facto and de lege symmetry breaking). While the Salam phase transition corresponds to parity violation of the second kind (de lege), the phase transition we discuss in this work corresponds to parity violation of what we may call a third kind. This is suggested by recent experimental phenomena which correlate chiral symmetry breaking and pattern formation (spontaneous symmetry breaking that separates an initial racemic mixture into right- and left-handed space domains by means of a substrate). Tentative comments are given on the eventual design of possible experiments that may test this new hypothesis. © 1994 Wiley-Liss, Inc.     

An oxorhenium complex bearing a chiral cyclohexane‐1‐olato‐2‐thiolato ligand: Synthesis,stereochemistry, and theoretical study of parity violation vibrational frequency shifts

In our effort towards measuring the parity violation energy difference between two enantiomers, a simple chiral oxorhenium complex 5 bearing enantiopure 2‐mercaptocyclohexan‐1‐ol has been prepared as a potential candidate species. Vibrational circular dichroism revealed a chiral environment surrounding the rhenium atom, even though the rhenium is not a stereogenic center itself, and enabled to assign the (1S,2S)‐(?) and (1R,2R)‐(+) absolute configuration for 5 . For both compound 5 and complex 4 , previously studied by us and bearing a propane‐2‐olato‐3‐thiolato ligand, relativistic calculations predict parity violating vibrational frequency differences of a few hundreds of millihertz, above the expected sensitivity attainable by a molecular beam Ramsey interferometer that we are constructing.     

Contribution of the trifluoroacetyl group in the thermodynamics of antigen–antibody binding

We analyzed the binding of the 7C8 antibody to the chloramphenicol phosphonate antigens—one containing a trifluoroacetyl group (CP‐F) and the other containing an acetyl group (CP‐H)—by using isothermal titration calorimetry (ITC). The thermodynamic difference due to the substitution of F by H was evaluated using free energy calculations based on molecular dynamics (MD) simulations. We have previously shown that another antibody, namely, 6D9, binds more weakly to CP‐H than to CP‐F, mainly due to the different hydration free energies of the dissociated state and not due to the unfavorable hydrophobic interactions with the antibody in the bound state. Unlike in the binding of the trifluoroacetyl group with 6D9, in its binding with 7C8, it is exposed to the solvent, as seen in the crystal structure of the complex of 7C8 with CP‐F. The thermodynamic analysis performed in this study showed that the binding affinity of 7C8 for CP‐H is similar to that for CP‐F, but this binding to CP‐H is accompanied with less favorable enthalpy and more favorable entropy changes. The free energy calculations indicated that, upon the substitution of F by H, enthalpy and entropy changes in the associated and dissociated states were decreased, but the magnitude of enthalpy and entropy changes in the dissociated state was larger than that in the associated state. The differences in binding free energy, enthalpy, and entropy changes determined by the free energy calculations for the substitution of F by H are in good agreement with the experimental results. Copyright © 2009 John Wiley & Sons, Ltd.     

Thermodynamic descriptors, profiles and driving forces in membrane receptor-ligand interactions

Extension of the (isothermal) Gibbs-Helmholtz equation for the heat capacity terms (ΔC(p)) allows formulating a temperature function of the free (Gibbs) energy change (ΔG). An approximation of the virtually unknown ΔC(p) temperature function enables then to determine and numerically solve temperature functions of thermodynamic parameters ΔH and ΔS (enthalpy and entropy change, respectively). Analytical solutions and respective numeric procedures for several such approximation formulas are suggested in the presented paper. Agreement between results obtained by this analysis with direct microcalorimetric measurements of ΔH (and ΔC(p) derived from them) was approved on selected cases of biochemical interactions presented in the literature. Analysis of several ligand-membrane receptor systems indicates that temperature profiles of ΔH and ΔS are parallel, largely not monotonic, and frequently attain both positive and negative values within the current temperature range of biochemical reactions. Their course is determined by the reaction change of heat capacity: temperature extremes (maximum or minimum) of both ΔH and ΔS occur at ΔC(p)=0, for most of these systems at roughly 285-305 K. Thus, the driving forces of these interactions may change from enthalpy-, entropy-, or enthalpy-entropy-driven in a narrow temperature interval. In contrast, thermodynamic parameters of ligand-macromolecule interactions in solutions (not bound to a membrane) mostly display a monotonic course. In the case of membrane receptors, thermodynamic discrimination between pharmacologically defined groups-agonists, partial agonists, antagonists-is in general not specified and can be achieved, in the best, solely within single receptor groups.     

Beta-gamma decay and stereoselective molecular breaking

Summary Firm empirical evidence demonstrating that optical activity in biological molecules could result from the violation of parity in the weak interactions has yet to be produced. This effect, when mediated by beta particles or by their subsequent circularly polarized external bremsstrahlung, seems to exist, but it is too small for experimental verification. In this paper, I suggest another mechanism to transfer the dissymmetry from the nucleus to the molecule: a compound beta (or electron capture)-gamma disintegration process which, carrying also the information of parity violation, might be more efficient at the time of a selective breaking of one of the two enantiomers in a racemic mixture. A new type of experiment to test this idea is suggested.     

Thermodynamic properties of two-dimensional nonideal structures with isotropic pair potential

An approximation is proposed for energy density in two-dimensional nonideal systems for a wide class of isotropic repulsive pair interparticle potentials. The approximation allows one to determine the main thermodynamic functions and characteristics of the system by using well-known thermodynamic formulas. The results obtained with the help of this approximation are compared with the data of numerical simulations of thermodynamic properties of the structures under study. The simulations were performed in a wide range of parameters typical of laboratory dusty plasmas. Main attention was paid to the screened Coulomb potential.     

Gas-chromatographic separation of tri(hetero)halogenomethane enantiomers

Five-atomic tri(hetero)halogenomethanes represent the simplest class of non-isotopic small chiral molecules suitable for the study of fundamental aspects of chirality. The analytical gas-chromatographic separation of the enantiomers of bromochlorofluoromethane 1 and of chlorofluoroiodomethane 2 on the immobilized chiral stationary phase octakis(3-O-butanoyl-2,6-di-O-n-pentyl)-gamma-cyclodextrin 3, chemically linked to polydimethylsiloxane, is described. By temperature-dependent thermodynamic measurements very low isoenantioselective temperatures T(iso) are found and for optimum enantiomeric separations cryogenic temperatures are required. The ee values of enantiomerically enriched tri(hetero)halogenomethanes 1 and 2 are determined and relative configurations are correlated with the chromatographic elution order of 1 and 2 on 3.     

Effect of competitive terminal electron acceptor processes on dechlorination of cis-1,2-dichloroethene and 1,2-dichloroethane in constructed wetland soils

Thermodynamic calculations were coupled with time-series measurements of chemical species (parent and daughter chlorinated solvents, H(2), sulfite, sulfate and methane) to predict the anaerobic transformation of cis-1,2-dichloroethene (cis-1,2-DCE) and 1,2-dichloroethane (1,2-DCA) in constructed wetland soil microcosms inoculated with a dehalorespiring culture. For cis-1,2-DCE, dechlorination occurred simultaneously with sulfite and sulfate reduction but competitive exclusion of methanogenesis was observed due to the rapid H(2) drawdown by the dehalorespiring bacteria. Rates of cis-1,2-DCE dechlorination decreased proportionally to the free energy yield of the competing electron acceptor and proportionally to the rate of H(2) drawdown, suggesting that H(2) competition between dehalorespirers and other populations was occurring, affecting the dechlorination rate. For 1,2-DCA, dechlorination occurred simultaneously with methanogenesis and sulfate reduction but occurred only after sulfite was completely depleted. Rates of 1,2-DCA dechlorination were unaffected by the presence of competing electron-accepting processes. The absence of a low H(2) threshold suggests that 1,2-DCA dechlorination is a cometabolic transformation, occurring at a higher H(2) threshold, despite the high free energy yields available for dehalorespiration of 1,2-DCA. We demonstrate the utility of kinetic and thermodynamic calculations to understand the complex, H(2)-utilizing reactions occurring in the wetland bed and their effect on rates of dechlorination of priority pollutants.     

Thermodynamic efficiency of bacterial growth calculated from growth yield of Pseudomonas oxalaticus OX1 in the chemostat

In order to determine the thermodynamic efficiency of bacterial growth, Pseudomonas oxalaticus OX1 was grown in carbon-limited continuous cultures. 11 different carbon sources, ranging from oxalate (most oxidised component) to ethanol (most reduced component), were used as limiting substrate in these experiments. From the experimental yield values (expressed as C-mol dry weight produced per C-mol carbon substrate consumed) the thermodynamic efficiencies were calculated. On substrates more reduced than biomass (such as ethanol and glycerol) the thermodynamic efficiency of growth of P. oxalaticus was negative but it reached a maximum of 23 +/- 3% with substrates with a degree of reduction of 3 (citrate) and lower. The actual concentrations of the components involved were incorporated into the calculations but this affected the overall thermodynamic efficiency only to a small extent. This result strengthens the conclusion of Westerhoff et al. (Westerhoff, H.V., Hellingwerf, K.J. and Van Dam, K. (1983) Proc. Natl. Acad. Sci. 80, 305-309) that bacteria have been optimised towards a theoretical thermodynamic efficiency of 24%, corresponding with maximisation of growth rate at optimal efficiency, with highly oxidised substrates.     

Quantitative determination of the main glucose metabolic fluxes in human erythrocytes by 13C- and 1H-MR spectroscopy.

Information displayed by homonuclear and heteronuclear spin-coupling patterns in 13C- and 1H-MR spectra allowed us to identify the major lactate isotopomers produced either from [1-(13)C]-glucose or from [2-(13)C]-glucose by human erythrocytes. Relative concentrations of detectable isotopomers were determined by integrating the corresponding MR signals. The interpretation of these data in terms of the fractional glucose metabolised through glycolysis and pentose phosphate pathway was performed by a computer simulation of the metabolism that took into account metabolic schemes pertaining to glycolysis and to the F-type of pentose phosphate pathway. The simulation was organised in a way to anticipate the populations of the isotopomers produced from any precursor at a priori established metabolic steady state. By the simulation, isotopomer populations were determined according to different values of pentose cycle, defined as the flux of glyceraldehyde 3-phosphate originating from pentose phosphate pathway at unitary glucose uptake. The populations of the isotopomers originating from [2-(13)C]-glucose were described by polynomials, and ratios between the polynomials were used in conjunction with 13C- and 1H-MR data to determine pentose cycle values. The knowledge of glucose uptake and of pentose cycle value allowed us to perform accurate measurement of the pentose phosphate pathway flux, of the hexokinase and phosphofructokinase fluxes as well as, indirectly, of the carbon dioxide production.     

Cumulative bondomers: a new concept in flux analysis from 2D [13C,1H] COSY NMR data

A well-established way of determining metabolic fluxes is to measure 2D [(13)C,(1)H] COSY NMR spectra of components of biomass grown on uniformly (13)C-labeled carbon sources. When using the entire set of measured data to simultaneously determine all fluxes in a proposed metabolic network model, the (13)C-labeling distribution in all measured compounds has to be simulated. This requires very large sets of isotopomer or cumomer balances. This article introduces the new concept of bondomers; entities that only vary in the numbers and positions of C-C bonds that have remained intact since the medium substrate molecule entered the metabolism. Bondomers are shown to have many analogies to isotopomers. One of these is that bondomers can be transformed to cumulative bondomers, just like isotopomers can be transformed to cumomers. Similarly to cumomers, cumulative bondomers allow an analytical solution of the entire set of balances describing a metabolic network. The main difference is that cumulative bondomer models are considerably smaller than corresponding cumomer models. This saves computational time, allows easier identifiability analysis, and yields new insights in the information content of 2D [(13)C,(1)H] COSY NMR data. We illustrate the theoretical concepts by means of a realistic example of the glycolytic and pentose phosphate pathways. The combinations of 2D [(13)C,(1)H] COSY NMR data that allow identification of all metabolic fluxes in these pathways are analyzed, and it is found that the NMR data contain less information than was previously expected.     

Chiroptical signatures of life and fundamental physics

This paper aims to inspire experimentalists to carry out proposed new chiroptical experiments springing from the theoretical study of the role of parity violation in the origin of biomolecular homochirality and to provide a brief update on the current status of calculations of the electroweak parity-violating energy difference (PVED) between enantiomers. If the PVED did select life's handedness, we would expect to find life on other planets consistently using the same hand as terrestrial biochemistry. Much more importantly, even finding the "wrong" hand (rather than a racemic mixture) on another planet could be the homochiral signature of life, and we discuss our proposal for chiroptical detection of life on extra-solar planets. The PVED may also have an exciting future as a "molecular footprint" of fundamental physics: comparison of calculated PVEDs with measured values could one day allow chemists to do "table-top particle physics" more cheaply with improved chiroptical techniques instead of ever larger particle accelerators. We discuss our proposed chiroptical method to measure the PVED by using molecular beams. To our knowledge, optical rotation has not yet been measured in molecular beams, but the rewards of doing so include a host of other "first ever" results in addition to measurement of the PVED. Chirality 24:764-769, 2012. ? 2012 Wiley Periodicals, Inc.     

The heat capacity paradox of ligand binding proteins: reconciling the microscopic and macroscopic world

Differential scanning microcalorimetry (DSC) is a superb method for the analysis of protein energetics. However, the relative simplicity of application has led astray many to assume that a proper analysis of the data was possible without a sound knowledge of the underlying statistical thermodynamic principles. In this study, the question is addressed of how to calculate properly the heat capacity signal of a protein in the presence of high affinity ligands. It is shown that the signal corresponds neither to grand canonic nor to canonic heat capacity. Statistical thermodynamic model calculations result only in the observed macroscopic heat capacity signal, if the protein in the calorimetric cell is assumed to form a grand canonic ensemble (T, p, mu controlled) which is, however, heated under constraints typical for a canonic ensemble (T, p, N controlled). As a consequence, the microscopic statistical thermodynamic heat capacity must be carefully distinguished from the macroscopically observable thermodynamic heat capacity in those cases where proteins unfold in the presence of high affinity ligands.     

Alveolar progenitor cells develop in mouse mammary glands independent of pregnancy and lactation

We have previously described pluripotent, parity-induced mammary epithelial cells (PI-MEC) marked by Rosa26-lacZ expression in the mammary glands of parous females. PI-MEC act as lobule-limited epithelial stem/progenitor cells. To determine whether parity is necessary to generate PI-MEC, we incubated mammary explant cultures from virgin mice in vitro with insulin alone (I), hydrocortisone alone (H), prolactin alone (Prl), or a combination of these lactogenic hormones (IHPrl). Insulin alone activated the WAP-Cre gene. Hydrocortisone and prolactin alone did not. Any combination of hormones that included insulin was effective. Only I, H and Prl together were able to induce secretory differentiation and milk protein synthesis. In addition, EGF, IGF-2 and IGF-1 added individually produced activated (lacZ(+)) PI-MEC in explant cultures. Neither estrogen nor progesterone induced WAP-Cre expression in the explants. None of these positive initiators of WAP-Cre expression in PI-MEC were effective in mammospheres or two-dimensional cultures of mammary epithelium, indicating the indispensability of epithelial-stromal interaction in PI-MEC activation. Like PI-MEC, lacZ(+) cells from virgin explants proliferated and contributed progeny to mammospheres in vitro and to epithelial outgrowths in vivo after transplantation. LacZ(+) cells induced in virgin mouse mammary explants were multipotent (like PI-MEC) in impregnated hosts producing lacZ(+) mammary alveolar structures comprised of both myoepithelial and luminal progeny. These data demonstrate PI-MEC, a mammary epithelial sub-population of lobule-limited progenitor cells, are present in nulliparous female mice before parity and, like the PI-MEC observed following parity, are capable of proliferation, self-renewal and the capacity to produce progeny of diverse epithelial cell fates.     

Epidemic Models with Heterogeneous Mixing and Treatment

We consider a two-group epidemic model with treatment and establish a final size relation that gives the extent of the epidemic. This relation can be established with arbitrary mixing between the groups even though it may not be feasible to determine the reproduction number for the model. If the mixing of the two groups is proportionate, there is an explicit expression for the reproductive number and the final size relation is expressible in terms of the components of the reproduction number. We also extend the results to a two-group influenza model with proportionate mixing. Some numerical simulations suggest that (i) the assumption of no disease deaths is a good approximation if the disease death rate is small and (ii) a one-group model is a close approximation to a two-group model but a two-group model is necessary for comparing targeted treatment strategies. This research was supported by MITACS and an NSERC Research Grant.     

Can contact potentials reliably predict stability of proteins?

The simplest approximation of interaction potential between amino acid residues in proteins is the contact potential, which defines the effective free energy of a protein conformation by a set of amino acid contacts formed in this conformation. Finding a contact potential capable of predicting free energies of protein states across a variety of protein families will aid protein folding and engineering in silico on a computationally tractable time-scale. We test the ability of contact potentials to accurately and transferably (across various protein families) predict stability changes of proteins upon mutations. We develop a new methodology to determine the contact potentials in proteins from experimental measurements of changes in protein's thermodynamic stabilities (DeltaDeltaG) upon mutations. We apply our methodology to derive sets of contact interaction parameters for a hierarchy of interaction models including solvation and multi-body contact parameters. We test how well our models reproduce experimental measurements by statistical tests. We evaluate the maximum accuracy of predictions obtained by using contact potentials and the correlation between parameters derived from different data-sets of experimental (DeltaDeltaG) values. We argue that it is impossible to reach experimental accuracy and derive fully transferable contact parameters using the contact models of potentials. However, contact parameters may yield reliable predictions of DeltaDeltaG for datasets of mutations confined to the same amino acid positions in the sequence of a single protein.