Mathematical modelling of competitive LDL/VLDL binding and uptake by hepatocytes

Elevated levels of low-density-lipoprotein cholesterol (LDL-C) in the plasma are a well-established risk factor for the development of coronary heart disease. Plasma LDL-C levels are in part determined by the rate at which LDL particles are removed from the bloodstream by hepatic uptake. The uptake of LDL by mammalian liver cells occurs mainly via receptor-mediated endocytosis, a process which entails the binding of these particles to specific receptors in specialised areas of the cell surface, the subsequent internalization of the receptor-lipoprotein complex, and ultimately the degradation and release of the ingested lipoproteins' constituent parts. We formulate a mathematical model to study the binding and internalization (endocytosis) of LDL and VLDL particles by hepatocytes in culture. The system of ordinary differential equations, which includes a cholesterol-dependent pit production term representing feedback regulation of surface receptors in response to intracellular cholesterol levels, is analysed using numerical simulations and steady-state analysis. Our numerical results show good agreement with in vitro experimental data describing LDL uptake by cultured hepatocytes following delivery of a single bolus of lipoprotein. Our model is adapted in order to reflect the in vivo situation, in which lipoproteins are continuously delivered to the hepatocyte. In this case, our model suggests that the competition between the LDL and VLDL particles for binding to the pits on the cell surface affects the intracellular cholesterol concentration. In particular, we predict that when there is continuous delivery of low levels of lipoproteins to the cell surface, more VLDL than LDL occupies the pit, since VLDL are better competitors for receptor binding. VLDL have a cholesterol content comparable to LDL particles; however, due to the larger size of VLDL, one pit-bound VLDL particle blocks binding of several LDLs, and there is a resultant drop in the intracellular cholesterol level. When there is continuous delivery of lipoprotein at high levels to the hepatocytes, VLDL particles still out-compete LDL particles for receptor binding, and consequently more VLDL than LDL particles occupy the pit. Although the maximum intracellular cholesterol level is similar for high and low levels of lipoprotein delivery, the maximum is reached more rapidly when the lipoprotein delivery rates are high. The implications of these results for the design of in vitro experiments is discussed.

New mathematical modeling approach for predicting microbial inactivation by high hydrostatic pressure

A new primary model based on a thermodynamically consistent first-order kinetic approach was constructed to describe non-log-linear inactivation kinetics of pressure-treated bacteria. The model assumes a first-order process in which the specific inactivation rate changes inversely with the square root of time. The model gave reasonable fits to experimental data over six to seven orders of magnitude. It was also tested on 138 published data sets and provided good fits in about 70% of cases in which the shape of the curve followed the typical convex upward form. In the remainder of published examples, curves contained additional shoulder regions or extended tail regions. Curves with shoulders could be accommodated by including an additional time delay parameter and curves with tails shoulders could be accommodated by omitting points in the tail beyond the point at which survival levels remained more or less constant. The model parameters varied regularly with pressure, which may reflect a genuine mechanistic basis for the model. This property also allowed the calculation of (a) parameters analogous to the decimal reduction time D and z, the temperature increase needed to change the D value by a factor of 10, in thermal processing, and hence the processing conditions needed to attain a desired level of inactivation; and (b) the apparent thermodynamic volumes of activation associated with the lethal events. The hypothesis that inactivation rates changed as a function of the square root of time would be consistent with a diffusion-limited process.

Modelling the uptake and growth kinetics of Penicillium glabrum in a tannic acid-containing solid-state fermentation for tannase production

A mathematical growth model for the batch solid-state fermentation process for fungal tannase production was developed and tested experimentally. The unstructured model describes the uptake and growth kinetics of Penicillium glabrum in an impregnated polyurethane foam substrate system. In general, good agreement between the experimental data and model simulations was obtained. Biomass, tannase and spore production are described by logistic kinetics with a time delay between biomass production and tannase and spore formation. Possible induction mechanisms for the latter are proposed. Hydrolysis of tannic acid, the main carbon source in the substrate system, is reasonably well described with Michaelis-Menten kinetics with time-varying enzyme concentration but a more complex reaction mechanism is suspected. The metabolism of gallic acid, a tannase-hydrolysis product of tannic acid, was shown to be growth limiting during the main growth phase. (c) 2004 Elsevier Ltd. All rights reserved.

Modelling and validation of Lactobacillus plantarum fermentations in cereal-based media with different sugar concentrations and buffering capacities

An unstructured mathematical model is proposed to describe the fermentation kinetics of growth, lactic acid production, pH and sugar consumption by Lactobacillus plantarum as a function of the buffering capacity and initial glucose concentration of the culture media. Initially the experimental data of L plantarum fermentations in synthetic media with different buffering capacity and glucose were fitted to a set of primary models. Later the parameters obtained from these models were used to establish mathematical relationships with the independent variables tested. The models were validated with 6 fermentations of L. plantarum in different cereal-based media. In most cases the proposed models adequately describe the biochemical changes taking place during fermentation and are a promising approach for the formulation of cereal-based probiotic foods. (C) 2008 Elsevier B.V. All rights reserved.

Kinetic modeling of the generation of 2-and 3-methylbutanal in a heated extract of beef liver

Quantitative control of aroma generation during the Maillard reaction presents great scientific and industrial interest. Although there have been many studies conducted in simplified model systems, the results are difficult to apply to complex food systems, where the presence of other components can have a significant impact. In this work, an aqueous extract of defatted beef liver was chosen as a simplified food matrix for studying the kinetics of the Mallard reaction. Aliquots of the extract were heated under different time and temperature conditions and analyzed for sugars, amino acids, and methylbutanals, which are important Maillard-derived aroma compounds formed in cooked meat. Multiresponse kinetic modeling, based on a simplified mechanistic pathway, gave a good fit with the experimental data, but only when additional steps were introduced to take into account the interactions of glucose and glucose-derived intermediates with protein and other amino compounds. This emphasizes the significant role of the food matrix in controlling the Maillard reaction.

Phonological similarity effects without a phonological store: an individual differences model

The feature model of immediate memory (Nairne, 1990) is applied to an experiment testing individual differences in phonological confusions amongst a group (N=100) of participants performing a verbal memory test. By simulating the performance of an equivalent number of “pseudo-participants” the model fits both the mean performance and the variability within the group. Experimental data show that high-performing individuals are significantly more likely to demonstrate phonological confusions than low performance individuals and this is also true of the model, despite the model’s lack of either an explicit phonological store or a performance-linked strategy shift away from phonological storage. It is concluded that a dedicated phonological store is not necessary to explain the basic phonological confusion effect, and the reduction in such an effect can also be explained without requiring a change in encoding or rehearsal strategy or the deployment of a different storage buffer.

Alternative explanations for rising dissolved organic carbon export from organic soils

Since 1988, there has been, on average, a 91% increase in dissolved organic carbon (DOC) concentrations of UK lakes and streams in the Acid Waters Monitoring Network (AWMN). Similar DOC increases have been observed in surface waters across much of Europe and North America. Much of the debate about the causes of rising DOC has, as in other studies relating to the carbon cycle, focused on factors related to climate change. Data from our peat-core experiments support an influence of climate on DOC, notably an increase in production with temperature under aerobic, and to a lesser extent anaerobic, conditions. However, we argue that climatic factors may not be the dominant drivers of DOC change. DOC solubility is suppressed by high soil water acidity and ionic strength, both of which have decreased as a result of declining sulphur deposition since the 1980s, augmented during the 1990s in the United Kingdom by a cyclical decline in sea-salt deposition. Our observational and experimental data demonstrate a clear, inverse and quantitatively important link between DOC and sulphate concentrations in soil solution. Statistical analysis of 11 AWMN lakes suggests that rising temperature, declining sulphur deposition and changing sea-salt loading can account for the majority of the observed DOC trend. This combination of evidence points to the changing chemical composition of atmospheric deposition, particularly the substantial reduction in anthropogenic sulphur emissions during the last 20 years, as a key cause of rising DOC. The implications of rising DOC export for the carbon cycle will be very different if linked primarily to decreasing acid deposition, rather than to changes in climate, suggesting that these systems may be recovering rather than destabilising.

Ocean acoustic models for low frequency propagation in 2D and 3D environments

In this paper we are mainly concerned with the development of efficient computer models capable of accurately predicting the propagation of low-to-middle frequency sound in the sea, in axially symmetric (2D) and in fully 3D environments. The major physical features of the problem, i.e. a variable bottom topography, elastic properties of the subbottom structure, volume attenuation and other range inhomogeneities are efficiently treated. The computer models presented are based on normal mode solutions of the Helmholtz equation on the one hand, and on various types of numerical schemes for parabolic approximations of the Helmholtz equation on the other. A new coupled mode code is introduced to model sound propagation in range-dependent ocean environments with variable bottom topography, where the effects of an elastic bottom, of volume attenuation, surface and bottom roughness are taken into account. New computer models based on finite difference and finite element techniques for the numerical solution of parabolic approximations are also presented. They include an efficient modeling of the bottom influence via impedance boundary conditions, they cover wide angle propagation, elastic bottom effects, variable bottom topography and reverberation effects. All the models are validated on several benchmark problems and versus experimental data. Results thus obtained were compared with analogous results from standard codes in the literature.

Enantiospecific adsorption of alanine on the chiral Cu{531} surface

We have studied enantiospecific differences in the adsorption of (S)- and (R)-alanine on Cu{531}R using low-energy electron diffraction (LEED), X-ray photoelectron spectroscopy, and near edge X-ray absorption fine structure (NEXAFS) spectroscopy. At saturation coverage, alanine adsorbs as alaninate forming a p(1 4) superstructure. LEED shows a significantly higher degree of long-range order for the S than for the R enantiomer. Also carbon K-edge NEXAFS spectra show differences between (S)- and (R)-alanine in the variations of the ð resonance when the linear polarization vector is rotated within the surface plane. This indicates differences in the local adsorption geometries of the molecules, most likely caused by the interaction between the methyl group and the metal surface and/or intermolecular hydrogen bonds. Comparison with model calculations and additional information from LEED and photoelectron spectroscopy suggest that both enantiomers of alaninate adsorb in two different orientations associated with triangular adsorption sites on {110} and {311} microfacets of the Cu{531} surface. The experimental data are ambiguous as to the exact difference between the local geometries of the two enantiomers. In one of two models that fit the data equally well, significantly more (R)-alaninate molecules are adsorbed on {110} sites than on {311} sites whereas for (S)-alaninate the numbers are equal. The enantiospecific differences found in these experiments are much more pronounced than those reported from other ultrahigh vacuum techniques applied to similar systems.

Computer simulations of structures, energetics and dynamics of myoglobin center dot center dot center dot ligand complexes

Myoglobin has been studied in considerable detail using different experimental and computational techniques over the past decades. Recent developments in time-resolved spectroscopy have provided experimental data amenable to detailed atomistic simulations. The main theme of the present review are results on the structures, energetics and dynamics of ligands ( CO, NO) interacting with myoglobin from computer simulations. Modern computational methods including free energy simulations, mixed quantum mechanics/molecular mechanics simulations, and reactive molecular dynamics simulations provide insight into the dynamics of ligand dynamics in confined spaces complementary to experiment. Application of these methods to calculate and understand experimental observations for myoglobin interacting with CO and NO are presented and discussed.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

Structure–activity relations (SARs) for gas-phase reactions of NO3, OH and O3 with alkenes: an update

We present updated structure-activity relations (SARs) for the prediction of rate coefficients for gas-phase reactions with alkenes of the major atmospheric oxidants NO3, OH and O-3. Such SARs provide one way of incorporating essential information about reactivity into atmospheric models. Rate coefficients obtained from correlations relating the logarithms of the rate coefficients to the energies of the highest occupied molecular orbitals (HOMOs) of the alkenes were used to refine the SARs. SARs have an advantage for the user over the direct application of the correlations in that knowledge of the structure of the alkene of interest is sufficient to estimate rate coefficients, and no quantum-mechanical calculations need to be performed. A comparison of the values predicted by the SARs with experimental data where they exist allowed us to assess the reliability of our method.

Correlations for gas-phase reactions of NO3, OH and O3 with alkenes: an update

Methods are developed for predicting rate coefficients for reactions of initiators of tropospheric oxidation with unsaturated compounds that are abundant in the atmosphere; prognostic tools of this kind are essential for atmospheric chemists and modellers. To pursue the aim of exploring such tools, the kinetics of reactions of NO3, OH and O-3 with a series of alkenes are examined for correlations relating the logarithms of the rate coefficients to the energies of the highest occupied molecular orbitals (HOMOs) of the alkenes. A comparison of the values predicted by the correlations with experimental data (where the latter exist) allowed us to assess the reliability of our method. We used a series of theoretical methods to calculate the HOMO energies, and found that higher computational effort improves the agreement of the predicted rate coefficients with experimental values, especially for reactions of NO3 with alkenes that possess vinyllic halogen substituents. As a consequence, it is expedient to suggest new correlations to replace those presented by us and others that were based on the lower level of theory. We propose the following correlations for the reactions of NO3, OH and O-3 with alkenes: ln(k(NO3)/cm(3) molecule(-1) s(-1)) = 6.40(E-HOMO/eV) + 31.69, ln(k(OH)/cm(3) molecule(-1) s(-1)) = 1.21 (E-HOMO/eV)-12.34 and ln(k(O3)/cm(3) molecule(-1) s(-1)) = 3.28(E-HOMO/eV)-6.78. These new correlations have been developed using the larger experimental data sets now available, and the impact of the extended data on the quality of the correlations is examined in the paper. Atmospheric lifetimes have been calculated from both experimental and estimated rate coefficients to provide an overview of removal efficiencies for different classes of alkenes with respect to oxidative processes initiated by NO3, OH and O-3. A figure is presented to show the spatial scales over which alkenes may survive transport in competition with attack by NO3, OH and O-3. Removal by NO3 or OH is always more important than removal by O-3, and reactions with NO3 dominate for scales up to a few hundred metres.

Gas-phase rate coefficients for the reactions of O(3P), S(3P), Se(3P), and Te(3P) with alkenes: application of perturbation frontier molecular orbital theory, correlations, and structure-activity relations (SARs)

The kinetics of the reactions of the atoms O(P-3), S(P-3), Se(P-3), and Te((3)p) with a series of alkenes are examined for correlations relating the logarithms of the rate coefficients to the energies of the highest occupied molecular orbitals (HOMOs) of the alkenes. These correlations may be employed to predict rate coefficients from the calculated HOMO energy of any other alkene of interest. The rate coefficients obtained from the correlations were used to formulate structure-activity relations (SARs) for reactions of O((3)p), S(P-3), Se (P-3), and Te((3)p) with alkenes. A comparison of the values predicted by both the correlations and the SARs with experimental data where they exist allowed us to assess the reliability of our method. We demonstrate the applicability of perturbation frontier molecular orbital theory to gas-phase reactions of these atoms with alkenes. The correlations are apparently not applicable to reactions of C(P-3), Si(P-3), N(S-4), and Al(P-2) atoms with alkenes, a conclusion that could be explained in terms of a different mechanism for reaction of these atoms.

A square-planar nickel(II) monoradical complex with a bis(salicylidene)diamine ligand

The new square-planar Ni-II-N2O2 complex [Ni(L-Me)] (1(Me)), where L-Me, stands for the dianionic phenolato form of N,N'bis(3,5-di-tert-butyl-salicylidene)-4,5-dimethyl-1,2-phenyl- enediamine ((LH2)-L-Me), has been synthesised and fully characterised. X-ray crystallography was also used for the characterisation. The electrochemical one-electron oxidation of 1(Me) produces the thermally stable (within the temperature range 10-295 K) cationic species (1(Me))(+). The UV/Vis and X-band EPR experimental data, supported by DFT calculations, indicate that (1(Me))(+), is best described as a Ni-II monoradical complex and, thus, does NOT exist in a Ni-III ground state, in contrast to its demethylated counterpart [Ni(L-H)](+) (1(H))(+) below 170 K.