Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

Litter quality and its response to water level drawdown in boreal peatlands at plant species and community level

Changes in the structure of plant communities may have much more impact on ecosystem carbon (C) cycling than any phenotypic responses to environmental changes. We studied these impacts via the response of plant litter quality, at the level of species and community, to persistent water-level (WL) drawdown in peatlands. We studied three sites with different nutrient regimes, and water-level manipulations at two time scales. The parameters used to characterize litter quality included extractable substances, cellulose, holocellulose, composition of hemicellulose (neutral sugars, uronic acids), Klason lignin, CuO oxidation phenolic products, and concentrations of C and several nutrients. The litters formed four chemically distinct groups: non-graminoid foliar litters, graminoids, mosses and woody litters. Direct effects of WL drawdown on litter quality at the species level were overruled by indirect effects via changes in litter type composition. The pristine conditions were characterized by Sphagnum moss and graminoid litters. Short-term (years) responses of the litter inputs to WL drawdown were small. In longterm (decades), total litter inputs increased, due to increased tree litter inputs. Simultaneously, the litter type composition and its chemical quality at the community level greatly changed. The changes that we documented will strongly affect soil properties and C cycle of peatlands.Peer reviewe

Shedding light on plant litter decomposition: Advances, implications and new directions in understanding the role of photodegradation

Litter decomposition contributes to one of the largest fluxes of carbon (C) in the terrestrial biosphere and is a primary control on nutrient cycling. The inability of models using climate and litter chemistry to predict decomposition in dry environments has stimulated investigation of non-traditional drivers of decomposition, including photodegradation, the abiotic decomposition of organic matter via exposure to solar radiation. Recent work in this developing field shows that photodegradation may substantially influence terrestrial C fluxes, including abiotic production of carbon dioxide, carbon monoxide and methane, especially in arid and semi-arid regions. Research has also produced contradictory results regarding controls on photodegradation. Here we summarize the state of knowledge about the role of photodegradation in litter decomposition and C cycling and investigate drivers of photodegradation across experiments using a meta-analysis. Overall, increasing litter exposure to solar radiation increased mass loss by 23% with large variation in photodegradation rates among and within ecosystems. This variation was tied to both litter and environmental characteristics. Photodegradation increased with litter C to nitrogen (N) ratio, but not with lignin content, suggesting that we do not yet fully understand the underlying mechanisms. Photodegradation also increased with factors that increased solar radiation exposure (latitude and litter area to mass ratio) and decreased with mean annual precipitation. The impact of photodegradation on C (and potentially N) cycling fundamentally reshapes our thinking of decomposition as a solely biological process and requires that we define the mechanisms driving photodegradation before we can accurately represent photodegradation in global C and N models. © 2012 US Government

Quantitative Analysis of Lipid Droplet Fusion: Inefficient Steady State Fusion but Rapid Stimulation by Chemical Fusogens

Lipid droplets (LDs) are dynamic cytoplasmic organelles containing neutral lipids and bounded by a phospholipid monolayer. Previous studies have suggested that LDs can undergo constitutive homotypic fusion, a process linked to the inhibitory effects of fatty acids on glucose transporter trafficking. Using strict quantitative criteria for LD fusion together with refined light microscopic methods and real-time analysis, we now show that LDs in diverse cell types show low constitutive fusogenic activity under normal growth conditions. To investigate the possible modulation of LD fusion, we screened for agents that can trigger fusion. A number of pharmacological agents caused homotypic fusion of lipid droplets in a variety of cell types. This provided a novel cell system to study rapid regulated fusion between homotypic phospholipid monolayers. LD fusion involved an initial step in which the two adjacent membranes became continuous (<10 s), followed by the slower merging (100 s) of the neutral lipid cores to produce a single spherical LD. These fusion events were accompanied by changes to the LD surface organization. Measurements of LDs undergoing homotypic fusion showed that fused LDs maintained their initial volume, with a corresponding decrease in surface area suggesting rapid removal of membrane from the fused LD. This study provides estimates for the level of constitutive LD fusion in cells and questions the role of LD fusion in vivo. In addition, it highlights the extent of LD restructuring which occurs when homotypic LD fusion is triggered in a variety of cell types

Bordetella Adenylate Cyclase Toxin Mobilizes Its β2 Integrin Receptor into Lipid Rafts to Accomplish Translocation across Target Cell Membrane in Two Steps

Bordetella adenylate cyclase toxin (CyaA) binds the αMβ2 integrin (CD11b/CD18, Mac-1, or CR3) of myeloid phagocytes and delivers into their cytosol an adenylate cyclase (AC) enzyme that converts ATP into the key signaling molecule cAMP. We show that penetration of the AC domain across cell membrane proceeds in two steps. It starts by membrane insertion of a toxin ‘translocation intermediate’, which can be ‘locked’ in the membrane by the 3D1 antibody blocking AC domain translocation. Insertion of the ‘intermediate’ permeabilizes cells for influx of extracellular calcium ions and thus activates calpain-mediated cleavage of the talin tether. Recruitment of the integrin-CyaA complex into lipid rafts follows and the cholesterol-rich lipid environment promotes translocation of the AC domain across cell membrane. AC translocation into cells was inhibited upon raft disruption by cholesterol depletion, or when CyaA mobilization into rafts was blocked by inhibition of talin processing. Furthermore, CyaA mutants unable to mobilize calcium into cells failed to relocate into lipid rafts, and failed to translocate the AC domain across cell membrane, unless rescued by Ca2+ influx promoted in trans by ionomycin or another CyaA protein. Hence, by mobilizing calcium ions into phagocytes, the ‘translocation intermediate’ promotes toxin piggybacking on integrin into lipid rafts and enables AC enzyme delivery into host cytosol

Regulation of the V-ATPase along the Endocytic Pathway Occurs through Reversible Subunit Association and Membrane Localization

The lumen of endosomal organelles becomes increasingly acidic when going from the cell surface to lysosomes. Luminal pH thereby regulates important processes such as the release of internalized ligands from their receptor or the activation of lysosomal enzymes. The main player in endosomal acidification is the vacuolar ATPase (V-ATPase), a multi-subunit transmembrane complex that pumps protons from the cytoplasm to the lumen of organelles, or to the outside of the cell. The active V-ATPase is composed of two multi-subunit domains, the transmembrane V0 and the cytoplasmic V1. Here we found that the ratio of membrane associated V1/Vo varies along the endocytic pathway, the relative abundance of V1 being higher on late endosomes than on early endosomes, providing an explanation for the higher acidity of late endosomes. We also found that all membrane-bound V-ATPase subunits were associated with detergent resistant membranes (DRM) isolated from late endosomes, raising the possibility that association with lipid-raft like domains also plays a role in regulating the activity of the proton pump. In support of this, we found that treatment of cells with U18666A, a drug that leads to the accumulation of cholesterol in late endosomes, affected acidification of late endosome. Altogether our findings indicate that the activity of the vATPase in the endocytic pathway is regulated both by reversible association/dissociation and the interaction with specific lipid environments

Integrative omics identifies conserved and pathogen-specific responses of sepsis-causing bacteria

Even in the setting of optimal resuscitation in high-income countries severe sepsis and septic shock have a mortality of 20–40%, with antibiotic resistance dramatically increasing this mortality risk. To develop a reference dataset enabling the identification of common bacterial targets for therapeutic intervention, we applied a standardized genomic, transcriptomic, proteomic and metabolomic technological framework to multiple clinical isolates of four sepsis-causing pathogens: Escherichia coli, Klebsiella pneumoniae species complex, Staphylococcus aureus and Streptococcus pyogenes. Exposure to human serum generated a sepsis molecular signature containing global increases in fatty acid and lipid biosynthesis and metabolism, consistent with cell envelope remodelling and nutrient adaptation for osmoprotection. In addition, acquisition of cholesterol was identified across the bacterial species. This detailed reference dataset has been established as an open resource to support discovery and translational research

The ELBA Force Field for Coarse-Grain Modeling of Lipid Membranes

A new coarse-grain model for molecular dynamics simulation of lipid membranes is presented. Following a simple and conventional approach, lipid molecules are modeled by spherical sites, each representing a group of several atoms. In contrast to common coarse-grain methods, two original (interdependent) features are here adopted. First, the main electrostatics are modeled explicitly by charges and dipoles, which interact realistically through a relative dielectric constant of unity (). Second, water molecules are represented individually through a new parametrization of the simple Stockmayer potential for polar fluids; each water molecule is therefore described by a single spherical site embedded with a point dipole. The force field is shown to accurately reproduce the main physical properties of single-species phospholipid bilayers comprising dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylethanolamine (DOPE) in the liquid crystal phase, as well as distearoylphosphatidylcholine (DSPC) in the liquid crystal and gel phases. Insights are presented into fundamental properties and phenomena that can be difficult or impossible to study with alternative computational or experimental methods. For example, we investigate the internal pressure distribution, dipole potential, lipid diffusion, and spontaneous self-assembly. Simulations lasting up to 1.5 microseconds were conducted for systems of different sizes (128, 512 and 1058 lipids); this also allowed us to identify size-dependent artifacts that are expected to affect membrane simulations in general. Future extensions and applications are discussed, particularly in relation to the methodology's inherent multiscale capabilities

Sensing the fuels: glucose and lipid signaling in the CNS controlling energy homeostasis

The central nervous system (CNS) is capable of gathering information on the body’s nutritional state and it implements appropriate behavioral and metabolic responses to changes in fuel availability. This feedback signaling of peripheral tissues ensures the maintenance of energy homeostasis. The hypothalamus is a primary site of convergence and integration for these nutrient-related feedback signals, which include central and peripheral neuronal inputs as well as hormonal signals. Increasing evidence indicates that glucose and lipids are detected by specialized fuel-sensing neurons that are integrated in these hypothalamic neuronal circuits. The purpose of this review is to outline the current understanding of fuel-sensing mechanisms in the hypothalamus, to integrate the recent findings in this field, and to address the potential role of dysregulation in these pathways in the development of obesity and type 2 diabetes mellitus