Pressure-induced protein-folding/unfolding kinetics

We use an off-lattice minimalist model to describe the effects of pressure in slowing down the folding/unfolding kinetics of proteins when subjected to increasingly larger pressures. The potential energy function used to describe the interactions between beads in the model includes the effects of pressure on the pairwise interaction of hydrophobic groups in water. We show that pressure affects the participation of contacts in the transition state. More significantly, pressure exponentially decreases the chain reconfigurational diffusion coefficient. These results are consistent with experimental results on the kinetics of pressure-denaturation of staphylococcal nuclease.

How optimization of potential functions affects protein folding.

The relationship between the optimization of the potential function and the foldability of theoretical protein models is studied based on investigations of a 27-mer cubic-lattice protein model and a more realistic lattice model for the protein crambin. In both the simple and the more complicated systems, optimization of the energy parameters achieves significant improvements in the statistical-mechanical characteristics of the systems and leads to foldable protein models in simulation experiments. The foldability of the protein models is characterized by their statistical-mechanical properties--e.g., by the density of states and by Monte Carlo folding simulations of the models. With optimized energy parameters, a high level of consistency exists among different interactions in the native structures of the protein models, as revealed by a correlation function between the optimized energy parameters and the native structure of the model proteins. The results of this work are relevant to the design of a general potential function for folding proteins by theoretical simulations.

On the magnitude of the electrostatic contribution to ligand-DNA interactions.

A model based on the nonlinear Poisson-Boltzmann equation is used to study the electrostatic contribution to the binding free energy of a simple intercalating ligand, 3,8-diamino-6-phenylphenanthridine, to DNA. We find that the nonlinear Poisson-Boltzmann model accurately describes both the absolute magnitude of the pKa shift of 3,8-diamino-6-phenylphenanthridine observed upon intercalation and its variation with bulk salt concentration. Since the pKa shift is directly related to the total electrostatic binding free energy of the charged and neutral forms of the ligand, the accuracy of the calculations implies that the electrostatic contributions to binding are accurately predicted as well. Based on our results, we have developed a general physical description of the electrostatic contribution to ligand-DNA binding in which the electrostatic binding free energy is described as a balance between the coulombic attraction of a ligand to DNA and the disruption of solvent upon binding. Long-range coulombic forces associated with highly charged nucleic acids provide a strong driving force for the interaction of cationic ligands with DNA. These favorable electrostatic interactions are, however, largely compensated for by unfavorable changes in the solvation of both the ligand and the DNA upon binding. The formation of a ligand-DNA complex removes both charged and polar groups at the binding interface from pure solvent while it displaces salt from around the nucleic acid. As a result, the total electrostatic binding free energy is quite small. Consequently, nonpolar interactions, such as tight packing and hydrophobic forces, must play a significant role in ligand-DNA stability.

Development of a novel pyrolysis-gas chromatography/mass spectrometry method for the analysis of poly(lactic acid) thermal degradation products

Thermal degradation of PLA is a complex process since it comprises many simultaneous reactions. The use of analytical techniques, such as differential scanning calorimetry (DSC) and thermogravimetry (TGA), yields useful information but a more sensitive analytical technique would be necessary to identify and quantify the PLA degradation products. In this work the thermal degradation of PLA at high temperatures was studied by using a pyrolyzer coupled to a gas chromatograph with mass spectrometry detection (Py-GC/MS). Pyrolysis conditions (temperature and time) were optimized in order to obtain an adequate chromatographic separation of the compounds formed during heating. The best resolution of chromatographic peaks was obtained by pyrolyzing the material from room temperature to 600 °C during 0.5 s. These conditions allowed identifying and quantifying the major compounds produced during the PLA thermal degradation in inert atmosphere. The strategy followed to select these operation parameters was by using sequential pyrolysis based on the adaptation of mathematical models. By application of this strategy it was demonstrated that PLA is degraded at high temperatures by following a non-linear behaviour. The application of logistic and Boltzmann models leads to good fittings to the experimental results, despite the Boltzmann model provided the best approach to calculate the time at which 50% of PLA was degraded. In conclusion, the Boltzmann method can be applied as a tool for simulating the PLA thermal degradation.

Surface growth for molten silicon infiltration into carbon millimeter-sized channels: Lattice–Boltzmann simulations, experiments and models

The process of liquid silicon infiltration is investigated for channels with radii from 0.25 to 0.75 [mm] drilled in compact carbon preforms. The advantage of this setup is that the study of the phenomenon results to be simplified. For comparison purposes, attempts are made in order to work out a framework for evaluating the accuracy of simulations. The approach relies on dimensionless numbers involving the properties of the surface reaction. It turns out that complex hydrodynamic behavior derived from second Newton law can be made consistent with Lattice-Boltzmann simulations. The experiments give clear evidence that the growth of silicon carbide proceeds in two different stages and basic mechanisms are highlighted. Lattice-Boltzmann simulations prove to be an effective tool for the description of the growing phase. Namely, essential experimental constraints can be implemented. As a result, the existing models are useful to gain more insight on the process of reactive infiltration into porous media in the first stage of penetration, i.e. up to pore closure because of surface growth. A way allowing to implement the resistance from chemical reaction in Darcy law is also proposed.

ACcESS: Australia's contribution to the iSERVO Institute's development

The new Australian Computational Earth Systems Simulator research facility provides a virtual laboratory for studying the solid earth and its complex system behavior. The facility's capabilities complement those developed by overseas groups, thereby creating the infrastructure for an international computational solid earth research virtual observatory.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

Dynamic rupture in a 3-D particle-based simulation of a rough planar fault

An appreciation of the physical mechanisms which cause observed seismicity complexity is fundamental to the understanding of the temporal behaviour of faults and single slip events. Numerical simulation of fault slip can provide insights into fault processes by allowing exploration of parameter spaces which influence microscopic and macroscopic physics of processes which may lead towards an answer to those questions. Particle-based models such as the Lattice Solid Model have been used previously for the simulation of stick-slip dynamics of faults, although mainly in two dimensions. Recent increases in the power of computers and the ability to use the power of parallel computer systems have made it possible to extend particle-based fault simulations to three dimensions. In this paper a particle-based numerical model of a rough planar fault embedded between two elastic blocks in three dimensions is presented. A very simple friction law without any rate dependency and no spatial heterogeneity in the intrinsic coefficient of friction is used in the model. To simulate earthquake dynamics the model is sheared in a direction parallel to the fault plane with a constant velocity at the driving edges. Spontaneous slip occurs on the fault when the shear stress is large enough to overcome the frictional forces on the fault. Slip events with a wide range of event sizes are observed. Investigation of the temporal evolution and spatial distribution of slip during each event shows a high degree of variability between the events. In some of the larger events highly complex slip patterns are observed.

Assessment of physical and hydrological properties of biological soil crusts using x-ray microtomography and modelling

Biological soil crusts (BSCs) are formed by aggregates of soil particles and communities of microbial organisms and are common in all drylands. The role of BSCs on infiltration remains uncertain due to the lack of data on their role in affecting soil physical properties such as porosity and structure. Quantitative assessment of these properties is primarily hindered by the fragile nature of the crusts. Here we show how the use of a combination of non-destructive imaging X-ray microtomography (XMT) and Lattice Boltzmann method (LBM) enables quantification of key soil physical parameters and the modeling of water flow through BSCs samples from Kalahari Sands, Botswana. We quantify porosity and flow changes as a result of mechanical disturbance of such a fragile cyanobacteria-dominated crust. Results show significant variations in porosity between different types of crusts and how they affect the flow and that disturbance of a cyanobacteria-dominated crust results in the breakdown of larger pore spaces and reduces flow rates through the surface layer. We conclude that the XMT–LBM approach is well suited for study of fragile surface crust samples where physical and hydraulic properties cannot be easily quantified using conventional methods.

Simulations of charged droplet collisions in shear flow

Acknowledgments This research has been enabled by the use of computing resources provided by WestGrid, the Shared Hierarchical Academic Research Computing Network (SHARCNET: www.sharcnet.ca), and Compute/Calcul Canada. O.S. thanks NSERC for an Alexander Graham Bell Canada Graduate Scholarship.

Finite difference calculations of permeability in large domains in a wide porosity range.

Determining effective hydraulic, thermal, mechanical and electrical properties of porous materials by means of classical physical experiments is often time-consuming and expensive. Thus, accurate numerical calculations of material properties are of increasing interest in geophysical, manufacturing, bio-mechanical and environmental applications, among other fields. Characteristic material properties (e.g. intrinsic permeability, thermal conductivity and elastic moduli) depend on morphological details on the porescale such as shape and size of pores and pore throats or cracks. To obtain reliable predictions of these properties it is necessary to perform numerical analyses of sufficiently large unit cells. Such representative volume elements require optimized numerical simulation techniques. Current state-of-the-art simulation tools to calculate effective permeabilities of porous materials are based on various methods, e.g. lattice Boltzmann, finite volumes or explicit jump Stokes methods. All approaches still have limitations in the maximum size of the simulation domain. In response to these deficits of the well-established methods we propose an efficient and reliable numerical method which allows to calculate intrinsic permeabilities directly from voxel-based data obtained from 3D imaging techniques like X-ray microtomography. We present a modelling framework based on a parallel finite differences solver, allowing the calculation of large domains with relative low computing requirements (i.e. desktop computers). The presented method is validated in a diverse selection of materials, obtaining accurate results for a large range of porosities, wider than the ranges previously reported. Ongoing work includes the estimation of other effective properties of porous media.

Efficient lattice (H)IBE in the standard model

We construct an efficient identity based encryption system based on the standard learning with errors (LWE) problem. Our security proof holds in the standard model. The key step in the construction is a family of lattices for which there are two distinct trapdoors for finding short vectors. One trapdoor enables the real system to generate short vectors in all lattices in the family. The other trapdoor enables the simulator to generate short vectors for all lattices in the family except for one. We extend this basic technique to an adaptively-secure IBE and a Hierarchical IBE.

Lattice-based certificateless public-key encryption in the standard model

The notion of certificateless public-key encryption (CL-PKE) was introduced by Al-Riyami and Paterson in 2003 that avoids the drawbacks of both traditional PKI-based public-key encryption (i.e., establishing public-key infrastructure) and identity-based encryption (i.e., key escrow). So CL-PKE like identity-based encryption is certificate-free, and unlike identity-based encryption is key escrow-free. In this paper, we introduce simple and efficient CCA-secure CL-PKE based on (hierarchical) identity-based encryption. Our construction has both theoretical and practical interests. First, our generic transformation gives a new way of constructing CCA-secure CL-PKE. Second, instantiating our transformation using lattice-based primitives results in a more efficient CCA-secure CL-PKE than its counterpart introduced by Dent in 2008.

Lattice-based completely non-malleable public-key encryption in the standard model

An encryption scheme is non-malleable if giving an encryption of a message to an adversary does not increase its chances of producing an encryption of a related message (under a given public key). Fischlin introduced a stronger notion, known as complete non-malleability, which requires attackers to have negligible advantage, even if they are allowed to transform the public key under which the related message is encrypted. Ventre and Visconti later proposed a comparison-based definition of this security notion, which is more in line with the well-studied definitions proposed by Bellare et al. The authors also provide additional feasibility results by proposing two constructions of completely non-malleable schemes, one in the common reference string model using non-interactive zero-knowledge proofs, and another using interactive encryption schemes. Therefore, the only previously known completely non-malleable (and non-interactive) scheme in the standard model, is quite inefficient as it relies on generic NIZK approach. They left the existence of efficient schemes in the common reference string model as an open problem. Recently, two efficient public-key encryption schemes have been proposed by Libert and Yung, and Barbosa and Farshim, both of them are based on pairing identity-based encryption. At ACISP 2011, Sepahi et al. proposed a method to achieve completely non-malleable encryption in the public-key setting using lattices but there is no security proof for the proposed scheme. In this paper we review the mentioned scheme and provide its security proof in the standard model. Our study shows that Sepahi’s scheme will remain secure even for post-quantum world since there are currently no known quantum algorithms for solving lattice problems that perform significantly better than the best known classical (i.e., non-quantum) algorithms.

Lattice-based completely non-malleable PKE in the standard model

This paper presents ongoing work toward constructing efficient completely non-malleable public-key encryption scheme based on lattices in the standard (common reference string) model. An encryption scheme is completely non-malleable if it requires attackers to have negligible advantage, even if they are allowed to transform the public key under which the related message is encrypted. Ventre and Visconti proposed two inefficient constructions of completely non-malleable schemes, one in the common reference string model using non-interactive zero-knowledge proofs, and another using interactive encryption schemes. Recently, two efficient public-key encryption schemes have been proposed, both of them are based on pairing identity-based encryption.