Structural and functional characteristics of homing endonucleases

Mobile genetic elements constitute a remarkably diverse group of nonessential selfish genes that provide no apparent function to the host. These selfish genes have been implicated in host extinction, speciation and architecture of genetic systems. Homing endonucleases, encoded by the open reading frames embedded in introns or inteins of mobile genetic elements, possess double-stranded DNA-specific endonuclease activity. They inflict sequence-specific double-strand breaks at or near the homing site in intron- or intein-less allele. Subsequently, through nonreciprocal exchange the insertion sequence (intron or intein) is transferred from an intein- or intron-containing allele to an intein- or intron-less allele. The components of host double-strand break repair pathway are thought to finish the "homing" process. Several lines of evidence suggest that homing endonucleases are capable of promoting transposition into ectopic sites within or across genomes for their survival as well as dispersal in natural populations. The occurrence of inteins at high frequencies serves as instructive models for understanding the mechanistic aspects of the process of homing and its evolution. This review focuses on genetic, biochemical, structural, and phylogenetic aspects of homing endonucleases, and their comparison with restriction endonucleases.

A multistep transversal linearization (MTL) method in non-linear dynamics through a Magnus characterization

A new form of a multi-step transversal linearization (MTL) method is developed and numerically explored in this study for a numeric-analytical integration of non-linear dynamical systems under deterministic excitations. As with other transversal linearization methods, the present version also requires that the linearized solution manifold transversally intersects the non-linear solution manifold at a chosen set of points or cross-section in the state space. However, a major point of departure of the present method is that it has the flexibility of treating non-linear damping and stiffness terms of the original system as damping and stiffness terms in the transversally linearized system, even though these linearized terms become explicit functions of time. From this perspective, the present development is closely related to the popular practice of tangent-space linearization adopted in finite element (FE) based solutions of non-linear problems in structural dynamics. The only difference is that the MTL method would require construction of transversal system matrices in lieu of the tangent system matrices needed within an FE framework. The resulting time-varying linearized system matrix is then treated as a Lie element using Magnus’ characterization [W. Magnus, On the exponential solution of differential equations for a linear operator, Commun. Pure Appl. Math., VII (1954) 649–673] and the associated fundamental solution matrix (FSM) is obtained through repeated Lie-bracket operations (or nested commutators). An advantage of this approach is that the underlying exponential transformation could preserve certain intrinsic structural properties of the solution of the non-linear problem. Yet another advantage of the transversal linearization lies in the non-unique representation of the linearized vector field – an aspect that has been specifically exploited in this study to enhance the spectral stability of the proposed family of methods and thus contain the temporal propagation of local errors. A simple analysis of the formal orders of accuracy is provided within a finite dimensional framework. Only a limited numerical exploration of the method is presently provided for a couple of popularly known non-linear oscillators, viz. a hardening Duffing oscillator, which has a non-linear stiffness term, and the van der Pol oscillator, which is self-excited and has a non-linear damping term.

A solution framework for scheduling a BPM with non-identical job dimensions

A Batch Processing Machine (BPM) is one which processes a number of jobs simultaneously as a batch with common beginning and ending times. Also, a BPM, once started cannot be interrupted in between (Pre-emption not allowed). This research is motivated by a BPM in steel casting industry. There are three main stages in any steel casting industry viz., pre-casting stage, casting stage and post-casting stage. A quick overview of the entire process, is shown in Figure 1. There are two BPMs : (1) Melting furnace in the pre-casting stage and (2) Heat Treatment Furnace (HTF) in the post casting stage of steel casting manufacturing process. This study focuses on scheduling the latter, namely HTF. Heat-treatment operation is one of the most important stages of steel casting industries. It determines the final properties that enable components to perform under demanding service conditions such as large mechanical load, high temperature and anti-corrosive processing. In general, different types of castings have to undergo more than one type of heat-treatment operations, where the total heat-treatment processing times change. To have a better control, castings are primarily classified into a number of job-families based on the alloy type such as low-alloy castings and high alloy castings. For technical reasons such as type of alloy, temperature level and the expected combination of heat-treatment operations, the castings from different families can not be processed together in the same batch.

Genetic analysis of structural brain connectivity using DICCCOL models of diffusion MRI in 522 twins

Genetic and environmental factors affect white matter connectivity in the normal brain, and they also influence diseases in which brain connectivity is altered. Little is known about genetic influences on brain connectivity, despite wide variations in the brain's neural pathways. Here we applied the 'DICCCOL' framework to analyze structural connectivity, in 261 twin pairs (522 participants, mean age: 21.8 y ± 2.7SD). We encoded connectivity patterns by projecting the white matter (WM) bundles of all 'DICCCOLs' as a tracemap (TM). Next we fitted an A/C/E structural equation model to estimate additive genetic (A), common environmental (C), and unique environmental/error (E) components of the observed variations in brain connectivity. We found 44 'heritable DICCCOLs' whose connectivity was genetically influenced (α2>1%); half of them showed significant heritability (α2>20%). Our analysis of genetic influences on WM structural connectivity suggests high heritability for some WM projection patterns, yielding new targets for genome-wide association studies.

Ion pairs in non-redundant protein structures

Ion pairs contribute to several functions including the activity of catalytic triads, fusion of viral membranes, stability in thermophilic proteins and solvent-protein interactions. Furthermore, they have the ability to affect the stability of protein structures and are also a part of the forces that act to hold monomers together. This paper deals with the possible ion pair combinations and networks in 25% and 90% non-redundant protein chains. Different types of ion pairs present in various secondary structural elements are analysed. The ion pairs existing between different subunits of multisubunit protein structures are also computed and the results of various analyses are presented in detail. The protein structures used in the analysis are solved using X-ray crystallography, whose resolution is better than or equal to 1.5 angstrom and R-factor better than or equal to 20%. This study can, therefore, be useful for analyses of many protein functions. It also provides insights into the better understanding of the architecture of protein structure.

Using a structural approach to identify relationships between soil and erosion in a semi-humid forested area, South India

Biogeochemical and hydrological cycles are currently studied on a small experimental forested watershed (4.5 km(2)) in the semi-humid South India. This paper presents one of the first data referring to the distribution and dynamics of a widespread red soil (Ferralsols and Chromic Luvisols) and black soil (Vertisols and Vertic intergrades) cover, and its possible relationship with the recent development of the erosion process. The soil map was established from the observation of isolated soil profiles and toposequences, and surveys of soil electromagnetic conductivity (EM31, Geonics Ltd), lithology and vegetation. The distribution of the different parts of the soil cover in relation to each other was used to establish the dynamics and chronological order of formation. Results indicate that both topography and lithology (gneiss and amphibolite) have influenced the distribution of the soils. At the downslope, the following parts of the soil covers were distinguished: i) red soil system, ii) black soil system, iii) bleached horizon at the top of the black soil and iv) bleached sandy saprolite at the base of the black soil. The red soil is currently transforming into black soil and the transformation front is moving upslope. In the bottom part of the slope, the chronology appears to be the following: black soil > bleached horizon at the top of the black soil > streambed > bleached horizon below the black soil. It appears that the development of the drainage network is a recent process, which was guided by the presence of thin black soil with a vertic horizon less than 2 in deep. Three distinctive types of erosional landforms have been identified: 1. rotational slips (Type 1); 2. a seepage erosion (Type 2) at the top of the black soil profile; 3. A combination of earthflow and sliding in the non-cohesive saprolite of the gneiss occurs at midslope (Type 3). Types 1 and 2 erosion are mainly occurring downslope and are always located at the intersection between the streambed and the red soil-black soil contact. Neutron probe monitoring, along an area vulnerable to erosion types 1 and 2, indicates that rotational slips are caused by a temporary watertable at the base of the black soil and within the sandy bleached saprolite, which behaves as a plane of weakness. The watertable is induced by the ephemeral watercourse. Erosion type 2 is caused by seepage of a perched watertable, which occurs after swelling and closing of the cracks of the vertic clay horizon and within a light textured and bleached horizon at the top of black soil. Type 3 erosion is not related to the red soil-black soil system but is caused by the seasonal seepage of saturated throughflow in the sandy saprolite of the gneiss occurring at midslope. (c) 2006 Elsevier B.V. All rights reserved.

Structural and electrical transport of carbon nanofibers derived from dihydro-2,5-furandione

Carbon nanofibers of 50–500 nm diameter and several micrometer length were synthesized by high-temperature pyrolysis of dihydro-2,5-furandione (C4H4O3) in the temperature range of 600–980 °C. The formation of both graphitic and non-graphitic structured carbon fibers was observed in high-resolution transmission electron microscope. The Raman spectra of the samples showed the presence of both the D and G bands of varying intensity and sharpness. The low-temperature electrical transport studies on the samples have shown interesting metal–insulator transitions. The films showed variable range hopping conduction in the insulating regime and power law behavior in the critical regime at low temperatures.

Structural studies of decaying fluid turbulence: Effect of initial conditions

We present results from a systematic numerical study of structural properties of an unforced, incompressible, homogeneous, and isotropic three-dimensional turbulent fluid with an initial energy spectrum that develops a cascade of kinetic energy to large wave numbers. The results are compared with those from a recently studied set of power-law initial energy spectra [C. Kalelkar and R. Pandit, Phys. Rev. E 69, 046304 (2004)] which do not exhibit such a cascade. Differences are exhibited in plots of vorticity isosurfaces, the temporal evolution of the kinetic energy-dissipation rate, and the rates of production of the mean enstrophy along the principal axes of the strain-rate tensor. A crossover between "non-cascade-type" and "cascade-type" behavior is shown numerically for a specific set of initial energy spectra.

Non-Protein Amino Acids in the Design of Secondary Structure Scaffolds

The use of stereochemically constrained amino acids permits the design of short peptides as models for protein secondary structures. Amino acid residues that are restrained to a limited range of backbone torsion angles (ϕ-ψ) may be used as folding nuclei in the design of helices and β-hairpins. α-Amino-isobutyric acid (Aib) and related Cαα dialkylated residues are strong promoters of helix formation, as exemplified by a large body of experimentally determined structures of helical peptides. DPro-Xxx sequences strongly favor type II’ turn conformations, which serve to nucleate registered β-hairpin formation. Appropriately positioned DPro-Xxx segments may be used to nucleate the formation of multistranded antiparallel β-sheet structures. Mixed (α/β) secondary structures can be generated by linking rigid modules of helices and β-hairpins. The approach of using stereochemically constrained residues promotes folding by limiting the local structural space at specific residues. Several aspects of secondary structure design are outlined in this chapter, along with commonly used methods of spectroscopic characterization.

Solvolytic disproportionation of non-stoichiometric oxides of praseodymium and terbiumstar, open

The solvolytic disproportionation of non-stoichiometric PrOχ and TbOχ in acid solutions to produce higher oxides has been investigated. Some new non-stoichiometric phases have been reported. A number of interesting features of the non-stoichiometric rare earth oxides have been discussed and the need for a satisfactory structural model has been pointed out.

Self-regularized pseudo time-marching schemes for structural system identification with static measurements

We explore the application of pseudo time marching schemes, involving either deterministic integration or stochastic filtering, to solve the inverse problem of parameter identification of large dimensional structural systems from partial and noisy measurements of strictly static response. Solutions of such non-linear inverse problems could provide useful local stiffness variations and do not have to confront modeling uncertainties in damping, an important, yet inadequately understood, aspect in dynamic system identification problems. The usual method of least-square solution is through a regularized Gauss-Newton method (GNM) whose results are known to be sensitively dependent on the regularization parameter and data noise intensity. Finite time,recursive integration of the pseudo-dynamical GNM (PD-GNM) update equation addresses the major numerical difficulty associated with the near-zero singular values of the linearized operator and gives results that are not sensitive to the time step of integration. Therefore, we also propose a pseudo-dynamic stochastic filtering approach for the same problem using a parsimonious representation of states and specifically solve the linearized filtering equations through a pseudo-dynamic ensemble Kalman filter (PD-EnKF). For multiple sets of measurements involving various load cases, we expedite the speed of thePD-EnKF by proposing an inner iteration within every time step. Results using the pseudo-dynamic strategy obtained through PD-EnKF and recursive integration are compared with those from the conventional GNM, which prove that the PD-EnKF is the best performer showing little sensitivity to process noise covariance and yielding reconstructions with less artifacts even when the ensemble size is small.

An automated framework for understanding structural variations in the binding grooves of MHC class II molecules

Background: MHC/HLA class II molecules are important components of the immune system and play a critical role in processes such as phagocytosis. Understanding peptide recognition properties of the hundreds of MHC class II alleles is essential to appreciate determinants of antigenicity and ultimately to predict epitopes. While there are several methods for epitope prediction, each differing in their success rates, there are no reports so far in the literature to systematically characterize the binding sites at the structural level and infer recognition profiles from them. Results: Here we report a new approach to compare the binding sites of MHC class II molecules using their three dimensional structures. We use a specifically tuned version of our recent algorithm, PocketMatch. We show that our methodology is useful for classification of MHC class II molecules based on similarities or differences among their binding sites. A new module has been used to define binding sites in MHC molecules. Comparison of binding sites of 103 MHC molecules, both at the whole groove and individual sub-pocket levels has been carried out, and their clustering patterns analyzed. While clusters largely agree with serotypic classification, deviations from it and several new insights are obtained from our study. We also present how differences in sub-pockets of molecules associated with a pair of autoimmune diseases, narcolepsy and rheumatoid arthritis, were captured by PocketMatch(13). Conclusion: The systematic framework for understanding structuralvariations in MHC class II molecules enables large scale comparison of binding grooves and sub-pockets, which is likely to have direct implications towards predicting epitopes and understanding peptide binding preferences.

Self-regularized pseudo time-marching schemes for structural system identification with static measurements

We explore the application of pseudo time marching schemes, involving either deterministic integration or stochastic filtering, to solve the inverse problem of parameter identification of large dimensional structural systems from partial and noisy measurements of strictly static response. Solutions of such non-linear inverse problems could provide useful local stiffness variations and do not have to confront modeling uncertainties in damping, an important, yet inadequately understood, aspect in dynamic system identification problems. The usual method of least-square solution is through a regularized Gauss-Newton method (GNM) whose results are known to be sensitively dependent on the regularization parameter and data noise intensity. Finite time, recursive integration of the pseudo-dynamical GNM (PD-GNM) update equation addresses the major numerical difficulty associated with the near-zero singular values of the linearized operator and gives results that are not sensitive to the time step of integration. Therefore, we also propose a pseudo-dynamic stochastic filtering approach for the same problem using a parsimonious representation of states and specifically solve the linearized filtering equations through apseudo-dynamic ensemble Kalman filter (PD-EnKF). For multiple sets ofmeasurements involving various load cases, we expedite the speed of the PD-EnKF by proposing an inner iteration within every time step. Results using the pseudo-dynamic strategy obtained through PD-EnKF and recursive integration are compared with those from the conventional GNM, which prove that the PD-EnKF is the best performer showing little sensitivity to process noise covariance and yielding reconstructions with less artifacts even when the ensemble size is small. Copyright (C) 2009 John Wiley & Sons, Ltd.

An automated framework for understanding structural variations in the binding grooves of MHC class II molecules

Background: MHC/HLA class II molecules are important components of the immune system and play a critical role in processes such as phagocytosis. Understanding peptide recognition properties of the hundreds of MHC class II alleles is essential to appreciate determinants of antigenicity and ultimately to predict epitopes. While there are several methods for epitope prediction, each differing in their success rates, there are no reports so far in the literature to systematically characterize the binding sites at the structural level and infer recognition profiles from them. Results: Here we report a new approach to compare the binding sites of MHC class II molecules using their three dimensional structures. We use a specifically tuned version of our recent algorithm, PocketMatch. We show that our methodology is useful for classification of MHC class II molecules based on similarities or differences among their binding sites. A new module has been used to define binding sites in MHC molecules. Comparison of binding sites of 103 MHC molecules, both at the whole groove and individual sub-pocket levels has been carried out, and their clustering patterns analyzed. While clusters largely agree with serotypic classification, deviations from it and several new insights are obtained from our study. We also present how differences in sub-pockets of molecules associated with a pair of autoimmune diseases, narcolepsy and rheumatoid arthritis, were captured by PocketMatch(13). Conclusion: The systematic framework for understanding structural variations in MHC class II molecules enables large scale comparison of binding grooves and sub-pockets, which is likely to have direct implications towards predicting epitopes and understanding peptide binding preferences.

USING FLUORESCENT AND NON-FLUORESCENT STEROLS TO STUDY RAFTS AND INTRACELLULAR CHOLESTEROL TRANSPORT IN MAMMALIAN CELLS

Cholesterol is an essential component in the membranes of most eukaryotic cells, in which it mediates many functions including membrane fluidity, permeability and the formation of ordered membrane domains. In this work a fluorescent and a non-fluorescent cholesterol analog were characterized as tools to study cholesterol. Next, these analogs were used to study two specific cell biological processes that involve cholesterol, i.e. the structure and function of ordered membrane domains/rafts and intracellular cholesterol transport. The most common method for studying ordered membrane domains is by disrupting them by cholesterol depletion. Because cholesterol depletion affects many cellular functions besides those mediated by membrane domains, this procedure is highly unspecific. The cellular exchange of cholesterol by desmosterol as a tool to study ordered membrane domains was characterized. It turned out that the ability of desmosterol to form and stabilize membrane domains in vitro was weaker compared to cholesterol. This result was reinforced by atomistic scale simulations that indicated that desmosterol has a lower ordering effect on phospholipid acyl chains. Three procedures were established for exchanging cellular cholesterol by desmosterol. In cells in which desmosterol was the main sterol, insulin signaling was attenuated. The results suggest that this was caused by desmosterol destabilizing membrane rafts. Contrary to its effect on ordered membrane domains it was found that replacing cholesterol by desmosterol does not change cell growth/viability, subcellular sterol distribution, Golgi integrity, secretory pathway, phospholipid composition and membrane fluidity. Together these results suggest that exchanging cellular cholesterol by desmosterol provides a selective tool for perturbing rafts. Next, the importance of cholesterol for the structure and function of caveolae was analyzed by exchanging the cellular cholesterol by desmosterol. The sterol exchange reduced the stability of caveolae as determined by detergent resistance of caveolin-1 and heat resistance of caveolin-1 oligomers. Also the sterol exchange led to aberrations in the caveolar structure; the morphology of caveolae was altered and there was a larger variation in the amount of caveolin-1 molecules per caveola. These results demonstrate that cholesterol is important for caveolar stability and structural homogeneity. In the second part of this work a fluorescent cholesterol analog was characterized as a tool to study cholesterol transport. Tight control of the intracellular cholesterol distribution is essential for many cellular processes. An important mechanism by which cells regulate their membrane cholesterol content is by cholesterol traffic, mostly from the plasma membrane to lipid droplets. The fluorescent sterol probe BODIPY-cholesterol was characterized as a tool to analyze cholesterol transport between the plasma membrane, the endoplasmic reticulum (ER) and lipid droplets. The behavior of BODIPY-cholesterol was compared to that of natural sterols, using both biochemical and live-cell microcopy assays. The results show that the transport kinetics of BODIPY-cholesterol between the plasma membrane, the ER and lipid droplets is similar to that of unesterified cholesterol. Next, BODIPY-cholesterol was utilized to analyze the importance of oxysterol binding protein related proteins (ORPs) for cholesterol transport between the plasma membrane, the ER, and lipid droplets in mammalian cells. By overexpressing all human ORPs it turned out that especially ORP1S and ORP2 enhanced sterol transport from the plasma membrane to lipid droplets. Our results suggest that the increased sterol transport takes place between the plasma membrane and ER and not between the ER and lipid droplets. Simultaneous knockdown of ORP1S and ORP2 resulted in a moderate but significant inhibition of sterol traffic from the plasma membrane to ER and lipid droplets, suggesting a physiological role for these ORPs in this process. The two phenylalanines in an acidic tract (FFAT) motif in ORPs, which mediates interaction with vesicle associated membrane protein associated proteins (VAPs) in the ER, was not necessary for mediating sterol transport. However, VAP silencing slowed down sterol transport, most likely by destabilizing ORPs containing a FFAT motif.