A Modified Valley Restrained Monte Carlo Method to Efficiently Search the Low Energy Structures of Peptides

Abstract

A new Monte Carlo sampling scheme, namely the Modified Valley Restrained Monte Carlo procedure, is used to obtain the global energy minimum conformations for polypeptides, such as Met-enkephalin and Melittin. For each peptide, we found close agreement with previous results from both theoretical and experimental studies. The simple idea for controlling the step size according to the Valley Function, provides useful suggestions in searching the global energy minimum structures, and furthermore helps solve the multiple minima problem.     

A combined steepest descent and genetic algorithm (SD/GA) approach for the optimization of solvation parameters

An accurate solvation model is essential for computer modeling of protein folding and other biomolecular self-assembly processes. Compared to explicit solvent models, implicit solvent models, such as the Poisson-Boltzmann (PB) with solvent accessible surface area model (PB/SA), offer a much faster speed—the most compelling reason for the popularity of these implicit solvent models. Since these implicit solvent models typically use empirical parameters, such as atomic radii and the surface tensions, an optimal fit of these parameters is crucial for the final accuracy of properties such as solvation free energy and folding free energy. In this paper, we proposed a combined approach, namely SD/GA, which takes the advantage of both local optimization with the steepest descent (SD), and global optimization with the genetic algorithm (GA), for parameters optimization in multi-dimensional space. The SD/GA method is then applied to the optimization of solvation parameters in the non-polar cavity term of the PB/SA model. The results show that the newly optimized parameters from SD/GA not only increase the accuracy in the solvation free energies for ～200 organic molecules, but also significantly improve the free energy landscape of a β-hairpin folding. The current SD/GA method can be readily applied to other multi-dimensional parameter space optimization as well.     

New potential model for molecular dynamic simulation of liquid HF. I—Parameter optimization for charge equilibration method

In order to build a complete potential model to perform classical molecular dynamic simulations of liquid HF, a new optimization method is proposed to obtain transferable parameters for charge equilibration method on the basis of ab initio reference data. The optimized parameters (the electronegativity χ and the Slater orbital exponent ζ for H and F atoms) appear to be able to reproduce the variations of the electrostatic potential calculated from an ab initio method in a liquid phase of HF molecules for different thermodynamic conditions. It is concluded that the proposed method is general, precise and efficient to obtain transferable and realistic parameters.     

BetaDock: Shape-Priority Docking Method Based on Beta-Complex

Abstract

This paper presents an approach and a software, BetaDock, to the docking problem by putting the priority on shape complementarity between a receptor and a ligand. The approach is based on the theory of the β-complex. Given the Voronoi diagram of the receptor whose topology is stored in the quasi-triangulation, the β-complex corresponding to water molecule is computed. Then, the boundary of the β-complex defines the β-shape which has the complete proximity information among all atoms on the receptor boundary. From the β-shape, we first compute pockets where the ligand may bind. Then, we quickly place the ligand within each pocket by solving the singular value decomposition problem and the assignment problem. Using the conformations of the ligands within the pockets as the initial solutions, we run the genetic algorithm to find the optimal solution for the docking problem. The performance of the proposed algorithm was verified through a benchmark test and showed that BetaDock is superior to a popular docking software AutoDock 4.     

Melatonin in the multi‐oscillatory mammalian circadian world

In mammals, the complex interaction of neural, hormonal, and behavioral outputs from the suprachiasmatic nucleus (SCN) drives circadian expression of events, either directly or through coordination of the timing of peripheral oscillators. Melatonin, one of the endocrine output signals of the clock, provides the organism with circadian information and can be considered as an endogenous synchronizer, able to stabilize and reinforce circadian rhythms and to maintain their mutual phase‐relationship at the different levels of the circadian network. Moreover, exogenous melatonin, through an action on the circadian clock, affects all levels of the circadian network. The molecular mechanisms underlying this chronobiotic effect have also been investigated in rats. REV‐ERB α seems to be the initial molecular target.     

High production of cold-tolerant chitinases on shrimp wastes in bench-top bioreactor by the Antarctic fungus Lecanicillium muscarium CCFEE 5003: Bioprocess optimization and characterization of two main enzymes

The Antarctic fungus Lecanicillium muscarium CCFEE-5003 was preliminary cultivated in shaken flasks to check its chitinase production on rough shrimp and crab wastes. Production on shrimp shells was much higher than that on crab shells (104.6 ± 9.3 and 48.6 ± 3.1 U/L, respectively). For possible industrial applications, bioprocess optimization was studied on shrimp shells in bioreactor using RSM to state best conditions of pH and substrate concentration. Optimization improved the production by 137% (243.6 ± 17.3). Two chitinolytic enzymes (CHI1 and CHI2) were purified and characterized. CHI1 (MW ca. 61 kDa) showed optima at pH 5.5 and 45 °C while CHI2 (MW ca. 25 kDa) optima were at pH 4.5 and 40 °C. Both enzymes maintained high activity levels at 5 °C and were inhibited by Fe⁺⁺, Hg⁺⁺ and Cu⁺⁺. CHI2 was strongly allosamidin-sensitive. Both proteins were N-acetyl-hexosaminidases (E.C. 3.2.1.52) but showed different roles in chitin hydrolysis: CHI1 could be defined as “chitobiase” while CHI2 revealed a main “eso-chitinase” activity.     

Optimal stomatal conductance in relation to photosynthesis in climatically contrasting Eucalyptus species under drought

Models of stomatal conductance (g_s) are based on coupling between g_s and CO₂ assimilation (A_net), and it is often assumed that the slope of this relationship (‘g₁’) is constant across species. However, if different plant species have adapted to different access costs of water, then there will be differences in g₁ among species. We hypothesized that g₁ should vary among species adapted to different climates, and tested the theory and its linkage to plant hydraulics using four Eucalyptus species from different climatic origins in a common garden. Optimal stomatal theory predicts that species from sub‐humid zones have a lower marginal water cost of C gain, hence lower g₁ than humid‐zone species. In agreement with the theory that g₁ is related to tissue carbon costs for water supply, we found a relationship between wood density and g₁ across Eucalyptus species of contrasting climatic origins. There were significant reductions in the parameter g₁ during drought in humid but not sub‐humid species, with the latter group maintaining g₁ in drought. There are strong differences in stomatal behaviour among related tree species in agreement with optimal stomatal theory, and these differences are consistent with the economics involved in water uptake and transport for carbon gain.     

Parameters Effects Evaluation of Microbial Strengthening of Sandy Soils in Mixing Experiments Using Taguchi Methodology

The effects of experimental parameters including soil type, curing duration, inoculum size, and biomass and nutrients concentration on soil strengthening due to calcite precipitation by Sporosarcina pasteurii PTCC 1645 were investigated. The laboratory-scale mixing experiments on remolded samples were designed by the Taguchi method. Soil type proved to be the most incorporating factor, followed by curing time and nutrient concentration. The main effect and the interactions of the parameters were presented and the optimal conditions were obtained. This suggests the importance of local conditions including soil type on any future large-scale, in situ application.     

The Trypsin Inhibitor Panulirin Regulates the Prophenoloxidase-activating System in the Spiny Lobster Panulirus argus

The melanization reaction promoted by the prophenoloxidase-activating system is an essential defense response in invertebrates subjected to regulatory mechanisms that are still not fully understood. We report here the finding and characterization of a novel trypsin inhibitor, named panulirin, isolated from the hemocytes of the spiny lobster Panulirus argus with regulatory functions on the melanization cascade. Panulirin is a cationic peptide (pI 9.5) composed of 48 amino acid residues (5.3 kDa), with six cysteine residues forming disulfide bridges. Its primary sequence was determined by combining Edman degradation/N-terminal sequencing and electrospray ionization-MS/MS spectrometry. The low amino acid sequence similarity with known proteins indicates that it represents a new family of peptidase inhibitors. Panulirin is a competitive and reversible tight-binding inhibitor of trypsin (K_i = 8.6 nm) with a notable specificity because it does not inhibit serine peptidases such as subtilisin, elastase, chymotrypsin, thrombin, and plasmin. The removal of panulirin from the lobster hemocyte lysate leads to an increase in phenoloxidase response to LPS. Likewise, the addition of increasing concentrations of panulirin to a lobster hemocyte lysate, previously depleted of trypsin-inhibitory activity, decreased the phenoloxidase response to LPS in a concentration-dependent fashion. These results indicate that panulirin is implicated in the regulation of the melanization cascade in P. argus by inhibiting peptidase(s) in the pathway toward the activation of the prophenoloxidase enzyme.     

New Insights into the Assembly of Bacterial Secretins: STRUCTURAL STUDIES OF THE PERIPLASMIC DOMAIN OF XcpQ FROM PSEUDOMONAS AERUGINOSA*

The type II secretion system is a multiprotein assembly spanning the inner and outer membranes in Gram-negative bacteria. It is found in almost all pathogenic bacteria where it contributes to virulence, host tissue colonization, and infection. The exoproteins are secreted across the outer membrane via a large translocation channel, the secretin, which typically adopts a dodecameric structure. These secretin channels have large periplasmic N-terminal domains that reach out into the periplasm for communication with the inner membrane platform and with a pseudopilus structure that spans the periplasm. Here we report the crystal structure of the N-terminal periplasmic domain of the secretin XcpQ from Pseudomonas aeruginosa, revealing a two-lobe dimeric assembly featuring parallel subunits engaging in well defined interactions at the tips of each lobe. We have employed structure-based engineering of disulfide bridges and native mass spectrometry to show that the periplasmic domain of XcpQ dimerizes in a concentration-dependent manner. Validation of these insights in the context of cellular full-length XcpQ and further evaluation of the functionality of disulfide-linked XcpQ establishes that the basic oligomerization unit of XcpQ is a dimer. This is consistent with the notion that the dodecameric secretin assembles as a hexamer of dimers to ensure correct projection of the N-terminal domains into the periplasm. Therefore, our studies provide a key conceptual advancement in understanding the assembly principles and dynamic function of type II secretion system secretins and challenge recent studies reporting monomers as the basic subunit of the secretin oligomer.