Concentration Dependent Ion-Protein Interaction Patterns Underlying Protein Oligomerization Behaviours.

Salts and proteins comprise two of the basic molecular components of biological materials. Kosmotropic/chaotropic co-solvation and matching ion water affinities explain basic ionic effects on protein aggregation observed in simple solutions. However, it is unclear how these theories apply to proteins in complex biological environments and what the underlying ionic binding patterns are. Using the positive ion Ca(2+) and the negatively charged membrane protein SNAP25, we studied ion effects on protein oligomerization in solution, in native membranes and in molecular dynamics (MD) simulations. We find that concentration-dependent ion-induced protein oligomerization is a fundamental chemico-physical principle applying not only to soluble but also to membrane-anchored proteins in their native environment. Oligomerization is driven by the interaction of Ca(2+) ions with the carboxylate groups of aspartate and glutamate. From low up to middle concentrations, salt bridges between Ca(2+) ions and two or more protein residues lead to increasingly larger oligomers, while at high concentrations oligomers disperse due to overcharging effects. The insights provide a conceptual framework at the interface of physics, chemistry and biology to explain binding of ions to charged protein surfaces on an atomistic scale, as occurring during protein solubilisation, aggregation and oligomerization both in simple solutions and membrane systems.

Computer simulations suggest direct and stable tip to tip interaction between the outer membrane channel TolC and the isolated docking domain of the multidrug RND efflux transporter AcrB.

One way by which bacteria achieve antibiotics resistance is preventing drug access to its target molecule for example through an overproduction of multi-drug efflux pumps of the resistance nodulation division (RND) protein super family of which AcrAB-TolC in Escherichia coli is a prominent example. Although representing one of the best studied efflux systems, the question of how AcrB and TolC interact is still unclear as the available experimental data suggest that either both proteins interact in a tip to tip manner or do not interact at all but are instead connected by a hexamer of AcrA molecules. Addressing the question of TolC-AcrB interaction, we performed a series of 100 ns - 1 µs-molecular dynamics simulations of membrane-embedded TolC in presence of the isolated AcrB docking domain (AcrB(DD)). In 5/6 simulations we observe direct TolC-AcrB(DD) interaction that is only stable on the simulated time scale when both proteins engage in a tip to tip manner. At the same time we find TolC opening and closing freely on extracellular side while remaining closed at the inner periplasmic bottleneck region, suggesting that either the simulated time is too short or additional components are required to unlock TolC.

Application of room-temperature aprotic and protic ionic liquids for oxidative folding of cysteine-rich peptides.

The oxidation of the conotoxin µ-SIIIA in different ionic liquids was investigated, and the results were compared with those obtained in [C2 mim][OAc]. Conversion of the reduced precursor into the oxidized product was observed in the protic ILs methyl- and ethylammonium formate (MAF and EAf, respectively), whereas choline dihydrogenphosphate and Ammoeng 110 failed to yield folded peptide. However, the quality and yield of the peptide obtained in MAF and EAF were lower than in the case of the product from [C2 mim][OAc]. Reaction conditions (temperature, water content) also had an impact on peptide conversion. A closer look at the activities of µ-SIIIA versions derived from an up-scaled synthesis in [C2 mim][OAc] revealed a significant loss of the effect on ion channel NaV 1.4 relative to the buffer-oxidized peptide, whereas digestion of either µ-SIIIA product by trypsin was unaffected. This was attributed to adherence of ions from the IL to the peptide, because the disulfide connectivity is basically the same for the differentially oxidized µ-SIIIA versions.

Structural dynamics of the cell wall precursor lipid II in the presence and absence of the lantibiotic nisin.

Representing a physiological "Achilles' heel", the cell wall precursor lipid II (LII) is a prime target for various classes of antibiotics. Over the years LII-binding agents have been recognized as promising candidates and templates in the search for new antibacterial compounds to complement or replace existing drugs. To elucidate the molecular structural basis underlying LII functional mechanism and to better understand if and how lantibiotic binding alters the molecular behavior of LII, we performed molecular dynamics (MD) simulations of phospholipid membrane-embedded LII in the absence and presence of the LII-binding lantibiotic nisin. In a series of 2×4 independent, unbiased 100ns MD simulations we sampled the conformational dynamics of nine LII as well as nine LII-nisin complexes embedded in an aqueous 150mM NaCl/POPC phospholipid membrane environment. We found that nisin binding to LII induces a reduction of LII mobility and flexibility, an outward shift of the LII pentapeptide, an inward movement of the LII disaccharide section, and an overall deeper insertion of the LII tail group into the membrane. The latter effect might indicate an initial step in adopting a stabilizing, scaffold-like structure in the process of nisin-induced membrane leakage. At the same time nisin conformation and LII interaction remain similar to the 1WCO LII-nisin NMR solution structure.

Efflux pump-mediated antibiotics resistance: insights from computational structural biology.

The continuous rise of bacterial resistance against formerly effective pharmaceuticals is a major challenge for biomedical research. Since the first computational studies published seven years ago the simulation-based investigation of antibiotics resistance mediated by multidrug efflux pumps of the resistance nodulation division (RND) protein super family has grown into a vivid field of research. Here we review the employment of molecular dynamics computer simulations to investigate RND efflux pumps focusing on our group's recent contributions to this field studying questions of energy conversion and substrate transport in the inner membrane antiporter AcrB in Escherichia coli as well as access regulation and gating mechanism in the outer membrane efflux ducts TolC and OprM in E. coli and Pseudomonas aeruginosa.

Porter domain opening and closing motions in the multi-drug efflux transporter AcrB.

Acriflavine resistance protein B acts as the active transporter in the multi-drug efflux pump Acriflavine resistance proteins A / B - Tolerance to colicins protein in Escherichia coli. Within the same reaction cycle intermediate all Acriflavine resistance protein B X-ray structures display highly similar conformations of the substrate-recruiting and transporting porter domain. To assess if this structural homogeneity is an intrinsic feature of Acriflavine resistance protein B or stems from other causes we performed a series of six independent, unbiased 100 ns molecular dynamics simulations of membrane-embedded, asymmetric, substrate-free wild type Acriflavine resistance protein B in a 150 mM NaCl solution. We find the porter domain more flexible than previously assumed displaying clear opening and closing motions of the proximal binding pocket (L and T-state) and the exit of the drug transport channels (O-intermediate). Concurrently the hydrophobic binding pocket favors a closed conformation in all three protomers. Our findings suggest that the conformational homogeneity seen in the crystal structures is likely an effect of bound but structurally unresolved substrate. Our simulations further imply that each of the known three reaction cycle intermediates occurs in at least two variants, the Thr676 loop independently regulates porter domain access and likely plays a key role in substrate transport. On a 100 ns time scale we find no evidence supporting the proposed LLL resting state in the absence of substrate. If the proximal binding pocket dynamics have an inhibiting effect on Acriflavine resistance protein B pump activity lowering the life time of substrate-accessible conformations, the observed dynamics could provide a structural explanation for the Acriflavine resistance protein B activity-enhancing effect of the adaptor protein Acriflavine resistance protein A stabilizing PC1 and PC2 subdomain orientations.

Unilateral access regulation: ground state dynamics of the Pseudomonas aeruginosa outer membrane efflux duct OprM.

Acting as an efflux duct in the MexA-MexB-OprM multidrug efflux pump, OprM plays a major role in the antibiotic resistance capability of Pseudomonas aeruginosa, trafficking substrates through the outer cell membrane. Whereas the available crystal structures showed restricted OprM access on both ends, the underlying gating mechanism is not yet fully understood. To gain insight into the functional mechanism of OprM access regulation, we conducted a series of five independent, unbiased molecular dynamics simulations, computing 200 ns dynamics samples of the wild-type protein in a phospholipid membrane/150 mM NaCl water environment. On the extracellular side, OprM opens and closes freely under the simulated conditions, suggesting the absence of a gating mechanism on this side of the isolated protein. On the periplasmic side, we observe an opening of the tip regions at Val408 and to a lesser degree Asp416 located 1.5 nm further into the channel, leading to OprM end conformations being up to 3 and 1.4 times, respectively, more open than the asymmetric crystal structure. If our simulations are correct, our findings imply that periplasmic gating involves only the Asp416 region and that in vivo additional components, absent in our simulation, might be required for periplasmic gating if the observed opening trend near Asp416 is not negligible. In addition to that ,we identified in each monomer a previously unreported sodium binding site in the channel interior coordinated by Asp171 and Asp230 whose functional role remains to be investigated.

Molecular Dynamics Computer Simulations of Multidrug RND Efflux Pumps.

Over-expression of multidrug efflux pumps of the Resistance Nodulation Division (RND) protein super family counts among the main causes for microbial resistance against pharmaceuticals. Understanding the molecular basis of this process is one of the major challenges of modern biomedical research, involving a broad range of experimental and computational techniques. Here we review the current state of RND transporter investigation employing molecular dynamics simulations providing conformational samples of transporter components to obtain insights into the functional mechanism underlying efflux pump-mediated antibiotics resistance in Escherichia coli and Pseudomonas aeruginosa.

LAMBADA and InflateGRO2: efficient membrane alignment and insertion of membrane proteins for molecular dynamics simulations.

At the beginning of each molecular dynamics membrane simulation stands the generation of a suitable starting structure which includes the working steps of aligning membrane and protein and seamlessly accommodating the protein in the membrane. Here we introduce two efficient and complementary methods based on pre-equilibrated membrane patches, automating these steps. Using a voxel-based cast of the coarse-grained protein, LAMBADA computes a hydrophilicity profile-derived scoring function based on which the optimal rotation and translation operations are determined to align protein and membrane. Employing an entirely geometrical approach, LAMBADA is independent from any precalculated data and aligns even large membrane proteins within minutes on a regular workstation. LAMBADA is the first tool performing the entire alignment process automatically while providing the user with the explicit 3D coordinates of the aligned protein and membrane. The second tool is an extension of the InflateGRO method addressing the shortcomings of its predecessor in a fully automated workflow. Determining the exact number of overlapping lipids based on the area occupied by the protein and restricting expansion, compression and energy minimization steps to a subset of relevant lipids through automatically calculated and system-optimized operation parameters, InflateGRO2 yields optimal lipid packing and reduces lipid vacuum exposure to a minimum preserving as much of the equilibrated membrane structure as possible. Applicable to atomistic and coarse grain structures in MARTINI format, InflateGRO2 offers high accuracy, fast performance, and increased application flexibility permitting the easy preparation of systems exhibiting heterogeneous lipid composition as well as embedding proteins into multiple membranes. Both tools can be used separately, in combination with other methods, or in tandem permitting a fully automated workflow while retaining a maximum level of usage control and flexibility. To assess the performance of both methods, we carried out test runs using 22 membrane proteins of different size and transmembrane structure.

Locked on one side only: ground state dynamics of the outer membrane efflux duct TolC.

Playing a major role in the expulsion of antibiotics and the secretion of cell toxins in conjunction with inner membrane transporters of three protein superfamilies, the outer membrane channel TolC occurs in at least two states blocking or permitting the passage of substrates. The details of the underlying gating mechanism are not fully understood. Addressing the questions of extracellular access control and periplasmic gating mechanism, we conducted a series of independent, unbiased 150-300 ns molecular dynamics simulations of wild-type TolC in a phospholipid membrane/150 mM NaCl water environment. We find that TolC opens and closes freely on the extracellular side, suggesting the absence of a gating mechanism on this side in the isolated protein. On the periplasmic side, we observe the outer periplasmic bottleneck region adopting in all simulations a conformation more open than the TolC wild-type crystal structures until in one run the successive binding of two sodium ions induces the transition to a conformation more closed than any of the available TolC X-ray structures. Concurrent with a heightened sodium residence probability near Asp374, the inner periplasmic bottleneck region at Asp374 remains closed throughout the simulations unless all NaCl is removed from the system, inducing a reopening of the outer and inner bottleneck. Our findings suggest that TolC is locked only on the periplasmic side in a sodium-dependent manner.

Three ways in, one way out: water dynamics in the trans-membrane domains of the inner membrane translocase AcrB.

Powered by proton-motive force, the inner membrane translocase AcrB is the engine of the AcrAB-TolC efflux pump in Escherichia coli. As proton conduction in proteins occurs along hydrogen-bonded networks of polar residues and water molecules, knowledge of the protein-internal water distribution and water-interacting residues allows drawing conclusions to possible pathways of proton conduction. Here, we report a series of 6× 50 ns independent molecular dynamics simulations of asymmetric AcrB embedded in a phospholipid/water environment. Simulating each monomer in its proposed protonation state, we calculated for each trans-membrane domain the average water distribution, identified residues interacting with these waters and quantified each residue's frequency of water hydrogen bond contact. Combining this information we find three possible routes of proton transfer connecting a continuously hydrated region of known key residues in the TMD interior to bulk water by one cytoplasmic and up to three periplasm water channels in monomer B and A. We find that water access of the trans-membrane domains is regulated by four groups of residues in a combination of side chain re-orientations and shifts of trans-membrane helices. Our findings support a proton release event via Arg971 during the C intermediate or in the transition to A, and proton uptake occurring in the A or B state or during a so far unknown intermediate in between B and C where cytoplasmic water access is still possible. Our simulations suggest experimentally testable hypotheses, which have not been investigated so far.

dxTuber: detecting protein cavities, tunnels and clefts based on protein and solvent dynamics.

Empty space in a protein structure can provide valuable insight into protein properties such as internal hydration, structure stabilization, substrate translocation, storage compartments or binding sites. This information can be visualized by means of cavity analysis. Numerous tools are available depicting cavities directly or identifying lining residues. So far, all available techniques base on a single conformation neglecting any form of protein and cavity dynamics. Here we report a novel, grid-based cavity detection method that uses protein and solvent residence probabilities derived from molecular dynamics simulations to identify (I) internal cavities, (II) tunnels or (III) clefts on the protein surface. Driven by a graphical user interface, output can be exported in PDB format where cavities are described as individually selectable groups of adjacent voxels representing regions of high solvent residence probability. Cavities can be analyzed in terms of solvent density, cavity volume and cross-sectional area along a principal axis. To assess dxTuber performance we performed test runs on a set of six example proteins representing the three main classes of protein cavities and compared our findings to results obtained with SURFNET, CAVER and PyMol.

Membrane protein dynamics from femtoseconds to seconds.

Membrane proteins play a key role in energy conversion, transport, signal recognition, transduction, and other fundamental biological processes. Despite considerable progress in experimental techniques, the determination of structure and dynamics of membrane proteins still represents a great challenge. Computer simulation methods are becoming an increasingly important tool not only in the interpretation of experiments but also in the prediction of membrane protein dynamics. In the present review, we give a brief introduction to molecular modeling techniques currently used to explore protein dynamics on time scales ranging from femtoseconds to microseconds. We then describe a few recent example applications of these techniques to membrane proteins. In conclusion, we also discuss some of the newest developments in simulation methodology that have the potential to further extend the time scale accessible to explore (membrane) protein dynamics.

Holo-BtuF stabilizes the open conformation of the vitamin B12 ABC transporter BtuCD.

While there is evidence that other ABC transporters can tell between empty and loaded substrate binding protein, reconstitution experiments suggest otherwise for the Escherichia coli vitamin B12 importer BtuCD-F. Here, we address the question of BtuCD-F substrate sensitivity in a combined protein-protein docking and molecular dynamics simulation approach. Starting from the BtuCD and holo-BtuF crystal structures, we model two holo-BtuCD-F docking complexes differing by a 180 degrees orientation of BtuF. One of these is similar to the apo-BtuCD-F crystal structure. Both docking complexes were embedded in a lipid/water environment to investigate their dynamics and BtuCD's conformational response to the presence and absence of BtuF, vitamin B12, and Mg-ATP in a series of 28 independent MD simulations. We find holo-BtuF stabilizing the open conformation of BtuCD, whereas the transporter begins to close again when BtuF or vitamin B12 is removed-suggesting BtuCD-F is capable of substrate sensitivity. We identified BtuC transmembrane helices 3 and 5, the L-loops and the adjacent helices comprised of BtuC residues 170-180 as hotspots of conformational change. We propose the latter to act as substrate sensors. BtuF-Trp44 appears to act as a lid on the vitamin B12 binding cleft in BtuF X-ray structures and protrudes into the BtuCD transport channel in one of our simulations, which might represent an initial step in vitamin B12 uptake. On an average, we observe subunit motions where the nucleotide binding domains approach each other while the transmembrane domains display an opening trend toward the periplasm.

ATP-binding cassette transporters in Escherichia coli.

ATP-binding cassette (ABC) transporters are integral membrane proteins that actively transport molecules across cell membranes. In Escherichia coli they consist primarily of import systems that involve in addition to the ABC transporter itself a substrate binding protein and outer membrane receptors or porins, and a number of transporters with varied functions. Recent crystal structures of a number of ATPase domains, substrate binding proteins, and full-length transporters have given new insight in the molecular basis of transport. Bioinformatics approaches allow an approximate identification of all ABC transporters in E. coli and their relation to other known transporters. Computational approaches involving modeling and simulation are beginning to yield insight into the dynamics of the transporters. We summarize the function of the known ABC transporters in E. coli and mechanistic insights from structural and computational studies.

Computer simulation of antimicrobial peptides.

Naturally occurring and synthetic peptides may be a novel class of clinically useful antibiotics. A large body of experimental data on structure function relationships for such peptides is available, but the molecular mechanism of their action remains elusive in most cases. Computer simulations can give detailed insights into the interactions between peptides and lipid bilayers, at least one crucial step in the antimicrobial mechanism. Here we review recent simulations of antimicrobial peptides and discuss potential future contributions of computer simulations in understanding and ultimately designing antimicrobial peptides.

The distribution and conformation of very long-chain plant wax components in a lipid bilayer.

Plant wax contains long-chain alkanes and related components which are transported to the surface of the plant by specialized ABC transporter proteins. Here, we determine the distribution and conformation of three wax components, nonacosane, nonacosan-15-one, and nonacosan-15-ol, using unbiased and umbrella sampling molecular dynamics simulations. The molecules all partitioned to the center of the bilayer, with a free-energy difference of -70 kJ/mol between bulk water and the center of the bilayer for the alkane and -55 kJ/mol for the two more-polar molecules. All of the wax molecules were highly mobile in the bilayer, freely moving between opposite leaflets on a time scale of a few nanoseconds. Nonacosan-15-one and nonacosan-15-ol folded double to expose their hydrophilic group to the solvent, whereas nonacosane alternated between orientations spanning the full bilayer and orientations in the center of the bilayer.

Setting up and running molecular dynamics simulations of membrane proteins.

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.

Simulation of the coupling between nucleotide binding and transmembrane domains in the ATP binding cassette transporter BtuCD.

The nucleotide-induced structural rearrangements in ATP binding cassette (ABC) transporters, leading to substrate translocation, are largely unknown. We have modeled nucleotide binding and release in the vitamin B(12) importer BtuCD using perturbed elastic network calculations and biased molecular dynamics simulations. Both models predict that nucleotide release decreases the tilt between the two transmembrane domains and opens the cytoplasmic gate. Nucleotide binding has the opposite effect. The observed coupling may be relevant for all ABC transporters because of the conservation of nucleotide binding domains and the shared role of ATP in ABC transporters. The rearrangements in the cytoplasmic gate region do not provide enough space for B(12) to diffuse from the transporter pore into the cytoplasm, which could suggest that peristaltic forces are needed to exclude B(12) from the transporter pore.

Opening and closing motions in the periplasmic vitamin B12 binding protein BtuF.

BtuF is the periplasmic binding protein (PBP) in the vitamin B(12) uptake system in Escherichia coli where it is associated with the ABC transporter BtuCD. When the ligand binds, PBPs generally display large conformational changes, commonly termed the Venus flytrap mechanism. BtuF belongs to a group of PBPs that, on the basis of crystal structures, does not appear to display such behavior. Using 480 ns multicopy molecular dynamics simulations of apo and holo forms of the protein, we investigate the dynamics of BtuF. We find BtuF to be more flexible than previously assumed, displaying clear opening and closing motions which are more pronounced in the apo form. The protein behavior is compatible with a PBP functional model that postulates a closed conformation for the ligand-bound state, whereas the empty form fluctuates between open and closed conformations. Elastic network normal-mode analysis suggests that all BtuF-like PBPs are capable of similar opening and closing motions. It also makes the typical Venus flytrap domain motions a likely common means of how PBP-ABC transporter interaction could occur.