Molecular modelling of human CYP1B1 substrate interactions and investigation of allelic variant effects on metabolism

Molecular modelling of human CYP1B1 based on homology with the mammalian P450, CYP2C5, of known three-dimensional structure is reported. The enzyme model has been used to investigate the likely mode of binding for selected CYP1B1 substrates, particularly with regard to the possible effects of allelic variants of CYP1B1 on metabolism. In general, it appears that the CYP1B1 model is consistent with known substrate selectivity for the enzyme, and the sites of metabolism can be rationalized in terms of specific contacts with key amino acid residues within the CYP1B1 heme locus. Further-more, a mode of binding interaction for the inhibitor, a-naphthoflavone, is presented which accords with currently available information. The current paper shows that a combination of molecular modelling and experimental determinations on the substrate metabolism for CYP1B1 allelic variants can aid in the understanding of structure-function relationships within P450 enzymes. (C) 2003 Elsevier Science Ireland Ltd. All rights reserved.

Molecular dynamics study of MscL interactions with a curved lipid bilayer

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.

Stimulated Raman adiabatic passage from an atomic to a molecular Bose-Einstein condensate

The process of stimulated Raman adiabatic passage (STIRAP) provides a possible route for the generation of a coherent molecular Bose-Einstein condensate (BEC) from an atomic BEC. We analyze this process in a three-dimensional mean-field theory, including atom-atom interactions and nonresonant intermediate levels. We find that the process is feasible, but at larger Rabi frequencies than anticipated from a crude single-mode lossless analysis, due to two-photon dephasing caused by the atomic interactions. We then identify optimal strategies in STIRAP allowing one to maintain high conversion efficiencies with smaller Rabi frequencies and under experimentally less demanding conditions.

Wolbachia infections of tephritid fruit flies: molecular evidence for five distinct strains in a single host species

Endosymbiotic bacteria of the genus Wolbachia are widespread among arthropods and can induce cytoplasmic incompatibility, thelytokous parthenogenesis, male-killing or feminization in their hosts. Here, we report phylogenetic relationships of Wolbachia in tephritid fruit flies based on wsp gene sequences. We also report, for the first time, five distinct strains of Wolbachia in Bactrocera ascita sp. B. Four of the five Wolbachia strains found in this species were in the same groups as those found in other tephritid fruit flies, suggesting possible horizontal transmission of Wolbachia from other fruit flies into B. ascita sp. B. The unreliability of wsp-specific group primers demonstrated in this study suggests that these primers might be useful only for preliminary identification of Wolbachia. Final determination of group affiliation needs to be verified with wsp sequence data.

Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

Molecular recognition of macrocyclic peptidomimetic inhibitors by HIV-1 protease

High-resolution crystal structures are described for seven macrocycles complexed with HIV-1 protease (HIVPR). The macrocycles possess two amides and an aromatic group within 15-17 membered rings designed to replace N- or C-terminal tripeptides from peptidic inhibitors of HIVPR. Appended to each macrocycle is a transition state isostere and either an acyclic peptide, nonpeptide, or another macrocycle. These cyclic analogues are potent inhibitors of HIVPR, and the crystal structures show them to be structural mimics of acyclic peptides, binding in the active site of HIVPR via the same interactions. Each macrocycle is restrained to adopt a P-strand conformation which is preorganized for protease binding. An unusual feature of the binding of C-terminal macrocyclic inhibitors is the interaction between a positively charged secondary amine and a catalytic aspartate of HIVPR. A bicyclic inhibitor binds similarly through its secondary amine that lies between its component N-terminal and C-terminal macrocycles. In contrast, the corresponding tertiary amine of the N-terminal macrocycles does not interact with the catalytic aspartates. The amine-aspartate interaction induces a 1.5 Angstrom N-terminal translation of the inhibitors in the active site and is accompanied by weakened interactions with a water molecule that bridges the ligand to the enzyme, as well as static disorder in enzyme flap residues. This flexibility may facilitate peptide cleavage and product dissociation during catalysis. Proteases [Aba(67,95)]HIVPR and [Lys(7),Ile(33),Aba(67,95)]- HIVPR used in this work were shown to have very similar crystal structures.

Characterization of the sequential non-covalent and covalent interactions of the antitumour antibiotic hedamycin with double stranded DNA by NMR spectroscopy

Hedamycin, a member of the pluramycin class of antitumour antibiotics, consists of a planar anthrapyrantrione chromophore to which is attached two aminosugar rings at one end and a bisepoxide-containing sidechain at the other end, Binding to double-stranded DNA is known to involve both reversible and non-reversible modes of interaction. As a part of studies directed towards elucidating the structural basis for the observed 5'-pyGT-3' sequence selectivity of hedamycin, we conducted one-dimensional NMR titration experiments at low temperature using the hexadeoxyribonucleotide duplexes d(CACGTG)(2) and d(CGTACG)(2). Spectral changes which occurred during these titrations are consistent with hedamycin initially forming a reversible complex in slow exchange on the NMR timescale and binding through intercalation of the chromophore. Monitoring of this reversible complex over a period of hours revealed a second type of spectral change which corresponds with formation of a non-reversible complex. Copyright (C) 1999 John Wiley & Sons, Ltd.

Multidimensional quantum solitons with nondegenerate parametric interactions: Photonic and Bose-Einstein condensate environments

We consider the quantum theory of three fields interacting via parametric and repulsive quartic couplings. This can be applied to treat photonic chi((2)) and chi((3)) interactions, and interactions in atomic Bose-Einstein condensates or quantum Fermi gases, describing coherent molecule formation together with a-wave scattering. The simplest two-particle quantum solitons or bound-state solutions of the idealized Hamiltonian, without a momentum cutoff, are obtained exactly. They have a pointlike structure in two and three dimensions-even though the corresponding classical theory is nonsingular. We show that the solutions can be regularized with a momentum cutoff. The parametric quantum solitons have much more realistic length scales and binding energies than chi((3)) quantum solitons, and the resulting effects could potentially be experimentally tested in highly nonlinear optical parametric media or interacting matter-wave systems. N-particle quantum solitons and the ground state energy are analyzed using a variational approach. Applications to atomic/molecular Bose-Einstein condensates (BEC's) are given, where we predict the possibility of forming coupled BEC solitons in three space dimensions, and analyze superchemistry dynamics.

Superhyperfine interactions in inhomogeneously broadened paramagnetic centers observed via a hole-burned free induction decay

Superhyperfine interactions in inhomogeneously broadened paramagnetic centers are observed using a single high-turn-angle microwave pulse. The free induction signal that follows the hole-burning pulse exhibits oscillations that are distinct from the oscillatory free induction decay observable in some inhomogeneously broadened systems. It contains frequencies characteristic of the superhyperfine splittings, together with a zero frequency component. Experimental examples of the effect in both orientationally disordered (powdered) and structurally disordered (glassy) systems are presented and compared with the conceptually similar Fourier transform electron paramagnetic resonance detected nuclear magnetic resonance experiment, together with numerical simulations. (C) 2003 American Institute of Physics.

Radiation of the Australian flora: What can comparisons of molecular phylogenies across multiple taxa tell us about the evolution of diversity in present-day communities?

The Australian fossil record shows that from ca. 25 Myr ago, the aseasonal-wet biome (rainforest and wet heath) gave way to the unique Australian sclerophyll biomes dominated by eucalypts, acacias and casuarinas. This transition coincided with tectonic isolation of Australia, leading to cooler, drier, more seasonal climates. From 3 Myr ago, aridification caused rapid opening of the central Australian and zone. Molecular phylogenies with dated nodes have provided new perspectives on how these events could have affected the evolution of the Australian flora. During the Mid-Cenozoic (25-10 Myr ago) period of climatic change, there were rapid radiations in sclerophyll taxa, such as Banksia, eucalypts, pea-flowered legumes and Allocasuarina. At the same time, taxa restricted to the aseasonal-wet biome (Nothofagus, Podocarpaceae and Araucariaceae) did not radiate or were depleted by extinction. During the Pliocene aridification, two Eremean biome taxa (Lepidium and Chenopodiaceae) radiated rapidly after dispersing into Australia from overseas. It is clear that the biomes have different histories. Lineages in the aseasonal-wet biome are species poor, with sister taxa that are species rich, either outside Australia or in the sclerophyll biomes. In conjunction with the fossil record, this indicates depletion of the Australian aseasonal-wet biome from the Mid-Cenozoic. In the sclerophyll biomes, there have been multiple exchanges between the southwest and southeast, rather than single large endemic radiations after a vicariance event. There is need for rigorous molecular phylogenetic studies so that additional questions can be addressed, such as how interactions between biomes may have driven the speciation process during radiations. New studies should include the hither-to neglected monsoonal tropics.

Interfacial interactions and structure of organic-inorganic nanohybrids

Understanding the interfacial interactions and structure is important to better design and application of organic-inorganic nanohybrids. This paper presents our recent molecular dynamic studies on organoclays and polymer nanocomposites, including the layering behavior of organoclays, structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays, and interfacial interactions and structure of polyurethane nanocomposites. The results demonstrate that the layering behaviors of organoclays are closely related to the chain length of quaternary alkyl ammoniums and cation exchangeable capacity of clays. In addition to typical layered structures such as monolayer, bilayer and pseudo-trilayer, a pseudo-quadrilayer structure was also observed in organoclays modified with dioctadecyldimethyl ammoniums (DODDMA). In such a structure, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer or even the next nearest layer. Moreover, the diffusion constants of nitrogen and methylene atoms increase with the temperature and methelene towards the tail groups. For polyurethane nanocomposite, the van der Waals interaction between apolar alkyl chains and soft segments of polyurethane predominates the interactions between organoclay and polyurethane. Different from most bulk polyurethane systems, there is no distinct phase-separated structure for the polyurethane.

Giant protein kinases: Domain interactions and structural basis of autoregulation

The myosin-associated giant protein kinases twitchin and titin are composed predominantly of fibronectin- and immunoglobulin-like modules, We report the crystal structures of two autoinhibited twitchin kinase fragments, one from Aplysia and a larger fragment from Caenorhabditis elegans containing an additional C-terminal immunoglobulin-like domain, The structure of the longer fragment shoes that the immunoglobulin domain contacts the protein kinase domain on the opposite side from the catalytic cleft, laterally exposing potential myosin binding residues, Together, the structures reveal the cooperative interactions between the autoregulatory region and the residues from the catalytic domain involved in protein substrate binding, ATP binding, catalysis and the activation loop, and explain the differences between the observed autoinhibitory mechanism and the one found in the structure of calmodulin-dependent kinase I.

Molecular motion in miscible polymer blends .1. Motion in blends of PEO and PVPh studied by solid-state C-13 T-1 rho measurements

Changes in molecular motion in blends of PEO-PVPh have been studied using measurements of C-13 T-1 rho relaxation times. C-13 T-1 rho relaxation has been confirmed as arising from spin-lattice interactions by observation of the variation in T-1 rho with rf field strength and temperature. In the pure homopolymers a minimum in T-1 rho is observed at ca. 50 K above the glass transition temperatures detected by DSC. After blending, the temperature of the minimum in T-1 rho for PEO increased, while that for PVPh decreased, however, the minima, which correspond to the temperatures where the average correlation times for reorientation are close to 3.1 mu s, are separated by 45 K (in a 45% PEO-PVPh blend). These phenomena are explained in terms of the local nature of T-1 rho measurements. The motions of the individual homopolymer chains are only partially coupled in the blend. A short T-1 rho has been observed for protonated aromatic carbons, and assigned to phenyl rings undergoing large-angle oscillatory motion, The effects of blending, and temperature, on the proportion of rings undergoing oscillatory motion are analyzed.

Concentration dependence of sodium permeation and sodium ion interactions in the cyclic AMP-gated channels of mammalian olfactory receptor neurons

The dependence of currents through the cyclic nucleotide-gated (CNG) channels of mammalian olfactory receptor neurons (ORNs) on the concentration of NaCl was studied in excised inside-out patches from their dendritic knobs using the patch-clamp technique. With a saturating concentration (100 mu M) of adenosine 3', 5'-cyclic monophosphate (cAMP), the changes in the reversal potential of macroscopic currents were studied at NaCl concentrations from 25 to 300 mM. In symmetrical NaCl solutions without the addition of divalent cations, the current-voltage relations were almost linear, reversing close to O mV. When the external NaCl concentration was maintained at 150 mM and the internal concentrations were varied, the reversal potentials of the cAMP-activated currents closely followed the Na+ equilibrium potential indicating that P-Cl/P-Na approximate to 0. However, at low external NaCl concentrations (less than or equal to 100 mM) there was some significant chloride permeability. Our results further indicated that Na+ currents through these channels: (i) did not obey the independence principle; (ii) showed saturation kinetics with K(m)s in the range of 100-150 mM and (iii) displayed a lack of voltage dependence of conductance in asymmetric solutions that suggested that ion-binding sites were situated midway along the channel. Together, these characteristics indicate that the permeation properties of the olfactory CNG channels are significantly different from those of photoreceptor CNG channels.