Non Equilibrium Molecular Dynamics Simulation of Coupled Heat- and Mass Transport Across a Liquid/Vapor Interface

Abstract

The transport properties of bulk liquid, gas and at the gas/liquid interface were studied for two binary Lennard-Jones/spline mixtures by use of nonequilibrium molecular dynamics. One of the mixtures was an ideal isotope mixture, the other a non-ideal mixture. The simulations gave the thermal conductivity, mutual diffusion coefficient, heat flux, mass flux, and the changes in these quantities across the interface. The local entropy production was expressed in terms of fluxes and thermodynamic forces, and numercial estimates are given. It was shown that the largest contribution to the total entropy production occurs in the vapor phase under the chosen conditions. We expect, however that if the mass flux were larger, the major contribution to the entropy production would come from the liquid phase.     

Marker-based adjustment of the additive relationship matrix for estimation of genetic parameters��an example using Eucalyptus cladocalyx

The effects of adjusting additive (numerator) relationship matrices (A) for inbreeding estimates taken from molecular markers were investigated using a small, model population of Eucalyptus cladocalyx. A number of individual-tree, mixed-models were compared, incorporating estimates of population- and family-level selfing and ancestral inbreeding applied either as average values to the entire population or as variable estimates for subpopulation and family groups. The consequences of ignoring inbreeding were inflated additive genetic variance estimates and underestimation of residual variance, with resulting inflation of heritability. We found models that correct for differential inbreeding at the subpopulation level give similar results to more complex ones including family-level estimates. Our analysis indicates that the commonly applied coefficient of relationship for first-generation eucalypt progeny of ρ = 1/2.5 appears to be quite suitable for correcting variance component and heritability estimates. However, if inbreeding is not specifically corrected for by adjustment of A, some minor rank changes of individual breeding values can occur, especially where levels of inbreeding vary among families, and some suboptimal selections and loss of genetic gain may ensue.     

Effects of coating roughness and biofouling on ship resistance and powering

Abstract

Predictions of full-scale ship resistance and powering are made for antifouling coating systems with a range of roughness and fouling conditions. The estimates are based on results from laboratory-scale drag measurements and boundary layer similarity law analysis. In the present work, predictions are made for a mid-sized naval surface combatant at cruising speed and near maximum speed. The results indicate that slime films can lead to significant increases in resistance and powering, and heavy calcareous fouling results in powering penalties up to 86% at cruising speed. The present estimates show good agreement with results from full-scale ship power trials.     

Above-ground biomass references for urban trees from terrestrial laser scanning data

Background and AimsWithin extending urban areas, trees serve a multitude of functions (e.g. carbon storage, suppression of air pollution, mitigation of the ‘heat island’ effect, oxygen, shade and recreation). Many of these services are positively correlated with tree size and structure. The quantification of above-ground biomass (AGB) is of especial importance to assess its carbon storage potential. However, quantification of AGB is difficult and the allometries applied are often based on forest trees, which are subject to very different growing conditions, competition and form. In this article we highlight the potential of terrestrial laser scanning (TLS) techniques to extract highly detailed information on urban tree structure and AGB.MethodsFifty-five urban trees distributed over seven cities in Switzerland were measured using TLS and traditional forest inventory techniques before they were felled and weighed. Tree structure, volume and AGB from the TLS point clouds were extracted using quantitative structure modelling. TLS-derived AGB estimates were compared with AGB estimates based on forest tree allometries dependent on diameter at breast height only. The correlations of various tree metrics as AGB predictors were assessed.Key ResultsEstimates of AGB derived by TLS showed good performance when compared with destructively harvested references, with an R² of 0.954 (RMSE = 556 kg) compared with 0.837 (RMSE = 1159 kg) for allometrically derived AGB estimates. A correlation analysis showed that different TLS-derived wood volume estimates as well as trunk diameters and tree crown metrics show high correlation in describing total wood AGB, outperforming tree height.ConclusionsWood volume estimates based on TLS show high potential to estimate tree AGB independent of tree species, size and form. This allows us to retrieve highly accurate non-destructive AGB estimates that could be used to establish new allometric equations without the need for extensive destructive harvesting.     

Bayesian inference for parameter estimation in lactoferrin-mediated iron transport across blood-brain barrier

BackgroundIn neurodegenerative diseases such as Alzheimer's and Parkinson's, excessive irons as well as lactoferrin (Lf), but not transferrin (Tf), have been found in and around the affected regions of the brain. These evidences suggest that lactoferrin plays a critical role during neurodegenerative diseases, although Lf-mediated iron transport across blood-brain barrier (BBB) is negligible compared to that of transferrin in normal condition. However, the kinetics of lactoferrins and lactoferrin-mediated iron transport are still unknown.MethodTo determine the kinetic rate constants of lactoferrin-mediated iron transport through BBB, a mass-action based ordinary differential equation model has been presented. A Bayesian framework is developed to estimate the kinetic rate parameters from posterior probability density functions. The iron transport across BBB is studied by considering both Lf- and Tf-mediated pathways for both normal and pathologic conditions.ResultsUsing the point estimates of kinetic parameters, our model can effectively reproduce the experimental data of iron transport through BBB endothelial cells. The robustness of the model and parameter estimation process are further verified by perturbation of kinetic parameters. Our results show that surge in high-affinity receptor density increases lactoferrin as well as iron in the brain.ConclusionsDue to the lack of a feedback loop such as iron regulatory proteins (IRPs) for lactoferrin, iron can transport to the brain continuously, which might increase brain iron to pathological levels and can contribute to neurodegeneration.General significanceThis study provides an improved understanding of presence of lactoferrin and iron in the brain during neurodegenerative diseases.     

Molecular Dynamics as a Mathematical Mapping. I. Differentiable Force Functions

Abstract

The molecular dynamics technique can be viewed as a deterministic mathematical mapping between, on one side, the force field parameters that describe the potential energy interactions and the input macroscopic conditions, and, on the other, the calculated macroscopic properties of the bulk molecular system.

The differentiability of such a mapping in the conventional molecular dynamics calculations is affected by the discontinuities in particle positions introduced by the periodic boundary conditions and the discontinuities introduced by the minimum image convention and other methods commonly employed to approximate the calculation of interparticle potential and force.

This paper proposes an alternative molecular dynamics framework based on modified force functions which are almost everywhere continuous and differentiable, and exhibit a natural periodicity. These characteristics obviate the need for both the periodic boundary conditions and the minimum image convention, as well as for any corrections for long-range interactions. They also make it possible to apply standard methods of variational calculus for the computation of partial derivatives of the molecular dynamics mapping.

The modified framework is first introduced for the case of simple monoatomic fluids where the nature of the forces exerted between any pair of two particles is identical. A more general model describing the interactions of flexible molecules is then developed. We describe the application of this approach to mixtures of alkane molecules interacting via the NERD force field.     

Nucleocytoplasmic transport under stress conditions and its role in HSP70 chaperone systems

BackgroundIn eukaryotic cells, molecular trafficking between the nucleus and cytoplasm is a highly regulated process related to cellular homeostasis and cellular signaling. However, various cellular stresses induce the perturbation of conventional nucleocytoplasmic transport pathways, resulting in the nucleocytoplasmic redistribution of many functional proteins.Scope of reviewWe describe the recent insights into the mechanism and functions of nuclear import of cytosolic chaperone HSP70 under stress conditions and the cellular distribution and functions of its co-chaperones.Major conclusionsHikeshi mediates the nuclear import of the molecular chaperone HSP70. A few of the regulators of the HSP70 chaperone system also accumulate in the nucleus under heat stress conditions. These proteins function collaboratively to protect cells from stress-induced damage and aid in the recovery of cells from stress.General significanceStudies on the regulation of nucleocytoplasmic transport under several cellular stresses should provide new insights into the fundamental principles of protein homeostasis (proteostasis) in both compartments, the nucleus and cytoplasm.     

Robustness of individual and marginal model-based estimates: A sensitivity analysis of flexible parametric models

BackgroundFlexible parametric survival models (FPMs) are commonly used in epidemiology. These are preferred as a wide range of hazard shapes can be captured using splines to model the log-cumulative hazard function and can include time-dependent effects for more flexibility. An important issue is the number of knots used for splines. The reliability of estimates are assessed using English data for 10 cancer types and the use of online interactive graphs to enable a more comprehensive sensitivity analysis at the control of the user is demonstrated.MethodsSixty FPMs were fitted to each cancer type with varying degrees of freedom to model the baseline excess hazard and the main and time-dependent effect of age. For each model, we obtained age-specific, age-group and internally age-standardised relative survival estimates. The Akaike Information Criterion and Bayesian Information Criterion were also calculated and comparative estimates were obtained using the Ederer II and Pohar Perme methods. Web-based interactive graphs were developed to present results.ResultsAge-standardised estimates were very insensitive to the exact number of knots for the splines. Age-group survival is also stable with negligible differences between models. Age-specific estimates are less stable especially for the youngest and oldest patients, of whom there are very few, but for most scenarios perform well.ConclusionAlthough estimates do not depend heavily on the number of knots, too few knots should be avoided, as they can result in a poor fit. Interactive graphs engage researchers in assessing model sensitivity to a wide range of scenarios and their use is highly encouraged.     

Multi-scale mechanical characterization of prostate cancer cell lines: Relevant biological markers to evaluate the cell metastatic potential

BackgroundConsidering the importance of cellular mechanics in the birth and evolution of cancer towards increasingly aggressive stages, we compared nano-mechanical properties of non-tumoral (WPMY-1) and highly aggressive metastatic (PC-3) prostate cell lines both on cell aggregates, single cells, and membrane lipids.MethodsCell aggregate rheological properties were analyzed during dynamic compression stress performed on a homemade rheometer. Single cell visco-elasticity measurements were performed by Atomic Force Microscopy using a cantilever with round tip on surface-attached cells. At a molecular level, the lateral diffusion coefficient of total extracted lipids deposited as a Langmuir monolayer on an air-water interface was measured by the FRAP technique.ResultsAt cellular pellet scale, and at single cell scale, PC-3 cells were less stiff, less viscous, and thus more prone to deformation than the WPMY-1 control. Interestingly, stress-relaxation curves indicated a two-step response, which we attributed to a differential response coming from two cell elements, successively stressed. Both responses are faster for PC-3 cells. At a molecular scale, the dynamics of the PC-3 lipid extracts are also faster than that of WPMY-1 lipid extracts.ConclusionsAs the evolution of cancer towards increasingly aggressive stages is accompanied by alterations both in membrane composition and in cytoskeleton dynamical properties, we attribute differences in viscoelasticity between PC-3 and WPMY-1 cells to modifications of both elements.General significanceA decrease in stiffness and a less viscous behavior may be one of the diverse mechanisms that cancer cells adopt to cope with the various physiological conditions that they encounter.     

The efficacy of respondent-driven sampling for the health assessment of minority populations

BackgroundRespondent driven sampling (RDS) is a relatively new network sampling technique typically employed for hard-to-reach populations. Like snowball sampling, initial respondents or “seeds” recruit additional respondents from their network of friends. Under certain assumptions, the method promises to produce a sample independent from the biases that may have been introduced by the non-random choice of “seeds.” We conducted a survey on health communication in Guam’s general population using the RDS method, the first survey that has utilized this methodology in Guam. It was conducted in hopes of identifying a cost-efficient non-probability sampling strategy that could generate reasonable population estimates for both minority and general populations.MethodsRDS data was collected in Guam in 2013 (n = 511) and population estimates were compared with 2012 BRFSS data (n = 2031) and the 2010 census data. The estimates were calculated using the unweighted RDS sample and the weighted sample using RDS inference methods and compared with known population characteristics.ResultsThe sample size was reached in 23 days, providing evidence that the RDS method is a viable, cost-effective data collection method, which can provide reasonable population estimates. However, the results also suggest that the RDS inference methods used to reduce bias, based on self-reported estimates of network sizes, may not always work. Caution is needed when interpreting RDS study findings.ConclusionsFor a more diverse sample, data collection should not be conducted in just one location. Fewer questions about network estimates should be asked, and more careful consideration should be given to the kind of incentives offered to participants.     

Nonlinear Regression Analysis of the time Course of Ligand Binding Experiments

Abstract

The binding and release of hormones and growth factors are often relatively slow processes under biological conditions. Consequently, a knowledge of the underlying rate constants may be of greater physiological relevance than the equilibrium constant. Here we show how, by following a single time course of binding, the rate constants for both binding and release can be determined. The ratio of these rate constants allows the binding constant to be calculated. A nonlinear regression computer program is described which facilitates these calculations and which provides estimates and standard errors of the constants determined. The method is illustrated by the binding of human growth hormone to the human growth hormone binding protein, and the binding of ovine prolactin to the rabbit prolactin receptor.     

Cross-genera amplification of Cajanus spp. specific SSR markers in Clitoria ternatea (L.) and their application in genetic diversity studies

Clitoria ternatea (L.) is a medicinal leguminous plant and is cultivated to cater the need of herbal industries and asthetic purposes. The unavailability of steady molecular marker impedes the genetic improvement of C. ternatea. In the present study, transferability of 98 pairs of Cajanus spp. specific SSR primers were assessed among 14 genotypes of C. ternatea, varied for their flower color, floral architecture and bio-metabolite (taraxerol and delphinidin) content, and out of them 43 had successfully amplified the fragments. Among them, 36 pairs of primers showed 100% transferability, whereas rest seven varied from 42.86 to 92.85% transferability. The transferable 43 pairs of SSR primers generated 196 alleles across the 14 genotypes and the AMOVA analysis showed moderate genetic variation (55.1%) among the genotypes of C. ternatea, which was also reinforced by Nei’s genetic distance and gene identity estimates derived haplotype matrix. Similarly, both the principal coordinate analysis and dendrogram grouped these 14 genotypes of C. ternatea into two major clusters based on SSR allele distribution and frequency, and the clustering pattern is in accordance with petal color but in contrast to floral architecture. MCheza based outlier analysis revealed 16 alleles for balancing selection, which are putatively involved in the maintenance of genetic polymorphism in C. ternatea. Moreover, the estimates of molecular diversity and bio-metabolite content revealed the possible use of these genotypes in future breeding programme of this species.Electronic supplementary materialThe online version of this article (10.1007/s12298-020-00907-x) contains supplementary material, which is available to authorized users.     

The organization and repair of DNA in the mammalian chromosome. I. Calibration procedures and errors in the determination of the molecular weight of a native DNA

The semiautomated sucrose gradient system of Lange and Liberman (1974) (Anal. Biochem. 59 , 129–145) has been used to determine an accurate value of the exponent k of the Burgi and Hershey equation (1963) (Biophys. J. 3 , 309–321) under conditions of 1M salt. A complete analysis of the variance of k was done for data from both systems. The results lead to the conclusion that the value k = 0.401 (±0.012 for reproducibility; ±0.019 for total error) obtained here is considerably more accurate than that of Burgi and Hershey (even when corrected for DNA concentration effects and the improved size estimates of T₄ DNA size). The use of the Burgi and Hershey relationship and a T₄ DNA metric to determine molecular weights of mammalian DNA requires the utmost precision in the exponent if useful estimates are to be obtained. This also has implications for double-strand breakage and repair estimates.     

Estimating the proportion cured of cancer: Some practical advice for users

BackgroundCure models can provide improved possibilities for inference if used appropriately, but there is potential for misleading results if care is not taken. In this study, we compared five commonly used approaches for modelling cure in a relative survival framework and provide some practical advice on the use of these approaches.Patients and methodsData for colon, female breast, and ovarian cancers were used to illustrate these approaches. The proportion cured was estimated for each of these three cancers within each of three age groups. We then graphically assessed the assumption of cure and the model fit, by comparing the predicted relative survival from the cure models to empirical life table estimates.ResultsWhere both cure and distributional assumptions are appropriate (e.g., for colon or ovarian cancer patients aged <75 years), all five approaches led to similar estimates of the proportion cured. The estimates varied slightly when cure was a reasonable assumption but the distributional assumption was not (e.g., for colon cancer patients ≥75 years). Greater variability in the estimates was observed when the cure assumption was not supported by the data (breast cancer).ConclusionsIf the data suggest cure is not a reasonable assumption then we advise against fitting cure models. In the scenarios where cure was reasonable, we found that flexible parametric cure models performed at least as well, or better, than the other modelling approaches. We recommend that, regardless of the model used, the underlying assumptions for cure and model fit should always be graphically assessed.     

Breast and colorectal cancer recurrence-free survival estimates in the US: Modeling versus active data collection

BackgroundA modeling method was developed to estimate recurrence-free survival using cancer registry survival data. This study aims to validate the modeled recurrence-free survival against “gold-standard” estimates from data collected by the National Program of Cancer Registries (NPCR) Patient-Centered Outcomes Research (PCOR) project.MethodsWe compared 5-year metastatic recurrence-free survival using modeling and empirical estimates from the PCOR project that collected disease-free status, tumor progression and recurrence for colorectal and female breast cancer cases diagnosed in 2011 in 5 U.S. state registries. To estimate empirical recurrence-free survival, we developed an algorithm that combined disease-free, recurrence, progression, and date information from NPCR-PCOR data. We applied the modeling method to relative survival for patients diagnosed with female breast and colorectal cancer in 2000–2015 in the SEER-18 areas.ResultsWhen grouping patients with stages I-III, the 5-year metastatic recurrence-free modeled and NPCR-PCOR estimates are very similar being respectively, 90.2 % and 88.6 % for female breast cancer, 74.6 % and 75.3 % for colon cancer, and 68.8 % and 68.5 % for rectum cancer. In general, the 5-year recurrence-free NPCR-PCOR and modeled estimates are still similar when controlling by stage. The modeled estimates, however, are not as accurate for recurrence-free survival in years 1–3 from diagnosis.ConclusionsThe alignment between NPCR-PCOR and modeled estimates supports their validity and provides robust population-based estimates of 5-year metastatic recurrence-free survival for female breast, colon, and rectum cancers. The modeling approach can in principle be extended to other cancer sites to provide provisional population-based estimates of 5-year recurrence free survival.     

Population differences in vision acuity: A review,with speculative notes on selection relaxation

Abstract

The contrast between the lower frequencies of vision defects reported among several different samples of hunting populations and the higher frequencies found repeatedly among populations with relatively long histories of agriculture and settled habitations points strongly to the relaxation of natural selection. Two rough estimates of the heritability of that conglomeration of genetic entities comprising “vision defects” are found in studies of MZ versus DZ twins and in studies of offspring of consanguineous parents versus those of unrelated parents. The hypothesis of relaxation of natural selection is plausible as man developed culture habitats more tolerant of persons with slight refractive aberrations. One can reasonably imagine that hereditary deficiencies accumulated sufficiently under such conditions to make possible quantitative estimates of differential population frequencies.     

Model Testing in Radioligand/Receptor Interaction by Monte Carlo Simulation

Abstract

Monte Carlo simulations have been applied for evaluating the reliability of parameter estimates as well as for testing models in radioligand saturation binding experiments. Scatchard analysis was compared to the nonlinear least-square curve fitting method for one-site saturation binding curves. It was found that linear regression analysis from the transformed data in the Scatchard plot yielded generally less accurate parameter estimates than nonlinear regression analysis of untransformed data. The advantage of the nonlinear least-squares curve fitting method was especially pronounced in cases where the scatter and number of data points, as well as the radioligand concentration range, were chosen similar to less optimal experimental conditions. Under such circumstances, several K_D and B_max values derived by Scatchard analysis led to physically impossible negative values whereas the same data analyzed by nonlinear regression yielded reasonable parameter estimates. Furthermore, it was found that for both means of analysis, K_D and B_max correlated positively. In another set of Monte Carlo experiments, saturation binding curves involving two receptor sites were generated and subsequently analyzed according to both a one-site and a two-site model. The confidence with which one is able to distinguish the two-site model from nonlinear least-squares curve fitting was then estimated for optimal, as well as for, less ideal experimental condigions.     

On the re-calibration process in radiochromic film dosimetry

PurposeThe accuracy and precision of the dose estimates obtained with radiochromic film dosimetry are investigated in a clinical environment. The improvement in the accuracy of dose estimates reached with corrective methods is analyzed. Two novel re-calibration algorithms for radiochromic film dosimetry are presented.MethodsTwo different EBT3 lots are evaluated in two different centres. They are calibrated in Varian linacs and read in two different EPSON scaners. Once the lots are calibrated, three films per lot are considered and divided into stripes that are exposed to known doses. Several dosimetry protocols usually employed in radiochromic film dosimetry are used to convert film responses to absorbed doses. These protocols are characterized by different choices of the film responses or different sensitometric curves. Finally, the accuracy and reproducibility of the dose estimates is investigated with and without the corrective methods.Results and ConclusionsThe variabilities that affect radiochromic film dosimetry, such as intra-lot variability, inter-scan variability, post-exposure time and film autodevelopment may give rise to inaccuracies in the dose estimates. However, the implementation of re-calibration methods leads to more accurate dose estimates. All the investigated protocols showed more accurate and reproducible results when the re-calibrated methods were employed. So, the novel re-calibration methods may be applied in order to improve the accuracy and reproducibility of radiochromic film dosimetry.     

CRISPR/Cas9编辑大肠杆菌工程菌生产氨基葡萄糖

【背景】氨基葡萄糖(glucosamine, GlcN)及其衍生物N-乙酰氨基葡萄糖(N-acetylglucosamine,GlcNAc)是合成糖胺聚糖的重要前体物质,在医药、化妆品和保健品领域具有广泛的应用价值。传统的生产方式存在诸多弊端,如环境污染、原料限制、不适于海鲜易过敏人群等问题,因此利用微生物发酵法生产GlcN和GlcNAc越来越受到青睐。【目的】利用微生物发酵生产并提高N-乙酰氨基葡萄糖的产量,探索分子改造及发酵条件优化策略。【方法】以大肠杆菌MG1655为出发菌株,首先利用表达载体共表达大肠杆菌来源的glmS和酿酒酵母来源的gna1,构建GlcNAc的生物合成路径,然后利用CRISPR/Cas9技术敲除GlcNAc的分解代谢与转运途径,以提高GlcNAc的产量,最后结合发酵条件优化使GlcNAc的产量得到进一步提升。【结果】通过分子改造得到一株产GlcNAc菌株RY-5,发酵20 h后GlcNAc的产量达到了2.36 g/L,相较于初始构建的菌株RY-1提高了29倍,进一步对装液量和诱导剂IPTG的添加时间等条件进行发酵优化,GlcNAc产量达到了7.74g/L,与优...