Transport properties of nonideal systems with isotropic pair interactions between particles

Results are presented from numerical investigations of the mass transfer and pair correlation in systems of interacting grains for different types of isotropic interaction potentials. The parameters are determined that govern the transport properties of nonideal dissipative systems with a large variety of model potentials. An analytic approximation for the dust grain diffusion coefficient in strongly nonideal systems is obtained.     

Determination of pair interaction potential of particles in nonideal dissipative systems

A new method is proposed for reconstructing physical characteristics of nonideal systems by solving of an inverse problem in which the motion of interacting particles is described by a set of Langevin equations. The method allows one to simultaneously determine the interaction potential between dust grains, the friction coefficients, and the external confining potential. The method was verified by solving the problem numerically in a wide range of parameters typical of laboratory dusty plasmas.     

Thermodynamic properties of bulk knitted structures

The mechanism of interaction of an external magnetic field with liquid was proposed. The statistical integral and configurational contributions for a free energy, entropy and specific heat for the soliton model of bulk knitted structures in a magnetic field were calculated. It was shown that the concentration of solitons depends on the effect of external fields. In the specific case of bulk knitted structures (liquid water without magnetic field), the theoretical data are consistent with experimental. The memory effects in systems with hydrogen bonds in magnetic field was explained in the framework of the continuum soliton concept.     

Formation of chain structures in systems of charged grains interacting via isotropic pair potentials

Conditions for the formation of chain structures of charged grains confined in the gravitational field by external electric fields are studied analytically and numerically. The relationships between the parameters of the pair interaction potential, the number of grains, and the electric field gradient in the trap are found. A criterion for the violation of stable equilibrium in a quasi-one-dimensional chain of grains and the formation of a new configuration in the system is proposed.     

A two-dimensional representation of protein structures

A two-dimensional chart of distances between residues of a protein and between their side-chains provides a concise and convenient representation of steric information derived from crystallographic data. Such a chart shows characteristic features corresponding to details of secondary structures of a polypeptide and also reveals which side-chains are interacting in the tertiary structure of the macromolecule.     

Thermodynamic stability and statistical significance of potential stem-loop structures situated at the frameshift sites of retroviruses.

RNA stem-loop structures situated just 3' to the frameshift sites of the retroviral gag-pol or gag-pro and pro-pol regions may make important contributions to frame-shifting in retroviruses. In this study, the thermodynamic stability and statistical significance of such secondary structural features relative to others in the sequence have been assessed using a newly developed method that combines calculations of the lowest free energy of formation of RNA secondary structures and the Monte Carlo simulations. Our results show that stem-loop structures situated just 3' to the frameshift sites are both highly stable and statistically significant relative to others in the gag-pol or gag-pro and pro-pol junction domains (both 300 nucleotides upstream and downstream from the possible frameshift sites are included) of Rous sarcoma virus (RSV), human immunodeficiency virus (HIV-1), bovine leukemia virus (BLV), human T-cell leukemia virus type II (HTLV-II), and mouse mammary tumor virus (MMTV). No other more stable, or significant folding regions are predicted in these domains.     

Thermodynamic properties of purple membrane

We measured the density, expansivity, specific heat at constant pressure, and sound velocity of suspensions of purple membrane from Halobacterium halobium and their constituent buffers. From these quantities we calculated the apparent values for the density, expansivity, adiabatic compressibility, isothermal compressibility, specific heat at constant pressure, and specific heat at constant volume for the purple membrane. These results are discussed with respect to previously reported measurements on globular proteins and lipids. Our data suggest a simple additive model in which the protein and lipid molecules expand and compress independently of each other. However, this simple model seems to fail to describe the specific heat data. Our compressibility data suggest that bacteriorhodopsin in native purple membrane binds less water than many globular proteins in neutral aqueous solution, a finding consistent with the lipid surround of bacteriorhodopsin in purple membrane.     

Decomposing the space of protein quaternary structures with the interface fragment pair library

Background

The physical interactions between proteins constitute the basis of protein quaternary structures. They dominate many biological processes in living cells. Deciphering the structural features of interacting proteins is essential to understand their cellular functions. Similar to the space of protein tertiary structures in which discrete patterns are clearly observed on fold or sub-fold motif levels, it has been found that the space of protein quaternary structures is highly degenerate due to the packing of compact secondary structure elements at interfaces. Therefore, it is necessary to further decompose the protein quaternary structural space into a more local representation.

Results

Here we constructed an interface fragment pair library from the current structure database of protein complexes. After structural-based clustering, we found that more than 90% of these interface fragment pairs can be represented by a limited number of highly abundant motifs. These motifs were further used to guide complex assembly. A large-scale benchmark test shows that the native-like binding is highly likely in the structural ensemble of modeled protein complexes that were built through the library.

Conclusions

Our study therefore presents supportive evidences that the space of protein quaternary structures can be represented by the combination of a small set of secondary-structure-based packing at binding interfaces. Finally, after future improvements such as adding sequence profiles, we expect this new library will be useful to predict structures of unknown protein-protein interactions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0437-4) contains supplementary material, which is available to authorized users.     

Thermodynamic contributions of the reactions of DNA intramolecular structures with their complementary strands

One focus of our research is to further our understanding of the physico-chemical properties of unusual DNA structures and their interaction with complementary oligonucleotides. We have investigated three types of reactions involving the interaction of intramolecular DNA complexes with their complementary single strands of varied length. Specifically, we have used a combination of isothermal titration (ITC) and differential scanning (DSC) calorimetry and spectroscopy techniques to determine standard thermodynamic profiles for the reaction of an i-motif, G-quadruplex, and triplex with their complementary strands. The enthalpies for each reaction are measured directly in ITC titrations and compared with those obtained indirectly from Hess cycles using DSC unfolding data. All reactions investigated yielded favorable free energy contributions, indicating that each single strand is able to invade and disrupt the corresponding intramolecular DNA complex. These favorable free energy terms are enthalpy driven, which result from a compensation of exothermic contributions, due to the formation of additional base-pair stacks (or base-triplet stacks) in the duplex product (or triplex product), immobilization of electrostricted water by the base-pair and base-triplet stacks, and the removal of structural water from the reactant single strands; and endothermic contributions from the disruption of base-base stacking interactions of the reactant single strands. This investigation of nucleic acid reactions has provided new methodology, based on physico-chemical principles, to determine the molecular forces involved in the interactions between DNA nucleic acid structures. This methodology may be used in targeting reactions for the control of gene expression.     

Thermodynamic properties of nucleotide reductase reactions

Recent thermodynamic measurements have made it possible to calculate the apparent equilibrium constants of the ribonucleoside diphosphate reductase reaction and the ribonucleoside triphosphate reductase reaction with various reducing agents. Third law heat capacity measurements on crystals of d-ribose and other calorimetric measurements make it possible to calculate Delta(f)G degrees for D-ribose and two species of D-ribose 5-phosphate. The experimental value of the apparent equilibrium constant K' for the deoxyribose-phosphate aldolase reaction makes it possible to calculate the standard Gibbs energies of formation Delta(f)G degrees for two protonation states of 2'-deoxy-D-ribose 5-phosphate. This shows that Delta(f)G degrees (2'-deoxy-D-ribose 5-phosphate(2)(-)) - Delta(f)G degrees (D-ribose 5-phosphate(2)(-)) = 147.86 kJ mol(-1) at 298.15 K and zero ionic strength in dilute aqueous solutions. This difference between reduced and oxidized forms is expected to apply to D-ribose, D-ribose 1-phosphate, ribonucleosides, and ribonucleotides in general. This expectation is supported by two other enzyme-catalyzed reactions for which apparent equilibrium constants have been determined. The availability of Delta(f)G degrees values for the species of 2'-deoxy-D-ribose and its derivatives makes it possible to calculate standard transformed Gibbs energies of formation of these reactants, apparent equilibrium constants for their reactions, changes in the binding of hydrogen ions in these reactions, and standard apparent reduction potentials of the half reactions involved as a function of pH and ionic strength at 298.15 K. The apparent equilibrium constant for ADP + thioredoxin(red) = 2'-deoxyADP + H(2)O + thioredoxin(ox) is 1.4 x 10(11) at 298.15 K, pH 7, and 0.25 M ionic strength.     

Thermodynamic description of oligonucleotide self-association in DNA concatamer structures

A scheme of equilibrium formation of concatamers by two different oligonucleotides has been considered. It is shown that in the general case, the dependence of the concentration of oligonucleotide components on temperature cannot be found in analytical form. Therefore, it is impossible to find the thermodynamic parameters of concatamer formation (ΔH ⁰, ΔS ⁰) and melting temperatures by analyzing the profiles of thermal denaturation of oligonucleotide complexes. An algorithm for numerical solution of implicit dependences has been developed. A number of approaches are considered that simplify the analysis of heat denaturation curves for concatamer complexes. It is shown that the dependence of the efficiency of concatamerization on temperature can be described analytically when duplex fragments have close stability and there is no cooperativity at the oligonucleotide junction. In this case, the dependence of melting temperature on thermodynamic parameters and oligonucleotide concentration has the same form as in the case of the duplex structure formed by a pair of non-self-complementary oligonucleotides. The ability of various model approaches to describe the experimental curves of concatamer heat denaturation is evaluated. For concatamer structures used as signal amplifiers in DNA hybridization analysis, a function is introduced that shows the relative contribution of a concatamer of given length to the magnitude of signal amplification. The dependence of the maximum of this function on the concentration of oligonucleotides, the thermodynamic characteristics of their complexes, and temperature has been determined. It is shown by the gel retardation assay that the function of the length distribution of concatamers qualitatively correlates with the experimental dependences.     

Thermodynamic properties of BPTI variants with highly simplified amino acid sequences

We report the first detailed thermodynamic analysis of simplified proteins by differential scanning calorimetry (DSC). The experiments were carried out with five simplified BPTI variants, whose structures and activities have been reported, in which several residues not essential for specifying the tertiary structure were replaced by alanine. In most aspects, the thermodynamics of simplified proteins were very similar to, if not essentially identical with, those of natural proteins. In particular, they undergo a highly cooperative two-state thermal unfolding process with a large enthalpy change, which is a thermodynamic hallmark of the native state of natural globular proteins. Furthermore, the specific enthalpy and entropy changes upon unfolding at 110 degrees C were close to values invariably observed for small natural globular proteins (55 J g(-1) and ~16 J K(-1) g(-1), respectively). On the other hand, two simplified BPTI variants, BPTI-21 and BPTI-22 (containing 21 and 22 alanine residues), were enthalpically stabilized while entropically destabilized with respect to the reference BPTI-[5,55] molecule. This peculiar type of entropy-enthalpy compensation is in sharp contrast to the usual enthalpy destabilization/entropy stabilization observed in mutational studies of natural proteins. Overall, we conclude that a thermodynamic native state can be achieved by proteins encoded with extensively simplified sequences.     

Comparative ab initio prediction of gene structures using pair HMMs

We present a novel comparative method for the ab initio prediction of protein coding genes in eukaryotic genomes. The method simultaneously predicts the gene structures of two un-annotated input DNA sequences which are homologous to each other and retrieves the subsequences which are conserved between the two DNA sequences. It is capable of predicting partial, complete and multiple genes and can align pairs of genes which differ by events of exon-fusion or exon-splitting. The method employs a probabilistic pair hidden Markov model. We generate annotations using our model with two different algorithms: the Viterbi algorithm in its linear memory implementation and a new heuristic algorithm, called the stepping stone, for which both memory and time requirements scale linearly with the sequence length. We have implemented the model in a computer program called DOUBLESCAN. In this article, we introduce the method and confirm the validity of the approach on a test set of 80 pairs of orthologous DNA sequences from mouse and human. More information can be found at: http://www.sanger.ac.uk/Software/analysis/doublescan/     

Nonlinear dynamics of two-dimensional cardiac action potential duration mapping model with memory

The aim of this work is the analysis of the nonlinear dynamics of two-dimensional mapping model of cardiac action potential duration (2D-map APD) with memory derived from one dimensional map (1D-map). Action potential duration (APD) restitution, which relates APD to the preceding diastolic interval (DI), is a useful tool for predicting cardiac arrhythmias. For a constant rate of stimulation the short action potential during alternans is followed by a longer DI and inversely. It has been suggested that these differences in DI are responsible for the occurrence and maintenance of APD alternans. We focus our attention on the observed bifurcations produced by a change in the stimulation period and a fixed value of a particular parameter in the model. This parameter provides new information about the dynamics of the APD with memory, such as the occurrence of bistabilities not previously described in the literature, as well as the fact that synchronization rhythms occur in different ways and in a new fashion as the stimulation frequency increases. Moreover, we show that this model is flexible enough as to accurately reflect the chaotic dynamics properties of the APD: we have highlighted the fractal structure of the strange attractor of the 2D-map APD, and we have characterized chaos by tools such as the calculation of the Lyapunov exponents, the fractal dimension and the Kolmogorov entropy, with the next objective of refining the study of the nonlinear dynamics of the duration of the action potential and to apply methods of controlling chaos.

     

Thermodynamic properties of an intramolecular DNA four-way junction

We have investigated the thermodynamic properties of two homologous DNA four-way junctions, J4 and J4M, based on 46-mer linear DNA molecules. J4 and J4M have the same base sequence with the only difference that the latter contains an uncharged methylene-acetal linkage, -O3'-CH2-O5', instead of the phosphodiester linkage, -O3'-PO2-O5'-, between the residues T18 and C19. The comparison of the thermal unfolding of the J4 junction and J4M junction serves to investigate the effect of the uncharged methylene-acetal linkage on the stability of the junction. Our analysis is based on CD, UV absorbance spectroscopy, DSC, and chemical footprinting. The aim is to characterize in detail the structure and stability of the junctions. As demonstrated before by NMR, in the presence of 5 mM MgCl2 +/- 50 mM NaCl, both J4 and J4M form a complete four-way junction. This is now evidenced by protection from OsO4 cleavage (chemical footprinting). We can assume that full base pairing occurs throughout the arms even at the center of the junction. CD spectra suggest that the helices within the junctions adopt the regular B-DNA conformation. Almost identical melting temperatures and unfolding enthalpies are obtained for J4 and J4M both by UV and DSC. Furthermore, the Van't Hoff enthalpy (DeltaHVH) derived from UV melting equals the calorimetric enthalpy (DeltaHcal), which means that the melting process of the structures proceeds in a two-state manner. All results taken together support the conclusion that there are no major conformational and energetic differences between J4 and J4M. The inclusion of the uncharged methylene-acetal group into the junction has no effect on its stability.     

Thermodynamic and mechanical properties of model mitochondrial membranes

Cardiolipin is a unique four-tailed, doubly negatively charged lipid found predominantly within the inner mitochondrial membrane, and is thought to be influential in determining membrane potential and permeability. To determine the role of cardiolipin in modulating the properties of membranes, this study investigates the thermodynamics of mixed cardiolipin and phosphatidylcholine monolayers and bilayers. Gibbs free energy analysis of mixed monolayers indicates that at low cardiolipin concentrations (5-10 mol%), there is a positive deviation from ideality on a pure water subphase, while at physiological salt concentrations a negative deviation from ideality is observed. The mechanical properties of bilayers containing cardiolipin were measured using micropipette aspiration. Both apparent area compressibility modulus, as well as lysis tension, decrease with increasing cardiolipin content. This destabilization indicates a decrease in the cohesive energy of the membrane. This interplay between interactions of lipids in monolayers and bilayers, suggests cardiolipin plays a dual role in modulating membrane properties. Cardiolipin enhances lateral interactions between lipids within monolayer leaflets, while simultaneously decreasing the cohesive energy of membranes at physiologically relevant concentrations. Taken together, these findings correlate with the decreased permeability and creation of folds in the inner mitochondrial membrane.