Macromolecular structure modeling from 3D EM using VolRover 2.0

We review tools for structure identification and model-based refinement from three-dimensional electron microscopy implemented in our in-house software package, VOLROVER 2.0. For viral density maps with icosahedral symmetry, we segment the capsid, polymeric, and monomeric subunits using techniques based on automatic symmetry detection and multidomain fast marching. For large biomolecules without symmetry information, we again use our multidomain fast-marching method with manual or fit-based multiseeding to segment meaningful substructures. In either case, we subject the resulting segmented subunit to secondary structure detection when the EM resolution is sufficiently high, and rigid-body structure fitting when the corresponding X-ray structure is available. Secondary structure elements are identified by three techniques: our earlier volume-based and boundary-based skeletonization methods as well as a new method, currently in development, based on solving the grassfire flow equation. For rigid-body fitting, we adapt our earlier fast Fourier-based correlation scheme F2Dock. Our reported segmentation, secondary structure elements identification, and rigid-body fitting techniques, implemented in VOLROVER 2.0 are applied to the PSB 2011 cryo-EM modeling challenge data, and our results are briefly compared to similar results submitted from other research groups. The comparisons show that our techniques are equally capable of segmenting relatively accurate subunits from a viral or protein assembly, and that high segmentation quality leads in turn to higher-quality results of secondary structure elements identification and correlation-based rigid-body fitting. ? 2012 Wiley Periodicals, Inc. Biopolymers 97: 709-731, 2012.     

Using motion planning to study protein folding pathways.

We present a framework for studying protein folding pathways and potential landscapes which is based on techniques recently developed in the robotics motion planning community. Our focus in this work is to study the protein folding mechanism assuming we know the native fold. That is, instead of performing fold prediction, we aim to study issues related to the folding process, such as the formation of secondary and tertiary structure, and the dependence of the folding pathway on the initial denatured conformation. Our work uses probabilistic roadmap (PRM) motion planning techniques which have proven successful for problems involving high-dimensional configuration spaces. A strength of these methods is their efficiency in rapidly covering the planning space without becoming trapped in local minima. We have applied our PRM technique to several small proteins (~60 residues) and validated the pathways computed by comparing the secondary structure formation order on our paths to known hydrogen exchange experimental results. An advantage of the PRM framework over other simulation methods is that it enables one to easily and efficiently compute folding pathways from any denatured starting state to the (known) native fold. This aspect makes our approach ideal for studying global properties of the protein's potential landscape, most of which are difficult to simulate and study with other methods. For example, in the proteins we study, the folding pathways starting from different denatured states sometimes share common portions when they are close to the native fold, and moreover, the formation order of the secondary structure appears largely independent of the starting denatured conformation. Another feature of our technique is that the distribution of the sampled conformations is correlated with the formation of secondary structure and, in particular, appears to differentiate situations in which secondary structure clearly forms first and those in which the tertiary structure is obtained more directly. Overall, our results applying PRM techniques are very encouraging and indicate the promise of our approach for studying proteins for which experimental results are not available.     

Comparison of the Predicted Structures of Loops in the ras–SOS Protein Bound to a Single ras-p21 Protein with the Crystallographically Determined Structures in SOS Bound to Two ras-p21 Proteins

We have previously computed the structures of three loops, residues 591–596, 654–675 and 742–751, in the ras-p21 protein-binding domain (residues 568–1044) of the guanine nucleotide-exchange-promoting SOS protein that were crystallographically undefined when one molecule of ras-p21 (unbound to nucleotide) binds to SOS. Based on our computational results, we synthesized three peptides corresponding to sequences of each of these three loops and found that all three peptides strongly inhibit ras-p21 signaling. More recently, a new crystal structure of SOS has been determined in which this protein binds to two molecules of ras-p21, one unbound to GTP and one bound to GTP. In this structure, the 654–675 loop and residues 742–743 and 750–751 are now crystallographically defined. We have superimposed our energy-minimized structure of the ras-binding domain of SOS bound to one molecule of ras-p21 on the X-ray structure for SOS bound to two molecules of ras-p21. We find that, while the two structures are superimposable, there are large deviations of the residues 673 and 676 and 741 and 752, flanking the two loop segments. This suggests that the binding of the extra ras-p21 molecule, which is far from each of the three loops, induces conformational changes in these domains and further supports their role in signal transduction. In spite of these differences, we have superimposed our computed structures for the loop residues on those from the more recent X-ray structure. Our structure for the 654–675 segment is an anti-parallel beta-sheet with a reverse turn at residues 663–665; in the X-ray structure residues 655–662 adopt an alpha-helical conformation; on the other hand, our computed structure for residues 663–675 superimpose on the X-ray structure for these residues. We further find that our computed structures for residues 742–743 and 750–751 are superimposable on the X-ray structure for these residues.     

NMR solution structure and topological orientation of monomeric phospholamban in dodecylphosphocholine micelles

Phospholamban is an integral membrane protein that regulates the contractility of cardiac muscle by maintaining cardiomyocyte calcium homeostasis. Abnormalities in association of protein kinase A with PLB have recently been linked to human heart failure, where a single mutation is responsible for dilated cardiomyopathy. To date, a high-resolution structure of phospholamban in a lipid environment has been elusive. Here, we describe the first structure of recombinant, monomeric, biologically active phospholamban in lipid-mimicking dodecylphosphocholine micelles as determined by multidimensional NMR experiments. The overall structure of phospholamban is "L-shaped" with the hydrophobic domain approximately perpendicular to the cytoplasmic portion. This is in agreement with our previously published solid-state NMR data. In addition, there are two striking discrepancies between our structure and those reported previously for synthetic phospholamban in organic solvents: a), in our structure, the orientation of the cytoplasmic helix is consistent with the amphipathic nature of these residues; and b), within the hydrophobic helix, residues are positioned on two discrete faces of the helix as consistent with their functional roles ascribed by mutagenesis. This topology renders the two phosphorylation sites, Ser-16 and Thr-17, more accessible to kinases.     

A coarse graining approach to determine nucleic acid structures from small angle neutron scattering profiles in solution

We present a theoretical method to calculate the small angle neutron scattering profile of nucleic acid structures in solution. Our approach is sensitive to the sequence and the structure of the nucleic acid. In order to test our approach, we apply this method to the calculation of the experimental scattered intensity of the decamer d(CCAACGTTGG)₂ in H₂O. This sequence was specifically chosen for this study as it is believed to adopt a canonical B-form structure in 0.3 M NaCl. We find that not only will our methodology reproduce the experimental scattered intensity for this sequence, but our method will also discriminate between B-, A- and Z-form DNA. By studying the scattering profile of this structure in 0.5 and 1.0 M NaCl, we are also able to identify tetraplex and other similar oligomers formation and to model the complex using the experimental scattering data in conjunction with our methodology.     

Vibrational analysis of crystalline tri-L-alanine

We have found that tri-L-alanine (Ala3) can crystallize in a parallel-chain beta structure in addition to the previously known antiparallel-chain beta structure. Although the chain conformations in each structure are essentially similar, the ir and Raman spectra are distinctively different. We have calculated the normal modes of each structure, and can account in significant detail for these differences. This demonstrates the essential validity of our empirically refined force fields, as well as showing that deeper insights into polypeptide and protein structure can be achieved through the rigorous analyses of normal mode calculations.     

Statistical Computations Underlying the Dynamics of Memory Updating

Psychophysical and neurophysiological studies have suggested that memory is not simply a carbon copy of our experience: Memories are modified or new memories are formed depending on the dynamic structure of our experience, and specifically, on how gradually or abruptly the world changes. We present a statistical theory of memory formation in a dynamic environment, based on a nonparametric generalization of the switching Kalman filter. We show that this theory can qualitatively account for several psychophysical and neural phenomena, and present results of a new visual memory experiment aimed at testing the theory directly. Our experimental findings suggest that humans can use temporal discontinuities in the structure of the environment to determine when to form new memory traces. The statistical perspective we offer provides a coherent account of the conditions under which new experience is integrated into an old memory versus forming a new memory, and shows that memory formation depends on inferences about the underlying structure of our experience.     

Homology-modeled structure of the yeast mitochondrial citrate transport protein

We have used homology modeling to construct a three-dimensional model of the yeast mitochondrial citrate transport protein (CTP), based on the recently published x-ray crystal structure of another mitochondrial transport protein, the ADP/ATP carrier. Superposition of the backbone traces of the homology-modeled CTP onto the crystallographically determined ADP carrier structure indicates that the CTP transmembrane domains are well modeled (i.e., root mean square deviation of 0.94 A), whereas the loops facing the intermembrane space and the mitochondrial matrix are less certain (i.e., root mean square deviation values of 0.72-2.06 A). The homology-modeled CTP is consistent with our earlier de novo models of the transporter's transmembrane domains, with respect to residues which face into the transport path. Importantly, the resulting model is consistent with our previous experimental data obtained from measuring reactivity of 34 single cysteine mutants in transmembrane domains 3 and 4 with methanethiosulfonate reagents. The model also points to a likely dimer interface region. In conclusion, our data help to define the substrate translocation pathway in both the modeled CTP structure and the crystallographic ADP carrier structure.     

Trinodular structure of fibrinogen. Confirmation by both shadowing and negative stain electron microscopy

Fibrinogen molecules sprayed on mica, dried in vacuo and shadowed with platinum appear as trinodular rods. We describe here the flotation technique of negative staining by which we were able to visualize individual fibrinogen molecules. The fibrinogen molecules in our negatively stained preparations have a trinodular structure identical to that of shadowed molecules. We believe that the 20 nm diameter globules seen in previous negative staining studies are aggregates of these molecules. This is the first study in which both shadowed and negatively stained preparations of fibrinogen support a single, consistent model for the structure of fibrinogen. The molecules in our preparation are 45 nm long (±2.5 nm, our overall estimate of accuracy), the two outer nodules are 6 to 7 nm in diameter and the middle nodule is 4 to 5 nm in diameter.     

Assessing scales of variability in benthic diatom community structure

Sources of variability such as sampling method, sample preparation, and sample analysis (taxonomy) might affect our ability to measure differences in community structure between sites in environmental effects studies. It is therefore important to ensure that changes in community structure related to the physical or chemical differences between sites are not hidden by other sources of variability within a site. The goal of this study was to quantify the amount of variability in benthic diatom community structure related to sampling and laboratory procedures. Our results showed that variability in community structure was minimal among replicate microscope slides (< 1%) and among samples collected within a site (1.8%). Variability in community structure was substantially higher between sites located in one stream (16.6%), and highest across streams (59.6%). This suggests that field sampling and laboratory methods do not contribute a large amount of variation into our analyses of benthic diatom community structure across sites.     

Structural change in a B-DNA helix with hydrostatic pressure

Study of the effects of pressure on macromolecular structure improves our understanding of the forces governing structure, provides details on the relevance of cavities and packing in structure, increases our understanding of hydration and provides a basis to understand the biology of high-pressure organisms. A study of DNA, in particular, helps us to understand how pressure can affect gene activity. Here we present the first high-resolution experimental study of B-DNA structure at high pressure, using NMR data acquired at pressures up to 200 MPa (2 kbar). The structure of DNA compresses very little, but is distorted so as to widen the minor groove, and to compress hydrogen bonds, with AT pairs compressing more than GC pairs. The minor groove changes are suggested to lead to a compression of the hydration water in the minor groove.     

Local interactions drive the formation of nonnative structure in the denatured state of human alpha-lactalbumin: a high resolution structural characterization of a peptide model in aqueous solution.

There are a small number of peptides derived from proteins that have a propensity to adopt structure in aqueous solution which is similar to the structure they possess in the parent protein. There are far fewer examples of protein fragments which adopt stable nonnative structures in isolation. Understanding how nonnative interactions are involved in protein folding is crucial to our understanding of the topic. Here we show that a small, 11 amino acid peptide corresponding to residues 101-111 of the protein alpha-lactalbumin is remarkably structured in isolation in aqueous solution. The peptide has been characterized by 1H NMR, and 170 ROE-derived constraints were used to calculate a structure. The calculations yielded a single, high-resolution structure for residues 101-107 that is nonnative in both the backbone and side-chain conformations. In the pH 6.5 crystal structure, residues 101-105 are in an irregular turn-like conformation and residues 106-111 form an alpha-helix. In the pH 4.2 crystal structure, residues 101-105 form an alpha-helix, and residues 106-111 form a loopike structure. Both of these structures are significantly different from the conformation adopted by our peptide. The structure in the peptide model is primarily the result of local side-chain interactions that force the backbone to adopt a nonnative 310/turn-like structure in residues 103-106. The structure in aqueous solution was compared to the structure in 30% trifluoroethanol (TFE), and clear differences were observed. In particular, one of the side-chain interactions, a hydrophobic cluster involving residues 101-105, is different in the two solvents and residues 107-111 are considerably more ordered in 30% TFE. The implications of the nonnative structure for the folding of alpha-lactalbumin is discussed.     

Mutagenesis data in the automated prediction of transmembrane helix dimers

We present a molecular modeling protocol that selects modeled protein structures based on experimental mutagenesis results. The computed effect of a point mutation should be consistent with its experimental effect for correct models; mutations that do not affect protein stability and function should not affect the computed energy of a correct model while destabilizing mutations should have unfavorable computed energies. On the other hand, an incorrect model will likely display computed energies that are inconsistent with experimental results. We added terms to our energy function which penalize models that are inconsistent with experimental results. This creates a selective advantage for models that are consistent with experimental results in the Monte Carlo simulated annealing protocol we use to search conformational space. We calibrated our protocol to predict the structure of transmembrane helix dimers using glycophorin A as a model system. Inclusion of mutational data in this protocol compensates for the limitations of our force field and the limitations of our conformational search. We demonstrate an application of this structure prediction protocol by modeling the transmembrane region of the BNIP3 apoptosis factor.     

Bayesian association mapping of multiple quantitative trait loci and its application to the analysis of genetic variation among <Emphasis Type="Italic">Oryza sativa</Emphasis> L. germplasms

One way to use a crop germplasm collection directly to map QTLs without using line-crossing experiments is the whole genome association mapping. A major problem with association mapping is the presence of population structure, which can lead to both false positives and failure to detect genuine associations (i.e., false negatives). Particularly in highly selfing species such as Asian cultivated rice, high levels of population structure are expected and therefore the efficiency of association mapping remains almost unknown. Here, we propose an approach that combines a Bayesian method for mapping multiple QTLs with a regression method that directly incorporates estimates of population structure. That is, the effects due to both multiple QTLs and population structure were included in our statistical model. We evaluated the efficiency of our approach in simulated- and real-trait analyses of a rice germplasm collection. Simulation analyses based on real marker data showed that our model could suppress both false-positive and false-negative rates and the error of estimation of genetic effects over single QTL models, indicating that our model has statistically desirable attributes over single QTL models. As real traits, we analyzed the size and shape of milled rice grains and found significant markers that may be linked to QTLs reported previously. Association mapping should have good prospects in highly selfing species such as rice if proper methods are adopted. Our approach will be useful for the whole genome association mapping of various selfing crop species.     

Constructing templates for protein structure prediction by simulation of protein folding pathways

How a one-dimensional protein sequence folds into a specific 3D structure remains a difficult challenge in structural biology. Many computational methods have been developed in an attempt to predict the tertiary structure of the protein; most of these employ approaches that are based on the accumulated knowledge of solved protein structures. Here we introduce a novel and fully automated approach for predicting the 3D structure of a protein that is based on the well accepted notion that protein folding is a hierarchical process. Our algorithm follows the hierarchical model by employing two stages: the first aims to find a match between the sequences of short independently-folding structural entities and parts of the target sequence and assigns the respective structures. The second assembles these local structural parts into a complete 3D structure, allowing for long-range interactions between them. We present the results of applying our method to a subset of the targets from CASP6 and CASP7. Our results indicate that for targets with a significant sequence similarity to known structures we are often able to provide predictions that are better than those achieved by two leading servers, and that the most significant improvements in comparison with these methods occur in regions of a gapped structural alignment between the native structure and the closest available structural template. We conclude that in addition to performing well for targets with known homologous structures, our method shows great promise for addressing the more general category of comparative modeling targets, which is our next goal.     

Complex formation with kinesin motor domains affects the structure of microtubules

Microtubules are highly dynamic components of the cytoskeleton. They are important for cell movement and they are involved in a variety of transport processes together with motor proteins, such as kinesin. The exact mechanism of these transport processes is not known and so far the focus has been on structural changes within the motor domains, but not within the underlying microtubule structure.Here we investigated the interaction between kinesin and tubulin and our experimental data show that microtubules themselves are changing structure during that process. We studied unstained, vitrified samples of microtubules composed of 15 protofilaments using cryo electron microscopy and helical image analysis. 3D maps of plain microtubules and microtubules decorated with kinesin have been reconstructed to approximately 17A resolution. The alphabeta-tubulin dimer could be identified and, according to our data, alpha- and beta-tubulin adopt different conformations in plain microtubules. Significant differences were detected between maps of plain microtubules and microtubule-kinesin complexes. Most pronounced is the continuous axial inter-dimer contact in the microtubule-kinesin complex, suggesting stabilized protofilaments along the microtubule axis. It seems, that mainly structural changes within alpha-tubulin are responsible for this observation. Lateral effects are less pronounced. Following our data, we believe, that microtubules play an active role in intracellular transport processes through modulations of their core structure.     

Extracting information on folding from the amino acid sequence: consensus regions with preferred conformation in homologous proteins.

It is investigated whether protein segments predicted to have a well-defined conformational preference in the absence of tertiary interactions are conserved in families of homologous proteins. The prediction method follows the procedures of Rooman, M., Kocher, J.-P., and Wodak, S. (preceding paper in this issue). It uses a knowledge-based force field that incorporates only local interactions along the sequence and identifies segments whose lowest energy structure displays a sizable energy gap relative to other computed conformations. In 13 of the protein families and subfamilies considered that are sufficiently homologous to have similar 3D structures, at least one region is consistently predicted as having the same preferred conformation in virtually all family members. These regions are between 4 and 26 residues long. They are often located at chain ends and correspond primarily to segments of secondary structure heavily involved in interactions with the rest of the protein, suggesting that they could act as nuclei around which other parts of the structure would assemble. Experimental data on early folding intermediates or on protein fragments with appreciable structure in aqueous solution are available for more than half of the protein families. Comparison of our results with these data is quite favorable. They reveal that each of the experimentally identified early formed, or independently stable, substructures harbors at least one of the segments consistently predicted as having a preferred conformation by our procedure. The implications of our findings for the conservation of folding pathways in homologous proteins are discussed.     

浙江天童常绿阔叶林林冠结构与群落物种组成的关系

森林林冠结构能改变林下微气候条件, 可能会形成独立于地面生境的空间结构, 进而影响群落物种组成差异。该研究利用机载激光雷达获取浙江天童20 hm²常绿阔叶林样地的高精度林冠结构信息, 初步探讨了林冠结构与群落物种组成差异的关系, 结果表明: (1)未考虑林冠结构时, 独立于地面生境的空间结构是天童样地群落物种组成差异的重要影响因子, 在100 m²、400 m²、2 500 m²样方尺度上, 其对群落物种组成差异的解释率分别为25.2%、28.1%、8.0%。(2)考虑林冠结构后, 林冠结构使独立于地面生境的空间结构对群落物种组成差异的解释率降低了约1/3 (26.2%-36.0%)。(3)林冠结构因子中, 林冠高度对群落物种组成差异影响最大, 其次为林冠内部结构; 随样方尺度增大, 林冠高度对群落物种组成差异的影响降低, 林冠内部结构的影响逐渐增加。该研究结果证明了林冠结构是独立于地面生境的空间结构的主要驱动因子, 对天童植物群落物种组成差异具有不可忽视的重要作用。这些结果明晰了林冠结构因子中林冠高度和内部结构的重要性。     

Improving a consensus approach for protein structure selection by removing redundancy

In protein tertiary structure prediction, a crucial step is to select near-native structures from a large number of predicted structural models. Over the years, extensive research has been conducted for the protein structure selection problem with most approaches focusing on developing more accurate energy or scoring functions. Despite significant advances in this area, the discerning power of current approaches is still unsatisfactory. In this paper, we propose a novel consensus-based algorithm for the selection of predicted protein structures. Given a set of predicted models, our method first removes redundant structures to derive a subset of reference models. Then, a structure is ranked based on its average pairwise similarity to the reference models. Using the CASP8 data set containing a large collection of predicted models for 122 targets, we compared our method with the best CASP8 quality assessment (QA) servers, which are all consensus based, and showed that our QA scores correlate better with the GDT-TSs than those of the CASP8 QA servers. We also compared our method with the state-of-the-art scoring functions and showed its improved performance for near-native model selection. The GDT-TSs of the top models picked by our method are on average more than 8 percent better than the ones selected by the best performing scoring function.