A molecular dynamics investigation of mono and dimeric states of the outer membrane enzyme OMPLA

OMPLA is a phospholipase found in the outer membranes of many Gram-negative bacteria. Enzyme activation requires calcium-induced dimerisation plus bilayer perturbation. As the conformation of OMPLA in the different crystal forms (monomer versus dimer; with/without bound Ca(2+)) is remarkably similar we have used multi-nanosecond molecular dynamics (MD) simulations to probe possible differences in conformational dynamics that may be related to enzyme activation. Simulations of calcium-free monomeric OMPLA, of the Ca(2+)-bound dimer, and of the Ca(2+)-bound dimer with a substrate analogue covalently linked to the active site serine have been performed, all with the protein embedded in a phospholipid (POPC) bilayer. All simulations were stable, but differences in the dynamic behaviour of the protein between the various states were observed. In particular, the stability of the active site and the hydrophobic substrate-binding cleft varied. Dimeric OMPLA is less flexible than monomeric OMPLA, especially around the active site. In the absence of bound substrate analogue, the hydrophobic substrate-binding cleft of dimeric OMPLA collapses. A model is proposed whereby the increased stability of the active site in dimeric OMPLA is a consequence of the local ordering of water around the nearby calcium ion. The observed collapse of the substrate-binding cleft may explain the experimentally observed occurrence of multiple dimer conformations of OMPLA, one of which is fully active while the other shows significantly reduced activity.     

Immunogenicity of peptide-vaccine candidates predicted by molecular dynamics simulations

We present an in silico, structure-based approach for design and evaluation of conformationally restricted peptide-vaccines. In particular, we designed four cyclic peptides of ten or 11 residues mimicking the crystallographically observed beta-turn conformation of a predicted immunodominant loop of PorA from Neisseria meningitidis. Conformational correctness and stability of the peptide designs, as evaluated by molecular dynamics simulations, correctly predicted the immunogenicity of the peptides. We observed a peptide-induced functional antibody response that, remarkably, exceeded the response induced by the native protein in outer membrane vesicles, without losing specificity for related strains. The presented approach offers tools for a priori design and selection of peptide-vaccine candidates with full biological activity. This approach could be widely applicable: to outer membrane proteins of Gram-negative bacteria, and to other epitopes in a large range of pathogens.     

The tertiary structure and backbone dynamics of human prolactin

Human prolactin is a 199-residue (23 kDa) protein closely related to growth hormone and placental lactogen with properties and functions resembling both a hormone and a cytokine. As a traditional hormone, prolactin is produced by lactotrophic cells in the pituitary and secreted into the bloodstream where it acts distally to regulate reproduction and promote lactation. Pituitary cells store prolactin in secretory granules organized around large prolactin aggregates, which are produced within the trans layer of the Golgi complex. Extrapituitary prolactin is synthesized by a wide variety of cells but is not stored in secretory granules. Extrapituitary prolactin displays immunomodulatory activities and acts as a growth factor for cancers of the breast, prostate and tissues of the female reproductive system. We have determined the tertiary structure of human prolactin using three-dimensional (3D) and four-dimensional (4D) heteronuclear NMR spectroscopy. As expected, prolactin adopts an "up-up-down-down" four-helical bundle topology and resembles other members of the family of hematopoietic cytokines. Prolactin displays three discrete structural differences from growth hormone: (1) a structured N-terminal loop in contact with the first helix, (2) a missing mini-helix in the loop between the first and second helices, and (3) a shorter loop between the second and third helices lacking the perpendicular mini-helix observed in growth hormone. Residues necessary for functional binding to the prolactin receptor are clustered on the prolactin surface in a position similar to growth hormone. The backbone dynamics of prolactin were investigated by analysis of 15N NMR relaxation phenomena and demonstrated a rigid four-helical bundle with relatively mobile interconnecting loops. Comparison of global macromolecular tumbling at 0.1mM and 1.0mM prolactin revealed reversible oligomerization, which was correlated to dynamic light scattering experiments. The existence of a reversible oligomerization reaction in solution provides insight into previous results describing the in vitro and in vivo aggregation properties of human prolactin.     

Optimal conservation effort for a population in a stochastic environment

We study the optimal conservation effort for a population in a fluctuating environment. The survivorship of a population is affected by unpredictable environmental fluctuation (noise) and can be improved by conservation effort accompanied by a cost. The optimal effort level is the one that minimizes the total cost, defined as the weighted sum of the population extinction risk and the economic cost of conservation effort. The optimal effort depends on the variance and the probability distribution of the noise, the relative importance of the population's survival vs. the economic cost, the effectiveness of conservation effort, and the time scope over which we optimize. The analysis of dynamic programming illustrates that the choice of extinction risk function greatly affects the optimal effort level. The conservation effort level that is the best solution of a multiple-year optimization may be higher than that for the corresponding single-year optimization, if the population is relatively safe. However, the conservation level for the multiple-year optimization becomes lower than for the single-year optimization if the population is endangered. In a similar manner, the optimal conservation effort level for the problem with a short time scope is either higher or lower than that for the problem with a long time scope, depending on the extinction risk of the population. Next, for each parameter of the model, we define five different sensitivities of extinction probability or of the total cost. We then study the mean increase in the total cost caused by the uncertainty of parameters. To achieve the best conservation result, we need to invest the limited research effort to the parameter with the largest effect to the optimal effort level, rather than to those with large impacts on the extinction probability or on the total cost. The recommended policy should depend critically on the choice of the criterion to optimize, which shows the importance of theoretical study of the relationship in performing proper decision making in conservation practice.     

Neurophysical theory of coherence and correlations of electroencephalographic and electrocorticographic signals

Correlation and coherence functions of electroencephalographic signals are calculated using a recent continuum theory that has previously yielded excellent agreement with observations of electroencephalographic spectra. The predicted properties of these functions are found to be in semiquantitative agreement with observations for parameters consistent with those used in previous studies of spectra. The corresponding results for electrocorticographic signals point to an additional contribution at characteristic scales of around 6mm, as has been previously inferred, and are consistent with a crossover to the long-range behavior seen in scalp data. Analysis within the framework of the model enables constraints on the relative strengths of the two contributions to be inferred, and makes it plausible that the short-range component reflects the point-spread function of external stimuli and corticothalamic feedback reaching the cortex.     

Sensitivity analysis of stoichiometric networks: an extension of metabolic control analysis to non-steady state trajectories

A sensitivity analysis of general stoichiometric networks is considered. The results are presented as a generalization of Metabolic Control Analysis, which has been concerned primarily with system sensitivities at steady state. An expression for time-varying sensitivity coefficients is given and the Summation and Connectivity Theorems are generalized. The results are compared to previous treatments. The analysis is accompanied by a discussion of the computation of the sensitivity coefficients and an application to a model of phototransduction.     

The inner membrane protein Mdm33 controls mitochondrial morphology in yeast

Mitochondrial distribution and morphology depend on MDM33, a Saccharomyces cerevisiae gene encoding a novel protein of the mitochondrial inner membrane. Cells lacking Mdm33 contain ring-shaped, mostly interconnected mitochondria, which are able to form large hollow spheres. On the ultrastructural level, these aberrant organelles display extremely elongated stretches of outer and inner membranes enclosing a very narrow matrix space. Dilated parts of Delta mdm33 mitochondria contain well-developed cristae. Overexpression of Mdm33 leads to growth arrest, aggregation of mitochondria, and generation of aberrant inner membrane structures, including septa, inner membrane fragments, and loss of inner membrane cristae. The MDM33 gene is required for the formation of net-like mitochondria in mutants lacking components of the outer membrane fission machinery, and mitochondrial fusion is required for the formation of extended ring-like mitochondria in cells lacking the MDM33 gene. The Mdm33 protein assembles into an oligomeric complex in the inner membrane where it performs homotypic protein-protein interactions. Our results indicate that Mdm33 plays a distinct role in the mitochondrial inner membrane to control mitochondrial morphology. We propose that Mdm33 is involved in fission of the mitochondrial inner membrane.     

RhoA/ROCK regulation of neuritogenesis via profilin IIa-mediated control of actin stability

Neuritogenesis, the first step of neuronal differentiation, takes place as nascent neurites bud from the immediate postmitotic neuronal soma. Little is known about the mechanisms underlying the dramatic morphological changes that characterize this event. Here, we show that RhoA activity plays a decisive role during neuritogenesis of cultured hippocampal neurons by recruiting and activating its specific kinase ROCK, which, in turn, complexes with profilin IIa. We establish that this previously uncharacterized brain-specific actin-binding protein controls neurite sprouting by modifying actin stability, a function regulated by ROCK-mediated phosphorylation. Furthermore, we determine that this novel cascade is switched on or off by physiological stimuli. We propose that RhoA/ROCK/PIIa-mediated regulation of actin stability, shown to be essential for neuritogenesis, may constitute a central mechanism throughout neuronal differentiation.     

Fighting a virus with a virus: a dynamic model for HIV-1 therapy

A mathematical model examined a potential therapy for controlling viral infections using genetically modified viruses. The control of the infection is an indirect effect of the selective elimination by an engineered virus of infected cells that are the source of the pathogens. Therefore, this engineered virus could greatly compensate for a dysfunctional immune system compromised by AIDS. In vitro studies using engineered viruses have been shown to decrease the HIV-1 load about 1000-fold. However, the efficacy of this potential treatment for reducing the viral load in AIDS patients is unknown. The present model studied the interactions among the HIV-1 virus, its main host cell (activated CD4+ T cells), and a therapeutic engineered virus in an in vivo context; and it examined the conditions for controlling the pathogen. This model predicted a significant drop in the HIV-1 load, but the treatment does not eradicate HIV. A basic estimation using a currently engineered virus indicated an HIV-1 load reduction of 92% and a recovery of host cells to 17% of their normal level. Greater success (98% HIV reduction, 44% host cells recovery) is expected as more competent engineered viruses are designed. These results suggest that therapy using viruses could be an alternative to extend the survival of AIDS patients.     

Rapid Axillary Bud Proliferation and ex vitro Rooting of Eupatorium triplinerve

Effective protocol was established for micropropagation of the medicinal plant Eupatorium triplinerve Vahl through rapid axillary bud proliferation and ex vitro rooting. Murashige and Skoog (MS) medium fortified with 8.87 M benzylaminopurine (BAP) and 2.46 M indole-3-butyric acid (IBA) was the best for axillary bud proliferation and developed a mean of 8.1 shoots per node. Excision and culture of the node segments of the in vitro shoots on medium supplemented with the same concentration of growth regulators developed more than 30 shoots within 40 d. Shoot multiplication did not exhibit decrease in the number of shoots even at 7^th subculture. Dipping of the basal end of shoots in 2.46 M IBA solution for 10 d induced roots and its transfer to small pots facilitated the survival of all rooted shoots (100 %). Ex vitro rooting by direct transfer of the shoots from multiplication medium showed 92 % survival.