Functional diversity among cardiolipin binding sites on the mitochondrial ADP/ATP carrier.

Lipid-protein interactions play a multitude of essential roles in membrane homeostasis. Mitochondrial membranes have a unique lipid-protein environment that ensures bioenergetic efficiency. Cardiolipin (CL), the signature mitochondrial lipid, plays multiple roles in promoting oxidative phosphorylation (OXPHOS). In the inner mitochondrial membrane, the ADP/ATP carrier (AAC in yeast; adenine nucleotide translocator, ANT in mammals) exchanges ADP and ATP, enabling OXPHOS. AAC/ANT contains three tightly bound CLs, and these interactions are evolutionarily conserved. Here, we investigated the role of these buried CLs in AAC/ANT using a combination of biochemical approaches, native mass spectrometry, and molecular dynamics simulations. We introduced negatively charged mutations into each CL-binding site of yeast Aac2 and established experimentally that the mutations disrupted the CL interactions. While all mutations destabilized Aac2 tertiary structure, transport activity was impaired in a binding site-specific manner. Additionally, we determined that a disease-associated missense mutation in one CL-binding site in human ANT1 compromised its structure and transport activity, resulting in OXPHOS defects. Our findings highlight the conserved significance of CL in AAC/ANT structure and function, directly tied to specific lipid-protein interactions.

Inside the capsid: Revealing viral genome organization through multiscale simulations.

In a recent issue of Nature, Coshic et al. employ a computational multiscale approach to package the complete HK97 viral genome into its capsid. They find both good agreement with experimental observations and shed new light on the heterogeneity of genome structures and the mechanism by which they package.

The α C helix is a central regulator of PKR activation.

Protein kinase R (PKR) functions in the eukaryotic innate immune system as a first-line defense against viral infections. PKR binds viral dsRNA, leading to autophosphorylation and activation. In its active state, PKR can phosphorylate its primary substrate, eIF2 α , which blocks initiation of translation in the infected cell. It has been established that PKR activation occurs when the kinase domain dimerizes in a back-to-back configuration. However, the mechanism by which dimerization leads to enzymatic activation is not fully understood. Here, we investigate the structural mechanistic basis and energy landscape for PKR activation, with a focus on the α C helix - a kinase activation and signal integration hub - using all-atom equilibrium and enhanced sampling molecular dynamics simulations. By employing window-exchange umbrella sampling, we compute free energy profiles of activation which show that back-to-back dimerization stabilizes a catalytically competent conformation of PKR. Key hydrophobic residues in the homodimer interface contribute to stabilization of the α C helix in an active conformation and the position of its glutamate residue. Using linear mutual information analysis, we analyze allosteric communication connecting the protomers' N-lobes and the α C helix dimer interface with the α C helix.

Curvature sensing lipid dynamics in a mitochondrial inner membrane model.

Membrane curvature is essential for many cellular structures and processes, and factors such as leaflet asymmetry, lipid composition, and proteins all play important roles. Cardiolipin is the signature lipid of mitochondrial membranes and is essential for maintaining the highly curved shapes of the inner mitochondrial membrane (IMM) and the spatial arrangement of membrane proteins. In this study, we investigate the partitioning behavior of various lipids present in the IMM using coarse-grained molecular dynamics simulations. This study explores curved bilayer systems containing phosphatidylcholine (PC), phosphatidylethanolamine (PE), and cardiolipin (CDL) in binary and ternary component mixtures. Curvature properties such as mean and Gaussian curvatures, as well as the distribution of lipids into the various curved regions of the cristae models, are quantified. Overall, this work represents an advance beyond previous studies on lipid curvature sensing by simulating these systems in a geometry that has the morphological features and scales of curvature consistent with regions of the IMM. We find that CDL has a stronger preference for accumulating in regions of negative curvature than PE lipids, in agreement with previous results. Furthermore, we find lipid partitioning propensity is dominated by sensitivity to mean curvature, while there is a weaker correlation with Gaussian curvature.

Peptide Power: Mechanistic Insights into the Effect of Mitochondria-Targeted Tetrapeptides on Membrane Electrostatics from Molecular Simulations.

Mitochondrial dysfunction is implicated in nine of the ten leading causes of death in the US, yet there are no FDA-approved therapeutics to treat it. Synthetic mitochondria-targeted peptides (MTPs), including the lead compound SS-31, offer promise, as they have been shown to restore healthy mitochondrial function and treat a variety of common diseases. At the cellular level, research has shown that MTPs accumulate strongly at the inner mitochondrial membrane (IMM), slow energy sinks (e.g., proton leaks), and improve ATP production. Modulation of electrostatic fields around the IMM has been implicated as a key aspect in the mechanism of action (MoA) of these peptides; however, molecular and mechanistic details have remained elusive. In this study, we employed all-atom molecular dynamics simulations (MD) to investigate the interactions of four MTPs with lipid bilayers and calculate their effect on structural and electrostatic properties. In agreement with previous experimental findings, we observed the modulation of the membrane surface and dipole potentials by MTPs. The simulations reveal that the MTPs achieve a reduction in the dipole potential by acting to disorder both lipid head groups and water layers proximal to the bilayer surface. We also find that MTPs decrease the bilayer thickness and increase the membrane's capacitance. These changes suggest that MTPs may enhance how much potential energy can be stored across the IMM at a given transmembrane potential difference. The MTPs also displace cations away from the bilayer surface, modulating the surface potential and offering an alternative mechanism for how these MTPs reduce mitochondrial energy sinks like proton leaks and mitigate Ca2+ accumulation stress. In conclusion, this study highlights the therapeutic potential of MTPs and underlines how interactions of MTPs with lipid bilayers serve as a fundamental component of their MoA.

Conserved cardiolipin-mitochondrial ADP/ATP carrier interactions assume distinct structural and functional roles that are clinically relevant.

The mitochondrial phospholipid cardiolipin (CL) promotes bioenergetics via oxidative phosphorylation (OXPHOS). Three tightly bound CLs are evolutionarily conserved in the ADP/ATP carrier (AAC in yeast; adenine nucleotide translocator, ANT in mammals) which resides in the inner mitochondrial membrane and exchanges ADP and ATP to enable OXPHOS. Here, we investigated the role of these buried CLs in the carrier using yeast Aac2 as a model. We introduced negatively charged mutations into each CL-binding site of Aac2 to disrupt the CL interactions via electrostatic repulsion. While all mutations disturbing the CL-protein interaction destabilized Aac2 monomeric structure, transport activity was impaired in a pocket-specific manner. Finally, we determined that a disease-associated missense mutation in one CL-binding site in ANT1 compromised its structure and transport activity, resulting in OXPHOS defects. Our findings highlight the conserved significance of CL in AAC/ANT structure and function, directly tied to specific lipid-protein interactions.

Non-equilibrium virus particle dynamics: Microsecond MD simulations of the complete Flock House virus capsid under different conditions.

Flock House virus (FHV) is an animal virus and considered a model system for non-enveloped viruses. It has a small, icosahedral capsid (T=3) and a bipartite positive-sense RNA genome. We present an extensive study of the FHV capsid dynamics from all-atom molecular dynamics simulations of the complete capsid. The simulations explore different biologically relevant conditions (neutral/low pH, with/without RNA in the capsid) using the CHARMM force field. The results show that low pH destabilizes the capsid, causing radial expansion, and RNA stabilizes the capsid. The finding of low pH destabilization is biologically relevant because the capsid is exposed to low pH in the endosome, where conformational changes occur leading to genome release. We also observe structural changes at the fivefold and twofold symmetry axes that likely relate to the externalization of membrane active γ peptides through the fivefold vertex and extrusion of RNA at the twofold axis. Simulations using the Amber force field at neutral pH are also performed and display similar characteristics to the CHARMM simulations.

Structure-activity relationships of mitochondria-targeted tetrapeptide pharmacological compounds.

Mitochondria play a central role in metabolic homeostasis, and dysfunction of this organelle underpins the etiology of many heritable and aging-related diseases. Tetrapeptides with alternating cationic and aromatic residues such as SS-31 (elamipretide) show promise as therapeutic compounds for mitochondrial disorders. In this study, we conducted a quantitative structure-activity analysis of three alternative tetrapeptide analogs, benchmarked against SS-31, that differ with respect to aromatic side chain composition and sequence register. We present the first structural models for this class of compounds, obtained with Nuclear Magnetic Resonance (NMR) and molecular dynamics approaches, showing that all analogs except for SS-31 form compact reverse turn conformations in the membrane-bound state. All peptide analogs bound cardiolipin-containing membranes, yet they had significant differences in equilibrium binding behavior and membrane interactions. Notably, analogs had markedly different effects on membrane surface charge, supporting a mechanism in which modulation of membrane electrostatics is a key feature of their mechanism of action. The peptides had no strict requirement for side chain composition or sequence register to permeate cells and target mitochondria in mammalian cell culture assays. All four peptides were pharmacologically active in serum withdrawal cell stress models yet showed significant differences in their abilities to restore mitochondrial membrane potential, preserve ATP content, and promote cell survival. Within our peptide set, the analog containing tryptophan side chains, SPN10, had the strongest impact on most membrane properties and showed greatest efficacy in cell culture studies. Taken together, these results show that side chain composition and register influence the activity of these mitochondria-targeted peptides, helping provide a framework for the rational design of next-generation therapeutics with enhanced potency.

A favorable path to domain separation in the orange carotenoid protein.

The orange carotenoid protein (OCP) is responsible for nonphotochemical quenching (NPQ) in cyanobacteria, a defense mechanism against potentially damaging effects of excess light conditions. This soluble two-domain protein undergoes profound conformational changes upon photoactivation, involving translocation of the ketocarotenoid inside the cavity followed by domain separation. Domain separation is a critical step in the photocycle of OCP because it exposes the N-terminal domain (NTD) to perform quenching of the phycobilisomes. Many details regarding the mechanism and energetics of OCP domain separation remain unknown. In this work, we apply metadynamics to elucidate the protein rearrangements that lead to the active, domain-separated, form of OCP. We find that translocation of the ketocarotenoid canthaxanthin has a profound effect on the energetic landscape and that domain separation only becomes favorable following translocation. We further explore, characterize, and validate the free energy surface (FES) using equilibrium simulations initiated from different states on the FES. Through pathway optimization methods, we characterize the most probable path to domain separation and reveal the barriers along that pathway. We find that the free energy barriers are relatively small (<5 kcal/mol), but the overall estimated kinetic rate is consistent with experimental measurements (>1 ms). Overall, our results provide detailed information on the requirement for canthaxanthin translocation to precede domain separation and an energetically feasible pathway to dissociation.

Atomistic dynamics of a viral infection process: Release of membrane lytic peptides from a non-enveloped virus.

Molecular simulations have played an instrumental role in uncovering the structural dynamics and physical properties of virus capsids. In this work, we move beyond equilibrium physicochemical characterization of a virus system to study a stage of the infection process that is required for viral proliferation. Despite many biochemical and functional studies, the molecular mechanism of host cell entry by non-enveloped viruses remains largely unresolved. Flock House virus (FHV) is a model system for non-enveloped viruses and is the subject of the current study. FHV infects through the acid-dependent endocytic pathway, where low pH triggers externalization of membrane-disrupting (γ) peptides from the capsid interior. Using all-atom equilibrium and enhanced sampling simulations, the mechanism and energetics of γ peptide liberation and the effect of pH on this process are investigated. Our computations agree with experimental findings and reveal nanoscopic details regarding the pH control mechanism, which are not readily accessible in experiments.

Structural and dynamic asymmetry in icosahedrally symmetric virus capsids.

A common characteristic of virus capsids is icosahedral symmetry, yet these highly symmetric structures can display asymmetric features within their virions and undergo asymmetric dynamics. The fields of structural and computational biology have entered a new realm in the investigation of virus infection mechanisms, with the ability to observe symmetry-breaking features. This review will cover important studies on icosahedral virus structure and dynamics, covering both symmetric and asymmetric conformational changes. However, the main emphasis of the review will be towards recent studies employing cryo-electron microscopy or molecular dynamics simulations, which can uncover asymmetric aspects of these systems relevant to understanding viral physical-chemical properties and their biological impact.

The mitochondria-targeted peptide SS-31 binds lipid bilayers and modulates surface electrostatics as a key component of its mechanism of action.

Mitochondrial dysfunction underlies many heritable diseases, acquired pathologies, and aging-related declines in health. Szeto-Schiller (SS) peptides comprise a class of amphipathic tetrapeptides that are efficacious toward a wide array of mitochondrial disorders and are believed to target mitochondrial membranes because they are enriched in the anionic phospholipid cardiolipin (CL). However, little is known regarding how SS peptides interact with or alter the physical properties of lipid bilayers. In this study, using biophysical and computational approaches, we have analyzed the interactions of the lead compound SS-31 (elamipretide) with model and mitochondrial membranes. Our results show that this polybasic peptide partitions into the membrane interfacial region with an affinity and a lipid binding density that are directly related to surface charge. We found that SS-31 binding does not destabilize lamellar bilayers even at the highest binding concentrations; however, it did cause saturable alterations in lipid packing. Most notably, SS-31 modulated the surface electrostatics of both model and mitochondrial membranes. We propose nonexclusive mechanisms by which the tuning of surface charge could underpin the mitoprotective properties of SS-31, including alteration of the distribution of ions and basic proteins at the interface, and/or modulation of bilayer physical properties. As a proof of concept, we show that SS-31 alters divalent cation (calcium) distribution within the interfacial region and reduces the energetic burden of calcium stress in mitochondria. The mechanistic details of SS-31 revealed in this study will help inform the development of future compound variants with enhanced efficacy and bioavailability.

Markov State Model of Lassa Virus Nucleoprotein Reveals Large Structural Changes during the Trimer to Monomer Transition.

Lassa virus contains a nucleoprotein (NP) that encapsulates the viral genomic RNA forming the ribonucleoprotein (RNP). The NP forms trimers that do not bind RNA, but a structure of only the NP N-terminal domain was co-crystallized with RNA bound. These structures suggested a model in which the NP forms a trimer to keep the RNA gate closed, but then is triggered to undergo a change to a form competent for RNA binding. Here, we investigate the scenario in which the trimer is disrupted to observe whether monomeric NP undergoes significant conformational changes. From multi-microsecond molecular dynamics simulations and an adaptive sampling scheme to sample the conformational space, a Markov state model (MSM) is constructed. The MSM reveals an energetically favorable conformational change, with the most significant changes occurring at the domain interface. These results support a model in which significant structural reorganization of the NP is required for RNP formation.

A Potent Host Defense Peptide Triggers DNA Damage and Is Active against Multidrug-Resistant Gram-Negative Pathogens.

Gram-negative bacteria are some of the biggest threats to public health due to a large prevalence of antibiotic resistance. The difficulty in treating bacterial infections, stemming from their double membrane structure combined with efflux pumps in the outer membrane, has resulted in a much greater need for antimicrobials with activity against these pathogens. Tunicate host defense peptide (HDP), Clavanin A, is capable of not only inhibiting Gram-negative growth but also potentiating activity in the presence of Zn(II). Here, we provide evidence that the improvements of Clavanin A activity in the presence of Zn(II) are due to its novel mechanism of action. We employed E. coli TD172 (ΔrecA::kan) and the terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) assay to show in cellulae that DNA damage occurs upon treatment with Clavanin A. In vitro assays demonstrated that Zn(II) ions are required for the nuclease activity of the peptide. The quantum mechanics/molecular mechanics (QM/MM) calculations were used to investigate the mechanism of DNA damage. In the rate-determining step of the proposed mechanism, due to its Lewis acidity, the Zn(II) ion activates the scissile P-O bond of DNA and creates a hydroxyl nucleophile from a water molecule. A subsequent attack by this group to the electrophilic phosphorus cleaves the scissile phosphoester bond. Additionally, we utilized bacterial cytological profiling (BCP), circular dichroism (CD) spectroscopy in the presence of lipid vesicles, and surface plasmon resonance combined with electrical impedance spectroscopy in order to address the apparent discrepancies between our results and the previous studies regarding the mechanism of action of Clavanin A. Finally, our approach may lead to the identification of additional Clavanin A like HDPs and promote the development of antimicrobial peptide based therapeutics.

Molecular dynamics study of membrane permeabilization by wild-type and mutant lytic peptides from the non-enveloped Flock House virus.

Flock House virus (FHV) serves as a model system for understanding infection mechanisms utilized by non-enveloped viruses to transport across cellular membranes. During the infection cycle of FHV, a fundamental stage involves disruption of the endosomal membrane by membrane active peptides, following externalization of the peptides from the capsid interior. The FHV lytic agents are the 44 C-terminal amino acids residues of the capsid protein, which are auto-catalytically cleaved during the capsid maturation process. The cleaved peptides are termed γ peptides. In this study, we perform multi-scale molecular dynamics simulations including 40 µs all-atom molecular dynamics simulations to study the behavior of pre-inserted transmembrane lytic peptides at a high concentration in a neutral membrane. We study the dynamical organization among peptides to form oligomeric bundles in four systems including the wild-type γ peptide and three mutant forms; namely, a truncation mutant in which the 23 C-terminal residues are deleted (γ1), a construct where the 8 C-terminal residues of γ are fused to γ1 (Δ385-399 γ) and a single-point mutant (F402A γ), all of which have been experimentally shown to drastically affect infectivity and lytic activity compared to the wild-type γ. Our results shed light on the actions of varied forms of the FHV lytic peptide including membrane insertion, trans-membrane stability, peptide oligomerization, water permeation activity and dynamic pore formation. Findings from this study provide detailed structural information and rationale for the differences in lytic activity among variants of FHV γ.

Functional and structural basis of E. coli enolase inhibition by SF2312: a mimic of the carbanion intermediate.

Many years ago, the natural secondary metabolite SF2312, produced by the actinomycete Micromonospora, was reported to display broad spectrum antibacterial properties against both Gram-positive and Gram-negative bacteria. Recent studies have revealed that SF2312, a natural phosphonic acid, functions as a potent inhibitor of human enolase. The mechanism of SF2312 inhibition of bacterial enolase and its role in bacterial growth and reproduction, however, have remained elusive. In this work, we detail a structural analysis of E. coli enolase bound to both SF2312 and its oxidized imide-form. Our studies support a model in which SF2312 acts as an analog of a high energy intermediate formed during the catalytic process. Biochemical, biophysical, computational and kinetic characterization of these compounds confirm that altering features characteristic of a putative carbanion (enolate) intermediate significantly reduces the potency of enzyme inhibition. When SF2312 is combined with fosfomycin in the presence of glucose-6 phosphate, significant synergy is observed. This suggests the two agents could be used as a potent combination, targeting distinct cellular mechanism for the treatment of bacterial infections. Together, our studies rationalize the structure-activity relationships for these phosphonates and validate enolase as a promising target for antibiotic discovery.

Computational Studies of Catalytic Loop Dynamics in Yersinia Protein Tyrosine Phosphatase Using Pathway Optimization Methods.

Yersinia Protein Tyrosine Phosphatase (YopH) is the most efficient enzyme among all known PTPases and relies on its catalytic loop movements for substrate binding and catalysis. Fluorescence, NMR, and UV resonance Raman (UVRR) techniques have been used to study the thermodynamic and dynamic properties of the loop motions. In this study, a computational approach based on the pathway refinement methods nudged elastic band (NEB) and harmonic Fourier beads (HFB) has been developed to provide structural interpretations for the experimentally observed kinetic processes. In this approach, the minimum potential energy pathways for the loop open/closure conformational changes were determined by NEB using a one-dimensional global coordinate. Two dimensional data analyses of the NEB results were performed as an efficient method to qualitatively evaluate the energetics of transitions along several specific physical coordinates. The free energy barriers for these transitions were then determined more precisely using the HFB method. Kinetic parameters were estimated from the energy barriers using transition state theory and compared against experimentally determined kinetic parameters. When the calculated energy barriers are calibrated by a simple "scaling factor", as have been done in our previous vibrational frequency calculations to explain the ligand frequency shift upon its binding to protein, it is possible to make structural interpretations of several observed enzyme dynamic rates. For example, the nanosecond kinetics observed by fluorescence anisotropy may be assigned to the translational motion of the catalytic loop and microsecond kinetics observed in fluorescence T-jump can be assigned to the loop backbone dihedral angle flipping. Furthermore, we can predict that a Trp354 conformational conversion associated with the loop movements would occur on the tens of nanoseconds time scale, to be verified by future UVRR T-jump studies.

Structural Basis of Protein Kinase R Autophosphorylation.

The RNA-activated protein kinase, PKR, is a key mediator of the innate immunity response to viral infection. Viral double-stranded RNAs induce PKR dimerization and autophosphorylation. The PKR kinase domain forms a back-to-back dimer. However, intermolecular ( trans) autophosphorylation is not feasible in this arrangement. We have obtained PKR kinase structures that resolves this dilemma. The kinase protomers interact via the known back-to-back interface as well as a front-to-front interface that is formed by exchange of activation segments. Mutational analysis of the front-to-front interface support a functional role in PKR activation. Molecular dynamics simulations reveal that the activation segment is highly dynamic in the front-to-front dimer and can adopt conformations conducive to phosphoryl transfer. We propose a mechanism where back-to-back dimerization induces a conformational change that activates PKR to phosphorylate a "substrate" kinase docked in a front-to-front geometry. This mechanism may be relevant to related kinases that phosphorylate the eukaryotic initiation factor eIF2α.

Molecular Dynamics Investigation into the Effect of Zinc(II) on the Structure and Membrane Interactions of the Antimicrobial Peptide Clavanin A.

Clavanin A (ClavA) is an antimicrobial peptide (AMP) whose antimicrobial activity is enhanced in the presence of Zn(II) ions. The antimicrobial activity of ClavA has been shown to increase 16-fold in the presence of Zn(II) ions. In this study, we investigate the potential sources of this enhancement, namely, the effect of Zn(II) binding on the helical conformation of ClavA and on the ClavA interaction with a model for gram-negative bacterial membranes. In addition, we investigate the effect of Zn(II) on the membrane mechanical properties. We employed all-atom equilibrium molecular dynamics simulations initiated from both fully helical and random coil structures of ClavA. We observe that Zn(II) can stabilize an existing helical conformation in the Zn(II)-binding region, but we do not observe induction of helical conformations in systems initiated in random coil configurations. Zn(II) binding to ClavA provides more favorable electrostatics for membrane association in the C-terminal region. This is evidenced by longer and stronger C-terminal-lipid interactions. Zn(II) is also capable of modulating the membrane properties in a manner which favors ClavA insertion and the potential for enhanced translocation into the cell. This work provides insights into the role of divalent metal cations in the antimicrobial activity of ClavA. This information can be used for the development of synthetic AMPs containing motifs that can bind metals (metalloAMPs) for therapeutic and medical purposes.