A Permeability Study of O2 and the Trace Amine p-Tyramine through Model Phosphatidylcholine Bilayers.

We study here the permeability of the hydrophobic O2 molecule through a model DPPC bilayer at 323K and 350K, and of the trace amine p-tyramine through PC bilayers at 310K. The tyramine results are compared to previous experimental work at 298K. Nonequilibrium work methods were used in conjunction to simultaneously obtain both the potential of mean force (PMF) and the position dependent transmembrane diffusion coefficient, D(z), from the simulations. These in turn were used to calculate the permeability coefficient, P, through the inhomogeneous solubility-diffusion model. The results for O2 are consistent with previous simulations, and agree with experimentally measured P values for PC bilayers. A temperature dependence in the permeability of O2 through DPPC was obtained, with P decreasing at higher temperatures. Two relevant species of p-tyramine were simulated, from which the PMF and D(z) were calculated. The charged species had a large energetic barrier to crossing the bilayer of ~ 21 kcal/mol, while the uncharged, deprotonated species had a much lower barrier of ~ 7 kcal/mol. The effective in silico permeability for p-tyramine was calculated by applying three approximations, all of which gave nearly identical results (presented here as a function of the pKa). As the permeability value calculated from simulation was highly dependent on the pKa of the amine group, a further pKa study was performed that also varied the fraction of the uncharged and zwitterionic p-tyramine species. Using the experimental P value together with the simulated results, we were able to label the phenolic group as responsible for the pKa1 and the amine for the pKa2, that together represent all of the experimentally measured pKa values for p-tyramine. This agrees with older experimental results, in contrast to more recent work that has suggested there is a strong ambiguity in the pKa values.

Structural basis for autoinhibition of CTP:phosphocholine cytidylyltransferase (CCT), the regulatory enzyme in phosphatidylcholine synthesis, by its membrane-binding amphipathic helix.

CTP:phosphocholine cytidylyltransferase (CCT) interconverts between an inactive soluble and active membrane-bound form in response to changes in membrane lipid composition. Activation involves disruption of an inhibitory interaction between the αE helices at the base of the active site and an autoinhibitory (AI) segment in the regulatory M domain and membrane insertion of the M domain as an amphipathic helix. We show that in the CCT soluble form the AI segment functions to suppress kcat and elevate the Km for CTP. The crystal structure of a CCT dimer composed of the catalytic and AI segments reveals an AI-αE interaction as a cluster of four amphipathic helices (two αE and two AI helices) at the base of the active sites. This interaction corroborates mutagenesis implicating multiple hydrophobic residues within the AI segment that contribute to its silencing function. The AI-αE interaction directs the turn at the C-terminal end of the AI helix into backbone-to-backbone contact with a loop (L2) at the opening to the active site, which houses the key catalytic residue, lysine 122. Molecular dynamics simulations suggest that lysine 122 side-chain orientations are constrained by contacts with the AI helix-turn, which could obstruct its engagement with substrates. This work deciphers how the CCT regulatory amphipathic helix functions as a silencing device.

Calculating diffusion and permeability coefficients with the oscillating forward-reverse method.

The forward-reverse or FR method is an efficient bidirectional work method for determining the potential of mean force w(z) and also supposedly gives in principle the position-dependent diffusion coefficient D(z). Results from a variation called the OFR (oscillating FR) method suggest inconsistencies in the D(z) values when calculated as prescribed by the FR method. A new steering protocol has thus been developed and applied to the OFR method for the accurate determination of D(z) and also provides greater convergence for w(z) in molecular dynamics simulations. The bulk diffusion coefficient for water was found to be (6.03±0.16)×10(-5) cm2/s at 350 K with system size dependence within the statistical error bars. Using this steering protocol, D(z) and w(z) for water permeating a dipalmitoylphosphatidylcholine (DPPC) bilayer were determined. The potential of mean force is shown to have a barrier of peak height, wmax/(kBT)=8.4, with a width of about 10 Å on either side from the membrane center. The diffusion constant is shown to be highest in the core region of the membrane [peak value ∼(8.0±0.8)×10(-5) cm2/s], lowest in the head-group region [minimum value ∼(2.0±0.3)×10(-5) cm2/s], and to tend toward the bulk value as the water molecule leaves the membrane. The permeability coefficient P for H2O in DPPC was determined using the simulated D(z) and w(z) to give values of (0.129±0.075) cm/s at 323 K and (0.141±0.043) cm/s at 350 K. The results show more spatial detail than results presented in previous work while reducing the computational and user effort.

Scattering density profile model of POPG bilayers as determined by molecular dynamics simulations and small-angle neutron and X-ray scattering experiments.

We combine molecular dynamics (MD) simulations and experiment, both small-angle neutron (SANS) and small-angle X-ray scattering (SAXS), to determine the precise structure of bilayers composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG), a lipid commonly encountered in bacterial membranes. Experiment and simulation are used to develop a one-dimensional scattering density profile (SDP) model suitable for the analysis of experimental data. The joint refinement of such data (i.e., SANS and SAXS) results in the area per lipid that is then used in the fixed-area simulations. In the final step, the direct comparison of simulated-to-experimental data gives rise to the detailed structure of POPG bilayers. From these studies we conclude that POPG's molecular area is 66.0 ± 1.3 Å(2), its overall bilayer thickness is 36.7 ± 0.7 Å, and its hydrocarbon region thickness is 27.9 ± 0.6 Å, assuming a simulated value of 1203 Å(3) for the total lipid volume.

Drift-oscillatory steering with the forward-reverse method for calculating the potential of mean force.

We present a method that enables the use of the forward-reverse (FR) method of Kosztin et al. on a broader range of problems in soft matter physics. Our method, which we call the oscillating forward-reverse (OFR) method, adds an oscillatory steering potential to the constant velocity steering potential of the FR method. This enables the calculation of the potential of mean force (PMF) in a single unidirectional oscillatory drift, rather than multiple drifts in both directions as required by the FR method. By following small forward perturbations with small reverse perturbations, the OFR method is able to generate a piecewise reverse path that follows the piecewise forward path much more closely than any practical set of paths used in the FR method. We calculate the PMF for four different systems: a dragged Brownian oscillator, a pair of atoms in a Lennard-Jones liquid, a Na(+)-Cl⁻ ion pair in an aqueous solution, and a deca-alanine molecule being stretched in an implicit solvent. In all cases, the PMF results are in good agreement with those published previously using various other methods, and, to our knowledge, we give for the first time PMFs calculated by nonequilibrium methods for the Lennard-Jones and Na(+)-Cl⁻ systems.