Search

The Eighth Central European Conference “Chemistry towards Biology” was held in Brno, Czech Republic, on 28 August–1 September 2016 to bring together experts in biology, chemistry and design of bioactive compounds promote the exchange of scientific results, methods and ideas and encourage cooperation between researchers from all over the world. The topics of the conference covered “Chemistry towards Biology”, meaning that the event welcomed chemists working on biology-related problems, biologists using chemical methods, and students and other researchers of the respective areas that fall within the common scope of chemistry and biology. The authors of this manuscript are plenary speakers and other participants of the symposium and members of their research teams. The following summary highlights the major points/topics of the meeting.

Membrane type 1-matrix metalloproteinase (MT1-MMP or MMP-14) is a zinc-transmembrane metalloprotease involved in the degradation of extracellular matrix and tumor invasion. While changes in solvation of MT1-MMP have been recently studied, little is known about the structural and energetic changes associated with MT1-MMP while interacting with substrates. Steady-state kinetic and thermodynamic data (including activation energies and activation volumes) were measured over a wide range of temperatures and pressures by means of a stopped-flow fluorescence technique. Complementary temperature- and pressure-dependent Fourier-transform infrared measurements provided corresponding structural information of the protein. MT1-MMP is stable and active over a wide range of temperatures (10–55°C). A small conformational change was detected at 37°C, which is responsible for the change in activity observed at the same temperature. Pressure decreases the enzymatic activity until complete inactivation occurs at 2 kbar. The inactivation is associated with changes in the rate-limiting step of the reaction caused by additional hydration of the active site upon compression and/or minor conformational changes in the active site region. Based on these data, an energy and volume diagram could be established for the various steps of the enzymatic reaction.

The increasing awareness of RNA’s central role in biology calls for a new understanding of how RNAs, like proteins, recognize biological partners. Because RNA is inherently flexible, it assumes a variety of conformations. This conformational flexibility can be a critical aspect of how RNA attracts and binds molecular partners. Structurally, RNA consists of rigid basepaired duplexes, separated by flexible non-basepaired regions. Here, using an RNA system consisting of two short helices, connected by a single-stranded (non-basepaired) junction, we explore the role of helix length and junction sequence in determining the range of conformations available to a model RNA. Single-molecule Förster resonance energy transfer reports on the RNA conformation as a function of either mono- or divalent ion concentration. Electrostatic repulsion between helices dominates at low salt concentration, whereas junction sequence effects determine the conformations at high salt concentration. Near physiological salt concentrations, RNA conformation is sensitive to both helix length and junction sequence, suggesting a means for sensitively tuning RNA conformations.

Eukaryotes express three or more multisubunit nuclear RNA polymerases (Pols) referred to as Pols I, II, and III, each of which synthesizes a specific subset of RNAs. Consistent with the diversity of their target genes, eukaryotic cells have evolved divergent cohorts of transcription factors and enzymatic properties for each RNA polymerase system. Over the years, many trans-acting factors that orchestrate transcription by the individual Pols have been described; however, little effort has been devoted to characterizing the molecular mechanisms of Pol I activity. To begin to address this gap in our understanding of eukaryotic gene expression, here we establish transient-state kinetic approaches to characterize the nucleotide incorporation mechanism of Pol I. We collected time courses for single turnover nucleotide incorporation reactions over a range of substrate ATP concentrations that provide information on both Pol I’s nucleotide addition and nuclease activities. The data were analyzed by model-independent and model-dependent approaches, resulting in, to our knowledge, the first minimal model for the nucleotide addition pathway for Pol I. Using a grid searching approach we provide rigorous bounds on estimated values of the individual elementary rate constants within the proposed model. This work reports the most detailed analysis of Pol I mechanism to date. Furthermore, in addition to their use in transient state kinetic analyses, the computational approaches applied here are broadly applicable to global optimization problems.

Several recent experiments suggest that sharply bent DNA has a surprisingly high bending flexibility, but the cause of this flexibility is poorly understood. Although excitation of flexible defects can explain these results, whether such excitation can occur with the level of DNA bending in these experiments remains unclear. Intriguingly, the DNA contained preexisting nicks in most of these experiments but whether nicks might play a role in flexibility has never been considered in the interpretation of experimental results. Here, using full-atom molecular dynamics simulations, we show that nicks promote DNA basepair disruption at the nicked sites, which drastically reduces DNA bending energy. In addition, lower temperatures suppress the nick-dependent basepair disruption. In the absence of nicks, basepair disruption can also occur but requires a higher level of DNA bending. Therefore, basepair disruption inside B-form DNA can be suppressed if the DNA contains preexisting nicks. Overall, our results suggest that the reported mechanical anomaly of sharply bent DNA is likely dependent on preexisting nicks, therefore the intrinsic mechanisms of sharply bent nick-free DNA remain an open question., Singapore. National Research Foundation (Singapore-MIT Alliance for Research and Technology)

A full understanding of RNA-mediated biology would require the knowledge of three-dimensional (3D) structures, structural flexibility and stability of RNAs. To predict RNA 3D structures and stability, we have previously proposed a three-bead coarse-grained predictive model with implicit salt/solvent potentials. In this study, we will further develop the model by improving the implicit-salt electrostatic potential and involving a sequence-dependent coaxial stacking potential to enable the model to simulate RNA 3D structure folding in divalent/monovalent ion solutions. As compared with the experimental data, the present model can predict 3D structures of RNA hairpins with bulge/internal loops (, 24 pages,8 figures

Ras proteins are small GTPases that act as signal transducers between cell surface receptors and several intracellular signaling cascades. They contain highly homologous catalytic domains and flexible C-terminal hypervariable regions (HVRs) that differ across Ras isoforms. KRAS is among the most frequently mutated oncogenes in human tumors. Surprisingly, we found that the C-terminal HVR of K-Ras4B, thought to minimally impact the catalytic domain, directly interacts with the active site of the protein. The interaction is almost 100-fold tighter with the GDP-bound than the GTP-bound protein. HVR binding interferes with Ras-Raf interaction, modulates binding to phospholipids, and slightly slows down nucleotide exchange. The data indicate that contrary to previously suggested models of K-Ras4B signaling, HVR plays essential roles in regulation of signaling. High affinity binding of short peptide analogs of HVR to K-Ras active site suggests that targeting this surface with inhibitory synthetic molecules for the therapy of KRAS-dependent tumors is feasible.

The number of fluorophores within a molecule complex can be revealed by single-molecule photobleaching imaging. A widely applied strategy to analyze intensity traces over time is the quantification of photobleaching step counts. However, several factors can limit and bias the detection of photobleaching steps, including noise, high numbers of fluorophores, and the possibility that several photobleaching events occur almost simultaneously. In this study, we propose a new approach, to our knowledge, to determine the fluorophore number that correlates the intensity decay of a population of molecule complexes with the decay of the number of visible complexes. We validated our approach using single and fourfold Atto-labeled DNA strands. As an example we estimated the subunit stoichiometry of soluble CD95L using GFP fusion proteins. To assess the precision of our method we performed in silico experiments showing that the estimates are not biased for experimentally observed intensity fluctuations and that the relative precision remains constant with increasing number of fluorophores. In case of fractional fluorescent labeling, our simulations predicted that the fluorophore number estimate corresponds to the product of the true fluorophore number with the labeling fraction. Our method, denoted by spot number and intensity correlation (SONIC), is fully automated, robust to noise, and does not require the counting of photobleaching events.

Metabolite-dependent conformational switching in RNA riboswitches is now widely accepted as a critical regulatory mechanism for gene expression in bacterial systems. More recently, similar gene regulation mechanisms have been found to be important for viral systems as well. One of the most abundant and best-studied systems is the tRNA-like structure (TLS) domain, which has been found to occur in many plant viruses spread across numerous genera. In this work, folding dynamics for the TLS domain of Brome Mosaic Virus have been investigated using single-molecule fluorescence resonance energy transfer techniques. In particular, burst fluorescence methods are exploited to observe metal-ion ([Mn+])-induced folding in freely diffusing RNA constructs resembling the minimal TLS element of brome mosaic virus RNA3. The results of these experiments reveal a complex equilibrium of at least three distinct populations. A stepwise, or consecutive, thermodynamic model for TLS folding is developed, which is in good agreement with the [Mn+]-dependent evolution of conformational populations and existing structural information in the literature. Specifically, this folding pathway explains the metal-ion dependent formation of a functional TLS domain from unfolded RNAs via two consecutive steps: 1) hybridization of a long-range stem interaction, followed by 2) formation of a 3′-terminal pseudoknot. These two conformational transitions are well described by stepwise dissociation constants for [Mg2+] (K1 = 328 ± 30 μM and K2 = 1092 ± 183 μM) and [Na+] (K1 = 74 ± 6 mM and K2 = 243 ± 52 mM)-induced folding. The proposed thermodynamic model is further supported by inhibition studies of the long-range stem interaction using a complementary DNA oligomer, which effectively shifts the dynamic equilibrium toward the unfolded conformation. Implications of this multistep conformational folding mechanism are discussed with regard to regulation of virus replication.

The β-defensins are a class of small cationic proteins that serve as components of numerous systems in vertebrate biology, including the immune and melanocortin systems. Human β-defensin 3 (HBD3), which is produced in the skin, has been found to bind to melanocortin receptors 1 and 4 through complementary electrostatics, a unique mechanism of ligand-receptor interaction. This finding indicates that electrostatics alone, and not specific amino acid contact points, could be sufficient for function in this ligand-receptor system, and further suggests that other small peptide ligands could interact with these receptors in a similar fashion. Here, we conducted molecular-similarity analyses and functional studies of additional members of the human β-defensin family, examining their potential as ligands of melanocortin-1 receptor, through selection based on their electrostatic similarity to HBD3. Using Poisson-Boltzmann electrostatic calculations and molecular-similarity analysis, we identified members of the human β-defensin family that are both similar and dissimilar to HBD3 in terms of electrostatic potential. Synthesis and functional testing of a subset of these β-defensins showed that peptides with an HBD3-like electrostatic character bound to melanocortin receptors with high affinity, whereas those that were anticorrelated to HBD3 showed no binding affinity. These findings expand on the central role of electrostatics in the control of this ligand-receptor system and further demonstrate the utility of employing molecular-similarity analysis. Additionally, we identified several new potential ligands of melanocortin-1 receptor, which may have implications for our understanding of the role defensins play in melanocortin physiology.

Plexins are single-pass transmembrane receptors that bind the axon guidance molecules semaphorins. Single-pass transmembrane proteins are an important class of receptors that display a wide variety of activation mechanisms, often involving ligand-dependent dimerization or conformational changes. Resolving the activation mechanism and dimerization state of these receptors is extremely challenging, especially in a live-cell environment. Here, we report on the dimerization state of PlexinA4 and its response to activation by semaphorin binding. Semaphorins are dimeric molecules that activate plexin by binding two copies of plexin simultaneously and inducing formation of a specific active dimer of plexin. An open question is whether there are preexisting plexin dimers that could act as autoinhibitory complexes. We address these questions with pulsed interleaved excitation fluorescence cross-correlation spectroscopy (PIE-FCCS). PIE-FCCS is a two-color fluorescence microscopy method that is directly sensitive to protein dimerization in a live-cell environment. With PIE-FCCS, we show that inactive PlexinA4 is dimerized in the live-cell plasma membrane. By comparing the cross correlation of full-length PlexinA4 to control proteins and plexin mutants, we show that dimerization of inactive PlexinA4 requires the Sema domain, but not the cytoplasmic domain. Ligand stimulation with Sema6A does not change the degree of cross correlation, indicating that plexin activation does not lead to higher-order oligomerization. Together, the results suggest that semaphorin activates plexin by disrupting an inhibitory plexin dimer and inducing the active dimer.

The observation of biological processes at the molecular scale in real time requires high spatial and temporal resolution. Magnetic tweezers are straightforward to implement, free of radiation or photodamage, and provide ample multiplexing capability, but their spatiotemporal resolution has lagged behind that of other single-molecule manipulation techniques, notably optical tweezers and AFM. Here, we present, to our knowledge, a new high-resolution magnetic tweezers apparatus. We systematically characterize the achievable spatiotemporal resolution for both incoherent and coherent light sources, different types and sizes of beads, and different types and lengths of tethered molecules. Using a bright coherent laser source for illumination and tracking at 6 kHz, we resolve 3 Å steps with a 1 s period for surface-melted beads and 5 Å steps with a 0.5 s period for double-stranded-dsDNA-tethered beads, in good agreement with a model of stochastic bead motion in the magnetic tweezers. We demonstrate how this instrument can be used to monitor the opening and closing of a DNA hairpin on millisecond timescales in real time, together with attendant changes in the hairpin dynamics upon the addition of deoxythymidine triphosphate. Our approach opens up the possibility of observing biological events at submillisecond timescales with subnanometer resolution using camera-based detection.

Proteins exhibit a variety of dense phases ranging from gels, aggregates, and precipitates to crystalline phases and dense liquids. Although the structure of the crystalline phase is known in atomistic detail, little attention has been paid to noncrystalline protein dense phases, and in many cases the structures of these phases are assumed to be fully amorphous. In this work, we used small-angle neutron scattering, electron microscopy, and electron tomography to measure the structure of ovalbumin precipitate particles salted out with ammonium sulfate. We found that the ovalbumin phase-separates into core-shell particles with a core radius of ∼2 μm and shell thickness of ∼0.5 μm. Within this shell region, nanostructures comprised of crystallites of ovalbumin self-assemble into a well-defined bicontinuous network with branches ∼12 nm thick. These results demonstrate that the protein gel is comprised in part of nanocrystalline protein.

Intrinsically disordered proteins and intrinsically disordered regions are implicated in many biological functions and associated with many diseases, but their conformational characterizations are challenging. The disordered β6/β7 loop of Staphylococcus aureus sortase A is involved in the binding of both sorting signals and calcium. Calcium binding allosterically activates the enzyme, but the detailed mechanism has been unclear. Here we adapted the replica exchange with solute tempering method to sample the conformations of the β6/β7 loop, in apo form and in three liganded forms (bound with a sorting signal or calcium or both). The extensive molecular dynamics simulations yield atomic details of the disordered-to-order transition of the loop and suggest a mechanism for allosteric activation: calcium binding results in partial closure and ordering of the loop, thereby leading to preorganization of the binding pocket for the sorting signal. The approach has general applicability to the study of intrinsically disordered regions.

The Engrailed Homeodomain (EnHD) transcription factor of Drosophila melanogaster was fused to the enhanced green fluorescent protein (eGFP) either at its C- or N-terminus via three- or ten-residue flexible linkers. Here, we show that EnHD undergoes destabilization upon fusing it to eGFP regardless of the linker length used and whether the tethering is to its N- or C-terminus. The destabilization is reflected in melting points that are lower by up to 9°C. Thermodynamic analysis and coarse-grained molecular dynamic simulations indicate that this destabilization is due to eGFP-promoted entropic stabilization of the denatured state ensemble of EnHD. Our results provide, therefore, an example for destabilizing interdomain allostery. They are also important given the widespread use of eGFP tagging in cell biology, as they indicate that such tagging can cause unintended protein destabilization and concomitant effects.

The glucocorticoid receptor (GR) is a steroid-hormone-activated transcription factor that modulates gene expression. Transcriptional regulation by the GR requires dynamic receptor binding to specific target sites located across the genome. This binding remodels the chromatin structure to allow interaction with other transcription factors. Thus, chromatin remodeling is an essential component of GR-mediated transcriptional regulation, and understanding the interactions between these molecules at the structural level provides insights into the mechanisms of how GR and chromatin remodeling cooperate to regulate gene expression. This study suggests models for the assembly of the SWI/SNF-A (SWItch/Sucrose-NonFermentable) complex and its interaction with the GR. We used the PRISM algorithm (PRotein Interactions by Structural Matching) to predict the three-dimensional complex structures of the target proteins. The structural models indicate that BAF57 and/or BAF250 mediate the interaction between the GR and the SWI/SNF-A complex, corroborating experimental data. They further suggest that a BAF60a/BAF155 and/or BAF60a/BAF170 interaction is critical for association between the core and variant subunits. Further, we model the interaction between GR and CCAAT-enhancer-binding proteins (C/EBPs), since the GR can regulate gene expression indirectly by interacting with other transcription factors like C/EBPs. We observe that GR can bind to bZip domains of the C/EBPα homodimer as both a monomer and dimer of the DNA-binding domain. In silico mutagenesis of the predicted interface residues confirm the importance of these residues in binding. In vivo analysis of the computationally suggested mutations reveals that double mutations of the leucine residues (L317D+L335D) may disrupt the interaction between GR and C/EBPα. Determination of the complex structures of the GR is of fundamental relevance to understanding its interactions and functions, since the function of a protein or a complex is dictated by its structure. In addition, it may help us estimate the effects of mutations on GR interactions and signaling. Fil: Muratcioglu, Serena. Koc University; Turquía Fil: Presman, Diego Martin. National Institutes of Health; Estados Unidos. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Fisiología, Biología Molecular y Neurociencias. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Fisiología, Biología Molecular y Neurociencias; Argentina Fil: Pooley, John R.. University Of Bristol; Reino Unido Fil: Grøntved, Lars. University Of Southern Denmark; Dinamarca Fil: Hager, Gordon L.. National Institutes of Health; Estados Unidos Fil: Nussinov, Ruth. Tel Aviv University; Israel Fil: Keskin, Ozlem. Koc University; Turquía Fil: Gursoy, Attila. Koc University; Turquía

The tunnel region at triosephosphate isomerase (TIM)’s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM’s catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand’s effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop’s dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.

We examine the contrast between mechanisms for allosteric signaling that involve structural change, and those that do not, from the perspective of allosteric pathways. In particular we treat in detail the case of fluctuation-allostery by which amplitude modulation of the thermal fluctuations of the elastic normal modes conveys the allosteric signal, and address the question of what an allosteric pathway means in this case. We find that a perturbation theory of thermal elastic solids and nonperturbative approach (by super-coarse-graining elasticity into internal bending modes) have opposite signatures in their structure of correlated pathways. We illustrate the effect from analysis of previous results from GlxR of Corynebacterium glutamicum, an example of the CRP/FNR transcription family of allosteric homodimers. We find that the visibility of both correlated pathways and disconnected sites of correlated motion in this protein suggests that mechanisms of local elastic stretch and bend are recruited for the purpose of creating and controlling allosteric cooperativity.

Proteins have a highly dynamic nature and there is a complex interrelation between their structural dynamics and binding behavior. By assuming various conformational ensembles, they perform both local and global fluctuations to interact with other proteins in a dynamic infrastructure adapted to functional motion. Here, we show that there is a significant association between allosteric mutations, which lead to high-binding-affinity changes, and the hinge positions of global modes, as revealed by a large-scale statistical analysis of data in the Structural Kinetic and Energetic Database of Mutant Protein Interactions (SKEMPI). We further examined the mechanism of allosteric dynamics by conducting studies on human growth hormone (hGH) and pyrin domain (PYD), and the results show how mutations at the hinge regions could allosterically affect the binding-site dynamics or induce alternative binding modes by modifying the ensemble of accessible conformations. The long-range dissemination of perturbations in local chemistry or physical interactions through an impact on global dynamics can restore the allosteric dynamics. Our findings suggest a mechanism for the coupling of structural dynamics to the modulation of protein interactions, which remains a critical phenomenon in understanding the effect of mutations that lead to functional changes in proteins.

The glucagon-like peptide 1 (GLP-1) analog, liraglutide, is a GLP-1 agonist and is used in the treatment of type-2 diabetes mellitus and obesity. From a pharmaceutical perspective, it is important to know the oligomerization state of liraglutide with respect to stability. Compared to GLP-1, liraglutide has an added fatty acid (FA) moiety that causes oligomerization of liraglutide as suggested by small-angle x-ray scattering (SAXS) and multiangle static light scattering (MALS) results. SAXS data suggested a global shape of a hollow elliptical cylinder of size hexa-, hepta-, or octamer, whereas MALS data indicate a hexamer. To elaborate further on the stability of these oligomers and the role of the FA chains, a series of molecular-dynamics simulations were carried out on 11 different hexa-, hepta-, and octameric systems. Our results indicate that interactions of the fatty acid chains contribute noticeably to the stabilization. The simulation results indicate that the heptamer with paired FA chains is the most stable oligomer when compared to the 10 other investigated structures. Theoretical SAXS curves extracted from the simulations qualitatively agree with the experimentally determined SAXS curves supporting the view that liraglutide forms heptamers in solution. In agreement with the SAXS data, the heptamer forms a water-filled oligomer of elliptical cylindrical shape.

PDEδ is a small protein that binds and controls the trafficking of RAS subfamily proteins. Its inhibition protects initiation of RAS signaling, and it is one of the common targets considered for oncological drug development. In this study, we used solved x-ray structures of inhibitor-bound PDEδ targets to investigate mechanisms of action of six independent all-atom MD simulations. An analysis of atomic simulations combined with the molecular mechanic-Poisson-Boltzmann solvent accessible surface area/generalized Born solvent accessible surface area calculations led to the identification of action mechanisms for a panel of novel PDEδ inhibitors. To the best of our knowledge, this study is one of the first in silico investigations on co-crystallized PDEδ protein. A detailed atomic-scale understanding of the molecular mechanism of PDEδ inhibition may assist in the design of novel PDEδ inhibitors. One of the most common side effects for diverse small molecules/kinase inhibitors is their off-target interactions with cardiac ion channels and human-ether-a-go-go channel specifically. Thus, all of the studied PDEδ inhibitors are also screened in silico at the central cavities of hERG1 potassium channels.

A balance of van der Waals, electrostatic, and hydrophobic forces drive the folding and packing of protein side chains. Although such interactions between residues are often approximated as being pairwise additive, in reality, higher-order many-body contributions that depend on environment drive hydrophobic collapse and cooperative electrostatics. Beginning from dead-end elimination, we derive the first algorithm, to our knowledge, capable of deterministic global repacking of side chains compatible with many-body energy functions. The approach is applied to seven PCNA x-ray crystallographic data sets with resolutions 2.5–3.8 Å (mean 3.0 Å) using an open-source software. While PDB_REDO models average an Rfree value of 29.5% and MOLPROBITY score of 2.71 Å (77th percentile), dead-end elimination with the polarizable AMOEBA force field lowered Rfree by 2.8–26.7% and improved mean MOLPROBITY score to atomic resolution at 1.25 Å (100th percentile). For structural biology applications that depend on side-chain repacking, including x-ray refinement, homology modeling, and protein design, the accuracy limitations of pairwise additivity can now be eliminated via polarizable or quantum mechanical potentials.

Sedimentation equilibrium (SE) analytical ultracentrifugation (AUC) is a gold standard for the rigorous determination of macromolecular buoyant molar masses and the thermodynamic study of reversible interactions in solution. A significant experimental drawback is the long time required to attain SE, which is usually on the order of days. We have developed a method for time-optimized SE (toSE) with defined time-varying centrifugal fields that allow SE to be attained in a significantly (up to 10-fold) shorter time than is usually required. To achieve this, numerical Lamm equation solutions for sedimentation in time-varying fields are computed based on initial estimates of macromolecular transport properties. A parameterized rotor-speed schedule is optimized with the goal of achieving a minimal time to equilibrium while limiting transient sample preconcentration at the base of the solution column. The resulting rotor-speed schedule may include multiple over- and underspeeding phases, balancing the formation of gradients from strong sedimentation fluxes with periods of high diffusional transport. The computation is carried out in a new software program called TOSE, which also facilitates convenient experimental implementation. Further, we extend AUC data analysis to sedimentation processes in such time-varying centrifugal fields. Due to the initially high centrifugal fields in toSE and the resulting strong migration, it is possible to extract sedimentation coefficient distributions from the early data. This can provide better estimates of the size of macromolecular complexes and report on sample homogeneity early on, which may be used to further refine the prediction of the rotor-speed schedule. In this manner, the toSE experiment can be adapted in real time to the system under study, maximizing both the information content and the time efficiency of SE experiments.

G protein-coupled receptors are vital membrane proteins that allosterically transduce biomolecular signals across the cell membrane. However, the process by which ligand binding induces protein conformation changes is not well understood biophysically. Rhodopsin, the mammalian dim-light receptor, is a unique test case for understanding these processes because of its switch-like activity; the ligand, retinal, is bound throughout the activation cycle, switching from inverse agonist to agonist after absorbing a photon. By contrast, the ligand-free opsin is outside the activation cycle and may behave differently. We find that retinal influences rhodopsin dynamics using an ensemble of all-atom molecular dynamics simulations that in aggregate contain 100 μs of sampling. Active retinal destabilizes the inactive state of the receptor, whereas the active ensemble was more structurally homogenous. By contrast, simulations of an active-like receptor without retinal present were much more heterogeneous than those containing retinal. These results suggest allosteric processes are more complicated than a ligand inducing protein conformational changes or simply capturing a shifted ensemble as outlined in classic models of allostery.

Spectral diversity of fluorescent proteins, crucial for multiparameter imaging, is based mainly on chemical diversity of their chromophores. Recently we have reported, to our knowledge, a new green fluorescent protein WasCFP—the first fluorescent protein with a tryptophan-based chromophore in the anionic state. However, only a small portion of WasCFP molecules exists in the anionic state at physiological conditions. In this study we report on an improved variant of WasCFP, named NowGFP, with the anionic form dominating at 37°C and neutral pH. It is 30% brighter than enhanced green fluorescent protein (EGFP) and exhibits a fluorescence lifetime of 5.1 ns. We demonstrated that signals of NowGFP and EGFP can be clearly distinguished by fluorescence lifetime in various models, including mammalian cells, mouse tumor xenograft, and Drosophila larvae. NowGFP thus provides an additional channel for multiparameter fluorescence lifetime imaging microscopy of green fluorescent proteins.

Supercoiled DNA polymer models for which the torsional energy depends on the total twist of molecules (Tw) are a priori well suited for thermodynamic analysis of long molecules. So far, nevertheless, the exact determination of Tw in these models has been based on a computation of the writhe of the molecules (Wr) by exploiting the conservation of the linking number, Lk=Tw+Wr, which reflects topological constraints coming from the helical nature of DNA. Because Wr is equal to the number of times the main axis of a DNA molecule winds around itself, current Monte Carlo algorithms have a quadratic time complexity, O(L2), with respect to the contour length (L) of the molecules. Here, we present an efficient method to compute Tw exactly, leading in principle to algorithms with a linear complexity, which in practice is O(L1.2). Specifically, we use a discrete wormlike chain that includes the explicit double-helix structure of DNA and where the linking number is conserved by continuously preventing the generation of twist between any two consecutive cylinders of the discretized chain. As an application, we show that long (up to 21 kbp) linear molecules stretched by mechanical forces akin to magnetic tweezers contain, in the buckling regime, multiple and branched plectonemes that often coexist with curls and helices, and whose length and number are in good agreement with experiments. By attaching the ends of the molecules to a reservoir of twists with which these can exchange helix turns, we also show how to compute the torques in these models. As an example, we report values that are in good agreement with experiments and that concern the longest molecules that have been studied so far (16 kbp).

Biopolymer homeostasis underlies the health of organisms, and protective osmolytes have emerged as one strat- egy used by Nature to preserve biopolymer homeostasis. However, a great deal remains unknown about the mechanism of action of osmolytes. Trehalose, as a prominent example, stabilizes proteins against denaturation by extreme temperature and denaturants, preserves membrane integrity upon freezing or in dry conditions, inhibits polyQ-mediated protein aggregation, and suppresses the aggregation of denatured proteins. The underlying thermodynamic mechanisms of such diverse effects of trehalose remain unclear or controversial. In this study, we applied the surface-additive method developed in the Record labo- ratory to attack this issue. We characterized the key features of trehalose-biopolymer preferential interactions and found that trehalose has strong unfavorable interactions with aliphatic carbon and significant favorable interactions with amide/anionic oxygen. This dissection has allowed us to elucidate the diverse effects of trehalose and to identify the crucial functional group(s) responsible for its effects. With (semi)quantitative thermodynamic analysis, we discovered that 1) the unfavorable interaction of trehalose with hydrophobic surfaces is the dominant factor in its effect on protein stability, 2) the favorable interaction of treha- lose with polar amides enables it to inhibit polyQ-mediated protein aggregation and the aggregation of denatured protein in gen- eral, and 3) the favorable interaction of trehalose with phosphate oxygens, together with its unfavorable interaction with aliphatic carbons, enables trehalose to preserve membrane integrity in aqueous solution. These results provide a basis for a full under- standing of the role of trehalose in biopolymer homeostasis and the reason behind its evolutionary selection as an osmolyte, as well as for a better application of trehalose as a chemical chaperone.

Neuronal calcium sensor-1 (NCS-1) is the primordial member of a family of proteins responsible primarily for sensing changes in neuronal Ca(2+) concentration. NCS-1 is a multispecific protein interacting with a number of binding partners in both calcium-dependent and independent manners, and acting in a variety of cellular processes in which it has been linked to a number of disorders such as schizophrenia and autism. Despite extensive studies on the Ca(2+)-activated state of NCS proteins, little is known about the conformational dynamics of the Mg(2+)-bound and apo states, both of which are populated, at least transiently, at resting Ca(2+) conditions. Here, we used optical tweezers to study the folding behavior of individual NCS-1 molecules in the presence of Mg(2+) and in the absence of divalent ions. Under tension, the Mg(2+)-bound state of NCS-1 unfolds and refolds in a three-state process by populating one intermediate state consisting of a folded C-domain and an unfolded N-domain. The interconversion at equilibrium between the different molecular states populated by NCS-1 was monitored in real time through constant-force measurements and the energy landscapes underlying the observed transitions were reconstructed through hidden Markov model analysis. Unlike what has been observed with the Ca(2+)-bound state, the presence of Mg(2+) allows both the N- and C-domain to fold through all-or-none transitions with similar refolding rates. In the absence of divalent ions, NCS-1 unfolds and refolds reversibly in a two-state reaction involving only the C-domain, whereas the N-domain has no detectable transitions. Overall, the results allowed us to trace the progression of NCS-1 folding along its energy landscapes and provided a solid platform for understanding the conformational dynamics of similar EF-hand proteins.

In Alzheimer’s disease (AD), the amyloid β (Aβ) peptide aggregates in the brain to form progressively larger oligomers, fibrils, and plaques. The aggregation process is strongly influenced by the presence of other macromolecular species, called crowders, that can exert forces on the proteins. One very common attribute of macromolecular crowders is their hydrophobicity. We examined the effect of hydrophobic crowders on protein aggregation by using discontinuous molecular dynamics (DMD) simulations in combination with an intermediate resolution protein model, PRIME20. The systems considered contained 48 Aβ (16–22) peptides and crowders with diameters of 5 Å, 20 Å, and 40 Å, represented by hard spheres or spheres with square-well/square-shoulder interactions, at a crowder volume fraction of ϕ = 0.10. Results show that low levels of crowder hydrophobicity are capable of increasing the fibrillation lag time and high levels of crowder hydrophobicity can fully prevent the formation of fibrils. The types of structures that remain during the final stages of the simulations are summarized in a global phase diagram that shows fibril, disordered oligomer, or β-sheet phases in the space spanned by crowder size and crowder hydrophobicity. In particular, at high levels of hydrophobicity, simulations with 5 Å crowders result in only disordered oligomers and simulations with 40 Å crowders result in only β-sheets. The presence of hydrophobic crowders reduces the antiparallel β-sheet content of fibrils, whereas hard sphere crowders increase it. Finally, strong hydrophobic crowders alter the secondary structure of the Aβ (16–22) monomers, bending them into a shape that is incapable of forming ordered β-sheets or fibrils. These results qualitatively agree with previous theoretical and experimental work.

The antibody immunoglobulin (Ig) 2D1 is effective against the 1918 hemagglutinin (HA) and also known to cross-neutralize the 2009 pandemic H1N1 influenza HA through a similar epitope. However, the detailed mechanism of neutralization remains unclear. We conducted molecular dynamics (MD) simulations to study the interactions between Ig-2D1 and the HAs from the 1918 pandemic flu (A/South Carolina/1/1918, 18HA), the 2009 pandemic flu (A/California/04/2009, 09HA), a 2009 pandemic flu mutant (A/California/04/2009, 09HA_mut), and the 2006 seasonal flu (A/Solomon Islands/3/2006, 06HA). MM-PBSA analyses suggest the approximate free energy of binding (ΔG) between Ig-2D1 and 18HA is −74.4 kcal/mol. In comparison with 18HA, 09HA and 06HA bind Ig-2D1 ∼6 kcal/mol (ΔΔG) weaker, and the 09HA_mut bind Ig-2D1 only half as strong. We also analyzed the contributions of individual epitope residues using the free-energy decomposition method. Two important salt bridges are found between the HAs and Ig-2D1. In 09HA, a serine-to-asparagine mutation coincided with a salt bridge destabilization, hydrogen bond losses, and a water pocket formation between 09HA and Ig-2D1. In 09HA_mut, a lysine-to-glutamic-acid mutation leads to the loss of both salt bridges and destabilizes interactions with Ig-2D1. Even though 06HA has a similar ΔG to 09HA, it is not recognized by Ig-2D1 in vivo. Because 06HA contains two potential glycosylation sites that could mask the epitope, our results suggest that Ig-2D1 may be active against 06HA only in the absence of glycosylation. Overall, our simulation results are in good agreement with observations from biological experiments and offer novel mechanistic insights, to our knowledge, into the immune escape of the influenza virus.

Recently various pathways of human telomere (ht) DNA folding into G-quadruplexes and of ligand binding to these structures have been proposed. However, the key issue as to the nature of forces driving the folding and recognition processes remains unanswered. In this study, structural changes of 22-mer ht-DNA fragment (Tel22), induced by binding of ions (K+, Na+) and specific bisquinolinium ligands, were monitored by calorimetric and spectroscopic methods and by gel electrophoresis. Using the global model analysis of a wide variety of experimental data, we were able to characterize the thermodynamic forces that govern the formation of stable Tel22 G-quadruplexes, folding intermediates, and ligand-quadruplex complexes, and then predict Tel22 behavior in aqueous solutions as a function of temperature, salt concentration, and ligand concentration. On the basis of the above, we believe that our work sets the framework for better understanding the heterogeneity of ht-DNA folding and binding pathways, and its structural polymorphism.

Septins comprise a family of proteins involved in a variety of cellular processes and related to several human pathologies. They are constituted by three structural domains: the N- and C-terminal domains, highly variable in length and composition, and the central domain, involved in the guanine nucleotide (GTP) binding. Thirteen different human septins are known to form heterogeneous complexes or homofilaments, which are stabilized by specific interactions between the different interfaces present in the domains. In this work, we have investigated by in-solution small-angle x-ray scattering the structural and thermodynamic properties of a human septin 3 construct, SEPT3-GC, which contains both of both interfaces (G and NC) responsible for septin-septin interactions. In order to shed light on the role of these interactions, small-angle x-ray scattering measurements were performed in a wide range of temperatures, from 2 up to 56°C, both with and without a nonhydrolysable form of GTP (GTPγS). The acquired data show a temperature-dependent coexistence of monomers, dimers, and higher-order aggregates that were analyzed using a global fitting approach, taking into account the crystallographic structure of the recently reported SEPT3 dimer, PDB:3SOP. As a result, the enthalpy, entropy, and heat capacity variations that control the dimer-monomer dissociation equilibrium in solution were derived and GTPγS was detected to increase the enthalpic stability of the dimeric species. Moreover, a temperature increase was observed to induce dissociation of SEPT3-GC dimers into monomers just preceding their reassembling into amyloid aggregates, as revealed by the Thioflavin-T fluorescence assays.

Both in vivo and in vitro, specific sequences in double-stranded DNA can adopt the left-handed Z-form when underwound. Recently, the B-Z transition of DNA has been studied in detail in magnetic tweezers experiments by several groups. We present a theoretical description of this transition, based on an annealed random copolymer model. The transition of a switchable sequence is discussed as a function of energetic and geometric parameters of the B- and Z-forms, of the applied boundary conditions, and of the characteristics of the B-Z interface. We address a possible torsional softening upon the B-Z transition. The model can be also applied to other biofilaments with annealed torsional/flexural degrees of freedom.

We have captured the binding of a peptide to a PDZ domain by unbiased molecular dynamics simulations. Analysis of the trajectories reveals on-pathway encounter complex formation, which is driven by electrostatic interactions between negatively charged carboxylate groups in the peptide and positively charged side chains surrounding the binding site. In contrast, the final stereospecific complex, which matches the crystal structure, features completely different interactions, namely the burial of the hydrophobic side chain of the peptide C-terminal residue and backbone hydrogen bonds. The simulations show that nonnative salt bridges stabilize kinetically the encounter complex during binding. Unbinding follows the inverse sequence of events with the same nonnative salt bridges in the encounter complex. Thus, in contrast to protein folding, which is driven by native interactions, the binding of charged peptides can be steered by nonnative interactions, which might be a general mechanism, e.g., in the recognition of histone tails by bromodomains.

Nanopore sensing involves an electrophoretic transport of analytes through a nanoscale pore, permitting label-free sensing at the single-molecule level. However, to date, the detection of individual small proteins has been challenging, primarily due to the poor signal/noise ratio that these molecules produce during passage through the pore. Here, we show that fine adjustment of the buffer pH, close to the isoelectric point, can be used to slow down the translocation speed of the analytes, hence permitting sensing and characterization of small globular proteins. Ubiquitin (Ub) is a small protein of 8.5 kDa, which is well conserved in all eukaryotes. Ub conjugates to proteins as a posttranslational modification called ubiquitination. The immense diversity of Ub substrates, as well as the complexity of Ub modification types and the numerous physiological consequences of these modifications, make Ub and Ub chains an interesting and challenging subject of study. The ability to detect Ub and to identify Ub linkage type at the single-molecule level may provide a novel tool for investigation in the Ub field. This is especially adequate because, for most ubiquitinated substrates, Ub modifies only a few molecules in the cell at a given time. Applying our method to the detection of mono- and poly-Ub molecules, we show that we can analyze their characteristics using nanopores. Of particular importance is that two Ub dimers that are equal in molecular weight but differ in 3D structure due to their different linkage types can be readily discriminated. Thus, to our knowledge, our method offers a novel approach for analyzing proteins in unprecedented detail using solid-state nanopores. Specifically, it provides the basis for development of single-molecule sensing of differently ubiquitinated substrates with different biological significance. Finally, our study serves as a proof of concept for approaching nanopore detection of sub-10-kDa proteins and demonstrates the ability of this method to differentiate among native and untethered proteins of the same mass.

Specific details concerning the spatial organization of nucleosomes in 30 nm fibers remain unknown. To investigate this, we analyzed all stereochemically possible configurations of two-start nucleosome fibers with short DNA linkers L = 13–37 bp (nucleosome repeat length (NRL) = 160–184 bp). Four superhelical parameters—inclination of nucleosomes, twist, rise, and diameter—uniquely describe a regular symmetric fiber. The energy of a fiber is defined as the sum of four terms: elastic energy of the linker DNA, steric repulsion, electrostatics, and a phenomenological (H4 tail–acidic patch) interaction between two stacked nucleosomes. By optimizing the fiber energy with respect to the superhelical parameters, we found two types of topological transition in fibers (associated with the change in inclination angle): one caused by an abrupt 360° change in the linker DNA twisting (change in the DNA linking number, ΔLk = 1), and another caused by overcrossing of the linkers ( ΔLk = 2). To the best of our knowledge, this topological polymorphism of the two-start fibers was not reported in the computations published earlier. Importantly, the optimal configurations of the fibers with linkers L = 10 n and 10 n + 5 bp are characterized by different values of the DNA linking number—that is, they are topologically different. Our results are consistent with experimental observations, such as the inclination 60° to 70° (the angle between the nucleosomal disks and the fiber axis), helical rise, diameter, and left-handedness of the fibers. In addition, we make several testable predictions, among them different degrees of DNA supercoiling in fibers with L = 10 n and 10 n + 5 bp, different flexibility of the two types of fibers, and a correlation between the local NRL and the level of transcription in different parts of the yeast genome.

The aggregation of α-synuclein (α-Syn) is linked to Parkinson’s disease. The mechanism of early aggregation steps and the effect of pathogenic single-point mutations remain elusive. We report here a single-molecule fluorescence study of α-Syn dimerization and the effect of mutations. Specific interactions between tethered fluorophore-free α-Syn monomers on a substrate and fluorophore-labeled monomers diffusing freely in solution were observed using total internal reflection fluorescence microscopy. The results showed that wild-type (WT) α-Syn dimers adopt two types of dimers. The lifetimes of type 1 and type 2 dimers were determined to be 197 ± 3 ms and 3334 ± 145 ms, respectively. All three of the mutations used, A30P, E46K, and A53T, increased the lifetime of type 1 dimer and enhanced the relative population of type 2 dimer, with type 1 dimer constituting the major fraction. The kinetic stability of type 1 dimers (expressed in terms of lifetime) followed the order A30P (693 ± 14 ms) > E46K (292 ± 5 ms) > A53T (226 ± 6 ms) > WT (197 ± 3 ms). Type 2 dimers, which are more stable, had lifetimes in the range of several seconds. The strongest effect, observed for the A30P mutant, resulted in a lifetime 3.5 times higher than observed for the WT type 1 dimer. This mutation also doubled the relative fraction of type 2 dimer. These data show that single-point mutations promote dimerization, and they suggest that the structural heterogeneity of α-Syn dimers could lead to different aggregation pathways.

The impact of flow disturbances on platelet adhesion is complex and incompletely understood. At the molecular scale, platelet glycoprotein Ibα (GPIbα) must associate with the von Willebrand factor A1 domain (VWF-A1) with a rapid on-rate under high hemodynamic forces, as occurs in arterial thrombosis, where various transport mechanisms are at work. Here, we theoretically modeled the coupled transport-reaction process of the two-dimensional (2D) receptor-ligand association kinetics in a biomembrane force probe to explicitly account for the effects of molecular length, confinement stiffness, medium viscosity, surface curvature, and separation distance. We experimentally verified the theoretical approach by visualizing association and dissociation of individual VWF-A1-GPIbα bonds in a real-time thermal fluctuation assay. The apparent on-rate, reciprocal of the average time intervals between sequential bonds, decreased with the increasing gap distance between A1- and GPIbα-bearing surfaces with an 80-nm threshold (beyond which bond formation became prohibitive) identified as the combined contour length of the receptor and ligand molecules. The biomembrane force probe spring constant and diffusivity of the protein-bearing beads also significantly influenced the apparent on-rate, in accordance with the proposed transport mechanisms. The global agreement between the experimental data and the model predictions supports the hypothesis that receptor-ligand association behaves distinctly in the transport- and reaction-limited scenarios. To our knowledge, our results represent the first detailed quantification of physical regulation of the 2D on-rate that allows platelets to sense and respond to local changes in their hemodynamic environment. In addition, they provide an approach for determining the intrinsic kinetic parameters that employs simultaneous experimental measurements and theoretical modeling of bond association in a single assay. The 2D intrinsic forward rate for VWF-A1-GPIbα association was determined from the measurements to be (3.5 ± 0.67) × 10−4 μm2 s−1.

The RNA-dependent RNA polymerase from the Hepatitis C Virus (gene product NS5B) is a validated drug target because of its critical role in genome replication. There are at least four distinct allosteric sites on the polymerase to which several small molecule inhibitors bind. In addition, numerous crystal structures have been solved with different allosteric inhibitors bound to the polymerase. However, the molecular mechanisms by which these small molecules inhibit the enzyme have not been fully elucidated. There is evidence that allosteric inhibitors alter the intrinsic motions and distribution of conformations sampled by the enzyme. In this study we use molecular dynamics simulations to understand the structural and dynamic changes that result when inhibitors are bound at three different allosteric binding sites on the enzyme. We observe that ligand binding at each site alters the structure and dynamics of NS5B in a distinct manner. Nonetheless, our studies also highlight commonalities in the mechanisms of action of the different inhibitors. Each inhibitor alters the conformational states sampled by the enzyme, either by rigidifying the enzyme and preventing transitions between functional conformational states or by destabilizing the enzyme and preventing functionally relevant conformations from being adequately sampled. By illuminating the molecular mechanisms of allosteric inhibition, these studies delineate the intrinsic functional properties of the enzyme and pave the way for designing novel and more effective polymerase inhibitors. This information may also be important to understand how allosteric regulation occurs in related viral polymerases and other enzymes.

PAP248–286 is a 39-residue fragment (residues 248 to 286) derived from protease cleavage of prostatic acidic phosphatase in semen. The amyloid fibrils formed in vitro by PAP248–286 can dramatically enhance human immunodeficiency virus (HIV) infection. To our knowledge, we present the first report that the HIV-enhancing potency of fibrils formed by PAP248–286 is morphology dependent. We identified pleomorphic fibrils by transmission electron microscopy in two buffer conditions. Our solid-state NMR data showed that these fibrils consist of molecules in distinct conformations. In agreement with NMR, fluorescence measurements confirmed that they are assembled along different pathways, with distinct molecular structures. Furthermore, our cell-based infectivity tests detected distinct HIV-enhancing potencies for fibrils in distinct morphologies. In addition, our transmission electron microscopy and NMR results showed that semen-derived enhancer of viral infection fibrils formed in sodium bicarbonate buffer remain stable over time, but semen-derived enhancer of viral infection fibrils formed in phosphate buffered saline keep evolving after the initial 7 days incubation period. Given time, most of the assemblies in phosphate buffered saline will turn into elongated thin fibrils. They have similar secondary structure but different packing than thin fibrils formed initially after 7 days incubation.

One and the same protein can self-assemble into amyloid fibrils with different morphologies. The phenomenon of fibril polymorphism is relevant biologically because different fibril polymorphs can have different toxicity, but there is no tool for predicting which polymorph forms and under what conditions. Here, we consider the nucleation of polymorphic amyloid fibrils occurring by direct polymerization of monomeric proteins into fibrils. We treat this process within the framework of our newly developed nonstandard nucleation theory, which allows prediction of the concentration dependence of the nucleation rate for different fibril polymorphs. The results highlight that the concentration dependence of the nucleation rate is closely linked with the protein solubility and a threshold monomer concentration below which fibril formation becomes biologically irrelevant. The relation between the nucleation rate, the fibril solubility, the threshold concentration, and the binding energies of the fibril building blocks within fibrils might prove a valuable tool for designing new experiments to control the formation of particular fibril polymorphs.

Although the magnitude of a protein’s net charge (Z) can control its rate of self-assembly into amyloid, and its interactions with cellular membranes, the net charge of a protein is not viewed as a druggable parameter. This article demonstrates that aspirin (the quintessential acylating pharmacon) can inhibit the amyloidogenesis of superoxide dismutase (SOD1) by increasing the intrinsic net negative charge of the polypeptide, i.e., by acetylation (neutralization) of multiple lysines. The protective effects of acetylation were diminished (but not abolished) in 100 mM NaCl and were statistically significant: a total of 432 thioflavin-T amyloid assays were performed for all studied proteins. The acetylation of as few as three lysines by aspirin in A4V apo-SOD1—a variant that causes familial amyotrophic lateral sclerosis (ALS)—delayed amyloid nucleation by 38% and slowed amyloid propagation by twofold. Lysines in wild-type- and ALS-variant apo-SOD1 could also be peracetylated with aspirin after fibrillization, resulting in supercharged fibrils, with increases in formal net charge of ∼2 million units. Peracetylated SOD1 amyloid defibrillized at temperatures below unacetylated fibrils, and below the melting temperature of native Cu2,Zn2-SOD1 (e.g., fibril Tm = 84.49°C for acetylated D90A apo-SOD1 fibrils). Targeting the net charge of native or misfolded proteins with small molecules—analogous to how an enzyme’s Km or Vmax are medicinally targeted—holds promise as a strategy in the design of therapies for diseases linked to protein self-assembly.

International audience; Diflavin reductases are bidomain electron transfer proteins in which structural reorientation is necessary to account for the various intramolecular and intermolecular electron transfer steps. Using small-angle x-ray scattering and nuclear magnetic resonance data, we describe the conformational free-energy landscape of the NADPH-cytochrome P450 reductase (CPR), a typical bidomain redox enzyme composed of two covalently-bound flavin domains, under various experimental conditions. The CPR enzyme exists in a salt- and pH-dependent rapid equilibrium between a previously described rigid, locked state and a newly characterized, highly flexible, unlocked state. We further establish that maximal electron flux through CPR is conditioned by adjustable stability of the locked-state domain interface under resting conditions. This is rationalized by a kinetic scheme coupling rapid conformational sampling and slow chemical reaction rates. Regulated domain interface stability associated with fast stochastic domain contacts during the catalytic cycle thus provides, to our knowledge, a new paradigm for improving our understanding of multidomain enzyme function.

Single-particle cryo-electron microscopy is widely used to study the structure of macromolecular assemblies. Tens of thousands of noisy two-dimensional images of the macromolecular assembly viewed from different directions are used to infer its three-dimensional structure. The first step is to estimate a low-resolution initial model and initial image orientations. This is a challenging global optimization problem with many unknowns, including an unknown orientation for each two-dimensional image. Obtaining a good initial model is crucial for the success of the subsequent refinement step. We introduce a probabilistic algorithm for estimating an initial model. The algorithm is fast, has very few algorithmic parameters, and yields information about the precision of estimated model parameters in addition to the parameters themselves. Our algorithm uses a pseudo-atomic model to represent the low-resolution three-dimensional structure, with isotropic Gaussian components as moveable pseudo-atoms. This leads to a significant reduction in the number of parameters needed to represent the three-dimensional structure, and a simplified way of computing two-dimensional projections. It also contributes to the speed of the algorithm. We combine the estimation of the unknown three-dimensional structure and image orientations in a Bayesian framework. This ensures that there are very few parameters to set, and specifies how to combine different types of prior information about the structure with the given data in a systematic way. To estimate the model parameters we use Markov chain Monte Carlo sampling. The advantage is that instead of just obtaining point estimates of model parameters, we obtain an ensemble of models revealing the precision of the estimated parameters. We demonstrate the algorithm on both simulated and real data.

Nitroxide spin labels are used for double electron-electron resonance (DEER) measurements of distances between sites in biomolecules. Rotation of gem-dimethyls in commonly used nitroxides causes spin echo dephasing times (Tm) to be too short to perform DEER measurements at temperatures between ∼80 and 295 K, even in immobilized samples. A spirocyclohexyl spin label has been prepared that has longer Tm between 80 and 295 K in immobilized samples than conventional labels. Two of the spirocyclohexyl labels were attached to sites on T4 lysozyme introduced by site-directed spin labeling. Interspin distances up to ∼4 nm were measured by DEER at temperatures up to 160 K in water/glycerol glasses. In a glassy trehalose matrix the Tm for the doubly labeled T4 lysozyme was long enough to measure an interspin distance of 3.2 nm at 295 K, which could not be measured for the same protein labeled with the conventional 1-oxyl-2,2,5,5-tetramethyl-3-pyrroline-3-(methyl)methanethio-sulfonate label.

The formation of linear protein fibrils has previously been shown to be enhanced by volume exclusion or crowding in the presence of a high concentration of chemically inert protein or polymer, and by adsorption to membrane surfaces. An equilibrium mesoscopic model for the combined effect of both crowding and adsorption upon the fibrillation of a dilute tracer protein is presented. The model exhibits behavior that differs qualitatively from that observed in the presence of crowding or adsorption alone. The model predicts that in a crowded solution, at critical values of the volume fraction of crowder or intrinsic energy of the tracer-wall interaction, the tracer protein will undergo an extremely cooperative transition—approaching a step function—from existence as a slightly self-associated species in solution to existence as a highly self-associated and completely adsorbed species. Criteria for a valid experimental test of these predictions are presented.

Single-molecule fluorescence microscopy is a powerful tool for observing biomolecular interactions with high spatial and temporal resolution. Detecting fluorescent signals from individual labeled proteins above high levels of background fluorescence remains challenging, however. For this reason, the concentrations of labeled proteins in in vitro assays are often kept low compared to their in vivo concentrations. Here, we present a new fluorescence imaging technique by which single fluorescent molecules can be observed in real time at high, physiologically relevant concentrations. The technique requires a protein and its macromolecular substrate to be labeled each with a different fluorophore. Making use of short-distance energy-transfer mechanisms, only the fluorescence from those proteins that bind to their substrate is activated. This approach is demonstrated by labeling a DNA substrate with an intercalating stain, exciting the stain, and using energy transfer from the stain to activate the fluorescence of only those labeled DNA-binding proteins bound to the DNA. Such an experimental design allowed us to observe the sequence-independent interaction of Cy5-labeled interferon-inducible protein 16 with DNA and the sliding via one-dimensional diffusion of Cy5-labeled adenovirus protease on DNA in the presence of a background of hundreds of nanomolar Cy5 fluorophore.

Protein structural analysis demonstrates that water molecules are commonly found in the internal cavities of proteins. Analysis of experimental data on the entropies of inorganic crystals suggests that the entropic cost of transferring such a water molecule to a protein cavity will not typically be greater than 7.0 cal/mol/K per water molecule, corresponding to a contribution of approximately +2.0 kcal/mol to the free energy. In this study, we employ the statistical mechanical method of inhomogeneous fluid solvation theory to quantify the enthalpic and entropic contributions of individual water molecules in 19 protein cavities across five different proteins. We utilize information theory to develop a rigorous estimate of the total two-particle entropy, yielding a complete framework to calculate hydration free energies. We show that predictions from inhomogeneous fluid solvation theory are in excellent agreement with predictions from free energy perturbation (FEP) and that these predictions are consistent with experimental estimates. However, the results suggest that water molecules in protein cavities containing charged residues may be subject to entropy changes that contribute more than +2.0 kcal/mol to the free energy. In all cases, these unfavorable entropy changes are predicted to be dominated by highly favorable enthalpy changes. These findings are relevant to the study of bridging water molecules at protein-protein interfaces as well as in complexes with cognate ligands and small-molecule inhibitors.

The bacterial outer membrane comprises two main classes of components, lipids and membrane proteins. These nonsoluble compounds are conveyed across the aqueous periplasm along specific molecular transport routes: the lipid lipopolysaccharide (LPS) is shuttled by the Lpt system, whereas outer membrane proteins (Omps) are transported by chaperones, including the periplasmic Skp. In this study, we revisit the specificity of the chaperone-lipid interaction of Skp and LPS. High-resolution NMR spectroscopy measurements indicate that LPS interacts with Skp nonspecifically, accompanied by destabilization of the Skp trimer and similar to denaturation by the nonnatural detergent lauryldimethylamine-N-oxide (LDAO). Bioinformatic analysis of amino acid conservation, structural analysis of LPS-binding proteins, and MD simulations further confirm the absence of a specific LPS binding site on Skp, making a biological relevance of the interaction unlikely. Instead, our analysis reveals a highly conserved salt-bridge network, which likely has a role for Skp function.

Measurements from hydroxyl radical footprinting (HRF) provide rich information about the solvent accessibility of amino acid side chains of a protein. Traditional HRF data analyses focus on comparing the difference in the modification/footprinting rate of a specific site to infer structural changes across two protein states, e.g., between a free and ligand-bound state. However, the rate information itself is not fully used for the purpose of comparing different protein sites within a protein on an absolute scale. To provide such a cross-site comparison, we present a new, to our knowledge, data analysis algorithm to convert the measured footprinting rate constant to a protection factor (PF) by taking into account the known intrinsic reactivity of amino acid side chain. To examine the extent to which PFs can be used for structural interpretation, this PF analysis is applied to three model systems where radiolytic footprinting data are reported in the literature. By visualizing structures colored with the PF values for individual peptides, a rational view of the structural features of various protein sites regarding their solvent accessibility is revealed, where high-PF regions are buried and low-PF regions are more exposed to the solvent. Furthermore, a detailed analysis correlating solvent accessibility and local structural contacts for gelsolin shows a statistically significant agreement between PF values and various structure measures, demonstrating that the PFs derived from this PF analysis readily explain fundamental HRF rate measurements. We also tested this PF analysis on alternative, chemical-based HRF data, showing improved correlations of structural properties of a model protein barstar compared to examining HRF rate data alone. Together, this PF analysis not only permits a novel, to our knowledge, approach of mapping protein structures by using footprinting data, but also elevates the use of HRF measurements from a qualitative, cross-state comparison to a quantitative, cross-site assessment of protein structures in the context of individual conformational states of interest.