A new set of relativistic screening constants for the screened hydrogenic model

AnewRelativisticScreenedHydrogenicModel has been developed to calculate atomic data needed to compute the optical and thermodynamic properties of high energy density plasmas. The model is based on anewset of universal screeningconstants, including nlj-splitting that has been obtained by fitting to a large database of ionization potentials and excitation energies. This database was built with energies compiled from the National Institute of Standards and Technology (NIST) database of experimental atomic energy levels, and energies calculated with the Flexible Atomic Code (FAC). The screeningconstants have been computed up to the 5p3/2 subshell using a Genetic Algorithm technique with an objective function designed to minimize both the relative error and the maximum error. To select the best set of screeningconstants some additional physical criteria has been applied, which are based on the reproduction of the filling order of the shells and on obtaining the best ground state configuration. A statistical error analysis has been performed to test the model, which indicated that approximately 88% of the data lie within a ±10% error interval. We validate the model by comparing the results with ionization energies, transition energies, and wave functions computed using sophisticated self-consistent codes and experimental data.

Optimal fits of diffusion constants from single-time data points of Brownian trajectories

Experimental methods based on single particle tracking (SPT) are being increasingly employed in the physical and biological sciences, where nanoscale objects are visualized with high temporal and spatial resolution. SPT can probe interactions between a particle and its environment but the price to be paid is the absence of ensemble averaging and a consequent lack of statistics. Here we address the benchmark question of how to accurately extract the diffusion constant of one single Brownian trajectory. We analyze a class of estimators based on weighted functionals of the square displacement. For a certain choice of the weight function these functionals provide the true ensemble averaged diffusion coefficient, with a precision that increases with the trajectory resolution.

Determination of the JWL constants for ANFO and emulsion explosives from cylinder test data

The Jones-Wilkins-Lee (JWL) equation of state parameters for ANFO and emulsion-type explosives have been obtained from cylinder test expansion measurements. The calculation method comprises a new radial expansion function, with a non-zero initial velocity at the onset of the expansion in order to comply with a positive Gurney energy at unit relative volume, as the isentropic expansion from the CJ state predicts. The equations reflecting the CJ state conditions and the measured expansion energy were solved for the JWL parameters by a non-linear least squares scheme. The JWL parameters of thirteen ANFO and emulsion type explosives have been determined in this way from their cylinder test expansion data. The results were evaluated through numerical modelling of the tests with the LS-DYNA hydrocode; the expansion histories from the modelling were compared with the measured ones, and excellent agreement was found.

Application of the homogenization techniques to the determination of the effective flexure constants of plate structural systems

A Mindlin plate with periodically distributed ribs patterns is analyzed by using homogenization techniques based on asymptotic expansion methods. The stiffness matrix of the homogenized plate is found to be dependent on the geometrical characteristics of the periodical cell, i.e. its skewness, plan shape, thickness variation etc. and on the plate material elastic constants. The computation of this plate stiffness matrix is carried out by averaging over the cell domain some solutions of different periodical boundary value problems. These boundary value problems are defined in variational form by linear first order differential operators on the cell domain and the boundary conditions of the variational equation correspond to a periodic structural problem. The elements of the stiffness matrix of homogenized plate are obtained by linear combinations of the averaged solution functions of the above mentioned boundary value problems. Finally, an illustrative example of application of this homogenization technique to hollowed plates and plate structures with ribs patterns regularly arranged over its area is shown. The possibility of using in the profesional practice the present procedure to the actual analysis of floors of typical buildings is also emphasized.

Simultaneous determination of helical unwinding angles and intrinsic association constants in ligand–DNA complexes: The interaction between DNA and calichearubicin B

We present a helical unwinding assay for reversibly binding DNA ligands that uses closed circular DNA, topoisomerase I (Topo I), and two-dimensional agarose gel electrophoresis. Serially diluted Topo I relaxation reactions at constant DNA/ligand ratio are performed, and the resulting apparent unwinding of the closed circular DNA is used to calculate both ligand unwinding angle (φ) and intrinsic association constant (Ka). Mathematical treatment of apparent unwinding is formally analogous to that of apparent extinction coefficient data for optical binding titrations. Extrapolation to infinite DNA concentration yields the true unwinding angle of a given ligand and its association constant under Topo I relaxation conditions. Thus this assay delivers simultaneous structural and thermodynamic information describing the ligand–DNA complex. The utility of this assay has been demonstrated by using calichearubicin B (CRB), a synthetic hybrid molecule containing the anthraquinone chromophore of (DA) and the carbohydrate domain of calicheamicin γ1I. The unwinding angle for CRB calculated by this method is −5.3 ± 0.5°. Its Ka value is 0.20 × 106 M−1. For comparison, the unwinding angles of ethidium bromide and DA have been independently calculated, and the results are in agreement with canonical values for these compounds. Although a stronger binder to selected sites, CRB is a less potent unwinder than its parent compound DA. The assay requires only small amounts of ligand and offers an attractive option for analysis of DNA binding by synthetic and natural compounds.

Effective intercellular communication distances are determined by the relative time constants for cyto/chemokine secretion and diffusion

A cell’s ability to effectively communicate with a neighboring cell is essential for tissue function and ultimately for the organism to which it belongs. One important mode of intercellular communication is the release of soluble cyto- and chemokines. Once secreted, these signaling molecules diffuse through the surrounding medium and eventually bind to neighboring cell’s receptors whereby the signal is received. This mode of communication is governed both by physicochemical transport processes and cellular secretion rates, which in turn are determined by genetic and biochemical processes. The characteristics of transport processes have been known for some time, and information on the genetic and biochemical determinants of cellular function is rapidly growing. Simultaneous quantitative analysis of the two is required to systematically evaluate the nature and limitations of intercellular signaling. The present study uses a solitary cell model to estimate effective communication distances over which a single cell can meaningfully propagate a soluble signal. The analysis reveals that: (i) this process is governed by a single, key, dimensionless group that is a ratio of biological parameters and physicochemical determinants; (ii) this ratio has a maximal value; (iii) for realistic values of the parameters contained in this dimensionless group, it is estimated that the domain that a single cell can effectively communicate in is ≈250 μm in size; and (iv) the communication within this domain takes place in 10–30 minutes. These results have fundamental implications for interpretation of organ physiology and for engineering tissue function ex vivo.

Equilibrium constants and free energies in unfolding of proteins in urea solutions

A novel thermodynamic approach to the reversible unfolding of proteins in aqueous urea solutions has been developed based on the premise that urea ligands are bound cooperatively to the macromolecule. When successive stoichiometric binding constants have values larger than expected from statistical effects, an equation for moles of bound urea can be derived that contains imaginary terms. For a very steep unfolding curve, one can then show that the fraction of protein unfolded, B̄, depends on the square of the urea concentration, U, and is given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\bar {B}=\frac{{\mathit{A}}^{{\mathit{2}}}_{{\mathit{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}n\bar {B}}{\mathit{U}}^{{\mathit{2}}}}{{\mathrm{1\hspace{.167em}+\hspace{.167em}}}{\mathit{A}}^{{\mathrm{2}}}_{{\mathrm{1}}}{\mathit{e}}^{{\mathrm{{\lambda}}}\bar {B}}{\mathit{U}}^{{\mathrm{2}}}}{\mathrm{.}}\end{equation*}\end{document} A12 is the binding constant as B̄→ 0, and λ is a parameter that reflects the augmentation in affinities of protein for urea as the moles bound increases to the saturation number, n. This equation provides an analytic expression that reproduces the unfolding curve with good precision, suggests a simple linear graphical procedure for evaluating A12 and λ, and leads to the appropriate standard free energy changes. The calculated ΔG° values reflect the coupling of urea binding with unfolding of the protein. Some possible implications of this analysis to protein folding in vivo are described.

Determination of the binding constants of the centromere protein Cbf1 to all 16 centromere DNAs of Saccharomyces cerevisiae

Cbf1p is a Saccharomyces cerevisiae chromatin protein belonging to the basic region helix–loop–helix leucine zipper (bHLHzip) family of DNA binding proteins. Cbf1p binds to a conserved element in the 5′-flanking region of methionine biosynthetic genes and to centromere DNA element I (CDEI) of S.cerevisiae centromeric DNA. We have determined the apparent equilibrium dissociation constants of Cbf1p binding to all 16 CDEI DNAs in gel retardation assays. Binding constants of full-length Cbf1p vary between 1.7 and 3.8 nM. However, the dissociation constants of a Cbf1p deletion variant that has been shown to be fully sufficient for Cbf1p function in vivo vary in a range between 3.2 and 12 nM. In addition, native polyacrylamide gel electrophoresis revealed distinct changes in the 3D structure of the Cbf1p/CEN complexes. We also show that the previously reported DNA binding stimulation activity of the centromere protein p64 functions on both the Cbf1 full-length protein and a deletion variant containing only the bHLHzip domain of Cbf1p. Our results suggest that centromeric DNA outside the consensus CDEI sequence and interaction of Cbf1p with adjacent centromere proteins contribute to the complex formation between Cbf1p and CEN DNA.

On the nucleation and growth of amyloid beta-protein fibrils: detection of nuclei and quantitation of rate constants.

We have studied the fibrillogenesis of synthetic amyloid beta-protein-(1-40) fragment (A beta) in 0.1 M HCl. At low pH, A beta formed fibrils at a rate amenable to detailed monitoring by quasi-elastic light-scattering spectroscopy. Examination of the fibrils with circular dichroism spectroscopy and electron microscopy showed them to be highly similar to those found in amyloid plaques. We determined the hydrodynamic radii of A beta aggregates during the entire process of fibril nucleation and growth. Above an A beta concentration of approximately 0.1 mM, the initial rate of elongation and the final size of fibrils were independent of A beta concentration. Below an A beta concentration of 0.1 mM, the initial elongation rate was proportional to the peptide concentration, and the resulting fibrils were significantly longer than those formed at higher concentration. We also found that the surfactant n-dodecylhexaoxyethylene glycol monoether (C12E6) slowed nucleation and elongation of fibrils in a concentration-dependent manner. Our observations are consistent with a model of A beta fibrillogenesis that includes the following key steps: (i) peptide micelles form above a certain critical A beta concentration, (ii) fibrils nucleate within these micelles or on heterogeneous nuclei (seeds), and (iii) fibrils grow by irreversible binding of monomers to fibril ends. Interpretation of our data enabled us to determine the sizes of fibril nuclei and A beta micelles and the rates of fibril nucleation (from micelles) and fibril elongation. Our approach provides a powerful means for the quantitative assay of A beta fibrillogenesis.

Improvements in DFT Calculations of Spin–Spin Coupling Constants

Different types of spin–spin coupling constants (SSCCs) for several representative small molecules are evaluated and analyzed using a combination of 10 exchange functionals with 12 correlation functionals. For comparison, calculations performed using MCSCF, SOPPA, other common DFT methods, and also experimental data are considered. A detailed study of the percentage of Hartree–Fock exchange energy in SSCCs and in its four contributions is carried out. From the above analysis, a combined functional formed with local Slater (34%), Hartree–Fock exchange (66%), and P86 correlation functional (S66P86) is proposed in this paper. The accuracy of the values obtained with this hybrid functional (mean absolute deviation of 4.5 Hz) is similar to that of the SOPPA method (mean absolute deviation of 4.6 Hz).

Intramolecular surfaces for vicinal proton–proton coupling constants 3JHH

Equations for the intramolecular surfaces of the 3JHH coupling constants in ethane, ethylene, and acetylene are formulated, and the corresponding coefficients are estimated from calculations at the DFT/B3LYP level. The chosen variables are changes in bond lengths, in the torsion angle φ between the coupled protons Ha and Hb, in bond angles, and in dihedral angles. The 3JHH surface of ethane is formulated as an extended Karplus equation with the coefficients of a truncated Fourier series on the torsion angle φ expanded as second-order Taylor series in the chosen variables taking into account the invariance of 3JHH under reflections and rotations of nuclear coordinates. Partial vibrational contributions from linear and square terms corresponding to changes in the geometry of the Ha − Ca − Cb − Hb fragment are important while those from cross terms are small with a few exceptions. The 3JHH surface of ethane is useful to predict contributions to 3JHH from changes in local geometry of derivatives but vibrational contributions are predicted less satisfactorily. The predicted values at the B3LYP/BS2 level of the 3JHH couplings (vibrational contributions at 300 K) from equilibrium geometries are 9.79 (−0.17) for acetylene, and 17.08 (1.93) and 10.73(0.93) for the trans and cis couplings of ethylene.

Table of conversion factors and physical constants.

Tables.