Binding of polar and hydrophobic molecules at the LiCoO2 (001)-water interface: force field development and molecular dynamics simulations.

A classical model in the framework of the INTERFACE force field has been developed for treating the LiCoO2 (LCO) (001)/water interface. In comparison to ab initio molecular dynamics (MD) simulations based on density functional theory, MD simulations using the classical model lead to generally reliable descriptions of interfacial properties, such as the density distribution of water molecules. Water molecules in close contact with the LCO surface form a strongly adsorbed layer, which leads to a free energy barrier for the adsorption of polar or charged molecules to the LCO surface. Moreover, due to the strong hydrogen bonding interactions with the LCO surface, the first water layer forms an interface that exhibits hydrophobic characters, leading to favorable adsorption of non-polar molecules to the interface. Therefore, despite its highly polar nature, the LCO (001) surface binds not only polar/charged but also non-polar solutes. As an application, the model is used to analyze the adsorption of reduced nicotinamide adenine dinucleotide (NADH) and its molecular components to the LCO (001) surface in water. The results suggest that recently observed redox activity of NADH at the LCO/water interface was due to the co-operativity between the ribose component, which drives binding to the LCO surface, and the nicotinamide moiety, which undergoes oxidation.

Anionic nanoparticle-induced perturbation to phospholipid membranes affects ion channel function.

Understanding the mechanisms of nanoparticle interaction with cell membranes is essential for designing materials for applications such as bioimaging and drug delivery, as well as for assessing engineered nanomaterial safety. Much attention has focused on nanoparticles that bind strongly to biological membranes or induce membrane damage, leading to adverse impacts on cells. More subtle effects on membrane function mediated via changes in biophysical properties of the phospholipid bilayer have received little study. Here, we combine electrophysiology measurements, infrared spectroscopy, and molecular dynamics simulations to obtain insight into a mode of nanoparticle-mediated modulation of membrane protein function that was previously only hinted at in prior work. Electrophysiology measurements on gramicidin A (gA) ion channels embedded in planar suspended lipid bilayers demonstrate that anionic gold nanoparticles (AuNPs) reduce channel activity and extend channel lifetimes without disrupting membrane integrity, in a manner consistent with changes in membrane mechanical properties. Vibrational spectroscopy indicates that AuNP interaction with the bilayer does not perturb the conformation of membrane-embedded gA. Molecular dynamics simulations reinforce the experimental findings, showing that anionic AuNPs do not directly interact with embedded gA channels but perturb the local properties of lipid bilayers. Our results are most consistent with a mechanism in which anionic AuNPs disrupt ion channel function in an indirect manner by altering the mechanical properties of the surrounding bilayer. Alteration of membrane mechanical properties represents a potentially important mechanism by which nanoparticles induce biological effects, as the function of many embedded membrane proteins depends on phospholipid bilayer biophysical properties.

Interfacial water and ion distribution determine ζ potential and binding affinity of nanoparticles to biomolecules.

The molecular features that dictate interactions between functionalized nanoparticles and biomolecules are not well understood. This is in part because for highly charged nanoparticles in solution, establishing a clear connection between the molecular features of surface ligands and common experimental observables such as ζ potential requires going beyond the classical models based on continuum and mean field models. Motivated by these considerations, molecular dynamics simulations are used to probe the electrostatic properties of functionalized gold nanoparticles and their interaction with a charged peptide in salt solutions. Counterions are observed to screen the bare ligand charge to a significant degree even at the moderate salt concentration of 50 mM. As a result, the apparent charge density and ζ potential are largely insensitive to the bare ligand charge densities, which fall in the range of ligand densities typically measured experimentally for gold nanoparticles. While this screening effect was predicted by classical models such as the Manning condensation theory, the magnitudes of the apparent surface charge from microscopic simulations and mean-field models are significantly different. Moreover, our simulations found that the chemical features of the surface ligand (e.g., primary vs. quaternary amines, heterogeneous ligand lengths) modulate the interfacial ion and water distributions and therefore the interfacial potential. The importance of interfacial water is further highlighted by the observation that introducing a fraction of hydrophobic ligands enhances the strength of electrostatic binding of the charged peptide. Finally, the simulations highlight that the electric double layer is perturbed upon binding interactions. As a result, it is the bare charge density rather than the apparent charge density or ζ potential that better correlates with binding affinity of the nanoparticle to a charged peptide. Overall, our study highlights the importance of molecular features of the nanoparticle/water interface and underscores a set of design rules for the modulation of electrostatic driven interactions at nano/bio interfaces.

Molecular Dynamics Simulation of Interaction between Functionalized Nanoparticles with Lipid Membranes: Analysis of Coarse-Grained Models.

We compare atomistic and two popular coarse-grained (POL- and BMW-MARTINI) models by studying the interaction between a cationic gold nanoparticle functionalized with primary alkane amines and a lipid bilayer that consists of either zwitterionic lipids or a mixture of zwitterionic and anionic lipids. In the atomistic simulations, the nanoparticle does not exhibit any notable affinity to the zwitterionic bilayer but readily binds to the 9:1 zwitterionic:anionic bilayer, and nanoparticle adsorption leads to local segregation of anionic lipids and slowing down of their diffusion. At the coarse-grained level, the POL-MARTINI model does not lead to nanoparticle-membrane binding for either bilayer system, while the BMW-MARTINI model leads to nanoparticle binding to both bilayers; with the BMW-MARTINI model, nanoparticle binding leads to much less demixing of zwitterionic and anionic lipids and moderately higher rates of lipid diffusion. Analysis of nanoparticle properties in solution reveals notable discrepancies in the interfacial charge and water distributions at the coarse-grained level that are likely relevant to their limitations in describing binding interactions with other (bio)molecules.

Counting charges on membrane-bound peptides.

Quantifying the number of charges on peptides bound to interfaces requires reliable estimates of (i) surface coverage and (ii) surface charge, both of which are notoriously difficult parameters to obtain, especially at solid/water interfaces. Here, we report the thermodynamics and electrostatics governing the interactions of l-lysine and l-arginine octamers (Lys8 and Arg8) with supported lipid bilayers prepared from a 9 : 1 mixture of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dimyristoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (sodium salt) (DMPG) from second harmonic generation (SHG) spectroscopy, quartz crystal microbalance with dissipation monitoring (QCM-D) and nanoplasmonic sensing (NPS) mass measurements, and atomistic simulations. The combined SHG/QCM-D/NPS approach provides interfacial charge density estimates from mean field theory for the attached peptides that are smaller by a factor of approximately two (0.12 ± 0.03 C m-2 for Lys8 and 0.10 ± 0.02 C m-2 for Arg8) relative to poly-l-lysine and poly-l-arginine. These results, along with atomistic simulations, indicate that the surface charge density of the supported lipid bilayer is neutralized by the attached cationic peptides. Moreover, the number of charges associated with each attached peptide is commensurate with those found in solution; that is, Lys8 and Arg8 are fully ionized when attached to the bilayer. Computer simulations indicate Lys8 is more likely than Arg8 to "stand-up" on the surface, interacting with lipid headgroups through one or two sidechains while Arg8 is more likely to assume a "buried" conformation, interacting with the bilayer through up to six sidechains. Analysis of electrostatic potential and charge distribution from atomistic simulations suggests that the Gouy-Chapman model, which is widely used for mapping surface potential to surface charge, is semi-quantitatively valid; despite considerable orientational preference of interfacial water, the apparent dielectric constant for the interfacial solvent is about 30, due to the thermal fluctuation of the lipid-water interface.

Hydrogen-Bond Networks near Supported Lipid Bilayers from Vibrational Sum Frequency Generation Experiments and Atomistic Simulations.

We report vibrational sum frequency generation (SFG) spectra in which the C-H stretches of lipid alkyl tails in fully hydrogenated single- and dual-component supported lipid bilayers are detected along with the O-H stretching continuum above the bilayer. As the salt concentration is increased from ∼10 µM to 0.1 M, the SFG intensities in the O-H stretching region decrease by a factor of 2, consistent with significant absorptive-dispersive mixing between χ(2) and χ(3) contributions to the SFG signal generation process from charged interfaces. A method for estimating the surface potential from the second-order spectral lineshapes (in the OH stretching region) is presented and discussed in the context of choosing truly zero-potential reference states. Aided by atomistic simulations, we find that the strength and orientation distribution of the hydrogen bonds over the purely zwitterionic bilayers are largely invariant between submicromolar and hundreds of millimolar concentrations. However, specific interactions between water molecules and lipid headgroups are observed upon replacing phosphocholine (PC) lipids with negatively charged phosphoglycerol (PG) lipids, which coincides with SFG signal intensity reductions in the 3100-3200 cm-1 frequency region. The atomistic simulations show that this outcome is consistent with a small, albeit statistically significant, decrease in the number of water molecules adjacent to both the lipid phosphate and choline moieties per unit area, supporting the SFG observations. Ultimately, the ability to probe hydrogen-bond networks over lipid bilayers holds the promise of opening paths for understanding, controlling, and predicting specific and nonspecific interactions between membranes and ions, small molecules, peptides, polycations, proteins, and coated and uncoated nanomaterials.

Analysis of the conformational properties of amine ligands at the gold/water interface with QM, MM and QM/MM simulations.

We describe a strategy of integrating quantum mechanical (QM), hybrid quantum mechanical/molecular mechanical (QM/MM) and MM simulations to analyze the physical properties of a solid/water interface. This protocol involves using a correlated ab initio (CCSD(T)) method to first calibrate Density Functional Theory (DFT) as the QM approach, which is then used in QM/MM simulations to compute relevant free energy quantities at the solid/water interface using a mean-field approximation of Yang et al. that decouples QM and MM thermal fluctuations; gas-phase QM/MM and periodic DFT calculations are used to determine the proper QM size in the QM/MM simulations. Finally, the QM/MM free energy results are compared with those obtained from MM simulations to directly calibrate the force field model for the solid/water interface. This protocol is illustrated by examining the orientations of an alkyl amine ligand at the gold/water interface, since the ligand conformation is expected to impact the chemical properties (e.g., charge) of the solid surface. DFT/MM and MM simulations using the INTERFACE force field lead to consistent results, suggesting that the effective gold/ligand interactions can be adequately described by a van der Waals model, while electrostatic and induction effects are largely quenched by solvation. The observed differences among periodic DFT, QM/MM and MM simulations, nevertheless, suggest that explicitly including electronic polarization and potentially charge transfer in the MM model can be important to the quantitative accuracy. The strategy of integrating multiple computational methods to cross-validate each other for complex interfaces is applicable to many problems that involve both inorganic/metallic and organic/biomolecular components, such as functionalized nanoparticles.

Quantifying the Electrostatics of Polycation-Lipid Bilayer Interactions.

Mechanistic insight into how polycations disrupt and cross cell membranes is needed for understanding and controlling polycation-membrane interactions, yet such information is surprisingly difficult to obtain at the molecular level. We use second harmonic and vibrational sum frequency generation spectroscopies along with quartz crystal microbalance with dissipation monitoring and computer simulations to quantify the interaction of poly(allylamine) hydrochloride (PAH) and its monomeric precursor allylamine hydrochloride (AH) with lipid bilayers. We find PAH adsorption to be reversible and nondisruptive to the bilayer under the conditions of our experiments. With an observed free adsorption energy of -52.7 ± 0.6 kJ/mol, PAH adsorption was found to be surprisingly less favorable relative to AH (-14.6 ± 0.4 kJ/mol) when considering a simple additive model. By experimentally quantifying the number of adsorbates and the average amount of charge carried by each adsorbate, we find that the PAH is associated with only 70% of the positive charges it could hold while the AH remains mostly charged while attached to the membrane. Simulations indicate that PAH pulls in condensed counterions from solution to avoid charge-repulsion along its backbone and with other PAH molecules to attach to, and completely cover, the bilayer surface. In addition, computations indicate that the amine groups shift their pKa values due to the confined environment upon adsorption to the surface. Our results provide experimental constraints for theoretical calculations, which yield atomistic views of the structures that are formed when polycations interact with lipid membranes that will be important for predicting polycation-membrane interactions.

[Acupotomy and acupuncture in the treatment of avascular necrosis of femoral head at the early and middle stages:a clinical randomized controlled trial].

OBJECTIVE: To compare the efficacy difference between acupotomy and acupuncture in the treatment of avascular necrosis of femoral head at the early and middle stages. METHODS: The randomized controlled prospective study method was adopted. Sixty cases of avascular necrosis of femoral head at Ficat-ArletⅠto Ⅱ stages were randomized into an acupotomy group (32 cases) and an acupuncture group (28 cases) by the third part. In the acupotomy group, the acupotomy was adopted for the loose solution at the treatment sites of hip joint, once every two weeks, totally for 3 times. In the acupuncture group, ashi points around the hip joint were selected and stimulated with warm acupuncture therapy, once every day, for 6 weeks. Harris hip score was observed before and after treatment. The efficacy was evaluated in the two groups. RESULTS: Harris hip score was improved significantly after treatment in the two groups (both P<0.05). The result in acupotomy group was better than that in the acupuncture group (P<0.05). The effective rate was 90.6% (29/32) in the acupotomy group, better than 75.0% (21/28) in the acupuncture group after treatment (P<0.05). CONCLUSIONS: Harris hip score and the effective rate in the acupotomy group are better than those in the treatment with routine acupuncture for avascular necrosis of femoral head at the early and middle stages.