Search

Free fatty acids in biomembranes have been proposed to be a central component in several cellular control and regulatory mechanisms. To elucidate some fundamental elements underlying this, we have applied molecular dynamics simulations and experimental density measurements to study the molecular packing and structure of oleic acid (HOA) and stearic acid (HSA) in fluid bilayers of dimyristoylphosphatidylcholine (DMPC). The experimental data show a small but consistent positive excess volume for fatty acid concentrations below 10 mol %. At higher concentrations the fatty acids mix ideally with fluid DMPC. The simulations, which were benchmarked against the densitometric data, revealed interesting differences in the structure and location of the fatty acids depending on their protonation status. Thus, the protonated (uncharged) acid is located rather deeply in the membrane with an average position of the carboxy group near the second carbon segment of the lipid chains with a typical end-to-end distance of 16−18 Å. This structure of the fatty acid brings about a rather tight lateral packing in the mixed membrane and a moderate ordering and hence stretching of the lipid chains. Deprotonation of the fatty acids is associated with a pronounced movement of their carboxy group to a more hydrated position at the membrane interface and a lateral expansion driven by the mutual repulsion of the anions. These changes increase both the disorder and the degree of interdigitation of the lipid chains, and they make the membrane thinner by 2−3 Å. [ABSTRACT FROM AUTHOR]

To study the flexibility of the substrate-binding site and in particular of Gln262, we have performed adiabatic conformational search and molecular dynamics simulations on the crystal structure of the catalytic domain of wild-type protein-tyrosine phosphatase (PTP) 1 B, a mutant PTP1BR47v,D48N,M258c,G259Q, and a model of the catalytically active form of PTPα. For each molecule two cases were modeled: the Michaelis-Menten complex with the substrate analogue p-nitrophenyl phosphate (p-PNPP) bound to the active site and the cysteine-phosphor complex, each corresponding to the first and second step of the phosphate hydrolysis. Analyses of the trajectories revealed that in the cysteine-phosphor complex of PTP1B, G1n262 oscillates freely between the bound phosphate group and G1y259 frequently forming, as observed in the crystal structure, a hydrogen bond with the backbone oxygen of G1y259. In contrast, the movement of G1n262 is restricted in PTPα and the mutant due to interactions with G1n259 reducing the frequency of the oscillation of G1n262 and thereby delaying the positioning of this residue for the second step in the catalysis, as reflected experimentally by a reduction in kcat. Additionally, in the simulation with the Michaelis-Menten complexes, we found that a glutamine in position 259 induces steric hindrance by pushing the G1n262 side chain further toward the substrate and thereby negatively affecting Km as indicated by kinetic studies. Detailed analysis of the water structure around G1n262 and the active site Cys215 reveals that the probability of finding a water molecule correctly positioned for catalysis is much larger in PTP1B than in PTP1BR47V,D48N,M258c,G259Q and PTPα, in accordance with experiments. [ABSTRACT FROM AUTHOR]

Secretory phospholipase A2 (sPLA2) is an interesting enzyme for triggered liposomal drug delivery to tumor tissue due the overexpression of sPLA2 in cancerous tissue. A drug delivery system based on the triggered release of therapeutics from sPLA2-sensitive liposomes constituted of pro anticancer ether lipids, which become cytotoxic upon sPLA2-catalyzed hydrolysis has previously been established. To optimize the hydrolysis rate of the lipids and thereby optimizing the release profile of the drugs from the liposomes, we have synthesized a thio-ester pro anticancer ether lipid. Liposomes constituted of this lipid showed an altered rate of hydrolysis by sPLA2. We have tested the cytotoxicity of the thio-ester pro anticancer ether lipids toward cancer cells, and the results showed that the cytotoxicity is indeed maintained upon sPLA2 exposure. To further understand the origin for the observed different hydrolysis rates for the esters, we have applied molecular dynamics simulations and density functional theory. The combination of these theoretical methods has given valuable insight into the molecular mechanism for sPLA2 action on sulfur-containing phospholipids. It appears that the enzyme-catalyzed hydrolysis of thio-esters follow a different pathway compared to the hydrolysis pathway of the free thio-ester. [ABSTRACT FROM AUTHOR]

Therapeutic proteins can be challenging to develop due to their complexity and the requirement of an acceptable formulation to ensure patient safety and efficacy. To date, there is no universal formulation development strategy that can identify optimal formulation conditions for all types of proteins in a fast and reliable manner. In this work, high-throughput characterization, employing a toolbox of five techniques, was performed on 14 structurally different proteins formulated in 6 different buffer conditions and in the presence of 4 different excipients. Multivariate data analysis and chemometrics were used to analyze the data in an unbiased way. First, observed changes in stability were primarily determined by the individual protein. Second, pH and ionic strength are the two most important factors determining the physical stability of proteins, where there exists a significant statistical interaction between protein and pH/ionic strength. Additionally, we developed prediction methods by partial least-squares regression. Colloidal stability indicators are important for prediction of real-time stability, while conformational stability indicators are important for prediction of stability under accelerated stress conditions at 40 °C. In order to predict real-time storage stability, protein-protein repulsion and the initial monomer fraction are the most important properties to monitor.

Secretory phospholipases A 2 (sPLA 2 s) are a subclass of enzymes that catalyze the hydrolysis at the sn-2 position of glycerophospholipids, producing free fatty acids and lysophospholipids. In this study, different phospholipids with structural modifications close to the scissile sn-2 ester bond were studied to determine the effect of the structural changes on the formation of the Michaelis-Menten complex and the water entry/exit pathways using molecular dynamics simulations and the computational tracking tool AQUA-DUCT. Structural modifications include methylation, dehydrogenation, and polarization close to the sn-2 scissile bond. We found that all water molecules reaching the active site of sPLA 2 -IIA pass by the aromatic residues Phe 5 and Tyr 51 and enter the active site through an active-site cleft. The relative amount of water available for the enzymatic reaction of the different phospholipid-sPLA 2 complexes was determined together with the distance between key atoms in the catalytic machinery. The results showed that ( Z )-unsaturated phospholipid is a good substrate for sPLA 2 -IIA. The computational results are in good agreement with previously reported experimental data on the ability of sPLA 2 -IIA to hydrolyze liposomes made from the different phospholipids, and the results provide new insights into the necessary active-site solvation of the Michaelis-Menten complex and can pave the road for rational design in engineering applications.

Water has a critical role in the stability of the higher-order structure of proteins. In addition, it is considered to be a major destabilization factor for the physical and chemical stability of freeze-dried proteins and peptides. Physical and chemical aspects of protein/water relationships are commonly studied with the use of water vapor sorption isotherms for amorphous lyophilized proteins, which, in turn, are commonly analyzed using the Brunauer-Emmett-Teller (BET) equation to obtain the parameters, W m and C B . The parameter W m is generally referred to as the "monolayer limit of adsorption" and has a narrow range of 6-8% for most proteins. In this study, the water distribution on an IgG1 surface is investigated by molecular dynamics (MD) simulations at different water contents. The monolayer of water molecules was found to have limited coverage of the protein surface, and the true monolayer coverage of the protein globule actually occurs at a hydration level above 30%. The distribution of water molecules on the IgG1 surface is also highly heterogeneous, and the heterogeneity is not considered in the BET theory. In this study, a mechanistic model has been developed to describe the water vapor sorption isotherm. This model is based on the analysis of the hydrogen bonding network extracted from the MD simulations. The model is consistent with the experimental Type-II isotherm, which is usually observed for proteins. The physical meaning of the BET monolayer was redefined as the onset of water cluster formation. A simple model to calculate the onset water level, W m , is proposed based on the hydration of different amino acids, as determined from the MD simulations.

Therapeutic protein candidates should exhibit favorable properties that render them suitable to become drugs. Nevertheless, there are no well-established guidelines for the efficient selection of proteinaceous molecules with desired features during early stage development. Such guidelines can emerge only from a large body of published research that employs orthogonal techniques to characterize therapeutic proteins in different formulations. In this work, we share a study on a diverse group of proteins, including their primary sequences, purity data, and computational and biophysical characterization at different pH and ionic strength. We report weak linear correlations between many of the biophysical parameters. We suggest that a stability comparison of diverse therapeutic protein candidates should be based on a computational and biophysical characterization in multiple formulation conditions, as the latter can largely determine whether a protein is above or below a certain stability threshold. We use the presented data set to calculate several stability risk scores obtained with an increasing level of analytical effort and show how they correlate with protein aggregation during storage. Our work highlights the importance of developing combined risk scores that can be used for early stage developability assessment. We suggest that such scores can have high prediction accuracy only when they are based on protein stability characterization in different solution conditions.

Plectasin is a small, cysteine-rich peptide antibiotic which belongs to the class of antimicrobial peptides and has potential antibacterial activity against various Gram-positive bacteria. In the current study, the effect of pH and ionic strength (NaCl) on the conformational stability of plectasin variants has been investigated. At all physiochemical conditions, peptide secondary structures are intact throughout simulations. However, flexibility increases with pH because of the change in electrostatics around the distinct anionic tetrapeptide (9-12) stretch. Furthermore, plectasin interactions with NaCl were measured by determining the preferential interaction coefficients, Γ 23 . Generally, wild-type plectasin has higher preference for sodium ions as 9ASP is mutated in other variants. Overall, the Γ 23 trend with pH for the two salt conditions remain the same for all variants predominately having accumulation of sodium ions around 10GLU and 12ASP. Insignificant changes in the overall peptide conformational stability are in agreement with the fact that plectasin has three cystines. Thermodynamic integration molecular dynamics simulations supplemented with nuclear magnetic resonance were employed to determine the degree of involvement of three different cystines to the overall structural integrity of the peptide. Both methods show the same order of cystine reduction and complete unfolding is observed only upon reduction of all cystines.

In the present study, small-angle X-ray scattering (SAXS) and static light scattering (SLS) have been used to study the solution properties and self-interaction of recombinant human serum albumin (rHSA) molecules in three pharmaceutically relevant buffer systems. Measurements are carried out up to high protein concentrations and as a function of ionic strength by adding sodium chloride to probe the role of electrostatic interactions. The effective structure factors ( S eff ) as a function of the scattering vector magnitude q have been extracted from the scattering profiles and fit to the solution of the Ornstein-Zernike equation using a screened Yukawa potential to describe the double-layer force. Although only a limited q range is used, accurate fits required including an electrostatic repulsion element in the model at low ionic strength, while only a hard sphere model with a tunable diameter is necessary for fitting to high-ionic-strength data. The fit values of net charge agree with available data from potentiometric titrations. Osmotic compressibility data obtained by extrapolating the SAXS profiles or directly from SLS measurements has been fit to a 10-term virial expansion for hard spheres and an equation of state for hard biaxial ellipsoids. We show that modeling rHSA as an ellipsoid, rather than a sphere, provides a much more accurate fit for the thermodynamic data over the entire concentration range. Osmotic virial coefficient data, derived at low protein concentration, can be used to parameterize the model for predicting the behavior up to concentrations as high as 450 g/L. The findings are especially important for the biopharmaceutical sector, which require approaches for predicting concentrated protein solution behavior using minimal sample consumption., Competing Interests: The authors declare no competing financial interest.

Membrane mechanical elastic properties regulate a variety of cellular processes involving local membrane deformation, such as ion channel function and vesicle fusion. In this work, we used molecular dynamics simulations to estimate the local elastic properties of a membrane. For this, we calculated the energy needed to extract a DOPE lipid molecule, modified with a linker chain, from a POPC bilayer membrane using the umbrella sampling technique. Although the extraction energy entails several contributions related not only to elastic deformation but also to solvation, careful analysis of the potential of mean force (PMF) allowed us to dissect the elastic contribution. With this information, we calculated an effective linear spring constant of 44 ± 4 kJ·nm(-2)·mol(-1) for the DOPC membrane, in agreement with experimental estimates. The membrane deformation profile was determined independently during the stretching process in molecular detail, allowing us to fit this profile to a previously proposed continuum elastic model. Through this approach, we calculated an effective membrane spring constant of 42 kJ·nm(-2)·mol(-1), which is in good agreement with the PMF calculation. Furthermore, the solvation energy we derived from the data is shown to match the solvation energy estimated from critical micelle formation constants. This methodology can be used to determine how changes in lipid composition or the presence of membrane modifiers can affect the elastic properties of a membrane at a local level.

The inclusion complexes of glycoconjugated bile salts with β-cyclodextrin (β-CD) and 2-hydroxypropyl-β-cyclodextrins (HP-β-CD) in aqueous solution were investigated by molecular dynamics simulations to provide a molecular explanation of the experimentally observed destabilizing effect of the HP substituents. Good agreement with experimental data was found with respect to penetration depths of CDs. An increased degree of HP substitution (DS) resulted in an increased probability of blocking the cavity opening, thereby hindering the bile salt from entering CD. Further, the residence time of water molecules in the cavity increased with the DS. Release of water from the cavity resulted in a positive enthalpy change, which correlates qualitatively with the experimentally determined increase in complexation enthalpy and contributes to the enthalpy-entropy compensation. The positive change in complexation entropy with DS was not able to compensate for this unfavorable change in enthalpy induced by the HP substituents, resulting in a destabilizing effect. This was found to originate from fixation of the HP substituents and decreased free rotation of the bile salts within the CD cavities.

Serotonin (5-hydroxytryptamine, 5-HT) is a prevalent neurotransmitter throughout the animal kingdom. It exerts its effect through the specific binding to the serotonin receptor, but recent research has suggested that neural transmission may also be affected by its nonspecific interactions with the lipid matrix of the synaptic membrane. However, membrane-5-HT interactions remain controversial and superficially investigated. Fundamental knowledge of this interaction appears vital in discussions of putative roles of 5-HT, and we have addressed this by thermodynamic measurements and molecular dynamics (MD) simulations. 5-HT was found to interact strongly with lipid bilayers (partitioning coefficient ~1200 in mole fraction units), and this is highly unusual for a hydrophilic solute like 5-HT which has a bulk, oil-water partitioning coefficient well below unity. It follows that membrane affinity must rely on specific interactions, and the MD simulations identified the salt-bridge between the primary amine of 5-HT and the lipid phosphate group as the most important interaction. This interaction anchored cationic 5-HT in the membrane interface with the aromatic ring system pointing inward and a prevailing residence between the phosphate and the carbonyl groups of the lipid. The unprotonated form of 5-HT shows the opposite orientation, with the primary amine pointing toward the membrane core. Partitioning of 5-HT was found to decrease lipid chain order. These distinctive interactions of 5-HT and model membranes could be related to nonspecific effects of this neurotransmitter.

In nonaqueous enzymology, control of enzyme hydration is commonly approached by fixing the thermodynamic water activity of the medium. In this work, we present a strategy for evaluating the water activity in molecular dynamics simulations of proteins in water/organic solvent mixtures. The method relies on determining the water content of the bulk phase and uses a combination of Kirkwood-Buff theory and free energy calculations to determine corresponding activity coefficients. We apply the method in a molecular dynamics study of Candida antarctica lipase B in pure water and the organic solvents methanol, tert-butyl alcohol, methyl tert-butyl ether, and hexane, each mixture at five different water activities. It is shown that similar water activity yields similar enzyme hydration in the different solvents. However, both solvent and water activity are shown to have profound effects on enzyme structure and flexibility.

We have studied secretory phospholipase A(2)-IIA (sPLA(2)) activity toward different phospholipid analogues by performing biophysical characterizations and molecular dynamics simulations. The phospholipids were natural substrates, triple alkyl phospholipids, a prodrug anticancer etherlipid, and an inverted ester. The latter were included to study head group-enzyme interactions. Our simulation results show that the lipids are optimally placed into the binding cleft and that water molecules can freely reach the active site through a well-defined pathway; both are indicative that these substrates are efficiently hydrolyzed, which is in good agreement with our experimental data. The phospholipid analogue with three alkyl side chains forms aggregates of different shapes with no well-defined sizes due to its cone-shape structure. Phosphatidylglycerol and phosphatidylcholine head groups interact with specific charged residues, but relatively large fluctuations are observed, suggesting that these interactions are not necessarily important for stabilizing substrate binding to the enzyme.

Weak interactions of neurotransmitters and the lipid matrix in the synaptic membrane have been hypothesized to play a role in synaptic transmission of nerve signals, particularly with respect to receptor desensitization (Cantor, R. S. Biochemistry 2003, 42, 11891). The strength of such interactions, however, was not measured, and this is an obvious impediment for further evaluation and understanding of a possible role for desensitization. We have used dialysis equilibrium to directly measure the net affinity of selected neurotransmitters for lipid membranes and analyzed this affinity data with respect to calorimetric measurements and molecular dynamics simulations. We studied an anionic (glutamate), a cationic (acetylcholine), and two zwitterionic (γ-aminobutyric acid and glycine) neurotransmitters, and membranes of pure dimyristoyl phosphatidylcholine (DMPC), DMPC doped with 10% anionic lipid (dimyristoyl phosphatidylglycerol, DMPG, or dimyristoyl phosphatidylserine, DMPS), or 1:1 mixtures of dipalmitoyl phosphatidylcholine (DPPC) and dilauroyl phosphatidylcholine (DLPC). The results showed a remarkable variability among the investigated systems. For example, the chloride salt of acetylcholine interacts unfavorably with DMPC and is thus preferentially excluded from the membrane's hydration layer. Conversely, the zwitterionic neurotransmitters are attracted to membranes with 10% anionic lipid and their local concentration at the interface is 5-10 times larger than in the aqueous bulk. The simulations suggest that this attraction mainly relies on electrostatic interactions of the amino group of the neurotransmitter and the lipid phosphate. We conclude that moderate attraction to lipid membranes occurs for some polar neurotransmitters and hence that one premise for a theory of bilayer-mediated modulation of nerve transmission seems to be fulfilled. However, the strong variability in interaction strengths also shows that this attraction is not an inherent property of all neurotransmitters.

Although the Escherichia coli ammonia transporter B (AmtB) protein has been the focus of several recent studies, there are still many questions and controversies regarding substrate binding and recognition. Specifically, how and where AmtB differentiates between substrates is not yet fully understood. The present computational study addresses the importance of intermolecular interactions with respect to substrate recruitment and recognition by means of ab initio QM/MM simulations. On the basis of calculations with substrates NH(3), NH(4)(+), Na(+), and K(+) positioned at the periplasmic binding site (Am1) and NH(3) and NH(4)(+) at intraluminal binding sites (Am1a/b), we conclude that D160 is the single most important residue for substrate recruitment, whereas cation-π interactions to W148 and F107 are found to be less important. Regarding substrate recruitment and recognition, we find that only NH(4)(+) and K(+) reach the Am1 site. However, NH(4)(+) has the largest affinity for this site due to its better dehydration compensation, while charge stabilization effects favor the binding of NH(4)(+) over NH(3) (i.e., if NH(3) would enter the Am1 site, it is likely to be protonated). Therefore, we conclude that the Am1 site selects NH(4)(+) over Na(+), K(+) and NH(3). Our calculations also suggest that translocation of NH(4)(+) from Am1 into the channel lumen is driven by rotation of the A162-G163 peptide bond, which coordinates NH(4)(+) but not NH(3) at both Am1 and Am1a/b sites.

The design of retinoid phospholipid prodrugs is described based on molecular dynamics simulations and cytotoxicity studies of synthetic retinoid esters. The prodrugs are degradable by secretory phospholipase A(2) IIA and have potential in liposomal drug delivery targeting tumors. We have synthesized four different retinoid phospholipid prodrugs and shown that they form particles in the liposome size region with average diameters of 94-118 nm. Upon subjection to phospholipase A(2), the lipid prodrugs were hydrolyzed, releasing cytotoxic retinoids and lysolipids. The formulated lipid prodrugs displayed IC(50) values in the range of 3-19 microM toward HT-29 and Colo205 colon cancer cells in the presence of phospholipase A(2), while no significant cell death was observed in the absence of the enzyme.

This paper reports all-atom computer simulations of five phospholipid membranes (DMPC, DPPC, DMPG, DMPS, and DMPSH) with focus on the thermal equilibrium fluctuations of volume, energy, area, thickness, and chain order. At constant temperature and pressure, volume and energy exhibit strong correlations of their slow fluctuations (defined by averaging over 0.5 ns). These quantities, on the other hand, do not correlate significantly with area, thickness, or chain order. The correlations are mainly reported for the fluid phase, but we also give some results for the ordered (gel) phase of two membranes, showing a similar picture. The cause of the observed strong correlations is identified by splitting volume and energy into contributions from tails, heads, and water, and showing that the slow volume-energy fluctuations derive from van der Waals interactions of the tail region; they are thus analogous to the similar strong correlations recently observed in computer simulations of the Lennard-Jones and other simple van der Waals type liquids (U. R. Pedersen et al., Phys. Rev. Lett. 2008, 100, 015701). The strong correlations reported here confirm one crucial assumption of a recent theory for nerve signal propagation proposed by Heimburg and Jackson (T. Heimburg and A. D. Jackson, Proc. Natl. Acad. Sci. 2005, 102, 9790-9795).

Free fatty acids in biomembranes have been proposed to be a central component in several cellular control and regulatory mechanisms. To elucidate some fundamental elements underlying this, we have applied molecular dynamics simulations and experimental density measurements to study the molecular packing and structure of oleic acid (HOA) and stearic acid (HSA) in fluid bilayers of dimyristoylphosphatidylcholine (DMPC). The experimental data show a small but consistent positive excess volume for fatty acid concentrations below 10 mol %. At higher concentrations the fatty acids mix ideally with fluid DMPC. The simulations, which were benchmarked against the densitometric data, revealed interesting differences in the structure and location of the fatty acids depending on their protonation status. Thus, the protonated (uncharged) acid is located rather deeply in the membrane with an average position of the carboxy group near the second carbon segment of the lipid chains with a typical end-to-end distance of 16-18 A. This structure of the fatty acid brings about a rather tight lateral packing in the mixed membrane and a moderate ordering and hence stretching of the lipid chains. Deprotonation of the fatty acids is associated with a pronounced movement of their carboxy group to a more hydrated position at the membrane interface and a lateral expansion driven by the mutual repulsion of the anions. These changes increase both the disorder and the degree of interdigitation of the lipid chains, and they make the membrane thinner by 2-3 A.

A new and unnatural type of phospholipids with the head group attached to the 2-position of the glycerol backbone has been synthesized and shown to be a good substrate for secretory phospholipase A2 (sPLA2). To investigate the unexpected sPLA2 activity, we have compared three different phospholipids by using fluorescence techniques and HPLC, namely: (R)-1,2-dipalmitoyl-glycero-3-phosphocholine (hereafter referred to as 1R), (R)-1-O-hexadecyl-2-palmitoyl-glycero-3-phoshocholine (2R), and (S)-1-O-hexadecyl-3-palmitoyl-glycero-2-phosphocholine (3S). Furthermore, to understand the underlying mechanisms for the observed differences, we have performed molecular dynamics simulations to clarify on a structural level the substrate specificity of sPLA2 toward phospholipid analogues with their head groups in the 2-position of the glycerol backbone. We have studied the lipids above 1R, 2R, and 3S as well as their enantiomers 1S, 2S, and 3R. In the simulations of sPLA2-1S and sPLA2-3R, structural distortion in the binding cleft induced by the phospholipids showed that these are not substrates for sPLA2. In the case of the phospholipids 1R, 2R, and 3S, our simulations revealed that the difference observed experimentally in sPLA2 activity might be caused by reduced access of water molecules to the active site. We have monitored the number of water molecules that enter the active site region for the different sPLA2-phospholipid complexes and found that the probability of a water molecule reaching the correct position such that hydrolysis can occur is reduced for the unnatural lipids. The relative water count follows 1R > 2R > 3S. This is in good agreement with experimental data that indicate the same trend for sPLA2 activity: 1R > 2R > 3S.