Hydrogen bonding heterogeneity correlates with protein folding transition state passage time as revealed by data sonification.

Protein-protein and protein-water hydrogen bonding interactions play essential roles in the way a protein passes through the transition state during folding or unfolding, but the large number of these interactions in molecular dynamics (MD) simulations makes them difficult to analyze. Here, we introduce a state space representation and associated "rarity" measure to identify and quantify transition state passage (transit) events. Applying this representation to a long MD simulation trajectory that captured multiple folding and unfolding events of the GTT WW domain, a small protein often used as a model for the folding process, we identified three transition categories: Highway (faster), Meander (slower), and Ambiguous (intermediate). We developed data sonification and visualization tools to analyze hydrogen bond dynamics before, during, and after these transition events. By means of these tools, we were able to identify characteristic hydrogen bonding patterns associated with "Highway" versus "Meander" versus "Ambiguous" transitions and to design algorithms that can identify these same folding pathways and critical protein-water interactions directly from the data. Highly cooperative hydrogen bonding can either slow down or speed up transit. Furthermore, an analysis of protein-water hydrogen bond dynamics at the surface of WW domain shows an increase in hydrogen bond lifetime from folded to unfolded conformations with Ambiguous transitions as an outlier. In summary, hydrogen bond dynamics provide a direct window into the heterogeneity of transits, which can vary widely in duration (by a factor of 10) due to a complex energy landscape.

Metastable States in the Hinge-Bending Landscape of an Enzyme in an Atomistic Cytoplasm Simulation.

Many enzymes undergo major conformational changes to function in cells, particularly when they bind to more than one substrate. We quantify the large-amplitude hinge-bending landscape of human phosphoglycerate kinase (PGK) in a human cytoplasm. Approximately 70 µs of all-atom simulations, upon coarse graining, reveal three metastable states of PGK with different hinge angle distributions and additional substates. The "open" state was more populated than the "semi-open" or "closed" states. In addition to free energies and barriers within the landscape, we characterized the average transition state passage time of ≈0.3 µs and reversible substrate and product binding. Human PGK in a dilute solution simulation shows a transition directly from the open to closed states, in agreement with previous SAXS experiments, suggesting that the cell-like model environment promotes stability of the human PGK semi-open state. Yeast PGK also sampled three metastable states within the cytoplasm model, with the closed state favored in our simulation.

Tuning sterol extraction kinetics yields a renal-sparing polyene antifungal.

Decades of previous efforts to develop renal-sparing polyene antifungals were misguided by the classic membrane permeabilization model1. Recently, the clinically vital but also highly renal-toxic small-molecule natural product amphotericin B was instead found to kill fungi primarily by forming extramembraneous sponge-like aggregates that extract ergosterol from lipid bilayers2-6. Here we show that rapid and selective extraction of fungal ergosterol can yield potent and renal-sparing polyene antifungals. Cholesterol extraction was found to drive the toxicity of amphotericin B to human renal cells. Our examination of high-resolution structures of amphotericin B sponges in sterol-free and sterol-bound states guided us to a promising structural derivative that does not bind cholesterol and is thus renal sparing. This derivative was also less potent because it extracts ergosterol more slowly. Selective acceleration of ergosterol extraction with a second structural modification yielded a new polyene, AM-2-19, that is renal sparing in mice and primary human renal cells, potent against hundreds of pathogenic fungal strains, resistance evasive following serial passage in vitro and highly efficacious in animal models of invasive fungal infections. Thus, rational tuning of the dynamics of interactions between small molecules may lead to better treatments for fungal infections that still kill millions of people annually7,8 and potentially other resistance-evasive antimicrobials, including those that have recently been shown to operate through supramolecular structures that target specific lipids9.

In-Cell Dynamics: The Next Focus of All-Atom Simulations.

The cell is a crowded space where large biomolecules and metabolites are in continuous motion. Great strides have been made in in vitro studies of protein dynamics, folding, and protein-protein interactions, and much new data are emerging of how they differ in the cell. In this Perspective, we highlight the current progress in atomistic modeling of in-cell environments, both bacteria and mammals, with emphasis on classical all-atom molecular dynamics simulations. These simulations have been recently used to capture and characterize functional and non-functional protein-protein interactions, protein folding dynamics of small proteins with varied topologies, and dynamics of metabolites. We further discuss the challenges and efforts for updating modern force fields critical to the progress of cellular environment simulations. We also briefly summarize developments in relevant state-of-the-art experimental techniques. As computational and experimental methodologies continue to progress and produce more directly comparable data, we are poised to capture the complex atomistic picture of the cell.

Lactadherin's Multistate Binding Predicts Stable Membrane-Bound Conformations of Factors V and VIII's C Domains.

Protein binding to negatively charged lipids is essential for maintaining numerous vital cellular processes where its dysfunction can lead to various diseases. One such protein that plays a crucial role in this process is lactadherin, which competes with coagulation factors for membrane binding sites to regulate blood clotting. Despite identifying key binding regions of these proteins through structural and biochemical studies, models incorporating membrane dynamics are still lacking. In this study, we report on the multimodal binding of lactadherin and use it to gain insight into the binding mechanisms of its C domain homologs, factor V and factor VIII. Molecular dynamics simulations enhanced with the highly mobile mimetic model enabled the determination of lactadherin's multimodal binding on membranes that revealed critical interacting residues consistent with prior NMR and mutagenesis data. The binding occurred primarily via two dynamic structural ensembles: an inserted state and an unreported, highly conserved side-lying state driven by a cationic patch. We utilized these findings to analyze the membrane binding domains of coagulation factors V and VIII and identified their preferred membrane-bound conformations. Specifically, factor V's C domains maintained an inserted state, while factor VIII preferred a tilted, side-lying state that permitted antibody binding. Insight into lactadherin's atomistically resolved membrane interactions from a multistate perspective can guide new therapeutic opportunities in treating diseases related to blood coagulation.

In silico protein dynamics in the human cytoplasm: Partial folding, misfolding, fold switching, and non-native interactions.

We examine the influence of cellular interactions in all-atom models of a section of the Homo sapiens cytoplasm on the early folding events of the three-helix bundle protein B (PB). While genetically engineered PB is known to fold in dilute water box simulations in three microseconds, the three initially unfolded PB copies in our two cytoplasm models using a similar force field did not reach the native state during 30-microsecond simulations. We did however capture the formation of all three helices in a compact native-like topology. Folding in vivo is delayed because intramolecular contact formation within PB is in direct competition with intermolecular contacts between PB and surrounding macromolecules. In extreme cases, intermolecular beta-sheets are formed. Interactions with other macromolecules are also observed to promote structure formation, for example when a PB helix in our simulations is shielded from solvent by macromolecular crowding. Sticking and crowding in our models initiate sampling of helix/sheet structural plasticity of PB. Relatedly, in past in vitro experiments, similar GA domains were shown to switch between two different folds. Finally, we also observed that stickiness between PB and the cellular environment can be modulated in our simulations through the reduction in protein hydrophobicity when we reversed PB back to the wild-type sequence. This study demonstrates that even fast-folding proteins can get stuck in non-native states in the cell, making them useful models for protein-chaperone interactions and early stages of aggregate formation relevant to cellular disease.

Utility of FRET in studies of membrane protein oligomerization: The concept of the effective dissociation constant.

The activity of many membrane receptors is controlled through their lateral association into dimers or higher-order oligomers. Although Förster resonance energy transfer (FRET) measurements have been used extensively to characterize the stability of receptor dimers, the utility of FRET in studies of larger oligomers has been limited. Here we introduce an effective equilibrium dissociation constant that can be extracted from FRET measurements for EphA2, a receptor tyrosine kinase (RTK) known to form active oligomers of heterogeneous distributions in response to its ligand ephrinA1-Fc. The newly introduced effective equilibrium dissociation constant has a well-defined physical meaning and biological significance. It denotes the receptor concentration for which half of the receptors are monomeric and inactive, and the other half are associated into oligomers and are active, irrespective of the exact oligomer size. This work introduces a new dimension to the utility of FRET in studies of membrane receptor association and signaling in the plasma membrane.

Segmental Dynamics of Membranous Cholesterol are Coupled.

Cholesterol promotes the structural integrity of the fluid cell membrane and interacts dynamically with many membrane proteins to regulate function. Understanding site-resolved cholesterol structural dynamics is thus important. This long-standing challenge has thus far been addressed, in part, by selective isotopic labeling approaches. Here we present a new 3D solid-state NMR (SSNMR) experiment utilizing scalar 13C-13C polarization transfer and recoupling of the 1H-13C interactions in order to determine average dipolar couplings for all 1H-13C vectors in uniformly 13C-enriched cholesterol. The experimentally determined order parameters (OP) agree exceptionally well with molecular dynamics (MD) trajectories and reveal coupling among several conformational degrees of freedom in cholesterol molecules. Quantum chemistry shielding calculations further support this conclusion and specifically demonstrate that ring tilt and rotation are coupled to changes in tail conformation and that these coupled segmental dynamics dictate the orientation of cholesterol. These findings advance our understanding of physiologically relevant dynamics of cholesterol, and the methods that revealed them have broader potential to characterize how structural dynamics of other small molecules impact their biological functions.

TMAO: Protecting proteins from feeling the heat.

Osmolytes are ubiquitous in the cell and play an important role in controlling protein stability under stress. The natural osmolyte trimethylamine N-oxide (TMAO) is used by marine animals to counteract the effect of pressure denaturation at large depths. The molecular mechanism of TMAO stabilization against pressure and urea denaturation has been extensively studied, but unlike the case of other osmolytes, the ability of TMAO to protect proteins from high temperature has not been quantified. To reveal the effect of TMAO on folded and unfolded protein ensembles and the hydration shell at different temperatures, we study a mutant of the well-characterized, fast-folding model protein B (PRB). We carried out, in total, >190 µs all-atom simulations of thermal folding/unfolding of PRB at multiple temperatures and concentrations of TMAO. The simulations show increased thermal stability of PRB in the presence of TMAO. Partly structured, compact ensembles are favored over the unfolded state. TMAO forms two shells near the protein: an outer shell away from the protein surface has altered H-bond lifetimes of water molecules and increases hydration of the protein to help stabilize it; a less-populated inner shell with an opposite TMAO orientation closer to the protein surface binds exclusively to basic side chains. The cooperative cosolute effect of the inner and outer shell TMAO has a small number of TMAO molecules "herding" water molecules into two hydration shells at or near the protein surface. The stabilizing effect of TMAO on our protein saturates at 1 M despite higher TMAO solubility, so there may be little evolutionary pressure for extremophiles to produce higher intracellular TMAO concentrations, if true in general.

Enzymatic Pyridine Aromatization during Thiopeptide Biosynthesis.

Thiazole-containing pyritides (thiopeptides) are ribosomally synthesized and post-translationally modified peptides (RiPPs) that have attracted interest owing to their potent biological activities and structural complexity. The class-defining feature of a thiopeptide is a six-membered, nitrogenous heterocycle formed by an enzymatic [4 + 2]-cycloaddition. In rare cases, piperidine or dehydropiperidine (DHP) is present; however, the aromatized pyridine is considerably more common. Despite significant effort, the mechanism by which the central pyridine is formed remains poorly understood. Building on our recent observation of the Bycroft-Gowland intermediate (i.e., the direct product of the [4 + 2]-cycloaddition), we interrogated thiopeptide pyridine synthases using a combination of targeted mutagenesis, kinetic assays, substrate analogs, enzyme-substrate cross-linking, and chemical rescue experiments. Collectively, our data delineate roles for several conserved residues in thiopeptide pyridine synthases. A critical tyrosine facilitates the final aromatization step of pyridine formation. This work provides a foundation for further exploration of the [4 + 2]-cycloaddition reaction and future customization of pyridine-containing macrocyclic peptides.

Impact of the Cellular Environment on Adenosine Triphosphate Conformations.

The cytoplasm is an environment crowded by macromolecules and filled with metabolites and ions. Recent experimental and computational studies have addressed how this environment affects protein stability, folding kinetics, and protein-protein and protein-nucleic acid interactions, though its impact on metabolites remains largely unknown. Here we show how a simulated cytoplasm affects the conformation of adenosine triphosphate (ATP), a key energy source and regulatory metabolite present at high concentrations in cells. Analysis of our all-atom model of a small volume of the Escherichia coli cytoplasm when contrasted with ATP modeled in vitro or resolved with protein structures deposited in the Protein Data Bank reveals that ATP molecules bound to proteins in cell form specific pitched conformations that are not observed at significant concentrations in the other environments. We hypothesize that these interactions evolved to fulfill functional roles when ATP interacts with protein surfaces.

Fungicidal amphotericin B sponges are assemblies of staggered asymmetric homodimers encasing large void volumes.

Amphotericin B (AmB) is a powerful but toxic fungicide that operates via enigmatic small molecule-small molecule interactions. This mechanism has challenged the frontiers of structural biology for half a century. We recently showed AmB primarily forms extramembranous aggregates that kill yeast by extracting ergosterol from membranes. Here, we report key structural features of these antifungal 'sponges' illuminated by high-resolution magic-angle spinning solid-state NMR, in concert with simulated annealing and molecular dynamics computations. The minimal unit of assembly is an asymmetric head-to-tail homodimer: one molecule adopts an all-trans C1-C13 motif, the other a C6-C7-gauche conformation. These homodimers are staggered in a clathrate-like lattice with large void volumes similar to the size of sterols. These results illuminate the atomistic interactions that underlie fungicidal assemblies of AmB and suggest this natural product may form biologically active clathrates that host sterol guests.

Bioinformatics-Guided Expansion and Discovery of Graspetides.

Graspetides are a class of ribosomally synthesized and post-translationally modified peptide natural products featuring ATP-grasp ligase-dependent formation of macrolactones/macrolactams. These modifications arise from serine, threonine, or lysine donor residues linked to aspartate or glutamate acceptor residues. Characterized graspetides include serine protease inhibitors such as the microviridins and plesiocin. Here, we report an update to Rapid ORF Description and Evaluation Online (RODEO) for the automated detection of graspetides, which identified 3,923 high-confidence graspetide biosynthetic gene clusters. Sequence and co-occurrence analyses doubled the number of graspetide groups from 12 to 24, defined based on core consensus sequence and putative secondary modification. Bioinformatic analyses of the ATP-grasp ligase superfamily suggest that extant graspetide synthetases diverged once from an ancestral ATP-grasp ligase and later evolved to introduce a variety of ring connectivities. Furthermore, we characterized thatisin and iso-thatisin, two graspetides related by conformational stereoisomerism from Lysobacter antibioticus. Derived from a newly identified graspetide group, thatisin and iso-thatisin feature two interlocking macrolactones with identical ring connectivity, as determined by a combination of tandem mass spectrometry (MS/MS), methanolytic, and mutational analyses. NMR spectroscopy of thatisin revealed a cis conformation for a key proline residue, while molecular dynamics simulations, solvent-accessible surface area calculations, and partial methanolytic analysis coupled with MS/MS support a trans conformation for iso-thatisin at the same position. Overall, this work provides a comprehensive overview of the graspetide landscape, and the improved RODEO algorithm will accelerate future graspetide discoveries by enabling open-access analysis of existing and emerging genomes.

A cancer mutation promotes EphA4 oligomerization and signaling by altering the conformation of the SAM domain.

The Eph receptor tyrosine kinases and their ephrin ligands regulate many physiological and pathological processes. EphA4 plays important roles in nervous system development and adult homeostasis, while aberrant EphA4 signaling has been implicated in neurodegeneration. EphA4 may also affect cancer malignancy, but the regulation and effects of EphA4 signaling in cancer are poorly understood. A correlation between decreased patient survival and high EphA4 mRNA expression in melanoma tumors that also highly express ephrinA ligands suggests that enhanced EphA4 signaling may contribute to melanoma progression. A search for EphA4 gain-of-function mutations in melanoma uncovered a mutation of the highly conserved leucine 920 in the EphA4 sterile alpha motif (SAM) domain. We found that mutation of L920 to phenylalanine (L920F) potentiates EphA4 autophosphorylation and signaling, making it the first documented EphA4 cancer mutation that increases kinase activity. Quantitative Föster resonance energy transfer and fluorescence intensity fluctuation (FIF) analyses revealed that the L920F mutation induces a switch in EphA4 oligomer size, from a dimer to a trimer. We propose this switch in oligomer size as a novel mechanism underlying EphA4-linked tumorigenesis. Molecular dynamics simulations suggest that the L920F mutation alters EphA4 SAM domain conformation, leading to the formation of EphA4 trimers that assemble through two aberrant SAM domain interfaces. Accordingly, EphA4 wild-type and the L920F mutant are affected differently by the SAM domain and are differentially regulated by ephrin ligand stimulation. The increased EphA4 activation induced by the L920F mutation, through the novel mechanism we uncovered, supports a functional role for EphA4 in promoting pathogenesis.

Discovery of Prenyltransferase Inhibitors with In Vitro and In Vivo Antibacterial Activity.

Cis-prenyltransferases such as undecaprenyl diphosphate synthase (UPPS) and decaprenyl diphosphate synthase (DPPS) are essential enzymes in bacteria and are involved in cell wall biosynthesis. UPPS and DPPS are absent in the human genome, so they are of interest as targets for antibiotic development. Here, we screened a library of 750 compounds from National Cancer Institute Diversity Set V for the inhibition of Mycobacterium tuberculosis DPPS and found 17 hits, and then IC50s were determined using dose-response curves. Compounds were tested for growth inhibition against a panel of bacteria, for in vivo activity in a Staphylococcus aureus/Caenorhabditis elegans model, and for mammalian cell toxicity. The most active DPPS inhibitor was the dicarboxylic acid redoxal (compound 10), which also inhibited undecaprenyl diphosphate synthase (UPPS) as well as farnesyl diphosphate synthase. 10 was active against S. aureus, Clostridiodes difficile, Bacillus anthracis Sterne, and Bacillus subtilis, and there was a 3.4-fold increase in IC50 on addition of a rescue agent, undecaprenyl monophosphate. We found that 10 was also a weak protonophore uncoupler, leading to the idea that it targets both isoprenoid biosynthesis and the proton motive force. In an S. aureus/C. elegans in vivo model, 10 reduced the S. aureus burden 3 times more effectively than did ampicillin.

Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism.

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

Coaction of Electrostatic and Hydrophobic Interactions: Dynamic Constraints on Disordered TrkA Juxtamembrane Domain.

In the receptor tyrosine kinase family, conformational change induced by ligand binding is transmitted across the membrane via a single transmembrane helix and a flexible juxtamembrane domain (JMD). Membrane dynamics makes it challenging to study the structural mechanism of receptor activation experimentally. In this study, we employ all-atom molecular dynamics with highly mobile membrane mimetic (HMMM) to capture the native conformation of the JMD in tropomyosin receptor kinase A (TrkA). We find that phosphatidylinositol 4,5-bisphosphate (PIP2) lipids engage in stable binding with multiple basic residues. Anionic lipids can compete with salt bridges within the peptide and alter TrkA-JMD conformation. We discover three-residue insertion into the membrane and are able to either enhance or reduce the level of insertion through computationally-designed point mutations. The vesicle-binding experiment supports computational results and indicates that hydrophobic insertion is comparable to electrostatic binding for membrane anchoring. Biochemical assays on cell lines with mutated TrkA show that enhanced TrkA-JMD insertion promotes receptor degradation but does not affect the short-term signaling capacity. Our joint work points to a scenario where lipid headgroups and tails interact with basic and hydrophobic residues on disordered domain, respectively, to restrain flexibility and potentially modulate protein function.

In-Cell Protein-Protein Contacts: Transient Interactions in the Crowd.

Proteins in vivo are immersed in a crowded environment of water, ions, metabolites, and macromolecules. In-cell experiments highlight how transient weak protein-protein interactions promote (via functional "quinary structure") or hinder (via competitive binding or "sticking") complex formation. Computational models of the cytoplasm are expensive. We tackle this challenge with an all-atom model of a small volume of the E. coli cytoplasm to simulate protein-protein contacts up to the 5 µs time scale on the special-purpose supercomputer Anton 2. We use three CHARMM-derived force fields: C22*, C36m, and C36mCU (with CUFIX corrections). We find that both C36m and C36mCU form smaller contact surfaces than C22*. Although CUFIX was developed to reduce protein-protein sticking, larger contacts are observed with C36mCU than C36m. We show that the lifespan Δt of protein-protein contacts obeys a power law distribution between 0.03 and 3 µs, with ∼90% of all contacts lasting <1 µs (similar to the time scale for downhill folding).

Bioinformatic Mapping of Radical S-Adenosylmethionine-Dependent Ribosomally Synthesized and Post-Translationally Modified Peptides Identifies New Cα, Cß, and Cγ-Linked Thioether-Containing Peptides.

Recently developed bioinformatic tools have bolstered the discovery of ribosomally synthesized and post-translationally modified peptides (RiPPs). Using an improved version of Rapid ORF Description and Evaluation Online (RODEO 2.0), a biosynthetic gene cluster mining algorithm, we bioinformatically mapped the sactipeptide RiPP class via the radical S-adenosylmethionine (SAM) enzymes that form the characteristic sactionine (sulfur-to-α carbon) cross-links between cysteine and acceptor residues. Hundreds of new sactipeptide biosynthetic gene clusters were uncovered, and a novel sactipeptide "huazacin" with growth-suppressive activity against Listeria monocytogenes was characterized. Bioinformatic analysis further suggested that a group of sactipeptide-like peptides heretofore referred to as six cysteines in forty-five residues (SCIFFs) might not be sactipeptides as previously thought. Indeed, the bioinformatically identified SCIFF peptide "freyrasin" was demonstrated to contain six thioethers linking the ß carbons of six aspartate residues. Another SCIFF, thermocellin, was shown to contain a thioether cross-linked to the γ carbon of threonine. SCIFFs feature a different paradigm of non-α carbon thioether linkages, and they are exclusively formed by radical SAM enzymes, as opposed to the polar chemistry employed during lanthipeptide biosynthesis. Therefore, we propose the renaming of the SCIFF family as radical non-α thioether peptides (ranthipeptides) to better distinguish them from the sactipeptide and lanthipeptide RiPP classes.