Monte Carlo simulation algorithm for B‐DNA

Understanding the structure–function relationship of biomolecules containing DNA has motivated experiments aimed at determining molecular structure using methods such as small‐angle X‐ray and neutron scattering (SAXS and SANS). SAXS and SANS are useful for determining macromolecular shape in solution, a process which benefits by using atomistic models that reproduce the scattering data. The variety of algorithms available for creating and modifying model DNA structures lack the ability to rapidly modify all‐atom models to generate structure ensembles. This article describes a Monte Carlo algorithm for simulating DNA, not with the goal of predicting an equilibrium structure, but rather to generate an ensemble of plausible structures which can be filtered using experimental results to identify a sub‐ensemble of conformations that reproduce the solution scattering of DNA macromolecules. The algorithm generates an ensemble of atomic structures through an iterative cycle in which B‐DNA is represented using a wormlike bead–rod model, new configurations are generated by sampling bend and twist moves, then atomic detail is recovered by back mapping from the final coarse‐grained configuration. Using this algorithm on commodity computing hardware, one can rapidly generate an ensemble of atomic level models, each model representing a physically realistic configuration that could be further studied using molecular dynamics. © 2016 Wiley Periodicals, Inc.     

Ab initio computational modeling of long loops in G-protein coupled receptors

A newly developed approach for predicting the structure of segments that connect known elements of secondary structure in proteins has been applied to some of the longer loops in the G-protein coupled receptors (GPCRs) rhodopsin and the dopamine receptor D2R. The algorithm uses Monte Carlo (MC) simulation in a temperature annealing protocol combined with a scaled collective variables (SCV) technique to search conformation space for loop structures that could belong to the native ensemble. Except for rhodopsin, structural information is only available for the transmembrane helices (TMHs), and therefore the usual approach of finding a single conformation of lowest energy has to be abandoned. Instead the MC search aims to find the ensemble located at the absolute minimum free energy, i.e., the native ensemble. It is assumed that structures in the native ensemble can be found by an MC search starting from any conformation in the native funnel. The hypothesis is that native structures are trapped in this part of conformational space because of the high-energy barriers that surround the native funnel. In this work it is shown that the crystal structure of the second extracellular loop (e2) of rhodopsin is a member of this loop’s native ensemble. In contrast, the crystal structure of the third intracellular loop is quite different in the different crystal structures that have been reported. Our calculations indicate, that of three crystal structures examined, two show features characteristic of native ensembles while the other one does not. Finally the protocol is used to calculate the structure of the e2 loop in D2R. Here, the crystal structure is not known, but it is shown that several side chains that are involved in interaction with a class of substituted benzamides assume conformations that point into the active site. Thus, they are poised to interact with the incoming ligand.     

Folding transition-state and denatured-state ensembles of FSD-1 from folding and unfolding simulations

Characterization of the folding transition-state ensemble and the denatured-state ensemble is an important step toward a full elucidation of protein folding mechanisms. We report herein an investigation of the free-energy landscape of FSD-1 protein by a total of four sets of folding and unfolding molecular dynamics simulations with explicit solvent. The transition-state ensemble was initially identified from unfolding simulations at 500 K and was verified by simulations at 300 K starting from the ensemble structures. The denatured-state ensemble and the early-stage folding were studied by a combination of unfolding simulations at 500 K and folding simulations at 300 K starting from the extended conformation. A common feature of the transition-state ensemble was the substantial formation of the native secondary structures, including both the alpha-helix and beta-sheet, with partial exposure of the hydrophobic core in the solvent. Both the native and non-native secondary structures were observed in the denatured-state ensemble and early-stage folding, consistent with the smooth experimental melting curve. Interestingly, the contact orders of the transition-state ensemble structures were similar to that of the native structure and were notably lower than those of the compact structures found in early-stage folding, implying that chain and topological entropy might play significant roles in protein folding. Implications for FSD-1 folding mechanisms and the rate-limiting step are discussed. Analyses further revealed interesting non-native interactions in the denatured-state ensemble and early-stage folding and the possibility that destabilization of these interactions could help to enhance the stability and folding rate of the protein.     

Structural characterization of flexible proteins using small-angle X-ray scattering

Structural analysis of flexible macromolecular systems such as intrinsically disordered or multidomain proteins with flexible linkers is a difficult task as high-resolution techniques are barely applicable. A new approach, ensemble optimization method (EOM), is proposed to quantitatively characterize flexible proteins in solution using small-angle X-ray scattering (SAXS). The flexibility is taken into account by allowing for the coexistence of different conformations of the protein contributing to the experimental scattering pattern. These conformers are selected using a genetic algorithm from a pool containing a large number of randomly generated models covering the protein configurational space. Quantitative criteria are developed to analyze the EOM selected models and to determine the optimum number of conformers in the ensemble. Simultaneous fitting of multiple scattering patterns from deletion mutants, if available, provides yet more detailed local information about the structure. The efficiency of EOM is demonstrated in model and practical examples on completely or partially unfolded proteins and on multidomain proteins interconnected by linkers. In the latter case, EOM is able to distinguish between rigid and flexible proteins and to directly assess the interdomain contacts.     

Structure of Cytochromec and Its Platinum-modified Derivatives by Fourier Transform Infrared Spectroscopy

ＳｔｒｕｃｔｕｒｅｏｆＣｙｔｏｃｈｒｏｍｅｃａｎｄＩｔｓＰｌａｔｉｎｕｍ┐ｍｏｄｉｆｉｅｄＤｅｒｉｖａｔｉｖｅｓｂｙＦｏｕｒｉｅｒＴｒａｎｓｆｏｒｍＩｎｆｒａｒｅｄＳｐｅｃｔｒｏｓｃｏｐｙＪＩＡＮＧＬｉ－ｊｕａｎａ，ＳＵＮＷｅｉ－ｙｉｎａ，ＦＡＮＧＪ...     

Explaining the structural plasticity of α-synuclein

Given that α-synuclein has been implicated in the pathogenesis of several neurodegenerative disorders, deciphering the structure of this protein is of particular importance. While monomeric α-synuclein is disordered in solution, it can form aggregates rich in cross-β structure, relatively long helical segments when bound to micelles or lipid vesicles, and a relatively ordered helical tetramer within the native cell environment. To understand the physical basis underlying this structural plasticity, we generated an ensemble for monomeric α-synuclein using a Bayesian formalism that combines data from NMR chemical shifts, RDCs, and SAXS with molecular simulations. An analysis of the resulting ensemble suggests that a non-negligible fraction of the ensemble (0.08, 95% confidence interval 0.03-0.12) places the minimal toxic aggregation-prone segment in α-synuclein, NAC(8-18), in a solvent exposed and extended conformation that can form cross-β structure. Our data also suggest that a sizable fraction of structures in the ensemble (0.14, 95% confidence interval 0.04-0.23) contains long-range contacts between the N- and C-termini. Moreover, a significant fraction of structures that contain these long-range contacts also place the NAC(8-18) segment in a solvent exposed orientation, a finding in contrast to the theory that such long-range contacts help to prevent aggregation. Lastly, our data suggest that α-synuclein samples structures with amphipathic helices that can self-associate via hydrophobic contacts to form tetrameric structures. Overall, these observations represent a comprehensive view of the unfolded ensemble of monomeric α-synuclein and explain how different conformations can arise from the monomeric protein.     

The implication of chemical extraction treatments on the cell wall nanostructure of softwood

Spruce wood was subjected to well-defined extraction treatments with sodium chlorite (NaClO₂) for delignification, as well as with sodium hydroxide (NaOH) at different concentrations for extraction of hemicelluloses. The corresponding changes of the macromolecular polymer assembly were investigated by small-angle X-ray scattering (SAXS). Measurements with Fourier-transform infrared (FTIR) spectroscopy and wide-angle scattering (WAXS) gave qualitative information about the effectiveness of the extraction process, while the scattering experiments provided information about the regularity and typical dimensions of the molecular structures. The scattering data indicated that delignification had only a moderate effect on the structural organisation of the cell wall, while further extraction with NaOH induced considerable nanostructural changes.     

间规1,2-聚丁二烯结晶结构研究

以间规1,2-聚丁二烯(s-PB)为研究对象,通过原位同步辐射小角X射线散射(SR-SAXS)和广角X射线衍射(WAXD)研究其结晶结构的变化过程.SR-SAXS曲线中存在明显的散射峰,表明在等温结晶过程中形成有序结晶结构;在等温结晶后间规1,2-聚丁二烯的片晶厚度、微晶尺寸均正比于1/Tc∞-T,根据高分子结晶中介相机理可以做出合理的解释.     

Local structure of pseudoboehmites

Summary The structure of pseudoboehmite was studied by the method of wide-angle X-ray scattering (WAXS). The samples of pseudoboehmite with various characteristics produced by different methods were examined in order to obtain the most complete information about their structure. In contrast to reference data, it was found that modifications in the structure of pseudoboehmite layers occur viathe insertion of additional water molecules.     

Structure of amorphous poly-(ethylmethacrylate): a wide-angle x-ray scattering study

The local, temperature dependent structure of poly-(ethylmethacrylate) was studied with wide-angle x-ray scattering. The results were set into context to recent wide-angle neutron scattering results as well as to the results from a multitude of studies of the dynamics of this polymer. The temperature dependence of the wide angle x-ray results point to the development of local order which is backed by the neutron scattering results and which is connected to characteristic temperatures of the relaxation dynamic T(g) and T(c). The poly-(ethylmethacrylate) was studied in its predominantly syndiotactic as well as predominantly isotactic state displaying vastly different local structures as manifested in the x-ray results.     

Contrast-Matched Small-Angle X-ray Scattering from a Heavy-Atom-Labeled Protein in Structure Determination: Application to a Lead-Substituted Calmodulin-Peptide Complex

The information content in 1-D solution X-ray scattering profiles is generally restricted to low-resolution shape and size information that, on its own, cannot lead to unique 3-D structures of biological macromolecules comparable to all-atom models derived from X-ray crystallography or NMR spectroscopy. Here we show that contrast-matched X-ray scattering data collected on a protein incorporating specific heavy-atom labels in 65% aqueous sucrose buffer can dramatically enhance the power of conventional small- and wide-angle X-ray scattering (SAXS/WAXS) measurements. Under contrast-matching conditions the protein is effectively invisible and the main contribution to the X-ray scattering intensity arises from the heavy atoms, allowing direct extraction of pairwise distances between them. In combination with conventional aqueous SAXS/WAXS data, supplemented by NMR-derived residual dipolar couplings (RDCs) measured in a weakly aligning medium, we show that it is possible to position protein domains relative to one another within a precision of 1 ?. We demonstrate this approach with respect to the determination of domain positions in a complex between calmodulin, in which the four Ca(2+) ions have been substituted by Pb(2+), and a target peptide. The uniqueness of the resulting solution is established by an exhaustive search over all models compatible with the experimental data, and could not have been achieved using aqueous SAXS and RDC data alone. Moreover, we show that the correct structural solution can be recovered using only contrast-matched SAXS and aqueous SAXS/WAXS data.     

Resting and active states of the ERK2:HePTP complex

The MAP kinase ERK2 (ERK2, extracellular signal-regulated kinase 2) is regulated by numerous phosphatases that tightly control its activity. For example, the hematopoietic tyrosine phosphatase (HePTP) negatively regulates T cell activation in lymphocytes via ERK2 dephosphorylation. However, only very limited structural information is available for these biologically important complexes. Here, we use small-angle X-ray scattering combined with EROS ensemble refinement to characterize the structures of the resting and active states of ERK2:HePTP complexes. Our data show that the resting state ERK2:HePTP complex adopts a highly extended, dynamic conformation that becomes compact and ordered in the active state complex. This work experimentally demonstrates that these complexes undergo significant dynamic structural changes in solution and provides the first structural insight into an active state MAPK complex.     

Solution-phase structural characterization of supramolecular assemblies by molecular diffraction

Structures of four molecular squares based on rhenium coordination chemistry have been characterized in the solution phase using pair distribution function (PDF) analyses of wide-angle X-ray scattering measured to better than 1 A spatial resolution. In this report we have focused, in particular, on a comparison of structures for pyrazine- and bipyridine-edged squares measured in solution with structures determined for these molecules in the solid state using X-ray crystallography and models derived from geometry optimization and molecular dynamics simulations using a classical force field. The wide-angle scattering for these assemblies is dominated by pair correlations involving one or more rhenium atoms, with both edge and diagonal Re-Re interactions appearing prominently in PDF plots. The pyrazine square is characterized by a relatively rigid structure in solution, with PDF peak positions and linewidths corresponding closely to those calculated from crystal structure data. For the bipyridine-edged square, the experimental PDF peaks measured along the molecular sides match the positions and linewidths of the PDF peaks calculated from static models. In contrast, PDF peaks measured across the diagonal distances of the molecular square deviate significantly from those calculated from the static crystallographic and energy minimized models. The experimental data are instead indicative of configurational broadening of the diagonal distances. In this respect, molecular dynamics simulations point to the importance of butterfly type motions that modulate the Re-Re diagonal distance. Indeed, the experimental data are reasonably well fit by assuming a bimodal distribution of butterfly conformers differing by approximately 25 degrees in the Re-Re-Re-Re torsion angle. Additionally, the measurements provide evidence for solvent ordering by the supramolecular assemblies detected as regions of solvent association and exclusion.     

RNACluster: An integrated tool for RNA secondary structure comparison and clustering

RNA structure comparison is a fundamental problem in structural biology, structural chemistry, and bioinformatics. It can be used for analysis of RNA energy landscapes, conformational switches, and facilitating RNA structure prediction. The purpose of our integrated tool RNACluster is twofold: to provide a platform for computing and comparison of different distances between RNA secondary structures, and to perform cluster identification to derive useful information of RNA structure ensembles, using a minimum spanning tree (MST) based clustering algorithm. RNACluster employs a cluster identification approach based on a MST representation of the RNA ensemble data and currently supports six distance measures between RNA secondary structures. RNACluster provides a user-friendly graphical interface to allow a user to compare different structural distances, analyze the structure ensembles, and visualize predicted structural clusters.     

Structure of triplex DNA in the gas phase

Extensive (more than 90 microseconds) molecular dynamics simulations complemented with ion-mobility mass spectrometry experiments have been used to characterize the conformational ensemble of DNA triplexes in the gas phase. Our results suggest that the ensemble of DNA triplex structures in the gas phase is well-defined over the experimental time scale, with the three strands tightly bound, and for the most abundant charge states it samples conformations only slightly more compact than the solution structure. The degree of structural alteration is however very significant, mimicking that found in duplex and much larger than that suggested for G-quadruplexes. Our data strongly supports that the gas phase triplex maintains an excellent memory of the solution structure, well-preserved helicity, and a significant number of native contacts. Once again, a linear, flexible, and charged polymer as DNA surprises us for its ability to retain three-dimensional structure in the absence of solvent. Results argue against the generally assumed roles of the different physical interactions (solvent screening of phosphate repulsion, hydrophobic effect, and solvation of accessible polar groups) in modulating the stability of DNA structures.     

The elastic net algorithm and protein structure prediction

Predicting protein structures from their amino acid sequences is a problem of global optimization. Global optima (native structures) are often sought using stochastic sampling methods such as Monte Carlo or molecular dynamics, but these methods are slow. In contrast, there are fast deterministic methods that find near-optimal solutions of well-known global optimization problems such as the traveling salesman problem (TSP). But fast TSP strategies have yet to be applied to protein folding, because of fundamental differences in the two types of problems. Here, we show how protein folding can be framed in terms of the TSP, to which we apply a variation of the Durbin-Willshaw elastic net optimization strategy. We illustrate using a simple model of proteins with database-derived statistical potentials and predicted secondary structure restraints. This optimization strategy can be applied to many different models and potential functions, and can readily incorporate experimental restraint information. It is also fast; with the simple model used here, the method finds structures that are within 5-6 A all-Calpha-atom RMSD of the known native structures for 40-mers in about 8 s on a PC; 100-mers take about 20 s. The computer time tau scales as tau approximately n, where n is the number of amino acids. This method may prove to be useful for structure refinement and prediction.     

Structure and reactivity of bis(iodozincio)methane solution

Bis(iodozincio)methane, which has been shown to be an efficient reagent for organic synthesis, is obtained as THF solution. The structural information about the reagent as THF solution was corrected by small angle neutron scattering and by anomalous X-ray scattering. Those scattering experiments implied that the prepared bis(iodozincio)methane exists without forming any oligomer or aggregate. A coordination of tetrahydrothiophene to bis(iodozincio)methane enhances the nucleophilicity of the reagent and stabilizes its monomeric structure in the solution.     

A Multifaceted Secondary Structure Mimic Based On Piperidine‐piperidinones

Minimalist secondary structure mimics are typically made to resemble one interface in a protein–protein interaction (PPI), and thus perturb it. We recently proposed suitable chemotypes can be matched with interface regions directly, without regard for secondary structures. Here we describe a modular synthesis of a new chemotype 1 , simulation of its solution‐state conformational ensemble, and correlation of that with ideal secondary structures and real interface regions in PPIs. Scaffold 1 presents amino acid side‐chains that are quite separated from each other, in orientations that closely resemble ideal sheet or helical structures, similar non‐ideal structures at PPI interfaces, and regions of other PPI interfaces where the mimic conformation does not resemble any secondary structure. 68 different PPIs where conformations of 1 matched well were identified. A new method is also presented to determine the relevance of a minimalist mimic crystal structure to its solution conformations. Thus dld ‐ 1 faf crystallized in a conformation that is estimated to be 0.91 kcal mol^?1 above the minimum energy solution state.     

Tertiary interactions determine the accuracy of RNA folding

RNAs must fold into unique three-dimensional structures to function in the cell, but how each polynucleotide finds its native structure is not understood. To investigate whether the stability of the tertiary structure determines the speed and accuracy of RNA folding, docking of a tetraloop with its receptor in a bacterial group I ribozyme was perturbed by site-directed mutagenesis. Disruption of the tetraloop or its receptor destabilizes tertiary interactions throughout the ribozyme by 2-3 kcal/mol, demonstrating that tertiary interactions form cooperatively in the transition from a native-like intermediate to the native state. Nondenaturing PAGE and RNase T1 digestion showed that base pairs form less homogeneously in the mutant RNAs during the transition from the unfolded state to the intermediate. Thus, tertiary interactions between helices bias the ensemble of secondary structures toward native-like conformations. Time-resolved hydroxyl radical footprinting showed that the wild-type ribozyme folds completely within 5-20 ms. By contrast, only 40-60% of a tetraloop mutant ribozyme folds in 30-40 ms, with the remainder folding in 30-200 s via nonnative intermediates. Therefore, destabilization of tetraloop-receptor docking introduces an alternate folding pathway in the otherwise smooth energy landscape of the wild-type ribozyme. Our results show that stable tertiary structure increases the flux through folding pathways that lead directly and rapidly to the native structure.     

Stable DOPG/Glycyrrhizin Vesicles with a Wide Range of Mixing Ratios: Structure and Stability as Seen by Scattering Experiments and Cryo-TEM

Phosphatidylglycerols represent a large share of the lipids in the plasmamembrane of procaryotes. Therefore, this study investigates the role of charged lipids in the plasma membrane with respect to the interaction of the antiviral saponin glycyrrhizin with such membranes. Glycyrrhizin is a natural triterpenic-based surfactant found in licorice. Vesicles made of 1,2-dioleoyl-sn-glycero-3-phospho-rac-(1’-glycerol) (DOPG)/glycyrrhizin are characterized by small-angle scattering with neutrons and X-rays (SANS and SAXS). Small-angle scattering data are first evaluated by the model-independent modified Kratky–Porod method and afterwards fitted by a model describing the shape of small unilamellar vesicles (SUV) with an internal head-tail contrast. Complete miscibility of DOPG and glycyrrhizin was revealed even at a ratio of lipid:saponin of 1:1. Additional information about the chain-chain correlation distance of the lipid/saponin mixtures in the SUV structures is obtained from wide-angle X-ray scattering (WAXS).