We have enhanced our parallel molecular dynamics (MD) simulation software LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator, lammps.sandia.gov) to include many new features for accelerated simulation including articulated rigid body dynamics via coupling to the Rensselaer Polytechnic Institute code POEMS (Parallelizable Open-source Efficient Multibody Software). We use new features of the LAMMPS software package to investigate rhodopsin photoisomerization, and water model surface tension and capillary waves at the vapor-liquid interface. Finally, we motivate the recipes of MD for practitioners and researchers in numerical analysis and computational mechanics.
There are many important biological processes involving lipid bilayers on times scales beyond that accessible by atomistic simulations. We have developed coarse-grained, bead-spring models of lipid molecules to treat membrane fusion, domain formation and the general physical characteristics of lipid bilayers. A key aspect of these coarse-grained models is that the liquid nature of a bilayer is explicitly present in the simulations; the lipids diffuse far beyond their neighbors in contrast to atomistic simulations. With these models self-assembly into a bilayer starting from a random configuration of lipids and solvent is readily simulated. We have performed extensive simulations to characterize these lipid models in single component lipid bilayers. For a variety of tail lengths, the area per lipid as a function of temperature has been calculated; the liquid-gel transition has been characterized. Models have been developed for a variety of lipids including double bonds in the lipid tails. Simulation results will be presented for fusion and domain formation.
Understanding the dynamics of the membrane protein rhodopsin will have broad implications for other membrane proteins and cellular signaling processes. Rhodopsin (Rho) is a light activated G-protein coupled receptor (GPCR). When activated by ligands, GPCRs bind and activate G-proteins residing within the cell and begin a signaling cascade that results in the cell's response to external stimuli. More than 50% of all current drugs are targeted toward G-proteins. Rho is the prototypical member of the class A GPCR superfamily. Understanding the activation of Rho and its interaction with its Gprotein can therefore lead to a wider understanding of the mechanisms of GPCR activation and G-protein activation. Understanding the dark to light transition of Rho is fully analogous to the general ligand binding and activation problem for GPCRs. This transition is dependent on the lipid environment. The effect of lipids on membrane protein activity in general has had little attention, but evidence is beginning to show a significant role for lipids in membrane protein activity. Using the LAMMPS program and simulation methods benchmarked under the IBIG program, we perform a variety of allatom molecular dynamics simulations of membrane proteins.
Using molecular dynamics simulations we examine the effective interactions between two like-charged rods as a function of angle and separation. In particular, we determine how the competing electrostatic repulsions and multivalent-ion-induced attractions depend upon concentrations of simple and multivalent salts. We find that with increasing multivalent salt, the stable configuration of two rods evolves from isolated rods to aggregated perpendicular rods to aggregated parallel rods; at sufficiently high concentration, additional multivalent salt reduces the attraction. Monovalent salt enhances the attraction near the onset of aggregation and reduces it at a higher concentration of multivalent salt.
The effect of cross-linker functionality and interfacial bond density on the fracture behavior of highly cross-linked polymer networks bonded to a solid surface is studied using large-scale molecular dynamics simulations. Three different cross-linker functionalities (f = 3, 4, and 6) are considered. The polymer networks are created between two solid surfaces with the number of bonds to the surfaces varying from zero to full bonding to the network. Stress?strain curves are determined for each system from tensile pull and shear deformations. At full interfacial bond density the failure mode is cohesive. The cohesive failure stress is almost identical for shear and tensile modes. The simulations directly show that cohesive failure occurs when the number of interfacial bonds is greater than in the bulk. Decreasing the number of interfacial bonds results in cohesive to adhesive transition consistent with recent experimental results. The correspondence between the stress?strain curves at different f and the sequence of molecular deformations is obtained. The failure stress decreases with smaller f while failure strain increases with smaller f.
This LDRD project has involved the development and application of Sandia's massively parallel materials modeling software to several significant biophysical systems. They have been successful in applying the molecular dynamics code LAMMPS to modeling DNA, unstructured proteins, and lipid membranes. They have developed and applied a coupled transport-molecular theory code (Tramonto) to study ion channel proteins with gramicidin A as a prototype. they have used the Towhee configurational bias Monte-Carlo code to perform rigorous tests of biological force fields. they have also applied the MP-Sala reacting-diffusion code to model cellular systems. Electroporation of cell membranes has also been studied, and detailed quantum mechanical studies of ion solvation have been performed. In addition, new molecular theory algorithms have been developed (in FasTram) that may ultimately make protein solvation calculations feasible on workstations. Finally, they have begun implementation of a combined molecular theory and configurational bias Monte-Carlo code. They note that this LDRD has provided a basis for several new internal (e.g. several new LDRD) and external (e.g. 4 NIH proposals and a DOE/Genomes to Life) proposals.
For highly cross-linked polymer networks bonded to a solid surface, the effect of interfacial bond density and system size on interfacial fracture is studied using molecular dynamics simulations. Results for tensile and shear mode simulations are given. The correspondence between the stress-strain curve and the sequence of molecular deformations is obtained. The failure strain for a fully bonded surface is equal to the strain necessary to make taut the average of the minimal paths through the network from a bonded site on the bottom solid surface to a bonded site on the top surface. At fractional interfacial bond densities, cavities form above the nonbonded surface, yielding an inhomogeneous strain profile and a smaller failure strain. The failure strain and stress are linearly proportional to the number of bonds at the interface except in the tensile mode when number of bonds is so few that van der Waals interactions dominate. The failure mode is successfully constructed to be interracial by limiting the interfacial bond density to be less than the bulk bond density.
Molecular dynamics simulations of a simple, bead-spring model of semiflexible polyelectrolytes such as DNA are performed. All charges are explicitly treated. Starting from extended, noncondensed conformations, condensed structures form in the simulations with tetravalent or trivalent counterions. No condensates form or are stable for divalent counterions. The mechanism by which condensates form is described. Briefly, condensation occurs because electrostatic interactions dominate entropy, and the favored Coulombic structure is a charge ordered state. Condensation is a generic phenomena and occurs for a variety of polyelectrolyte parameters. Toroids and rods are the condensate structures. Toroids form preferentially when the molecular stiffness is sufficiently strong.
The fracture of highly-crosslinked networks is investigated by molecular dynamics simulations. The network is modeled as a bead-spring polymer network between two solid surfaces. The network is dynamically formed by crosslinking an equilibrated liquid mixture. Tensile pull fracture is simulated as a function of the number of interracial bonds. The sequence of molecular structural deformations that lead to failure are determined, and the connectivity is found to strongly control the stress-strain response and failure modes. The failure strain is related to the minimal paths in the network that connect the two solid surfaces. The failure stress is a fraction of the ideal stress required to fracture all the interracial bonds, and is linearly proportional to the number of interracial bonds. By allowing only a single bond between a crosslinker and the surface, interracial failure always occurs. Allowing up to half of the crosslinker's bonds to occur with the surface, cohesive failure can occur.
