An extended aqueous solvation model based on atom-weighted solvent accessible surface areas: SAWSA v2.0 model

A new method is proposed for calculating aqueous solvation free energy based on atom-weighted solvent accessible surface areas. The method, SAWSA v2.0, gives the aqueous solvation free energy by summing the contributions of component atoms and a correction factor. We applied two different sets of atom typing rules and fitting processes for small organic molecules and proteins, respectively. For small organic molecules, the model classified the atoms in organic molecules into 65 basic types and additionally. For small organic molecules we proposed a correction factor of hydrophobic carbon to account for the aggregation of hydrocarbons and compounds with long hydrophobic aliphatic chains. The contributions for each atom type and correction factor were derived by multivariate regression analysis of 379 neutral molecules and 39 ions with known experimental aqueous solvation free energies. Based on the new atom typing rules, the correlation coefficient (r) for fitting the whole neutral organic molecules is 0.984, and the absolute mean error is 0.40 kcal mol^–1, which is much better than those of the model proposed by Wang et al. and the SAWSA model previously proposed by us. Furthermore, the SAWSA v2.0 model was compared with the simple atom-additive model based on the number of atom types (NA). The calculated results show that for small organic molecules, the predictions from the SAWSA v2.0 model are slightly better than those from the atom-additive model based on NA. However, for macromolecules such as proteins, due to the connection between their molecular conformation and their molecular surface area, the atom-additive model based on the number of atom types has little predictive power. In order to investigate the predictive power of our model, a systematic comparison was performed on seven solvation models including SAWSA v2.0, GB/SA_1, GB/SA_2, PB/SA_1, PB/SA_2, AM1/SM5.2R and SM5.0R. The results showed that for organic molecules the SAWSA v2.0 model is better than the other six solvation models. For proteins, the model classified the atoms into 20 basic types and the predicted aqueous free energies of solvation by PB/SA were used for fitting. The solvation model based on the new parameters was employed to predict the solvation free energies of 38 proteins. The predicted values from our model were in good agreement with those from the PB/SA model and were much better than those given by the other four models developed for proteins.Figure The definition of hydrophobic carbons. Here CA, CB and CD are three carbon atoms; X represents a heteroatom. According to our definition, CB is a hydrophobic carbon, CA is not a hydrophobic carbon because a heteroatom is within four atoms and CD is not a hydrophobic carbon because CD is sp²- hydridized and in a six-member ring.Electronic Supplementary Material Supplementary material is available for this article at     

A new computational model for protein folding based on atomic solvation.

A new model for calculating the solvation energy of proteins is developed and tested for its ability to identify the native conformation as the global energy minimum among a group of thousands of computationally generated compact non-native conformations for a series of globular proteins. In the model (called the WZS model), solvation preferences for a set of 17 chemically derived molecular fragments of the 20 amino acids are learned by a training algorithm based on maximizing the solvation energy difference between native and non-native conformations for a training set of proteins. The performance of the WZS model confirms the success of this learning approach; the WZS model misrecognizes (as more stable than native) only 7 of 8,200 non-native structures. Possible applications of this model to the prediction of protein structure from sequence are discussed.     

Applicability of a continuum solvation model to the octanol water transfer: CFF91 based model for amino acids

A continuum hydration model based upon the atomic charges provided with the CFF91 force field [A. B. Schmidt and R. M. Fine (1994) Molecular Simulation, 13. 347–365] has been extended to the octanol–water transfer. The electrostatic component of the transfer free energy is calculated using the finite-difference solution to the Poisson–Boltzmann equation while the nonpolar contributions are assumed to be proportional to the solute-excluded volume in water. All atomic charges and radii besides the aromatic carbon radius are equal in both solvents. The octanol dielectric constant and the probe radius are the main fitting parameters defining the octanol phase. The model has been tested for 38 organic molecules related to the amino acid residues and generally provides a high accuracy. In particular, the mean unsigned error for N-acetyl amino acid amides is 0.5 kcal/mol. © 1995 John Wiley & Sons, Inc.     

A model of bacterial DNA segregation based upon helical geometry

A new mechanism to segregate daughter genomes in bacterial cells is suggested that is based upon the rules of geometry governing the helix clock (Mendelson, 1982a). The reorientation of cell surface string arrays used as a timing reference in the helix clock is capable of drawing apart the initial products of DNA replication. Physically linking the sister DNA replication origins to the ends of the initial cell surface string inserted into the cell surface at the start of a helix clock cycle, and linking the DNA terminus to a point along the length of the same string provides a means to mark the locations to which the genomes will segregate as well as the place where cell division will occur. The parallel packing of additional cell surface strings into an array which includes the string to which DNA is attached provides the necessary spatial rearrangements. The helical segregation model can account for the precise registration of cell divisions with the completion of replication forks in a multifork replication system, provides a basis for determining the relationship of sister cell sizes at division, and can also accommodate the asymmetrical divisions associated with minicell production and sporulation. Examination of the helical segregation theory under multifork DNA replication conditions moreover reveals that adjacent helical clocks are physically linked to one another although totally independent in terms of their progression through the clock cycle. A relationship between the initiation of DNA replication forks and the insertion of the first cell surface string associated with the start of a helix clock cycle is predicted by the model.     

An articular cartilage contact model based on real surface geometry

Abnormal, excessive stresses acting on articular joint surfaces are speculated to be one of the causes for joint degeneration. However, articular surface stresses have not been studied systematically, since it is technically difficult to measure in vivo contact areas and pressures in dynamic situations. Therefore, we implemented a numerical model of articular surface contact using accurate surface geometries. The model was developed for the cat patellofemoral joint. We demonstrated that small misalignments of the patella relative to the femur change the joint contact mechanics substantially for a given external load. These results suggest that misalignment might be studied as one of the factors causing articular cartilage disorder and joint degeneration.     

ASPDock: protein-protein docking algorithm using atomic solvation parameters model

Background  

Atomic Solvation Parameters (ASP) model has been proven to be a very successful method of calculating the binding free energy of protein complexes. This suggests that incorporating it into docking algorithms should improve the accuracy of prediction. In this paper we propose an FFT-based algorithm to calculate ASP scores of protein complexes and develop an ASP-based protein-protein docking method (ASPDock).     

Efficient electrostatic solvation model for protein-fragment docking

A method is presented for the fast evaluation of the binding energy of a protein-small molecule complex with electrostatic solvation. It makes use of a fast preprocessing step based on the assumption that the main contribution to electrostatic desolvation upon ligand binding originates from the displacement of the first shell of water molecules. For a rigid protein, the precomputation of the energy contributions on a set of grids allows the estimation of the energy in solution of about 300 protein-fragment binding modes per second on a personal computer. The docking procedure is applied to five rigid binding sites whose size ranges from 17 residues to a whole protein of 107 amino acids. Using a library of 70 mainly rigid molecules, known micromolar inhibitors or close analogs are docked and prioritized correctly. The docking based rank-ordering of the library requires about 5 h and is proposed as a complementary approach to structure-activity relationships by nuclear magnetic resonance. Proteins 2001;42:256-268.     

A predictive model of cell traction forces based on cell geometry

Recent work has indicated that the shape and size of a cell can influence how a cell spreads, develops focal adhesions, and exerts forces on the substrate. However, it is unclear how cell shape regulates these events. Here we present a computational model that uses cell shape to predict the magnitude and direction of forces generated by cells. The predicted results are compared to experimentally measured traction forces, and show that the model can predict traction force direction, relative magnitude, and force distribution within the cell using only cell shape as an input. Analysis of the model shows that the magnitude and direction of the traction force at a given point is proportional to the first moment of area about that point in the cell, suggesting that contractile forces within the cell act on the entire cytoskeletal network as a single cohesive unit. Through this model, we demonstrate that intrinsic properties of cell shape can facilitate changes in traction force patterns, independently of heterogeneous mechanical properties or signaling events within the cell.     

Differential growth dynamics control aerial organ geometry

1. Download : Download high-res image (76KB)
2. Download : Download full-size image

     

An augmented Lagrangian finite element formulation for 3D contact of biphasic tissues

Biphasic contact analysis is essential to obtain a complete understanding of soft tissue biomechanics, and the importance of physiological structure on the joint biomechanics has long been recognised; however, up to date, there are no successful developments of biphasic finite element contact analysis for three-dimensional (3D) geometries of physiological joints. The aim of this study was to develop a finite element formulation for biphasic contact of 3D physiological joints. The augmented Lagrangian method was used to enforce the continuity of contact traction and fluid pressure across the contact interface. The biphasic contact method was implemented in the commercial software COMSOL Multiphysics 4.2^® (COMSOL, Inc., Burlington, MA). The accuracy of the implementation was verified using 3D biphasic contact problems, including indentation with a flat-ended indenter and contact of glenohumeral cartilage layers. The ability of the method to model multibody biphasic contact of physiological joints was proved by a 3D knee model. The 3D biphasic finite element contact method developed in this study can be used to study the biphasic behaviours of the physiological joints.     

Landscape as a model: the importance of geometry

In all models, but especially in those used to predict uncertain processes (e.g., climate change and nonnative species establishment), it is important to identify and remove any sources of bias that may confound results. This is critical in models designed to help support decisionmaking. The geometry used to represent virtual landscapes in spatially explicit models is a potential source of bias. The majority of spatial models use regular square geometry, although regular hexagonal landscapes have also been used. However, there are other ways in which space can be represented in spatially explicit models. For the first time, we explicitly compare the range of alternative geometries available to the modeller, and present a mechanism by which uncertainty in the representation of landscapes can be incorporated. We test how geometry can affect cell-to-cell movement across homogeneous virtual landscapes and compare regular geometries with a suite of irregular mosaics. We show that regular geometries have the potential to systematically bias the direction and distance of movement, whereas even individual instances of landscapes with irregular geometry do not. We also examine how geometry can affect the gross representation of real-world landscapes, and again show that individual instances of regular geometries will always create qualitative and quantitative errors. These can be reduced by the use of multiple randomized instances, though this still creates scale-dependent biases. In contrast, virtual landscapes formed using irregular geometries can represent complex real-world landscapes without error. We found that the potential for bias caused by regular geometries can be effectively eliminated by subdividing virtual landscapes using irregular geometry. The use of irregular geometry appears to offer spatial modellers other potential advantages, which are as yet underdeveloped. We recommend their use in all spatially explicit models, but especially for predictive models that are used in decisionmaking.     

Photobioreactors for microalgal cultures: A Lagrangian model coupling hydrodynamics and kinetics

Closed photobioreactors have to be optimized in terms of light utilization and overall photosynthesis rate. A simple model coupling the hydrodynamics and the photosynthesis kinetics has been proposed to analyze the photosynthesis dynamics due to the continuous shuttle of microalgae between dark and lighted zones of the photobioreactor. Microalgal motion has been described according to a stochastic Lagrangian approach adopting the turbulence model suitable for the photobioreactor configuration (single vs. two‐phase flows). Effects of light path, biomass concentration, turbulence level and irradiance have been reported in terms of overall photosynthesis rate. Different irradiation strategies (internal, lateral and rounding) and several photobioreactor configurations (flat, tubular, bubble column, airlift) have been investigated. Photobioreactor configurations and the operating conditions to maximize the photosynthesis rate have been pointed out. Results confirmed and explained the common experimental observation that high concentrated cultures are not photoinhibited at high irradiance level. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:1259–1272, 2015     

Differential interaction of equinatoxin II with model membranes in response to lipid composition

Equinatoxin II is a 179-amino-acid pore-forming protein isolated from the venom of the sea anemone Actinia equina. Large unilamellar vesicles and lipid monolayers of different lipid compositions have been used to study its interaction with membranes. The critical pressure for insertion is the same in monolayers made of phosphatidylcholine or sphingomyelin (approximately 26 mN m(-1)) and explains why the permeabilization of large unilamellar vesicles by equinatoxin II with these lipid compositions is null or moderate. In phosphatidylcholine-sphingomyelin (1:1) monolayers, the critical pressure is higher (approximately 33 mN m(-1)), thus permitting the insertion of equinatoxin II in large unilamellar vesicles, a process that is accompanied by major conformational changes. In the presence of vesicles made of phosphatidylcholine, a fraction of the protein molecules remains associated with the membranes. This interaction is fully reversible, does not involve major conformational changes, and is governed by the high affinity for membrane interfaces of the protein region comprising amino acids 101-120. We conclude that although the presence of sphingomyelin within the membrane creates conditions for irreversible insertion and pore formation, this lipid is not essential for the initial partitioning event, and its role as a specific receptor for the toxin is not so clear-cut.     

A mathematical model of tumor growth. II. effects of geometry and spatial nonuniformity on stability

A theoretical account of mitotic inhibition in one-, two-, and three-dimensional configurations is presented. Based on an earlier paper, the inhibitor production rate is taken to be nonuniform throughout the tissue, resulting in significant deviations from the prediction of uniform models. Geometry affects the stability of growth also. The analysis presented here represents a detailed study of the properties of highly nonuniform inhibition, from which information on intermediate inhibition models can be readily deduced. This information is used to compare such a model with experimental results in the third paper in this series.     

Novel flexible vesicles based topical formulation of levocetirizine: in vivo evaluation using oxazolone-induced atopic dermatitis in murine model

Background: Levocetirizine, an active enantiomer of cetirizine is third generation antihistaminic agent used for treating various allergies like atopic dermatitis, chronic idiopathic urticaria and allergic rhinitis.

Objective: Development of novel topical formulation of levocetirizine based on flexible vesicles (FVs) with an aim to have targeted peripheral antihistaminic effect.

Materials and methods: The FVs were prepared by thin film hydration method and characterized for drug content, entrapment efficiency, pH, vesicular size, spreadability, morphological characteristics and drug leakage studies. Franz diffusion cell assembly was used to carry out the ex vivo permeation studies through mice skin and the permeation profile of the developed FV formulation was compared with conventional formulations of levocetirizine.

Results and discussion: The ex vivo permeation studies revealed 1.78-fold increase in percent permeation of levocetirizine from FV formulation as compared to conventional formulations of levocetirizine in 8?h. Further, oxazolone induced atopic dermatitis murine model was selected to study the in vivo pharmacodynamic activity. The developed formulation was evaluated for scratching score, erythema score and histological evaluation. There was marked reduction in scratching score from 15.25 scratches/20?min with conventional levocetirizine cream to 6.75 scratches/20?min with application of levocetirizine FV formulation. Also, there was significant reduction in erythema score as well as dermal eosinophil count. Results of skin sensitivity and toxicity studies suggest that the developed formulation was dermally safe and nontoxic.

Conclusion: A novel FVs based topical formulation of levocetirizine was successfully developed for treatment of atopic dermatitis.     

Modeling individual effects in the Cormack-Jolly-Seber model: a state-space formulation

Summary .   In population and evolutionary biology, there exists considerable interest in individual heterogeneity in parameters of demographic models for open populations. However, flexible and practical solutions to the development of such models have proven to be elusive. In this article, I provide a state-space formulation of open population capture–recapture models with individual effects. The state-space formulation provides a generic and flexible framework for modeling and inference in models with individual effects, and it yields a practical means of estimation in these complex problems via contemporary methods of Markov chain Monte Carlo. A straightforward implementation can be achieved in the software package WinBUGS . I provide an analysis of a simple model with constant parameter detection and survival probability parameters. A second example is based on data from a 7-year study of European dippers, in which a model with year and individual effects is fitted.