Search

The spinach plant aquaporin SoPIP2;1 is a gated water channel, which switches between open and closed states depending on the conformation of a 20-residue cytoplasmic loop, the D-loop. Using fully atomistic molecular dynamics simulations, we have investigated the possibility of driving the conformational equilibrium of the protein toward a constitutively open state. We introduce two separate mutations in the D-loop, while being in the closed conformation. We show that the single channel permeability of both mutants is comparable to that of the open conformation. This Article provides new molecular insight into the gating mechanism of SoPIP2;1. It is proposed that residues Arg190, Asp191, and Ser36 might play important roles in the gating of the protein.

Escherichia coli lactose permease (LacY) transports sugar across the inner membrane of the bacterium using the proton motive force to accumulate sugar in the cytosol. We have probed lactose conduction across LacY using steered molecular dynamics, permitting us to follow molecular and energetic details of lactose interaction with the lumen of LacY during its permeation. Lactose induces a widening of the narrowest parts of the channel during permeation, the widening being largest within the periplasmic half-channel. During permeation, the water-filled lumen of LacY only partially hydrates lactose, forcing it to interact with channel lining residues. Lactose forms a multitude of direct sugar-channel hydrogen bonds, predominantly with residues of the flexible N-domain, which is known to contribute a major part of LacY's affinity for lactose. In the periplasmic half-channel lactose predominantly interacts with hydrophobic channel lining residues, whereas in the cytoplasmic half-channel key protein-substrate interactions are mediated by ionic residues. A major energy barrier against transport is found within a tight segment of the periplasmic half-channel where sugar hydration is minimal and protein-sugar interaction maximal. Upon unbinding from the binding pocket, lactose undergoes a rotation to permeate either half-channel with its long axis aligned parallel to the channel axis. The results hint at the possibility of a transport mechanism, in which lactose permeates LacY through a narrow periplasmic half-channel and a wide cytoplasmic half-channel, the opening of which is controlled by changes in protonation states of key protein side groups.

Molecular dynamics simulations of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers using the CHARMM27 force field in the tensionless isothermal-isobaric (NPT) ensemble give highly ordered, gel-like bilayers with an area per lipid of approximately 48 A(2). To obtain fluid (L(alpha)) phase properties of DPPC bilayers represented by the CHARMM energy function in this ensemble, we reparameterized the atomic partial charges in the lipid headgroup and upper parts of the acyl chains. The new charges were determined from the electron structure using both the Mulliken method and the restricted electrostatic potential fitting method. We tested the derived charges in molecular dynamics simulations of a fully hydrated DPPC bilayer. Only the simulation with the new restricted electrostatic potential charges shows significant improvements compared with simulations using the original CHARMM27 force field resulting in an area per lipid of 60.4 +/- 0.1 A(2). Compared to the 48 A(2), the new value of 60.4 A(2) is in fair agreement with the experimental value of 64 A(2). In addition, the simulated order parameter profile and electron density profile are in satisfactory agreement with experimental data. Thus, the biologically more interesting fluid phase of DPPC bilayers can now be simulated in all-atom simulations in the NPT ensemble by employing our modified CHARMM27 force field.

To investigate substrate recruitment and transport across the Escherichia coli Ammonia transporter B (AmtB) protein, we performed molecular dynamics simulations of the AmtB trimer. We have identified residues important in recruitment of ammonium and intraluminal binding sites selective of ammonium, which provide a means of cation selectivity. Our results indicate that A162 guides translocation of an extraluminal ammonium into the pore lumen. We propose a mechanism for transporting the intraluminally recruited proton back to periplasm. Our mechanism conforms to net transport of ammonia and can explain why ammonia conduction is lost upon mutation of the conserved residue D160. We unify previous suggestions of D160 having either a structural or an ammonium binding function. Finally, our simulations show that the channel lumen is hydrated from the cytoplasmic side via the formation of single file water, while the F107/F215 stack at the inner-most part of the periplasmic vestibule constitutes a hydrophobic filter preventing AmtB from conducting water.

Integral membrane proteins are characterized by having a preference for aromatic residues, e.g., tryptophan (W), at the interface between the lipid bilayer core and the aqueous phase. The reason for this is not clear, but it seems that the preference is related to a complex interplay between steric and electrostatic forces. The flat rigid paddle-like structure of tryptophan, associated with a quadrupolar moment (aromaticity) arising from the pi-electron cloud of the indole, interacts primarily with moieties in the lipid headgroup region hardly penetrating into the bilayer core. We have studied the interaction between the nitrogen moiety of lipid molecule headgroups and the pi-electron distribution of gramicidin (gA) tryptophan residues (W9, W11, W13, and W15) using molecular dynamics (MD) simulations of gA embedded in two hydrated lipid bilayers composed of 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) and 1-palmitoyl-2-oleoylphosphatidyl-choline (POPC), respectively. We use a force field model for tryptophan in which polarizability is only implicit, but we believe that classical molecular dynamics force fields are sufficient to capture the most prominent features of the cation-pi interaction. Our criteria for cation-pi interactions are based on distance and angular requirements, and the results from our model suggest that cation-pi interactions are relevant for W(PE)1), W(PE)13, W(PE)15, and, to some extent, W(PC)11 and W(PC)13. In our model, W9 does not seem to engage in cation-pi interactions with lipids, neither in POPE nor POPC. The criteria for the cation-pi effect are satisfied more often in POPE than in POPC, whereas the H-bonding ability between the indole donor and the carbonyl acceptor is similar in POPE and POPC. This suggests an increased affinity for lipids with ethanolamine headgroups to transmembrane proteins enriched in interfacial tryptophans.

The thermotropic phase behavior and lateral structure of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers containing an acylated peptide has been characterized by differential scanning calorimetry (DSC) on vesicles and atomic force microscopy (AFM) on mica-supported bilayers. The acylated peptide, which is a synthetic decapeptide N-terminally linked to a C14 acyl chain (C14-peptide), is incorporated into DPPC bilayers in amounts ranging from 0-20 mol %. The calorimetric scans of the two-component system demonstrate a distinct influence of the C14-peptide on the lipid bilayer thermodynamics. This is manifested as a concentration-dependent downshift of both the main phase transition and the pretransition. In addition, the main phase transition peak is significantly broadened, indicating phase coexistence. In the AFM imaging scans we found that the C14-peptide, when added to supported gel phase DPPC bilayers, inserts preferentially into preexisting defect regions and has a noticeable influence on the organization of the surrounding lipids. The presence of the C14-peptide gives rise to a laterally heterogeneous bilayer structure with coexisting lipid domains characterized by a 10 A height difference. The AFM images also show that the appearance of the ripple phase of the DPPC lipid bilayers is unaffected by the C14-peptide. The experimental results are supported by molecular dynamics simulations, which show that the C14-peptide has a disordering effect on the lipid acyl chains and causes a lateral expansion of the lipid bilayer. These effects are most pronounced for gel-like bilayer structures and support the observed downshift in the phase-transition temperature. Moreover, the molecular dynamics data indicate a tendency of a tryptophan residue in the peptide sequence to position itself in the bilayer headgroup region.

Water permeation and electrostatic interactions between water and channel are investigated in the Escherichia coli glycerol uptake facilitator GlpF, a member of the aquaporin water channel family, by molecular dynamics simulations. A tetrameric model of the channel embedded in a 16:0/18:1c9-palmitoyloleylphosphatidylethanolamine membrane was used for the simulations. During the simulations, water molecules pass through the channel in single file. The movement of the single file water molecules through the channel is concerted, and we show that it can be described by a continuous-time random-walk model. The integrity of the single file remains intact during the permeation, indicating that a disrupted water chain is unlikely to be the mechanism of proton exclusion in aquaporins. Specific hydrogen bonds between permeating water and protein at the channel center (at two conserved Asp-Pro-Ala "NPA" motifs), together with the protein electrostatic fields enforce a bipolar water configuration inside the channel with dipole inversion at the NPA motifs. At the NPA motifs water-protein electrostatic interactions facilitate this inversion. Furthermore, water-water electrostatic interactions are in all regions inside the channel stronger than water-protein interactions, except near a conserved, positively charged Arg residue. We find that variations of the protein electrostatic field through the channel, owing to preserved structural features, completely explain the bipolar orientation of water. This orientation persists despite water translocation in single file and blocks proton transport. Furthermore, we find that for permeation of a cation, ion-protein electrostatic interactions are more unfavorable at the conserved NPA motifs than at the conserved Arg, suggesting that the major barrier against proton transport in aquaporins is faced at the NPA motifs.