Molecular basis for microbial adhesion to geochemical surfaces: computer simulation of Pseudomonas aeruginosa adhesion to goethite

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Multiscale Model Describing Bacterial Adhesion and Detachment

Bacterial surfaces are complex structures with nontrivial adhesive properties. The physics of bacterial adhesion deviates from that of ideal colloids as a result of cell-surface roughness and because of the mechanical properties of the polymers covering the cell surface. In the present study, we develop a simple multiscale model for the interplay between the potential energy functions that characterize the cell surface biopolymers and their interaction with the extracellular environment. We then use the model to study a discrete network of bonds in the presence of significant length heterogeneities in cell-surface polymers. The model we present is able to generate force curves (both approach and retraction) that closely resemble those measured experimentally. Our results show that even small-length-scale heterogeneities can lead to macroscopically nonlinear behavior that is qualitatively and quantitatively different from the homogeneous case. We also report on the energetic consequences of such structural heterogeneity.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Amphiphile disruption of pathogen attachment at the hematite (α-Fe2O3)-water interface

Prior studies have indicated that the subsurface transport of Cryptosporidium parvum oocysts is diminished in sediments containing iron oxides and that inner-sphere complexation of oocyst surficial carboxylate plays a role in the retardation. However, the impacts of natural organic matter (NOM) remain poorly understood. In this study, we used a model anionic surfactant, sodium dodecyl sulfate (SDS), as a surrogate for amphiphilic NOM components to examine the impacts of amphiphilic components on oocyst adhesion mechanisms. We employed in situ attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectroscopy to determine the effects of SDS on the molecular bonds that mediate interactions between oocyst surficial biomolecules and hematite (α-Fe(2)O(3)) surface functional groups over a wide range of solution pH. The results show that the presence of SDS significantly diminishes Fe-carboxylate complexation, as indicated by progressive decreases in intensity of asymmetric and symmetric stretching vibrations of carboxylate [ν(as)(COO(-)) and ν(s)(COO(-))] with reaction time. In addition, one of the ν(s)(COO(-)) bands shifted from 1370 to 1418 cm(-1) upon SDS introduction, suggesting that SDS also changed the complexation mode. The data indicate that competition from the sulfonate groups (OSO(3)(-)) of SDS at α-Fe(2)O(3) surface sites is a primary mechanism resulting in decreased Fe-carboxylate complexation. Sorptive competition from amphiphilic NOM components may therefore increase the mobility of C. parvum oocysts in the environment through disruption of interfacial pathogen-mineral surface bonds.

Molecular dynamics simulations of a bacterial autotransporter: NalP fromNeisseria meningitidis

NalP is an autotransporter secretory protein found in the outer membrane of Neisseria meningitidis. The crystal structure of the NalP translocator domain revealed a transmembrane beta-barrel containing a central alpha-helix. The role of this alpha-helix, and of the conformational dynamics of the beta-barrel pore have been studied via atomistic molecular dynamics simulations. Three simulations, each of 10 ns duration, of NalP embedded within a solvated DMPC bilayer were performed. The helix was removed from the barrel interior in one simulation. The conformational stability of the protein is similar to that of other outer membrane proteins, e.g., OmpA, in comparable simulations. The transmembrane beta-barrel is stable even in the absence of the alpha-helix. Removal of the helix results in an influx of water into the pore region, suggesting the helix acts as a 'plug'. Water molecules entering the resultant pore form hydrogen bonds with the barrel lining that compensate for the loss of helix-barrel hydrogen bonds. The dimensions of the pore fluctuate over the course of the simulation revealing it to be flexible, but only wide enough to allow transport of the passenger domain in an unfolded or extended conformation. The simulations help us to understand the role of the central helix in plugging the pore and in maintaining the width of the barrel, and show that the NalP monomer is sufficient for the transport of the passenger domain in an unfolded or extended conformation.

Characterization of the outer membrane protein OprF of Pseudomonas aeruginosa in a lipopolysaccharide membrane by computer simulation

The N-terminal domain of outer membrane protein OprF of Pseudomonas aeruginosa forms a membrane spanning eight-stranded antiparallel beta-barrel domain that folds into a membrane channel with low conductance. The structure of this protein has been modeled after the crystal structure of the homologous protein OmpA of Escherichia coli. A number of molecular dynamics simulations have been carried out for the homology modeled structure of OprF in an explicit molecular model for the rough lipopolysaccharide (LPS) outer membrane of P. aeruginosa. The structural stability of the outer membrane model as a result of the strong electrostatic interactions compared with simple lipid bilayers is restricting both the conformational flexibility and the lateral diffusion of the porin in the membrane. Constricting side-chain interactions within the pore are similar to those found in reported simulations of the protein in a solvated lipid bilayer membrane. Because of the strong interactions between the loop regions of OprF and functional groups in the saccharide core of the LPS, the entrance to the channel from the extracellular space is widened compared with the lipid bilayer simulations in which the loops are extruding in the solvent. The specific electrostatic signature of the LPS membrane, which results in a net intrinsic dipole across the membrane, is found to be altered by the presence of OprF, resulting in a small electrically positive patch at the position of the channel.

GLYCAM06: a generalizable biomolecular force field. Carbohydrates

A new derivation of the GLYCAM06 force field, which removes its previous specificity for carbohydrates, and its dependency on the AMBER force field and parameters, is presented. All pertinent force field terms have been explicitly specified and so no default or generic parameters are employed. The new GLYCAM is no longer limited to any particular class of biomolecules, but is extendible to all molecular classes in the spirit of a small-molecule force field. The torsion terms in the present work were all derived from quantum mechanical data from a collection of minimal molecular fragments and related small molecules. For carbohydrates, there is now a single parameter set applicable to both alpha- and beta-anomers and to all monosaccharide ring sizes and conformations. We demonstrate that deriving dihedral parameters by fitting to QM data for internal rotational energy curves for representative small molecules generally leads to correct rotamer populations in molecular dynamics simulations, and that this approach removes the need for phase corrections in the dihedral terms. However, we note that there are cases where this approach is inadequate. Reported here are the basic components of the new force field as well as an illustration of its extension to carbohydrates. In addition to reproducing the gas-phase properties of an array of small test molecules, condensed-phase simulations employing GLYCAM06 are shown to reproduce rotamer populations for key small molecules and representative biopolymer building blocks in explicit water, as well as crystalline lattice properties, such as unit cell dimensions, and vibrational frequencies.

Glycolipid Membranes through Atomistic Simulations: Effect of Glucose and Galactose Head Groups on Lipid Bilayer Properties

Though glycolipids are involved in a multitude of cellular functions, the understanding of their atom-scale properties in lipid membranes has remained very limited due to the lack of atomistic simulations. In this work, we employ extensive simulations to characterize one-component membranes comprised of glycoglycerolipids, focusing on two common glyco head groups, namely glucose and galactose. The properties of these two glycoglycerolipid bilayers are compared in a systematic manner with membranes consisting of phosphatidylcholine (PC) or phosphatidylethanolamine (PE) lipids, whose structures aside from the head group are identical with those of the two glycolipids. We find that the glycolipid systems are characterized by a substantial number of hydrogen bonds in the head group region, leading to membrane packing that is stronger than in a PC but less significant than that in a PE bilayer. The role played by the glyco head group is especially evident in the electrostatic membrane potential, which is particularly large in the glycolipid membranes. For the same reason, the interfacial forces near glycolipid bilayers are significantly different from those found in PC and PE bilayers, affecting, e.g., the ordering of water close to the membrane. These effects are particularly important for the case of galactose, an important component in thylacoids.

Outer membrane proteins: comparing X-ray and NMR structures by MD simulations in lipid bilayers

The structures of three bacterial outer membrane proteins (OmpA, OmpX and PagP) have been determined by both X-ray diffraction and NMR. We have used multiple (7 x 15 ns) MD simulations to compare the conformational dynamics resulting from the X-ray versus the NMR structures, each protein being simulated in a lipid (DMPC) bilayer. Conformational drift was assessed via calculation of the root mean square deviation as a function of time. On this basis the 'quality' of the starting structure seems mainly to influence the simulation stability of the transmembrane beta-barrel domain. Root mean square fluctuations were used to compare simulation mobility as a function of residue number. The resultant residue mobility profiles were qualitatively similar for the corresponding X-ray and NMR structure-based simulations. However, all three proteins were generally more mobile in the NMR-based than in the X-ray simulations. Principal components analysis was used to identify the dominant motions within each simulation. The first two eigenvectors (which account for >50% of the protein motion) reveal that such motions are concentrated in the extracellular loops and, in the case of PagP, in the N-terminal alpha-helix. Residue profiles of the magnitude of motions corresponding to the first two eigenvectors are similar for the corresponding X-ray and NMR simulations, but the directions of these motions correlate poorly reflecting incomplete sampling on a approximately 10 ns timescale.

Modeling and simulations of a bacterial outer membrane protein: OprF from Pseudomonas aeruginosa

OprF is a major outer membrane protein from Pseudomonas aeruginosa, a homolog of OmpA from Escherichia coli. The N-terminal domains of both proteins have been demonstrated to form low conductance channels in lipid bilayers. Homology models, consisting of an eight-stranded beta-barrel, of the N-terminal domain OprF have been constructed based on the crystal structure of the corresponding domain from E. coli OmpA. OprF homology models have been evaluated via a set (6 x 10 ns) of simulations of the beta-barrel embedded within a solvated dimyristoyl-phosphatidylcholine (DMPC) bilayer. The conformational stability of the models is similar to that of the crystal structure of OmpA in comparable simulations. There is a degree of water penetration into the pore-like center of the OprF barrel. The presence of an acidic/basic (E8/K121) side-chain interaction within the OprF barrel may form a "gate" able to close/open a central pore. Lipid-protein interactions within the simulations were analyzed and revealed that aromatic side-chains (Trp, Tyr) of OprF interact with lipid headgroups. Overall, the behavior of the OprF model in simulations supports the suggestion that this molecule is comparable to OmpA. The simulations help to explain the mechanism of formation of low conductance pores within the outer membrane.

Monitoring metal ion binding in single-layer Pseudomonas aeruginosa biofilms using ATR-IR spectroscopy

The binding of metal ions to Pseudomonas aeruginosa PAO1 cells attached to a ZnSe surface has been observed in this research through cation exchange experiments using ATR-IR spectroscopy. A biofilm consisting of a single layer of Pseudomonas aeruginosa PAO1 cells was formed on a ZnSe prism by flowing a bacterial suspension in a 0.03 mol L(-)(1) NaNO(3) solution at pH 5.0 across its surface. Exposure of the biofilm to chromium(III) nitrate solution resulted in increases in all band absorbances. This absorbance increase has been attributed to the binding of chromium(III) to the bacterial exopolymers associated with the prism surface. The chromium(III) binding causes the exopolymers to contract and move the bacterial cell closer to the ZnSe surface. Further study of chromium(III) ion exchange using a mutant P. aeruginosa with a truncated lipopolysaccharide (LPS) chain resulted in much smaller absorbance changes. This observation supports the view that the extension of bacterial exopolymers and hence the distance of the bacterial cell from the surface is strongly influenced by environmental factors such as the presence of metal cations. Following chromium(III) cation exchange, the bacterial band absorbances remained constant even when the bacteria were washed with a 0.03 mol L(-)(1) NaNO(3) solution, indicating that the chromium(III) was irreversibly bound. Ion exchange with nickel(II) and cobalt(II) nitrate solutions within identical biofilms showed that these cations caused relatively small increases in absorbances that were reversible, indicating that nickel(II) and cobalt(II) are less strongly bound than chromium(III) within P. aeruginosa biofilms. The absence of discernible IR spectral changes with metal binding appears to indicate a predominantly electrostatic mechanism for binding of Cr(III), Ni(II), and Co(II) ions by bacteria in the early stages of biofilm formation.