Investigating the Bioactive Conformation of Angiotensin II Using Markov State Modeling Revisited with Web-Scale Clustering.

Molecular dynamics simulation is a powerful technique for studying the structure and dynamics of biomolecules in atomic-level detail by sampling their various conformations in real time. Because of the long timescales that need to be sampled to study biomolecular processes and the big and complex nature of the corresponding data, relevant analyses of important biophysical phenomena are challenging. Clustering and Markov state models (MSMs) are efficient computational techniques that can be used to extract dominant conformational states and to connect those with kinetic information. In this work, we perform Molecular Dynamics simulations to investigate the free energy landscape of Angiotensin II (AngII) in order to unravel its bioactive conformations using different clustering techniques and Markov state modeling. AngII is an octapeptide hormone, which binds to the AT1 transmembrane receptor, and plays a vital role in the regulation of blood pressure, conservation of total blood volume, and salt homeostasis. To mimic the water-membrane interface as AngII approaches the AT1 receptor and to compare our findings with available experimental results, the simulations were performed in water as well as in water-ethanol mixtures. Our results show that in the water-ethanol environment, AngII adopts more compact U-shaped (folded) conformations than in water, which resembles its structure when bound to the AT1 receptor. For clustering of the conformations, we validate the efficiency of an inverted-quantized k-means algorithm, as a fast approximate clustering technique for web-scale data (millions of points into thousands or millions of clusters) compared to k-means, on data from trajectories of molecular dynamics simulations with reasonable trade-offs between time and accuracy. Finally, we extract MSMs using various clustering techniques for the generation of microstates and macrostates and for the selection of the macrostate representatives.

Insights into the mechanism of the PIK3CA E545K activating mutation using MD simulations.

Phosphoinositide 3-kinase alpha (PI3Kα) is involved in fundamental cellular processes including cell proliferation and differentiation and is frequently mutated in human malignancies. One of the most common mutations is E545K, which results in an amino acid substitution of opposite charge. It has been recently proposed that in this oncogenic charge-reversal mutation, the interactions between the protein catalytic and regulatory subunits are abrogated, resulting in loss of regulation and constitutive PI3Kα activity, which can lead to oncogenesis. To assess the mechanism of the PI3Kα E545K activating mutation, extensive Molecular Dynamics simulations were performed to examine conformational changes differing between the wild type (WT) and mutant proteins as they occur in microsecond simulations. In the E545K mutant PI3Kα, we observe a spontaneous detachment of the nSH2 PI3Kα domain (regulatory subunit, p85α) from the helical domain (catalytic subunit, p110α) causing significant loss of communication between the regulatory and catalytic subunits. We examine the allosteric network of the two proteins and show that a cluster of residues around the mutation is important for delivering communication signals between the catalytic and regulatory subunits. Our results demonstrate the dynamical and structural effects induced by the p110α E545K mutation in atomic level detail and indicate a possible mechanism for the loss of regulation that E545K confers on PI3Kα.

Toroidal pores formed by antimicrobial peptides show significant disorder.

A large variety of antimicrobial peptides have been shown to act, at least in vitro, by poration of the lipid membrane. The nanometre size of these pores, however, complicates their structural characterization by experimental techniques. Here we use molecular dynamics simulations, to study the interaction of a specific class of antimicrobial peptides, melittin, with a dipalmitoylphosphatidylcholine bilayer in atomic detail. We show that transmembrane pores spontaneously form above a critical peptide to lipid ratio. The lipid molecules bend inwards to form a toroidally shaped pore but with only one or two peptides lining the pore. This is in strong contrast to the traditional models of toroidal pores in which the peptides are assumed to adopt a transmembrane orientation. We find that peptide aggregation, either prior or after binding to the membrane surface, is a prerequisite to pore formation. The presence of a stable helical secondary structure of the peptide, however is not. Furthermore, results obtained with modified peptides point to the importance of electrostatic interactions in the poration process. Removing the charges of the basic amino-acid residues of melittin prevents pore formation. It was also found that in the absence of counter ions pores not only form more rapidly but lead to membrane rupture. The rupture process occurs via a novel recursive poration pathway, which we coin the Droste mechanism.

Ion transport across transmembrane pores.

To study the pore-mediated transport of ionic species across a lipid membrane, a series of molecular dynamics simulations have been performed of a dipalmitoyl-phosphatidyl-choline bilayer containing a preformed water pore in the presence of sodium and chloride ions. It is found that the stability of the transient water pores is greatly reduced in the presence of the ions. Specifically, the binding of sodium cations at the lipid/water interface increases the pore line tension, resulting in a destabilization of the pore. However, the application of mechanical stress opposes this effect. The flux of ions through these mechanically stabilized pores has been analyzed. Simulations indicate that the transport of the ions through the pores depends strongly on the size of the water channel. In the presence of small pores (radius <1.5 nm) permeation is slow, with both sodium and chloride permeating at similar rates. In the case in which the pores are larger (radius >1.5 nm), a crossover is observed to a regime where the anion flux is greatly enhanced. Based on these observations, a mechanism for the basal membrane permeability of ions is discussed.

Antimicrobial peptides in action.

Molecular dynamics simulations of the magainin MG-H2 peptide interacting with a model phospholipid membrane have been used to investigate the mechanism by which antimicrobial peptides act. Multiple copies of the peptide were randomly placed in solution close to the membrane. The peptide readily bound to the membrane, and above a certain concentration, the peptide was observed to cooperatively induce the formation of a nanometer-sized, toroidally shaped pore in the bilayer. In sharp contrast with the commonly accepted model of a toroidal pore, only one peptide was typically found near the center of the pore. The remaining peptides lay close to the edge of the pore, maintaining a predominantly parallel orientation with respect to the membrane.

Molecular dynamics simulations of hydrophilic pores in lipid bilayers.

Hydrophilic pores are formed in peptide free lipid bilayers under mechanical stress. It has been proposed that the transport of ionic species across such membranes is largely determined by the existence of such meta-stable hydrophilic pores. To study the properties of these structures and understand the mechanism by which pore expansion leads to membrane rupture, a series of molecular dynamics simulations of a dipalmitoylphosphatidylcholine (DPPC) bilayer have been conducted. The system was simulated in two different states; first, as a bilayer containing a meta-stable pore and second, as an equilibrated bilayer without a pore. Surface tension in both cases was applied to study the formation and stability of hydrophilic pores inside the bilayers. It is observed that below a critical threshold tension of approximately 38 mN/m the pores are stabilized. The minimum radius at which a pore can be stabilized is 0.7 nm. Based on the critical threshold tension the line tension of the bilayer was estimated to be approximately 3 x 10(-11) N, in good agreement with experimental measurements. The flux of water molecules through these stabilized pores was analyzed, and the structure and size of the pores characterized. When the lateral pressure exceeds the threshold tension, the pores become unstable and start to expand causing the rupture of the membrane. In the simulations the mechanical threshold tension necessary to cause rupture of the membrane on a nanosecond timescale is much higher in the case of the equilibrated bilayers, as compared with membranes containing preexisting pores.

Simulation of pore formation in lipid bilayers by mechanical stress and electric fields.

We describe computer simulations of pore formation and membrane rupture of phospholipid bilayers under mechanical and electrical stress. On the nanosecond simulation time scale, pores are induced by a lateral pressure exceeding -200 bar or by an applied electric field of 0.5 V/nm.