Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: A sequential molecular dynamics/quantum mechanical approach

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.

Electronic properties of liquid ammonia: A sequential molecular dynamics/quantum mechanics approach

The electronic properties of liquid ammonia are investigated by a sequential molecular dynamics/quantum mechanics approach. Quantum mechanics calculations for the liquid phase are based on a reparametrized hybrid exchange-correlation functional that reproduces the electronic properties of ammonia clusters [(NH(3))(n); n=1-5]. For these small clusters, electron binding energies based on Green's function or electron propagator theory, coupled cluster with single, double, and perturbative triple excitations, and density functional theory (DFT) are compared. Reparametrized DFT results for the dipole moment, electron binding energies, and electronic density of states of liquid ammonia are reported. The calculated average dipole moment of liquid ammonia (2.05 +/- 0.09 D) corresponds to an increase of 27% compared to the gas phase value and it is 0.23 D above a prediction based on a polarizable model of liquid ammonia [Deng , J. Chem. Phys. 100, 7590 (1994)]. Our estimate for the ionization potential of liquid ammonia is 9.74 +/- 0.73 eV, which is approximately 1.0 eV below the gas phase value for the isolated molecule. The theoretical vertical electron affinity of liquid ammonia is predicted as 0.16 +/- 0.22 eV, in good agreement with the experimental result for the location of the bottom of the conduction band (-V(0)=0.2 eV). Vertical ionization potentials and electron affinities correlate with the total dipole moment of ammonia aggregates. (c) 2008 American Institute of Physics.

Aspects of a noncommutative scalar/tensor duality

We study the noncommutative massless Kalb-Ramond gauge field coupled to a dynamical U(1) gauge field in the adjoint representation together with a compensating vector field. We derive the Seiberg-Witten map and obtain the corresponding mapped action to first order in theta. The (emergent) gravity structure found in other situations is not present here. The off-shell dual scalar theory is derived and it does not coincide with the Seiberg-Witten mapped scalar theory. Dispersion relations are also discussed. The p-form generalization of the Seiberg-Witten map to order theta is also derived.

Atomistic molecular dynamics study of interface formation: Al on poly(p-phenylene vinylene)

We model interface formation by metal deposition on the conjugated polymer poly-para-phenylene vinylene, studying direct aluminum and layered aluminum-calcium structures Al/PPV and Al/Ca/PPV. To do that we use classical molecular dynamics simulations, checked by ab initio density-functional theory calculations, for selected relevant configurations. We find that Al not only migrates easily into the film, with a strong charge transfer to the neighboring chains, but also promotes rearrangement of the polymer in the interfacial region to the hexagonal structure. On the other hand, our results indicate that a thin Ca layer is sufficient to protect the film and maintain a well-defined metal/polymer interface, and that also a thin Al capping layer may protect the whole from environmental degradation.

Amorphous HfO(2) and Hf(1-x)Si(x)O via a melt-and-quench scheme using ab initio molecular dynamics

We have performed ab initio molecular dynamics simulations to generate an atomic structure model of amorphous hafnium oxide (a-HfO(2)) via a melt-and-quench scheme. This structure is analyzed via bond-angle and partial pair distribution functions. These results give a Hf-O average nearest-neighbor distance of 2.2 angstrom, which should be compared to the bulk value, which ranges from 1.96 to 2.54 angstrom. We have also investigated the neutral O vacancy and a substitutional Si impurity for various sites, as well as the amorphous phase of Hf(1-x)Si(x)O(2) for x=0.25, 0375, and 0.5.

First-principles study of the adsorption of atomic and molecular hydrogen on BC(2)N nanotubes

The adsorption of atomic and molecular hydrogen on armchair and zigzag boron carbonitride nanotubes is investigated within the ab initio density functional theory. The adsorption of atomic H on the BC(2)N nanotubes presents properties which are promising for nanoelectronic applications. Depending on the adsorption site for the H, the Fermi energy moves toward the bottom of the conduction band or toward the top of the valence band, leading the system to exhibit donor or acceptor characteristics, respectively. The H(2) molecules are physisorbed on the BC(2)N surface for both chiralities. The binding energies for the H(2) molecules are slightly dependent on the adsorption site, and they are near to the range to work as a hydrogen storage medium.

Artificial molecular quantum rings under magnetic field influence

The ground states of a few electrons confined in two vertically coupled quantum rings in the presence of an external magnetic field are studied systematically within the current spin-density functional theory. Electron-electron interactions combined with inter-ring tunneling affect the electronic structure and the persistent current. For small values of the external magnetic field, we recover the zero magnetic field molecular quantum ring ground state configurations. Increasing the magnetic field many angular momentum, spin, and isospin transitions are predicted to occur in the ground state. We show that these transitions follow certain rules, which are governed by the parity of the number of electrons, the single-particle picture, Hund's rules, and many-body effects. (C) 2009 American Institute of Physics. [doi:10.1063/1.3223360]

Molecular determinants of improved cathepsin B inhibition by new cystatins obtained by DNA shuffling

Background: Cystatins are inhibitors of cysteine proteases. The majority are only weak inhibitors of human cathepsin B, which has been associated with cancer, Alzheimer's disease and arthritis. Results: Starting from the sequences of oryzacystatin-1 and canecystatin-1, a shuffling library was designed and a hybrid clone obtained, which presented higher inhibitory activity towards cathepsin B. This clone presented two unanticipated point mutations as well as an N-terminal deletion. Reversing each point mutation independently or both simultaneously abolishes the inhibitory activity towards cathepsin B. Homology modeling together with experimental studies of the reverse mutants revealed the likely molecular determinants of the improved inhibitory activity to be related to decreased protein stability. Conclusion: A combination of experimental approaches including gene shuffling, enzyme assays and reverse mutation allied to molecular modeling has shed light upon the unexpected inhibitory properties of certain cystatin mutants against Cathepsin B. We conclude that mutations disrupting the hydrophobic core of phytocystatins increase the flexibility of the N-terminus, leading to an increase in inhibitory activity. Such mutations need not affect the inhibitory site directly but may be observed distant from it and manifest their effects via an uncoupling of its three components as a result of increased protein flexibility.

3-Nitrophenol-4,4 '-bipyridyl N,N '-dioxide (2/1): a DFT study and CSD analysis of DPNO molecular complexes

The title 2:1 complex of 3-nitrophenol (MNP) and 4,4'-bipyridyl N, N'-dioxide (DPNO), 2C(6)H(5)NO(3)center dot C(10)H(8)N(2)O(2) or 2MNP center dot DPNO, crystallizes as a centrosymmetric three-component adduct with a dihedral angle of 59.40 (8)degrees between the planes of the benzene rings of MNP and DPNO (the DPNO moiety lies across a crystallographic inversion centre located at the mid-point of the C-C bond linking its aromatic rings). The complex owes its formation to O-H center dot center dot center dot O hydrogen bonds [O center dot center dot center dot O = 2.605 (3) angstrom]. Molecules are linked by intermolecular C-H center dot center dot center dot O and C-H center dot center dot center dot N interactions forming R(2)(1) (6) and R(2)(2) (10) rings, and R(6)(6) (34) and R(4)(4) (26) macro-rings, all of which are aligned along the [(1) over bar 01] direction, and R(2)(2) (10) and R(2)(1) (7) rings aligned along the [010] direction. The combination of chains of rings along the [(1) over bar 01] and [010] directions generates the three-dimensional structure. A total of 27 systems containing the DNPO molecule and forming molecular complexes of an organic nature were analysed and compared with the structural characteristics of the dioxide reported here. The N-O distance [1.325 (2) angstrom] depends not only on the interactions involving the O atom at the N-O group, but also on the structural ordering and additional three-dimensional interactions in the crystal structure. A density functional theory (DFT) optimized structure at the B3LYP/6-311G(d,p) level is compared with the molecular structure in the solid state.

Adenosine binding to low-molecular-weight purine nucleoside phosphorylase: the structural basis for recognition based on its complex with the enzyme from Schistosoma mansoni

Schistosomes are unable to synthesize purines de novo and depend exclusively on the salvage pathway for their purine requirements. It has been suggested that blockage of this pathway could lead to parasite death. The enzyme purine nucleoside phosphorylase (PNP) is one of its key components and molecules designed to inhibit the low-molecular-weight (LMW) PNPs, which include both the human and schistosome enzymes, are typically analogues of the natural substrates inosine and guanosine. Here, it is shown that adenosine both binds to Schistosoma mansoni PNP and behaves as a weak micromolar inhibitor of inosine phosphorolysis. Furthermore, the first crystal structures of complexes of an LMW PNP with adenosine and adenine are reported, together with those with inosine and hypoxanthine. These are used to propose a structural explanation for the selective binding of adenosine to some LMW PNPs but not to others. The results indicate that transition-state analogues based on adenosine or other 6-amino nucleosides should not be discounted as potential starting points for alternative inhibitors.

Evolutionary dynamics on rugged fitness landscapes: Exact dynamics and information theoretical aspects

The parallel mutation-selection evolutionary dynamics, in which mutation and replication are independent events, is solved exactly in the case that the Malthusian fitnesses associated to the genomes are described by the random energy model (REM) and by a ferromagnetic version of the REM. The solution method uses the mapping of the evolutionary dynamics into a quantum Ising chain in a transverse field and the Suzuki-Trotter formalism to calculate the transition probabilities between configurations at different times. We find that in the case of the REM landscape the dynamics can exhibit three distinct regimes: pure diffusion or stasis for short times, depending on the fitness of the initial configuration, and a spin-glass regime for large times. The dynamic transition between these dynamical regimes is marked by discontinuities in the mean-fitness as well as in the overlap with the initial reference sequence. The relaxation to equilibrium is described by an inverse time decay. In the ferromagnetic REM, we find in addition to these three regimes, a ferromagnetic regime where the overlap and the mean-fitness are frozen. In this case, the system relaxes to equilibrium in a finite time. The relevance of our results to information processing aspects of evolution is discussed.

Density functional theory and ab initio molecular dynamics study of NO adsorption on Pd(111) and Pt(111) surfaces

The origin of the unique geometry for nitric oxide (NO) adsorption on Pd(111) and Pt(111) surfaces as well as the effect of temperature were studied by density functional theory calculations and ab initio molecular dynamics at finite temperature. We found that at low coverage, the adsorption geometry is determined by electronic interactions, depending sensitively on the adsorption sites and coverages, and the effect of temperature on geometries is significant. At coverage of 0.25 monolayer (ML), adsorbed NO at hollow sites prefer an upright configuration, while NO adsorbed at top sites prefer a tilting configuration. With increase in the coverage up to 0.50 ML, the enhanced steric repulsion lead to the tilting of hollow NO. We found that the tilting was enhanced by the thermal effects. At coverage of 0.75 ML with p(2 x 2)-3NO(fcc+hcp+top) structure, we found that there was no preferential orientation for tilted top NO. The interplay of the orbital hybridization, thermal effects, steric repulsion, and their effects on the adsorption geometries were highlighted at the end.

Determination of the molecular weight of proteins in solution from a single small-angle X-ray scattering measurement on a relative scale

This paper describes a new and simple method to determine the molecular weight of proteins in dilute solution, with an error smaller than similar to 10%, by using the experimental data of a single small-angle X-ray scattering (SAXS) curve measured on a relative scale. This procedure does not require the measurement of SAXS intensity on an absolute scale and does not involve a comparison with another SAXS curve determined from a known standard protein. The proposed procedure can be applied to monodisperse systems of proteins in dilute solution, either in monomeric or multimeric state, and it has been successfully tested on SAXS data experimentally determined for proteins with known molecular weights. It is shown here that the molecular weights determined by this procedure deviate from the known values by less than 10% in each case and the average error for the test set of 21 proteins was 5.3%. Importantly, this method allows for an unambiguous determination of the multimeric state of proteins with known molecular weights.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Non-Cardiac Surgery in Developing Countries: Epidemiological Aspects and Economical Opportunities - The Case of Brazil

Background: Worldwide distribution of surgical interventions is unequal. Developed countries account for the majority of surgeries and information about non-cardiac operations in developing countries is scarce. The purpose of our study was to describe the epidemiological data of non-cardiac surgeries performed in Brazil in the last years. Methods and Findings: This is a retrospective cohort study that investigated the time window from 1995 to 2007. We collected information from DATASUS, a national public health system database. The following variables were studied: number of surgeries, in-hospital expenses, blood transfusion related costs, length of stay and case fatality rates. The results were presented as sum, average and percentage. The trend analysis was performed by linear regression model. There were 32,659,513 non-cardiac surgeries performed in Brazil in thirteen years. An increment of 20.42% was observed in the number of surgeries in this period and nowadays nearly 3 million operations are performed annually. The cost of these procedures has increased tremendously in the last years. The increment of surgical cost was almost 200%. The total expenses related to surgical hospitalizations were more than $10 billion in all these years. The yearly cost of surgical procedures to public health system was more than $1.27 billion for all surgical hospitalizations, and in average, U$445.24 per surgical procedure. The total cost of blood transfusion was near $98 million in all years and annually approximately $10 million were spent in perioperative transfusion. The surgical mortality had an increment of 31.11% in the period. Actually, in 2007, the surgical mortality in Brazil was 1.77%. All the variables had a significant increment along the studied period: r square (r(2)) = 0.447 for the number of surgeries (P = 0.012), r(2) = 0.439 for in-hospital expenses (P = 0.014) and r(2) = 0.907 for surgical mortality (P = 0.0055). Conclusion: The volume of surgical procedures has increased substantially in Brazil through the past years. The expenditure related to these procedures and its mortality has also increased as the number of operations. Better planning of public health resource and strategies of investment are needed to supply the crescent demand of surgery in Brazil.