Modeling the diffusion-erosion crossover dynamics in drug release.

A computational model is proposed to investigate drug delivery systems in which erosion and diffusion mechanisms are participating in the drug release process. Our approach allowed us to analytically estimate the crossover point between those mechanisms through the value of the parameter b (b_{c}=1) and the scaling behavior of parameter τ on the Weibull function, exp[-(t/τ)^{b}], used to adjust drug release data in pharmaceutical literature. Numerical investigations on the size dependence of the characteristic release time τ found it to satisfy either linear or quadratic scaling relations on either erosive or diffusive regimes.

Residual entropy and waterlike anomalies in the repulsive one dimensional lattice gas.

The thermodynamics and kinetics of the one dimensional lattice gas with repulsive interaction are investigated using transfer matrix technique and Monte Carlo simulations. This simple model is shown to exhibit waterlike anomalies in density, thermal expansion coefficient, and self-diffusion. An unified description for the thermodynamic anomalies in this model is achieved based on the ground state residual entropy which appears in the model due to mixing entropy in a ground state phase transition.

Thermodynamic and dynamic anomalies in a one-dimensional lattice model of liquid water.

We investigate the occurrence of waterlike thermodynamic and dynamic anomalous behavior in a one dimensional lattice gas model. The system thermodynamics is obtained using the transfer matrix technique and anomalies on density and thermodynamic response functions are found. When the hydrogen bond (molecules separated by holes) is more attractive than the van der Waals interaction (molecules in contact) a transition between two fluid structures is found at null temperature and high pressure.

Solution of an associating lattice-gas model with density anomaly on a Husimi lattice.

We study a model of a lattice gas with orientational degrees of freedom which resemble the formation of hydrogen bonds between the molecules. In this model, which is the simplified version of the Henriques-Barbosa model, no distinction is made between donors and acceptors in the bonding arms. We solve the model in the grand-canonical ensemble on a Husimi lattice built with hexagonal plaquettes with a central site.

Molecular correlations and solvation in simple fluids.

We study the molecular correlations in a lattice model of a solution of a low-solubility solute, with emphasis on how the thermodynamics is reflected in the correlation functions. The model is treated in the Bethe-Guggenheim approximation, which is exact on a Bethe lattice (Cayley tree). The solution properties are obtained in the limit of infinite dilution of the solute.

Frustration and anomalous behavior in the Bell-Lavis model of liquid water.

We have reconsidered the Bell-Lavis model of liquid water and investigated its relation to its isotropic version, the antiferromagnetic Blume-Emery-Griffiths model on the triangular lattice. Our study was carried out by means of an exact solution on the sequential Husimi cactus. We show that the ground states of both models share the same topology and that fluid phases (gas and low- and high-density liquids) can be mapped onto magnetic phases (paramagnetic, antiferromagnetic, and dense paramagnetic, respectively).

Alternative hydrogen bond implementations produce opposite effects on collapse cooperativity of lattice homopolypeptide models.

We use complete enumeration of self-avoiding chains of up to N=26 monomers in two-dimensional lattices to investigate the effect of alternative implementations of backbone hydrogen bonds on the cooperativity of homopolypeptide collapse. Following a recent study on protein folding models, we use the square lattice with z=3 local conformations per monomer and lattice extensions containing diagonal steps which result in z=5 or z=7 and assume that only a subset of zh View Article and Find Full Text PDF

Entropy reduction effect imposed by hydrogen bond formation on protein folding cooperativity: evidence from a hydrophobic minimalist model.

Conformational restrictions imposed by hydrogen bond formation during protein folding are investigated by Monte Carlo simulations of a non-native-centric, two-dimensional, hydrophobic model in which the formation of favorable contacts is coupled to an effective reduction in lattice coordination. This scheme is intended to mimic the requirement that polar backbone groups of real proteins must form hydrogen bonds concomitantly to their burial inside the apolar protein core. In addition to the square lattice, with z=3 conformations per monomer, we use extensions in which diagonal step vectors are allowed, resulting in z=5 and z=7.

Relevance of structural segregation and chain compaction for the thermodynamics of folding of a hydrophobic protein model.

The relevance of inside-outside segregation and chain compaction for the thermodynamics of folding of a hydrophobic protein model is probed by complete enumeration of two-dimensional chains of up to 18 monomers in the square lattice. The exact computation of Z scores for uniquely designed sequences confirms that Z tends to decrease linearly with sigma square root of N, as previously suggested by theoretical analysis and Monte Carlo simulations, where sigma, the standard deviation of the number of contacts made by different monomers in the target structure, is a measure of structural segregation and N is the chain length. The probability that the target conformation is indeed the unique global energy minimum of the designed sequence is found to increase dramatically with sigma, approaching unity at maximal segregation.

Non-native interactions, effective contact order, and protein folding: a mutational investigation with the energetically frustrated hydrophobic model.

By Monte Carlo simulations, we explored the effect of single mutations on the thermodynamics and kinetics of the folding of a two-dimensional, energetically frustrated, hydrophobic protein model. Phi-Value analysis, corroborated by simulations beginning from given sets of judiciously chosen initial contacts, suggests that the transition state of the model consists of a limited region of the native structure, that is, a folding nucleus. It seems that the most important contacts in the transition state (large and positive Phi) are not the ones with the highest contact order, because in this case the entropic cost of their formation would be too high, but exactly the ones that decrease the entropic cost of difficult contacts, reducing their effective contact order.