Free energy calculations of glycosaminoglycan - Protein interactions

Glycosaminoglycans (GAGs) are complex highly charged linear polysaccharides that have a variety of roles in biological processes. We report the first use of molecular dynamics (MD) free energy calculations using the MM/PBSA method to investigate the binding of GAGs to protein molecules, namely the platelet endothelial cell adhesion molecule 1 (PECAM-1) and annexin A2. Calculations of the free energy of the binding of heparin fragments of different sizes reveal the existence of a region of low GAG-binding affinity in domains 5-6 of PECAM-1 and a region of high affinity in domains 2-3, consistent with experimental data and ligand-protein docking studies. A conformational hinge movement between domains 2 and 3 was observed, which allows the binding of heparin fragments of increasing size (pentasaccharides to octasaccharides) with an increasingly higher binding affinity. Similar simulations of the binding of a heparin fragment to annexin A2 reveal the optimization of electrostatic and hydrogen bonding interactions with the protein and protein-bound calcium ions. In general, these free energy calculations reveal that the binding of heparin to protein surfaces is dominated by strong electrostatic interactions for longer fragments, with equally important contributions from van der Waals interactions and vibrational entropy changes, against a large unfavorable desolvation penalty due to the high charge density of these molecules.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

Calculations of the free energy of interaction of the c-Fos−c-Jun coiled coil: Effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein−protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson−Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos−c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos−c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.

Theoretical calculations of base-base interactions in nucleic acids: II Stacking interactions in polynucleotides

Base-base interactions were computed for single- and double- stranded polynucleotides, for all possible base sequences. In each case, both right and left stacking arrangements are energetically possible. The preference of one over the other depends upon the base-sequence and the orientation of the bases with respect to helix-axis. Inverted stacking arrangement is also energetically possible for both single- and double-stranded polynucleotides. Finally, interaction energies of a regular duplex and the alternative structures3 were compared. It was found that the type II model3 is energetically more favourable than the rest.

Theoretical calculations of base-base interactions in nucleic acids: I Stacking interactions in free bases

Stacking interactions in free bases were computed on the basis of molecular association. The results of the calculations were compared with the stacking patterns observed in a few single crystals of nucleic acid components as examples. The following are the conclusions: (i) there can be two types of stacking pattern classified as normal and inverted types for any two interacting bases and both can be energetically favourable (ii) in both the types the stacking interaction is a combined effect of the overlap of the interacting bases and relative positions and orientations of the atomic centres of the two bases (iii) crystal symmetry and H-bonding interaction may influence stacking patterns.

Simplified theory of optical nonlinearities in spin-polarized semiconductors

The spin degree of freedom is largely disregarded in existing theories of the density-dependent optical properties of an interacting electron-hole plasma in quasiequilibrium. Here, we extended the pair equation, which is applicable to a bulk semiconductor at elevated temperatures, to calculate optical nonlinearities due to a spin-polarized plasma. We obtained agreement with recent circular dichroism data in laser-excited GaAs by using the plasma density alone as the fitting parameter. The simplicity of our theory, based on the analytical pair-equation formula, makes it ideal for conveniently modelling absorption in a carrier spin-polarized semiconductor.

Observation of the intense de-pairing effect in YBa2Cu3O7-delta due to the spin injection from La0.5Sr0.5CoO3

We report experimental evidence for a huge pair breaking effect induced by spin polarized quasiparticles in a YBa2Cu3O7-delta/La0.5Sr0.5CoO3 bi-layer fabricated by pulsed laser deposition. The temperature dependent magnetization measurements show evidence for the presence of both ferromagnetic and diamagnetic phases in the bi-layer. The current dependent electrical transport studies in the bi-layer exhibit a significant reduction in the superconducting transition temperature with the increase in applied current as compared to a single YBa2Cu3O7-delta layer and it follows a I-2/3 dependence in accordance with the pair breaking effect. Here, we find that the current driven from a ferromagnetic electrode with low spin polarization, such as La0.5Sr0.5CoO3 (-11%), into the superconductor can act as a strong pair breaker. This indicates that the spin polarization of the injecting electrode is not the only criterion in determining the pair breaking effect, rather the transparency of the interface for the spin polarization may also be significant. More interestingly, the spin diffusion length for YBa2Cu3O7-delta has a much longer length scale than that reported earlier in the study of ferromagnetic/superconducting heterostructures.

Spin-Probe ESR Studies on Nanocomposite Polymer Electrolytes

The possibility of using spin-probe electron spin resonance (ESR) as a tool to study glass transition temperature, T g, of polymer electrolytes is explored in 4 hydroxy 2,2,6,6 tetramethylpiperidine N oxyl (TEMPOL) doped composite polymer electrolyte (PEG)46LiClO4 dispersed with nanoparticles of hydrotalcite. The T g is estimated from the measured values of T 50G, the temperature at which the extrema separation 2A zz of the broad powder spectrum decreases to 50 G. In another method, the correlation time τc for the spin probe dynamics was determined by computer simulation of the ESR spectra and T g has been identified as the temperature at which τc begins to show temperature dependence. While both methods give values of T g close to those obtained from differential scanning calorimetry, it is concluded that more work is required to establish spin-probe ESR as a reliable technique for the determination of T g.

Vitrification, relaxation and free volume in glycerol-water binary liquid mixture:Spin probe ESR studies

Glass transition and relaxation of the glycerol-water (G-W) binary mixture system have been studied over the glycerol concentration range of 5-85 mol% by using the highly sensitive technique of electron spin resonance (ESR). For the water rich mixture the glass transition,sensed by the dissolved spin probe, arises from the vitrified mesoscopic portion of the binary system. The concentration dependence of the glass transition temperature manifests a closely related molecular level cooperativity in the system. A drastic change in the mesoscopic structure of the system at the critical concentration of 40 mol is confirmed by an estimation of the spin probe effective volume in a temperature range where the tracer reorientation is strongly coupled to the system dynamics.

Nonspecifically bound proteins spin while diffusing along DNA

It is known that DNA-binding proteins can slide along the DNA helix while searching for specific binding sites, but their path of motion remains obscure. Do these proteins undergo simple one-dimensional (1D) translational diffusion, or do they rotate to maintain a specific orientation with respect to the DNA helix? We measured 1D diffusion constants as a function of protein size while maintaining the DNA-protein interface. Using bootstrap analysis of single-molecule diffusion data, we compared the results to theoretical predictions for pure translational motion and rotation-coupled sliding along the DNA. The data indicate that DNA-binding proteins undergo rotation-coupled sliding along the DNA helix and can be described by a model of diffusion along the DNA helix on a rugged free-energy landscape. A similar analysis including the 1D diffusion constants of eight proteins of varying size shows that rotation-coupled sliding is a general phenomenon. The average free-energy barrier for sliding along the DNA was 1.1 +/- 0.2 k(B)T. Such small barriers facilitate rapid search for binding sites.

Analyzing changes in macro-level drivers of country-specific emissions : The role of consumption, technology and trade

Industrial ecology is an important field of sustainability science. It can be applied to study environmental problems in a policy relevant manner. Industrial ecology uses ecosystem analogy; it aims at closing the loop of materials and substances and at the same time reducing resource consumption and environmental emissions. Emissions from human activities are related to human interference in material cycles. Carbon (C), nitrogen (N) and phosphorus (P) are essential elements for all living organisms, but in excess have negative environmental impacts, such as climate change (CO2, CH4 N2O), acidification (NOx) and eutrophication (N, P). Several indirect macro-level drivers affect emissions change. Population and affluence (GDP/capita) often act as upward drivers for emissions. Technology, as emissions per service used, and consumption, as economic intensity of use, may act as drivers resulting in a reduction in emissions. In addition, the development of country-specific emissions is affected by international trade. The aim of this study was to analyse changes in emissions as affected by macro-level drivers in different European case studies. ImPACT decomposition analysis (IPAT identity) was applied as a method in papers I III. The macro-level perspective was applied to evaluate CO2 emission reduction targets (paper II) and the sharing of greenhouse gas emission reduction targets (paper IV) in the European Union (EU27) up to the year 2020. Data for the study were mainly gathered from official statistics. In all cases, the results were discussed from an environmental policy perspective. The development of nitrogen oxide (NOx) emissions was analysed in the Finnish energy sector during a long time period, 1950 2003 (paper I). Finnish emissions of NOx began to decrease in the 1980s as the progress in technology in terms of NOx/energy curbed the impact of the growth in affluence and population. Carbon dioxide (CO2) emissions related to energy use during 1993 2004 (paper II) were analysed by country and region within the European Union. Considering energy-based CO2 emissions in the European Union, dematerialization and decarbonisation did occur, but not sufficiently to offset population growth and the rapidly increasing affluence during 1993 2004. The development of nitrogen and phosphorus load from aquaculture in relation to salmonid consumption in Finland during 1980 2007 was examined, including international trade in the analysis (paper III). A regional environmental issue, eutrophication of the Baltic Sea, and a marginal, yet locally important source of nutrients was used as a case. Nutrient emissions from Finnish aquaculture decreased from the 1990s onwards: although population, affluence and salmonid consumption steadily increased, aquaculture technology improved and the relative share of imported salmonids increased. According to the sustainability challenge in industrial ecology, the environmental impact of the growing population size and affluence should be compensated by improvements in technology (emissions/service used) and with dematerialisation. In the studied cases, the emission intensity of energy production could be lowered for NOx by cleaning the exhaust gases. Reorganization of the structure of energy production as well as technological innovations will be essential in lowering the emissions of both CO2 and NOx. Regarding the intensity of energy use, making the combustion of fuels more efficient and reducing energy use are essential. In reducing nutrient emissions from Finnish aquaculture to the Baltic Sea (paper III) through technology, limits of biological and physical properties of cultured fish, among others, will eventually be faced. Regarding consumption, salmonids are preferred to many other protein sources. Regarding trade, increasing the proportion of imports will outsource the impacts. Besides improving technology and dematerialization, other viewpoints may also be needed. Reducing the total amount of nutrients cycling in energy systems and eventually contributing to NOx emissions needs to be emphasized. Considering aquaculture emissions, nutrient cycles can be partly closed through using local fish as feed replacing imported feed. In particular, the reduction of CO2 emissions in the future is a very challenging task when considering the necessary rates of dematerialisation and decarbonisation (paper II). Climate change mitigation may have to focus on other greenhouse gases than CO2 and on the potential role of biomass as a carbon sink, among others. The global population is growing and scaling up the environmental impact. Population issues and growing affluence must be considered when discussing emission reductions. Climate policy has only very recently had an influence on emissions, and strong actions are now called for climate change mitigation. Environmental policies in general must cover all the regions related to production and impacts in order to avoid outsourcing of emissions and leakage effects. The macro-level drivers affecting changes in emissions can be identified with the ImPACT framework. Statistics for generally known macro-indicators are currently relatively well available for different countries, and the method is transparent. In the papers included in this study, a similar method was successfully applied in different types of case studies. Using transparent macro-level figures and a simple top-down approach are also appropriate in evaluating and setting international emission reduction targets, as demonstrated in papers II and IV. The projected rates of population and affluence growth are especially worth consideration in setting targets. However, sensitivities in calculations must be carefully acknowledged. In the basic form of the ImPACT model, the economic intensity of consumption and emission intensity of use are included. In seeking to examine consumption but also international trade in more detail, imports were included in paper III. This example demonstrates well how outsourcing of production influences domestic emissions. Country-specific production-based emissions have often been used in similar decomposition analyses. Nevertheless, trade-related issues must not be ignored.

Highly efficient spin-conversion effect leading to energy up-converted electroluminescence in singlet fission photovoltaics

Free charge generation in donor-acceptor (D-A) based organic photovoltaic diodes (OPV) progresses through formation of charge-transfer (CT) and charge-separated (CS) states and excitation decay to the triplet level is considered as a terminal loss. On the other hand a direct excitation decay to the triplet state is beneficial for multiexciton harvesting in singlet fission photovoltaics (SF-PV) and the formation of CT-state is considered as a limiting factor for multiple triplet harvesting. These two extremes when present in a D-A system are expected to provide important insights into the mechanism of free charge generation and spin-character of bimolecular recombination in OPVs. Herein, we present the complete cycle of events linked to spin conversion in the model OPV system of rubrene/C60. By tracking the spectral evolution of photocurrent generation at short-circuit and close to open-circuit conditions we are able to capture spectral changes to photocurrent that reveal the triplet character of CT-state. Furthermore, we unveil an energy up-conversion effect that sets in as a consequence of triplet population build-up where triplet-triplet annihilation (TTA) process effectively regenerates the singlet excitation. This detailed balance is shown to enable a rare event of photon emission just above the open-circuit voltage (VOC) in OPVs.

The incongruous observation of magnetic phase separation in La0.85Sr0.15CoO3 spin glass system

Phase separation (PS) in hole-doped cobaltites (La1-xSrxCoxO3) is drawing renewed interest recently. In particular, the magnetic behavior of La0.85Sr0.15CoO3 has been subjected to a controversial debate for the past several years; while some groups show evidence for magnetic PS, others show spin glass (SG) behavior. Here, an attempt is made to resolve the controversy related to ``PS versus SG'' behavior in this compound. We present the results of a comprehensive investigation of the dc magnetization, ac susceptibility, and the magnetotransport properties of La0.85Sr0.15CoO3 samples. We contemplate that the magnetic PS in La0.85Sr0.15CoO3 is neither intrinsic nor inherent, but it is a consequence of the preparation conditions. It is realized that a low temperature annealed (LTA) sample shows PS whereas the high temperature annealed (HTA) sample shows SG behavior. The Brillouin-like behavior of field cooled dc magnetization and apparently no frequency dependent peak shift in ac susceptibility for the LTA sample characterize it to be of ferromagneticlike whereas a kink in field cooled dc magnetization and a considerable amount (similar to 3 K) of frequency dependent peak shift in the ac susceptibility for the HTA sample characterize it to be of SG state. The magnetotransport properties show that the HTA sample is more semiconducting as compared to the LTA sample. This is interpreted in terms of the presence of isolated as well as coalescing metallic ferromagnetic clusters in the case of LTA sample. The magnetoresistance (MR) at 10 K for the HTA sample exhibits a huge value (similar to 65%) as compared to the LTA sample, and it monotonically decreases with the rise in temperature. Such a high value of MR in the case of HTA sample is strongly believed to be due to the spin dependent part of random potential distribution. Further, the slow decay of remnant magnetization with progress of time and the existence of hysteresis at higher temperatures (up to 200 K) in the case of LTA sample as compared to the HTA sample clearly unveil different magnetic states associated with them.