3D multi-agent models for protein release from PLGA spherical particles with complex inner morphologies

In order to better understand and predict the release of proteins from bioerodible micro- or nanospheres, it is important to know the influences of different initial factors on the release mechanisms. Often though it is difficult to assess what exactly is at the origin of a certain dissolution profile. We propose here a new class of fine-grained multi-agent models built to incorporate increasing complexity, permitting the exploration of the role of different parameters, especially that of the internal morphology of the spheres, in the exhibited release profile. This approach, based on Monte-Carlo (MC) and Cellular Automata (CA) techniques, has permitted the testing of various assumptions and hypotheses about several experimental systems of nanospheres encapsulating proteins. Results have confirmed that this modelling approach has increased the resolution over the complexity involved, opening promising perspectives for future developments, especially complementing in vitro experimentation

Influence of a knot on the strength of a polymer strand

Many experiments have been done to determine the relative strength of different knots, and these show that the break in a knotted rope almost invariably occurs at a point just outside the `entrance' to the knot. The influence of knots on the properties of polymers has become of great interest, in part because of their effect on mechanical properties. Knot theory applied to the topology of macromolecules indicates that the simple trefoil or `overhand' knot is likely to be present with high probability in any long polymer strand. Fragments of DNA have been observed to contain such knots in experiments and computer simulations. Here we use {\it ab initio} computational methods to investigate the effect of a trefoil knot on the breaking strength of a polymer strand. We find that the knot weakens the strand significantly, and that, like a knotted rope, it breaks under tension at the entrance to the knot.Comment: 3 pages, 4 figure

Report of first recurrent glioma patients examined with PET-MRI prior to re-irradiation

Background and purpose The advantage of combined PET-MRI over sequential PET and MRI is the high spatial conformity and the absence of time delay between the examinations. The benefit of this technique for planning of re-irradiation (re-RT) treatment is unkown yet. Imaging data from a phase 1 trial of re-RT for recurrent glioma was analysed to assess whether planning target volumes and treatment margins in glioma re-RT can be adjusted by PET-MRI with rater independent PET based biological tumour volumes (BTVs). Patients and methods Combined PET-MRI with the tracer O-(2-F-18-fluoroethyl)-1-tyrosine (F-18-FET) prior to re-RT was performed in recurrent glioma patients in a phase I trial. GTVs including all regions suspicious of tumour on contrast enhanced MRI were delineated by three experienced radiation oncologists and included into MRI based consensus GTVs (mRGTVs). BTVs were semiautomatically delineated with a fixed threshold of 1.6 x background activity. Corresponding BTVs and mRGTVs were fused into union volume RET-NARGIVs. The Sorensen Dice coefficient and the conformity index were used to assess the geometric overlap of the BTVs with the mRGTVs. A recurrence pattern analysis was performed based on the original planning target volumes (PTVs = GTV + 10 mm margin or 5 mm in one case) and the RET-NARGTVs with margins of 10, 8, 5 and 3 mm. Results Seven recurrent glioma patients, who received PET-MRI prior to re-RT, were included into the present planning study. At the time of re-RT, patients were in median 54 years old and had a median Karnofsky Performance Status (KPS) score of 80. Median post-recurrence survival after the beginning of re-RT was 13 months. Concomitant bevacizumab therapy was applied in six patients and one patient received chemoradiation with temozolomide. Median GTV volumes of the three radiation oncologists were 35.0, 37.5 and 40.5 cubic centimeters (cc) and median (MR)GTV volume 41.8 cc. Median BTV volume was 36.6 cc and median (PET-MR)GTV volume 59.3 cc. The median Sorensen-Dice coefficient for the comparison between (MR)GTV and BTV was 0.61 and the median conformity index 0.44. Recurrence pattern analysis revealed two central, two in-field and one distant recurrence within both, the original PTV, as well as the (PER-MR)GTV with a reduced margin of 3 mm. Conclusion PET-MRI provides radiation treatment planning imaging with high spatial and timely conformity for high-grade glioma patients treated with re-RT with potential advancements for target volume delineation. Prospective randomised trials are warranted to further investigate the treatment benefits of PET-MRI based re-RT planning

Parallel Excluded Volume Tempering for Polymer Melts

We have developed a technique to accelerate the acquisition of effectively uncorrelated configurations for off-lattice models of dense polymer melts which makes use of both parallel tempering and large scale Monte Carlo moves. The method is based upon simulating a set of systems in parallel, each of which has a slightly different repulsive core potential, such that a thermodynamic path from full excluded volume to an ideal gas of random walks is generated. While each system is run with standard stochastic dynamics, resulting in an NVT ensemble, we implement the parallel tempering through stochastic swaps between the configurations of adjacent potentials, and the large scale Monte Carlo moves through attempted pivot and translation moves which reach a realistic acceptance probability as the limit of the ideal gas of random walks is approached. Compared to pure stochastic dynamics, this results in an increased efficiency even for a system of chains as short as $N = 60$ monomers, however at this chain length the large scale Monte Carlo moves were ineffective. For even longer chains the speedup becomes substantial, as observed from preliminary data for $N = 200$

The electric double layer has a life of its own

Using molecular dynamics simulations with recently developed importance sampling methods, we show that the differential capacitance of a model ionic liquid based double-layer capacitor exhibits an anomalous dependence on the applied electrical potential. Such behavior is qualitatively incompatible with standard mean-field theories of the electrical double layer, but is consistent with observations made in experiment. The anomalous response results from structural changes induced in the interfacial region of the ionic liquid as it develops a charge density to screen the charge induced on the electrode surface. These structural changes are strongly influenced by the out-of-plane layering of the electrolyte and are multifaceted, including an abrupt local ordering of the ions adsorbed in the plane of the electrode surface, reorientation of molecular ions, and the spontaneous exchange of ions between different layers of the electrolyte close to the electrode surface. The local ordering exhibits signatures of a first-order phase transition, which would indicate a singular charge-density transition in a macroscopic limit

Mathematical models for estimating effective diffusion parameters of spherical drug delivery devices

Mathematical modeling of drug delivery is of increasing academic and industrial importance in manyaspects. In this paper, we propose an optimization approach for the estimation of the parameters characterizing the diffusion process of a drug from a spherical porous polymer device to an external finite volume. The approach is based on a nonlinear least-squares method and a novel mathematical model which takes into consideration both boundary layer effect and initial burst phenomenon. Ananalytical solution to the model is derived and a formula for the ratio of the mass released in a given time interval and the total mass released in infinite time is also obtained. The approach has been tested using experimental data of the diffusion of prednisolone 21-hemisuccinate sodium saltfrom spherical devices made of porous poly(2-hydroxyethyl methacrylate) hydrogels. The effectiveness and accuracy of the method are well demonstrated by the numerical results. The model was used to determine the diffusion parameters including the effective diffusion coefficient of the drug from a series of devices that vary in both the porous structure and the drug loading levels. The computed diffusion parameters are discussed in relation to the physical properties of the devices

Modulation of functional pendant chains within poly(ethylene glycol) hydrogels for refined control of protein release

Hydrogels are highly attractive delivery vehicles for therapeutic proteins. Their innate biocompatibility, hydrophilicity and aqueous permeability allow stable encapsulation and release of proteins. The release rates also can be controlled simply by altering the crosslinking density of the polymeric network. However, the crosslinking density also influences the mechanical properties of hydrogels, generally opposite to the permeability. In addition, the release of larger proteins may be hindered below critically diminished porosity determined by the crosslinking density. Herein, the physical properties of the hydrogels are tuned by presenting functional pendant chains, independent of crosslinking density. Heterobifunctional poly(ethylene glycol) monomethacrylate (PEGMA) with various end functional groups is synthesized and copolymerized with PEG dimethacrylate (PEGDA) to engineer PEG hydrogels with pendant PEG chains. The pendant chains of the PEG hydrogels consisting of sulfonate, trimethylammonium chloride, and phenyl groups are utilized to provide negative charge, positive charge and hydrophobicity, respectively, to the hydrogels. The release rates of proteins with different isoelectric points are controlled in a wide range by the type and the density of functional pendant chains via electrostatic and hydrophobic interactions

Bioadhesive Controlled Metronidazole Release Matrix Based on Chitosan and Xanthan Gum

Metronidazole, a common antibacterial drug, was incorporated into a hydrophilic polymer matrix composed of chitosan xanthan gum mixture. Hydrogel formation of this binary chitosan-xanthan gum combination was tested for its ability to control the release of metronidazole as a drug model. This preparation (MZ-CR) was characterized by in vitro, ex vivo bioadhesion and in vivo bioavailability study. For comparison purposes a commercial extended release formulation of metronidazole (CMZ) was used as a reference. The in vitro drug-release profiles of metronidazole preparation and CMZ were similar in 0.1 M HCl and phosphate buffer pH 6.8. Moreover, metronidazole preparation and CMZ showed a similar detachment force to sheep stomach mucosa, while the bioadhesion of the metronidazole preparation was higher three times than CMZ to sheep duodenum. The results of in vivo study indicated that the absorption of metronidazole from the preparation was faster than that of CMZ. Also, MZ-CR leads to higher metronidazole Cmax and AUC relative to that of the CMZ. This increase in bioavailability might be explained by the bioadhesion of the preparation at the upper part of the small intestine that could result in an increase in the overall intestinal transit time. As a conclusion, formulating chitosan-xanthan gum mixture as a hydrophilic polymer matrix resulted in a superior pharmacokinetic parameters translated by better rate and extent of absorption of metronidazole

Use-Exposure Relationships of Pesticides for Aquatic Risk Assessment

Field-scale environmental models have been widely used in aquatic exposure assessments of pesticides. Those models usually require a large set of input parameters and separate simulations for each pesticide in evaluation. In this study, a simple use-exposure relationship is developed based on regression analysis of stochastic simulation results generated from the Pesticide Root-Zone Model (PRZM). The developed mathematical relationship estimates edge-of-field peak concentrations of pesticides from aerobic soil metabolism half-life (AERO), organic carbon-normalized soil sorption coefficient (KOC), and application rate (RATE). In a case study of California crop scenarios, the relationships explained 90–95% of the variances in the peak concentrations of dissolved pesticides as predicted by PRZM simulations for a 30-year period. KOC was identified as the governing parameter in determining the relative magnitudes of pesticide exposures in a given crop scenario. The results of model application also indicated that the effects of chemical fate processes such as partitioning and degradation on pesticide exposure were similar among crop scenarios, while the cross-scenario variations were mainly associated with the landscape characteristics, such as organic carbon contents and curve numbers. With a minimum set of input data, the use-exposure relationships proposed in this study could be used in screening procedures for potential water quality impacts from the off-site movement of pesticides