A semiempirical extension of polyelectrolyte theory to the treatment of oligoelectrolytes: Application to oligonucleotide helix-coil transitions

The interaction of counterions with a suitably long, charged oligomer appears susceptible to treatment in the context of polyelectrolyte theory by the introduction of an end-effect parameter that reflects the reduced association of counterions with the terminal regions of the oligo-ion. Use of a physically reasonable value for the end-effect parameter provides excellent agreement between theory and the experimental data of Elson, Scheffler, and Baldwin [J. Mol. Biol. 54 , 401–415 (1970)] on the dependences of melting temperature on salt concentration and chain length for a series of hairpin helices formed by d(TA) oligomers. The differences in behavior expected for hairpin, dimer, and oligomer-polymer helices are discussed. The salt dependence of the end-joining equilibrium investigated for λ DNA by Wang and Davidson [Cold Spring Harbor Symp. Quant. Biol. 33 , 409–415 (1968)] is treated as an oligomer–polymer interconversion. The dependence of equilibrium constant for this reaction on counterion concentration is in good agreement with that predicted by theory for an end-region totalling 24 nucleotides, the known length of the λ ends.     

Nucleic acid base-pairing and N-methylacetamide self-association in chloroform: affinity and conformation.

A recently developed computational method, 'mining minima', is used to examine the hydrogen-bonding interactions of nucleic acid base-pairs and of the N-methylacetamide homodimer in chloroform. The mining minima algorithm aggressively samples molecular conformations, identifies the most important local minima, and computes their contributions to the overall free energy of the system. Here, the CHARMM 98 parameter set is used for the potential energy and the generalized Born/surface area solvent model is used to account for the influence of the solvent. Good agreement with experiment is obtained for the non-covalent binding affinities of a series of complexes. The computational approach used here is applicable to a range of molecular systems.     

Yeast基因组编码区特征参数的研究

以碱基成分偏移量D值为基本参数定义参数d,以d为Yeast编码区的特征参数,对Yeast的第1、2、3类ORF（open reading frame)进行了统计,得到d的特征参数区间,并且,以此区间为标准为Yeast的6类ORF,以及5′帽、3′尾、内含子、组分随机序列等非编码序列进行了检验。结果表明,d作编码区的特征参数是可行的,它可以很好地区分编码序列和非编码序列。别外,又讨论了参数d与基因表达水平（用CAI值来衡量）的关系。发现,参数d与基因表达水平成很好的正相关关系。发现密码子的第1位点和第2位点的某些碱基分布与基因表达水平有关。     

桔青霉摇瓶发酵生产核酸酶P1的动力学研究

研究了桔青霉发酵生产核酸酶P1的发酵动力学特性：以Logistic方程和Luedeking—Piret方程为基础,进行最优参数估计和非线性拟合,得到了描述整个发酵过程中的茵体生长、产物合成和基质消耗的动力学模型。对实验数据与模型预测值进行比较,发现模型预测值与实验数据能较好地拟合,基本上反映了桔青霉发酵过程的动力学特征,为以后进一步研究和预测核酸酶P1发酵过程奠定了理论基础。     

Heat Generation/Absorption Effects in a Boundary Layer Stretched Flow of Maxwell Nanofluid: Analytic and Numeric Solutions

Analysis has been done to investigate the heat generation/absorption effects in a steady flow of non-Newtonian nanofluid over a surface which is stretching linearly in its own plane. An upper convected Maxwell model (UCM) has been utilized as the non-Newtonian fluid model in view of the fact that it can predict relaxation time phenomenon which the Newtonian model cannot. Behavior of the relaxations phenomenon has been presented in terms of Deborah number. Transport phenomenon with convective cooling process has been analyzed. Brownian motion “D_b” and thermophoresis effects “D_t” occur in the transport equations. The momentum, energy and nanoparticle concentration profiles are examined with respect to the involved rheological parameters namely the Deborah number, source/sink parameter, the Brownian motion parameters, thermophoresis parameter and Biot number. Both numerical and analytic solutions are presented and found in nice agreement. Comparison with the published data is also made to ensure the validity. Stream lines for Maxwell and Newtonian fluid models are presented in the analysis.     

Kinetics of denitritification and denitratification in anoxic filters

Denitritification and denitratification in anoxic filters were performed to generate experimental data. Also, a kinetic model of denitratification that accounts for intrinsic biokinetics and hydrodynamic behavior of the biofilter is proposed. In denitritification, the simulated results are in good agreement with the experimental data; and a higher nitrite influent concentration gives a higher nitrite reduction efficiency if the denitrifying loading is kept the same. In denitratification, the intermediate nitrite tends to accumulate, and a higher denitrifying loading results in a higher nitrite effluent concentration. By inserting biological and physical parameter values into the kinetic model, the variations in distributed fractions of nitrate-reductase (f) and nitrite-reductase (1-f) with different denitrifying loadings can be estimated by fitting in experimental data. The estimated f increased with an increase in denitrifying loading, implying that a higher denitrifying loading results in a higher nitrite effluent concentration. From parametric sensitivity analyses, the parameter f is more sensitive than other biological and physical parameters. Accordingly, the proposed kinetic model of denitratification can be used to predict the treatment performance of anoxic filters appropriately. Copyright 1998 John Wiley & Sons, Inc.     

Modelling root-induced solubilization of nutrients

Theory for coupled diffusion processes in soil is briefly described and three examples of its application to understand root-induced solubilization of nutrients given. The examples are: (1) solubilization of P through root-induced pH changes in the rhizosphere of rice plants growing in flooded soil; (2) solubilization of P through excretion of organic chelating agents from rice roots growing in aerobic soil; and (3) the effects of root geometry on P solubilization, particularly cylindrical versus planar geometry and the effect of excretion of a solubilizing agent being localized along the root axis. The theory is tested by comparing measured concentration profiles of P near roots with the predictions of the theory made using independently measured parameter values. In the examples given, the agreement between the observed and predicted concentration profiles is very good, indicating that the theory is sound and the processes involved well understood.     

Relationship between preferential interaction of a protein in an aqueous mixed solvent and its solubility

The present paper is devoted to the derivation of a relation between the preferential solvation of a protein in a binary aqueous solution and its solubility. The preferential binding parameter, which is a measure of the preferential solvation (or preferential hydration) is expressed in terms of the derivative of the protein activity coefficient with respect to the water mole fraction, the partial molar volume of protein at infinite dilution and some characteristics of the protein-free mixed solvent. This expression is used as the starting point in the derivation of a relationship between the preferential binding parameter and the solubility of a protein in a binary aqueous solution. The obtained expression is used in two different ways: (1) to produce a simple criterion for the salting-in or salting-out by various cosolvents on the protein solubility in water, (2) to derive equations which predict the solubility of a protein in a binary aqueous solution in terms of the preferential binding parameter. The solubilities of lysozyme in aqueous sodium chloride solutions (pH=4.5 and 7.0), in aqueous sodium acetate (pH=8.3) and in aqueous magnesium chloride (pH=4.1) solutions are predicted in terms of the preferential binding parameter without any adjustable parameter. The results are compared with experiment, and for aqueous sodium chloride mixtures the agreement is excellent, for aqueous sodium acetate and magnesium chloride mixtures the agreement is only satisfactory.     

Preferential hydration and solubility of proteins in aqueous solutions of polyethylene glycol

This paper is focused on the local composition around a protein molecule in aqueous mixtures containing polyethylene glycol (PEG) and the solubility of proteins in water + PEG mixed solvents. Experimental data from literature regarding the preferential binding parameter were used to calculate the excesses (or deficits) of water and PEG in the vicinity of β-lactoglobulin, bovine serum albumin, lysozyme, chymotrypsinogen and ribonuclease A. It was concluded that the protein molecule is preferentially hydrated in all cases (for all proteins and PEGs investigated). The excesses of water and deficits of PEG in the vicinity of a protein molecule could be explained by a steric exclusion mechanism, i.e. the large difference in the sizes of water and PEG molecules.

The solubility of different proteins in water + PEG mixed solvent was expressed in terms of the preferential binding parameter. The slope of the logarithm of protein (lysozyme, β-lactoglobulin and bovine serum albumin) solubility versus the PEG concentration could be predicted on the basis of experimental data regarding the preferential binding parameter. For all the cases considered (various proteins, various PEGs molecular weights and various pHs), our theory predicted that PEG acts as a salting-out agent, conclusion in full agreement with experimental observations. The predicted slopes were compared with experimental values and while in some cases good agreement was found, in other cases the agreement was less satisfactory. Because the established equation is a rigorous thermodynamic one, the disagreement might occur because the experimental results used for the solubility and/or the preferential binding parameter do not correspond to thermodynamic equilibrium.     

Interaction of anticancer drug cisplatin with guanine: density functional theory and surface-enhanced Raman spectroscopy study

Surface-enhanced Raman spectroscopy (SERS) in a silver sol assisted by density functional theory (DFT) calculations is shown to be a promising tool in the characterization of platinum complexes and their interaction with nucleic acid bases. This is demonstrated using cisplatin and guanine as a model. The energies and geometric parameters of cisplatin, guanine, and their reaction products are calculated at Becke's nonlocal three parameter exchange and correlation functional and the Lee-Yang-Parr correlation functional level using the 6-31++G(d,p) basis set on the light elements and the effective core potential by Hay and Wadt on platinum. Available X-ray crystallography data are mostly in agreement with predictions within the experimental precision level, although Pt-N bond lengths tend to be systematically overestimated. The normal Raman spectrum of cisplatin is assigned. The SERS spectra of cisplatin and its reaction product with guanine are measured from 10(-6) M aqueous solution. The observed spectral changes in the SERS spectrum of guanine upon cisplatin binding are modeled by DFT calculations. The best agreement between theory and experiment is achieved when the adsorbed reaction product is assumed to be the 1:1 adduct cis-Pt(NH3)2ClG in which Pt is bound to N7 and guanine is deprotonated at N9.     

Measurement of interrater agreement with adjustment for covariates

The kappa coefficient measures chance-corrected agreement between two observers in the dichotomous classification of subjects. The marginal probability of classification by each rater may depend on one or more confounding variables, however. Failure to account for these confounders may lead to inflated estimates of agreement. A multinomial model is used that assumes both raters have the same marginal probability of classification, but this probability may depend on one or more covariates. The model may be fit using software for conditional logistic regression. Additionally, likelihood-based confidence intervals for the parameter representing agreement may be computed. A simple example is discussed to illustrate model-fitting and application of the technique.     

Biochemical characterization and membrane fluidity of membranous vesicles isolated from boar seminal plasma

Mammalian seminal plasma contains membranous vesicles (MV), which differ in composition and origin. Among these particles, human prostasomes and equine prostasome-like MV have been the most studied. The aim of the present work is to characterize the biochemical composition and membrane fluidity of MV isolated from boar seminal plasma. The MV from boar seminal plasma were isolated by ultracentrifugation and further purification by gel filtration on Sephadex G-200. The MV were examined by electron microscopy (EM), amount of cholesterol, total phospholipid, protein content, and phospholipid composition were analyzed. Membrane fluidity of MV and spermatozoa were estimated from the electron spin resonance (ESR) spectra of the 5-doxilstearic acid incorporated into the vesicle membranes by the order parameter (S). The S parameter gives a measure of degree of structural order in the membrane and is defined as the ratio of the spectral anisotropy in the membranes to the maximum anisotropy obtained in a rigidly oriented system. The S parameter takes into consideration that S = 1 for a rapid spin-label motion of about only one axis and S = 0 for a rapid isotropic motion. Intermediate S values between S = 0 and S = 1 represents the consequence of decreased membrane fluidity. The EM revealed the presence of bilaminar and multilaminar electron-dense vesicles. Cholesterol to phospholipid molar ratio from the isolated MV was 1.8. Phospholipid composition showed a predominance of sphingomyelin. The S parameter for porcine MV and for boar spermatozoa was 0.73 +/- 0.02 and 0.644 +/- 0.008, respectively, with the S for MV being greater (p < 0.001) than the S for spermatozoa. The high order for S found for boar MV was in agreement with the greater cholesterol/phospholipids ratio and the lesser ratio for phosphatidylcholine/sphingomyelin. Results obtained in the present work indicate that MV isolated from boar semen share many biochemical and morphological characteristics with equine prostasome-like MV and human prostasomes. The characteristics of the porcine MV of the seminal plasma, however, differed from those of boar sperm plasma membranes.     

Generalized ordering parameter S2 for N-H peptide bonding as measure of the conformational flexibility of proteins: comparison of algorithms of S2 calculation from molecular dynamics simulation data

NMR spectroscopy in combination with molecular dynamics (MD) simulation technique is a powerful tool for investigation of biomolecule conformations. In the present study we compare the two algorithms of calculation of the generalized order parameter S2 for the peptide N-H bond from MD simulation data. The program based on this algorithm has been written and the effect of an adjustible parameter on results of calculation was studied. HIV-1 protease MD simulation was performed during 7.36 ns using GROMACS 3.1.4. software package. The calculated order parameter values are in good agreement with those obtained by other authors using NMR.     

Influence of sequence on the conformation of the B-DNA helix.

We have tried to ascertain whether the variability found in the conformational features of the 10 base steps in B-DNA is mainly due to the flanking sequences or to interactions with the environment. From an analysis of the twist parameter of the base-pair steps available from crystals of oligonucleotides and protein/oligonucleotide complexes, we conclude that in most cases the flanking sequences show little influence: the conformation of a DNA region results from the combination of the independent intrinsic features of each base step (average conformation and intrinsic variability), modulated by their interactions with the environment. Only in some cases (YR steps, in particular CG and CA/TG) does it appear that flanking sequences have an influence on the conformation of the central base step. The values obtained allow an approximation to the parameters expected for repetitive DNA sequences. In particular, it is found that poly[d(AG/CT)] should have a strongly alternating conformation, in agreement with recently reported oligonucleotide structures.     

Helix–coil transition in poly(γ-benzyl-L-glutamate) and poly-ε-carbobenzoxy-L-lysine

In this paper two points are considered: the methods of evaluating the helical content θ and the calculation of the parameters of the transition from experimental data and its interpretation. The parameter ΔH obtained is in good agreement with the calorimetric one and v is found to be independent of temperature and solvent and in agreement with the ordinarily accepted value for poly(γ-benzyl-L -glutamate). The different methods of estimating θ are discussed for both polypeptides.     

Bacterial chemotaxis toward a NAPL source within a pore-scale microfluidic chamber

Chemotaxis toward chemical pollutants provides a mechanism for bacteria to migrate to locations of high contamination, which may improve the effectiveness of bioremediation. A microfluidic device was designed to mimic the dissolution of an organic-phase contaminant from a single pore into a larger macropore representing a preferred pathway for microorganisms that are carried along by groundwater flow. The glass windows of the microfluidic device allowed direct image analysis of bacterial distributions within the vicinity of the organic contaminant. Concentrations of chemotactic bacteria P. putida F1 near the organic/aqueous interface were 25% greater than those of a nonchemotactic mutant in the vicinity of toluene for a fluid velocity of 0.5 m/d. For E. coli responding to phenol, the bacterial concentrations were 60% greater than the controls, also at a velocity of 0.5 m/d. Velocities in the macropore were varied over a range from 0.5 to 10 m/d, the lower end of which is typical of groundwater velocities. The accumulation of chemotactic bacteria near the NAPL chemoattractant source decreased as the fluid velocity increased. Good agreement between computer-based simulations, generated using reasonable values of the model parameters, and the experimental data for P. putida strains confirmed the contribution due to chemotaxis. The experimental data for E. coli required a larger chemotactic sensitivity coefficient than that for P. putida, which was consistent with parameter values reported in the literature.     

Computer modeling of actinomycin D interactions with double-helical DNA

We have performed molecular mechanical calculations on intercalation complexes of actinomycin D with a series of base-paired hexanucleoside pentaphosphates; d(GCGCGC)2, d(GCCGGC)2, d(GCATGC)2, d(GCTAGC)2 and d(ATGCAT)2. Our results are in good agreement with previous experimental work on sequence selectivity. The results provide a rationalization for the strong preference of actinomycin D to intercalate on the 3' side of guanine residues, consistent with previously proposed models. Finally, the computed structures for d(ATGCAT)2-actinomycin D complexes have been compared with two-dimensional nuclear magnetic resonance nuclear Overhauser effect experimental results. To our knowledge, this is the first extensive comparison of molecular mechanical model structures for a drug-DNA complex with experimental solution phase data. We find generally good agreement between our computational models and the experimental solution phase structures.     

SBaddon: high performance simulation for the Systems Biology Toolbox for MATLAB

We present the SBaddon package as an extension to the Systems Biology Toolbox for MATLAB (SBtoolbox). The goal of this extension is to provide the users of the SBtoolbox with important functionality that is needed for parameter estimation applications. While simulation in the SBtoolbox relies on the MATLAB ODE solvers, the SBaddon package provides considerably increased simulation performance through automatic generation of compiled simulation functions. Furthermore, the package contains improved optimization algorithms, forward parameter sensitivity analysis and basic numeric parameter identifiability analysis. Availability: The SBaddon package is open source and freely available for non-commercial use. Commercial use of SBaddon is only possible through a specific licensing agreement (contact sbaddon@sbtoolbox.org). SBaddon can be obtained from http://www.sbtoolbox.org/SBaddon. The website also contains extensive documentation, and examples.     

Computational and analytical modeling of cationic lipid-DNA complexes

We present a theoretical study of the physical properties of cationic lipid-DNA (CL-DNA) complexes--a promising synthetically based nonviral carrier of DNA for gene therapy. The study is based on a coarse-grained molecular model, which is used in Monte Carlo simulations of mesoscopically large systems over timescales long enough to address experimental reality. In the present work, we focus on the statistical-mechanical behavior of lamellar complexes, which in Monte Carlo simulations self-assemble spontaneously from a disordered random initial state. We measure the DNA-interaxial spacing, d(DNA), and the local cationic area charge density, sigma(M), for a wide range of values of the parameter (c) representing the fraction of cationic lipids. For weakly charged complexes (low values of (c)), we find that d(DNA) has a linear dependence on (c)(-1), which is in excellent agreement with x-ray diffraction experimental data. We also observe, in qualitative agreement with previous Poisson-Boltzmann calculations of the system, large fluctuations in the local area charge density with a pronounced minimum of sigma(M) halfway between adjacent DNA molecules. For highly-charged complexes (large (c)), we find moderate charge density fluctuations and observe deviations from linear dependence of d(DNA) on (c)(-1). This last result, together with other findings such as the decrease in the effective stretching modulus of the complex and the increased rate at which pores are formed in the complex membranes, are indicative of the gradual loss of mechanical stability of the complex, which occurs when (c) becomes large. We suggest that this may be the origin of the recently observed enhanced transfection efficiency of lamellar CL-DNA complexes at high charge densities, because the completion of the transfection process requires the disassembly of the complex and the release of the DNA into the cytoplasm. Some of the structural properties of the system are also predicted by a continuum free energy minimization. The analysis, which semiquantitatively agrees with the computational results, shows that that mesoscale physical behavior of CL-DNA complexes is governed by interplay among electrostatic, elastic, and mixing free energies.     

Three-Dimensional Flow of an Oldroyd-B Nanofluid towards Stretching Surface with Heat Generation/Absorption

This article addresses the steady three-dimensional flow of an Oldroyd-B nanofluid over a bidirectional stretching surface with heat generation/absorption effects. Suitable similarity transformations are employed to reduce the governing partial differential equations into coupled nonlinear ordinary differential equations. These nonlinear ordinary differential equations are then solved analytically by using the homotpy analysis method (HAM). Graphically results are presented and discussed for various parameters, namely, Deborah numbers and , heat generation/absorption parameter Prandtl parameter , Brownian motion parameters, thermophoresis parameter and Lewis number . We have seen that the increasing values of the Brownian motion parameter and thermophoresis parameter leads to an increase in the temperature field and thermal boundary layer thickness while the opposite behavior is observed for concentration field and concentration boundary layer thickness. To see the validity of the present work, the numerical results are compared with the analytical solutions obtained by Homotopy analysis method and noted an excellent agreement for the limiting cases.