Are Current Molecular Dynamics Force Fields too Helical?

Accurate force fields are essential for the success of molecular dynamics simulations. In apparent contrast to the conformational preferences of most force fields, recent NMR experiments suggest that short polyalanine peptides in water populate the polyproline II structure almost exclusively. To investigate this apparent contradiction, with its ramifications for the assessment of molecular force fields and the structure of unfolded proteins, we performed extensive simulations of Ala₅ in water (∼5 μs total time), using twelve different force fields and three different peptide terminal groups. Using either empirical or density-functional-based Karplus relations for the J-couplings, we find that most current force fields do overpopulate the α-region, with quantitative results depending on the choice of Karplus relation and on the peptide termini. Even after reweighting to match experiment, we find that Ala₅ retains significant α- and β-populations. In fact, several force fields match the experimental data well before reweighting and have a significant helical population. We conclude that radical changes to the best current force fields are not necessary, based on the NMR data. Nevertheless, experiments on short peptides open the way toward the systematic improvement of current simulation models.     

Is it possible to produce succinic acid at a low pH?

Bio-based succinate is still a matter of special emphasis in biotechnology and adjacent research areas. The vast majority of natural and engineered producers are bacterial strains that accumulate succinate under anaerobic conditions. Recently, we succeeded in obtaining an aerobic yeast strain capable of producing succinic acid at low pH. Herein, we discuss some difficulties and advantages of microbial pathways producing "succinic acid" rather than "succinate." It was concluded that the peculiar properties of the constructed yeast strain could be clarified in view of a distorted energy balance. There is evidence that in an acidic environment, the majority of the cellular energy available as ATP will be spent for proton and anion efflux. The decreased ATP:ADP ratio could essentially reduce the growth rate or even completely inhibit growth. In the same way, the preference of this elaborated strain for certain carbon sources could be explained in terms of energy balance. Nevertheless, the opportunity to exclude alkali and mineral acid waste from microbial succinate production seems environmentally friendly and cost-effective.     

The binding of inhibitors to α-chymotrypsin at alkaline pH

1. The binding of the competitive inhibitor N-acetyl-d-tryptophan amide to alpha-chymotrypsin has now been studied at pH values up to 10.6, by the technique of equilibrium dialysis. 2. This binding depends on the ionization of a group on the free enzyme with apparent pK(a) 9.3 at 5 degrees . 3. This group is tentatively identified as that responsible for an enzyme conformation change at high pH values, on which the catalytic activity of the enzyme also depends.     

Molecular systems biology at the crossroads: to know less about more,or to know more about less?

Systems biology is a rapidly evolving discipline that endeavours to understand the detailed coordinated workings of entire organisms, with the ultimate goal to detect differences between health and disease, or to understand how cells or entire organisms react to the environment. The editorial provides a critical evaluation of what molecular systems analysis can and cannot accomplish with existing methodologies, and how systems biology needs to merge with reductionism to yield a more comprehensive and mechanistically insightful model of a cell or organism.     

A Novel Implicit Solvent Model for Simulating the Molecular Dynamics of?RNA

Although molecular dynamics simulations can be accelerated by more than an order of magnitude by implicitly describing the influence of the solvent with a continuum model, most currently available implicit solvent simulations cannot robustly simulate the structure and dynamics of nucleic acids. The difficulties become exacerbated especially for RNAs, suggesting the presence of serious physical flaws in the prior continuum models for the influence of the solvent and counter ions on the nucleic acids. We present a novel, to our knowledge, implicit solvent model for simulating nucleic acids by combining the Langevin–Debye model and the Poisson–Boltzmann equation to provide a better estimate of the electrostatic screening of both the water and counter ions. Tests of the model involve comparisons of implicit and explicit solvent simulations for three RNA targets with 20, 29, and 75 nucleotides. The model provides reasonable agreement with explicit solvent simulations, and directions for future improvement are noted.     

Do Insect Populations Die at Constant Rates as They Become Older? Contrasting Demographic Failure Kinetics with Respect to Temperature According to the Weibull Model

Temperature implies contrasting biological causes of demographic aging in poikilotherms. In this work, we used the reliability theory to describe the consistency of mortality with age in moth populations and to show that differentiation in hazard rates is related to extrinsic environmental causes such as temperature. Moreover, experiments that manipulate extrinsic mortality were used to distinguish temperature-related death rates and the pertinence of the Weibull aging model. The Newton-Raphson optimization method was applied to calculate parameters for small samples of ages at death by estimating the maximum likelihoods surfaces using scored gradient vectors and the Hessian matrix. The study reveals for the first time that the Weibull function is able to describe contrasting biological causes of demographic aging for moth populations maintained at different temperature regimes. We demonstrate that at favourable conditions the insect death rate accelerates as age advances, in contrast to the extreme temperatures in which each individual drifts toward death in a linear fashion and has a constant chance of passing away. Moreover, slope of hazard rates shifts towards a constant initial rate which is a pattern demonstrated by systems which are not wearing out (e.g. non-aging) since the failure, or death, is a random event independent of time. This finding may appear surprising, because, traditionally, it was mostly thought as rule that in aging population force of mortality increases exponentially until all individuals have died. Moreover, in relation to other studies, we have not observed any typical decelerating aging patterns at late life (mortality leveling-off), but rather, accelerated hazard rates at optimum temperatures and a stabilized increase at the extremes.In most cases, the increase in aging-related mortality was simulated reasonably well according to the Weibull survivorship model that is applied. Moreover, semi log- probability hazard rate model illustrations and maximum likelihoods may be usefully in defining periods of mortality leveling off and provide clear evidence that environmental variability may affect parameter estimates and insect population failure rate. From a reliability theory standpoint, failure rates vary according to a linear function of age at the extremes indicating that the life system (i.e., population) is able to eliminate earlier failure and/or to keep later failure rates constant. The applied model was able to identify the major correlates of extended longevity and to suggest new ideas for using demographic concepts in both basic and applied population biology and aging.     

GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics Simulations

We introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged density maps from molecular dynamics simulations. GROmaρs efficiently computes density maps by fast multi-Gaussian spreading of atomic densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of atomic average localization during the simulations. Most importantly, it allows the comparison between computed and reference maps (e.g., experimental) through calculation of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quantitatively contrast perturbed and control simulation data sets and to examine how much biomolecular systems resemble both synthetic and experimental density maps. This was especially advantageous for multimolecule systems in which standard comparisons like RMSDs are difficult to compute. In addition, GROmaρs incorporates absolute and relative spatial free-energy estimates to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the density calculation. Furthermore, masking of regions was implemented to speed up calculations and to facilitate the comparison with experimental maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and average charge distribution in biomolecular systems. We employed all these functionalities to inspect the localization of lipid and water molecules in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, we anticipate a high applicability of GROmaρs for the analysis of molecular dynamics simulations and their comparison with experimentally determined densities.     

Molecular Dynamics Study of Binding of μ-Conotoxin GIIIA to the Voltage-Gated Sodium Channel Nav1.4

Homology models of mammalian voltage-gated sodium (Na_V) channels based on the crystal structures of the bacterial counterparts are needed to interpret the functional data on sodium channels and understand how they operate. Such models would also be invaluable in structure-based design of therapeutics for diseases involving sodium channels such as chronic pain and heart diseases. Here we construct a homology model for the pore domain of the Na_V1.4 channel and use the functional data for the binding of µ-conotoxin GIIIA to Na_V1.4 to validate the model. The initial poses for the Na_V1.4–GIIIA complex are obtained using the HADDOCK protein docking program, which are then refined in molecular dynamics simulations. The binding mode for the final complex is shown to be in broad agreement with the available mutagenesis data. The standard binding free energy, determined from the potential of mean force calculations, is also in good agreement with the experimental value. Because the pore domains of Na_V1 channels are highly homologous, the model constructed for Na_V1.4 will provide an excellent template for other Na_V1 channels.     

Molecular Dynamics Simulation of Limiting Conductance for Na2+, Cl2−, Na°, and Cl° in Supercritical Water

Abstract

We report results of molecular dynamics simulations of the limiting conductance of Na²⁺, Cl^2?, Na°, and Cl° in supercritical water using the SPC/E model for water in conjuction with our previous study (Lee et al., Chem. Phys. Lett. 293, 289 (1998)). The behavior of the limiting conductances of Na²⁺ and Cl^2? in the whole range of water density shows almost the same trend as those of Na⁺ and Cl^?, but the deviation from the assumed linear dependence of limiting conductances of Na²⁺ and Cl^2? on the water density is smaller than that of Na⁺ and Cl^?. The ratio of the limiting conductance of the divalentions to that of the corresponding monovalentions over the whole range of water density is almost constant. In the cases of Na²⁺ and Cl^2?, the dominating factor of the number of hydration water molecules around ions in the higher-density region and the dominating factor of the interaction strength between the ions and the hydration water molecules in the lower-density region are also found as was the cases for Na⁺ and Cl^?. These factors, however, are not so strong as for the corresponding monovalent ions because the change in the energetics, structure, and dynamics are very small mainly due to the strong Coulomb interaction of the divalent ions with the hydration water molecules. The diffusion coefficient of Na° and Cl° monotonically increases with decreasing water density over the whole range of water density. The increase of the diffusion coefficient with decreasing water density is attributed only to the dramatic decrease of the hydration number of water in the first solvation shell around the uncharged species. Among the two important competing factors in the limiting conductance of Na⁺ and Cl^?, the effect of the number of hydration water molecules around the uncharged species is the only existing factor over the whole range of water density since the interaction strength between the uncharged species and the hydration water molecules very small through the LJ interaction. This result has confirmed the dominating factor of the number of hydration water molecules around ions in the higher-density region in the explanation of the limiting conductance of Na⁺ and Cl^? in supercritical water at 673 K.     

Introducing weights to life cycle sustainability assessment—how do decision-makers weight sustainability dimensions?

Purpose

Decisions based on life cycle sustainability assessment (LCSA) pose a multi-criteria decision issue, as impacts on the three different sustainability dimensions have to be considered which themselves are often measured through several indicators. To support decision-making at companies, a method to interpret multi-criteria assessment and emerging trade-offs would be beneficial. This research aims at enabling decision-making within LCSA by introducing weights to the sustainability dimensions.

Methods

To derive weights, 54 decision-makers of different functions at a German automotive company were asked via limit conjoint analysis how they ranked the economic, environmental, and social performance of a vehicle component. Results were evaluated for the entire sample and by functional clusters. Additionally, sustainability respondents, i.e., respondents that dealt with sustainability in their daily business, were contrasted with non-sustainability respondents. As a last step, the impact of outliers was determined. From this analysis, practical implications for ensuring company-optimal decision-making in regard to product sustainability were derived.

Results and discussion

The results showed a large spread in weighting without clear clustering. On average, all sustainability dimensions were considered almost equally important: the economic dimension tallied at 33.5%, the environmental at 35.2%, and the social at 31.2%. Results were robust as adjusting for outliers changed weights on average by less than 10%. Results by function showed low consistency within clusters hinting that weighting was more of a personal than a functional issue. Sustainability respondents weighted the social before the environmental and economic dimension while non-sustainability respondents put the economic before the other two dimensions. Provided that the results of this research could be generalized, the retrieved weighting set was seen as a good way to introduce weights into an operationalized LCSA framework as it represented the quantification of the already existing decision process. Therefore, the acceptance of this weighting set within the respective company was expected to be increased.

Conclusions

It could be shown that conjoint analysis enabled decision-making within LCSA by introducing weights to solve a multi-criteria decision issue. Furthermore, implications for practitioners could be derived to ensure company-optimal decision-making related to product sustainability. Future research should look at expanding the sample size and geographical scope as well as investigating the weighting of indicators within sustainability dimensions and the drivers that influence personal decision-making in regard to weighting sustainability dimensions.

     

Disturb or Stabilize? A Molecular Dynamics Study of the Effects of Resorcinolic Lipids on Phospholipid Bilayers

Resorcinolic lipids, or resorcinols, are commonly found in plant membranes. They consist of a substituted benzene ring forming the hydrophilic lipid head, attached to an alkyl chain forming the hydrophobic tail. Experimental results show alternative effects of resorcinols on lipid membranes. Depending on whether they are added to lipid solutions before or after the formation of the liposomes, they either stabilize or destabilize these liposomes. Here we use atomistic molecular dynamics simulations to elucidate the molecular nature of this dual effect. Systems composed of either one of three resorcinol homologs, differing in the alkyl tail length, interacting with dimyristoylphosphatidylcholine lipid bilayers were studied. It is shown that resorcinols preincorporated into bilayers induce order within the lipid acyl chains, decrease the hydration of the lipid headgroups, and make the bilayers less permeable to water. In contrast, simulations in which the resorcinols are incorporated from the aqueous solution into a preformed phospholipid bilayer induce local disruption, leading to either transient pore formation or even complete rupture of the membrane. In line with the experimental data, our simulations thus demonstrate that resorcinols can either disturb or stabilize the membrane structure, and offer a detailed view of the underlying molecular mechanism.     

Molecular Dynamics Simulations to Gain Insights into the Stability and Morphologies of K3 Oligomers from β2-microglobulin

Abstract

β2-Microglobulin (β2-m) forms amyloid fibrils in patients undergoing long-term hemodialysis. K3 peptide, a Ser20-Lys41 fragment of β2-m, has been known to form fibrils over a wide range of pH and solvent conditions. Recent solid-state NMR has revealed that K3 oligomer adopts a parallel U-shaped β-strand-turn-β-strand motif. In order to investigate the stability and morphologies of K3 oligomers with different sizes (dimer, trimer, and tetrameri and organizations (single and double layers), several all-atom molecular dynamics simulations were conducted at 310 K and pH 2 in water and 2,2,2-trifluoroethanol (TFE). For single-layered organizations, our results show that TFE destabilizes the stacking of K3 peptides due to the fact that TFE weakens the intermolecular hydrophobic interactions of K3 oligomers. In addition, we also identified that the loop region is stabilized by the hydrophobic cluster involving resides Y7, Fll, and I16. Our results further suggest that K3 tetramer is a potential minimal nucleus seed for the formation of K3 protofibrils. For dou-ble-layered organizations in water, our data demonstrate that K3 peptides can form various stable assemblies through different interfacial arrangements, such as NN, NC, and CC, by different driving forces. We further propose that the stacking of different interfaces between two facing β-sheets of K3 peptides could be related to different fibril morphologies, which is in good agreement with the previous experimental results, showing that K3 protofibrils associated to formed mature fibrils with a wide range of diameters from 4 to 15 nm when they were transferred from 20% (v/v) TFE to aqueous solution.     

Rutland Water — from conception to operation

Rutland Water was planned in the 1960s to meet the water needs of the expanding East Midland region of England. Chosen from 64 possible sites, it was approved in 1970 and construction started in 1971. The dam was closed and the reservoir started to fill from its natural catchment in 1975; pumping started from the rivers Welland and Nene in 1976 and the reservoir was full by March 1979.Details of the stages in its development are given together with an indication of the ways in which the recreational uses of the reservoir have been integrated with its water supply function.     

Molecular Dynamics Simulation Reveals the Selective Binding of Human Leukocyte Antigen Alleles Associated with Beh?et's Disease

Behçet’s disease (BD), a multi-organ inflammatory disorder, is associated with the presence of the human leukocyte antigen (HLA) HLA-B*51 allele in many ethnic groups. The possible antigen involvement of the major histocompatibility complex class I chain related gene A transmembrane (MICA-TM) nonapeptide (AAAAAIFVI) has been reported in BD symptomatic patients. This peptide has also been detected in HLA-A*26:01 positive patients. To investigate the link of BD with these two specific HLA alleles, molecular dynamics (MD) simulations were applied on the MICA-TM nonapeptide binding to the two BD-associated HLA alleles in comparison with the two non-BD-associated HLA alleles (B*35:01 and A*11:01). The MD simulations were applied on the four HLA/MICA-TM peptide complexes in aqueous solution. As a result, stabilization for the incoming MICA-TM was found to be predominantly contributed from van der Waals interactions. The P2/P3 residue close to the N-terminal and the P9 residue at the C-terminal of the MICA-TM nonapeptide served as the anchor for the peptide accommodated at the binding groove of the BD associated HLAs. The MM/PBSA free energy calculation predicted a stronger binding of the HLA/peptide complexes for the BD-associated HLA alleles than for the non-BD-associated ones, with a ranked binding strength of B*51:01 > B*35:01 and A*26:01 > A*11:01. Thus, the HLAs associated with BD pathogenesis expose the binding efficiency with the MICA-TM nonapeptide tighter than the non-associated HLA alleles. In addition, the residues 70, 73, 99, 146, 147 and 159 of the two BD-associated HLAs provided the conserved interaction for the MICA-TM peptide binding.     

Hygiene,Sanitation, and Water: What Needs to Be Done?

In the final article in a four-part PLoS Medicine series on water and sanitation, Sandy Cairncross and colleagues outline what needs to be done to make significant progress in providing more and better hygiene, sanitation, and water for all.     

Molecular biology of cell nonpermissiveness to retroviruses. Has the time come?

The problem of host cell nonpermissiveness to retrovirus infection is characterized and illustrated on several retroviral models, including the role of viral receptors, cell fusion, and endogenous retroviral genomes as modifiers of the outcome of retroviral infection. Special attention is paid to different barriers against the infection of mammalian cells with avian leukosis/sarcoma viruses (ALV/ASV). Even when avian retroviruses become integrated in mammalian cells, several blocks at the level of provirus expression, processing of viral RNAs, and posttranslational modification prevent virus production in such virogenic cells. The significance of these blocks and new strategies making it possible to overcome some of them are discussed in relation to the development of ALV/ASV-based vectors suitable for gene therapy in mammals.     

Applying Molecular Crowding Models to Simulations of Virus Capsid Assembly In?Vitro

Virus capsid assembly has been widely studied as a biophysical system, both for its biological and medical significance and as an important model for complex self-assembly processes. No current technology can monitor assembly in detail and what information we have on assembly kinetics comes exclusively from in vitro studies. There are many differences between the intracellular environment and that of an in vitro assembly assay, however, that might be expected to alter assembly pathways. Here, we explore one specific feature characteristic of the intracellular environment and known to have large effects on macromolecular assembly processes: molecular crowding. We combine prior particle simulation methods for estimating crowding effects with coarse-grained stochastic models of capsid assembly, using the crowding models to adjust kinetics of capsid simulations to examine possible effects of crowding on assembly pathways. Simulations suggest a striking difference depending on whether or not a system uses nucleation-limited assembly, with crowding tending to promote off-pathway growth in a nonnucleation-limited model but often enhancing assembly efficiency at high crowding levels even while impeding it at lower crowding levels in a nucleation-limited model. These models may help us understand how complicated assembly systems may have evolved to function with high efficiency and fidelity in the densely crowded environment of the cell.     

Molecular Determinants of Mechanical Properties of V.?cholerae Biofilms at?the Air-Liquid Interface

Biofilm formation increases both the survival and infectivity of Vibrio cholerae, the causative agent of cholera. V. cholerae is capable of forming biofilms on solid surfaces and at the air-liquid interface, termed pellicles. Known components of the extracellular matrix include the matrix proteins Bap1, RbmA, and RbmC, an exopolysaccharide termed Vibrio polysaccharide, and DNA. In this work, we examined a rugose strain of V. cholerae and its mutants unable to produce matrix proteins by interfacial rheology to compare the evolution of pellicle elasticity in real time to understand the molecular basis of matrix protein contributions to pellicle integrity and elasticity. Together with electron micrographs, visual inspection, and contact angle measurements of the pellicles, we defined distinct contributions of the matrix proteins to pellicle morphology, microscale architecture, and mechanical properties. Furthermore, we discovered that Bap1 is uniquely required for the maintenance of the mechanical strength of the pellicle over time and contributes to the hydrophobicity of the pellicle. Thus, Bap1 presents an important matrix component to target in the prevention and dispersal of V. cholerae biofilms.