Enhanced Lipid Diffusion and Mixing in Accelerated Molecular Dynamics

Accelerated molecular dynamics (aMD) is an enhanced sampling technique that expedites conformational space sampling by reducing the barriers separating various low-energy states of a system. Here, we present the first application of the aMD method on lipid membranes. Altogether, ∼1.5 μs simulations were performed on three systems: a pure POPC bilayer, a pure DMPC bilayer, and a mixed POPC:DMPC bilayer. Overall, the aMD simulations are found to produce significant speedup in trans–gauche isomerization and lipid lateral diffusion versus those in conventional MD (cMD) simulations. Further comparison of a 70-ns aMD run and a 300-ns cMD run of the mixed POPC:DMPC bilayer shows that the two simulations yield similar lipid mixing behaviors, with aMD generating a 2–3-fold speedup compared to cMD. Our results demonstrate that the aMD method is an efficient approach for the study of bilayer structural and dynamic properties. On the basis of simulations of the three bilayer systems, we also discuss the impact of aMD parameters on various lipid properties, which can be used as a guideline for future aMD simulations of membrane systems

Performance of ab initio and density functional methods for conformational equilibria of CnH2n+2 alkane isomers (n=2-8)

Conformational energies of n-butane, n-pentane, and n-hexane have been calculated at the CCSD(T) level and at or near the basis set limit. Post-CCSD(T) contribution were considered and found to be unimportant. The data thus obtained were used to assess the performance of a variety of density functional methods. Double-hybrid functionals like B2GP-PLYP and B2K-PLYP, especially with a small Grimme-type empirical dispersion correction, are capable of rendering conformational energies of CCSD(T) quality. These were then used as a `secondary standard' for a larger sample of alkanes, including isopentane and the branched hexanes as well as key isomers of heptane and octane. Popular DFT functionals like B3LYP, B3PW91, BLYP, PBE, and PBE0 tend to overestimate conformer energies without dispersion correction, while the M06 family severely underestimates GG interaction energies. Grimme-type dispersion corrections for these overcorrect and lead to qualitatively wrong conformer orderings. All of these functionals also exhibit deficiencies in the conformer geometries, particularly the backbone torsion angles. The PW6B95 and, to a lesser extent, BMK functionals are relatively free of these deficiencies. Performance of these methods is further investigated to derive conformer ensemble corrections to the enthalpy function, $H_{298}-H_0$, and the Gibbs energy function, ${\rm gef}(T)\equiv - [G(T)-H_0]/T$, of these alkanes. While $H_{298}-H_0$ is only moderately sensitive to the level of theory, ${\rm gef}(T)$ exhibits more pronounced sensitivity. Once again, double hybrids acquit themselves very well.Comment: J. Phys. Chem. A, revised [Walter Thiel festschrift

Promover a fluência em leitura: um estudo com alunos do 2º ano de escolaridade

Neste artigo são analisados os resultados de um Programa de Promoção da Fluência em Leitura (PPFL), implementado junto de 74 alunos do 2º ano de escolaridade de um Agrupamento Escolar do norte de Portugal. O PPFL é constituído por 22 sequências didáticas, elaboradas a partir de 22 textos (9 narrativos, 4 informativos e 9 poemas). Cada sequência foi operacionalizada em sessões de 10 a 15 minutos, durante 22 semanas, em ciclos de cinco dias (quinta-feira a quarta-feira). Em cada semana foi trabalhado apenas uma sequência didáctica (i.e. um texto). Foi utilizado um design quase experimental, com grupo experimental e grupo de controlo e com pré e pós-teste. Os sujeitos foram avaliados através de um teste de fluência de leitura – o Teste de Fluência em Leitura (TFL) –, considerando as variáveis velocidade e precisão. Os resultados evidenciam diferenças significativas a favor do grupo experimental, que superou as diferenças iniciais que se registavam no pré-teste. A análise qualitativa do impacto do PPFL aponta para um incremento na motivação para ler, no envolvimento da família e para a mudança de práticas por parte dos professores envolvidos.CIEC - Centro de Investigação em Estudos da Criança, IE, UMinho (UI 317 da FCT), PortugalFundos Nacionais através da FCT (Fundação para a Ciência e a Tecnologia) e cofinanciado pelo Fundo Europeu de Desenvolvimento Regional (FEDER) através do COMPETE 2020 – Programa Operacional Competitividade e Internacionalização (POCI) no âmbito do CIEC (Centro de Investigação em Estudos da Criança, da Universidade do Minho) com a referência POCI-01-0145-FEDER-00756

Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations

Molecular dynamics simulations of membrane proteins have provided deeper insights into their functions and interactions with surrounding environments at the atomic level. However, compared to solvation of globular proteins, building a realistic protein/membrane complex is still challenging and requires considerable experience with simulation software. Membrane Builder in the CHARMM-GUI website (http://www.charmm-gui.org) helps users to build such a complex system using a web browser with a graphical user interface. Through a generalized and automated building process including system size determination as well as generation of lipid bilayer, pore water, bulk water, and ions, a realistic membrane system with virtually any kinds and shapes of membrane proteins can be generated in 5 minutes to 2 hours depending on the system size. Default values that were elaborated and tested extensively are given in each step to provide reasonable options and starting points for both non-expert and expert users. The efficacy of Membrane Builder is illustrated by its applications to 12 transmembrane and 3 interfacial membrane proteins, whose fully equilibrated systems with three different types of lipid molecules (DMPC, DPPC, and POPC) and two types of system shapes (rectangular and hexagonal) are freely available on the CHARMM-GUI website. One of the most significant advantages of using the web environment is that, if a problem is found, users can go back and re-generate the whole system again before quitting the browser. Therefore, Membrane Builder provides the intuitive and easy way to build and simulate the biologically important membrane system

All-d-Enantiomer of β-Amyloid Peptide Forms Ion Channels in Lipid Bilayers

Alzheimer’s disease (AD) is the most common type of senile dementia in aging populations. Amyloid β (Aβ)-mediated dysregulation of ionic homeostasis is the prevailing underlying mechanism leading to synaptic degeneration and neuronal death. Aβ-dependent ionic dysregulation most likely occurs either directly via unregulated ionic transport through the membrane or indirectly via Aβ binding to cell membrane receptors and subsequent opening of existing ion channels or transporters. Receptor binding is expected to involve a high degree of stereospecificity. Here, we investigated whether an Aβ peptide enantiomer, whose entire sequence consists of d-amino acids, can form ion-conducting channels; these channels can directly mediate Aβ effects even in the absence of receptor–peptide interactions. Using complementary approaches of planar lipid bilayer (PLB) electrophysiological recordings and molecular dynamics (MD) simulations, we show that the d-Aβ isomer exhibits ion conductance behavior in the bilayer indistinguishable from that described earlier for the l-Aβ isomer. The d isomer forms channel-like pores with heterogeneous ionic conductance similar to the l-Aβ isomer channels, and the d-isomer channel conductance is blocked by Zn2+, a known blocker of l-Aβ isomer channels. MD simulations further verify formation of β-barrel-like Aβ channels with d- and l-isomers, illustrating that both d- and l-Aβ barrels can conduct cations. The calculated values of the single-channel conductance are approximately in the range of the experimental values. These findings are in agreement with amyloids forming Ca2+ leaking, unregulated channels in AD, and suggest that Aβ toxicity is mediated through a receptor-independent, nonstereoselective mechanism

Recovery of a US Endangered Fish

BACKGROUND: More fish have been afforded US Endangered Species Act protection than any other vertebrate taxonomic group, and none has been designated as recovered. Shortnose sturgeon (Acipenser brevirostrum) occupy large rivers and estuaries along the Atlantic coast of North America, and the species has been protected by the US Endangered Species Act since its enactment. METHODOLOGY/PRINCIPAL FINDINGS: Data on the shortnose sturgeon in the Hudson River (New York to Albany, NY, USA) were obtained from a 1970s population study, a population and fish distribution study we conducted in the late 1990s, and a fish monitoring program during the 1980s and 1990s. Population estimates indicate a late 1990s abundance of about 60,000 fish, dominated by adults. The Hudson River population has increased by more than 400% since the 1970s, appears healthy, and has attributes typical for a long-lived species. Our population estimates exceed the government and scientific population recovery criteria by more than 500%, we found a positive trend in population abundance, and key habitats have remained intact despite heavy human river use. CONCLUSIONS/SIGNIFICANCE: Scientists and legislators have called for changes in the US Endangered Species Act, the Act is being debated in the US Congress, and the Act has been characterized as failing to recover species. Recovery of the Hudson River population of shortnose sturgeon suggests the combination of species and habitat protection with patience can yield successful species recovery, even near one of the world's largest human population centers

Outer membrane protein size and LPS O-antigen define protective antibody targeting to the Salmonella surface

Lipopolysaccharide (LPS) O-antigen (O-Ag) is known to limit antibody binding to surface antigens, although the relationship between antibody, O-Ag and other outer-membrane antigens is poorly understood. Here we report, immunization with the trimeric porin OmpD from Salmonella Typhimurium (STmOmpD) protects against infection. Atomistic molecular dynamics simulations indicate this is because OmpD trimers generate footprints within the O-Ag layer sufficiently sized for a single IgG Fab to access. While STmOmpD differs from its orthologue in S. Enteritidis (SEn) by a single amino-acid residue, immunization with STmOmpD confers minimal protection to SEn. This is due to the OmpD-O-Ag interplay restricting IgG binding, with the pairing of OmpD with its native O-Ag being essential for optimal protection after immunization. Thus, both the chemical and physical structure of O-Ag are key for the presentation of specific epitopes within proteinaceous surface-antigens. This enhances combinatorial antigenic diversity in Gram-negative bacteria, while reducing associated fitness costs

Structure and Dynamics of a Fusion Peptide Helical Hairpin on the Membrane Surface: Comparison of Molecular Simulations and NMR

The conserved N-terminal residues of the HA2 subunit of influenza hemagglutinin (fusion peptide) are essential for membrane fusion and viral entry. Recent NMR studies showed that the 23-residue fusion peptide forms a helical hairpin that undergoes rocking motion relative to the membrane surface on a nanosecond time scale. To compare with NMR and to obtain a detailed molecular picture of the peptide–membrane interaction, we performed molecular dynamics simulations of the fusion peptide in explicit dimyristoylphosphatidylcholine and with the IMM1 implicit membrane model. To account for low and neutral pH conditions, simulations were performed with acidic groups (E11 and D19) protonated and unprotonated, respectively. The hairpin structure was stable in the simulations, with the N-terminal helix buried more deeply into the hydrophobic membrane interior than the C-terminal helix. Interactions between the tryptophans in the fusion peptide and phospholipid residues contribute to peptide orientation. Higher flexibility of the hairpin was observed in the implicit membrane simulations. Internal correlation functions of backbone N–H vectors were fit to the extended Lipari–Szabo model-free approach to obtain order parameters and correlation times. Good agreement with the NMR results was obtained for orientational fluctuations around the hairpin axis (rotation), but those around the perpendicular axis (tilting) were more limited in the simulations than inferred from the NMR experiments

Membrainy: a ‘smart’, unified membrane analysis tool

BACKGROUND: The study of biological membranes using Molecular Dynamics has become an increasingly popular means by which to investigate the interactions of proteins, peptides and potentials with lipid bilayers. These interactions often result in changes to the properties of the lipids which can modify the behaviour of the membrane. Membrainy is a unified membrane analysis tool that contains a broad spectrum of analytical techniques to enable: measurement of acyl chain order parameters; presentation of 2D surface and thickness maps; determination of lateral and axial headgroup orientations; measurement of bilayer and leaflet thickness; analysis of the annular shell surrounding membrane-embedded objects; quantification of gel percentage; time evolution of the transmembrane voltage; area per lipid calculations; and quantification of lipid mixing/demixing entropy. RESULTS: Each analytical component within Membrainy has been tested on a variety of lipid bilayer systems and was found to be either comparable to or an improvement upon existing software. For the analytical techniques that have no direct comparable software, our results were confirmed with experimental data. CONCLUSIONS: Membrainy is a user-friendly, intelligent membrane analysis tool that automatically interprets a variety of input formats and force fields, is compatible with both single and double bilayers, and capable of handling asymmetric bilayers and lipid flip-flopping. Membrainy has been designed for ease of use, requiring no installation or configuration and minimal user-input to operate