Molecular dynamics simulation of organic-inorganic nanocomposites: Layering behavior and interlayer structure of organoclays

Isothermal-isobaric (NPT) molecular dynamics simulation has been performed to investigate the layering behavior and structure of nanoconfined quaternary alkylammoniums in organoclays. This work is focused on systems consisting of two clay layers and a number of alkylammoniums, and involves the use of modified Dreiding force field. The simulated basal spacings of organoclays agree satisfactorily with the experimental results in the literature. The atomic density profiles in the direction normal to the clay surface indicate that the alkyl chains within the interlayer space of montmorillonite exhibit an obvious layering behavior. The headgroups of long alkyl chains are distributed within two layers close to the clay surface, whereas the distributions of methyl and methylene groups are strongly dependent on the alkyl chain length and clay layer charge. Monolayer, bilayer, and pseudo-trilayer structures are found in organoclays modified with single long alkyl chains, which are identical to the structural models based on the measured basal spacings. A pseudo-quadrilayer structure, for the first time to our knowledge, is also identified in organoclays with double long alkyl chains. In the mixture structure of paraffin-type and multilayer, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer in pseudo-trilayer as well as next nearest layer in pseudo-quadrilayer.

Optimal design for model discrimination and parameter estimation for itraconazole population pharmacokinetics in cystic fibrosis patients

Optimal sampling times are found for a study in which one of the primary purposes is to develop a model of the pharmacokinetics of itraconazole in patients with cystic fibrosis for both capsule and solution doses. The optimal design is expected to produce reliable estimates of population parameters for two different structural PK models. Data collected at these sampling times are also expected to provide the researchers with sufficient information to reasonably discriminate between the two competing structural models.

Analysis of population pharmacokinetic data using NONMEM and WinBUGS

The aim of this report is to describe the use of WinBUGS for two datasets that arise from typical population pharmacokinetic studies. The first dataset relates to gentamicin concentration-time data that arose as part of routine clinical care of 55 neonates. The second dataset incorporated data from 96 patients receiving enoxaparin. Both datasets were originally analyzed by using NONMEM. In the first instance, although NONMEM provided reasonable estimates of the fixed effects parameters it was unable to provide satisfactory estimates of the between-subject variance. In the second instance, the use of NONMEM resulted in the development of a successful model, albeit with limited available information on the between-subject variability of the pharmacokinetic parameters. WinBUGS was used to develop a model for both of these datasets. Model comparison for the enoxaparin dataset was performed by using the posterior distribution of the log-likelihood and a posterior predictive check. The use of WinBUGS supported the same structural models tried in NONMEM. For the gentamicin dataset a one-compartment model with intravenous infusion was developed, and the population parameters including the full between-subject variance-covariance matrix were available. Analysis of the enoxaparin dataset supported a two compartment model as superior to the one-compartment model, based on the posterior predictive check. Again, the full between-subject variance-covariance matrix parameters were available. Fully Bayesian approaches using MCMC methods, via WinBUGS, can offer added value for analysis of population pharmacokinetic data.

A quantum yield map for synthetic eumelanin

The quantum yield of synthetic eumelanin is known to be extremely low and it has recently been reported to be dependent on excitation wavelength. In this paper, we present quantum yield as a function of excitation wavelength between 250 and 500 nm, showing it to be a factor of 4 higher at 250 nm than at 500 nm. In addition, we present a definitive map of the steady-state fluorescence as a function of excitation and emission wavelengths, and significantly, a three-dimensional map of the specific quantum yield: the fraction of photons absorbed at each wavelength that are subsequently radiated at each emission wavelength. This map contains clear features, which we attribute to certain structural models, and shows that radiative emission and specific quantum yield are negligible at emission wavelengths outside the range of 585 and 385 nm (2.2 and 3.2 eV), regardless of excitation wavelength. This information is important in the context of understanding melanin biofunctionality, and the quantum molecular biophysics therein. (c) 2005 American Institute of Physics.

Quantitative fluorescence excitation spectra of synthetic eumelanin

Previously reported excitation spectra for eumelanin are sparse and inconsistent. Moreover, these studies have failed to account for probe beam attenuation and emission reabsorption within the samples, making them qualitative at best. We report for the first time quantitative excitation spectra for synthetic eumelanin, acquired for a range of solution concentrations and emission wavelengths. Our data indicate that probe beam attenuation and emission reabsorption significantly affect the spectra even in low-concentration eumelanin solutions and that previously published data do not reflect the true excitation profile. We apply a correction procedure (previously applied to emission spectra) to account for these effects. Application of this procedure reconstructs the expected relationship of signal intensity with concentration, and the normalized spectra show a similarity in form to the absorption profiles. These spectra reveal valuable information regarding the photophysics and photochemistry of eumelanin. Most notably, an excitation peak at 365 urn (3.40 eV), whose position is independent of emission wavelength, is possibly attributable to a 5,6-dihydroxyindole-2-carboxylic acid (DHICA) component singly linked to a polymeric structure.

Mimicking the action of GroEL in molecular dynamics simulations: Application to the refinement of protein structures

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results Sur port the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

The physical and chemical properties of eumelanin

In this article, we review the current state of knowledge concerning the physical and chemical properties of the eumelanin pigment. We examine properties related to its photoprotective functionality, and draw the crucial link between fundamental molecular structure and observable macroscopic behaviour. Where necessary, we also briefly review certain aspects of the pheomelanin literature to draw relevant comparison. A full understanding of melanin function, and indeed its role in retarding or promoting the disease state, can only be obtained through a full mapping of key structure-property relationships in the main pigment types. We are engaged in such an endeavor for the case of eumelanin.

Radiative relaxation in synthetic pheomelanin

We report a detailed photoluminescence study of cysteinyldopa-melanin ( CDM), the synthetic analogue of pheomelanin. Emission spectra are shown to be a far more sensitive probe of CDM's spectroscopic behavior than are absorption spectra. Although CDM and dopa-melanin ( DM, the synthetic analogue of eumelanin) have very similar absorption spectra, we find that they have very different excitation and emission characteristics; CDM has two distinct photoluminescence peaks that do not shift with excitation wavelength. Additionally, our data suggest that the radiative quantum yield of CDM is excitation energy dependent, an unusual property among biomolecules that is indicative of a chemically disordered system. Finally, we find that the radiative quantum yield for CDM is similar to 0.2%, twice that of DM, although still extremely low. This means that 99.8% of the energy absorbed by CDM is dissipated via nonradiative pathways, consistent with its role as a pigmentary photoprotectant.

Quantitative scattering of melanin solutions

The optical scattering coefficient of a dilute, well-solubilized eumelanin solution has been accurately measured as a function of incident wavelength, and found to contribute < 6% of the total optical attenuation between 210 and 325 nm. At longer wavelengths (325-800 nm), the scattering was less than the minimum sensitivity of our instrument. This indicates that ultraviolet and visible optical density spectra can be interpreted as true absorption with a high degree of confidence. The scattering coefficient versus wavelength was found to be consistent with Rayleigh theory for a particle radius of 38 6 1 nm. Our results shed important light on the role of melanins as photoprotectants.

Using the canonical modelling approach to simplify the simulation of function in functional-structural plant models

Functional-structural plant models that include detailed mechanistic representation of underlying physiological processes can be expensive to construct and the resulting models can also be extremely complicated. On the other hand, purely empirical models are not able to simulate plant adaptability and response to different conditions. In this paper, we present an intermediate approach to modelling plant function that can simulate plant response without requiring detailed knowledge of underlying physiology. Plant function is modelled using a 'canonical' modelling approach, which uses compartment models with flux functions of a standard mathematical form, while plant structure is modelled using L-systems. Two modelling examples are used to demonstrate that canonical modelling can be used in conjunction with L-systems to create functional-structural plant models where function is represented either in an accurate and descriptive way, or in a more mechanistic and explanatory way. We conclude that canonical modelling provides a useful, flexible and relatively simple approach to modelling plant function at an intermediate level of abstraction.

Validation and refinement of survival models for liver retransplantation

Orthotopic liver retransplantation (re-OLT) is highly controversial. The objectives of this study were to determine the validity of a recently developed United Network for Organ Sharing (UNOS) multivariate model using an independent cohort of patients undergoing re-OLT outside the United States, to determine whether incorporation of other variables that were incomplete in the UNOS registry would provide additional prognostic information, to develop new models combining data sets from both cohorts, and to evaluate the validity of the model for end-stage liver disease (MELD) in patients undergoing re-OLT. Two hundred eighty-one adult patients undergoing re-OLT (between 1986 and 1999) at 6 foreign transplant centers comprised the validation cohort. We found good agreement between actual survival and predicted survival in the validation cohort; 1-year patient survival rates in the low-, intermediate-, and high-risk groups (as assigned by the original UNOS model) were 72%, 68%, and 36%, respectively (P < .0001). In the patients for whom the international normalized ratio (INR) of prothrombin time was available, MELD correlated with outcome following re-OLT; the median MELD scores for patients surviving at least 90 days compared with those dying within 90 days were 20.75 versus 25.9, respectively (P = .004). Utilizing both patient cohorts (n = 979), a new model, based on recipient age, total serum bilirubin, creatinine, and interval to re-OLT, was constructed (whole model χ(2) = 105, P < .0001). Using the c-statistic with 30-day, 90-day, 1-year, and 3-year mortality as the end points, the area under the receiver operating characteristic (ROC) curves for 4 different models were compared. In conclusion, prospective validation and use of these models as adjuncts to clinical decision making in the management of patients being considered for re-OLT are warranted.

Hurt feelings in couple relationships: Towards integrative models of the negative effects of hurtful events

Previous research suggests that hurt feelings can have powerful effects on individual and relational outcomes. This study examined a typology of hurtful events in couple relationships, together with integrative models predicting ongoing effects on victims and relationships. Participants were 224 students from introductory and third-year psychology classes, who completed open-ended and structured measures concerning an event in which a partner had hurt their feelings. By tailoring Leary et al.'s (1998) typology to the context of romantic relationships, five categories of hurtful events were proposed: active disassociation, passive disassociation, criticism, infidelity, and deception. Analyses assessing similarities and differences among the categories confirmed the utility of the typology. Structural equation modeling showed that longer-term effects on the victim were predicted by relationship anxiety and by the victim's immediate reactions to the event (negative emotions and self-perceptions; feelings of rejection and powerlessness). In contrast, ongoing effects on the relationship were predicted by avoidance, the victim's attributions and perceptions of offender remorse, and the victim's own behavior. The results highlight the utility of an integrated approach to hurt, incorporating emotional, cognitive, and behavioral responses, and dimensions of attachment security.