Representation of Coarse-grained Potentials for Polymer Simulations

In order to be able to simulate long time and large space scale properties of polymer melts one has to resort to coarse grained models, for example by subdividing all polymers into parts and restricting attention to the center of mass positions and velocities of these parts. The dynamics of these variables is governed by Langevin equations in which the free energy obtained by integrating the remaining variables provides the potential of the conservative forces. In general this leads to many particle interactions on the coarse-grained level. Methods suggested in the literature to represent these many particle interactions by effective two body interactions are reviewed and a new method, based on the Gibbs-Bogoliubov inequality, is proposed. The reason why none of these methods is able to reproduce the pressure of the underlying atomistic model is discussed.     

The concentrations and thermodynamic activities of cations in intact Bryopsis chloroplasts

The internal cation levels of chloroplasts isolated from a green sea alga, Bryopsis maxima, were studied. Atomic absorption spectroscopy, combined with the determination of the sorbitol-impermeable and water-permeable spaces, revealed that chloroplasts contain an extremely high concentration of K⁺ and high levels of Na⁺, Mg²⁺ and Ca²⁺. A method was developed to estimate the thermodynamic activities of monovalent and divalent cations present in chloroplasts. pH changes induced by the addition of an ionophore (plus an H⁺ carrier), which makes the outer limiting membranes of chloroplasts permeable to both a cation and H⁺, were determined. Provided that the external pH was set equal to the internal pH, the internal concentration of the cation was estimated by determining the external cation concentration which gave rise to no electrochemical potential difference of the cation and hence no pH change on addition of the ionophore. The internal pH was determined by measuring distributions of radioactive methylamine and 5,5-dimethyloxazolidine-2,4-dione between the chloroplast and medium (Heldt, H.W., Werdan, K., Milovancev, M. and Geller, G. (1973) Biochim. Biophys. Acta 314, 224–241). The internal pH was also estimated by measuring pH changes caused by the disruption of the outer limiting membrane with Triton X-100. The results indicate that a significant part of the monovalent cations and most of the divalent cations are attracted into a diffuse layer adjacent to the negatively charged surfaces of membranes and proteins, or form complexes with organic and inorganic compounds present in the intact chloroplasts.     

The Influence of Nearest Neighbours on the Efficiency of Coaxial Stacking at Contiguous Stacking Hybridization of Oligodeoxyribonucleotides

Contiguous stacking hybridization of oligodeoxyribonucleotides with a stem of preformed minihairpin structure of a DNA template was studied with the use of UV‐melting technique. It was shown that the free‐energy of the coaxial stacking interaction (ΔG°_ST at 37°C, 1 M NaCl, pH 7.4) at the complementary interface XA^*pTY/ZATV (an asterisk stands for a nick) strongly depends on the type of nearest neighbor bases X and Y flanking the nicked dinucleotide step. The maximum efficiency of the coaxial stacking was observed for the PuA^*pTPy/PuATPy interface, whereas the minimum efficiency was obtained for the PyA^*pTPu/PyATPu interface. A 5′‐phosphate residue in the nick enhances the coaxial stacking. In dependence on duplex structure the observed efficiency of A^*T/AT coaxial stacking varied from (? 0.97 kcal/mol) for unphosphorylated TA^*TA/TATA interface to three‐fold higher value (? 2.78 kcal/mol) for GA^*pTT/AATC interface.     

Ago-Allosteric Modulation and Other Types of Allostery in Dimeric 7TM Receptors

Conventionally, an allosteric modulator is neutral in respect of efficacy and binds to a receptor site distant from the orthosteric site of the endogenous agonist. However, recently compounds being ago-allosteric modulators have been described i.e., compounds acting both as agonists on their own and as enhancers for the endogenous agonists in both increasing agonist potency and providing additive efficacy—superagonism. The additive efficacy can also be observed with agonists, which are neutral or even negative modulators of the potency of the endogenous ligand. Based on the prevailing dimeric concept for 7TM receptors, it is proposed that the ago-allosteric modulators bind in the orthosteric binding site, but–importantly–in the “other” or allosteric protomer of the dimer. Hereby, they can act both as additive co-agonists, and through intermolecular cooperative effects between the protomers, they may influence the potency of the endogenous agonist. It is of interest that at least some endogenous agonists can only occupy one protomer of a dimeric 7TM receptor complex at a time and thereby they leave the orthosteric binding site in the allosteric protomer free, potentially for binding of exogenous, allosteric modulators. If the allosteric modulator is an agonist, it is an ago-allosteric modulator; if it is neutral, it is a classical enhancer. Molecular mapping in hetero-dimeric class-C receptors, where the endogenous agonist clearly binds only in one protomer, supports the notion that allosteric modulators can act through binding in the “other” protomer. It is suggested that for the in vivo, clinical setting a positive ago-allosteric modulator should be the preferred agonist drug.     

Synthesis and Hybridization Properties of Oligonucleotides Containing (2 S ,3 R )-9-(2,3,4-Trihydroxybutyl)Adenine

The synthesis and properties of oligonucleotides (ONs) containing 9-(2,3,4-trihydroxybutyl)adenine, A ^C2 and A ^C3, are described. The ON containing A ^C2 involves the 3′ → 4′ and 3′ → 5′ phosphodiester linkages in the strand, whereas that containing A ^C3 possesses the 3′ → 4′ and 2′ → 5′ phosphodiester linkages. It was found that incorporation of the analogs, A ^C2 or A ^C3, into ONs significantly reduces the thermal and thermodynamic stabilities of the ON/DNA duplexes, but does not largely decrease the thermal and thermodynamic stabilities of the ON/RNA duplexes as compared with the case of the ON/DNA duplexes. It was revealed that the base recognition ability of A ^C2 is greater than that of A ^C3 in the ON/RNA duplexes.     

Depletion interaction in colloid/polymer mixtures: application of density functional theory

Insert-route density functional approach (IRDFT), modified fundamental measure theory (MFMT) and thermodynamic perturbation theory (TPT1 and TPT2) are combined to study the depletion force between colloidal particles in hard sphere/hard sphere chain mixtures which represent a model of systems containing colloids dispersed in an athermal polymer solution. The predicted results are compared to simulations showing the reliability of the method used which captures the main characteristics of depletion interaction between colloids induced by polymers. Results of TPT2 are slightly more repulsive and better than that of TPT1 especially when the inter-particle distance is small than the diameter of polymer segment indicating the essential influence of the three-body correlations. Effects of the polymer density, polymer chain length and size ratio of colloid to polymer segment on the depletion force are studied in detail. Due to a little deterioration of the prediction in the high density region, further improvement is anticipated to better balance the competition between the excluded-volume effect and the chain connectivity.     

Base Sequence Effects in Double Helical DNA. I. Potential Energy Estimates of Local Base Morphology

Abstract

A series of potential energy calculations have been carried out to estimate base sequence dependent structural differences in B-DNA. Attention has been focused on the simplest dimeric fragments that can be used to build long chains, computing the energy as a function of the orientation and displacement of the 16 possible base pair combinations within the double helix. Calculations have been performed, for simplicity, on free base pairs rather than complete nucleotide units. Conformational preferences and relative flexibilities are reported for various combinations of the roll, tilt, twist, lateral displacement, and propeller twist of individual residues. The predictions are compared with relevant experimental measures of conformation and flexibility, where available. The energy surfaces are found to fit into two distinct categories, some dimer duplexes preferring to bend in a symmetric fashion and others in a skewed manner. The effects of common chemical substitutions (uracil for thymine, 5-methyl cytosine for cytosine, and hypoxanthine for guanine) on the preferred arrangements of neighboring residues are also examined, and the interactions of the sugar-phosphate backbone are included in selected cases. As a first approximation, long range interactions between more distant neighbors, which may affect the local chain configuration, are ignored. A rotational isomeric state scheme is developed to describe the average configurations of individual dimers and is used to develop a static picture of overall double helical structure. The ability of the energetic scheme to account for documented examples of intrinsic B-DNA curvature is presented, and some new predictions of sequence directed chain bending are offered.     

Unusual Protonated Structure in the Homopurine · Homopyrimidine Tract of Supercoiled and Linearized Plasmids Recognized by Chemical Probes

Abstract

Plasmid pEJ4, which is a derivative of pUC19 containing an insert with 60-bp-long · homopurine homopyrimidine tract from sea urchin P. miliaris histone gene spacer, was studied by chemical probes of the DNA structure osmium tetroxide and glyoxal. The former probe reacts with pyrimidine bases, while the latter forms a stable product only with guanine residues. These probes can thus be applied as specific probes for the homopyrimidine and homopurine strands, respectively.

At pH 6.0 the site-specific modification of the homopurine · homopyrimidine tract by both probes was observed at native superhelical density of the plasmid. In the linear plasmid under the same conditions this modification was absent; it appeared, however, at more acid pH values. In supercoiled DNA the hypersensitivity of the homopurine homopyrimidine tract to osmium tetroxide did not substantially change when pH was decreased from 6.0 to 4.0. Changes in NaCl concentration at pH 4.5 did not influence the hypersensitivity to osmium tetroxide; at pH 6.0 this hypersensitivity decreased with increasing NaCl concentration. These results thus show that the chemical probes recognize an unusual protonated structure containing unpaired bases or non-Watson-Crick base pairs. At pH 5.6 the site-specific modification occurred at or near to the middle of the homopurine · homopyrimidine tract, suggesting that a hairpin may be involved in the unusual structure under the given conditions. From the models suggested so far for the unusual structure of homopurine · homopyrimidine tracts our results fit best the protonated triplex H form suggest by V.I. Lyamichev, S.M. Mirkin and M.D. Frank-Kamenetskii, J. Biomol. Struct. Dyn. 3, 667 (1986).     

A thermodynamic approach for assessing agroecosystem sustainability

By revisiting theoretical concepts in biogeography and the importance of thermodynamic laws in biosphere-atmosphere interactions, ecological sustainability in agricultural systems may be better defined. In this case study, we employed a multidisciplinary methodology for exploring agroecosystem sustainability by using eddy covariance (EC) data to compute thermodynamic entropy production (σ) and relate it to water, energy and carbon cycling in croplands and grasslands of the Central US. From 2002 to 2012, the biophysical metric of σ was compared across AmeriFlux sites, each with site-specific land management practices of irrigation, crop rotation, and tillage. Results show that σ is most correlated with net ecosystem exchange (NEE) of carbon, and when cropland and grassland sites are close to being carbon neutral, σ values range from 0.51–1.0 W K⁻¹ m⁻² for grasslands, 0.81–1.0 W K⁻¹ m⁻² for rainfed croplands, and 0.81–1.1 W K⁻¹ m⁻² for irrigated croplands. Irrigated maize stressed by hydrologic and high temperature anomalies associated with the 2012 drought exhibit the greatest increase in σ, indicating the possibility of decreased sustainability compared to rainfed croplands and grasslands. These results suggest that maximizing carbon uptake with irrigation and fertilizer use tends to move agroecosystems further away from thermodynamic equilibrium, which has implications for ecological sustainability and greenhouse gas (GHG) mitigation in climate-smart agriculture. The underlying theoretical concepts, multidisciplinary methodology, and use of eddy covariance data for biophysical indicators in this study contribute to a unique understanding of ecological sustainability in agricultural systems.     

Amino-acid substitutions at the fully exposed P1 site of bovine pancreatic trypsin inhibitor affect its stability

It is widely accepted that solvent-exposed sites in proteins play only a negligible role in determining protein energetics. In this paper we show that amino acid substitutions at the fully exposed Lys15 in bovine pancreatic trypsin inhibitor (BPTI) influenced the CD- and DSC-monitored stability: The T(den) difference between the least (P1 Trp) and the most stable (P1 His) mutant is 11.2 degrees C at pH 2.0. The DeltaH(den) versus T(den) plot for all the variants at three pH values (2.0, 2.5, 3.0) is linear (DeltaC(p,den) = 0.41 kcal* mole(-1) * K(-1); 1 cal = 4.18 J) leading to a DeltaG(den) difference of 2.1 kcal*mole(-1). Thermal denaturation of the variants monitored by CD signal at pH 2.0 in the presence of 6 M GdmCl again showed differences in their stability, albeit somewhat smaller (DeltaT(den) =7.1 degrees C). Selective reduction of the Cys14-Cys 38 disulfide bond, which is located in the vicinity of the P1 position did not eliminate the stability differences. A correlation analysis of the P1 stability with different properties of amino acids suggests that two mechanisms may be responsible for the observed stability differences: the reverse hydrophobic effect and amino acid propensities to occur in nonoptimal dihedral angles adopted by the P1 position. The former effect operates at the denatured state level and causes a drop in protein stability for hydrophobic side chains, due to their decreased exposure upon denaturation. The latter factor influences the native state energetics and results from intrinsic properties of amino acids in a way similar to those observed for secondary structure propensities. In conclusion, our results suggest that the protein-stability-derived secondary structure propensity scales should be taken with more caution.