Biomedical subjects

Mark J Stevens

Publications and source records attributed to Mark J Stevens.

12 recordsLinked to original sources

PubMed2007

How a small change in retinal leads to G-protein activation: initial events suggested by molecular dynamics calculations.

Rhodopsin is the prototypical G-protein coupled receptor, coupling light activation with high efficiency to signaling molecules. The dark-state X-ray structures of the protein provide a starting point for consideration of the relaxation from initial light activation to conformational changes that may lead to signaling. In this study we create an energetically unstable retinal in the light activated state and then use molecular dynamics simulations to examine the types of compensation, relaxation, and conformational changes that occur following the cis-trans light activation. The results suggest that changes occur throughout the protein, with changes in the orientation of Helices 5 and 6, a closer interaction between Ala 169 on Helix 4 and retinal, and a shift in the Schiff base counterion that also reflects changes in sidechain interactions with the retinal. Taken together, the simulation is suggestive of the types of changes that lead from local conformational change to light-activated signaling in this prototypical system.

Alanine↗

PubMed2006

Capillary waves at the liquid-vapor interface and the surface tension of water.

Capillary waves occurring at the liquid-vapor interface of water are studied using molecular dynamics simulations. In addition, the surface tension, determined thermodynamically from the difference in the normal and tangential pressure at the liquid-vapor interface, is compared for a number of standard three- and four-point water models. We study four three-point models (SPC/E, TIP3P, TIP3P-CHARMM, and TIP3P-Ew) and two four-point models (TIP4P and TIP4P-Ew). All of the models examined underestimate the surface tension; the TIP4P-Ew model comes closest to reproducing the experimental data. The surface tension can also be determined from the amplitude of capillary waves at the liquid-vapor interface by varying the surface area of the interface. The surface tensions determined from the amplitude of the logarithmic divergence of the capillary interfacial width and from the traditional thermodynamic method agree only if the density profile is fitted to an error function instead of a hyperbolic tangent function.

Journal Article↗

PubMed2005

Tribological properties of alkylsilane self-assembled monolayers.

In this study, we perform molecular dynamics simulations of adhesive contact and friction between alkylsilane Si(OH)(3)(CX(2))(10)CX(3) and alkoxylsilane Si(OH)(2)(CX(2))(10)CX(3) (where X = H or F) self-assembled monolayers (SAMs) on an amorphous silica substrate. The alkylsilane SAMs are primarily hydrogen-bonded or physisorbed to the surface. The alkoxylsilane SAMs are covalently bonded or chemisorbed to the surface. Previously, we studied the chemisorbed systems. In this work, we study the physisorbed systems and compare the tribological properties with the chemisorbed systems. Furthermore, we examine how water at the interface of the SAMs and substrate affects the tribological properties of the physisorbed systems. When less than a third of a monolayer is present, very little difference in the microscopic friction coefficient mu or shear stresses is observed. For increasing amounts of water, the values of mu and the shear stresses decrease; this effect is somewhat more pronounced for fluorocarbon alkylsilane SAMs than for the hydrocarbon SAMs. The observed decrease in friction is a consequence of a slip plane that occurs in the water as the amount of water is increased. We studied the frictional behavior using relative shear velocities ranging from v = 2 cm/s to 2 m/s. Similar to previously reported results for alkoxylsilane SAMs, the values of the measured stress and mu for the alkylsilane SAM systems decrease monotonically with v.

Journal Article↗

PubMed2005

Complementary matching in domain formation within lipid bilayers.

Domain structure and formation in lipid bilayers are investigated by molecular dynamics simulations using a coarse-grained lipid model. The lipid bilayers consist of two lipid types that are identical except for tail length. At a temperature intermediate to the two melting temperatures of the constituent lipid types, gel domains spontaneously form from an initial random structure. The simulations reveal that the gel domains consist of both lipid types in a complementary match. If a long lipid is in the top monolayer, then a short lipid is underneath and vice versa. The gel domains have a larger thickness than the surrounding liquid phase. The thickness of the gel domains is close to that of the pure long lipid gel phase bilayers. However, since in the mixed gel domains the lipids are not tilted and in the pure gel phase the lipids are tilted, the two thicknesses are similar, and the underlying structure is therefore not distinguishable solely by thickness measurements.

Gels↗

PubMed2005

Molecular dynamics simulations of grafted polyelectrolytes on two apposing walls.

Molecular dynamics simulations of polyelectrolytes grafted to two apposing surfaces were performed. Bead-spring polymer models are used to treat flexible chains [e.g., sodium poly(styrene sulfonate)] and stiff chains (double-stranded DNA). The counterions are explicitly treated. The effect of the surface density of the grafted polymer, the chain length, and the gap width on the structure and the pressure were studied. Results are compared to experimental measurements and to simulations of polyelectrolyte brushes on a single surface. The density profiles exhibit a maximum not found in single surface data. The maximum is due to the brushes shrinking to avoid interpenetration.

Computer Simulation↗

PubMed2004

Coarse-grained simulations of lipid bilayers.

A minimal model of lipid molecules consisting of bead-spring representation is developed. The basic interactions are hydrophobic and polar interactions. Essential physical features of lipid bilayers are maintained using this model, and relatively long times can be simulated in comparison to atomistic models. Self-assembly from a random starting configuration to a bilayer can readily be followed using molecular dynamics simulations. The diffusion of lipid molecules well beyond their nearest neighbors is attained. As a basis for description of the model, the area per lipid, the bending modulus, and the area compressibility as a function of temperature and tail length are calculated. A liquid to gel transition is observed and quantitatively characterized. Both saturated and unsaturated lipids are treated.

Computer Simulation↗

PubMed2004

Frictional dynamics of fluorine-terminated alkanethiol self-assembled monolayers.

The frictional dynamics of fluorine-terminated alkanethiol (S(CH2)8CF3) self-assembled monolayers (SAMs) on gold are studied using molecular dynamics simulations. The simulations treat the interactions between two SAMs on flat surfaces. The structure and frictional behavior are investigated as a function of applied pressure (200 MPa to 1 GPa) for a shear velocity of 2 m/s and compared to methyl-terminated alkanethiol SAMs. The maximum adhesive pressure between the SAMs is 220 MPa for both end groups. In agreement with experiments on the molecular scale, the shear stress and the coefficient of friction for CF3-terminated alkanethiols are larger than for CH3-terminated alkanethiols. The main source for the difference is primarily the tighter packing of the fluorinated terminal group resulting in a higher degree of order. The molecular scale coefficient of friction is correlated with the degree of order among all the systems.

Journal Article↗

PubMed2004

Effect of mono- and multivalent salts on angle-dependent attractions between charged rods.

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.

Journal Article↗

PubMed2004

How environment supports a state: molecular dynamics simulations of two states in bacteriorhodopsin suggest lipid and water compensation.

The light-driven proton pump bacteriorhodopsin (bR) is a transmembrane protein that uses large conformational changes for proton transfer from the cytoplasmic to the extracellular regions. Crystal structures, due to their solvent conditions, do not resolve the effect of lipid molecules on these protein conformational changes. To begin to understand the molecular details behind such large conformational changes, we simulated two conformations of wild-type bacteriorhodopsin, one of the dark-adapted state and the second of an intermediate (M(O)) state, each within an explicit dimyristoyl-phosphatidylcholine (DMPC) lipid bilayer. The simulations included all-hydrogen and all-atom representations of protein, lipid, and water and were performed for 20 ns. We investigate the equilibrium properties and the dynamic motions of the two conformations in the lipid setting. We note that the conformational state of the M(O) intermediate bR remains markedly different from the dark-adapted bR state in that the M(O) intermediate shows rearrangement of the cytoplasmic portions of helices C, F, and G, and nearby loops. This difference in the states remained throughout the simulations, and the results are stable on the molecular dynamics timescale and provide an illustration of the changes in both lipid and water that help to stabilize a particular state. Our analysis focuses on how the environment adjusts to these two states and on how the dynamics of the helices, loops, and water molecules can be related to the pump mechanism of bacteriorhodopsin. For example, water generally behaves in the same manner on the extracellular sides of both simulations but is decreased in the cytoplasmic region of the M(O) intermediate. We suspect that the different water behavior is closely related to the fluctuations of microcavities volume in the protein interior, which is strongly coupled to the collective motion of the protein. Our simulation result suggests that experimental observation can be useful to verify a decreased number of waters in the cytoplasmic regions of the late-intermediate stages by measuring the rate of water exchange with the interior of the protein.

Bacteriorhodopsins↗

PubMed2003

Insights into the molecular mechanism of membrane fusion from simulation: evidence for the association of splayed tails.

We present coarse-grained simulations of fusion between two liposomes from which a detailed picture of lipid movements emerges. In these simulations the bilayers dilate at the contact edge, and the resulting increase in the area per lipid produces a tilting of the individual molecules as predicted. Fusion is initiated when some of these tilted lipids splay their aliphatic tails, such that the molecules are shared between the opposing leaflets. Multiple splayed lipids subsequently associate with their aliphatic tails in contact, which produces a new hydrophobic core. As the tails extend into a more parallel conformation the two outer leaflets become contiguous to produce a hemifused structure. The results have interesting implications for biological membrane fusion and suggest new possibilities for designing molecules that control fusion.

Computer Simulation↗

PubMed2003

Molecular dynamics simulation of dark-adapted rhodopsin in an explicit membrane bilayer: coupling between local retinal and larger scale conformational change.

The light-driven photocycle of rhodopsin begins the photoreceptor cascade that underlies visual response. In a sequence of events, the retinal covalently attached to the rhodopsin protein undergoes a conformational change that communicates local changes to a global conformational change throughout the whole protein. In turn, the large-scale protein change then activates G-proteins and signal amplification throughout the cell. The nature of this change, involving a coupling between a local process and larger changes throughout the protein, may be important for many membrane proteins. In addition, functional work has shown that this coupling occurs with different efficiency in different lipid settings. To begin to understand the nature of the efficiency of this coupling in different lipid settings, we present a molecular dynamics study of rhodopsin in an explicit dioleoyl-phosphatidylcholine bilayer. Our system was simulated for 40 ns and provides insights into the very early events of the visual cascade, before the full transition and activation have occurred. In particular, we see an event near 10 ns that begins with a change in hydrogen bonding near the retinal and that leads through a series of coupled changes to a shift in helical tilt. This type of event, though rare on the molecular dynamics time-scale, could be an important clue to the types of coupling that occur between local and large-scale conformational change in many membrane proteins.

Animals↗

PubMed2003

Fracture behavior of Lennard-Jones glasses.

The fracture behavior of binary Lennard-Jones (LJ) glasses is studied by extensive molecular dynamics simulations. These LJ glasses represent a nonbond limit of polymer network glasses. We determine that the low strain behavior of the LJ and polymer glasses is similar. Two different LJ glasses are fractured under tensile strain without any preexisting crack. Void formation and resulting growth as strain increases is the mechanism through which the system fails. Void formation initiates at the yield strain of varepsilon (y) approximately 0.09, which is approximately the same strain at which the yielding behavior is first observed in cross-linked network models of polymer adhesives. The yield stress increased only by small amounts with increased strain rate and with increased system size (from N=30 000 atoms to 120 000 atoms). Within the ranges tested, the stress-strain behavior of these systems is independent of the temperature drop during quench and the initial molecular configuration.

Journal Article↗