Heterogeneous dynamics in partially disordered proteins

MDverse: Shedding Light on the Dark Matter of Molecular Dynamics Simulations

The rise of open science and the absence of a global dedicated data repository for molecular dynamics (MD) simulations has led to the accumulation of MD files in generalist data repositories, constituting the dark matter of MD - data that is technically accessible, but neither indexed, curated, or easily searchable. Leveraging an original search strategy, we found and indexed about 250,000 files and 2,000 datasets from Zenodo, Figshare and Open Science Framework. With a focus on files produced by the Gromacs MD software, we illustrate the potential offered by the mining of publicly available MD data. We identified systems with specific molecular composition and were able to characterize essential parameters of MD simulation such as temperature and simulation length, and could identify model resolution, such as all-atom and coarse-grain. Based on this analysis, we inferred metadata to propose a search engine prototype to explore the MD data. To continue in this direction, we call on the community to pursue the effort of sharing MD data, and to report and standardize metadata to reuse this valuable matter.

Functional in vitro diversity of an intrinsically disordered plant protein during freeze-thawing is encoded by its structural plasticity

Intrinsically disordered late embryogenesis abundant (LEA) proteins play a central role in the tolerance of plants and other organisms to dehydration brought upon, for example, by freezing temperatures, high salt concentration, drought or desiccation, and many LEA proteins have been found to stabilize dehydration-sensitive cellular structures. Their conformational ensembles are highly sensitive to the environment, allowing them to undergo conformational changes and adopt ordered secondary and quaternary structures and to participate in formation of membraneless organelles. In an interdisciplinary approach, we discovered how the functional diversity of the Arabidopsis thaliana LEA protein COR15A found in vitro is encoded in its structural repertoire, with the stabilization of membranes being achieved at the level of secondary structure and the stabilization of enzymes accomplished by the formation of oligomeric complexes. We provide molecular details on intra- and inter-monomeric helix-helix interactions, demonstrate how oligomerization is driven by an α-helical molecular recognition feature (α-MoRF) and provide a rationale that the formation of noncanonical, loosely packed, right-handed coiled-coils might be a recurring theme for homo- and hetero-oligomerization of LEA proteins.

Probing the dynamic landscape of peptides in molecular assemblies by synergized NMR experiments and MD simulations

Peptides or proteins containing small biomolecular aggregates, such as micelles, bicelles, droplets and nanodiscs, are pivotal in many fields ranging from structural biology to pharmaceutics. Monitoring dynamics of such systems has been limited by the lack of experimental methods that could directly detect their fast (picosecond to nanosecond) timescale dynamics. Spin relaxation times from NMR experiments are sensitive to such motions, but their interpretation for biomolecular aggregates is not straightforward. Here we show that the dynamic landscape of peptide-containing molecular assemblies can be determined by a synergistic combination of solution state NMR experiments and molecular dynamics (MD) simulations. Solution state NMR experiments are straightforward to implement without an excessive amount of sample, while direct combination of spin relaxation data to MD simulations enables interpretation of dynamic landscapes of peptides and other aggregated molecules. To demonstrate this, we interpret NMR data from transmembrane, peripheral, and tail anchored peptides embedded in micelles. Our results indicate that peptides and detergent molecules do not rotate together as a rigid body, but peptides rotate in a viscous medium composed of detergent micelle. Spin relaxation times also provide indirect information on peptide conformational ensembles. This work gives new perspectives on peptide dynamics in complex biomolecular assemblies.

Ton motor conformational switch and peptidoglycan role in bacterial nutrient uptake

Active nutrient uptake is fundamental for survival and pathogenicity of Gram-negative bacteria, which operate a multi-protein Ton system to transport essential nutrients like metals and vitamins. This system harnesses the proton motive force at the inner membrane to energize the import through the outer membrane, but the mechanism of energy transfer remains enigmatic. Here, we study the periplasmic domain of ExbD, a crucial component of the proton channel of the Ton system. We show that this domain is a dynamic dimer switching between two conformations representing the proton channel's open and closed states. By in vivo phenotypic assays we demonstrate that this conformational switch is essential for the nutrient uptake by bacteria. The open state of ExbD triggers a disorder to order transition of TonB, enabling TonB to supply energy to the nutrient transporter. We also reveal the anchoring role of the peptidoglycan layer in this mechanism. Herein, we propose a mechanistic model for the Ton system, emphasizing ExbD duality and the pivotal catalytic role of peptidoglycan. Sequence analysis suggests that this mechanism is conserved in other systems energizing gliding motility and membrane integrity. Our study fills important gaps in understanding bacterial motor mechanism and proposes novel antibacterial strategies.

Multiple TonB homologs are important for carbohydrate utilization by Bacteroides thetaiotaomicron

IMPORTANCE

The human gut microbiota, including Bacteroides, is required for the degradation of otherwise undigestible polysaccharides. The gut microbiota uses polysaccharides as an energy source, and fermentation products such as short-chain fatty acids are beneficial to the human host. This use of polysaccharides is dependent on the proper pairing of a TonB protein with polysaccharide-specific TonB-dependent transporters; however, the formation of these protein complexes is poorly understood. In this study, we examine the role of 11 predicted TonB homologs in polysaccharide uptake. We show that two proteins, TonB4 and TonB6, may be functionally redundant. This may allow for the development of drugs targeting Bacteroides species containing only a TonB4 homolog with limited impact on species encoding the redundant TonB6.

Ton Motor Conformational Switch and Peptidoglycan Role in Bacterial Nutrient Uptake

Active nutrient uptake is fundamental for survival and pathogenicity of Gram-negative bacteria, which operate a multi-protein Ton system to transport essential nutrients like metals and vitamins. This system harnesses the proton motive force at the inner membrane to energize the import through the outer membrane, but the mechanism of energy transfer remains enigmatic. Here, we study the periplasmic domain of ExbD, a crucial component of the proton channel of the Ton system. We show that this domain is a dynamic dimer switching between two conformations representing the proton channel's open and closed states. By in vivo phenotypic assays we demonstrate that this conformational switch is essential for the nutrient uptake by bacteria. The open state of ExbD triggers a disorder to order transition of TonB, enabling TonB to supply energy to the nutrient transporter. We also reveal the anchoring role of the peptidoglycan layer in this mechanism. Herein, we propose a mechanistic model for the Ton system, emphasizing ExbD duality and the pivotal catalytic role of peptidoglycan. Sequence analysis suggests that this mechanism is conserved in other systems energizing gliding motility and membrane integrity. Our study fills important gaps in understanding bacterial motor mechanism and proposes novel antibacterial strategies.

Multiple TonB Homologs are Important for Carbohydrate Utilization by Bacteroides thetaiotaomicron

The human gut microbiota is able to degrade otherwise undigestible polysaccharides, largely through the activity of the Bacteroides. Uptake of polysaccharides into Bacteroides is controlled by TonB-dependent transporters (TBDT) whose transport is energized by an inner membrane complex composed of the proteins TonB, ExbB, and ExbD. Bacteroides thetaiotaomicron (B. theta) encodes 11 TonB homologs which are predicted to be able to contact TBDTs to facilitate transport. However, it is not clear which TonBs are important for polysaccharide uptake. Using strains in which each of the 11 predicted tonB genes are deleted, we show that TonB4 (BT2059) is important but not essential for proper growth on starch. In the absence of TonB4, we observed an increase in abundance of TonB6 (BT2762) in the membrane of B. theta, suggesting functional redundancy of these TonB proteins. Growth of the single deletion strains on pectin galactan, chondroitin sulfate, arabinan, and levan suggests a similar functional redundancy of the TonB proteins. A search for highly homologous proteins across other Bacteroides species and recent work in B. fragilis suggests that TonB4 is widely conserved and may play a common role in polysaccharide uptake. However, proteins similar to TonB6 are found only in B. theta and closely related species suggesting that the functional redundancy of TonB4 and TonB6 may be limited across the Bacteroides. This study extends our understanding of the protein network required for polysaccharide utilization in B. theta and highlights differences in TonB complexes across Bacteroides species.

Accurate Simulations of Lipid Monolayers Require a Water Model with Correct Surface Tension

Lipid monolayers provide our lungs and eyes their functionality and serve as proxy systems in biomembrane research. Therefore, lipid monolayers have been studied intensively including using molecular dynamics simulations, which are able to probe their lateral structure and interactions with, e.g., pharmaceuticals or nanoparticles. However, such simulations have struggled in describing the forces at the air-water interface. Particularly, the surface tension of water and long-range van der Waals interactions have been considered critical, but their importance in monolayer simulations has been evaluated only separately. Here, we combine the recent C36/LJ-PME lipid force field that includes long-range van der Waals forces with water models that reproduce experimental surface tensions to elucidate the importance of these contributions in monolayer simulations. Our results suggest that a water model with correct surface tension is necessary to reproduce experimental surface pressure-area isotherms and monolayer phase behavior. The latter includes the liquid expanded and liquid condensed phases, their coexistence, and the opening of pores at the correct area per lipid upon expansion. Despite these improvements of the C36/LJ-PME with certain water models, the standard cutoff-based CHARMM36 lipid model with the 4-point OPC water model still provides the best agreement with experiments. Our results emphasize the importance of using high-quality water models in applications and parameter development in molecular dynamics simulations of biomolecules.

Energetics and J-coupling constants for Ala, Gly, and Val peptides demonstrated using ABEEM polarizable force field in vacuo and an aqueous solution

The development of an atom-bond electronegativity equalisation method at the σπ-level (ABEEM) polarisable force field (PFF) for peptides is presented. ABEEM PFF utilises a fluctuating charge model to explicitly describe...

The diversity of molecular interactions involving intrinsically disordered proteins: A molecular modeling perspective

Intrinsically Disordered Proteins and Regions (IDPs/IDRs) are key components of a multitude of biological processes. Conformational malleability enables IDPs/IDRs to perform very specialized functions that cannot be accomplished by globular proteins. The functional role for most of these proteins is related to the recognition of other biomolecules to regulate biological processes or as a part of signaling pathways. Depending on the extent of disorder, the number of interacting sites and the type of partner, very different architectures for the resulting assemblies are possible. More recently, molecular condensates with liquid-like properties composed of multiple copies of IDPs and nucleic acids have been proven to regulate key processes in eukaryotic cells. The structural and kinetic details of disordered biomolecular complexes are difficult to unveil experimentally due to their inherent conformational heterogeneity. Computational approaches, alone or in combination with experimental data, have emerged as unavoidable tools to understand the functional mechanisms of this elusive type of assemblies. The level of description used, all-atom or coarse-grained, strongly depends on the size of the molecular systems and on the timescale of the investigated mechanism. In this mini-review, we describe the most relevant architectures found for molecular interactions involving IDPs/IDRs and the computational strategies applied for their investigation.

How to assess the structural dynamics of transcription factors by integrating sparse NMR and EPR constraints with molecular dynamics simulations

We review recent advances in modeling structural ensembles of transcription factors from nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) spectroscopic data, integrated with molecular dynamics (MD) simulations. We focus on approaches that confirm computed conformational ensembles by sparse constraints obtained from magnetic resonance. This combination enables the deduction of functional and structural protein models even if nuclear Overhauser effects (NOEs) are too scarce for conventional structure determination. We highlight recent insights into the folding-upon-DNA binding transitions of intrinsically disordered transcription factors that could be assessed using such integrative approaches.

Ferreira T, Ollila OHS, Miettinen MS. Using Open Data to Rapidly Benchmark Biomolecular Simulations: Phospholipid Conformational Dynamics

Molecular dynamics (MD) simulations are widely used to monitor time-resolved motions of biomacromolecules, although it often remains unknown how closely the conformational dynamics correspond to those occurring in real life. Here, we used a large set of open-access MD trajectories of phosphatidylcholine (PC) lipid bilayers to benchmark the conformational dynamics in several contemporary MD models (force fields) against nuclear magnetic resonance (NMR) data available in the literature: effective correlation times and spin-lattice relaxation rates. We found none of the tested MD models to fully reproduce the conformational dynamics. That said, the dynamics in CHARMM36 and Slipids are more realistic than in the Amber Lipid14, OPLS-based MacRog, and GROMOS-based Berger force fields, whose sampling of the glycerol backbone conformations is too slow. The performance of CHARMM36 persists when cholesterol is added to the bilayer, and when the hydration level is reduced. However, for conformational dynamics of the PC headgroup, both with and without cholesterol, Slipids provides the most realistic description because CHARMM36 overestimates the relative weight of ∼1 ns processes in the headgroup dynamics. We stress that not a single new simulation was run for the present work. This demonstrates the worth of open-access MD trajectory databanks for the indispensable step of any serious MD study: benchmarking the available force fields. We believe this proof of principle will inspire other novel applications of MD trajectory databanks and thus aid in developing biomolecular MD simulations into a true computational microscope-not only for lipid membranes but for all biomacromolecular systems.