Studying the Effects of Dissolved Noble Gases and High Hydrostatic Pressure on the Spherical DOPC Bilayer Using Molecular Dynamic Simulations.

Fine-grained molecular dynamics simulations have been conducted to depict lipid objects enclosed in water and interacting with a series of noble gases dissolved in the medium. The simple point-charge (SPC) water system, featuring a boundary composed of 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) molecules, maintained stability throughout the simulation under standard conditions. This allowed for the accurate modeling of the effects of hydrostatic pressure at an ambient pressure of 25 bar. The chosen pressure references the 240 m depth of seawater: the horizon frequently used by commercial divers, who comprise the primary patient population of the neurological complication of inert gas narcosis and the consequences of high-pressure neurological syndrome. To quantify and validate the neurological effects of noble gases and discriminate them from high hydrostatic pressure, we reduced the dissolved gas molar concentration to 1.5%, three times smaller than what we previously tested for the planar bilayer (3.5%). The nucleation and growth of xenon, argon and neon nanobubbles proved consistent with the data from the planar bilayer simulations. On the other hand, hyperbaric helium induces only a residual distorting effect on the liposome, with no significant condensed gas fraction observed within the hydrophobic core. The bubbles were distributed over a large volume-both in the bulk solvent and in the lipid phase-thereby causing substantial membrane distortion. This finding serves as evidence of the validity of the multisite distortion hypothesis for the neurological effect of inert gases at high pressure.

Full-Atomistic Optimized Potentials for Liquid Simulations and Polymer Consistent Force Field Models for Biocompatible Shape-Memory Poly(ε-caprolactone).

Thermally induced shape memory poly(ε-caprolactone) (PCL)-based polymers are one of the most extensively researched families of biocompatible materials. They are degradable under physiological conditions and have high applicability in general biomedical engineering, with cross-linked PCL networks being particularly useful for tissue engineering. In this study, we used the optimized potentials for liquid simulations (OPLS) force field, which is well suited for describing intermolecular interactions in biomolecules, and the class II polymer consistent force field (PCFF) to investigate the properties of telechelic PCL with diacrylates as reactive functionalities on its end groups. PCFF has been specifically parameterized for simulating synthetic polymeric materials. We compare the findings of all-atom molecular dynamics simulations with known experimental data and theoretical assumptions to verify the applicability of both these force fields. We estimated the melt density, volume, transition temperatures, and mechanical characteristics of two-branched PCL diacrylates with a molecular weight of 2481 Da. Our findings point to the utility of the aforementioned force fields in predicting the properties of PCL-based polymers. It also opens avenues for developing PCL cross-linked polymer models and employing OPLS to investigate the interactions of synthetic polymers with biomolecules.

The effect of high pressure on the NMDA receptor: molecular dynamics simulations.

Professional divers exposed to ambient pressures above 11 bar develop the high pressure neurological syndrome (HPNS), manifesting as central nervous system (CNS) hyperexcitability, motor disturbances, sensory impairment, and cognitive deficits. The glutamate-type N-methyl-D-aspartate receptor (NMDAR) has been implicated in the CNS hyperexcitability of HPNS. NMDARs containing different subunits exhibited varying degrees of increased/decreased current at high pressure. The mechanisms underlying this phenomenon remain unclear. We performed 100 ns molecular dynamics (MD) simulations of the NMDAR structure embedded in a dioleoylphosphatidylcholine (DOPC) lipid bilayer solvated in water at 1 bar, hydrostatic 25 bar, and in helium at 25 bar. MD simulations showed that in contrast to hydrostatic pressure, high pressure helium causes substantial distortion of the DOPC membrane due to its accumulation between the two monolayers: reduction of the Sn-1 and Sn-2 DOPC chains and helium-dependent dehydration of the NMDAR pore. Further analysis of important regions of the NMDAR protein such as pore surface (M2 α-helix), Mg2+ binding site, and TMD-M4 α-helix revealed significant effects of helium. In contrast with previous models, these and our earlier results suggest that high pressure helium, not hydrostatic pressure per se, alters the receptor tertiary structure via protein-lipid interactions. Helium in divers' breathing mixtures may partially contribute to HPNS symptoms.

Modelling of noble anaesthetic gases and high hydrostatic pressure effects in lipid bilayers.

Our objective was to study molecular processes that might be responsible for inert gas narcosis and high-pressure nervous syndrome. The classical molecular dynamics trajectories (200 ns) of dioleoylphosphatidylcholine (DOPC) bilayers simulated by the Berger force field were evaluated for water and the atomic distribution of noble gases around DOPC molecules in the pressure range of 1-1000 bar and at a temperature of 310 K. Xenon and argon have been tested as model gases for general anaesthetics, and neon has been investigated for distortions that are potentially responsible for neurological tremors in hyperbaric conditions. The analysis of stacked radial pair distribution functions of DOPC headgroup atoms revealed the explicit solvation potential of the gas molecules, which correlates with their dimensions. The orientational dynamics of water molecules at the biomolecular interface should be considered as an influential factor, while excessive solvation effects appearing in the lumen of membrane-embedded ion channels could be a possible cause of inert gas narcosis. All the noble gases tested exhibit similar order parameter patterns for both DOPC acyl chains, which are opposite of the patterns found for the order parameter curve at high hydrostatic pressures in intact bilayers. This finding supports the 'critical volume' hypothesis of anaesthesia pressure reversal. The irregular lipid headgroup-water boundary observed in DOPC bilayers saturated with neon in the pressure range of 1-100 bar could be associated with the possible manifestation of neurological tremors at the atomic scale. The non-immobiliser neon also demonstrated the highest momentum impact on the normal component of the DOPC diffusion coefficient representing the monolayer undulation rate, which indicates that enhanced diffusivity rather than atomic size is the key factor.

Conformational changes of globular proteins upon adsorption on a hydrophobic surface.

This paper presents a study of protein adsorption and denaturation using coarse-grained Monte Carlo simulations with simulated annealing. Intermolecular interactions are modeled using the Miyazawa-Jernigan (MJ) knowledge-based potential for an implicit solvent. Three different hydrophobicity scales are tested for adsorption of fibronectin on a hydrophobic surface. The hydrophobic scale BULDG was chosen for further analysis due to its greater stability during heating and its partial regenerative ability upon slow cooling. Differences between helical and sheet structures are observed upon denaturation -α-helices undergo spreading of their native helical order to an elliptical perturbed shape, while ß-sheets transform into random coils and other more structured conformations. Electronic calculations carried out on rebuilt all-atom coordinates of adsorbed lysozymes revealed consistent destabilization of helices, while beta sheets show a greater variety of trends.

Thermal stability limits of proteins in solution and adsorbed on a hydrophobic surface.

A coarse-grained Monte Carlo simulation is used to study thermal denaturation of small proteins in an infinitely dilute solution and adsorbed on a flat hydrophobic surface. Intermolecular interactions are modeled using the Miyazawa-Jernigan (MJ) knowledge-based potential for implicit solvent with the BULDG hydrophobicity scale. We analyze the thermal behavior of lysozyme for its prevalence of α-helices, fibronectin for its prevalence of ß-sheets, and a short single helical peptide. Protein dimensions and contact maps are studied in detail before and during isothermal adsorption and heating. The MJ potential is shown to correctly predict the native conformation in solution under standard conditions, and the anticipated thermal stabilization of adsorbed proteins is observed when compared with heating in solution. The helix of the peptide is found to be much less stable thermally than the helices of lysozyme, reinforcing the importance of long-range forces in defining the protein structure. Contact map analysis of the adsorbed proteins shows correlation between the hydrophobicity of the secondary structure and their thermal stability on the surface.

Mean-field model of immobilized enzymes embedded in a grafted polymer layer.

Two-dimensional mean-field lattice theory is used to model immobilization and stabilization of an enzyme on a hydrophobic surface using grafted polymers. Although the enzyme affords biofunctionality, the grafted polymers stabilize the enzyme and impart biocompatibility. The protein is modeled as a compact hydrophobic-polar polymer, designed to have a specific bulk conformation reproducing the catalytic cleft of natural enzymes. Three scenarios are modeled that have medical or industrial importance: 1), It is shown that short hydrophilic grafted polymers, such as polyethylene glycol, which are often used to provide biocompatibility, can also serve to protect a surface-immobilized enzyme from adsorption and denaturation on a hydrophobic surface. 2), Screening of the enzyme from the surface and nonspecific interactions with biomaterial in bulk solution requires a grafted layer composed of short hydrophilic polymers and long triblock copolymers. 3), Hydrophilic polymers grafted on a hydrophobic surface in contact with an organic solvent form a dense hydrophilic nanoenvironment near the surface that effectively shields and stabilizes the enzyme against both surface and solvent.

Stabilization of surface-immobilized enzymes using grafted polymers.

We introduce a two-dimensional lattice model of immobilization and stabilization of proteinlike polymers using grafted polymers. The protein is designed to have a specific bulk conformation reproducing a catalytic cleft of natural enzymes. Our model predicts a first order denaturing adsorption transition of free proteins. On the other hand, for an immobilized protein we observe a more gradual disappearance of the hydrophobic centers accompanied by adsorption. We show that, using hydrophilic grafted polymers of proper length and grafting density, the conformation as well as the hydrophobic centers of the protein can be restored.