Synthesis of an erodible biomimetic hydrogel for drug delivery using native chemical ligation

Hydrogels are hydrophilic, three dimensional polymers that imbibe large quantities of water while remaining insoluble in aqueous solutions due to chemical or physical cross-linking. The polymers swell in water or biological fluids, immobilizing the bioactive agent, leading to drug release in a well-defined specific manner. Thus the hydrogels’ elastic properties, swellability and biocompatibility make them excellent formulations for drug delivery. Currently, many drug potencies and therapeutic effects are limited or otherwise reduced because of the partial degradation that occurs before the administered drug reaches the desired site of action. On the other hand, sustained release medications release drugs continually, rather than providing relief of symptoms and protection solely when necessary. In fact, it would be much better if drugs could be administered in a manner that precisely matches physiological needs at desired times and at the desired site (site specific targeting). There is therefore an unmet need to develop controlled drug delivery systems especially for delivery of peptide and protein bound drugs. The purpose of this project is to produce hydrogels for structural drug delivery and time-dependent sustained release of drugs (bioactive agents). We use an innovative polymerisation strategy based on native chemical ligation (NCL) to covalently cross-link polymers to form hydrogels. When mixed in aqueous solution, four armed (polyethylene glycol) amine (PEG-4A) end functionalised with thioester and four branched Nterminal cysteine peptide dendrimers spontaneously conjugated to produce biomimetic hydrogels. These hydrogels showed superior resistance to shear stress compared to an equivalent PEG macromonomer system and were shown to be proteolytically degradable with concomitant release of a model payload molecule. This is the first report of a peptide dendrimers/PEG macromonomer approach to hydrogel production and opens up the prospect of facile hydrogel synthesis together with tailored payload release.

Preliminary bond-slip model for CFRP sheets bonded to steel plates

Carbon fiber reinforced polymer (CFRP) sheets have established a strong position as an effective method for innovative structural rehabilitation. However, the use of externally bonded CFRP in the repair and rehabilitation of steel structures is a relatively new technique that has the potential to improve the way structures are repaired. An important step toward understanding bond behaviour is to have an estimation of local bond stress versus slip relationship. The current study aims to establish the bond-slip model for CFRP sheets bonded to steel plate. To obtain the shear stress versus slippage relationship, a series of double strap tension type bond tests were conducted. This paper reports on the findings of the experimental studies. The strain and stress distributions measured in the specimens for two different bond lengths. The results show a preliminary bi-linear bond-slip model may be adopted for CFRP sheet bonded with steel plate.

Computational fluid dynamic analysis of intracranial aneurysmal bleb formation

Background The management of unruptured aneurysms is controversial with the decision to treat influenced by aneurysm characteristics including size and morphology. Aneurysmal bleb formation is thought to be associated with an increased risk of rupture. Objective To correlate computational fluid dynamic (CFD) indices with bleb formation. Methods Anatomical models were constructed from three-dimensional rotational angiogram (3DRA) data in 27 patients with cerebral aneurysms harbouring single blebs. Additional models representing the aneurysm before bleb formation were constructed by digitally removing the bleb. We characterised haemodynamic features of models both with and without the bleb using CFDs. Flow structure, wall shear stress (WSS), pressure and oscillatory shear index (OSI) were analysed. Results There was a statistically significant association between bleb location at or adjacent to the point of maximal WSS (74.1%, p=0.019), irrespective of rupture status. Aneurysmal blebs were related to the inflow or outflow jet in 88.9% of cases (p<0.001) whilst 11.1% were unrelated. Maximal wall pressure and OSI were not significantly related to bleb location. The bleb region attained a lower WSS following its formation in 96.3% of cases (p<0.001) and was also lower than the average aneurysm WSS in 86% of cases (p<0.001). Conclusion Cerebral aneurysm blebs generally form at or adjacent to the point of maximal WSS and are aligned with major flow structures. Wall pressure and OSI do not contribute to determining bleb location. The measurement of WSS using CFD models may potentially predict bleb formation and thus improve the assessment of rupture risk in unruptured aneurysms.

Atomistic simulation of interfacial behaviour in graphene-polymer nanocomposites

Graphene–polymer nanocomposites have promising properties as new structural and functional materials. The remarkable mechanical property enhancement in these nanocomposites is generally attributed to exceptional mechanical property of graphene and possible load transfer between graphene and polymer matrix. However, the underlying strengthening and toughening mechanisms have not been well understood. In this work, the interfacial behavior of graphene-polyethylene (PE) was investigated using molecular dynamics (MD) method. The interfacial shear force (ISF) and interfacial shear stress (ISS) between graphene and PE matrix were evaluated, taking into account graphene size, the number of graphene layers and the structural defects in graphene. MD results show that the ISS at graphene-PE interface mainly distributes at each end of the graphene nanofiller within the range of 1 nm, and much larger than that at carbon nanotube (CNT)-PE interface. Moreover, it was found that the ISS at graphene-PE interface is sensitive to the layer number.

Effects of scaffold architecture on mechanical characteristics and osteoblast response to static and perfusion bioreactor cultures

Tissue engineering focuses on the repair and regeneration of tissues through the use of biodegradable scaffold systems that structurally support regions of injury whilst recruiting and/or stimulating cell populations to rebuild the target tissue. Within bone tissue engineering, the effects of scaffold architecture on cellular response have not been conclusively characterized in a controlled-density environment. We present a theoretical and practical assessment of the effects of polycaprolactone (PCL) scaffold architectural modifications on mechanical and flow characteristics as well as MC3T3-E1 preosteoblast cellular response in an in vitro static plate and custom-designed perfusion bioreactor model. Four scaffold architectures were contrasted, which varied in inter-layer lay-down angle and offset between layers, whilst maintaining a structural porosity of 60 ± 5%. We established that as layer angle was decreased (90° vs. 60°) and offset was introduced (0 vs. 0.5 between layers), structural stiffness, yield stress, strength, pore size and permeability decreased, whilst computational fluid dynamics-modeled wall shear stress was increased. Most significant effects were noted with layer offset. Seeding efficiencies in static culture were also dramatically increased due to offset (~45% to ~86%), with static culture exhibiting a much higher seeding efficiency than perfusion culture. Scaffold architecture had minimal effect on cell response in static culture. However, architecture influenced osteogenic differentiation in perfusion culture, likely by modifying the microfluidic environment.

Impinging jet simulation of stationary downburst flow over topography

A non-translating, long duration thunderstorm downburst has been simulated experimentally and numerically by modelling a spatially stationary steady flow impinging air jet. Velocity profiles were shown to compare well with an upper-bound of velocity measurements reported for full-scale microbursts. Velocity speed-up over a range of topographic features in simulated downburst flow was also tested with comparisons made to previous work in a similar flow, and also boundary layer wind tunnel experiments. It was found that the amplification measured above the crest of topographic features in simulated downburst flow was up to 35% less than that observed in boundary layer flow for all shapes tested. From the computational standpoint we conclude that the Shear Stress Transport (SST) model performs the best from amongst a range of eddy-viscosity and second moment closures tested for modelling the impinging jet flow.

Modelling of arc welding : the importance of including the arc plasma in the computational domain

We present results of computational simulations of tungsten-inert-gas and metal-inert-gas welding. The arc plasma and the electrodes (including the molten weld pool when necessary) are included self-consistently in the computational domain. It is shown, using three examples, that it would be impossible to accurately estimate the boundary conditions on the weld-pool surface without including the arc plasma in the computational domain. First, we show that the shielding gas composition strongly affects the properties of the arc that influence the weld pool: heat flux density, current density, shear stress and arc pressure at the weld-pool surface. Demixing is found to be important in some cases. Second, the vaporization of the weld-pool metal and the diffusion of the metal vapour into the arc plasma are found to decrease the heat flux density and current density to the weld pool. Finally, we show that the shape of the wire electrode in metal-inert-gas welding has a strong influence on flow velocities in the arc and the pressure and shear stress at the weld-pool surface. In each case, we present evidence that the geometry and depth of the weld pool depend strongly on the properties of the arc.

On the biomechanical significance of inter-lamellar interfaces in the intervertebral disc

The intervertebral disc (IVD) is a unique soft tissue structure which provides structural support and flexibility in the axial skeleton of vertebrates. From a structural perspective, the disc behaves somewhat like a thick walled pressure vessel, where the walls are comprised of a series of composite annular rings (lamellae). However, a prior study (Marchand and Ahmed, 1990) found a high proportion of circumferentially discontinuous lamellae in human lumbar IVDs. The presence of these discontinuities raises important structural questions, because discontinuous lamellae cannot withstand high nucleus pressures via the generation of circumferential (hoop) stress. A possible alternative mechanism may be that inter-lamellar cohesion allows shear stress transfer between adjacent annular layers. The aim of the present study was therefore to investigate the importance of inter-lamellar shear resistance in the intervertebral disc. This work found that inter-lamellar shear resistance has a strong influence on the compressive stiffness of the intervertebral disc, with a change in interface condition from tied (no slip) to frictionless (no shear resistance) reducing disc compressive stiffness by 40%. However, it appears that substantial inter-lamellar shear resistance is present in the bovine tail disc. Decreases in inter-lamellar shear resistance due to degradation of bridging collagenous or elastic fibre structures could therefore be an important part of the process of disc degeneration.

Flow evaluation and hemolysis analysis of BVAD centrifugal blood pump by computational fluids dynamics

Computational fluid dynamics (CFD) and particle image velocimetry (PIV) are commonly used techniques to evaluate the flow characteristics in the development stage of blood pumps. CFD technique allows rapid change to pump parameters to optimize the pump performance without having to construct a costly prototype model. These techniques are used in the construction of a bi-ventricular assist device (BVAD) which combines the functions of LVAD and RVAD in a compact unit. The BVAD construction consists of two separate chambers with similar impellers, volutes, inlet and output sections. To achieve the required flow characteristics of an average flow rate of 5 l/min and different pressure heads (left – 100mmHg and right – 20mmHg), the impellers were set at different rotating speeds. From the CFD results, a six-blade impeller design was adopted for the development of the BVAD. It was also observed that the fluid can flow smoothly through the pump with minimum shear stress and area of stagnation which are related to haemolysis and thrombosis. Based on the compatible Reynolds number the flow through the model was calculated for the left and the right pumps. As it was not possible to have both the left and right chambers in the experimental model, the left and right pumps were tested separately.

The effect of thermal radiation on the natural convection boundary layer flow over a wavy horizontal surface

In this article, natural convection boundary layer flow is investigated over a semi-infinite horizontal wavy surface. Such an irregular (wavy) surface is used to exchange heat with an external radiating fluid which obeys Rosseland diffusion approximation. The boundary layer equations are cast into dimensionless form by introducing appropriate scaling. Primitive variable formulations (PVF) and stream function formulations (SFF) are independently used to transform the boundary layer equations into convenient form. The equations obtained from the former formulations are integrated numerically via implicit finite difference iterative scheme whereas equations obtained from lateral formulations are simulated through Keller-box scheme. To validate the results, solutions produced by above two methods are compared graphically. The main parameters: thermal radiation parameter and amplitude of the wavy surface are discussed categorically in terms of shear stress and rate of heat transfer. It is found that wavy surface increases heat transfer rate compared to the smooth wall. Thus optimum heat transfer is accomplished when irregular surface is considered. It is also established that high amplitude of the wavy surface in the boundary layer leads to separation of fluid from the plate.

Unsteady flow over a stretching surface with a magnetic field in a rotating fluid

The effect of the magnetic field on the unsteady flow over a stretching surface in a rotating fluid has been studied. The unsteadiness in the flow field is due to the time-dependent variation of the velocity of the stretching surface and the angular velocity of the rotating fluid. The Navier-Stokes equations and the energy equation governing the flow and the heat transfer admit a self-similar solution if the velocity of the stretching surface and the angular velocity of the rotating fluid vary inversely as a linear function of time. The resulting system of ordinary differential equations is solved numerically using a shooting method. The rotation parameter causes flow reversal in the component of the velocity parallel to the strerching surface and the magnetic field tends to prevent or delay the flow reversal. The surface shear stresses dong the stretching surface and in the rotating direction increase with the rotation parameter, but the surface heat transfer decreases. On the other hand, the magnetic field increases the surface shear stress along the stretching surface, but reduces the surface shear stress in the rotating direction and the surface heat transfer. The effect of the unsteady parameter is more pronounced on the velocity profiles in the rotating direction and temperature profiles.

Size-dependent elasticity of nanocrystalline titania

Synchrotron-based high-pressure x-ray diffraction measurements indicate that compressibility, a fundamental materials property, can have a size-specific minimum value. The bulk modulus of nanocrystalline titania has a maximum at particle size of 15 nm. This can be explained by dislocation behavior because very high dislocation contents can be achieved when shear stress induced within nanoparticles counters the repulsion between dislocations. As particle size decreases, compression increasingly generates dislocation networks hardened by overlap of strain fields that shield intervening regions from external pressure. However, when particles become too small to sustain high dislocation concentrations, elastic stiffening declines. The compressibility has a minimum at intermediate sizes.

Cardiovascular diseases and vulnerable plaques: data, modeling, predictions and clinical applications PREFACE

Introduction Two symposia on “cardiovascular diseases and vulnerable plaques” Cardiovascular disease (CVD) is the leading cause of death worldwide. Huge effort has been made in many disciplines including medical imaging, computational modeling, bio- mechanics, bioengineering, medical devices, animal and clinical studies, population studies as well as genomic, molecular, cellular and organ-level studies seeking improved methods for early detection, diagnosis, prevention and treatment of these diseases [1-14]. However, the mechanisms governing the initiation, progression and the occurrence of final acute clinical CVD events are still poorly understood. A large number of victims of these dis- eases who are apparently healthy die suddenly without prior symptoms. Available screening and diagnostic methods are insufficient to identify the victims before the event occurs [8,9]. Most cardiovascular diseases are associated with vulnerable plaques. A grand challenge here is to develop new imaging techniques, predictive methods and patient screening tools to identify vulnerable plaques and patients who are more vulnerable to plaque rupture and associated clinical events such as stroke and heart attack, and recommend proper treatment plans to prevent those clinical events from happening. Articles in this special issue came from two symposia held recently focusing on “Cardio-vascular Diseases and Vulnerable Plaques: Data, Modeling, Predictions and Clinical Applications.” One was held at Worcester Polytechnic Institute (WPI), Worcester, MA, USA, July 13-14, 2014, right after the 7th World Congress of Biomechanics. This symposium was endorsed by the World Council of Biomechanics, and partially supported by a grant from NIH-National Institute of Biomedical Image and Bioengineering. The other was held at Southeast University (SEU), Nanjing, China, April 18-20, 2014.

Flow pattern analysis in a highly stenotic patient-specific carotid bifurcation model using a turbulence model

The aim of this study is to investigate the blood flow pattern in carotid bifurcation with a high degree of luminal stenosis, combining in vivo magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). A newly developed two-equation transitional model was employed to evaluate wall shear stress (WSS) distribution and pressure drop across the stenosis, which are closely related to plaque vulnerability. A patient with an 80% left carotid stenosis was imaged using high resolution MRI, from which a patient-specific geometry was reconstructed and flow boundary conditions were acquired for CFD simulation. A transitional model was implemented to investigate the flow velocity and WSS distribution in the patient-specific model. The peak time-averaged WSS value of approximately 73Pa was predicted by the transitional flow model, and the regions of high WSS occurred at the throat of the stenosis. High oscillatory shear index values up to 0.50 were present in a helical flow pattern from the outer wall of the internal carotid artery immediately after the throat. This study shows the potential suitability of a transitional turbulent flow model in capturing the flow phenomena in severely stenosed carotid arteries using patient-specific MRI data and provides the basis for further investigation of the links between haemodynamic variables and plaque vulnerability. It may be useful in the future for risk assessment of patients with carotid disease.