Computational chemistry for elucidating protein function: energetics and dynamics of myoglobin-ligand systems

State-of-the-art computational methodologies are used to investigate the energetics and dynamics of photodissociated CO and NO in myoglobin (Mb···CO and Mb···NO). This includes the combination of molecular dynamics, ab initio MD, free energy sampling, and effective dynamics methods to compare the results with studies using X-ray crystallography and ultrafast spectroscopy metho ds. It is shown that modern simulation techniques along with careful description of the intermolecular interactions can give quantitative agreement with experiments on complex molecular systems. Based on this agreement predictions for as yet uncharacterized species can be made.

Three-centre hydrogen bonds in triphenylphosphine oxide-hydroquinone (1/1)

The title cocrystal, C18H15OP center dot C6H6O2, belongs to a series of molecular systems based on triphenylphosphine P-oxide. The O atom of the oxide group acts as an acceptor for hydrogen bonds from OH groups of two hydroquinone molecules which lie on inversion centres [O center dot center dot center dot O = 2.7451 (17) and 2.681 (2) A S]. The crystal structure is stabilized by weak C-H center dot center dot center dot O hydrogen bonds, forming a C-2(1)(8) chain which runs parallel to the [100] direction.

Systemic multicompartmental effects of the gut microbiome on mouse metabolic phenotypes

To characterize the impact of gut microbiota on host metabolism, we investigated the multicompartmental metabolic profiles of a conventional mouse strain (C3H/HeJ) (n=5) and its germ-free (GF) equivalent (n=5). We confirm that the microbiome strongly impacts on the metabolism of bile acids through the enterohepatic cycle and gut metabolism (higher levels of phosphocholine and glycine in GF liver and marked higher levels of bile acids in three gut compartments). Furthermore we demonstrate that (1) well-defined metabolic differences exist in all examined compartments between the metabotypes of GF and conventional mice: bacterial co-metabolic products such as hippurate (urine) and 5-aminovalerate (colon epithelium) were found at reduced concentrations, whereas raffinose was only detected in GF colonic profiles. (2) The microbiome also influences kidney homeostasis with elevated levels of key cell volume regulators (betaine, choline, myo-inositol and so on) observed in GF kidneys. (3) Gut microbiota modulate metabotype expression at both local (gut) and global (biofluids, kidney, liver) system levels and hence influence the responses to a variety of dietary modulation and drug exposures relevant to personalized health-care investigations.

Conformational analysis of oriented non-crystalline polymers using wide angle X-ray scattering

A procedure is presented for obtaining conformational parameters from oriented but non-crystalline polymers. This is achieved by comparison of the experimental wide angle X-ray scattering with that calculated from models but in such a way that foreknowledge of the orientation distribution function is not required. X-ray scattering intensity values for glassy isotactic poly(methylmethacrylate) are analysed by these techniques. The method could be usefully applied to other oriented molecular systems such as liquid crystalline materials.

Electrospun supramolecular polymer fibres

The electrospinning of urethane based low molecular weight polymers differing only in the nature of the hydrogen bonding end-groups has been investigated. For the end-groups with the lowest binding constants at maximum solubility only droplets, are produced at the electrode; in contrast, increasing the binding constant of the end-group results in electrospun fibres being produced. The properties of the fibres produced are subject to changes in solvent, concentration and temperature. Typical diameters for these fibres were found to be some 10 s of μm, rather than the sub-micron dimensions often produced in electrospinning systems. Such diameters are related to the high initial concentrations required; this also may influence the rate of solvent removal and preferential surface solidification which feature in these examples. A simple theoretical model is used to relate the association constant to the molecular weight required for fibre formation; significantly lower levels of association are required for higher molecular weight macromonomers compared to smaller molecular systems.

Substitutional and orientational disorder in organic crystals: a symmetry-adapted ensemble model

Modelling of disorder in organic crystals is highly desirable since it would allow thermodynamic stabilities and other disorder-sensitive properties to be estimated for such systems. Two disordered organic molecular systems are modeled using a symmetry-adapted ensemble approach, in which the disordered system is treated as an ensemble of the configurations of a supercell with respect to substitution of one disorder component for another. Computation time is kept manageable by performing calculations only on the symmetrically inequivalent configurations. Calculations are presented on a substitutionally disordered system, the dichloro/dibromobenzene solid solution, and on an orientationally disordered system, eniluracil, and the resultant free energies, disorder patterns, and system properties are discussed. The results are found to be in agreement with experiment following manual removal of physically implausible configurations from ensemble averages, highlighting the dangers of a completely automated approach to organic crystal thermodynamics which ignores the barriers to equilibration once the crystal has been formed.

Elucidating the structure of chiral molecules by using amplified vibrational circular dichroism: from theory to experimental realization

Recent experimental observations of enhanced vibrational circular dichroism (VCD) in molecular systems with low-lying electronically excited states suggest interesting new applications of VCD spectroscopy. The theory describing VCD enhancement through vibronic coupling schemes was derived by Nafie in 1983, but only recently experimental evidence of VCD amplification has demonstrated the extent to which this effect can be exploited as a structure elucidation tool to probe local structure. In this Concept paper, we give an overview of the physics behind vibrational circular dichroism, in particular the equations governing the VCD amplification effect, and review the latest experimental developments with a prospective view on the application of amplified VCD to locally probe biomolecular structure.

Systems biology meets stress ecology: linking molecular and organismal stress responses in Daphnia magna

Background: Ibuprofen and other nonsteroidal anti-inflammatory drugs have been designed to interrupt eicosanoid metabolism in mammals, but little is known of how they affect nontarget organisms. Here we report a systems biology study that simultaneously describes the transcriptomic and phenotypic stress responses of the model crustacean Daphnia magna after exposure to ibuprofen. Results: Our findings reveal intriguing similarities in the mode of action of ibuprofen between vertebrates and invertebrates, and they suggest that ibuprofen has a targeted impact on reproduction at the molecular, organismal, and population level in daphnids. Microarray expression and temporal real-time quantitative PCR profiles of key genes suggest early ibuprofen interruption of crustacean eicosanoid metabolism, which appears to disrupt signal transduction affecting juvenile hormone metabolism and oogenesis. Conclusion: Combining molecular and organismal stress responses provides a guide to possible chronic consequences of environmental stress for population health. This could improve current environmental risk assessment by providing an early indication of the need for higher tier testing. Our study demonstrates the advantages of a systems approach to stress ecology, in which Daphnia will probably play a major role.

Linear viscoelasticity from molecular dynamics simulation of entangled polymers

The linear viscoelastic (LVE) spectrum is one of the primary fingerprints of polymer solutions and melts, carrying information about most relaxation processes in the system. Many single chain theories and models start with predicting the LVE spectrum to validate their assumptions. However, until now, no reliable linear stress relaxation data were available from simulations of multichain systems. In this work, we propose a new efficient way to calculate a wide variety of correlation functions and mean-square displacements during simulations without significant additional CPU cost. Using this method, we calculate stress−stress autocorrelation functions for a simple bead−spring model of polymer melt for a wide range of chain lengths, densities, temperatures, and chain stiffnesses. The obtained stress−stress autocorrelation functions were compared with the single chain slip−spring model in order to obtain entanglement related parameters, such as the plateau modulus or the molecular weight between entanglements. Then, the dependence of the plateau modulus on the packing length is discussed. We have also identified three different contributions to the stress relaxation:  bond length relaxation, colloidal and polymeric. Their dependence on the density and the temperature is demonstrated for short unentangled systems without inertia.

Vesicular systems for delivering conventional small organic molecules and larger macromolecules to and through human skin

The history of using vesicular systems for drug delivery to and through skin started nearly three decades ago with a study utilizing phospholipid liposomes to improve skin deposition and reduce systemic effects of triamcinolone acetonide. Subsequently, many researchers evaluated liposomes with respect to skin delivery, with the majority of them recording localized effects and relatively few studies showing transdermal delivery effects. Shortly after this, Transfersomes were developed with claims about their ability to deliver their payload into and through the skin with efficiencies similar to subcutaneous administration. Since these vesicles are ultradeformable, they were thought to penetrate intact skin deep enough to reach the systemic circulation. Their mechanisms of action remain controversial with diverse processes being reported. Parallel to this development, other classes of vesicles were produced with ethanol being included into the vesicles to provide flexibility (as in ethosomes) and vesicles were constructed from surfactants and cholesterol (as in niosomes). Thee ultradeformable vesicles showed variable efficiency in delivering low molecular weight and macromolecular drugs. This article will critically evaluate vesicular systems for dermal and transdermal delivery of drugs considering both their efficacy and potential mechanisms of action.

Dynamic load balancing for the distributed mining of molecular structures

In molecular biology, it is often desirable to find common properties in large numbers of drug candidates. One family of methods stems from the data mining community, where algorithms to find frequent graphs have received increasing attention over the past years. However, the computational complexity of the underlying problem and the large amount of data to be explored essentially render sequential algorithms useless. In this paper, we present a distributed approach to the frequent subgraph mining problem to discover interesting patterns in molecular compounds. This problem is characterized by a highly irregular search tree, whereby no reliable workload prediction is available. We describe the three main aspects of the proposed distributed algorithm, namely, a dynamic partitioning of the search space, a distribution process based on a peer-to-peer communication framework, and a novel receiverinitiated load balancing algorithm. The effectiveness of the distributed method has been evaluated on the well-known National Cancer Institute’s HIV-screening data set, where we were able to show close-to linear speedup in a network of workstations. The proposed approach also allows for dynamic resource aggregation in a non dedicated computational environment. These features make it suitable for large-scale, multi-domain, heterogeneous environments, such as computational grids.

Efficient mining of discriminative molecular fragments

Frequent pattern discovery in structured data is receiving an increasing attention in many application areas of sciences. However, the computational complexity and the large amount of data to be explored often make the sequential algorithms unsuitable. In this context high performance distributed computing becomes a very interesting and promising approach. In this paper we present a parallel formulation of the frequent subgraph mining problem to discover interesting patterns in molecular compounds. The application is characterized by a highly irregular tree-structured computation. No estimation is available for task workloads, which show a power-law distribution in a wide range. The proposed approach allows dynamic resource aggregation and provides fault and latency tolerance. These features make the distributed application suitable for multi-domain heterogeneous environments, such as computational Grids. The distributed application has been evaluated on the well known National Cancer Institute’s HIV-screening dataset.

Mapping and controlling molecular processes: The application of control theory and system identification algorithms to ASOPS and other fast pulse spectroscopies

Asynchronous Optical Sampling (ASOPS) [1,2] and frequency comb spectrometry [3] based on dual Ti:saphire resonators operated in a master/slave mode have the potential to improve signal to noise ratio in THz transient and IR sperctrometry. The multimode Brownian oscillator time-domain response function described by state-space models is a mathematically robust framework that can be used to describe the dispersive phenomena governed by Lorentzian, Debye and Drude responses. In addition, the optical properties of an arbitrary medium can be expressed as a linear combination of simple multimode Brownian oscillator functions. The suitability of a range of signal processing schemes adopted from the Systems Identification and Control Theory community for further processing the recorded THz transients in the time or frequency domain will be outlined [4,5]. Since a femtosecond duration pulse is capable of persistent excitation of the medium within which it propagates, such approach is perfectly justifiable. Several de-noising routines based on system identification will be shown. Furthermore, specifically developed apodization structures will be discussed. These are necessary because due to dispersion issues, the time-domain background and sample interferograms are non-symmetrical [6-8]. These procedures can lead to a more precise estimation of the complex insertion loss function. The algorithms are applicable to femtosecond spectroscopies across the EM spectrum. Finally, a methodology for femtosecond pulse shaping using genetic algorithms aiming to map and control molecular relaxation processes will be mentioned.

Recent advances in the application of in vitro systems to plant improvement

The horticultural industry was instrumental in the early development and exploitation of genetic techniques over a century ago. This review will describe recent advances in a range of in vitro methods and their application to plant breeding, with especial emphasis on horticultural crops. These methods include improvements in the efficiency of haploid breeding techniques in many fruit and vegetable species using either microspore-derived or ovule-derived plants. Significant molecular information is now available to supplement these essentially empirical approaches and this may enable the more predictable application of these technologies in previously intransigent crops. Similarly there are now improved techniques for isolation of somatic hybrids, by application of either in vitro fertilisation or the culture of excised ovules from interspecific crosses. In addition to examples taken from the traditional scientific literature, emphasis will also be given to the use of patent databases as a valuable source of information on recent novel technologies developed in the commercial world.