Exact calculations of survival probability for diffusion on growing lines, disks, and spheres: The role of dimension

Unlike standard applications of transport theory, the transport of molecules and cells during embryonic development often takes place within growing multidimensional tissues. In this work, we consider a model of diffusion on uniformly growing lines, disks, and spheres. An exact solution of the partial differential equation governing the diffusion of a population of individuals on the growing domain is derived. Using this solution, we study the survival probability, S(t). For the standard nongrowing case with an absorbing boundary, we observe that S(t) decays to zero in the long time limit. In contrast, when the domain grows linearly or exponentially with time, we show that S(t) decays to a constant, positive value, indicating that a proportion of the diffusing substance remains on the growing domain indefinitely. Comparing S(t) for diffusion on lines, disks, and spheres indicates that there are minimal differences in S(t) in the limit of zero growth and minimal differences in S(t) in the limit of fast growth. In contrast, for intermediate growth rates, we observe modest differences in S(t) between different geometries. These differences can be quantified by evaluating the exact expressions derived and presented here.

Calculations of helium separation via uniform pores of stanene-based membranes

The development of low energy cost membranes to separate He from noble gas mixtures is highly desired. In this work, we studied He purification using recently experimentally realized, two-dimensional stanene (2D Sn) and decorated 2D Sn (SnH and SnF) honeycomb lattices by density functional theory calculations. To increase the permeability of noble gases through pristine 2D Sn at room temperature (298 K), two practical strategies (i.e., the application of strain and functionalization) are proposed. With their high concentration of large pores, 2D Sn-based membrane materials demonstrate excellent helium purification and can serve as a superior membrane over traditionally used, porous materials. In addition, the separation performance of these 2D Sn-based membrane materials can be significantly tuned by application of strain to optimize the He purification properties by taking both diffusion and selectivity into account. Our results are the first calculations of He separation in a defect-free honeycomb lattice, highlighting new interesting materials for helium separation for future experimental validation.

Field measurement of beef pen manure methane and nitrous oxide reveals a surprise for inventory calculations

Few data exist on direct greenhouse gas emissions from pen manure at beef feedlots. However, emission inventories attempt to account for these emissions. This study used a large chamber to isolate N2O and CH4 emissions from pen manure at two Australian commercial beef feedlots (stocking densities, 13-27 m(2) head) and related these emissions to a range of potential emission control factors, including masses and concentrations of volatile solids, NO3-, total N, NH4+, and organic C (OC), and additional factors such as total manure mass, cattle numbers, manure pack depth and density, temperature, and moisture content. Mean measured pen N2O emissions were 0.428 kg ha(-1) d(-1) (95% confidence interval [CI], 0.252-0.691) and 0.00405 kg ha(-1) d(-1) (95% CI, 0.00114-0.0110) for the northern and southern feedlots, respectively. Mean measured CH4 emission was 0.236 kg ha(-1) d(-1) (95% CI, 0.163-0.332) for the northern feedlot and 3.93 kg ha(-1) d(-1) (95% CI, 2.58-5.81) for the southern feedlot. Nitrous oxide emission increased with density, pH, temperature, and manure mass, whereas negative relationships were evident with moisture and OC. Strong relationships were not evident between N2O emission and masses or concentrations of NO3- or total N in the manure. This is significant because many standard inventory calculation protocols predict N2O emissions using the mass of N excreted by the animal.

A new approach to thermochemical calculations of condensed fuel-oxidizer mixtures

A simple method of calculating the elemental stoichiometric coefficient, φe has been developed, which can easily be applied to multicomponent fuel-oxidizer compositions. The method correctly predicts whether a mixture is fuel lean, fuel rich, or stoichiometrically balanced. The total composition of oxidizing (or reducing) elements of the mixture appears to be related to the thermochemistry of the system. For the reaction of ammonium perchlorate and an organic fuel the heat of reaction varies linearly with the total composition of oxidizing elements. The physical significance of such a correlation based on thermochemical reasoning is highlighted in the paper.

The nonplanar peptide unit III. Quantum chemical calculations for related compounds and experimental X-ray diffraction data

The possible nonplanar distortions of the amide group in formamide, acetamide, N-methylacetamide, and N-ethylacetamide have been examined using CNDO/2 and INDO methods. The predictions from these methods are compared with the results obtained from X-ray and neutron diffraction studies on crystals of small open peptides, cyclic peptides, and amides. It is shown that the INDO results are in good agreement with observations, and that the dihedral angles N and defining the nonplanarity of the amide unit are correlated approximately by the relation N = -2, while C is small and uncorrelated with . The present study indicates that the nonplanar distortions at the nitrogen atom of the peptide unit may have to be taken into consideration, in addition to the variation in the dihedral angles (,), in working out polypeptide and protein structures.

Application of ab initio molecular orbital calculations to the structural moieties of carbohydrates. 3

Ab initio RHF/4-31G level molecular orbital calculations have been carried out on dimethoxymethane as a model compound for the acetal moiety in methyl pyranosides. The calculations are consistent with the predictions of the anomeric effect and the exo-anomeric effect. They reproduce very successfully the differences in molecular geometry observed by x-ray and neutron diffraction of single crystals of the methyl cy-D- and methyl 0-D-pyranosides. Calculations carried out at the 6-3 1G* level for methanediol confirm the earlier calculations at the 4-31G level, with smaller energy differences between the four staggered conformations.

Dissipative dynamics of a harmonic oscillator: A nonperturbative approach

Starting from a microscopic theory, we derive a master equation for a harmonic oscillator coupled to a bath of noninteracting oscillators. We follow a nonperturbative approach, proposed earlier by us for the free Brownian particle. The diffusion constants are calculated analytically and the positivity of the master equation is shown to hold above a critical temperature. We compare the long time behavior of the average kinetic and potential energies with known thermodynamic results. In the limit of vanishing oscillator frequency of the system, we recover the results of the free Brownian particle.

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.

Synthesis and Characterization of a Class of Donor-Acceptor Conjugated Molecules: Experiments and Theoretical Calculations

A class of conjugated molecules containing donor (thiophene) and acceptor (malononitrile) is synthesized by Knoevenagel condensation reaction between 2-(2,6-dimethy1-4H-pyran-4-ylidene) malononitrile and thiophene carbaldehyde containing two and three thiophene units. The resulting molecules are characterized by H-1 and C-13 NMR. We have performed UV-vis absorption, fluorescence, and cyclic voltammetry measurements on these materials. The spectroscopic and electrochemical measurements proved beyond doubt that these materials possess lowexcitation gap and are suitable for being an active material in various electronic devices. We have also performed electronic structure calculations using density functional theory (DFT) and INDO/SCI methods to characterize the ground and excited states of this class of molecules. These donor-acceptor molecules show a strong charge transfercharacter that increases with the increase in the number of thiophene rings coupled to the malononitrile acceptor moiety. We have also calculated the pi-coherence length, Stoke's shift, and effect of solvents on excited states for this class of molecules, Our theoretical values agree well with experimental results.