The development of retinotectal maps: A review of models based on molecular gradients

Information about the world is often represented in the brain in the form of topographic maps. A paradigm example is the topographic representation of the visual world in the optic tectum/superior colliculus. This map initially forms during neural development using activity-independent molecular cues, most notably some type of chemospecific matching between molecular gradients in the retina and corresponding gradients in the tectum/superior colliculus. Exactly how this process might work has been studied both experimentally and theoretically for several decades. This review discusses the experimental data briefly, and then in more detail the theoretical models proposed. The principal conclusions are that (1) theoretical models have helped clarify several important ideas in the field, (2) earlier models were often more sophisticated than more recent models, and (3) substantial revisions to current modelling approaches are probably required to account for more than isolated subsets of the experimental data.

Algebraic Bethe ansatz method for the exact calculation of energy spectra and form factors: applications to models of Bose-Einstein condensates and metallic nanograins

In this review we demonstrate how the algebraic Bethe ansatz is used for the calculation of the-energy spectra and form factors (operator matrix elements in the basis of Hamiltonian eigenstates) in exactly solvable quantum systems. As examples we apply the theory to several models of current interest in the study of Bose-Einstein condensates, which have been successfully created using ultracold dilute atomic gases. The first model we introduce describes Josephson tunnelling between two coupled Bose-Einstein condensates. It can be used not only for the study of tunnelling between condensates of atomic gases, but for solid state Josephson junctions and coupled Cooper pair boxes. The theory is also applicable to models of atomic-molecular Bose-Einstein condensates, with two examples given and analysed. Additionally, these same two models are relevant to studies in quantum optics; Finally, we discuss the model of Bardeen, Cooper and Schrieffer in this framework, which is appropriate for systems of ultracold fermionic atomic gases, as well as being applicable for the description of superconducting correlations in metallic grains with nanoscale dimensions.; In applying all the above models to. physical situations, the need for an exact analysis of small-scale systems is established due to large quantum fluctuations which render mean-field approaches inaccurate.

Reactivity of chars and carbons: New insights through molecular modeling

A new model proposed for the gasification of chars and carbons incorporates features of the turbostratic nanoscale structure that exists in such materials. The model also considers the effect of initial surface chemistry and different reactivities perpendicular to the edges and to the faces of the underlying crystallite planes comprising the turbostratic structure. It may be more realistic than earlier models based on pore or grain structure idealizations when the carbon contains large amounts of crystallite matter. Shrinkage of the carbon particles in the chemically controlled regime is also possible due to the random complete gasification of crystallitic planes. This mechanism can explain observations in the literature of particle size reduction. Based on the model predictions, both initial surface chemistry and the number of stacked planes in the crystallites strongly influence the reactivity and particle shrinkage. Its test results agree well with literature data on the air-oxidation of Spherocarb and show that it accurately predicts the variation of particle size with conversion. Model parameters are determined entirely from rate measurements.

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.

Hemichordates and deuterostome evolution: robust molecular phylogenetic support for a hemichordate plus echinoderm clade

Hemichordates were traditionally allied to the chordates, but recent molecular analyses have suggested that hemichordates are a sister group to the echinoderms, a relationship that has important consequences for the interpretation of the evolution of deuterostome body plans. However, the molecular phylogenetic analyses to date have not provided robust support for the hemichordate + echinoderm clade. We use a maximum likelihood framework, including the parametric bootstrap, to reanalyze DNA data from complete mitochondrial genomes and nuclear 18S rRNA. This approach provides the first statistically significant support for the hemichordate + echinoderm clade from molecular data. This grouping implies that the ancestral deuterostome had features that included an adult with a pharynx and a dorsal nerve cord and an indirectly developing dipleurula-like larva.

Efficient discovery of immune response targets by cyclical refinement of QSAR models of peptide binding

Peptides that induce and recall T-cell responses are called T-cell epitopes. T-cell epitopes may be useful in a subunit vaccine against malaria. Computer models that simulate peptide binding to MHC are useful for selecting candidate T-cell epitopes since they minimize the number of experiments required for their identification. We applied a combination of computational and immunological strategies to select candidate T-cell epitopes. A total of 86 experimental binding assays were performed in three rounds of identification of HLA-All binding peptides from the six preerythrocytic malaria antigens. Thirty-six peptides were experimentally confirmed as binders. We show that the cyclical refinement of the ANN models results in a significant improvement of the efficiency of identifying potential T-cell epitopes. (C) 2001 by Elsevier Science Inc.

Virtual models of the HLA class I antigen processing pathway

Antigen recognition by cytotoxic CD8 T cells is dependent upon a number of critical steps in MHC class I antigen processing including proteosomal cleavage, TAP transport into the endoplasmic reticulum, and MHC class 1 binding. Based on extensive experimental data relating to each of these steps there is now the capacity to model individual antigen processing steps with a high degree of accuracy. This paper demonstrates the potential to bring together models of individual antigen processing steps, for example proteosome cleavage, TAP transport, and MHC binding, to build highly informative models of functional pathways. In particular, we demonstrate how an artificial neural network model of TAP transport was used to mine a HLA-binding database so as to identify H LA-binding peptides transported by TAP. This integrated model of antigen processing provided the unique insight that HLA class I alleles apparently constitute two separate classes: those that are TAP-efficient for peptide loading (HLA-B27, -A3, and -A24) and those that are TAP-inefficient (HLA-A2, -B7, and -B8). Hence, using this integrated model we were able to generate novel hypotheses regarding antigen processing, and these hypotheses are now capable of being tested experimentally. This model confirms the feasibility of constructing a virtual immune system, whereby each additional step in antigen processing is incorporated into a single modular model. Accurate models of antigen processing have implications for the study of basic immunology as well as for the design of peptide-based vaccines and other immunotherapies. (C) 2004 Elsevier Inc. All rights reserved.

Recent advances in processing and characterization of periodic mesoporous MCM-41 silicate molecular sieves

The discovery of periodic mesoporous MCM-41 and related molecular sieves has attracted significant attention from a fundamental as well as applied perspective. They possess well-defined cylindrical/hexagonal mesopores with a simple geometry, tailored pore size, and reproducible surface properties. Hence, there is an ever-growing scientific interest in the challenges posed by their processing and characterization and by the refinement of various sorption models. Further, MCM-41-based materials are currently under intense investigation with respect to their utility as adsorbents, catalysts, supports, ion-exchangers, and molecular hosts. In this article, we provide a critical review of the developments in these areas with particular emphasis on adsorption characteristics, progress in controlling the pore sizes, and a comparison of pore size distributions using traditional and newer models. The model proposed by the authors for adsorption isotherms and criticalities in capillary condensation and hysteresis is found to explain unusual adsorption behavior in these materials while providing a convenient characterization tool.

Molecular genetics of schizophrenia

1. Schizophrenia is a chronic, disabling brain disease that affects approxmately 1% of the world's population. It is characterized by delusions, hallucinations and formal thought disorder, together with a decline in socio-occupational functioning. While the causes for schizophrenia remain unknown, evidence from family, twin and adoption studies clearly demonstrates that it aggregates in families, with this clustering largely attributable to genetic rather than cultural or environmental factors. Identifying the genes involved, however, has proven to be a difficult task because schizophrenia is a complex trait characterized by an imprecise phenotype, the existence of phenocopies and the presence of low disease penetrance, 2. The current working hypothesis for schizophrenia causation is that multiple genes of small to moderate effect confer compounding risk through interactions with each other and with non-genetic risk factors, The same genes may be commonly involved in conferring risk across populations or they may vary in number and strength between different populations. To search for evidence of such genetic loci, both candidate gene and genome-wide linkage studies have been used in clinical cohorts collected from a variety of populations. Collectively, these works provide some evidence for the involvement of a number of specific genes (e.g. the 5-hydroxytryptamine (5-HT) type 2a receptor (5-HT2a) gene and the dopamine D-3 receptor gene) and as yet unidentified factors localized to specific chromosomal regions, including 6p, 6q, 8p, 13q and 22q, These data provide suggestive, but no conclusive, evidence for causative genes. 3. To enable further progress there is a need to: (i) collect fine-grained clinical datasets while searching the schizophrenia phenotype for subgroups or dimensions that may provide a more direct route to causative genes; and (ii) integrate recent refinements in molecular genetic technology, including modern composite marker maps, DNA expression assays and relevant animal models, while using the latest analytical techniques to extract maximum information in order to help distinguish a true result from a false-positive finding.

Hemoglobin-degrading, aspartic proteases of blood-feeding parasites - Substrate specificity revealed by homology models

Blood-feeding parasites, including schistosomes, hookworms, and malaria parasites, employ aspartic proteases to make initial or early cleavages in ingested host hemoglobin. To better understand the substrate affinity of these aspartic proteases, sequences were aligned with and/or three-dimensional, molecular models were constructed of the cathepsin D-like aspartic proteases of schistosomes and hookworms and of plasmepsins of Plasmodium falciparum and Plasmodium vivax, using the structure of human cathepsin D bound to the inhibitor pepstatin as the template. The catalytic subsites S5 through S4' were determined for the modeled parasite proteases. Subsequently, the crystal structure of mouse renin complexed with the nonapeptidyl inhibitor t-butyl-CO-His-Pro-Phe-His-Leu [CHOHCH2]Leu-Tyr-Tyr-Ser-NH2 (CH-66) was used to build homology models of the hemoglobin-degrading peptidases docked with a series of octapeptide substrates. The modeled octapeptides included representative sites in hemoglobin known to be cleaved by both Schistosoma japonicum cathepsin D and human cathepsin D, as well as sites cleaved by one but not the other of these enzymes. The peptidase-octapeptide substrate models revealed that differences in cleavage sites were generally attributable to the influence of a single amino acid change among the P5 to P4' residues that would either enhance or diminish the enzymatic affinity. The difference in cleavage sites appeared to be more profound than might be expected from sequence differences in the enzymes and hemoglobins. The findings support the notion that selective inhibitors of the hemoglobin-degrading peptidases of blood-feeding parasites at large could be developed as novel anti-parasitic agents.

Reconciling paleodistribution models and comparative phylogeography in the Wet Tropics rainforest land snail Gnarosophia bellendenkerensis (Brazier 1875)

Comparative phylogeography has proved useful for investigating biological responses to past climate change and is strongest when combined with extrinsic hypotheses derived from the fossil record or geology. However, the rarity of species with sufficient, spatially explicit fossil evidence restricts the application of this method. Here, we develop an alternative approach in which spatial models of predicted species distributions under serial paleoclimates are compared with a molecular phylogeography, in this case for a snail endemic to the rainforests of North Queensland, Australia. We also compare the phylogeography of the snail to those from several endemic vertebrates and use consilience across all of these approaches to enhance biogeographical inference for this rainforest fauna. The snail mtDNA phylogeography is consistent with predictions from paleoclimate modeling in relation to the location and size of climatic refugia through the late Pleistocene-Holocene and broad patterns of extinction and recolonization. There is general agreement between quantitative estimates of population expansion from sequence data (using likelihood and coalescent methods) vs. distributional modeling. The snail phylogeography represents a composite of both common and idiosyncratic patterns seen among vertebrates, reflecting the geographically finer scale of persistence and subdivision in the snail. In general, this multifaceted approach, combining spatially explicit paleoclimatological models and comparative phylogeography, provides a powerful approach to locating historical refugia and understanding species' responses to them.

Effects of molecular crowding by saccharides on alpha-chymotrypsin dimerization

Given the importance of protein complexes as therapeutic targets, it is necessary to understand the physical chemistry of these interactions under the crowded conditions that exist in cells. We have used sedimentation equilibrium to quantify the enhancement of the reversible homodimerization of alpha-chymotrypsin by high concentrations of the osmolytes glucose, sucrose, and raffinose. In an attempt to rationalize the ostuolyte-mediated stabilization of the a-chymotrypsin homodimer, we have used models based on binding interactions (transfer-free energy analysis) and steric interactions (excluded volume theory) to predict the stabilization. Although transfer-free energy analysis predicts reasonably well the relatively small stabilization observed for complex formation between cytochrome c and cytochrome c peroxidase, as well as that between bobtail quail lysozyme and a monoclonal Fab fragment, it underestimates the sugar-mediated stabilization of the alpha-chymotrypsin dimer. Although predictions based on excluded volume theory overestimate the stabilization, it would seem that a major determinant in the observed stabilization of the a-chymotrypsin homodimer is the thermodynamic nonideality arising from molecular crowding by the three small sugars.

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.