Empirical prediction of peptide octanol-water partition coefficients

Peptides are of great therapeutic potential as vaccines and drugs. Knowledge of physicochemical descriptors, including the partition coefficient P (commonly expressed in logarithm form: logP), is useful for screening out unsuitable molecules and also for the development of predictive Quantitative Structure-Activity Relationships (QSARs). In this paper we develop a new approach to the prediction of LogP values for peptides based on an empirical relationship between global molecular properties and measured physical properties. Our method was successful in terms of peptide prediction (total r2 = 0.641). The final model consisted of 5 physicochemical descriptors (molecular weight, number of single bonds, 2D-VDW volume, 2D-VSA hydrophobic and 2D-VSA polar). The approach is peptide specific and its predictive accuracy was high. Overall, 67% of the peptides were able to be predicted within +/-0.5 log units from the experimental values. Our method thus represents a novel prediction method with proven predictive ability.

Monte Carlo and molecular dynamics simulations of methane in potassium montmorillonite clay hydrates at elevated pressures and temperatures

The structure and dynamics of methane in hydrated potassium montmorillonite clay have been studied under conditions encountered in sedimentary basin and compared to those of hydrated sodium montmorillonite clay using computer simulation techniques. The simulated systems contain two molecular layers of water and followed gradients of 150 barkm-1 and 30 Kkm-1 up to a maximum burial depth of 6 km. Methane particle is coordinated to about 19 oxygen atoms, with 6 of these coming from the clay surface oxygen. Potassium ions tend to move away from the center towards the clay surface, in contrast to the behavior observed with the hydrated sodium form. The clay surface affinity for methane was found to be higher in the hydrated K-form. Methane diffusion in the two-layer hydrated K-montmorillonite increases from 0.39×10-9 m2s-1 at 280 K to 3.27×10-9 m2s-1 at 460 K compared to 0.36×10-9 m2s-1 at 280 K to 4.26×10-9 m2s-1 at 460 K in Na-montmorillonite hydrate. The distributions of the potassium ions were found to vary in the hydrates when compared to those of sodium form. Water molecules were also found to be very mobile in the potassium clay hydrates compared to sodium clay hydrates. © 2004 Elsevier Inc. All All rights reserved.

The role of lipid geometry in designing liposomes for the solubilisation of poorly water soluble drugs

Liposomes are well recognised for their ability to improve the delivery of a range of drugs. More commonly they are applied for the delivery of water-soluble drugs, but given their structural attributes, they can also be employed as solubilising agents for low solubility drugs as well as drug targeting agents. To further explore the potential of liposomes as solubilising agents, we have investigated the role of bilayer packaging in promoting drug solubilisation in liposome bilayers. The effect of alkyl chain length and symmetry was investigated to consider if using 'mis-matched' phospholipids could create 'voids' within the bilayers, and enhance bilayer loading capacity. Lipid packing was investigated using Langmuir studies, which demonstrated that increasing the alkyl chain length enhanced lipid packing, with condensed monolayers forming, whilst asymmetric lipids formed less condensed monolayers. However, this more open packing did not translate into improved drug loading, with the longer chain, condensed bilayers formed from long-chain, saturated lipids offering higher drug loading capacity. These studies demonstrate that liposomes formulated from longer chain, saturated lipids offer enhanced solubilisation capacity. However the molecular size, rather than lipophilicity, of the drug to be incorporated was also a key factor dominating bilayer incorporation efficiency. © 2012 Elsevier B.V. All rights reserved.

Dynamic complexity of chaotic transitions in high-dimensional classical dynamics:Leu-Enkephalin folding

Leu-Enkephalin in explicit water is simulated using classical molecular dynamics. A ß-turn transition is investigated by calculating the topological complexity (in the "computational mechanics" framework [J. P. Crutchfield and K. Young, Phys. Rev. Lett., 63, 105 (1989)]) of the dynamics of both the peptide and the neighbouring water molecules. The complexity of the atomic trajectories of the (relatively short) simulations used in this study reflect the degree of phase space mixing in the system. It is demonstrated that the dynamic complexity of the hydrogen atoms of the peptide and almost all of the hydrogens of the neighbouring waters exhibit a minimum precisely at the moment of the ß-turn transition. This indicates the appearance of simplified periodic patterns in the atomic motion, which could correspond to high-dimensional tori in the phase space. It is hypothesized that this behaviour is the manifestation of the effect described in the approach to molecular transitions by Komatsuzaki and Berry [T. Komatsuzaki and R.S. Berry, Adv. Chem. Phys., 123, 79 (2002)], where a "quasi-regular" dynamics at the transition is suggested. Therefore, for the first time, the less chaotic character of the folding transition in a realistic molecular system is demonstrated. © Springer-Verlag Berlin Heidelberg 2006.

Complexity of classical dynamics of molecular systems. I. Methodology

Methods for the calculation of complexity have been investigated as a possible alternative for the analysis of the dynamics of molecular systems. “Computational mechanics” is the approach chosen to describe emergent behavior in molecular systems that evolve in time. A novel algorithm has been developed for symbolization of a continuous physical trajectory of a dynamic system. A method for calculating statistical complexity has been implemented and tested on representative systems. It is shown that the computational mechanics approach is suitable for analyzing the dynamic complexity of molecular systems and offers new insight into the process.

Molecular dynamics simulation of methane in sodium montmorillonite clay hydrates at elevated pressures and temperatures

Computer simulation has been used to study the structure and dynamics of methane in hydrated sodium montmorillonite clays under conditions encountered in sedimentary basins. Systems containing approximately one, two, three and four molecular layers of water have followed gradients of 150 bar km-1 and 30Kkm-1, to a maximum burial depth of 6 km (900 bar and 460 K). Methane is coordinated to approximately 19 oxygen atoms, of which typically 6 are provided by the clay surface. Only in the three- and four-layer hydrates is methane able to leave the clay surface. Diffusion depends strongly on the porosity (water content) and burial depth: self-diffusion coefficients are in the range 0.12 × 10-9m2s-1 for water and 0.04 × 10−9m2s−1 < D < 8.64 × 10−9m2s−1 for methane. Bearing in mind that porosity decreases with burial depth, it is estimated that maximum diffusion occurs at around 3 km. This is in good agreement with the known location of methane reservoirs in sedimentary basins.

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

A hybrid molecular dynamics/fluctuating hydrodynamics method for modelling liquids at multiple scales in space and time

A new 3D implementation of a hybrid model based on the analogy with two-phase hydrodynamics has been developed for the simulation of liquids at microscale. The idea of the method is to smoothly combine the atomistic description in the molecular dynamics zone with the Landau-Lifshitz fluctuating hydrodynamics representation in the rest of the system in the framework of macroscopic conservation laws through the use of a single "zoom-in" user-defined function s that has the meaning of a partial concentration in the two-phase analogy model. In comparison with our previous works, the implementation has been extended to full 3D simulations for a range of atomistic models in GROMACS from argon to water in equilibrium conditions with a constant or a spatially variable function s. Preliminary results of simulating the diffusion of a small peptide in water are also reported.

Light scattering study of human serum albumin in pre-denaturation:relation to dynamic transition in water at 42°C

Protein functional motions are ultimately connected to water dynamics. The goal of this study is to link the conformational dynamics of albumin to a dynamic transition taking place at ∼ 42°C in water. We report the results of dynamic light scattering measurements of albumin aqueous solution in the temperature interval 20-65°C. The processing of the experimental data produced the temperature dependence of the macromolecular hydrodynamic radius. We demonstrate that the growth of the macromolecular size in this temperature range can be divided into two stages that are connected to the dynamical properties of water. © 2012 Elsevier B.V. All rights reserved.

Exploratory catalyst screening studies on the base free conversion of glycerol to lactic acid and glyceric acid in water using bimetallic Au-Pt nanoparticles on acidic zeolites

The base free oxidation of glycerol with molecular oxygen in water using bimetallic Au-Pt catalysts on three different acidic zeolite supports (H-mordenite, H-β and H-USY) was explored in a batch setup. At temperatures between 140 and 180 °C, lactic acid formation was significant and highest selectivity (60 % lactic acid at 80 % glycerol conversion) was obtained using Au-Pt/USY-600 (180 °C). A selectivity switch to glyceric acid (GLYA) was observed when the reactions were performed at 100 °C. Highest conversion and selectivity towards GLYA were obtained with Au-Pt/H-β as the catalyst (68 % selectivity at 68 % conversion).

Organic matter dynamics in the Florida coastal Everglades: A molecular marker and isotopic approach

Everglades National Park (ENP) is about to undergo the world's largest wetland restoration with the aim of improving the quality, timing and distribution of water flow. The changes in water flow are hypothesized to alter the nutrient fluxes and organic matter (OM) dynamics within ENP, especially in the estuarine areas. This study used a multi-proxy approach of molecular markers and stable δ 13C isotope measurements, to determine the present day OM dynamics in ENP. ^ OM dynamics in wetland soils/sediments have proved to be difficult to understand using traditional geochemical approaches. These are often inadequate to describe the multitude of OM sources (e.g. higher land plant, emergent vegetation, submerged vegetation) to the soils/sediments and the complex diagenetic processes that can alter the OM characteristics. A multi-proxy approach, however, that incorporates both molecular level and bulk parameter information is ideal to comprehend complex OM dynamics in aquatic environments. Therefore, biomass-specific molecular markers or proxies can be useful in tracing the sources and processing of OM. This approach was used to examine the OM dynamics in the two major drainage basins, Shark River Slough and Taylor River Slough, of ENP. Freshwater to marine transects were sampled in both systems for soils/sediments and suspended particulate organic matter (SPOM) to be characterized through bulk OM analyses, lipid biomarker determinations (e.g. sterols, fatty acids, hydrocarbons and triterpenoids) and compound-specific stable carbon isotope (δ 13C) determinations. ^ One key accomplishment of the research was the assessment of a molecular marker proxy (Paq) to distinguish between emergent/higher plant vegetation from submerged vegetation within ENP. This proxy proved to be quite useful at tracing OM inputs to the soils/sediments of ENP. A second key accomplishment was the development of a 3-way model using vegetation specific molecular markers. This novel, descriptive model was successfully applied to the estuarine areas of Taylor and Shark River sloughs, providing clear evidence of mixing of freshwater, estuarine and marine derived OM in these areas. In addition, diagenetic transformations of OM in these estuaries were found to be quite different between Taylor and Shark Rivers, and are likely a result of OM quality and hydrological differences. ^

The thermodynamics and statistical mechanics of protein folding

The physics of self-organization and complexity is manifested on a variety of biological scales, from large ecosystems to the molecular level. Protein molecules exhibit characteristics of complex systems in terms of their structure, dynamics, and function. Proteins have the extraordinary ability to fold to a specific functional three-dimensional shape, starting from a random coil, in a biologically relevant time. How they accomplish this is one of the secrets of life. In this work, theoretical research into understanding this remarkable behavior is discussed. Thermodynamic and statistical mechanical tools are used in order to investigate the protein folding dynamics and stability. Theoretical analyses of the results from computer simulation of the dynamics of a four-helix bundle show that the excluded volume entropic effects are very important in protein dynamics and crucial for protein stability. The dramatic effects of changing the size of sidechains imply that a strategic placement of amino acid residues with a particular size may be an important consideration in protein engineering. Another investigation deals with modeling protein structural transitions as a phase transition. Using finite size scaling theory, the nature of unfolding transition of a four-helix bundle protein was investigated and critical exponents for the transition were calculated for various hydrophobic strengths in the core. It is found that the order of the transition changes from first to higher order as the strength of the hydrophobic interaction in the core region is significantly increased. Finally, a detailed kinetic and thermodynamic analysis was carried out in a model two-helix bundle. The connection between the structural free-energy landscape and folding kinetics was quantified. I show how simple protein engineering, by changing the hydropathy of a small number of amino acids, can enhance protein folding by significantly changing the free energy landscape so that kinetic traps are removed. The results have general applicability in protein engineering as well as understanding the underlying physical mechanisms of protein folding. ^

Theoretical study of chloroperoxidase catalyzed chlorination of beta-cyclopentanedione and role of water in the chlorination mechanism

Chloroperoxidase (CPO) is a potential biocatalyst for use in asymmetric synthesis. The mechanisms of CPO catalysis are therefore of interest. The halogenation reaction, one of several chemical reactions that CPO catalyzes, is not fully understood and is the subject of this dissertation. The mechanism by which CPO catalyzes halogenation is disputed. It has been postulated that halogenation of substrates occurs at the active site. Alternatively, it has been proposed that hypochlorous acid, produced at the active site via oxidation of chloride, is released prior to reaction, so that halogenation occurs in solution. The free-solution mechanism is supported by the observation that halogenation of most substrates often occurs non-stereospecifically. On the other hand, the enzyme-bound mechanism is supported by the observation that some large substrates undergo halogenation stereospecifically. The major purpose of this research is to compare chlorination of the substrate β-cyclopentanedione in the two environments. One study was of the reaction with limited hydration because such a level of hydration is typical of the active site. For this work, a purely quantum mechanical approach was used. To model the aqueous environment, the limited hydration environment approach is not appropriate. Instead, reaction precursor conformations were obtained from a solvated molecular dynamics simulation, and reaction of potentially reactive molecular encounters was modeled with a hybrid quantum mechanical/molecular mechanical approach. Extensive work developing parameters for small molecules was pre-requisite for the molecular dynamics simulation. It is observed that a limited and optimized (active-site-like) hydration environment leads to a lower energetic barrier than the fully solvated model representative of the aqueous environment at room temperature, suggesting that the stable water network near the active site is likely to facilitate the chlorination mechanism. The influence of the solvent environment on the reaction barrier is critical. It is observed that stabilization of the catalytic water by other solvent molecules lowers the barrier for keto-enol tautomerization. Placement of water molecules is more important than the number of water molecules in such studies. The fully-solvated model demonstrates that reaction proceeds when the instantaneous dynamical water environment is close to optimal for stabilizing the transition state.

Using soil profiles of seeds and molecular markers as proxies for sawgrass and wet prairie slough vegetation in Shark Slough, Everglades National Park

We measured the abundance of Cladium jamaicense (Crantz) seeds and three biomarkers in freshwater marsh soils in Shark River Slough (SRS), Everglades National Park (ENP) to determine the degree to which these paleoecological proxies reflect spatial and temporal variation in vegetation. We found that C. jamaicense seeds and the biomarkers Paq, total lignin phenols (TLP) and kaurenes analyzed from surface soils were all significantly correlated with extant aboveground C. jamaicense biomass quantified along a vegetation gradient from a C. jamaicense to a wet prairie/slough (WPS) community. Our results also suggest that these individual proxies may reflect vegetation over different spatial scales: Paq and kaurenes correlated most strongly (R 2 = 0.88 and 0.99, respectively) with vegetation within 1 m of a soil sample, while seeds and TLP reflected vegetation 0–20 m upstream of soil samples. These differences in the spatial scale depicted by the different proxies may be complementary in understanding aspects of historic landscape patterning. Soil profiles of short (25 cm) cores showed that downcore variation in C. jamaicense seeds was highly correlated with two of the three biomarkers (Paq, R 2 = 0.84, p<0.005; TLP, R 2 = 0.97, p<0.0001), and all four of the proxies indicated a recent increase in C. jamaicense biomass at the site. Using a preliminary depth-to-age relationship based on matching charcoal peaks with available ENP fire records (1980-present) specific to our coring site, we found that peak-depths in C. jamaicense seed concentration appeared to correspond to recent minimum water levels (e.g., 1989 and 2001), and low seed abundance corresponded to high water levels (e.g., 1995), consistent with the known autecology of C. jamaicense. In summary, the combination of C. jamaicense seeds and biomarkers may be useful for paleoecological reconstruction of vegetation change and ultimately in guaging the success of ongoing efforts to restore historic hydrologic conditions in the South Florida Everglades.

Functional genomic and molecular analyses of transcripts related to glycerol synthesis and trafficking in rainbow smelt (Osmerus mordax) using cold temperature-induced models of glycerol production

The rainbow smelt (Osmerus mordax) is an anadromous teleost that produces type II antifreeze protein (AFP) and accumulates modest urea and high glycerol levels in plasma and tissues as adaptive cryoprotectant mechanisms in sub-zero temperatures. It is known that glyceroneogenesis occurs in liver via a branch in glycolysis and gluconeogenesis and is activated by low temperature; however, the precise mechanisms of glycerol synthesis and trafficking in smelt remain to be elucidated. The objective of this thesis was to provide further insight using functional genomic techniques [e.g. suppression subtractive hybridization (SSH) cDNA library construction, microarray analyses] and molecular analyses [e.g. cloning, quantitative reverse transcription - polymerase chain reaction (QPCR)]. Novel molecular mechanisms related to glyceroneogenesis were deciphered by comparing the transcript expression profiles of glycerol (cold temperature) and non-glycerol (warm temperature) accumulating hepatocytes (Chapter 2) and livers from intact smelt (Chapter 3). Briefly, glycerol synthesis can be initiated from both amino acids and carbohydrate; however carbohydrate appears to be the preferred source when it is readily available. In glycerol accumulating hepatocytes, levels of the hepatic glucose transporter (GLUT2) plummeted and transcript levels of a suite of genes (PEPCK, MDH2, AAT2, GDH and AQP9) associated with the mobilization of amino acids to fuel glycerol synthesis were all transiently higher. In contrast, in glycerol accumulating livers from intact smelt, glycerol synthesis was primarily fuelled by glycogen degradation with higher PGM and PFK (glycolysis) transcript levels. Whether initiated from amino acids or carbohydrate, there were common metabolic underpinnings. Increased PDK2 (an inhibitor of PDH) transcript levels would direct pyruvate derived from amino acids and / or DHAP derived from G6P to glycerol as opposed to oxidation via the citric acid cycle. Robust LIPL (triglyceride catabolism) transcript levels would provide free fatty acids that could be oxidized to fuel ATP synthesis. Increased cGPDH (glyceroneogenesis) transcript levels were not required for increased glycerol production, suggesting that regulation is more likely by post-translational modification. Finally, levels of a transcript potentially encoding glycerol-3-phosphatase, an enzyme not yet characterized in any vertebrate species, were transiently higher. These comparisons also led to the novel discoveries that increased G6Pase (glucose synthesis) and increased GS (glutamine synthesis) transcript levels were part of the low temperature response in smelt. Glucose may provide increased colligative protection against freezing; whereas glutamine could serve to store nitrogen released from amino acid catabolism in a non-toxic form and / or be used to synthesize urea via purine synthesis-uricolysis. Novel key aspects of cryoprotectant osmolyte (glycerol and urea) trafficking were elucidated by cloning and characterizing three aquaglyceroporin (GLP)-encoding genes from smelt at the gene and cDNA levels in Chapter 4. GLPs are integral membrane proteins that facilitate passive movement of water, glycerol and urea across cellular membranes. The highlight was the discovery that AQP10ba transcript levels always increase in posterior kidney only at low temperature. This AQP10b gene paralogue may have evolved to aid in the reabsorption of urea from the proximal tubule. This research has contributed significantly to a general understanding of the cold adaptation response in smelt, and more specifically to the development of a working scenario for the mechanisms involved in glycerol synthesis and trafficking in this species.