The back and forth of energy transfer between carotenoids and chlorophylls and its role in the regulation of light harvesting

Many aspects in the regulation of photosynthetic light-harvesting of plants are still quite poorly understood. For example, it is still a matter of debate which physical mechanism(s) results in the regulation and dissipation of excess energy in high light. Many researchers agree that electronic interactions between chlorophylls (Chl) and certain states of carotenoids are involved in these mechanisms. However, in particular, the role of the first excited state of carotenoids (Car S₁) is not easily revealed, because of its optical forbidden character. The use of two-photon excitation is an elegant approach to address directly this state and to investigate the energy transfer in the direction Car S₁ → Chl. Meanwhile, it has been applied to a large variety of systems starting from simple carotenoid–tetrapyrrole model compounds up to entire plants. Here, we present a systematic summary of the observations obtained by two-photon excitation about Car S₁ → Chl energy transfer in systems with increasing complexity and the correlation to fluorescence quenching. We compare these observations directly with the energy transfer in the opposite direction, Chl → Car S₁, for the same systems as obtained in pump-probe studies. We discuss what surprising aspects of this comparison led us to the suggestion that quenching excitonic Car–Chl interactions could contribute to the regulation of light harvesting, and how this suggestion can be connected to other models proposed.     

Phase-transition properties of glycerol-monopalmitate lipid bilayers investigated by molecular dynamics simulation: influence of the system size and force-field parameters

Phase-transition properties of glycerol-1-monopalmitate (GMP) bilayers are investigated using explicit-solvent molecular dynamics (MD) simulations, initiated from structures appropriate for the gel (GL) or liquid crystal (LC) phases, and carried out at different hydration levels and temperatures. Building up on a previous study and based on 600 ns simulations, the influence of the system size and of the force field on the equilibrium thermodynamic and dynamic parameters of the bilayers in the GL and LC phases, as well as on the temperature T_m and properties of the GL ? LC phase transition, are analysed. Qualitatively speaking, the results agree with the available experimental data for the area per lipid in the two phases and for the phase-transition temperatures at the three hydration levels irrespective of the selected model parameters. They also suggest that the total number of hydrogen bonds formed between a lipid headgroup and its environment is essentially constant, amounting to about four in both the LC and the GL phases. Quantitatively speaking, the dependence of T_m on the hydration level is found to be non-systematic across the different combinations of model parameters. This results in part from a sensitivity of the results on the system size and force-field parameters but also from the limited accuracy of the bracketing approach employed here to estimate T_m. Finally, a simple kinetic model is proposed to account for the timescales of the transitions. This model involves enthalpy and entropy increases of about 26 kJ mol^? 1 and 83 J mol^? 1 K^? 1 per lipid, upon going from the GL to the LC phase. The transition state is associated with activation parameters corresponding to 13% and 11%, respectively, of these values along the GL → LC transition, resulting in an activation free energy of about 0.3 kJ mol^? 1 per lipid at T_m.     

Molecular dynamics and quantum chemical studies of solvent effects on cyclo glycylglycine and glycylalanine dipeptides

Hydrogen bond (H-bond) interactions between the two cyclo dipeptides, cyclo(glycyl-glycine) (CGG) and cyclo(glycyl-alanine) (CGA), and water have been studied using molecular dynamics (MD) and quantum chemical methods. The MD studies have been carried out on CGG and CGA in water using fixed charge force field (AMBER ff03) for over 10 ns with a MD time step of 2 fs. The results of this study show that the solvation pattern influences the conformations of the cyclo dipeptides. Following molecular simulations, post Hartree–Fock and density functional theory methods have been used to explore the molecular properties of the cyclo dipeptides in gaseous and aqueous phase environments. The self-consistent reaction field theory has been used to optimise the cyclopeptides in diethyl ether (? = 4.3) and water (? = 78.5), and the solvent effects have been analysed. A cluster of eight water molecules leads to the formation of first solvation shell of CGG and CGA and the strong H-bonding mainly contributes to the interaction energies. The H-bond interactions have been analysed by the calculation of electron density ρ(r) and its Laplacian ▽²ρ(r) at bond critical points using atoms in molecules theory. The natural bond orbital analysis was carried out to reveal the nature of H-bond interactions. In the solvated complexes, the keto carbons registered the maximum NMR chemical shifts.     

Density functional computations of Rh(I)-catalysed hydroacylation and hydrogenation of ethene using formic acid

The rhodium-catalysed hydroacylation of alkene is one of the most useful C–H bond activation processes. The C–C bond-forming reactions via C–H bond activation have extensively been the focus of study in the fields of organic and organometallic chemistry. In this work, density functional theory has been used to study Rh(I)-catalysed hydroacylation and hydrogenation of ethene with formic acid. All the intermediates and the transition states were optimised completely at the B3LYP/6-311++G(d,p) level (LANL2DZ(d) for Rh, P). Calculation results confirm that Rh(I)-catalysed hydroacylation of ethene is exothermic and the released Gibbs free energy is ? 60.39 kJ/mol. Rh(I)-catalysed hydrogenation of ethene is also exothermic and the released Gibbs free energy is ? 150.97 kJ/mol. Rh(I)-catalysed hydroacylation of ethene is the dominant reaction mode for Rh(I)-catalysed hydroacylation and hydrogenation of ethene with formic acid. In Rh(I)-catalysed hydroacylation of ethene, the H-transfer reaction is prior to the C–C bond-forming reaction. Therefore, the reaction mode ‘a’ (i.e. ca → M1 → TS1 → M2 → TS2a → M3a → TS3a → M4 → P1) is the dominant reaction pathway for Rh(I)-catalysed hydroacylation and hydrogenation of ethene. The theoretically predicted dominant product is propane acid.     

Structural phase transition of BeTe: an ab initio molecular dynamics study

Beryllium telluride (BeTe) with cubic zinc-blende (ZB) structure was studied using ab initio constant pressure method under high pressure. The ab initio molecular dynamics (MD) approach for constant pressure was studied and it was found that the first order phase transition occurs from the ZB structure to the nickel arsenide (NiAs) structure. It has been shown that the MD simulation predicts the transition pressure P _T more than the value obtained by the static enthalpy and experimental data. The structural pathway reveals MD simulation such as cubic → tetragonal → orthorhombic → monoclinic → orthorhombic → hexagonal, leading the ZB to NiAs phase. The phase transformation is accompanied by a 10% volume drop and at 80 GPa is likely to be around 35 GPa in the experiment. In the present study, our obtained values can be compared with the experimental and theoretical results.

Open image in new window

Graphical abstract The energy-volume relation and ZB phase for the BeTe

     

Biosynthesis of the carbohydrate moieties of arabinogalactan proteins by membrane-bound β-glucuronosyltransferases from radish primary roots

A membrane fraction from etiolated 6-day-old primary radish roots (Raphanus sativus L. var hortensis) contained β-glucuronosyltransferases (GlcATs) involved in the synthesis of the carbohydrate moieties of arabinogalactan proteins (AGPs). The GlcATs transferred [¹⁴C]GlcA from UDP-[¹⁴C]GlcA on to β-(1 → 3)-galactan as an exogenous acceptor substrate, giving a specific activity of 50–150 pmol min^?1 (mg protein)^?1. The enzyme specimen also catalyzed the transfer of [¹⁴C]GlcA on to an enzymatically modified AGP from mature radish root. Analysis of the transfer products revealed that the transfer of [¹⁴C]GlcA occurred preferentially on to consecutive (1 → 3)-linked β-Gal chains as well as single branched β-(1 → 6)-Gal residues through β-(1 → 6) linkages, producing branched acidic side chains. The enzymes also transferred [¹⁴C]GlcA residues on to several oligosaccharides, such as β-(1 → 6)- and β-(1 → 3)-galactotrioses. A trisaccharide, α-l-Araf-(1 → 3)-β-Gal-(1 → 6)-Gal, was a good acceptor, yielding a branched tetrasaccharide, α-l-Araf-(1 → 3)[β-GlcA-(1 → 6)]-β-Gal-(1 → 6)-Gal. We report the first in vitro assay system for β-GlcATs involved in the AG synthesis as a step toward full characterization and cloning.     

Density functional studies on copper-catalysed asymmetric aziridination of diazoacetate with imines

Density functional theory has been used to study copper(I)-catalysed aziridination of diazoacetate with imines. All the intermediates and the transition states were optimised completely at B3LYP/6-31G(d) level. Calculation results confirm that copper(I)-catalysed aziridination of diazoacetate with imines is exothermic, and the total released Gibbs free energy is about ? 170 kJ/mol. Copper(I)-catalysed aziridination has two reaction modes: I and II, and thus the reaction mode I is dominant. The formation of the copper(I)–carbene–imine complex M3 (i.e. the attack of imines on copper–carbon(carbene) of copper–carbene intermediate M2) is the rate-determining step and the chirality-limiting step for copper-catalysed asymmetric aziridination. The reaction channel CA2 → M1a → TS1a → M2 → TS2a2 → M3a2 → TS3a2 → M4a2 → P1 is the most favourable one. The dominant products predicted theoretically are of (R)-chirality.     

Anion–π interactions: CP-corrected vs. standard optimisation,AIM and NBO analyses

The anion–π interactions between Br^? , Cl^? , F^?  and H^?  anions and hexafluorobenzene (HFB), 1,2,4,5-tetracyanobenzene (TCB) and tetracyanopyrazine (TCP) have been studied by standard and counterpoise (CP) corrected methods at HF, B3LYP and MP2/6-31+ + G (d,p) levels of theory. The complexation energies were corrected for basis set superposition error (ΔE _BSSE) and zero point energy (ΔE _BSSE + ZPE). Also, the B3LYP results were corrected by single-point calculation at B3LYP/aug-cc-PVTZ level of theory. Although the CP-corrected method results in higher distances between anions and rings, the standard method gives lower complexation energies. TCP…X^?  series gives higher complexation energies in both CP-corrected and standard methods. Topological analysis of the charge density ρ(r) has been performed by the means of atoms in molecules method on the wave functions obtained at MP2/6-31+ + G (d,p) level of theory. The number and the nature of critical points depend on aromatic ring and anion. Natural bond orbital analysis indicates that n_X → π*_CC and n_X → π*_CN interactions are the most important interactions for TCB (and HFB)…X^?  and TCP…X^?  complexes, respectively.     

Glycosaminoglycans from fish swim bladder: isolation,structural characterization and bioactive potential

The swim bladder of fish is an internal gas-filled organ that allows fish to control their buoyancy and swimming depth. Fish maws (the dried swim bladders of fish) have been used over many centuries as traditional medicines, tonics and a luxurious gourmet food in China and Southeast Asia. Little is known about the structural information of polysaccharides comprising this important functional material of fish tissue. In the present study, the total glycosaminoglycan (GAG) from fish maw was characterized. Two GAGs were identified, chondroitin sulfate (CS, having a molecular weight of 18–40 kDa) and heparan sulfate (HS), corresponding to 95% and 5% of the total GAG, respectively. Chondroitinase digestion showed that the major CS GAG was composed of ΔUA-1 → 3-GalNAc4S (59.7%), ΔUA-1 → 3-GalNAc4,6S (36.5%), ΔUA-1 → 3-GalNAc6S (2.2%) and ΔUA-1 → 3-GalNAc (1.6%) disaccharide units. ¹H–NMR analysis and degradation with specific chondroitinases, both CS-type A/C and CS-type B were present in a ratio of 1.4:1. Analysis using surface plasmon resonance showed that fibroblast growth factor (FGF)-2 bound to the CS fraction (K_D = 136 nM). These results suggest that this CS may be involved in FGF-signal pathway, mediating tissue repair, regeneration and wound healing. The CS, as the major GAG in fish maw, may have potential pharmacological activity in accelerating wound healing.     

Impacts of terminal (4R)-fluoroproline and (4S)-fluoroproline residues on polyproline conformation

Many interests have been focused on prolyl cis–trans isomerization which is related to protein folding and isomer-specific biochemical recognition. Since polyproline can adopt either type I (PPI) helices with all cis amide bonds or type II (PPII) helices with all trans amide bonds, it has been a valuable model to study the prolyl isomerization. Recent studies have shown that stereoelectronic effects govern the stability of PPII structure and the rate of PPII → PPI conversion. To further explore the terminal stereoelectronic effects on polyproline conformation, herein we synthesized a series of host–guest peptides in which (2S,4S)-4-fluoroproline (flp) or (2S,4R)-4-fluoroproline (Flp) residues are incorporated into the C- or N-terminal end of a peptide and studied the thermodynamic and kinetic consequences on polyproline conformation. Circular dichroism measurements revealed that inserting 4-fluoroproline residues into the C terminus of a polyproline peptide induces a great stereoelectronic effect on PPII stability and PPII → PPI conversion rates. From the C terminus, a (Flp)₃ triplet stabilizes PPII structure and increases the transition barrier of PPII → PPI conversion by 1.53 kJ mol^?1 while a (flp)₃ triplet destabilizes PPII conformation and reduce the PPII → PPI transition barrier by 4.61 kJ mol^?1. In contrast, the 4-fluoroproline substitutions at the N terminus do not exhibit distinct stereoelectronic effects on PPII stability and PPII → PPI conversion rates. Our data demonstrate that the C-terminal stereoelectronic effects have a more dramatic impact on PPII stability and PPII → PPI conversion kinetics.     

Nanomilling mechanism on Cu surfaces investigated using atomistic simulation

The nanoscale milling and scratching processes of copper workpieces are studied using molecular dynamics simulations based on the tight-binding and Morse potentials. The effects of the rotation velocity of the tool and the workpiece temperature are evaluated in terms of atomic trajectories, slip vectors, flow field of chips, cutting forces and groove characteristics. The simulation shows that a slip system in the ?110? direction on the workpiece surface occurs for milling with a tool rotation velocity of ω = 0°/fs. However, no apparent slip system appears for ω = 0.005°/fs or higher; instead, the number of amorphous areas increases. At ω = 0°/fs (nanoscratching), most of the removed atoms pile up in front of the tool and some gradually backfill when the tool rotates due to the effects of rotational friction and adhesion between the tool and the removed atoms. The largest number of removed atoms that piled up in front of the tool were obtained for milling with ω = 0°/fs; the number of removed atoms that piled up in front of the tool decreased with the increasing ω value. The component forces corresponding to the feed direction of the tool are the largest for the nanodrilling and nanomilling processes. High-precision grooves can be obtained at a low workpiece temperature (e.g. room temperature) with ω = 0°/fs.     

The effects of codon usage on the formation of secondary structures of nucleocapsid protein of peste des petits ruminants virus

The nucleocapsid (N) protein of peste des petits ruminants virus (PPRV) with a conserved amino acid usage pattern plays an important role in viral replication. The primary objective of this study was to estimate roles of synonymous codon usages of PPRV N gene and tRNA abundances of host in the formation of secondary structure of N protein. The potential effects of synonymous codon usages of N gene and tRNA abundances of host on shaping different folding units (α-helix, β-strand and the coil) in N protein were estimated, based on the information about the modeling secondary structure of PPRV N protein. The synonymous codon usage bias was found in different folding units in PPRV N protein. To better understand the role of translation speed caused by variant tRNA abundances in shaping the specific folding unit in N protein, we modeled the changing trends of tRNA abundance at the transition boundaries from one folding unit to another folding unit (β-strand → coil, coil → β-strand, α-helix → coil, coil → α-helix). The obvious fluctuations of tRNA abundance were identified at the two transition boundaries (β-strand → coil and coil → β-strand) in PPRV N protein. Our findings suggested that viral synonymous codon usage bias and cellular tRNA abundance variation might have potential effects on the formation of secondary structure of PPRV N protein.     

GPU-accelerated replica exchange molecular simulation on solid–liquid phase transition study of Lennard-Jones fluids

Determining the solid–liquid phase transition point by conventional molecular dynamics (MD) simulations is difficult because of the tendency of the system to get trapped in local minimum energy states at low temperatures and hysteresis during cooling and heating cycles. The replica exchange method, used in performing many MD simulations of the system at different temperature conditions simultaneously and performs exchanges of these temperatures at certain intervals, has been introduced as a tool to overcome this local-minimum problem. However, around the phase transition temperature, a greater number of different temperatures are required to adequately find the phase transition point. In addition, the number of different temperature values increases when treating larger systems resulting in huge computation times. We propose a computational acceleration of the replica exchange MD simulation on graphics processing units (GPUs) in studying first-order solid–liquid phase transitions of Lennard-Jones (LJ) fluids. The phase transition temperature for a 108-atom LJ fluid has been calculated to validate our new code. The result corresponds with that of a previous study using multicanonical ensemble. The computational speed is measured for various GPU-cluster sizes. A peak performance of 196.3 GFlops with one GPU and 8.13 TFlops with 64 GPUs is achieved.     

Mesoscopic simulation of self-assembly of linear and dendritic copolymer poly(styrene)-b-poly(ethyleneglycol) in polar and non-polar solvents

In this work, we simulate the microphase separation of aqueous and non-aqueous solutions of diblock copolymer poly(styrene)-b-poly(ethyleneglycol) under different architectures (linear and linear–dendritic) by dissipative particle dynamics. The observed morphologies in water where poly(ethyleneglycol) (PEG) block is soluble are as follows: (1) at low concentrations spherical micelles, cylinders and bicontinuous structures are formed in dendritic structures and spheres, cylinders and perforated lamellas in linear structure. The architectures simulated at low–moderate concentrations show an evolution sphere → cylinder → bicontinuous or perforated lamellas as the concentration is increased. (2) At high concentrated solutions rich defect structures of the sponge type are formed. In a non-aqueous non-polar solution such as cyclohexane, which is a good solvent for the polystyrene block, the formation of well-defined aggregates at low concentrations is not observed; however, irregular structures are achieved in concentrated solutions. We compare these results with a polymeric chimera consisting of a mixture of linear poly(styrene) homopolymer and PEG homopolymer in the linear, G1 or G2 dendritic configurations. Our simulations are in agreement with the experimentally observed structures of these polymers.     

Genetic Diversity of the Melanocortin 4 Receptor (MC4R) Gene and its Association with Slaughter Traits in the Landes Goose

Melanocortin 4 receptor (MC4R) plays a crucial part in regulating feeding behavior in humans and rodents. We detected two single nucleotide polymorphisms (SNPs; c.108G → A and c.627C → T) in the goose MC4R gene and genotyped 94 Landes geese for association analysis with several carcass traits. Significant associations (P < 0.05) were obtained for c.108G → A with carcass weight, breast muscle percentage, and leg muscle percentage, and for c.627C → T with body weight, carcass weight, semi-eviscerated weight, and eviscerated weight. We re-constructed haplotypes based on the two SNPs and analyzed diplotypes in association with carcass traits, obtaining significant associations with several of the traits. These results suggest that polymorphisms in the MC4R gene could have effects on carcass traits in Landes geese. More study is required to confirm these results.     

The Critical Roles of Information and Nonequilibrium Thermodynamics in Evolution of Living Systems

Living cells are spatially bounded, low entropy systems that, although far from thermodynamic equilibrium, have persisted for billions of years. Schrödinger, Prigogine, and others explored the physical principles of living systems primarily in terms of the thermodynamics of order, energy, and entropy. This provided valuable insights, but not a comprehensive model. We propose the first principles of living systems must include: (1) Information dynamics, which permits conversion of energy to order through synthesis of specific and reproducible, structurally-ordered components; and (2) Nonequilibrium thermodynamics, which generate Darwinian forces that optimize the system. Living systems are fundamentally unstable because they exist far from thermodynamic equilibrium, but this apparently precarious state allows critical response that includes: (1) Feedback so that loss of order due to environmental perturbations generate information that initiates a corresponding response to restore baseline state. (2) Death due to a return to thermodynamic equilibrium to rapidly eliminate systems that cannot maintain order in local conditions. (3) Mitosis that rewards very successful systems, even when they attain order that is too high to be sustainable by environmental energy, by dividing so that each daughter cell has a much smaller energy requirement. Thus, nonequilibrium thermodynamics are ultimately responsible for Darwinian forces that optimize system dynamics, conferring robustness sufficient to allow continuous existence of living systems over billions of years.     

Assessment of the elastic properties of amorphous calcium silicates hydrates (I) and (II) structures by molecular dynamics simulation

Based on Hamid model of 11Å tobermorite, amorphous calcium silicates hydrates (or C-S-H) structures (Ca₄Si₆O₁₄(OH)₄?2H₂O as the C-S-H(I) and (CaO)_1.67(SiO₂)(H₂O)_1.75 as the C-S-H(II)) with the Ca/Si ratio of 0.67 and 1.7 are concerned. Then, as the representative ‘globule’ C-S-H, two amorphous C-S-H structures with the size of 5.352 × 4.434 × 4.556 nm³ during the stretch process are simulated at a certain strain rate of 10^?3 ps^?1 by LAMMPS program for molecular dynamics simulation, using ClayFF force field. The tensile stress–strain curves are obtained and analysed. Besides, elastic modulus of the ‘globule’ C-S-H is calculated to assess the elastic modulus of C-S-H phases (the low-density C-S-H – LD C-S-H – and the high-density C-S-H – HD C-S-H), where the porosity is a critical factor for explaining the relationship between ‘globule’ C-S-H at nanoscale and C-S-H phases at microscale. Results show that: (1) The C-S-H(I) structure has transformed from crystalline to amorphous during the annealing process, Young’s moduli in x, y and z directions are almost the same. Besides, the extent of aggregation and aggregation path for water molecules in the structure is different in three directions. (2) Young’s modulus of both amorphous C-S-H(I) and C-S-H(II) structures with a size of about 5 nm under strain rate of 10^?3 ps^?1 at 300 K in three directions is averaged to be equal, of which C-S-H(II) structure is about 60.95 GPa thus can be seen as the elastic modulus of the ‘globule’ C-S-H. (3) Based on the ‘globule’ C-S-H, the LD C-S-H and HD C-S-H can be assessed by using the Self-Consistent Scheme (separately 18.11 and 31.45 GPa) and using the Mori–Tanaka scheme (29.78 and 37.71 GPa), which are close to the nanoindentation experiments by Constantinides et al. (21.7 and 29.4 GPa).     

Investigating lithium ion conductivity in amorphous solids containing [V10O28]6− clusters by computer simulation

Computer simulation method was applied to investigate the migration of lithium ion in three amorphous solid systems containing polyoxovanadate (POV) clusters [V₁₀O₂₈]^6 ? . The cluster was adopted from a recently synthesized crystalline poly[octa-μ-aqua-octaaqua-μ-decavanadato-hexalithium] (POAODH). The simulated POV systems correspond to amorphous solid half-dehydrated solid and completely dehydrated solid doped with LiCl salt. The simulation results show large diffusion constants of lithium ions in all systems in spite of highly negatively charged [V₁₀O₂₈]^6 ?  clusters presented in the system. The estimated ionic conductivity due to the migration of lithium ions reaches a magnitude of 10^? 4 S/m. The conductivity increases as the water content in the system decreases. The analysis of moving trajectories shows the lithium ion moves around the oxygen sites of POV clusters and hops between them. The estimated displacement of lithium ion is about 4～5 Å, which is much larger than the corresponding displacement of lithium ion in a polymer matrix. Rapidly rotating clusters shown by orientation correlation function analysis, in conjunction with the large separation between clusters in the system, provides favorable conditions for the large amplitude migration of lithium ions.     

Relationships Between Single Nucleotide Polymorphisms of the H-FABP Gene and Slaughter and Meat Quality Traits in Chicken

Using PCR-SSCP with five primer pairs, we detected six single nucleotide polymorphisms of the H-FABP gene: 332G → A, 534G → A, 783C → T, 835C → T, 1198T → C, and 2329C → T. Chi-square results showed significant differences (P < 0.05) in genotype frequency among breeds in Fragment 1 and extremely significant differences (P < 0.01) in Fragments 2–4. We found a significant association between Fragment 2 genotype and muscle fiber number, Arg and Thr (P < 0.05); between Fragment 3 genotype and living weight, carcass weight, breast muscle weight, abdominal fat weight, and abdominal fat percentage (P < 0.05); between Fragment 4 genotype and Thr, Phe, and inosinic acid (P < 0.05). It was concluded that H-FABP was the major gene influencing slaughter performance and meat quality or was linked with the major gene in these strains and that the C783T mutation could be used as a candidate molecular genetic marker for breeding selection. The combination M1C2–B2B2–D1D1 is an ideal model for breeding in these strains because it can improve slaughter and meat quality traits.