PubMed Health⌕ Search

Biomedical subjects

K M Merz

Publications and source records attributed to K M Merz.

24 records · Page 2Linked to original sources

Interaction of small peptides with lipid bilayers.

Molecular dynamics simulations of the tripeptide Ala-Phe-Ala-O-tert-butyl interacting with dimyristoylphosphatidylcholine lipid bilayers have been carried out. The lipid and aqueous environments of the peptide, the alkyl chain order, and the lipid and peptide dynamics have been investigated with use of density profiles, radial distribution functions, alkyl chain order parameter profiles, and time correlation functions. It appears that the alkyl chain region accommodates the peptides in the bilayer with minimal perturbation to this region. The peptide dynamics in the bilayer bound form has been compared with that of the free peptide in water. The peptide structure does not vary on the simulation time scale (of the order of hundreds of picoseconds) compared with the solution structure in which a random structure is observed.

Amino Acid Sequence↗

A comparison of DMPC- and DLPE-based lipid bilayers.

A 250 ps molecular dynamics simulation of the dimyristoylphosphatidylcholine (DMPC)-based lipid bilayer, including explicit water molecules, is reported. The solvent environment of the head groups and other structural properties of the bilayer have been analyzed and compared with experimental results as well as our previous simulation of the dilauroylphosphatidylethanolamine (DLPE)-based bilayer. From this comparison we find that the solvent structure around the DMPC head group (clathrate shell) is significantly different than that around the DLPE head group (typical hydrogen bonding interactions). We have modeled the probable relationship between the different solvent environments around the R-N(CH3)3+ (DMPC) and R-NH3+ (DLPE) head groups and the different interlammelar distances in these systems by performing potential of mean force (PMF) simulations on two N(CH3)4+ and NH4+ ions in water. From the PMF simulations it appears that the differences in the hydration of the DMPC and DLPE head groups is not responsible for the differences in the hydration force observed for these systems. We also find that the orientational polarization of DLPE and DMPC is similar, which suggests that solvent polarization is not responsible for the differences in the hydration repulsion behavior observed in these systems. We also examined the order parameters for DMPC and found them to be in reasonable agreement with experiment. Given the different characteristics of the DLPE and DMPC head groups, we suggest an explanation of the differences in the interlammellar spacings of bilayers composed of these like-charged lipids. From our DLPE simulations we find that the R-NH3+ head groups can interact with the nonesterified oxygens of the phosphate group in an intraleaflet or an interleaflet manner. For the latter a "cross link" between two leaflets can be formed, which causes a stabilization of the interlamellar spacings at fairly short distances. Moreover, due to the strong intraleaflet interaction we find that the DLPE interface is relatively "flat" (as opposed to DMPC-based bilayers), which results in a surface that has regions of positive and negative charge that reside in the same plane along the bilayer normal. Based on this we propose that the DLPE bilayer interface can correlate itself with another DLPE interface by alignment of the regions of positive (or negative) charge on one leaflet with the opposite charges on the opposing leaflet.

Biophysical Phenomena↗

Binding of cyanide, cyanate, and thiocyanate to human carbonic anhydrase II.

Computer simulation techniques are used to address the question of how cyanide and related ions interact with human carbonic anhydrase II (HCAII). Spectroscopic results have suggested that cyanide is coordinated with the zinc ion, while recent X-ray results suggest that the cyanide ion is noncovalently associated with the zinc-water or zinc-hydroxide form of the enzyme. We have carried out simulations on three models in an attempt to shed light on why the spectroscopic and X-ray results differ. The first model we studied (Model I) has cyanide directly coordinated to the zinc ion, the second has it noncovalently interacting with the zinc-hydroxide (high pH) form of the enzyme (Model II), and the third has cyanide noncovalently interacting with the zinc-water (low pH) form of the enzyme (Model III). None of these models is satisfactory in explaining the available structural data obtained from X-ray crystallography. This leads us to propose an alternative model, in which HCAII hydrates HCN to form an OH-/HCN complex coordinated to the Zn ion. Ab initio calculations are consistent with this model. Based on these results we are able to explain the observed crystallographic behavior of cyanate and, by inference, thiocyanate.

Bicarbonates↗

Theoretical examination of the mechanism of aldose-ketose isomerization.

Two mechanisms for an aldose-ketose isomerization have been examined using high level ab initio and semiempirical molecular orbital methods. The proton transfer pathway via an enediol intermediate is shown to be favored in the absence of a metal ion, while the hydride transfer pathway becomes favored in the presence of a metal ion. Our calculations explain why the proton transfer pathway is operative in most aldose-ketose isomerization reactions. These calculations also provide further support for the previously proposed metal ion-mediated hydride transfer mechanism for xylose isomerase.

Aldehydes↗

Structure and dynamics of the dilauroylphosphatidylethanolamine lipid bilayer.

A 200-ps molecular dynamics (MD) simulation trajectory of a model dilauroylphosphatidylethanolamine (DLPE) bilayer in water at 315 K has been generated. Segmental order parameters, electron density profiles, and water pair distribution functions have been calculated. Comparison to experiment is made where possible. The dynamics of the system has been studied by analyzing the velocity autocorrelation functions (VAF) of both water and lipid atoms. Furthermore, the diffusive properties of water have been analyzed by computing the mean square displacement (MSD) and orientational correlation function (OCF) of water in two regions around the bilayer. The calculated order parameters show a behavior similar to the liquid crystalline phase of other bilayers, but the region around C1-C3 does not show the expected behavior. The electron density profile shows features that are characteristic of the liquid crystalline phase. The radial distribution functions suggest ordering of water near the charged head groups, which results in about 15 water molecules solvating each lipid molecule. We find from the VAF, MSD, and OCF calculation that the water molecules near the head groups of the lipid bilayer move more slowly than those further away. The VAF of the hydrocarbon chains have features of low-frequency motions that are probably cooperative nature in addition to the high-frequency motions associated with bond angle and torsional motions.

Lipid Bilayers↗

Insights into the function of the zinc hydroxide-Thr199-Glu106 hydrogen bonding network in carbonic anhydrases.

The exact functional role of the zinc hydroxide (water)-Thr199-Glu106 hydrogen bond network in the carbonic anhydrases is unknown. However, from the results of molecular dynamics simulations (MD) we are able to better define its function. From computer graphics analysis and MD simulations on the zinc hydroxide form of human carbonic anhydrase II we find that this interaction forces the hydroxide hydrogen atom to be in a "down" position relative to the deep water-binding pocket. From previous work we have found that this pocket is a high-affinity binding site for CO2. We also note that during the timescale of our simulation (126 ps) the hydrogen bonds between the hydroxide hydrogen atom and Thr199 and the one between Thr199 and Glu106 are not fluxional. We propose that the role of the zinc hydroxide (water)-Thr199-Glu106 hydrogen bond network is to lock the hydrogen atom in the down position in order to expose the CO2 molecule bound in the deep water pocket to a lone pair of the hydroxide oxygen atom. This would allow for the rapid reaction of the CO2 molecule around the zinc ion. Furthermore, if the hydroxide hydrogen atom were not locked in the down position the binding of CO2 to the deep water pocket could be interfered with by the unrestrained hydroxide hydrogen atom (e.g. the N-Zn-O-H torsion could undergo rotational transitions that would partially block the deep water pocket). In summary, the roles we ascribe to this hydrogen bonding network are (1) to allow for facile access of CO2 to the deep water pocket and (2) to allow for maximal exposure of a hydroxide oxygen lone pair to the CO2 carbon atom.

Binding Sites↗