DNA and estrogen receptor interaction revealed by fragment molecular orbital calculations

Molecular orbital calculations of the complex between DNA-ERE (estrogen response element) and ER (estrogen receptor)-DBD (DNA-binding domain) were performed using the fragment molecular orbital (FMO) method, which enables large-scale MO (molecular orbital) calculations by reducing the computational cost and by significantly increasing efficiency for parallel computation. Such a large system, which contains 3354 atoms, is impractical via conventional MO methods due to the immense computational cost. Details of the interaction between DNA-ERE and ER-DBD were revealed in this study as follows by using the FMO calculations to analyze the interfragment interaction energies (IFIEs) and the electrostatic potentials (ESPs). An area with a high positive ESP is identified on the DNA-binding side of ER-DBD and is the main driving force behind access to the DNA. The position of the ER-DBD monomer can be fixed on a phosphate group of DNA-ERE by the strong electrostatic interactions, whereas the rotation cannot be fixed. In contrast, both the position and rotation of the ER-DBD dimer can be fixed and can therefore form the stable (ER-DBD)2...DNA-ERE complex. Dimerization of the ER-DBD monomers, each of which have a charge of +5 , is mainly due to large attractive interaction energies of the second Zn fragments. The base pairs in the consensus sequence of DNA-ERE interact only with the recognition helix located in the major groove due to the large shielding effect of the phosphate groups of DNA. The recognition helix has weaker interactions with the base pairs than the electrostatic interactions with the phosphate groups. Thus, the DNA-binding machinery of the ER-DBD dimer, which can secure the recognition helix in the major groove of DNA, is crucial for interactions between the recognition helix and base pairs.     

The symmetry of two-photon excited states as determined by time-resolved fluorescence depolarization experiments

A new experimental and theoretical approach is presented for the quantitative determination and assignment of the two-photon absorption tensor of fluorophores dissolved in liquid solutions. Two linearly independent time-resolved fluorescence anisotropies and the two-photon polarization ratio were determined from experiments based on using the time-correlated single photon counting technique. The data were analyzed in a global manner under the assumption of prevailing diffusive molecular reorientations and when accounting for the influence of rapid unresolved reorientations. The method has been applied in fluorescence studies of perylene, two-photon excited at 800 nm. The analysis suggests that the two-photon transition is mediated via vibronic coupling including at least two vibrations of different symmetry, and also that the first singlet excited electronic state acts as a dominating intermediate state.     

Selective pressurized liquid extraction for multi-residue analysis of organochlorine pesticides in soil

A selective pressurized liquid extraction (SPLE) procedure capable of performing simultaneous extraction and clean-up has been demonstrated for multi-residue analysis of organochlorine pesticides (OCPs) in soil. The final method was performed at 100 degrees C for 3 x 10 min using acetone/n-heptane (1:1, v/v). Florisil was placed inside the extraction cell downstream the sample to remove interfering compounds. Extraction of two soil samples by SPLE gave a recovery which was over 80% for beta-endosulfan, endosulfan sulfate, p,p'-DDT and p,p'-DDE compared to PLE with off-line clean-up. The same trend was observed when applying the SPLE method to a certified reference soil sample (CRM 811-050) containing 13 OCPs, where the SPLE method gave 80-90% recovery vis-à-vis the PLE method with off-line clean-up. Feasibility of the SPLE method was further demonstrated by applying it to five real soil samples collected in Ethiopia.     

In silico prediction of drug solubility: 2. Free energy of solvation in pure melts

The solubility of drugs in water is investigated in a series of papers and in the current work. The free energy of solvation, DeltaG*(vl), of a drug molecule in its pure drug melt at 673.15 K (400 degrees C) has been obtained for 46 drug molecules using the free energy perturbation method. The simulations were performed in two steps where first the Coulomb and then the Lennard-Jones interactions were scaled down from full to no interaction. The results have been interpreted using a theory assuming that DeltaG*(vl) = DeltaG(cav) + E(LJ) + E(C)/2 where the free energy of cavity formation, DeltaG(cav), in these pure drug systems was obtained using hard body theories, and E(LJ) and E(C) are the Lennard-Jones and Coulomb interaction energies, respectively, of one molecule with the other ones. Since the main parameter in hard body theories is the volume fraction, an equation of state approach was used to estimate the molecular volume. Promising results were obtained using a theory for hard oblates, in which the oblate axial ratio was calculated from the molecular surface area and volume obtained from simulations. The Coulomb term, E(C)/2, is half of the Coulomb energy in accord with linear response, which showed good agreement with our simulation results. In comparison with our previous results on free energy of hydration, the Coulomb interactions in pure drug systems are weaker, and the van der Waals interactions play a more important role.     

In silico prediction of drug solubility. 3. Free energy of solvation in pure amorphous matter

The solubility of drugs in water is investigated in a series of papers. In this work, we address the process of bringing a drug molecule from the vapor into a pure drug amorphous phase. This step enables us to actually calculate the solubility of amorphous drugs in water. In our general approach, we, on one hand, perform rigorous free energy simulations using a combination of the free energy perturbation and thermodynamic integration methods. On the other hand, we develop an approximate theory containing parameters that are easily accessible from conventional Monte Carlo simulations, thereby reducing the computation time significantly. In the theory for solvation, we assume that DeltaG* = DeltaGcav + ELJ + EC/2, where the free energy of cavity formation, DeltaGcav, in pure drug systems is obtained using a theory for hard-oblate spheroids, and ELJ and EC are the Lennard-Jones and Coulomb interaction energies between the chosen molecule and the others in the fluid. The theoretical predictions for the free energy of solvation in pure amorphous matter are in good agreement with free energy simulation data for 46 different drug molecules. These results together with our previous studies support our theoretical approach. By using our previous data for the free energy of hydration, we compute the total free energy change of bringing a molecule from the amorphous phase into water. We obtain good agreement between the theory and simulations. It should be noted that to obtain accurate results for the total process, high precision data are needed for the individual subprocesses. Finally, for eight different substances, we compare the experimental amorphous and crystalline solubility in water with the results obtained by the proposed theory with reasonable success.     

Spectroscopic Ellipsometry Characterisation and Estimation of the Hamaker Constant of Cellulose

Calculations of Hamaker constants using Lifshitz theory require the availability of accurate dielectric data, especially in the visible-ultraviolet region. We present spectroscopic ellipsometry data on well-defined cellulose films of a limited thickness range (100–140 layers) deposited on an oxidised and hydrophobised silicon substrate. The spectral data, representing measurements from a perpendicular orientation to the fibre deposition direction, was used for estimates of the necessary spectral parameters, i.e. the oscillator strengths and characteristic frequencies in the UV-range. Our calculations show that cellulose has a relatively low Hamaker constant in air (58 zJ) and water (8.0 zJ). The implications for the surface energy estimates of cellulose and colloidal interactions between cellulose and various types of fillers and coating colours are indicated.     

Some practical aspects of free energy calculations from molecular dynamics simulation

In this work, we address two critical aspects of calculation of the free energy differences in molecular systems from molecular simulations. The first aspect involves checking whether the calculated free energy difference depends significantly on the extent of perturbation used for accomplishment of a given transformation. The second aspect of interest is to verify if the sampling errors in calculating the free energy differences between the wild-type molecule and a mutated one in its free state and in a complex are similar, or not, for a finite-length dynamic simulation. The reliability of the free energy estimates obtained from molecular simulations using thermodynamic cycles depends in part on this fact. For investigating these aspects, we use a self-transformation scheme in which a transformation of a part of a molecular system into itself is considered. We perform MD simulations of DNA fragments in which a part of a specific base is subjected to such a self-transformation. Results indicate that the estimated free energy differences do not depend significantly on the extent of perturbation used to achieve the transformation. Interestingly, the variation in the cumulative free energy difference, ΔA, with the coupling parameter, λ, depends significantly on the extent of perturbation. We examine the physical basis of the observed nature of the variation of the accumulated free energy difference, ΔA, against the λ value in the case of a self-transformation. In a thermodynamic cycle, the sampling errors due to the finite-length simulation for the molecular system are found to be similar to each other for the two perturbations (free and in a complex) justifying the use of such approach in calculating ΔΔA in molecular complexes. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 877–885, 1999     

Reconstruction of Solar Spectral Surface UV Irradiances Using Radiative Transfer Simulations

UV radiation exerts several effects concerning life on Earth, and spectral information on the prevailing UV radiation conditions is needed in order to study each of these effects. In this paper, we present a method for reconstruction of solar spectral UV irradiances at the Earth's surface. The method, which is a further development of an earlier published method for reconstruction of erythemally weighted UV, relies on radiative transfer simulations, and takes as input (1) the effective cloud optical depth as inferred from pyranometer measurements of global radiation (300–3000 nm); (2) the total ozone column; (3) the surface albedo as estimated from measurements of snow depth; (4) the total water vapor column; and (5) the altitude of the location. Reconstructed daily cumulative spectral irradiances at Jokioinen and Sodankylä in Finland are, in general, in good agreement with measurements. The mean percentage difference, for instance, is mostly within ±8%, and the root mean square of the percentage difference is around 10% or below for wavelengths over 310 nm and daily minimum solar zenith angles (SZA) less than 70°. In this study, we used pseudospherical radiative transfer simulations, which were shown to improve the performance of our method under large SZA (low Sun).     

Detecting unalloyed tin in LiSn2(PO4)3-based anodes with Mössbauer spectroscopy

The first discharge of the Li⁺ ion anode material LiSn₂(PO₄)₃ was investigated with Mössbauer spectroscopy and electrochemical techniques. Mössbauer spectroscopy provided insight into the structure of the tin atoms of the fully discharged anode materials. Spectra consist of overlapping peaks, which are assigned to noncrystalline β-Sn and Li–Sn alloy domains. An analysis of the relative intensities of the Mössbauer spectra shows the relative abundance of β-Sn increases at the expense of the Li–Sn alloy as the discharge rate increases. Cell polarization occurs at higher discharge rates, leading to inefficient electrode utilization and poor cycling performance. Sluggish Li⁺ ion diffusion through the amorphous Li₃PO₄ network that is formed early in the discharge process might be responsible for the poor electrochemical performance and the accumulation of unalloyed tin.     

Optimization of the CHARMM additive force field for DNA: Improved treatment of the BI/BII conformational equilibrium

The B-form of DNA can populate two different backbone conformations: BI and BII, defined by the difference between the torsion angles ε and ζ (BI = ε-ζ < 0 and BII = ε-ζ > 0). BI is the most populated state, but the population of the BII state, which is sequence dependent, is significant and accumulating evidence shows that BII affects the overall structure of DNA, and thus influences protein-DNA recognition. This work presents a reparametrization of the CHARMM27 additive nucleic acid force field to increase the sampling of the BII form in MD simulations of DNA. In addition, minor modifications of sugar puckering were introduced to facilitate sampling of the A form of DNA under the appropriate environmental conditions. Parameter optimization was guided by quantum mechanical data on model compounds, followed by calculations on several DNA duplexes in the condensed phase. The selected optimized parameters were then validated against a number of DNA duplexes, with the most extensive tests performed on the EcoRI dodecamer, including comparative calculations using the Amber Parm99bsc0 force field. The new CHARMM model better reproduces experimentally observed sampling of the BII conformation, including sampling as a function of sequence. In addition, the model reproduces the A form of the 1ZF1 duplex in 75 % ethanol, and yields a stable Z-DNA conformation of duplex (GTACGTAC) in its crystal environment. The resulting model, in combination with a recent reoptimization of the CHARMM27 force field for RNA, will be referred to as CHARMM36.