Deciphering complex patterns of class-I HLA-peptide cross-reactivity via hierarchical grouping

T-cell responses in humans are initiated by the binding of a peptide antigen to a human leukocyte antigen (HLA) molecule. The peptide-HLA complex then recruits an appropriate T cell, leading to cell-mediated immunity. More than 2000 HLA class-I alleles are known in humans, and they vary only in their peptide-binding grooves. The polymorphism they exhibit enables them to bind a wide range of peptide antigens from diverse sources. HLA molecules and peptides present a complex molecular recognition pattern, as many peptides bind to a given allele and a given peptide can be recognized by many alleles. A powerful grouping scheme that not only provides an insightful classification, but is also capable of dissecting the physicochemical basis of recognition specificity is necessary to address this complexity. We present a hierarchical classification of 2010 class-I alleles by using a systematic divisive clustering method. All-pair distances of alleles were obtained by comparing binding pockets in the structural models. By varying the similarity thresholds, a multilevel classification was obtained, with 7 supergroups, each further subclassifying to yield 72 groups. An independent clustering performed based only on similarities in their epitope pools correlated highly with pocket-based clustering. Physicochemical feature combinations that best explain the basis of clustering are identified. Mutual information calculated for the set of peptide ligands enables identification of binding site residues contributing to peptide specificity. The grouping of HLA molecules achieved here will be useful for rational vaccine design, understanding disease susceptibilities and predicting risk of organ transplants.

Estimation of Power Input to Complex Dielectric Barrier Discharge Reactor Geometries used in NOx Cleaning

Studies were carried out to estimate the power input to Dielectric Barrier Discharge (DBD) reactors powered by AC high voltage in the context of their application in non-thermal plasma cleaning of exhaust gases. Power input to the reactors was determined both theoretically and experimentally. Four different reactor geometries energized with 50 Hz and 1.5 kHz AC excitation were considered for the study. The theoretically estimated power using Manley's equation was found to agree with the experimental results. Results show that the analytically computed capacitance, without including the electrode edge effects, gives sufficiently good results that are matching with the measured values. For complex geometries where analytical calculation of capacitance is often difficult, a novel method of estimating the reactor capacitance, and hence the power input to the reactor, was introduced in this paper. The predicted results were validated with experiments.

Quick measurement of H-1-F-19 scalar couplings from the complex NMR spectra by real-time spin edition

Most often the measurement of VHF from the conventional 1D H-1 NMR spectrum is severely hindered consequent to similar magnitudes of JHF and JHH couplings and the spectral multiplicity pattern. The present study reports a new 1D NMR technique based on real time spin edition, which removes all JHF and JHH while retaining only VHF of a chosen fluorine. The obtained spectrum is significantly simplified and permits straightforward determination of all possible VHF values of a chosen fluorine. Due to one dimensional nature, the method is much faster compared to 2D GET-SERF by 1-2 orders of magnitude. (C) 2015 Elsevier B.V. All rights reserved.

A Dynamic Comprehensive Method for Landslide Control

A slope failure is developed due to progressive external loads and deteriorations of slope geomaterials, thus forming a progressive and dynamic development and occurrence of landslides. Site geological properties and other active factors such as hydrodynamic load and human activities are complex and usually unknown, thus this dynamic development and occurrence of landslides can only be understood through the progressive accumulation of knowledge on the landslides. For such a progressive process, this paper proposes a dynamic comprehensive control method for landslide control. This control method takes full advantage of updated monitoring data and site investigations of landslides, and emphasizes the implementation of possible measures for landslide control at reasonable stages and in different groups. These measures are to prevent the occurrence of a landslide disaster. As a case study, a landslide project at the Panluo open-pit iron mine is analyzed to illustrate this dynamic comprehensive control method.

Analysis of super compact finite difference method and application to simulation of vortex-shock interaction

Turbulence and aeroacoustic noise high-order accurate schemes are required, and preferred, for solving complex flow fields with multi-scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth-order traditional and compact finite difference approximation. The comparison shows that the sixth-order accurate super compact method has higher resolving efficiency. The sixth-order super compact method, with a three-stage Runge-Kutta method for approximation of the compressible Navier-Stokes equations, is used to solve the complex flow structures induced by vortex-shock interactions. The basic nature of the near-field sound generated by interaction is studied.

A high order accurate difference scheme for complex flow fields

A high order accurate finite difference method for direct numerical simulation of coherent structure in the mixing layers is presented. The reason for oscillation production in numerical solutions is analyzed, It is caused by a nonuniform group velocity of wavepackets. A method of group velocity control for the improvement of the shock resolution is presented. In numerical simulation the fifth-order accurate upwind compact difference relation is used to approximate the derivatives in the convection terms of the compressible N-S equations, a sixth-order accurate symmetric compact difference relation is used to approximate the viscous terms, and a three-stage R-K method is used to advance in time. In order to improve the shock resolution the scheme is reconstructed with the method of diffusion analogy which is used to control the group velocity of wavepackets. (C) 1997 Academic Press.

Complex stress functions and plastic zone sizes for notch cracks subjected to various loading conditions

In this paper, the conformal mapping method is used to solve the plane problem of an infinite plate containing a central lip-shaped notch subjected to biaxial loading at a remote boundary or a surface uniform pressure on the notch. The stress intensity factors KI and KII are obtained by the derived complex stress functions. The simple analytical expressions can be applied to the situation of cracks originating from a circular or an elliptical notch. The plastic zone sizes for such notch cracks are subsequently evaluated in light of the Dugdale strip yield concept. The results are consistent with available numerical data.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

An Integrated Assessment of the Introduction of Lionfish (Pterois volitans/miles complex) to the Western Atlantic Ocean.

Lionfish (Pterois volitans/miles complex) are venomous coral reef fishes from the Indian and western Pacific oceans that are now found in the western Atlantic Ocean. Adult lionfish have been observed from Miami, Florida to Cape Hatteras, North Carolina, and juvenile lionfish have been observed off North Carolina, New York, and Bermuda. The large number of adults observed and the occurrence of juveniles indicate that lionfish are established and reproducing along the southeast United States coast. Introductions of marine species occur in many ways. Ballast water discharge, a very common method of introduction for marine invertebrates, is responsible for many freshwater fish introductions. In contrast, most marine fish introductions result from intentional stocking for fishery purposes. Lionfish, however, likely were introduced via unintentional or intentional aquarium releases, and the introduction of lionfish into United States waters should lead to an assessment of the threat posed by the aquarium trade as a vector for fish introductions. Currently, no management actions are being taken to limit the effect of lionfish on the southeast United States continental shelf ecosystem. Further, only limited funds have been made available for research. Nevertheless, the extent of the introduction has been documented and a forecast of the maximum potential spread of lionfish is being developed. Under a scenario of no management actions and limited research, three predictions are made: ● With no action, the lionfish population will continue to grow along the southeast United States shelf. ● Effects on the marine ecosystem of the southeast United States will become more noticeable as the lionfish population grows. ● There will be incidents of lionfish envenomations of divers and/or fishers along the east coast of the United States. Removing lionfish from the southeast United States continental shelf ecosystem would be expensive and likely impossible. A bounty could be established that would encourage the removal of fish and provide specimens for research. However, the bounty would need to be lower than the price of fish in the aquarium trade (~$25-$50 each) to ensure that captured specimens were from the wild. Such a low bounty may not provide enough incentive for capturing lionfish in the wild. Further, such action would only increase the interaction between the public and lionfish, increasing the risk of lionfish envenomations. As the introduction of lionfish is very likely irreversible, future actions should focus on five areas. 1) The population of lionfish should be tracked. 2) Research should be conducted so that scientists can make better predictions regarding the status of the invasion and the effects on native species, ecosystem function, and ecosystem services. 3) Outreach and education efforts must be increased, both specifically toward lionfish and more generally toward the aquarium trade as a method of fish introductions. 4) Additional regulation should be considered to reduce the frequency of marine fish introduction into U.S. waters. However, the issue is more complicated than simply limiting the import of non-native species, and these complexities need to be considered simultaneously. 5) Health care providers along the east coast of the United States need to be notified that a venomous fish is now resident along the southeast United States. The introduction and spread of lionfish illustrates the difficulty inherent in managing introduced species in marine systems. Introduced species often spread via natural mechanisms after the initial introduction. Efforts to control the introduction of marine fish will fail if managers do not consider the natural dispersal of a species following an introduction. Thus, management strategies limiting marine fish introductions need to be applied over the scale of natural ecological dispersal to be effective, pointing to the need for a regional management approach defined by natural processes not by political boundaries. The introduction and success of lionfish along the east coast should change the long-held perception that marine fish invasions are a minimal threat to marine ecosystems. Research is needed to determine the effects of specific invasive fish species in specific ecosystems. More broadly, a cohesive plan is needed to manage, mitigate and minimize the effects of marine invasive fish species on ecosystems that are already compromised by other human activities. Presently, the magnitude of marine fish introductions as a stressor on marine ecosystems cannot be quantified, but can no longer be dismissed as negligible. (PDF contains 31 pages)

A new method to analyse the Peierls-Nabarro model of an edge dislocation in a rectangular plate

A new method is presented here to analyse the Peierls-Nabarro model of an edge dislocation in a rectangular plate. The analysis is based on the superposition scheme and series expansions of complex potentials. The stress field and dislocation density field on the slip plane can be expressed as the first and the second Chebyshev polynomial series respectively. Two sets of governing equations are obtained on the slip plane and outer boundary of the rectangular plate respectively. Three numerical methods are used to solve the governing equations.

A neural extracellular matrix-based method for in vitro hippocampal neuron culture and dopaminergic differentiation of neural stem cells

Background: The ability to recreate an optimal cellular microenvironment is critical to understand neuronal behavior and functionality in vitro. An organized neural extracellular matrix (nECM) promotes neural cell adhesion, proliferation and differentiation. Here, we expanded previous observations on the ability of nECM to support in vitro neuronal differentiation, with the following goals: (i) to recreate complex neuronal networks of embryonic rat hippocampal cells, and (ii) to achieve improved levels of dopaminergic differentiation of subventricular zone (SVZ) neural progenitor cells. Methods: Hippocampal cells from E18 rat embryos were seeded on PLL- and nECM-coated substrates. Neurosphere cultures were prepared from the SVZ of P4-P7 rat pups, and differentiation of neurospheres assayed on PLL- and nECM-coated substrates. Results: When seeded on nECM-coated substrates, both hippocampal cells and SVZ progenitor cells showed neural expression patterns that were similar to their poly-L-lysine-seeded counterparts. However, nECM-based cultures of both hippocampal neurons and SVZ progenitor cells could be maintained for longer times as compared to poly-L-lysine-based cultures. As a result, nECM-based cultures gave rise to a more branched neurite arborization of hippocampal neurons. Interestingly, the prolonged differentiation time of SVZ progenitor cells in nECM allowed us to obtain a purer population of dopaminergic neurons. Conclusions: We conclude that nECM-based coating is an efficient substrate to culture neural cells at different stages of differentiation. In addition, neural ECM-coated substrates increased neuronal survival and neuronal differentiation efficiency as compared to cationic polymers such as poly-L-lysine.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

A load simulation method of piezoelectric actuator in FEM for smart structures

More and more piezoelectric materials and structures have been used for structure control in aviation and aerospace industry. More efficient and convenient computation method for large complex structure with piezoelectric actuation devices is required. A load simulation method of piezoelectric actuation is presented in this paper. By this method, the freedom degree of finite element simulation is significantly reduced, the difficulty in defining in-plane voltage for multi-layers piezoelectric composite is overcome and the transfer computation between material main direction and the element main direction is simplified. The concept of simulation load is comprehensible and suitable for engineers of structure strength in shape and vibration control, thereby is valuable for promoting the application of piezoelectric material and structures in practical aviation and aerospace fields.

Iso-specific Ziegler-Natta polymerization of α-olefins with a single-component organoyttrium catalyst

The synthesis of a sterically tailored ligand array (M)_2((C_5H_2-2-Si(CH_3)_3-4-C(CH_3)_3)S_2i(CH_3)_2]("M_2Bp") (M = Li, 16; K, 19) is described. Transmetallation of Li_2Bp with YCl_3(THF)_3 affords exclusively the C_2 symmetric product rac-[BpY(µ_2-Cl)_2Li(THF)_2], 20. A X-ray crystal structure of 20 has been determined; triclinic, P1, a= 13.110 (8), b = 17.163 (15), c = 20.623 (14) Å, α = 104.02 (7), β = 99.38 (5), γ = 100.24 (6)° , Z = 4, R = 0.056. Transmetallation of K_2Bp with YCl_3(THF)_3 affords the halide free complex rac-BpYCl, 23. The corresponding rac-BpLaCl, 28, is prepared in an anlogous manner. In all cases the achiral meso isomer is not obtained since only for the racemic isomers are the unfavorable steric interactions between the Si(CH3)_3 groups in the narrow portion of the [Cp-M'-Cp] wedge avoided. Alkylation of 20 or 23 with LiCH(Si(CH_3)_3)_2 affords rac-BpYCH(Si(CH_3)_3)_2, 26 in good yield. Alkylation of 28 with LiCH(Si(CH_3)_3)_2 affords rac-BpLaCH(Si(CH_3)_3)_2 29. Hydrogenation of 26 cleanly affords the bridging hydride species [BpY(µ_2-H)]_2, 27, as the homochiral (R,R) and (S,S) dimeric pairs. 26 is an efficient initiator for the polymerization of ethylene to high molecular weight linear polyethylene. 27 catalyzes the polymerization of propylene (25% v/v in methylcyclohexane) and neat samples of 1-butene, 1-pentene, 1-hexene to moderately high molecular weight polymers: polypropylene (M_n = 4,200, PDI 2.32, T_m 157 °C); poly-1-butene (M_n = 8,500, PDI 3.44, T_m 105 °C); poly-1-pentene (M_n = 20,000, PDI 1.99, T_m 73 °C); poly-1-hexene (M_n = 24,000, PDI 1.75, T_m < 25 °C). ^(13)C NMR spectra at the pentad analysis level indicates that the degree of isotacticity is 99% mmmm for all polymer samples. 27 is the first single component iso-specific α-olefin polymerization catalyst. The presumed origins of the high isospecificity are presented.