Exploring the conformational dynamics of alanine dipeptide in solution subjected to an external electric field: A nonequilibrium molecular dynamics simulation

In this paper, we investigate the conformational dynamics of alanine dipeptide under an external electric field by nonequilibrium molecular dynamics simulation. We consider the case of a constant and of an oscillatory field. In this context we propose a procedure to implement the temperature control, which removes the irrelevant thermal effects of the field. For the constant field different time-scales are identified in the conformational, dipole moment, and orientational dynamics. Moreover, we prove that the solvent structure only marginally changes when the external field is switched on. In the case of oscillatory field, the conformational changes are shown to be as strong as in the previous case, and non-trivial nonequilibrium circular paths in the conformation space are revealed by calculating the integrated net probability fluxes.Comment: 23 pages, 12 figure

HUWE1 E3 ligase promotes PINK1/PARKINindependent mitophagy by regulating AMBRA1 activation via IKKa

The selective removal of undesired or damaged mitochondria by autophagy, known as mitophagy, is crucial for cellular homoeostasis, and prevents tumour diffusion, neurodegeneration and ageing. The pro-autophagic molecule AMBRA1 (autophagy/beclin-1 regulator-1) has been defined as a novel regulator of mitophagy in both PINK1/PARKIN-dependent and -independent systems. Here, we identified the E3 ubiquitin ligase HUWE1 as a key inducing factor in AMBRA1-mediated mitophagy, a process that takes place independently of the main mitophagy receptors. Furthermore, we show that mitophagy function of AMBRA1 is post-translationally controlled, upon HUWE1 activity, by a positive phosphorylation on its serine 1014. This modification is mediated by the IKKα kinase and induces structural changes in AMBRA1, thus promoting its interaction with LC3/GABARAP (mATG8) proteins and its mitophagic activity. Altogether, these results demonstrate that AMBRA1 regulates mitophagy through a novel pathway, in which HUWE1 and IKKα are key factors, shedding new lights on the regulation of mitochondrial quality control and homoeostasis in mammalian cells

PIP2-Binding Site in Kir Channels: Definition by Multiscale Biomolecular Simulations†

Phosphatidylinositol bisphosphate (PIP(2)) is an activator of mammalian inwardly rectifying potassium (Kir) channels. Multiscale simulations, via a sequential combination of coarse-grained and atomistic molecular dynamics, enabled exploration of the interactions of PIP(2) molecules within the inner leaflet of a lipid bilayer membrane with possible binding sites on Kir channels. Three Kir channel structures were investigated: X-ray structures of KirBac1.1 and of a Kir3.1-KirBac1.3 chimera and a homology model of Kir6.2. Coarse-grained simulations of the Kir channels in PIP(2)-containing lipid bilayers identified the PIP(2)-binding site on each channel. These models of the PIP(2)-channel complexes were refined by conversion to an atomistic representation followed by molecular dynamics simulation in a lipid bilayer. All three channels were revealed to contain a conserved binding site at the N-terminal end of the slide (M0) helix, at the interface between adjacent subunits of the channel. This binding site agrees with mutagenesis data and is in the proximity of the site occupied by a detergent molecule in the Kir chimera channel crystal. Polar contacts in the coarse-grained simulations corresponded to long-lived electrostatic and H-bonding interactions between the channel and PIP(2) in the atomistic simulations, enabling identification of key side chains

Atomic-Resolution Simulations Predict a Transition State for Vesicle Fusion Defined by Contact of a Few Lipid Tails

Membrane fusion is essential to both cellular vesicle trafficking and infection by enveloped viruses. While the fusion protein assemblies that catalyze fusion are readily identifiable, the specific activities of the proteins involved and nature of the membrane changes they induce remain unknown. Here, we use many atomic-resolution simulations of vesicle fusion to examine the molecular mechanisms for fusion in detail. We employ committor analysis for these million-atom vesicle fusion simulations to identify a transition state for fusion stalk formation. In our simulations, this transition state occurs when the bulk properties of each lipid bilayer remain in a lamellar state but a few hydrophobic tails bulge into the hydrophilic interface layer and make contact to nucleate a stalk. Additional simulations of influenza fusion peptides in lipid bilayers show that the peptides promote similar local protrusion of lipid tails. Comparing these two sets of simulations, we obtain a common set of structural changes between the transition state for stalk formation and the local environment of peptides known to catalyze fusion. Our results thus suggest that the specific molecular properties of individual lipids are highly important to vesicle fusion and yield an explicit structural model that could help explain the mechanism of catalysis by fusion proteins

Beyond the Binding Site: The Role of the β2 – β3 Loop and Extra-Domain Structures in PDZ Domains

A general paradigm to understand protein function is to look at properties of isolated well conserved domains, such as SH3 or PDZ domains. While common features of domain families are well understood, the role of subtle differences among members of these families is less clear. Here, molecular dynamics simulations indicate that the binding mechanism in PSD95-PDZ3 is critically regulated via interactions outside the canonical binding site, involving both the poorly conserved loop and an extra-domain helix. Using the CRIPT peptide as a prototypical ligand, our simulations suggest that a network of salt-bridges between the ligand and this loop is necessary for binding. These contacts interconvert between each other on a time scale of a few tens of nanoseconds, making them elusive to X-ray crystallography. The loop is stabilized by an extra-domain helix. The latter influences the global dynamics of the domain, considerably increasing binding affinity. We found that two key contacts between the helix and the domain, one involving the loop, provide an atomistic interpretation of the increased affinity. Our analysis indicates that both extra-domain segments and loosely conserved regions play critical roles in PDZ binding affinity and specificity

Conformational Changes and Slow Dynamics through Microsecond Polarized Atomistic Molecular Simulation of an Integral Kv1.2 Ion Channel

Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins—not to mention numerous applications in drug design. Here, we present a full 1 µs atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120° rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation (∼35°) of the extracellular end of all S4 segments is present also in a reference 0.5 µs simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 310 helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4–lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5–1 µs). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations