Fluorescence of nitrobenzoxadiazole (NBD)-labeled lipids in model membranes is connected not to lipid mobility but to probe location

Nitrobenzoxadiazole (NBD)-labeled lipids are popular fluorescent membrane probes. However, the understanding of important aspects of the photophysics of NBD remains incomplete, including the observed shift in the emission spectrum of NBD-lipids to longer wavelengths following excitation at the red edge of the absorption spectrum (red-edge excitation shift or REES). REES of NBD-lipids in membrane environments has been previously interpreted as reflecting restricted mobility of solvent surrounding the fluorophore. However, this requires a large change in the dipole moment (Dm) of NBD upon excitation. Previous calculations of the value of Dm of NBD in the literature have been carried out using outdated semi-empirical methods, leading to conflicting values. Using up-to-date density functional theory methods, we recalculated the value of Dm and verified that it is rather small (B2 D). Fluorescence measurements confirmed that the value of REES is B16 nm for 1,2-dioleoyl-sn-glycero-3- phospho-L-serine-N-(NBD) (NBD-PS) in dioleoylphosphatidylcholine vesicles. However, the observed shift is independent of both the temperature and the presence of cholesterol and is therefore insensitive to the mobility and hydration of the membrane. Moreover, red-edge excitation leads to an increased contribution of the decay component with a shorter lifetime, whereas time-resolved emission spectra of NBD-PS displayed an atypical blue shift following excitation. This excludes restrictions to solvent relaxation as the cause of the measured REES and TRES of NBD, pointing instead to the heterogeneous transverse location of probes as the origin of these effects. The latter hypothesis was confirmed by molecular dynamics simulations, from which the calculated heterogeneity of the hydration and location of NBD correlated with the measured fluorescence lifetimes/REES. Globally, our combination of theoretical and experiment-based techniques has led to a considerably improved understanding of the photophysics of NBD and a reinterpretation of its REES in particular.

A multiscale multisurface constitutive model for the thermo-plastic behavior of polyethylene

We present a multiscale model bridging length and time scales from molecular to continuum levels with the objective of predicting the yield behavior of amorphous glassy polyethylene (PE). Constitutive pa- rameters are obtained from molecular dynamics (MD) simulations, decreasing the requirement for ad- hoc experiments. Consequently, we achieve: (1) the identification of multisurface yield functions; (2) the high strain rate involved in MD simulations is upscaled to continuum via quasi-static simulations. Validation demonstrates that the entire multisurface yield functions can be scaled to quasi-static rates where the yield stresses are possibly predicted by a proposed scaling law; (3) a hierarchical multiscale model is constructed to predict temperature and strain rate dependent yield strength of the PE.

Mechanical Properties And Fracture Dynamics Of Silicene Membranes.

As graphene has become one of the most important materials, there is renewed interest in other similar structures. One example is silicene, the silicon analogue of graphene. It shares some of the remarkable graphene properties, such as the Dirac cone, but presents some distinct ones, such as a pronounced structural buckling. We have investigated, through density functional based tight-binding (DFTB), as well as reactive molecular dynamics (using ReaxFF), the mechanical properties of suspended single-layer silicene. We calculated the elastic constants, analyzed the fracture patterns and edge reconstructions. We also addressed the stress distributions, unbuckling mechanisms and the fracture dependence on the temperature. We analysed the differences due to distinct edge morphologies, namely zigzag and armchair.

Structural basis of GC-1 selectivity for thyroid hormone receptor isoforms

Background: Thyroid receptors, TRa and TR beta, are involved in important physiological functions such as metabolism, cholesterol level and heart activities. Whereas metabolism increase and cholesterol level lowering could be achieved by TR beta isoform activation, TRa activation affects heart rates. Therefore, beta-selective thyromimetics have been developed as promising drug-candidates for treatment of obesity and elevated cholesterol level. GC-1 [ 3,5-dimethyl-4-(4'-hydroxy- 3'-isopropylbenzyl)-phenoxy acetic acid] has ability to lower LDL cholesterol with 600-to 1400-fold more potency and approximately two-to threefold more efficacy than atorvastatin (Lipitor(C)) in studies in rats, mice and monkeys. Results: To investigate GC-1 specificity, we solved crystal structures and performed molecular dynamics simulations of both isoforms complexed with GC-1. Crystal structures reveal that, in TRa Arg228 is observed in multiple conformations, an effect triggered by the differences in the interactions between GC-1 and Ser277 or the corresponding asparagine (Asn331) of TR beta. The corresponding Arg282 of TR beta is observed in only one single stable conformation, interacting effectively with the ligand. Molecular dynamics support this model: our simulations show that the multiple conformations can be observed for the Arg228 in TR alpha, in which the ligand interacts either strongly with the ligand or with the Ser277 residue. In contrast, a single stable Arg282 conformation is observed for TR beta, in which it strongly interacts with both GC-1 and the Asn331. Conclusion: Our analysis suggests that the key factors for GC-1 selectivity are the presence of an oxyacetic acid ester oxygen and the absence of the amino group relative to T(3). These results shed light into the beta-selectivity of GC-1 and may assist the development of new compounds with potential as drug candidates to the treatment of hypercholesterolemia and obesity.

Hydrogen bond interactions between acetone and supercritical water

Hydrogen bond interactions between acetone and supercritical water are investigated using a combined and sequential Monte Carlo/quantum mechanics (S-MC/QM) approach. Simulation results show a dominant presence of con. gurations with one hydrogen bond for different supercritical states, indicating that this specific interaction plays an important role on the solvation properties of acetone in supercritical water. Using QM MP2/aug-cc-pVDZ the calculated average interaction energy reveals that the hydrogen-bonded acetone-water complex is energetically more stable under supercritical conditions than ambient conditions and its stability is little affected by variations of temperature and/or pressure. All average results reported here are statistically converged.

Adiabatic intramolecular movements for water systems

An effective treatment of the intramolecular degrees of freedom is presented for water, where these modes are decoupled from the intermolecular ones, ""adiabatically"" allowing these coordinates to be positioned at their local minimum of the potential energy surface. We perform ab initio Monte Carlo simulations with the configurational energies obtained via density functional theory. We study a water dimer as a prototype system, and even in this simple case the intramolecular relaxations are very important to properly describe properties such as the dipole moment. We show that rigid simulations do not correctly sample the phase space, resulting in an average dipole moment smaller than the one obtained with the adiabatic model, which is closer to the experimental result. (c) 2008 American Institute of Physics.

Intermediate motions and dipolar couplings as studied by Lee-Goldburg cross-polarization NMR: Hartmann-Hahn matching profiles

In this article, we evaluate the use of simple Lee-Goldburg cross-polarization (LG-CP) NMR experiments for obtaining quantitative information of molecular motion in the intermediate regime. In particular, we introduce the measurement of Hartmann-Hahn matching profiles for the assessment of heteronuclear dipolar couplings as well as dynamics as a reliable and robust alternative to the more common analysis of build-up curves. We have carried out dynamic spin dynamics simulations in order to test the method's sensitivity to intermediate motion and address its limitations concerning possible experimental imperfections. We further demonstrate the successful use of simple theoretical concepts, most prominently Anderson-Weiss (AW) theory, to analyze the data. We further propose an alternative way to estimate activation energies of molecular motions, based upon the acquisition of only two LG-CP spectra per temperature at different temperatures. As experimental tests, molecular jumps in imidazole methyl sulfonate, trimethylsulfoxonium iodide, and bisphenol A polycarbonate were investigated with the new method.

HFE gene mutations in patients with primary iron overload: Is there a significant improvement in molecular diagnosis yield with HFE sequencing?

Rare HFE variants have been shown to be associated with hereditary hemochromatosis (HH), an iron overload disease. The low frequency of the HFE p.C282Y mutation in HH-affected Brazilian patients may suggest that other HFE-related mutations may also be implicated in the pathogenesis of HH in this population. The main aim was to screen for new HFE mutations in Brazilian individuals with primary iron overload and to investigate their relationship with HH. Fifty Brazilian patients with primary iron overload (transferrin saturation >50% in females and 60% in males) were selected. Subsequent bidirectional sequencing for each HFE exon was performed. The effect of HFE mutations on protein structure were analyzed by molecular dynamics simulation and free binding energy calculations. p.C282Y in homozygosis or in heterozygosis with p.H63D were the most frequent genotypic combinations associated with HH in our sample population (present in 17 individuals, 34%). Thirty-six (72.0%) out of the 50 individuals presented at least one HFE mutation. The most frequent genotype associated with HH was the homozygous p.C282Y mutation (n = 11, 22.0%). One novel mutation (p.V256I) was indentified in heterozygosis with the p.H63D mutation. In silico modeling analysis of protein behavior indicated that the p.V256I mutation does not reduce the binding affinity between HFE and beta 2-microglobulin ((beta 2M) in the same way the p.C282Y mutation does compared with the native HFE protein. In conclusion, screening of HFE through direct sequencing, as compared to p.C282Y/p.H63D genotyping, was not able to increase the molecular diagnosis yield of HH. The novel p.V256I mutation could not be implicated in the molecular basis of the HH phenotype, although its role cannot be completely excluded in HH-phenotype development. Our molecular modeling analysis can help in the analysis of novel, previously undescribed, HFE mutations. (C) 2010 Elsevier Inc. All rights reserved.

QSAR Modeling of a Set of Pyrazinoate Esters as Antituberculosis Prodrugs

Tuberculosis is an infection caused mainly by Mycobacterium tuberculosis. A first-line antimycobacterial drug is pyrazinamide (PZA), which acts partially as a prodrug activated by a pyrazinamidase releasing the active agent, pyrazinoic acid (POA). As pyrazinoic acid presents some difficulty to cross the mycobacterial cell wall, and also the pyrazinamide-resistant strains do not express the pyrazinamidase, a set of pyrazinoic acid esters have been evaluated as antimycobacterial agents. In this work, a QSAR approach was applied to a set of forty-three pyrazinoates against M. tuberculosis ATCC 27294, using genetic algorithm function and partial least squares regression (WOLF 5.5 program). The independent variables selected were the Balaban index (I), calculated n-octanol/water partition coefficient (ClogP), van-der-Waals surface area, dipole moment, and stretching-energy contribution. The final QSAR model (N = 32, r(2) = 0.68, q(2) = 0.59, LOF = 0.25, and LSE = 0.19) was fully validated employing leave-N-out cross-validation and y-scrambling techniques. The test set (N = 11) presented an external prediction power of 73%. In conclusion, the QSAR model generated can be used as a valuable tool to optimize the activity of future pyrazinoic acid esters in the designing of new antituberculosis agents.

Molecular cloning and biochemical characterization of a myotoxin inhibitor from Bothrops alternatus snake plasma

Phospholipases A(2) (PLA(2)s) are important components of Bothrops snake venoms, that can induce several effects on envenomations such as myotoxicity, inhibition or induction of platelet aggregation and edema. It is known that venomous and non-venomous snakes present PLA(2) inhibitory proteins (PLIs) in their blood plasma. An inhibitory protein that neutralizes the enzymatic and toxic activities of several PLA2s from Bothrops venoms was isolated from Bothrops alternatus snake plasma by affinity chromatography using the immobilized myotoxin BthTX-I on CNBr-activated Sepharose. Biochemical characterization of this inhibitory protein, denominated alpha BaltMIP, showed it to be a glycoprotein with Mr of similar to 24,000 for the monomeric subunit. CD spectra of the PLA(2)/inhibitor complexes are considerably different from those corresponding to the individual proteins and data deconvolution suggests that the complexes had a relative gain of helical structure elements in comparison to the individual protomers, which may indicate a more compact structure upon complexation. Theoretical and experimental structural studies performed in order to obtain insights into the structural features of aBaltMIP indicated that this molecule may potentially trimerize in solution, thus strengthening the hypothesis previously raised by other authors about snake PLIs oligomerization. (C) 2010 Elsevier Masson SAS. All rights reserved.

Computational chemistry on Fujitsu vector-parallel processors: Development and performance of applications software

In this and a preceding paper, we provide an introduction to the Fujitsu VPP range of vector-parallel supercomputers and to some of the computational chemistry software available for the VPP. Here, we consider the implementation and performance of seven popular chemistry application packages. The codes discussed range from classical molecular dynamics to semiempirical and ab initio quantum chemistry. All have evolved from sequential codes, and have typically been parallelised using a replicated data approach. As such they are well suited to the large-memory/fast-processor architecture of the VPP. For one code, CASTEP, a distributed-memory data-driven parallelisation scheme is presented. (C) 2000 Published by Elsevier Science B.V. All rights reserved.

A model for adsorption of condensable vapors on nonporous materials

A new isotherm is proposed here for adsorption of condensable vapors and gases on nonporous materials having type II isotherms according to the Brunauer-Deming-Deming-Teller (BDDH) classification. The isotherm combines the recent molecular-continuum model in the multilayer region, with other widely used models for sub-monolayer coverage, some of which satisfy the requirement of a Henry's law asymptote. The model is successfully tested using isotherm data for nitrogen adsorption on nonporous silica, carbon and alumina, as well as benzene and hexane adsorption on nonporous carbon. Based on the data fits, out of several different alternative choices of model for the monolayer region, the Freundlich and the Unilan models are found to be the most successful when combined with the multilayer model to predict the whole isotherm. The hybrid model is consequently applicable over a wide pressure range. (C) 2000 Elsevier Science B.V. All rights reserved.

A landscape-scale test of the predictive ability of a spatially explicit model for population viability analysis

1. Although population viability analysis (PVA) is widely employed, forecasts from PVA models are rarely tested. This study in a fragmented forest in southern Australia contrasted field data on patch occupancy and abundance for the arboreal marsupial greater glider Petauroides volans with predictions from a generic spatially explicit PVA model. This work represents one of the first landscape-scale tests of its type. 2. Initially we contrasted field data from a set of eucalypt forest patches totalling 437 ha with a naive null model in which forecasts of patch occupancy were made, assuming no fragmentation effects and based simply on remnant area and measured densities derived from nearby unfragmented forest. The naive null model predicted an average total of approximately 170 greater gliders, considerably greater than the true count (n = 81). 3. Congruence was examined between field data and predictions from PVA under several metapopulation modelling scenarios. The metapopulation models performed better than the naive null model. Logistic regression showed highly significant positive relationships between predicted and actual patch occupancy for the four scenarios (P = 0.001-0.006). When the model-derived probability of patch occupancy was high (0.50-0.75, 0.75-1.00), there was greater congruence between actual patch occupancy and the predicted probability of occupancy. 4. For many patches, probability distribution functions indicated that model predictions for animal abundance in a given patch were not outside those expected by chance. However, for some patches the model either substantially over-predicted or under-predicted actual abundance. Some important processes, such as inter-patch dispersal, that influence the distribution and abundance of the greater glider may not have been adequately modelled. 5. Additional landscape-scale tests of PVA models, on a wider range of species, are required to assess further predictions made using these tools. This will help determine those taxa for which predictions are and are not accurate and give insights for improving models for applied conservation management.

Macromolecular assembly of polycystin-2 intracytosolic C-terminal domain

Mutations in PKD2 are responsible for approximately 15% of the autosomal dominant polycystic kidney disease cases. This gene encodes polycystin-2, a calcium-permeable cation channel whose C-terminal intracytosolic tail (PC2t) plays an important role in its interaction with a number of different proteins. In the present study, we have comprehensively evaluated the macromolecular assembly of PC2t homooligomer using a series of biophysical and biochemical analyses. Our studies, based on a new delimitation of PC2t, have revealed that it is capable of assembling as a homotetramer independently of any other portion of the molecule. Our data support this tetrameric arrangement in the presence and absence of calcium. Molecular dynamics simulations performed with a modified all-atoms structure-based model supported the PC2t tetrameric assembly, as well as how different populations are disposed in solution. The simulations demonstrated, indeed, that the best-scored structures are the ones compatible with a fourfold oligomeric state. These findings clarify the structural properties of PC2t domain and strongly support a homotetramer assembly of PC2.

Ring flexibility within tricyclic antidepressant drugs

The internal flexibility of the central seven-membered ring of a series of tricyclic antidepressant drugs (TCAs), imipramine {l}, amitriptyline {2}, doxepin {3}, and dothiepin {4}, has been investigated by H-1 and C-13 nuclear magnetic (NMR) techniques. Two dynamic processes were examined: ring inversion and bridge flexing. H-1 NMR lineshape analysis was used to obtain ring inversion barriers for 2-4. These studies yielded energy barriers of 14.3, 16.7, and 15.7 +/- 0.6 kcal/mol for the hydrochloride salts of doxepin, dothiepin, and amitriptyline, respectively. The barriers for the corresponding free bases were lower by 0.6 kcal/mol on average. (CT1)-C-13 relaxation measurements were used to determine the degree of bridge flexing associated with the central seven-membered ring for all four compounds. By fitting the T-1 data to a two-state jump model, lifetimes and amplitudes of rapid bridge flexing motions were determined. The results show that imipramine has the fastest rate of bridge flexing, followed by amitriptyline, doxepin, and dothiepin. The pharmacological profiles of the TCAs are complex and they interact with many receptor sites, resulting in numerous side effects and a general lack of understanding of their precise mode of action in different anxiety-related disorders. They all have similar three-dimensional structures, which makes it difficult to rationalize their differing relative potency in different assays/clinical settings. However, the clear finding here that there are significantly different degrees of internal mobility suggests that molecular dynamics should be an additional factor considered when trying to understand the mode of action of this clinically important family of molecules. (C) 2001 Wiley-Liss, Inc. and the American Pharmaceutical Association J Pharm Sci 90:713-721, 2001.