Exploring models of the influenza A M2 channel: MD simulations in a phospholipid bilayer

The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-selective channels. The native form of the protein is 97 residues long, although peptides representing the transmembrane section display ion channel activity, which (like the native channel) is blocked by the antiviral drug amantadine. As a small ion channel, M2 may provide useful insights into more complex channel systems. Models of tetrameric bundles of helices containing either 18 or 22 residues have been simulated while embedded in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylcholine bilayer. Several different starting models have been used. These suggest that the simulation results, at least on a nanosecond time scale, are sensitive to the exact starting structure. Electrostatics calculations carried out on a ring of four ionizable aspartate residues at the N-terminal mouth of the channel suggest that at any one time, only one will be in a charged state. Helix bundle models were mostly stable over the duration of the simulation, and their helices remained tilted relative to the bilayer normal. The M2 helix bundles form closed channels that undergo breathing motions, alternating between a tetramer and a dimer-of-dimers structure. Under these conditions either the channel forms a pocket of trapped waters or it contains a column of waters broken predominantly at the C-terminal mouth of the pore. These waters exhibit restricted motion in the pore and are effectively "frozen" in a way similar to those seen in previous simulations of a proton channel formed by a four-helix bundle of a synthetic leucine-serine peptide (, Biophys. J. 77:2400-2410).     

The TatA subunit of Escherichia coli twin-arginine translocase has an N-in topology

The twin-arginine translocase (Tat) system is used by many bacteria to translocate folded proteins across the cytoplasmic membrane. The TatA subunit is the predicted pore-forming subunit and has been shown to form a homo-oligomeric complex. Through accessibility experiments using the thiol-reactive reagents 4-acetamido-4'-maleimidylstilbene-2,2'-disulfonic acid and Nalpha-(3-maleimidylproprionyl)biocytin toward site-specific cysteine mutants in TatA, we show that the N-terminus of TatA is located in the cytoplasm rather than the previously assumed periplasm. We also confirm previous observations that the C-terminus has a dual topology. By treatment with the membrane uncoupler carbonyl cyanide-m-chlorophenyl hydrazone, we show that the topological state of the C-terminus is dependent on the membrane potential. These results suggest two architectures of TatA in the membrane: one with a single transmembrane helix and the other with two transmembrane helices. Molecular models of both topologies were used to develop and cartoon a homo-oligomeric complex as a channel with a diameter of approximately 50 A and suggest that the double transmembrane helix topology might be the building block for the translocation channel. Additionally, in vivo cross-linking experiments of Gly2Cys and Thr22Cys mutants showed that Gly2, at the beginning of transmembrane helix-1, is in close proximity with Gly2 of a neighboring TatA, as Cys2 cross-linked immediately upon the addition of copper phenanthroline. On the other hand, Cys22, at the other end of the transmembrane helix, took at least 10 min to cross-link, suggesting that a possible movement or reorientation is required to bring this residue into proximity with a neighboring TatA subunit.     

A molecular dynamics study of the pores formed by Escherichia coli OmpF porin in a fully hydrated palmitoyloleoylphosphatidylcholine bilayer.

In this paper we study the properties of pores formed by OmpF porin from Escherichia coli, based on a molecular dynamics simulation of the OmpF trimer, 318 palmitoyl-oleoyl-phosphatidylethanolamine lipids, 27 Na+ ions, and 12,992 water molecules. After equilibration and a nanosecond production run, the OmpF trimer exhibits a C-alpha root mean square deviation from the crystal structure of 0.23 nm and a stable secondary structure. No evidence is found for large-scale motions of the L3 loop. We investigate the pore dimensions, conductance, and the properties of water inside the pore. This water forms a complicated pattern, even when averaged over 1 ns of simulation time. Around the pore constriction zone the water dipoles are highly structured in the plane of the membrane, oriented by the strong transversal electric field. In addition, there is a net orientation along the pore axis pointing from the extracellular to the intracellular side of the bilayer. The diffusion coefficients of water inside the pore are greatly reduced compared to bulk. We compare our results to results from model pores (Breed et al., 1996. Biophys. J. 70:1 643-1 661; Sansom et al. 1997. Biophys. J. 73:2404-241 5) and discuss implications for further theoretical work.     

Cystic Fibrosis Transmembrane Conductance Regulator (CFTR): CLOSED AND OPEN STATE CHANNEL MODELS*

The cystic fibrosis transmembrane conductance regulator (CFTR) is a member of the ATP-binding cassette (ABC) transporter superfamily. CFTR controls the flow of anions through the apical membrane of epithelia. Dysfunctional CFTR causes the common lethal genetic disease cystic fibrosis. Transitions between open and closed states of CFTR are regulated by ATP binding and hydrolysis on the cytosolic nucleotide binding domains, which are coupled with the transmembrane (TM) domains forming the pathway for anion permeation. Lack of structural data hampers a global understanding of CFTR and thus the development of “rational” approaches directly targeting defective CFTR. In this work, we explored possible conformational states of the CFTR gating cycle by means of homology modeling. As templates, we used structures of homologous ABC transporters, namely TM(287–288), ABC-B10, McjD, and Sav1866. In the light of published experimental results, structural analysis of the transmembrane cavity suggests that the TM(287–288)-based CFTR model could correspond to a commonly occupied closed state, whereas the McjD-based model could represent an open state. The models capture the important role played by Phe-337 as a filter/gating residue and provide structural information on the conformational transition from closed to open channel.     

Estimation of divergence time between two sibling species of the <Emphasis Type="Italic">Anopheles (Kerteszia) cruzii</Emphasis>complex using a multilocus approach

Background  

Anopheles cruzii is the primary human Plasmodium vector in southern and southeastern Brazil. The distribution of this mosquito follows the coast of the Brazilian Atlantic Forest. Previous studies indicated that An. cruzii is a complex of cryptic species.     

Cytosolic region of TM6 in P-glycoprotein: topographical analysis and functional perturbation by site directed labeling

Reduced intracellular drug accumulation due to the activity of the drug efflux pump ABC (B1) is a major mechanism in the resistance of cancer cells to chemotherapy. ABC (B1) is a poly specific transporter, and the molecular mechanism of its complex translocation process remains to be elucidated. To understand the process will require information on the regions involved in drug binding and those that couple this event to nucleotide hydrolysis. The present investigation focuses on the cytosolic region of transmembrane helix 6 (TM6), which has been widely attributed with a central role in the translocation process. A series of ABC (B1) isoforms containing a unique cysteine within TM6 was constructed and the resultant proteins purified and reconstituted. Accessibility of the cysteines to covalent modification by maleimide reagents was measured for the basal, ATP bound and vanadate trapped conformations of each isoform. Residues at the two extremes of the TM6 region examined (amino acids 344 to 360) were considerably more accessible than the central segment, the latter of which also failed to undergo significant conformational changes during the catalytic cycle. Covalent modification of the cytosolic segment of TM6 did, however, attenuate drug stimulation of ATP hydrolysis and demonstrates an important role for this segment in coupling drug binding to ATP hydrolysis during translocation.     

Structure and dynamics of the antifungal molecules Syringotoxin-B and Syringopeptin-25A from molecular dynamics simulation

The phytopathogen Pseudomonas syringae pv. syringae produces toxic cyclic lipodepsipeptides (CLPs): nona-peptides and syringopeptins. All CLPs inhibit the growth of many fungal species, including human pathogens, although different fungi display different degrees of sensitivity. The best studied CLPs are Syringomycin-E (SR-E), Syringotoxin-B (ST-B) and Syringopeptin-25A (SP-25A). Their biological activity is affected by membrane composition and their structural differences. We previously (Matyus et al. in Eur Biophys J 35:459–467, 2006) reported the molecular features and structural preferences of SR-E in water and octane environments. Here we investigate in atomic detail the molecular features of the two other main CLP components, ST-B and SP-25A, in water and octane by 200 ns molecular dynamics simulations (MD), using distance restraints derived from NMR NOE data (Ballio et al. in Eur J Biochem 234:747–758, 1995). We have obtained three-dimensional models of ST-B and SP-25A CLPs in different environments. These models can now be used as a basis to investigate the interactions of ST-B and SP-25A with lipid membranes an important further step towards a better understanding of the antifungal and antibacterial activity of these peptides. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Distribution of amino acids in a lipid bilayer from computer simulations

We have calculated the distribution in a lipid bilayer of small molecules mimicking 17 natural amino acids in atomistic detail by molecular dynamics simulation. We considered both charged and uncharged forms for Lys, Arg, Glu, and Asp. The results give detailed insight in the molecular basis of the preferred location and orientation of each side chain as well the preferred charge state for ionizable residues. Partitioning of charged and polar side chains is accompanied by water defects connecting the side chains to bulk water. These water defects dominate the energetic of partitioning, rather than simple partitioning between water and a hydrophobic phase. Lys, Glu, and Asp become uncharged well before reaching the center of the membrane, but Arg may be either charged or uncharged at the center of the membrane. Phe has a broad distribution in the membrane but Trp and Tyr localize strongly to the interfacial region. The distributions are useful for the development of coarse-grained and implicit membrane potentials for simulation and structure prediction. We discuss the relationship between the distribution in membranes, bulk partitioning to cyclohexane, and several amino acid hydrophobicity scales.