Molecular phylogenetic analysis of evolutionary trends in stonefly wing structure and locomotor behavior

Insects in the order Plecoptera (stoneflies) use a form of two-dimensional aerodynamic locomotion called surface skimming to move across water surfaces. Because their weight is supported by water, skimmers can achieve effective aerodynamic locomotion even with small wings and weak flight muscles. These mechanical features stimulated the hypothesis that surface skimming may have been an intermediate stage in the evolution of insect flight, which has perhaps been retained in certain modern stoneflies. Here we present a phylogeny of Plecoptera based on nucleotide sequence data from the small subunit rRNA (18S) gene. By mapping locomotor behavior and wing structural data onto the phylogeny, we distinguish between the competing hypotheses that skimming is a retained ancestral trait or, alternatively, a relatively recent loss of flight. Our results show that basal stoneflies are surface skimmers, and that various forms of surface skimming are distributed widely across the plecopteran phylogeny. Stonefly wings show evolutionary trends in the number of cross veins and the thickness of the cuticle of the longitudinal veins that are consistent with elaboration and diversification of flight-related traits. These data support the hypothesis that the first stoneflies were surface skimmers, and that wing structures important for aerial flight have become elaborated and more diverse during the radiation of modern stoneflies.

Large vortex-like structure of dipole field in computer models of liquid water and dipole-bridge between biomolecules

We propose a framework to describe the cooperative orientational motions of water molecules in liquid water and around solute molecules in water solutions. From molecular dynamics (MD) simulation a new quantity “site-dipole field” is defined as the averaged orientation of water molecules that pass through each spatial position. In the site-dipole field of bulk water we found large vortex-like structures of more than 10 Å in size. Such coherent patterns persist more than 300 ps although the orientational memory of individual molecules is quickly lost. A 1-ns MD simulation of systems consisting of two amino acids shows that the fluctuations of site-dipole field of solvent are pinned around the amino acids, resulting in a stable dipole-bridge between side-chains of amino acids. The dipole-bridge is significantly formed even for the side-chain separation of 14 Å, which corresponds to five layers of water. The way that dipole-bridge forms sensitively depends on the side-chain orientations and thereby explains the specificity in the solvent-mediated interactions between biomolecules.

Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water

The hydrophobic interaction, the tendency for nonpolar molecules to aggregate in solution, is a major driving force in biology. In a direct approach to the physical basis of the hydrophobic effect, nanosecond molecular dynamics simulations were performed on increasing numbers of hydrocarbon solute molecules in water-filled boxes of different sizes. The intermittent formation of solute clusters gives a free energy that is proportional to the loss in exposed molecular surface area with a constant of proportionality of 45 ± 6 cal/mol⋅Å2. The molecular surface area is the envelope of the solute cluster that is impenetrable by solvent and is somewhat smaller than the more traditional solvent-accessible surface area, which is the area transcribed by the radius of a solvent molecule rolled over the surface of the cluster. When we apply a factor relating molecular surface area to solvent-accessible surface area, we obtain 24 cal/mol⋅Å2. Ours is the first direct calculation, to our knowledge, of the hydrophobic interaction from molecular dynamics simulations; the excellent qualitative and quantitative agreement with experiment proves that simple van der Waals interactions and atomic point-charge electrostatics account for the most important driving force in biology.

The arginine finger of RasGAP helps Gln-61 align the nucleophilic water in GAP-stimulated hydrolysis of GTP

The Ras family of GTPases is a collection of molecular switches that link receptors on the plasma membrane to signaling pathways that regulate cell proliferation and differentiation. The accessory GTPase-activating proteins (GAPs) negatively regulate the cell signaling by increasing the slow intrinsic GTP to GDP hydrolysis rate of Ras. Mutants of Ras are found in 25–30% of human tumors. The most dramatic property of these mutants is their insensitivity to the negative regulatory action of GAPs. All known oncogenic mutants of Ras map to a small subset of amino acids. Gln-61 is particularly important because virtually all mutations of this residue eliminate sensitivity to GAPs. Despite its obvious importance for carcinogenesis, the role of Gln-61 in the GAP-stimulated GTPase activity of Ras has remained a mystery. Our molecular dynamics simulations of the p21ras–p120GAP–GTP complex suggest that the local structure around the catalytic region can be different from that revealed by the x-ray crystal structure. We find that the carbonyl oxygen on the backbone of the arginine finger supplied in trans by p120GAP (Arg-789) interacts with a water molecule in the active site that is forming a bridge between the NH2 group of the Gln-61 and the γ-phosphate of GTP. Thus, Arg-789 may play a dual role in generating the nucleophile as well as stabilizing the transition state for P—O bond cleavage.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

The human AQP4 gene: definition of the locus encoding two water channel polypeptides in brain.

The aquaporin family of membrane water transport proteins are expressed in diverse tissues, and in brain the predominant water channel protein is AQP4. Here we report the isolation and characterization of the human AQP4 cDNAs and genomic DNA. Two cDNAs were isolated corresponding to the two initiating methionines (M1 in a 323-aa polypeptide and M23 in a 301-aa polypeptide) previously identified in rat [Jung, J.S., Bhat, R.V., Preston, G.M., Guggino, W.B. & Agre, P. (1994) Proc. Natl. Acad. Sci. USA 91, 13052-13056]. Similar to other aquaporins, the AQP4 gene is composed of four exons encoding 127, 55, 27, and 92 amino acids separated by introns of 0.8, 0.3, and 5.2 kb. Unlike other aquaporins, an alternative coding initiation sequence (designated exon 0) was located 2.7 kb upstream of exon 1. When spliced together, M1 and the subsequent 10 amino acids are encoded by exon 0; the next 11 amino acids and M23 are encoded by exon 1. Transcription initiation sites have been mapped in the proximal promoters of exons 0 and 1. RNase protection revealed distinct transcripts corresponding to M1 and M23 mRNAs, and AQP4 immunoblots of cerebellum demonstrated reactive polypeptides of 31 and 34 kDa. Using a P1 and a lambda EMBL subclone, the chromosomal site of the human AQP4 gene was mapped to chromosome 18 at the junction of q11.2 and q12.1 by fluorescence in situ hybridization. These studies may now permit molecular characterization of AQP4 during human development and in clinical disorders.

Glass transition in DNA from molecular dynamics simulations.

Molecular dynamics simulations of the oligonucleotide duplex d(CGCGCG)2 in aqueous solution are used to investigate the glass transition phenomenon. The simulations were performed at temperatures in the 20 K to 340 K range. The mean square atomic fluctuations showed that the behavior of the oligonucleotide duplex was harmonic at low temperatures. A glass transition temperature at 223 K to 234 K was inferred for the oligonucleotide duplex, which is in agreement with experimental observations. The largest number of hydrogen bounds between the polar atoms of the oligonucleotide duplex and the water molecules was obtained at the glass transition temperature. With increasing temperature we observed a decrease in the average lifetime of the hydrogen bonds to water molecules.

Crystal and molecular structure of paclitaxel (taxol).

Paclitaxel (formerly called taxol), an important anticancer drug, inhibits cell replication by binding to and stabilizing microtubule polymers. As drug-receptor interactions are governed by the three-dimensional stereochemistries of both participants, we have determined the crystal structure of paclitaxel to identify its conformational preferences that may be related to biological activity. The monoclinic crystals contain two independent paclitaxel molecules in the asymmetric unit plus several water and dioxane solvent molecules. Taxane ring conformation is very similar in both paclitaxel molecules and is similar to the taxane ring conformation found in the crystal structure of the paclitaxel analogue docetaxel (formerly called taxotere). The two paclitaxel molecules have carbon-13 side-chain conformations that differ from each other and from that of the corresponding side chain in the docetaxel crystal structure. The carbon-13 side-chain conformation of one paclitaxel molecule is similar to what was proposed from NMR studies done in polar solvents, while that of the other paclitaxel molecule is different and hitherto unobserved. The paclitaxel molecules interact with each other and with solvent atoms through an extensive network of hydrogen bonds. Analysis of the hydrogen-bonding network together with structure-activity studies may suggest which atoms of paclitaxel are important for binding to microtubule receptors.

Immunolocalization of the mercurial-insensitive water channel and glycerol intrinsic protein in epithelial cell plasma membranes.

Two water channel homologs were cloned recently from rat kidney, mercurial-insensitive water channel (MIWC) and glycerol intrinsic protein (GLIP). Polyclonal antibodies were raised against synthetic C-terminal peptides and purified by affinity chromatography. MIWC and GLIP antibodies recognized proteins in rat kidney with an apparent molecular mass of 30 and 27 kDa, respectively, and did not cross-react. By immunofluorescence, MIWC and GLIP were expressed together on the basolateral plasma membrane of collecting duct principal cells in kidney. By immunohistochemistry, MIWC and GLIP were expressed on tracheal epithelial cells with greater expression of GLIP on the basal plasma membrane and MIWC on the lateral membrane; only MIWC was expressed in bronchial epithelia. In eye, GLIP was expressed in conjunctival epithelium, whereas MIWC was found in iris, ciliary body, and neural cell layers in retina. MIWC and GLIP colocalized on the basolateral membrane of villus epithelial cells in colon and brain ependymal cells. Expression of MIWC and GLIP was not detected in small intestine, liver, spleen, endothelia, and cells that express water channels CHIP28 or WCH-CD. These studies suggest water/solute transporting roles for MIWC and GLIP in the urinary concentrating mechanism, cerebrospinal fluid absorption, ocular fluid balance, fecal dehydration, and airway humidification. The unexpected membrane colocalization of MIWC and GLIP in several tissues suggests an interaction at the molecular and/or functional levels.

Interfacial activation-based molecular bioimprinting of lipolytic enzymes.

Interfacial activation-based molecular (bio)-imprinting (IAMI) has been developed to rationally improve the performance of lipolytic enzymes in nonaqueous environments. The strategy combinedly exploits (i) the known dramatic enhancement of the protein conformational rigidity in a water-restricted milieu and (ii) the reported conformational changes associated with the activation of these enzymes at lipid-water interfaces, which basically involves an increased substrate accessibility to the active site and/or an induction of a more competent catalytic machinery. Six model enzymes have been assayed in several model reactions in nonaqueous media. The results, rationalized in light of the present biochemical and structural knowledge, show that the IAMI approach represents a straightforward, versatile method to generate manageable, activated (kinetically trapped) forms of lipolytic enzymes, providing under optimal conditions nonaqueous rate enhancements of up to two orders of magnitude. It is also shown that imprintability of lipolytic enzymes depends not only on the nature of the enzyme but also on the "quality" of the interface used as the template.

Heterogeneity in molecular recognition by restriction endonucleases: osmotic and hydrostatic pressure effects on BamHI, Pvu II, and EcoRV specificity.

The cleavage specificity of the Pvu II and BamHI restriction endonucleases is found to be dramatically reduced at elevated osmotic pressure. Relaxation in specificity of these otherwise highly accurate and specific enzymes, previously termed "star activity," is uniquely correlated with osmotic pressure between 0 and 100 atmospheres. No other colligative solvent property exhibits a uniform correlation with star activity for all of the compounds tested. Application of hydrostatic pressure counteracts the effects of osmotic pressure and restores the natural selectivity of the enzymes for their canonical recognition sequences. These results indicate that water solvation plays an important role in the site-specific recognition of DNA by many restriction enzymes. Osmotic pressure did not induce an analogous effect on the specificity of the EcoRV endonuclease, implying that selective hydration effects do not participate in DNA recognition in this system. Hydrostatic pressure was found to have little effect on the star activity induced by changes in ionic strength, pH, or divalent cation, suggesting that distinct mechanisms may exist for these observed alterations in specificity. Recent evidence has indicated that BamHI and EcoRI share similar structural motifs, while Pvu II and EcoRV belong to a different structural family. Evidently, the use of hydration water to assist in site-specific recognition is a motif neither limited to nor defined by structural families.

Comportamento tribo-mecânico e desgaste adesivo de materiais em nanoescala: análises por dinâmica molecular e mecânica do contínuo.

As formulações baseadas na mecânica do contínuo, embora precisas até certo ponto, por vezes não podem ser utilizadas, ou não são conceitualmente corretas para o entendimento de fenômenos em escalas reduzidas. Estas limitações podem aparecer no estudo dos fenômenos tribológicos em escala nanométrica, que passam a necessitar de novos métodos experimentais, teóricos e computacionais que permitam explorar estes fenômenos com a resolução necessária. Simulações atomísticas são capazes de descrever fenômenos em pequena escala, porém, o número necessário de átomos modelados e, portanto, o custo computacional - geralmente torna-se bastante elevado. Por outro lado, os métodos de simulação associados à mecânica do contínuo são mais interessantes em relação ao custo computacional, mas não são precisos na escala atômica. A combinação entre essas duas abordagens pode, então, permitir uma compreensão mais realista dos fenômenos da tribologia. Neste trabalho, discutem-se os conceitos básicos e modelos de atrito em escala atômica e apresentam-se estudos, por meio de simulação numérica, para a análise e compreensão dos mecanismos de atrito e desgaste no contato entre materiais. O problema é abordado em diferentes escalas, e propõe-se uma abordagem conjunta entre a Mecânica do Contínuo e a Dinâmica Molecular. Para tanto, foram executadas simulações numéricas, com complexidade crescente, do contato entre superfícies, partindo-se de um primeiro modelo que simula o efeito de defeitos cristalinos no fenômeno de escorregamento puro, considerando a Dinâmica Molecular. Posteriormente, inseriu-se, nos modelos da mecânica do contínuo, considerações sobre o fenômeno de adesão. A validação dos resultados é feita pela comparação entre as duas abordagens e com a literatura.

Novel silica membrane material for molecular sieve applications

Development of new silica membranes properties, e.g., molecular sieving properties, has been increasingly gaining importance in the last few years. A novel unsupported silica membrane, referred to as hydrophobic metal-doped silica, was developed by cobalt-doping within the organic templated silica matrix. The novel material was prepared by the acid-catalyzed hydrolysis and condensation process of tetraethylorthosilicate (TEOS) and methyltriethoxysilane (MTES), which is the precursor for methyl ligand covalently bounded to the silica matrix. The synthesis and surface properties of the novel unsupported silica membrane as well as the unsupported blank silica and modified silica membranes were revealed by surface and microstructural techniques, such as water contact angle measurement, FTIR, X-ray, Solid-state 29Si MAS NMR, TGA and N2 and CO2 adsorption measurements. The results showed that the thermal stability of the organic templated silica matrix was enhanced by cobalt-doping process. A hydrophobic microporous silica membrane material with high thermal stability up to ∼560 °C in oxidizing atmosphere and a narrow pore size distribution centered at 1.1 nm was obtained. Therefore, a novel precursor material for molecular sieve silica membranes applications has been achieved and developed.

Water adsorption in hydrophilic zeolites: experiment and simulation

We have measured experimental adsorption isotherms of water in zeolite LTA4A, and studied the regeneration process by performing subsequent adsorption cycles after degassing at different temperatures. We observed incomplete desorption at low temperatures, and cation rearrangement at successive adsorption cycles. We also developed a new molecular simulation force field able to reproduce experimental adsorption isotherms in the range of temperatures between 273 K and 374 K. Small deviations observed at high pressures are attributed to the change in the water dipole moment at high loadings. The force field correctly describes the preferential adsorption sites of water at different pressures. We tested the influence of the zeolite structure, framework flexibility, and cation mobility when considering adsorption and diffusion of water. Finally, we performed checks on force field transferability between different hydrophilic zeolite types, concluding that classical, non-polarizable water force fields are not transferable.

New Approaches for Teaching Soil and Rock Mechanics Using Information and Communication Technologies

Soil and rock mechanics are disciplines with a strong conceptual and methodological basis. Initially, when engineering students study these subjects, they have to understand new theoretical phenomena, which are explained through mathematical and/or physical laws (e.g. consolidation process, water flow through a porous media). In addition to the study of these phenomena, students have to learn how to carry out estimations of soil and rock parameters in laboratories according to standard tests. Nowadays, information and communication technologies (ICTs) provide a unique opportunity to improve the learning process of students studying the aforementioned subjects. In this paper, we describe our experience of the incorporation of ICTs into the classical teaching-learning process of soil and rock mechanics and explain in detail how we have successfully developed various initiatives which, in summary, are: (a) implementation of an online social networking and microblogging service (using Twitter) for gradually sending key concepts to students throughout the semester (gradual learning); (b) detailed online virtual laboratory tests for a delocalized development of lab practices (self-learning); (c) integration of different complementary learning resources (e.g. videos, free software, technical regulations, etc.) using an open webpage. The complementary use to the classical teaching-learning process of these ICT resources has been highly satisfactory for students, who have positively evaluated this new approach.