Search

All-atom simulations have become increasingly popular to study conformational and dynamical properties of nucleic acids as they are accurate and provide high spatial and time resolutions. This high resolution however comes at a heavy computational cost and within the time scales of simulations nucleic acids weakly fluctuate around their ideal structure exploring a limited set of conformations. We introduce the RBB-NA algorithm which is capable of controlling rigid base parameters in all-atom simulations of Nucleic Acids. With suitable biasing potentials this algorithm can "force" a DNA or RNA molecule to assume specific values of the six rotational (tilt, roll, twist, buckle, propeller, opening) and/or the six translational parameters (shift, slide, rise, shear, stretch, stagger). The algorithm enables the use of advanced sampling techniques to probe the structure and dynamics of locally strongly deformed Nucleic Acids. We illustrate its performance showing some examples in which DNA is strongly twisted, bent or locally buckled. In these examples RBB-NA reproduces well the unconstrained simulations data and other known features of DNA mechanics, but it also allows one to explore the anharmonic behavior characterizing the mechanics of nucleic acids in the high deformation regime., 12 pages, 6 figures

Flexibility is essential for many proteins to function, but can be difficult to characterize. Experiments lack resolution in space and time, while the time scales involved are prohibitively long for straightforward molecular dynamics simulations. In this work, we present a multiple state transition path sampling simulation study of a protein that has been notoriously difficult to characterize in its active state. The GTPase enzyme KRas is a signal transduction protein in pathways for cell differentiation, growth, and division. When active, KRas tightly binds guanosine triphosphate (GTP) in a rigid state. The protein-GTP complex can also visit more flexible states, in which it is not active. KRas mutations can affect the conversion between these rigid and flexible states, thus prolonging the activation of signal transduction pathways, which may result in tumor formation. In this work, we apply path sampling simulations to investigate the dynamic behavior of KRas-4B (wild type, WT) and the oncogenic mutant Q61L (Q61L). Our results show that KRas visits several conformational states, which are the same for WT and Q61L. The multiple state transition path sampling (MSTPS) method samples transitions between the different states in a single calculation. Tracking which transitions occur shows large differences between WT and Q61L. The MSTPS results further reveal that for Q61L, a route to a more flexible state is inaccessible, thus shifting the equilibrium to more rigid states. The methodology presented here enables a detailed characterization of protein flexibility on time scales not accessible with brute-force molecular dynamics simulations.

DNA predominantly contains Watson–Crick (WC) base pairs, but a non-negligible fraction of base pairs are in the Hoogsteen (HG) hydrogen bonding motif at any time. In HG, the purine is rotated ∼180° relative to the WC motif. The transitions between WC and HG may play a role in recognition and replication, but are difficult to investigate experimentally because they occur quickly, but only rarely. To gain insight into the mechanisms for this process, we performed transition path sampling simulations on a model nucleotide sequence in which an AT pair changes from WC to HG. This transition can occur in two ways, both starting with loss of hydrogen bonds in the base pair, followed by rotation around the glycosidic bond. In one route the adenine base converts from WC to HG geometry while remaining entirely within the double helix. The other route involves the adenine leaving the confines of the double helix and interacting with water. Our results indicate that this outside route is more probable. We used transition interface sampling to compute rate constants and relative free energies for the transitions between WC and HG. Our results agree with experiments, and provide highly detailed insights into the mechanisms of this important process.

Computational chemistry uses computer simulation to assist in solving chemical problems. Typical workflows of computational chemists include the use of dozens of utilities. 3D modeling programs are powerful tools that help researchers visualize their work and create illustrative graphics. In this review, we describe and highlight tools and visualization packages that are commonly used in the field of (bio-)chemistry and material science.

Bacillus subtilis forms dormant spores upon nutrient depletion. Germinant receptors (GRs) in spore’s inner membrane respond to ligands such as L-alanine, and trigger spore germination. In B. subtilis spores, GerA is the major GR, and has three subunits, GerAA, GerAB, and GerAC. L-Alanine activation of GerA requires all three subunits, but which binds L-alanine is unknown. To date, how GRs trigger germination is unknown, in particular due to lack of detailed structural information about B subunits. Using homology modelling with molecular dynamics (MD) simulations, we present structural predictions for the integral membrane protein GerAB. These predictions indicate that GerAB is an α-helical transmembrane protein containing a water channel. The MD simulations with free L-alanine show that alanine binds transiently to specific sites on GerAB. These results provide a starting point for unraveling the mechanism of L-alanine mediated signaling by GerAB, which may facilitate early events in spore germination.

The mechanical properties of DNA are typically described by elastic theories with purely local couplings (on-site models). We discuss and analyze coarse-grained (oxDNA) and all-atom simulations, which indicate that in DNA distal sites are coupled. Hence, off-site models provide a more realistic description of the mechanics of the double helix. We show that off-site interactions are responsible for a length scale dependence of the elasticity, and we develop an analytical framework to estimate bending and torsional persistence lengths in models including these interactions. Our simulations indicate that off-site couplings are particularly strong for certain degrees of freedom, while they are very weak for others. If stiffness parameters obtained from DNA data are used, the theory predicts large length scale dependent effects for torsional fluctuations and a modest effect in bending fluctuations, which is in agreement with experiments., 16 pages, 9 figures; video abstract at https://youtu.be/8FUdFw8IbuY

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 P. 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 strengths on electrostatic interactions between the appositively 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.

In the past decade, great progress has been made in the development of enhanced sampling methods, aimed at overcoming the time-scale limitations of molecular dynamics (MD) simulations. Many sampling schemes rely on adding an external bias to favor the sampling of transitions and to estimate the underlying free energy landscape. Nevertheless, sampling molecular processes described by many order parameters, or collective variables (CVs), such as complex biomolecular transitions, remains often very challenging. The computational cost has a prohibitive scaling with the dimensionality of the CV-space. Inspiration can be taken from methods that focus on localizing transition pathways: the CV-space can be projected onto a path-CV that connects two stable states, and a bias can be exerted onto a one-dimensional parameter that captures the progress of the transition along the path-CV. In principle, such a sampling scheme can handle an arbitrarily large number of CVs. A standard enhanced sampling technique combined with an adaptive path-CV can then locate the mean transition pathway and obtain the free energy profile along the path. In this chapter, we discuss the adaptive path-CV formalism and its numerical implementation. We apply the path-CV with several enhanced sampling methods-steered MD, metadynamics, and umbrella sampling-to a biologically relevant process: the Watson-Crick to Hoogsteen base-pairing transition in double-stranded DNA. A practical guide is provided on how to recognize and circumvent possible pitfalls during the calculation of a free energy landscape that contains multiple pathways. Examples are presented on how to perform enhanced sampling simulations using PLUMED, a versatile plugin that can work with many popular MD engines.

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 macromo-lecular 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 here to 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 contain several nucleoid-associated proteins that organize their genomic DNA into the nucleoid by bending, wrapping or bridging DNA. The Histone-like Nucleoid Structuring protein H-NS found in many Gram-negative bacteria is a DNA bridging protein and can structure DNA by binding to two separate DNA duplexes or to adjacent sites on the same duplex, depending on external conditions. Several nucleotide sequences have been identified to which H-NS binds with high affinity, indicating H-NS prefers AT-rich DNA. To date, highly detailed structural information of the H-NS DNA complex remains elusive. Molecular simulation can complement experiments by modelling structures and their time evolution in atomistic detail. In this paper we report an exploration of the different binding modes of H-NS to a high affinity nucleotide sequence and an estimate of the associated rate constant. By means of molecular dynamics simulations, we identified three types of binding for H-NS to AT-rich DNA. To further sample the transitions between these binding modes, we performed Replica Exchange Transition Interface Sampling, providing predictions of the mechanism and rate constant of H-NS binding to DNA. H-NS interacts with the DNA through a conserved QGR motif, aided by a conserved arginine at position 93. The QGR motif interacts first with phosphate groups, followed by the formation of hydrogen bonds between acceptors in the DNA minor groove and the sidechains of either Q112 or R114. After R114 inserts into the minor groove, the rest of the QGR motif follows. Full insertion of the QGR motif in the minor groove is stable over several tens of nanoseconds, and involves hydrogen bonds between the bases and both backbone and sidechains of the QGR motif. The rate constant for the process of H-NS binding to AT-rich DNA resulting in full insertion of the QGR motif is in the order of 106 M−1s−1, which is rate limiting compared to the non-specific association of H-NS to the DNA backbone at a rate of 108 M−1s−1., Author summary The Histone-like Nucleoid Structuring protein (H-NS) occurs in enterobacteria, such as Salmonella typhimurium and Escherichia coli, and structures DNA by forming filaments along DNA duplexes. Several nucleotide sequences have been identified to which H-NS binds with high affinity. Yet, obtaining highly detailed structural information of the H-NS DNA complex has proven to be a major challenge, which has not been yet resolved. By employing molecular dynamics simulations we were able to provide high resolution insights into the mechanism of DNA binding by H-NS. We identified various ways in which H-NS can bind to DNA. In all binding events, a conserved region in the protein initiates the association of H-NS to DNA. Our results show that H-NS binds in the minor groove of AT-rich DNA via a series of intermediate steps. Using advanced molecular simulation methods we predicted that the process of H-NS binding to the DNA backbone to full insertion into the minor groove occurs in the order of a million times per second, which is slower than the non-specific association of H-NS to the DNA backbone.

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 the past decade, great progress has been made in the development of enhanced sampling methods, aimed at overcoming the time-scale limitations of molecular dynamics (MD) simulations. Many sampling schemes rely on adding an external bias to favor the sampling of transitions and to estimate the underlying free energy landscape. Nevertheless, sampling molecular processes described by many order parameters, or collective variables (CVs), such as complex biomolecular transitions, remains often very challenging. The computational cost has a prohibitive scaling with the dimensionality of the CV-space. Inspiration can be taken from methods that focus on localizing transition pathways: the CV-space can be projected onto a path-CV that connects two stable states, and a bias can be exerted onto a one-dimensional parameter that captures the progress of the transition along the path-CV. In principle, such a sampling scheme can handle an arbitrarily large number of CVs. A standard enhanced sampling technique combined with an adaptive path-CV can then locate the mean transition pathway and obtain the free energy profile along the path. In this chapter, we discuss the adaptive path-CV formalism and its numerical implementation. We apply the path-CV with several enhanced sampling methods-steered MD, metadynamics, and umbrella sampling-to a biologically relevant process: the Watson-Crick to Hoogsteen base-pairing transition in double-stranded DNA. A practical guide is provided on how to recognize and circumvent possible pitfalls during the calculation of a free energy landscape that contains multiple pathways. Examples are presented on how to perform enhanced sampling simulations using PLUMED, a versatile plugin that can work with many popular MD engines.

Azobenzene-based molecular photoswitches are becoming increasingly important for the development of photoresponsive, functional soft-matter material systems. Upon illumination with light, fast interconversion between a more stable trans and a metastable cis configuration can be established resulting in pronounced changes in conformation, dipole moment or hydrophobicity. A rational design of functional photosensitive molecules with embedded azo moieties requires a thorough understanding of isomerization mechanisms and rates, especially the thermally activated relaxation. For small azo derivatives considered in the gas phase or simple solvents, Eyring’s classical transition state theory (TST) approach yields useful predictions for trends in activation energies or corresponding half-life times of the cis isomer. However, TST or improved theories cannot easily be applied when the azo moiety is part of a larger molecular complex or embedded into a heterogeneous environment, where a multitude of possible reaction pathways may exist. In these cases, only the sampling of an ensemble of dynamic reactive trajectories (transition path sampling, TPS) with explicit models of the environment may reveal the nature of the processes involved. In the present work we show how a TPS approach can conveniently be implemented for the phenomenon of relaxation–isomerization of azobenzenes starting with the simple examples of pure azobenzene and a push–pull derivative immersed in a polar (DMSO) and apolar (toluene) solvent. The latter are represented explicitly at a molecular mechanical (MM) and the azo moiety at a quantum mechanical (QM) level. We demonstrate for the push–pull azobenzene that path sampling in combination with the chosen QM/MM scheme produces the expected change in isomerization pathway from inversion to rotation in going from a low to a high permittivity (explicit) solvent model. We discuss the potential of the simulation procedure presented for comparative calculation of reaction rates and an improved understanding of activated states.

Biological signal transduction by photoactive yellow protein (PYP) in halophilic purple sulfur bacteria is initiated by trans-to-cis isomerization of the p-coumaric acid chromophore (pCa) of PYP. pCa is engaged in two short hydrogen bonds with protein residues E46 and Y42, and it is negatively charged at the phenolate oxygen. We investigated the role in the isomerization process of the E46 short hydrogen bond and that of the negative charge on the anionic phenolate moiety of the chromophore. We used wild-type PYP and the mutant E46A, in protonated and deprotonated states (referred to as pE46A and dpE46A, respectively), to reduce the number of hydrogen bond interactions between the pCa phenolate oxygen and the protein and to vary the negative charge density in the chromophore-binding pocket. Their effects on the yield and rate of chromophore isomerization were determined by ultrafast spectroscopy. Molecular dynamics simulations were used to relate these results to structural changes in the mutant protein. We found that deprotonated pCa in E46A has a slower isomerization rate as the main part of this reaction was associated with time constants of 1 and 6 ps, significantly slower than the 0.6 ps time constant in wild-type PYP. The quantum yield of isomerization in dpE46A was estimated to be 30 ± 2%, and that of pE46A was 32 ± 3%, very close to the value determined for wtPYP of 32 ± 2%. Relaxation of the isomerized product state I0 to I1 was faster in dpE46A. We conclude that the negative charge on pCa stabilized by the short hydrogen bonds with E46 and Y42 affects the rate of isomerization but not the quantum yield of isomerization. With this information, we propose a scheme for the potential energy surfaces involved in the isomerization and suggest protein motions near the pCa backbone as key events in successful isomerization.

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., eLife digest The genetic information every cell needs to work properly is encoded in molecules of DNA that are much longer than the cell itself. A key challenge in biology is to understand how DNA is organized to fit inside each cell, whilst still providing access to the information that it contains. Since the way DNA is organized can influence which genes are active, rearranging DNA plays an important role in controlling how cells behave. In Escherichia coli and many other bacteria, a protein called H-NS contributes to DNA reorganization by forming or rupturing loops in the DNA in response to changes in temperature, the levels of salt and other aspects of the cell’s surroundings. In controlling loop formation, it dictates whether specific genes are switched on or off. However, it remains unclear how H-NS detects the environmental changes. To address this question, van der Valk et al. used biochemical techniques to study the activity of H-NS from E. coli under different environmental conditions. The experiments show that changes in the environment cause structural changes to H-NS, altering its ability to form DNA loops. A previously unnoticed region of the protein acts as a switch to control these structural changes, and ultimately affects which genes are active in the cell. These findings shed new light on how bacteria organize their DNA and the strategies they have developed to adapt to different environments. The new protein region identified in H-NS may also be present in similar proteins found in other organisms. In the future, this knowledge may ultimately help to develop new antibiotic drugs that target H-NS proteins in bacteria.

Salt bridges are known to play an essential role in the thermodynamic stability of the folded conformation of many proteins, but their influence on the kinetics of folding remains largely unknown. Here, we investigate the effect of Glu-Arg salt bridges on the kinetics of α-helix folding using temperature-jump transient-infrared spectroscopy and steady-state UV circular dichroism. We find that geometrically optimized salt bridges (Glu– and Arg+ are spaced four peptide units apart, and the Glu/Arg order is such that the side-chain rotameric preferences favor salt-bridge formation) significantly speed up folding and slow down unfolding, whereas salt bridges with unfavorable geometry slow down folding and slightly speed up unfolding. Our observations suggest a possible explanation for the surprising fact that many biologically active proteins contain salt bridges that do not stabilize the native conformation: these salt bridges might have a kinetic rather than a thermodynamic function.

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 inGram-negative bacteria, and is believed to be a crucial player in bacterial chromatin organization via its DNA bridging activity. H-NS activityin vivois 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 ability of H-NS to bridge DNA. 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 investigate the molecular mechanism of monomer addition to a growing amyloid fibril composed of the main amyloidogenic region from the insulin peptide hormone, the LVEALYL heptapeptide. Applying transition path sampling in combination with reaction coordinate analysis reveals that the transition from a docked peptide to a locked, fully incorporated peptide can occur in two ways. Both routes involve the formation of backbone hydrogen bonds between the three central amino acids of the attaching peptide and the fibril, as well as a reorientation of the central Glu side chain of the locking peptide toward the interface between two β-sheets forming the fibril. The mechanisms differ in the sequence of events. We also conclude that proper docking is important for correct alignment of the peptide with the fibril, as alternative pathways result in misfolding.

The histone-like nucleoid structuring protein (H-NS) is a DNA-organizing protein in bacteria. It contains a DNA-binding domain and a dimerization domain, connected by a flexible linker region. Dimerization occurs through the formation of a helical bundle, including a coiled-coil interaction motif. Two conformations have been resolved, for different sequences of Escherichia coli H-NS, resulting in an antiparallel coiled-coil for the shorter wild-type sequence, and a parallel coiled-coil for the longer C21S mutant. Because H-NS functions as a thermo- and osmosensor, these conformations may both be functionally relevant. Molecular simulation can complement experiments by modeling the dynamical time evolution of biomolecular systems in atomistic detail. We performed a molecular-dynamics study of the H-NS dimerization domain, showing that the parallel complex is sensitive to changes in salt conditions: it is unstable in absence of NaCl, but stable at physiological salt concentrations. In contrast, the stability of the antiparallel complex is not salt-dependent. The stability of the parallel complex also appears to be affected by mutation of the critical but nonconserved cysteine residue at position 21, whereas the antiparallel complex is not. Together, our simulations suggest that osmoregulation could be mediated by changes in the ratio of parallel- and antiparallel-oriented H-NS dimers.

Conformational changes in proteins often take place on long time scales compared to the molecular timescale. These long time scales, related to high free energy barriers, make such processes difficult to access with a straightforward molecular dynamics approach. The transition path sampling technique has been developed to overcome such timescale differences without assuming a predefined reaction coordinate. We review the transition path sampling methodology with the application of protein conformational change in mind. Using three case studies, based on previous work, we elucidate the strengths and pitfalls of the method. First, the extensive work on the folding of Trp-cage reveals how to sample parallel pathways, how to obtain rate constants, and how to extract reaction coordinates. The second case-study, on the folding of Trpzip4 beta-hairpin, illustrates how to treat long-lived intermediates in a conformational change. The final example showcases the light-triggered conformational transition of Photo-active Yellow Protein into its signaling state, highlighting the wealth of insight that can be gathered from a transition path sampling approach, including new hypotheses for reaction mechanisms. We end with an outlook discussing future developments and application of the methodology. (C) 2011 Elsevier B. V. All rights reserved.

Chemokines (chemotactic cytokines) can have direct antimicrobial activity, which is apparently related to the presence of a distinct positively charged patch on the surface. However, chemokines can retain antimicrobial activity upon linearization despite the loss of their positive patch, thus questioning the importance of this patch for activity. Thrombocidin-1 (TC-1) is a microbicidal protein isolated from human blood platelets. TC-1 only differs from the chemokine NAP-2/CXCL7 by a two-amino acid C-terminal deletion, but this truncation is crucial for antimicrobial activity. We assessed the structure-activity relationship for antimicrobial activity of TC-1. Reduction of the charge of the TC-1-positive patch by replacing lysine 17 with alanine reduced the activity against bacteria and almost abolished activity against the yeast Candida albicans. Conversely, augmentation of the positive patch by increasing charge density or size resulted in a 2-3-fold increased activity against Staphylococcus aureus, Escherichia coli, and Bacillus subtilis but did not substantially affect activity against C. albicans. Reduction of TC-1 resulted in loss of the folded conformation, but this disruption of the positive patch did not affect antimicrobial activity. Using overlapping 15-mer synthetic peptides, we demonstrate peptides corresponding to the N-terminal part of TC-1 to have similar antimicrobial activity as intact TC-1. Although we demonstrate that the positive patch is essential for activity of folded TC-1, unfolded TC-1 retained antimicrobial activity despite the absence of a positive patch. This activity is probably exerted by a linear peptide stretch in the N-terminal part of the molecule. We conclude that intact TC-1 and unfolded TC-1 exert antimicrobial activity via distinct structural elements.

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. © 2011 American Chemical Society.

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.

Chitotriosidase (CHIT1) is a chitinase that is secreted by activated macrophages. Plasma chitotriosidase activity reflects the presence of lipid-laden macrophages in patients with Gaucher disease. CHIT1 activity can be conveniently measured using fluorogenic 4-methylumbelliferyl (4MU)-chitotrioside or 4MU-chitobioside as the substrate, however, nonsaturating concentrations have to be used because of apparent substrate inhibition. Saturating substrate concentrations can, however, be used with the newly designed substrate 4MU-deoxychitobioside. We studied the impact of a known polymorphism, G102S, on the catalytic properties of CHIT1. The G102S allele was found to be common in type I Gaucher disease patients in the Netherlands ( approximately 24% of alleles). The catalytic efficiency of recombinant Ser102 CHIT1 was approximately 70% that of wild-type Gly102 CHIT1 when measured with 4MU-chitotrioside at a nonsaturating concentration. However, the activity was normal with 4MU-deoxychitobioside as the substrate at saturating concentrations, consistent with predictions from molecular dynamics simulations. In conclusion, interpretation of CHIT1 activity measurements with 4MU-chitotrioside with respect to CHIT1 protein concentrations depends on the presence of Ser102 CHIT1 in an individual, complicating estimation of the body burden of storage macrophages. Use of the superior 4MU-deoxychitobioside substrate avoids such complications because activity towards this substrate under saturating conditions is not affected by the G102S substitution.

The general stress response of Bacillus subtilis can be activated by stimuli such as the addition of salt or ethanol and with blue light. In the latter response, YtvA activates sigma(B) through a cascade of Rsb proteins, organized in stressosomes. YtvA is composed of an N-terminal LOV (light, oxygen, and voltage) domain and a C-terminal STAS (sulfate transporter and anti-sigma factor) domain and shows light-modulated GTP binding in vitro. Here, we examine the mechanism of YtvA-mediated activation of sigma(B) in vivo with site-directed mutagenesis. Constitutive off and constitutive on mutations have been identified. Disruption of GTP binding in the STAS domain eliminates light activation of sigma(B). In contrast, modification of sites relevant for phosphorylation of STAS domains does not affect the stress response significantly. The data obtained are integrated into a model for the structure of full-length YtvA, which presumably functions as a dimer.

The bacterial sensor Photoactive Yellow Protein (PYP) signals the presence of blue light by undergoing a series of conformational changes. We present atomistic Parallel Tempering (Replica Exchange Molecular Dynamics) simulations of conformational changes occurring during the photo-cycle of PYP. First, we study the signaling state formation of PYP in detail. Our previous simulations have shown that the formation of the signaling state is characterized by the solvent exposure of both the chromophore and Glu46 (Vreede J, Crielaard W, Hellingwerf KJ, Bolhuis PG. Biophys J, 2005;8:3525-3535). Subsequent NMR results agreed with this prediction, but as these experiments were performed on an N-terminally truncated mutant, a simulation of this mutant would further substantiate our previous results. Here, we compare simulations of the truncated PYP to the NMR structures, as well as to the wild type predictions. This comparison also gives some insight into the role of the N-terminal domain of PYP, which restricts the movement of the chromophore binding pocket (CBP) in the wild type. Second, we report simulations of the recovery of the receptor state from the signaling state. While we did not observe complete refolding of the protein, we did observe transient interactions between residues of the CBP occurring when the chromophore is in a trans configuration. Using simulations that sample anomalous exposure of the chromophore in the receptor state, we were able to sample chromophore re-entry into its binding pocket. While the involved time scales prohibit drawing definitive conclusions even when using parallel tempering, we nevertheless propose that the formation of a helix in the CBP is essential for a successful recovery of the receptor state, and forms a kinetic barrier in this process.

The bacterium Agrobacterium tumefaciens mediates quorum sensing via TraR, a N-3-oxo-octanoyl-l-homoserine lactone (OOHL) signaling transcription factor. TraR consists of an N-terminal Per-Arnt-Sim domain (PAS-domain) containing OOHL and a C-terminal helix-turn-helix motif (HTF-motif) binding DNA, thus comprising a complete signaling pathway. Using Molecular Dynamics, we investigate the influence of OOHL on the dynamics of the protein. The OOHL binding pocket is in contact with water, leading to dynamic changes in OOHL-protein interactions. These rearrangements increase fluctuations in the PAS-domain, and, interestingly, also enhance fluctuations in the HTF-motif. Our results imply OOHL in inducing fluctuations in TraR that may facilitate DNA binding.

The effective volume occupied by the genomes of all forms of life far exceeds that of the cells in which they are contained. Therefore, all organisms have developed mechanisms for compactly folding and functionally organizing their genetic material. Through recent advances in fluorescent microscopy and 3C-based technologies, we finally have a first glimpse into the complex mechanisms governing the 3-D folding of genomes. A key feature of genome organization in all domains of life is the formation of DNA loops. Here, we describe the main players in DNA organization with a focus on DNA-bridging proteins. Specifically, we discuss the properties of the bacterial DNA-bridging protein H-NS. Via two different modes of binding to DNA, this protein is a key driver of bacterial genome organization and provides a link between 3-D organization and transcription regulation. Importantly, H-NS function is modulated in response to environmental cues, which are translated into adapted gene expression patterns. We delve into the mechanisms underlying DNA looping and explore the complex and subtle modulation of these diverse, yet difficult-to-study, structures. DNA looping is universal and a conserved mechanism of genome organization throughout all domains of life.

As a bacterial blue light sensor the photoactive yellow protein (PYP) undergoes conformational changes upon signal transduction. The absorption of a photon triggers a series of events that are initially localized around the protein chromophore, extends to encompass the whole protein within microseconds, and leads to the formation of the transient pB signaling state. We study the formation of this signaling state pB by molecular simulation and predict its solution structure. Conventional straightforward molecular dynamics is not able to address this formation process due to the long (microsecond) timescales involved, which are (partially) caused by the presence of free energy barriers between the metastable states. To overcome these barriers, we employed the parallel tempering (or replica exchange) method, thus enabling us to predict qualitatively the formation of the PYP signaling state pB. In contrast to the receptor state pG of PYP, the characteristics of this predicted pB structure include a wide open chromophore-binding pocket, with the chromophore and Glu46 fully solvent-exposed. In addition, loss of α-helical structure occurs, caused by the opening motion of the chromophore-binding pocket and the disruptive interaction of the negatively charged Glu46 with the backbone atoms in the hydrophobic core of the N-terminal cap. Recent NMR experiments agree very well with these predictions.

Biological photoreceptors are very suitable for studies of the structure-function relationship in proteins. Of the many possible candidates, from at least six different families, we discuss here photoactive yellow protein (PYP, a member of the Xanthopsin family). The excellent physical and photochemical stability of PYP has allowed recording of a large amount of data relevant for this characteristic. The timescäles relevant for functional transitions in PYP are most easily studied with transient spectroscopy. Application of the combination of UV/Vis and FTIR spectroscopy has revealed (besides multiple electronically-excited states) the subsequent involvement of the following transient intermediates: I0, pR1, pR 2, pB′, pB and pBdeprot. The structures of these transient intermediates have been characterized with general structural techniques like X-ray diffraction and NMR spectroscopy. Most extensive results have been obtained with X-ray diffraction, applied both in a time-resolved mode and by making use of low-temperature trapping. This has led to a detailed description of the reversible spatial changes in PYP that follow light activation. The structure of the transient intermediates, however, is exquisitely sensitive to constraints from the molecular environment of the protein, like the absence/presence of a crystalline lattice. Extrapolating these results to the signalling state of PYP in vivo is extremely complicated but also very challenging. Rather than measure, one can also simulate and/or calculate the spatial structures of transient intermediates of PYP. This challenge is most straightforwardly tackled with a detailed analysis of the intrinsic dynamics of the stable receptor state of PYP. Such studies have revealed striking similarities in the dynamics of all PAS domain proteins of which a spatial structure has been determined. The next challenge is to properly describe the change in configuration of the chromophore, induced by light, and to expand such simulations to the timescale relevant for completion of this photocycle.

PAS (PER-ARNT-SIM) domains are a family of sensor protein domains involved in signal transduction in a wide range of organisms. Recent structural studies have revealed that these domains contain a structurally conserved α/β-fold, whereas almost no conservation is observed at the amino acid sequence level. The photoactive yellow protein, a bacterial light sensor, has been proposed as the PAS structural prototype yet contains an N-terminal helix-turn-helix motif not found in other PAS domains. Here we describe the atomic resolution structure of a photoactive yellow protein deletion mutant lacking this motif, revealing that the PAS domain is indeed able to fold independently and is not affected by the removal of these residues. Computer simulations of currently known PAS domain structures reveal that these domains are not only structurally conserved but are also similar in their conformational flexibilities. The observed motions point to a possible common mechanism for communicating ligand binding/activation to downstream transducer proteins.

How chaperones interact with protein chains to assist in their folding is a central open question in biology. Obtaining atomistic insight is challenging in particular, given the transient nature of the chaperone-substrate complexes and the large system sizes. Recent single-molecule experiments have shown that the chaperone Trigger Factor (TF) not only binds unfolded protein chains, but can also guide protein chains to their native state by interacting with partially folded structures. Here, we used all-atom MD simulations to provide atomistic insights into how Trigger Factor achieves this chaperone function. Our results indicate a crucial role for the tips of the finger-like appendages of TF in the early interactions with both unfolded chains and partially folded structures. Unfolded chains are kinetically trapped when bound to TF, which suppresses the formation of transient, non-native end-to-end contacts. Mechanical flexibility allows TF to hold partially folded structures with two tips (in a pinching configuration), and to stabilize them by wrapping around its appendages. This encapsulation mechanism is distinct from that of chaperones such as GroEL, and allows folded structures of diverse size and composition to be protected from aggregation and misfolding interactions. The results suggest that an ATP cycle is not required to enable both encapsulation and liberation., Author Summary Trigger Factor (TF) is an ATP-independent chaperone protein that assists in folding and prevents misfolding. Up to now, it is a general unsolved question how chaperones assist in the folding of protein chains. Experimental methods that can probe at the length and timescales of inter-residue interactions are scarce, while the systems are too large—and the folding process too long—to be studied by computer simulations. To overcome these obstacles, the authors performed molecular dynamics simulations at key moments along the folding pathway, and address the changes in the folding and unfolding dynamics of protein chains while in contact with TF. This study provides the first detailed view on a chaperone-protein complex in different stages of folding and offers an explanation for the ability of TF to guide chains to their native state. Moreover, the results demonstrates the role of TF’s flexibility in interacting with a wide range of client states. Overall, it explains how TF can interact with many types of substrates in various stages of folding, without the need for an ATP cycle to switch between encapsulation and liberation of client proteins.

DNA duplexes predominantly contain Watson-Crick (WC) base pairs. Yet, a non-negligible number of base pairs converts to the Hoogsteen (HG) hydrogen bonding pattern, involving a 180° rotation of the purine base relative to Watson-Crick. These WC to HG conversions alter the conformation of DNA, and may play a role in several processes including recognition and replication. The transient nature of these processes hamper thorough experimental investigation. Molecular dynamics simulations can provide complementary insights to experiments at high spatial and temporal resolution. By using path sampling techniques, a framework that harvests molecular dynamics trajectories that undergo reactions of interest, we avoid long waiting times in stable states, thus focusing on the actual transition. Our results reveal that WC to HG conversion can proceed along different routes with a varying degree of exposure of the purine. Furthermore, we computed the rate constants for this transition using transition interface sampling.

The HAMP domain is a linker region in prokaryotic sensor proteins and relays input signals to the transmitter domain and vice versa. Functional as a dimer, the structure of HAMP shows a parallel coiled-coil motif comprising four helices. To date, it is unclear how HAMP can relay signals from one domain to another, although several models exist. In this work, we use molecular simulation to test the hypothesis that HAMP adopts different conformations, one of which represents an active, signal-relaying configuration, and another an inactive, resting state. We first performed molecular dynamics simulation on the prototype HAMP domain Af1503 from Archaeoglobus fulgidus. We explored its conformational space by taking the structure of the A291F mutant disabling HAMP activity as a starting point. These simulations revealed additional conformational states that differ in the tilt angles between the helices as well as the relative piston shifts of the helices relative to each other. By enhancing the sampling in a metadynamics set up, we investigated three mechanistic models for HAMP signal transduction. Our results indicate that HAMP can access additional conformational states characterized by piston motion. Furthermore, the piston motion of the N-terminal helix of one monomer is directly correlated with the opposite piston motion of the C-terminal helix of the other monomer. The change in piston motion is accompanied by a change in tilt angle between the monomers, thus revealing that HAMP exhibits a collective motion, i.e. a combination of changes in tilt angles and a piston-like displacement. Our results provide insights into the conformational changes that underlie the signaling mechanism involving HAMP., Author Summary For survival, bacteria must constantly monitor their environmental conditions and adapt to these by generating a response. Protein sensors enable bacteria to perceive their surroundings and are typically built from modular compounds that are connected by linker regions. The HAMP domain is such a linker region that relays signals between different modules in a sensory cascade. HAMP is a dimer comprising four helices in a parallel coiled-coil interaction motif. One of the hypotheses explaining the mechanism of signal communication by HAMP is that the domain can adopt different stable conformations. In this work, we used a molecular simulation approach to investigate this hypothesis at high atomic resolution. We found that HAMP can adopt different conformations and that, in doing so, the helices shift and tilt with respect to each other. Furthermore, we found that if one helix moves upward, the helix at the other end in the other monomer moves down.

Bacterial chromosomal DNA is organized within a structure called the nucleoid, which is distinctly different from the rest of the cytoplasm. Bacteria have a number of nucleoid-associated proteins that influence the organization of the nucleoid by bending, wrapping or bridging DNA. The Histone-like Nucleoid Structuring protein H-NS can bridge DNA by binding to two separate DNA duplexes, or shield the DNA by binding to distant sites on the same duplex, depending on external conditions. H-NS occurs in Gram-negative enterobacteria and silences genes involved in bacterial virulence and antibiotic resistance.

Chitotriosidase (CHIT1) is a chitinase that is secreted by activated macrophages. Plasma chitotriosidase activity reflects the presence of lipid-laden macrophages in patients with Gaucher disease. CHIT1 activity can be conveniently measured using fluorogenic 4-methylumbelliferyl (4MU)-chitotrioside or 4MU-chitobioside as the substrate, however, nonsaturating concentrations have to be used because of apparent substrate inhibition. Saturating substrate concentrations can, however, be used with the newly designed substrate 4MU-deoxychitobioside. We studied the impact of a known polymorphism, G102S, on the catalytic properties of CHIT1. The G102S allele was found to be common in type I Gaucher disease patients in the Netherlands ( approximately 24% of alleles). The catalytic efficiency of recombinant Ser102 CHIT1 was approximately 70% that of wild-type Gly102 CHIT1 when measured with 4MU-chitotrioside at a nonsaturating concentration. However, the activity was normal with 4MU-deoxychitobioside as the substrate at saturating concentrations, consistent with predictions from molecular dynamics simulations. In conclusion, interpretation of CHIT1 activity measurements with 4MU-chitotrioside with respect to CHIT1 protein concentrations depends on the presence of Ser102 CHIT1 in an individual, complicating estimation of the body burden of storage macrophages. Use of the superior 4MU-deoxychitobioside substrate avoids such complications because activity towards this substrate under saturating conditions is not affected by the G102S substitution.

The general stress response of Bacillus subtilis can be activated by stimuli such as the addition of salt or ethanol and with blue light. In the latter response, YtvA activates sigma(B) through a cascade of Rsb proteins, organized in stressosomes. YtvA is composed of an N-terminal LOV (light, oxygen, and voltage) domain and a C-terminal STAS (sulfate transporter and anti-sigma factor) domain and shows light-modulated GTP binding in vitro. Here, we examine the mechanism of YtvA-mediated activation of sigma(B) in vivo with site-directed mutagenesis. Constitutive off and constitutive on mutations have been identified. Disruption of GTP binding in the STAS domain eliminates light activation of sigma(B). In contrast, modification of sites relevant for phosphorylation of STAS domains does not affect the stress response significantly. The data obtained are integrated into a model for the structure of full-length YtvA, which presumably functions as a dimer.

When receiving a signal trigger, sensor proteins undergo conformational changes resulting in the formation of a signaling state. Photoactive Yellow Protein (PYP) is a bacterial blue light sensor, 125 amino acids in size, including para-coumaric acid as a chromophore. Upon absorbing a blue-light photon, PYP undergoes a series of rearrangements to form a signaling state. The last step in this process is partial unfolding of the protein, occurring on a sub-millisecond timescale.Molecular simulation can provide detailed insight into the mechanisms underlying protein conformational changes and is complementary to experiments. Studying a protein folding reaction at atomistic resolution with conventional atomistic Molecular Dynamics (MD) is unpractical due to the long time scales involved. These long time scales originate from the presence of local free energy minima from which it is not trivial to escape.Advanced simulations enabled us to investigate the equilibrium characteristics, as well as the dynamical pathways of conformational changes linked to the formation of the signaling state of PYP.Replica exchange MD resulted in the identification of several intermediates during the light induced unfolding. Using these state as input for transition path sampling and subsequent reaction coordinate analysis led to new mechanistic insights in this conformational change. The conformational change starts with the unfolding of a helix in the chromophore binding pocket, followed by the solvent exposure of either the chromophore or glutamate at position 46. Furthermore, our simulations indicate that it is more likely that Glu46 becomes solvent exposed first.To our knowledge this is the first simulation study of unbiased dynamical pathways of a sub-millisecond timescale process of a biologically relevant protein. This work opens up the way for investigating conformational changes in other interesting systems in high detail.

The bacterial sensor Photoactive Yellow Protein (PYP) signals the presence of blue light by undergoing a series of conformational changes. We present atomistic Parallel Tempering (Replica Exchange Molecular Dynamics) simulations of conformational changes occurring during the photo-cycle of PYP. First, we study the signaling state formation of PYP in detail. Our previous simulations have shown that the formation of the signaling state is characterized by the solvent exposure of both the chromophore and Glu46 (Vreede J, Crielaard W, Hellingwerf KJ, Bolhuis PG. Biophys J, 2005;8:3525-3535). Subsequent NMR results agreed with this prediction, but as these experiments were performed on an N-terminally truncated mutant, a simulation of this mutant would further substantiate our previous results. Here, we compare simulations of the truncated PYP to the NMR structures, as well as to the wild type predictions. This comparison also gives some insight into the role of the N-terminal domain of PYP, which restricts the movement of the chromophore binding pocket (CBP) in the wild type. Second, we report simulations of the recovery of the receptor state from the signaling state. While we did not observe complete refolding of the protein, we did observe transient interactions between residues of the CBP occurring when the chromophore is in a trans configuration. Using simulations that sample anomalous exposure of the chromophore in the receptor state, we were able to sample chromophore re-entry into its binding pocket. While the involved time scales prohibit drawing definitive conclusions even when using parallel tempering, we nevertheless propose that the formation of a helix in the CBP is essential for a successful recovery of the receptor state, and forms a kinetic barrier in this process.

Mammals express two active chitinases, chitotriosidase and AMCase. Only AMCase displays an extremely acidic pH optimum, consistent with its observed presence in the gastrointestinal tract. A structural model of AMCase reveals the presence of a conserved histidine residue in the active site. Mutational analyses and molecular dynamics simulations show that His187 is responsible for the acidic optimum and suggest pH dependent modulation of the reaction mechanism that is unique to AMCases. Concluding, His187 is a crucial structural component of the active site of AMCase and this unique feature may serve as a lead for the development of specific inhibitors. (C) 2008 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved

The photoactive yellow protein (PYP) acts as a light sensor to its bacterial host: it responds to light by changing shape. After excitation by blue light, PYP undergoes several transformations, to partially unfold into its signaling state. One of the crucial steps in this photocycle is the protonation of p-coumaric acid after excitation and isomerization of this chromophore. Experimentalists still debate on the nature of the proton donor and on whether it donates the hydrogen directly or indirectly. To obtain better knowledge of the mechanism, we studied this proton transfer using Car-Parrinello molecular dynamics, classical molecular dynamics, and computer simulations combining these two methods (quantum mechanics/molecular mechanics, QMMM). The simulations reproduce the chromophore structure and hydrogen-bond network of the protein measured by X-ray crystallography and NMR. When the chromophore is protonated, it leaves the assumed proton donor, glutamic acid 46, with a negative charge in a hydrophobic environment. We show that the stabilization of this charge is a very important factor in the mechanism of protonation. Protonation frequently occurs in simplified ab initio simulations of the chromophore binding pocket in vacuum, where amino acids can easily hydrogen bond to Glu46. When the complete protein environment is incorporated in a QMMM simulation on the complete protein, no proton transfer is observed within 14 ps. The hydrogen-bond rearrangements in this time span are not sufficient to stabilize the new protonation state. Force field molecular dynamics simulations on a much longer time scale have shown which internal rearrangements of the protein are needed. Combining these simulations with more QMMM calculations enabled us to check the stability of protonation states and clarify the initial requirements for the proton transfer in PYP.