A functionally defined model for the M2 proton channel of influenza A virus suggests a mechanism for its ion selectivity

The M2 protein from influenza A virus forms proton-selective channels that are essential to viral function and are the target of the drug amantadine. Cys scanning was used to generate a series of mutants with successive substitutions in the transmembrane segment of the protein, and the mutants were expressed in Xenopus laevis oocytes. The effect of the mutations on reversal potential, ion currents, and amantadine resistance were measured. Fourier analysis revealed a periodicity consistent with a four-stranded coiled coil or helical bundle. A three-dimensional model of this structure suggests a possible mechanism for the proton selectivity of the M2 channel of influenza virus.

On the well-posedness of a mathematical model for lithium-ion batteries.

We discuss the well-posedness of a mathematical model that is used in the literature for the simulation of lithium-ion batteries. First, a mathematical model based on a macrohomogeneous approach is presented, following previous work. Then it is shown, from a physical and a mathematical point of view, that a boundary condition widely used in the literature is not correct. Although the errors could be just sign typos (which can be explained as carelessness in the use of d/dx versus d/dn, with n the outward unit vector) and authors using this model probably use the correct boundary condition when they solve it in order to do simulations, readers should be aware of the right choice. Therefore, the deduction of the correct boundary condition is done here, and a mathematical study of the well-posedness of the corresponding problem is presented.

Finite element model of a small automobile impacting a rigid pole. Final report.

Federal Highway Administration, Office of Safety and Traffic Operations Research Development, McLean, Va.

Crash III model improvements: MRB-to-car side impact of a 90 [degree] moving rigid barrier to a 1983 Mitsubishi Tredia: test #1 - 16.7 mph, test #2 - 33.5 mph. Test report.

Mode of access: Internet.

CRASH III model improvements: MRB-to-car side impact test of a 90 [degree] moving rigid barrier to a 1981 Chevrolet Citation test speed 35.2 mph.

Mode of access: Internet.

A novel concentration dependent amino acid ion pair strategy to mediate drug permeation using indomethacin as a model insoluble drug

Assessment of oral drug bioavailability is an important parameter for new chemical entities (NCEs) in drug development cycle. After evaluating the pharmacological response of these new molecules, the following critical stage is to investigate their in vitro permeability. Despite the great success achieved by prodrugs, covalent linking the drug molecule with a hydrophobic moiety might result in a new entity that might be toxic or ineffective. Therefore, an alternative that would improve the drug uptake without affecting the efficacy of the drug molecule would be advantageous. The aim of the current study is to investigate the effect of ion-pairing on the permeability profile of a model drug: indomethacin (IND) to understand the mechanism behind the permeability improvement across Caco-2 monolayers. Arginine and lysine formed ion-pairs with IND at various molar ratios 1:1, 1:2, 1:4 and 1:8 as reflected by the double reciprocal graphs. The partitioning capacities of the IND were evaluated using octanol/water partitioning studies and the apparent permeabilities (P app) were measured across Caco-2 monolayers for the different formulations. Partitioning studies reflected the high hydrophobicity of IND (Log P = 3) which dropped upon increasing the concentrations of arginine/lysine in the ion pairs. Nevertheless, the prepared ion pairs improved IND permeability especially after 60 min of the start of the experiment. Coupling partitioning and permeability results suggest a decrease in the passive transcellular uptake due to the drop in IND portioning capacities and a possible involvement of active carriers. Future work will investigate which transport gene might be involved in the absorption of the ion paired formulations using molecular biology technologies. © 2014 Elsevier B.V. All rights reserved.

In situ transmission electron microscopy ion irradiations of model iron-chromium alloys

Iron-chromium alloys are used as a model to study the microstructural evolution of defects in irradiated structural steel components of a nuclear reactor. We examine the effects of temperature and chromium concentration on the defect evolution and segregation behavior in the early stages of damage. In situ irradiations are conducted in a transmission electron microscope (TEM) at 300°C and 450°C with 150keV iron ions in single crystal Fe14Cr and Fe19Cr bicrystal to doses of 2E15 ions/cm^2. The microstructures resulting from annealing and irradiation of the alloy are characterized by analysis of TEM micrographs and diffraction patterns and compared with those of irradiated pure iron. We found the irradiation temperature to have little effect on the microstructural development. We also found that the presence of chromium in the sample leads to defect populations with small average loop size and no extended or nested loop structures, in contrast to the populations of large extended loops seen in irradiated pure iron. A very weak dependence was found on the specific chromium content of the alloy. Chromium was shown to suppress defect growth by inhibiting defect mobility in the alloy. While defects in pure iron are highly mobile and able to grow, those in the FeCr alloys remained small and relatively motionless due to the pinning effect of the chromium.

Creation of a validated 3D finite element model of an ovine tibia

Introduction Ovine models are widely used in orthopaedic research. To better understand the impact of orthopaedic procedures computer simulations are necessary. 3D finite element (FE) models of bones allow implant designs to be investigated mechanically, thereby reducing mechanical testing. Hypothesis We present the development and validation of an ovine tibia FE model for use in the analysis of tibia fracture fixation plates. Material & Methods Mechanical testing of the tibia consisted of an offset 3-pt bend test with three repetitions of loading to 350N and return to 50N. Tri-axial stacked strain gauges were applied to the anterior and posterior surfaces of the bone and two rigid bodies – consisting of eight infrared active markers, were attached to the ends of the tibia. Positional measurements were taken with a FARO arm 3D digitiser. The FE model was constructed with both geometry and material properties derived from CT images of the bone. The elasticity-density relationship used for material property determination was validated separately using mechanical testing. This model was then transformed to the same coordinate system as the in vitro mechanical test and loads applied. Results Comparison between the mechanical testing and the FE model showed good correlation in surface strains (difference: anterior 2.3%, posterior 3.2%). Discussion & Conclusion This method of model creation provides a simple method for generating subject specific FE models from CT scans. The use of the CT data set for both the geometry and the material properties ensures a more accurate representation of the specific bone. This is reflected in the similarity of the surface strain results.

CYR61 (CCN1) protein expression during fracture healing in an ovine tibial model and its relation to the mechanical fixation stability

The formation of new blood vessels is a prerequisite for bone healing. CYR61 (CCN1), an extracellular matrix-associated signaling protein, is a potent stimulator of angiogenesis and mesenchymal stem cell expansion and differentiation. A recent study showed that CYR61 is expressed during fracture healing and suggested that CYR61 plays a significant role in cartilage and bone formation. The hypothesis of the present study was that decreased fixation stability, which leads to a delay in healing, would lead to reduced CYR61 protein expression in fracture callus. The aim of the study was to quantitatively analyze CYR61 protein expression, vascularization, and tissue differentiation in the osteotomy gap and relate to the mechanical fixation stability during the course of healing. A mid-shaft osteotomy of the tibia was performed in two groups of sheep and stabilized with either a rigid or semirigid external fixator, each allowing different amounts of interfragmentary movement. The sheep were sacrificed at 2, 3, 6, and 9 weeks postoperatively. The tibiae were tested biomechanically and histological sections from the callus were analyzed immunohistochemically with regard to CYR61 protein expression and vascularization. Expression of CYR61 protein was upregulated at the early phase of fracture healing (2 weeks), decreasing over the healing time. Decreased fixation stability was associated with a reduced upregulation of the CYR61 protein expression and a reduced vascularization at 2 weeks leading to a slower healing. The maximum cartilage callus fraction in both groups was reached at 3 weeks. However, the semirigid fixator group showed a significantly lower CYR61 immunoreactivity in cartilage than the rigid fixator group at this time point. The fraction of cartilage in the semirigid fixator group was not replaced by bone as quickly as in the rigid fixator group leading to an inferior histological and mechanical callus quality at 6 weeks and therefore to a slower healing. The results supply further evidence that CYR61 may serve as an important regulator of bone healing.

Stochastic models and simulation of ion channel dynamics

The behaviour of ion channels within cardiac and neuronal cells is intrinsically stochastic in nature. When the number of channels is small this stochastic noise is large and can have an impact on the dynamics of the system which is potentially an issue when modelling small neurons and drug block in cardiac cells. While exact methods correctly capture the stochastic dynamics of a system they are computationally expensive, restricting their inclusion into tissue level models and so approximations to exact methods are often used instead. The other issue in modelling ion channel dynamics is that the transition rates are voltage dependent, adding a level of complexity as the channel dynamics are coupled to the membrane potential. By assuming that such transition rates are constant over each time step, it is possible to derive a stochastic differential equation (SDE), in the same manner as for biochemical reaction networks, that describes the stochastic dynamics of ion channels. While such a model is more computationally efficient than exact methods we show that there are analytical problems with the resulting SDE as well as issues in using current numerical schemes to solve such an equation. We therefore make two contributions: develop a different model to describe the stochastic ion channel dynamics that analytically behaves in the correct manner and also discuss numerical methods that preserve the analytical properties of the model.

A subject-specific model of human buttocks and thighs in a seated posture

Finite element analyses of the human body in seated postures requires digital models capable of providing accurate and precise prediction of the tissue-level response of the body in the seated posture. To achieve such models, the human anatomy must be represented with high fidelity. This information can readily be defined using medical imaging techniques such as Magnetic Resonance Imaging (MRI) or Computed Tomography (CT). Current practices for constructing digital human models, based on the magnetic resonance (MR) images, in a lying down (supine) posture have reduced the error in the geometric representation of human anatomy relative to reconstructions based on data from cadaveric studies. Nonetheless, the significant differences between seated and supine postures in segment orientation, soft-tissue deformation and soft tissue strain create a need for data obtained in postures more similar to the application posture. In this study, we present a novel method for creating digital human models based on seated MR data. An adult-male volunteer was scanned in a simulated driving posture using a FONAR 0.6T upright MRI scanner with a T1 scanning protocol. To compensate for unavoidable image distortion near the edges of the study, images of the same anatomical structures were obtained in transverse and sagittal planes. Combinations of transverse and sagittal images were used to reconstruct the major anatomical features from the buttocks through the knees, including bone, muscle and fat tissue perimeters, using Solidworks® software. For each MR image, B-splines were created as contours for the anatomical structures of interest, and LOFT commands were used to interpolate between the generated Bsplines. The reconstruction of the pelvis, from MR data, was enhanced by the use of a template model generated in previous work CT images. A non-rigid registration algorithm was used to fit the pelvis template into the MR data. Additionally, MR image processing was conducted to both the left and the right sides of the model due to the intended asymmetric posture of the volunteer during the MR measurements. The presented subject-specific, three-dimensional model of the buttocks and thighs will add value to optimisation cycles in automotive seat development when used in simulating human interaction with automotive seats.

Modeling ion channel dynamics through reflected stochastic differential equations

Ion channels are membrane proteins that open and close at random and play a vital role in the electrical dynamics of excitable cells. The stochastic nature of the conformational changes these proteins undergo can be significant, however current stochastic modeling methodologies limit the ability to study such systems. Discrete-state Markov chain models are seen as the "gold standard," but are computationally intensive, restricting investigation of stochastic effects to the single-cell level. Continuous stochastic methods that use stochastic differential equations (SDEs) to model the system are more efficient but can lead to simulations that have no biological meaning. In this paper we show that modeling the behavior of ion channel dynamics by a reflected SDE ensures biologically realistic simulations, and we argue that this model follows from the continuous approximation of the discrete-state Markov chain model. Open channel and action potential statistics from simulations of ion channel dynamics using the reflected SDE are compared with those of a discrete-state Markov chain method. Results show that the reflected SDE simulations are in good agreement with the discrete-state approach. The reflected SDE model therefore provides a computationally efficient method to simulate ion channel dynamics while preserving the distributional properties of the discrete-state Markov chain model and also ensuring biologically realistic solutions. This framework could easily be extended to other biochemical reaction networks. © 2012 American Physical Society.

A boundary preserving numerical algorithm for the Wright-Fisher model with mutation

The Wright-Fisher model is an Itô stochastic differential equation that was originally introduced to model genetic drift within finite populations and has recently been used as an approximation to ion channel dynamics within cardiac and neuronal cells. While analytic solutions to this equation remain within the interval [0,1], current numerical methods are unable to preserve such boundaries in the approximation. We present a new numerical method that guarantees approximations to a form of Wright-Fisher model, which includes mutation, remain within [0,1] for all time with probability one. Strong convergence of the method is proved and numerical experiments suggest that this new scheme converges with strong order 1/2. Extending this method to a multidimensional case, numerical tests suggest that the algorithm still converges strongly with order 1/2. Finally, numerical solutions obtained using this new method are compared to those obtained using the Euler-Maruyama method where the Wiener increment is resampled to ensure solutions remain within [0,1].