Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used hereto provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

Bacteria frequently need to adapt to altered environmental conditions. Adaptation requires changes in gene expression, often mediated by global regulators of transcription. The nucleoid-associated protein H-NS is a key global regulator in Gram-negative bacteria and is believed to be a crucial player in bacterial chromatin organization via its DNA-bridging activity. H-NS activity in vivo is modulated by physico-chemical factors (osmolarity, pH, temperature) and interaction partners. Mechanistically, it is unclear how functional modulation of H-NS by such factors is achieved. Here, we show that a diverse spectrum of H-NS modulators alter the DNA-bridging activity of H-NS. Changes in monovalent and divalent ion concentrations drive an abrupt switch between a bridging and non-bridging DNA-binding mode. Similarly, synergistic and antagonistic co-regulators modulate the DNA-bridging efficiency. Structural studies suggest a conserved mechanism: H-NS switches between a ‘closed’ and an ‘open’, bridging competent, conformation driven by environmental cues and interaction partners.

We have used dissipative particle dynamics (DPD) to simulate surfactant monolayers on the interface between oil and water. With a simple surfactant model, we investigate how variations in size and structure of surfactants influence their ability to reduce the interfacial tension. In particular, we studied the effect of branching of the hydrophobic tail. We found that branched surfactants are more efficient at the interface than linear ones only if the head groups are sufficiently hydrophilic to prevent the molecules from staggering. By combining DPD with a Monte Carlo method, we have imposed constant surfactant chemical potential and (normal) pressure in separate simulations of bulk and interface. From this, we can determine the bulk concentration needed to obtain a given interfacial tension. We found that higher concentrations of branched surfactants are required to obtain the same reduction of the interfacial tension. We argue that the stronger excluded volume interactions which make branched surfactants more efficient at the interface compared to their linear isomers at the same time make them less inclined to adsorb at the interface.

Understanding the dynamics of large-scale conformational changes in proteins still poses a challenge for molecular simulations. We employ transition path sampling of explicit solvent molecular dynamics trajectories to obtain atomistic insight in the reaction network of the millisecond timescale partial unfolding transition in the photocycle of the bacterial sensor photoactive yellow protein. Likelihood maximization analysis predicts the best model for the reaction coordinates of each substep as well as tentative transition states, without further simulation. We find that the unfolding of the α-helical region 43-51 is followed by sequential solvent exposure of both Glu46 and the chromophore. Which of these two residues is exposed first is correlated with the presence of a salt bridge that is part of the N-terminal domain. Additional molecular dynamics simulations indicate that the exposure of the chromophore does not result in a productive pathway. We discuss several possibilities for experimental validation of these predictions. Our results open the way for studying millisecond conformational changes in other medium-sized (signaling) proteins. molecular dynamics | protein unfolding mechanism | transition path sampling | transition state ensemble | committor analysis I nvestigation of the dynamics of large conformational changes in proteins leads to a better understanding of their functioning. Molecular simulation can complement experiments by modeling the dynamical time evolution of biomolecular systems in atomistic detail. Yet, even though straightforward all-atom molecular dynamics (MD) can now access the microsecond regime, simulating the kinetics of folding and other large conformational changes in proteins at ambient conditions has been limited to relatively small (model) proteins (1). Here we show how the transition path sampling technique enables investigation of the atomistic dynamical pathways for the millisecond timescale partial unfolding transition in a biologically relevant signaling protein, the photoactive yellow protein (PYP). Moreover, analysis of the ensembles of unbiased molecular trajectories allows the extraction of the (possibly complex) reaction coordinates.PYP is a water soluble blue-light photo receptor from Halorhodospira halophila (2, 3) composed of 125 amino acids and a covalently bound chromophore (para-coumaric acid, pCA) (4). The protein folds into a Per-Arnt-Sim core capped by an Nterminal domain (5, 6). Upon the absorption of a blue-light photon, PYP undergoes many different conformational rearrangements, culminating in the formation of the signaling state. The first steps in signaling state formation are the isomerization of pCA, on a picosecond timescale (7), followed by a proton transfer, taking place within microseconds (8, 9). These reactions result in an intermediate pB 0 with pCA in a strained configuration and a negative charge on a buried glutamate (Glu46) (10, 11). As a consequence, the protein partially unfolds to expose pCA and Glu46 to bulk water (12-14), completing ...

H-NS proteins act as osmotic sensors translating changes in osmolarity into altered DNA binding properties, thus, regulating enterobacterial genome organization and genes transcription. The molecular mechanism underlying the switching process and its conservation among H-NS family members remains elusive. Here, we focus on the H-NS family protein MvaT from Pseudomonas aeruginosa and demonstrate experimentally that its protomer exists in two different conformations, corresponding to two different functional states. In the half-opened state (dominant at low salt) the protein forms filaments along DNA, in the fully opened state (dominant at high salt) the protein bridges DNA. This switching is a direct effect of ionic strength on electrostatic interactions between the oppositely charged DNA binding and N-terminal domains of MvaT. The asymmetric charge distribution and intramolecular interactions are conserved among the H-NS family of proteins. Therefore, our study establishes a general paradigm for the molecular mechanistic basis of the osmosensitivity of H-NS proteins.

The histone-like nucleoid structuring protein (H-NS) is a nucleoid-associated protein, which is involved in both gene regulation and DNA compaction. H-NS can bind to DNA in two different ways: in trans, by binding to two separate DNA duplexes, or in cis, by binding to different sites on the same duplex. Based on scanning force microscopy imaging and optical trap-driven unzipping assays, it has recently been suggested that DNA compaction may result from the antagonistic effects of H-NS binding to DNA in trans and cis configurations. To get more insight into the compaction mechanism, we constructed a coarse-grained model description of the compaction of bacterial DNA by H-NS. These simulations highlight the fact that DNA compaction indeed results from the subtle equilibrium between several competing factors, which include the deformation dynamics of the plasmid and the several binding modes of protein dimers to DNA, i.e., dangling configurations, cis- and trans-binding. In particular, the degree of compaction is extremely sensitive to the difference in binding energies of the cis and trans configurations. Our simulations also point out that the conformations of the DNA-protein complexes are significantly different in bulk and in planar conditions, suggesting that conformations observed on mica surfaces may differ significantly from those that prevail in living cells.

Using a combination of ultraviolet circular dichroism, temperature-jump transient-infrared spectroscopy, and molecular dynamics simulations, we investigate the effect of salt bridges between different types of charged amino-acid residue pairs on α-helix folding. We determine the stability and the folding and unfolding rates of 12 alanine-based α-helical peptides, each of which has a nearly identical composition containing three pairs of positively and negatively charged residues (either Glu(-)/Arg(+), Asp(-)/Arg(+), or Glu(-)/Lys(+)). Within each set of peptides, the distance and order of the oppositely charged residues in the peptide sequence differ, such that they have different capabilities of forming salt bridges. Our results indicate that stabilizing salt bridges (in which the interacting residues are spaced and ordered such that they favor helix formation) speed up α-helix formation by up to 50% and slow down the unfolding of the α-helix, whereas salt bridges with an unfavorable geometry have the opposite effect. Comparing the peptides with different types of charge pairs, we observe that salt bridges between side chains of Glu(-) and Arg(+) are most favorable for the speed of folding, probably because of the larger conformational space of the salt-bridging Glu(-)/Arg(+) rotamer pairs compared to Asp(-)/Arg(+) and Glu(-)/Lys(+). We speculate that the observed impact of salt bridges on the folding kinetics might explain why some proteins contain salt bridges that do not stabilize the final, folded conformation.

In proteins and enzymes, the local environment of an active cofactor plays an important role in controlling the outcome of a functional reaction. In photoactive yellow protein (PYP), it ensures photoisomerization of the chromophore, a prerequisite for formation of a signaling state. PYP is the prototype of a PAS domain, and the preferred model system for the studies of molecular mechanisms of biological light sensing. We investigated the effect of replacing proline-68, positioned near but not in direct contact with the chromophore, with other neutral amino acids (alanine, glycine, and valine), using ultrafast spectroscopy probing the visible and the mid-IR spectral regions, and molecular simulation to understand the interactions tuning the efficiency of light signaling. Transient absorption measurements indicate that the quantum yield of isomerization in the mutants is lower than the yield observed for the wild type. Subpicosecond mid-IR spectra and molecular dynamics simulations of the four proteins reveal that the hydrogen bond interactions around the chromophore and the access of water molecules in the active site of the protein determine the efficiency of photoisomerization. The mutants provide additional hydrogen bonds to the chromophore, directly and by allowing more water molecules access to its binding pocket. We conclude that proline-68 in the wild type protein optimizes the yield of photochemistry by maintaining a weak hydrogen bond with the chromophore, at the same time restraining the entrance of water molecules close to the alkylic part of pCa. This study provides a molecular basis for the structural optimization of biological light sensing.