Molecular Dynamics Study on Protein and it's Water Structure at High Pressure

Abstract

We have performed NPT molecular dynamics simulations (Langevin Piston Method) on two types of solvated proteins-‘denaturation-unfavorable’ protein (insulin) and ‘denaturation-favorable protein’ (ribonuclease A) at high pressure (from 1 bar up to 20 kbar). The method is based on the extended system formalism introduced by Andersen, where the deterministic equations of motion for the piston degree of freedom are replaced by Langevin equation. We report the structural changes of proteins (ribonuclease A and insulin) and water molecules through radius of gyration, solvent accessible surface area, hydrogen bond pattern, and the topology of water clusters connected by the hydrogen bonded circular network. The solvent accessibility of ribonuclease A is mainly decreased by hydrophilic residues rather than hydrophobic residues under high pressure. From the results of hydrogen bond analysis, we have found that α-helix is more stable than β-sheet under high pressure. In addition, from the analysis of the water cluster, we have observed that for ribonuclease A, 5-membered ring structure is more favorable than 6-membered ring at higher pressure. However, for insulin, the ratio of 5 to 6-ring is constant over the pressure ranges for which we have performed MD simulation. This indicates that the water structure around insulin does not change under high pressure.     

Computational structure characterisation tools in application to ordered and disordered porous materials

In this article, we present a set of computational tools for systematic characterisation of ordered and disordered porous materials. These tools include calculation of the accessible surface area and geometric pore size distribution, analysis of the structure connectivity and percolation analysis of the porous space. We briefly discuss the algorithms behind these calculations. To demonstrate the capabilities of the tools and the type of insights that can be gained from their application, we consider a series of case studies. These case studies include small molecular fragments, several crystalline metal-organic materials, and variants of these materials with induced defects and disorder in their structure. The simulation package is available upon request.     

Hybrid Oligomer of Cyclonucleotides and Deoxynucleotides. A High Anti Left-handed Double Helical DNA Structure

Abstract

It has been shown by us that oligonucleotides containing cyclonucleosides with a high anti glycosidic conformation take left-handed, single and double helical structures (S. Uesugi, J. Yano, E. Yano and M. Ikehara, J. Am. Chem. Soc. 99, 2313 (1977) and references therein). In order to see whether DNA can adopt the high anti left-handed double helical structure or not, a self-complementary hexanucleotide containing 6,2′-O-cyclocytidine (C^○), 8,2′-O-cyclo- guanosine (G^○), deoxycytidine and deoxyguanosine, C^○G^{○dCdGC ○}G^○, was synthesized. Corresponding hexanucleotide containing only cyclonucleosides, C^○G^○C^○G^○C^○G^○, was also synthesized. Their conformation was examined by UV, CD and ¹H NMR spectroscopy. C^○G^○C^○G^○C^○G^○ forms an unusually stable, left-handed duplex. Imino proton NMR spectra and the results of nuclear Overhauser effect experiments strongly suggest that C^○G^○dCdGC^○G^○ take a left-handed double helical structure where the deoxynucleoside residues are involved in hydrogen bonding and take a high anti glycosidic conformation. Thus it is revealed that DNA could form a high anti, left-handed double helix which is different from that of Z-DNA under some constrained conditions.     

Molecular architecture of protein-RNA recognition sites

The molecular architecture of protein-RNA interfaces are analyzed using a non-redundant dataset of 152 protein-RNA complexes. We find that an average protein-RNA interface is smaller than an average protein-DNA interface but larger than an average protein–protein interface. Among the different classes of protein-RNA complexes, interfaces with tRNA are the largest, while the interfaces with the single-stranded RNA are the smallest. Significantly, RNA contributes more to the interface area than its partner protein. Moreover, unlike protein–protein interfaces where the side chain contributes less to the interface area compared to the main chain, the main chain and side chain contributions flipped in protein-RNA interfaces. We find that the protein surface in contact with the RNA in protein-RNA complexes is better packed than that in contact with the DNA in protein-DNA complexes, but loosely packed than that in contact with the protein in protein–protein complexes. Shape complementarity and electrostatic potential are the two major factors that determine the specificity of the protein-RNA interaction. We find that the H-bond density at the protein-RNA interfaces is similar with that of protein-DNA interfaces but higher than the protein–protein interfaces. Unlike protein-DNA interfaces where the deoxyribose has little role in intermolecular H-bonds, due to the presence of an oxygen atom at the 2′ position, the ribose in RNA plays significant role in protein-RNA H-bonds. We find that besides H-bonds, salt bridges and stacking interactions also play significant role in stabilizing protein-nucleic acids interfaces; however, their contribution at the protein–protein interfaces is insignificant.     

Organization of cortical processing for facial movements during licking in cats

We proposed that cortical organization for the execution of adequate licking in cats was processed under the control of two kinds of affiliated groups for face and jaw & tongue movements (Hiraba H, Sato T. 2005A. Cerebral control of face, jaw, and tongue movements in awake cats: Changes in regional cerebral blood flow during lateral feeding Somatosens Mot Res 22:307–317). We assumed the cortical organization for face movements from changes in MRN (mastication-related neuron) activities recorded at area M (motor cortex) and orofacial behaviors after the lesion in the facial SI (facial region in the primary somatosensory cortex). Although we showed the relationship between facial SI (area 3b) and area M (area 4δ), the property of area C (area 3a) was not fully described. The aim of this present study is to investigate the functional role of area C (the anterior part of the coronal sulcus) that transfers somatosensory information in facial SI to area M, as shown in a previous paper (Hiraba H. 2004. The function of sensory information from the first somatosensory cortex for facial movements during ingestion in cats Somatosens Mot Res 21:87--97). We examined the properties of MRNs in area C and changes in orofacial behaviors after the area C or area M lesion. MRNs in area C had in common RFs in the lingual, perioral, and mandibular parts, and activity patterns of MRNs showed both post- and pre-movement types. Furthermore, cats with the area C lesion showed similar disorders to cats with the area M lesion, such as the dropping of food from the contralateral mouth, prolongation of the period of ingestion and mastication, and so on. From these results, we believe firmly the organization of unilateral cortical processing in facial SI, area C, and area M for face movements during licking.     

Thermotactile thresholds at the fingertip: Effect of contact area and contact location

Thresholds for the detection of changes in temperature are used to indicate neuropathy, but a variety of different contact areas and contact locations are used. This study was designed to determine the effects of variations in contact area and contact location on both warm and cool thresholds at the fingertip. With 20 healthy subjects (10 females and 10 males aged 20–30 years), warm thresholds and cool thresholds were determined in two separate sessions using the method of limits. In the first part of each session, thresholds were determined around the centre of the whorl using circular contactors with five different diameters (3, 6, 9, 12, and 55 mm). In the second part of each session, thresholds were determined using two contactors (6- and 12-mm diameter) at three locations along the fingertip: (i) distal (5 mm from the nail), (ii) middle (centre of whorl), and (iii) proximal (3 mm from the distal interphalangeal joint). With increasing contact area, the warm thresholds decreased, the cool thresholds increased, and the inter-subject variability in both warm and cool thresholds decreased. Using the 6-mm diameter contactor, warm thresholds were independent of location but cool thresholds increased from distal to proximal locations. It is concluded that temperature sensitivity at the fingertip increases with increasing area of contact, with the variability in thresholds consistent with the existence of warm and cool “insensitive fields”. The findings show that the influence of contact area and contact location should be considered when assessing thermotactile thresholds at the fingertip.     

Cortical control for mastication in cats: Changes in masticatory movements following lesions in the masticatory cortex

In a previous paper (Hiraba and Sato ) we reported that an accurate mastication might be executed by the cortical processing in bilateral masticatory area (MA)and motor cortices. The aim of this study was to determine if cats with lesion of either unilateral or bilateral MA showed changes in mastication. After exploring mechanoreceptive fields and motor effects of mastication-related neurons (MRNs) in MA using the single unit recording and intracortical microstimulation methods, we made various lesions in MAs with injections of kainic acid (0.1%, 2.0?µl). Since the MA was divided into facial (F) and intraoral (I) projection areas as reported in the previous paper, cats with the unilateral lesion in F or I, and with the bilateral lesion in F & F, I & I or F & I (F on one side and I on other side) were prepared. Cats with unilateral lesion in F or I and with bilateral lesion in F & I showed no changes in mastication except for prolongation of the food intake and masticatory periods. Cats with bilateral lesion into F & F, or I & I showed wider jaw-opening during mastication. Particularly, the latter group showed enormous jaw-opening, delay in the start of mastication and difficulty in manipulating food on the tongue. In all cats with lesions of each type, masticatory and swallowing rhythms remained normal. These findings suggest that accurate mastication is executed by the close integration between F & F and I & I of the bilateral MA.     

The phylogeographical and population genetic approach to the investigation of the genetic diversity patterns in self‐incompatible clonal and polyploid Linnaea borealis subsp. borealis

Surveys of genetic diversity patterns of self‐incompatible clonal polyploid plant species are still scarcer than those of diploid plant species. Therefore, I studied the phylogeographical history of Linnaea borealis subsp. borealis to shed light on the colonization history of this clonal self‐incompatible polyploid plant in Eurasia using selected regions of plastid DNA and genetic diversity patterns of 22 populations of this species employing AFLP markers. I also addressed the question of whether the genetic diversity patterns in L. borealis subsp. borealis in Eurasia are similar to those of earlier published studies of clonal self‐incompatible diploid or polyploid plants. This survey revealed that the shallow phylogeographical history (six plastid haplotypes forming one haplogroup, 100% bootstrap support) and moderate genome‐wide diversity estimated using AFLP markers (Frag_poly = 10.8–38.9%, I = 0.060–0.180, F_ST = 0.289) were general characteristics of L. borealis subsp. borealis in its Eurasian range. The sampling strategy, in most cases at 1–2‐m or even 3–5‐m intervals, showed that a balance between vegetative and sexual reproduction and limited pollen dispersal among compatible mates can be important for genetic diversity patterns in populations of this taxon. Despite the fact that one‐half of the investigated populations were strongly isolated, they still preserved similar levels of genetic diversity across the geographical range. I found no support for the hypothesis that a bottleneck and/or inbreeding had accompanied habitat fragmentation as factors shaping genetic diversity. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2013, 173 , 64–76.     

Smaller stomata require less severe leaf drying to close: A case study in Rosa hydrida

Stomata formed at high relative air humidity (RH) close less as leaf dries; an effect that varies depending on the genotype. We here quantified the contribution of each stomatal response characteristic to the higher water loss of high RH-grown plants, and assessed the relationship between response characteristics and intraspecific variation in stomatal size. Stomatal size (length multiplied by width), density and responsiveness to desiccation, as well as pore dimensions were analyzed in ten rose cultivars grown at moderate (60%) or high (85%) RH. Leaf morphological components and transpiration at growth conditions were also assessed. High growth RH resulted in thinner (11%) leaves with larger area. A strong positive genetic correlation of daytime and nighttime transpiration at either RH was observed. Stomatal size determined pore area (r = 0.7) and varied by a factor of two, as a result of proportional changes in length and width. Size and density of stomata were not related. Following desiccation, high RH resulted in a significantly lower (6–19%) decline of transpiration in three cultivars, whereas the relative water content (RWC) of high RH-expanded leaflets was lower (29–297%) in seven cultivars. The lower RWC of these leaflets was caused by (a) higher (33–72%) stable transpiration and/or (b) lower (12–143%) RWC at which this stable transpiration occurred, depending on the cultivar. Stomatal size was significantly correlated with both characteristics (r = 0.5 and -0.7, respectively). These results indicate that stomatal size explains much of the intraspecific variation in the regulation of transpiration upon water deprivation on rose.