Coarse‐ and fine‐grained models for proteins: Evaluation by decoy discrimination

Coarse‐grained models for protein structure are increasingly used in simulations and structural bioinformatics. In this study, we evaluated the effectiveness of three granularities of protein representation based on their ability to discriminate between correctly folded native structures and incorrectly folded decoy structures. The three levels of representation used one bead per amino acid (coarse), two beads per amino acid (medium), and all atoms (fine). Multiple structure features were compared at each representation level including two‐body interactions, three‐body interactions, solvent exposure, contact numbers, and angle bending. In most cases, the all‐atom level was most successful at discriminating decoys, but the two‐bead level provided a good compromise between the number of model parameters which must be estimated and the accuracy achieved. The most effective feature type appeared to be two‐body interactions. Considering three‐body interactions increased accuracy only marginally when all atoms were used and not at all in medium and coarse representations. Though two‐body interactions were most effective for the coarse representations, the accuracy loss for using only solvent exposure or contact number was proportionally less at these levels than in the all‐atom representation. We propose an optimization method capable of selecting bead types of different granularities to create a mixed representation of the protein. We illustrate its behavior on decoy discrimination and discuss implications for data‐driven protein model selection. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

The evolution of bat nucleic acid‐sensing Toll‐like receptors

We characterized the nucleic acid‐sensing Toll‐like receptors (TLR) of a New World bat species, the common vampire bat (Desmodus rotundus), and through a comparative molecular evolutionary approach searched for general adaptation patterns among the nucleic acid‐sensing TLRs of eight different bats species belonging to three families (Pteropodidae, Vespertilionidae and Phyllostomidae). We found that the bat TLRs are evolving slowly and mostly under purifying selection and that the divergence pattern of such receptors is overall congruent with the species tree, consistent with the evolution of many other mammalian nuclear genes. However, the chiropteran TLRs exhibited unique mutations fixed in ligand‐binding sites, some of which involved nonconservative amino acid changes and/or targets of positive selection. Such changes could potentially modify protein function and ligand‐binding properties, as some changes were predicted to alter nucleic acid binding motifs in TLR 9. Moreover, evidence for episodic diversifying selection acting specifically upon the bat lineage and sublineages was detected. Thus, the long‐term adaptation of chiropterans to a wide variety of environments and ecological niches with different pathogen profiles is likely to have shaped the evolution of the bat TLRs in an order‐specific manner. The observed evolutionary patterns provide evidence for potential functional differences between bat and other mammalian TLRs in terms of resistance to specific pathogens or recognition of nucleic acids in general.     

Low‐coherent optical diffraction tomography by angle‐scanning illumination

Temporally low‐coherent optical diffraction tomography (ODT) is proposed and demonstrated based on angle‐scanning Mach‐Zehnder interferometry. Using a digital micromirror device based on diffractive tilting, the full‐field interference of incoherent light is successfully maintained during every angle‐scanning sequences. Further, current ODT reconstruction principles for temporally incoherent illuminations are thoroughly reviewed and developed. Several limitations of incoherent illumination are also discussed, such as the nondispersive assumption, optical sectioning capacity and illumination angle limitation. Using the proposed setup and reconstruction algorithms, low‐coherent ODT imaging of plastic microspheres, human red blood cells and rat pheochromocytoma cells is experimentally demonstrated.      

Studies of the B‐Z transition of DNA: The temperature dependence of the free‐energy difference,the composition of the counterion sheath in mixed salt,and the preparation of a sample of the 5′‐d[T‐(m5C‐G)12‐T] duplex in pure B‐DNA or Z‐DNA form

It is often envisioned that cations might coordinate at specific sites of nucleic acids and play an important structural role, for instance in the transition between B‐DNA and Z‐DNA. However, nucleic acid models explicitly devoid of specific sites may also exhibit features previously considered as evidence for specific binding. Such is the case of the “composite cylinder” (or CC) model which spreads out localized features of DNA structure and charge by cylindrical averaging, while sustaining the main difference between the B and Z structures, namely the better immersion of the B‐DNA phosphodiester charges in the solution. Here, we analyze the non‐electrostatic component of the free‐energy difference between B‐DNA and Z‐DNA. We also compute the composition of the counterion sheath in a wide range of mixed‐salt solutions and of temperatures: in contrast with the large difference of composition between the B‐DNA and Z‐DNA forms, the temperature dependence of sheath composition, previously unknown, is very weak. In order to validate the model, the mixed‐salt predictions should be compared to experiment. We design a procedure for future measurements of the sheath composition based on Anomalous Small‐Angle X‐ray Scattering and complemented by ³¹P NMR. With due consideration for the kinetics of the B‐Z transition and for the capacity of generating at will the B or Z form in a single sample, the 5′‐d[T‐(m⁵C‐G)₁₂‐T] 26‐mer emerges as a most suitable oligonucleotide for this study. Finally, the application of the finite element method to the resolution of the Poisson‐Boltzmann equation is described in detail. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 369–384, 2016.     

Label‐free indicator‐free nucleic acid biosensors using carbon nanotubes

Carbon nanotubes (CNTs) are promising components for electrical biosensors due to their high surface‐to‐volume ratio and improved electron transfer properties. This review surveys CNT‐based label‐free indicator‐free biosensing strategies that have been demonstrated for the sensitive detection of nucleic acids. After an introduction to CNTs, the fabrication of biosensors and techniques for the immobilization of probe nucleic acids are outlined. Subsequently, two major label‐free strategies namely electrochemical transduction and field‐effect detection are presented. The focus is on direct detection methods that avoid labels, indicators, intercalating agents, mediators, and even secondary receptors. The review concludes with a comparison between the various biosensors and presents ways of engineering them so that they can be deployed in realistic diagnostic applications.     

Selenium Derivatization of Nucleic Acids for X‐Ray Crystal‐Structure and Function Studies

It is estimated that over two thirds of all new crystal structures of proteins are determined via the protein selenium derivatization (selenomethionine (Se‐Met) strategy). This selenium derivatization strategy via MAD (multi‐wavelength anomalous dispersion) phasing has revolutionized protein X‐ray crystallography. Through our pioneer research, similarly, Se has also been successfully incorporated into nucleic acids to facilitate the X‐ray crystal‐structure and function studies of nucleic acids. Currently, Se has been stably introduced into nucleic acids by replacing nucleotide O‐atom at the positions 2′, 4′, 5′, and in nucleobases and non‐bridging phosphates. The Se derivatization of nucleic acids can be achieved through solid‐phase chemical synthesis and enzymatic methods, and the Se‐derivatized nucleic acids (SeNA) can be easily purified by HPLC, FPLC, and gel electrophoresis to obtain high purity. It has also been demonstrated that the Se derivatization of nucleic acids facilitates the phase determination via MAD phasing without significant perturbation. A growing number of structures of DNAs, RNAs, and protein–nucleic acid complexes have been determined by the Se derivatization and MAD phasing. Furthermore, it was observed that the Se derivatization can facilitate crystallization, especially when it is introduced to the 2′‐position. In addition, this novel derivatization strategy has many advantages over the conventional halogen derivatization, such as more choices of the modification sites via the atom‐specific substitution of the nucleotide O‐atom, better stability under X‐ray radiation, and structure isomorphism. Therefore, our Se‐derivatization strategy has great potentials to provide rational solutions for both phase determination and high‐quality crystal growth in nucleic‐acid crystallography. Moreover, the Se derivatization generates the nucleic acids with many new properties and creates a new paradigm of nucleic acids. This review summarizes the recent developments of the atomic site‐specific Se derivatization of nucleic acids for structure determination and function study. Several applications of this Se‐derivatization strategy in nucleic acid and protein research are also described in this review.     

Folding and stability of helical bundle proteins from coarse‐grained models

We develop a coarse‐grained model where solvent is considered implicitly, electrostatics are included as short‐range interactions, and side‐chains are coarse‐grained to a single bead. The model depends on three main parameters: hydrophobic, electrostatic, and side‐chain hydrogen bond strength. The parameters are determined by considering three level of approximations and characterizing the folding for three selected proteins (training set). Nine additional proteins (containing up to 126 residues) as well as mutated versions (test set) are folded with the given parameters. In all folding simulations, the initial state is a random coil configuration. Besides the native state, some proteins fold into an additional state differing in the topology (structure of the helical bundle). We discuss the stability of the native states, and compare the dynamics of our model to all atom molecular dynamics simulations as well as some general properties on the interactions governing folding dynamics. Proteins 2013; 81:1200–1211. © 2013 Wiley Periodicals, Inc.     

Covid‐19.bioreproducibility.org: A web resource for SARS‐CoV‐2‐related structural models

The COVID‐19 pandemic has triggered numerous scientific activities aimed at understanding the SARS‐CoV‐2 virus and ultimately developing treatments. Structural biologists have already determined hundreds of experimental X‐ray, cryo‐EM, and NMR structures of proteins and nucleic acids related to this coronavirus, and this number is still growing. To help biomedical researchers, who may not necessarily be experts in structural biology, navigate through the flood of structural models, we have created an online resource, covid19.bioreproducibility.org, that aggregates expert‐verified information about SARS‐CoV‐2‐related macromolecular models. In this article, we describe this web resource along with the suite of tools and methodologies used for assessing the structures presented therein.     

Breast Cancer‐Derived Exosomes Reflect the Cell‐of‐Origin Phenotype

A manner in which cells can communicate with each other is via secreted nanoparticles termed exosomes. These vesicles contain lipids, nucleic acids, and proteins, and are said to reflect the cell‐of‐origin. However, for the exosomal protein content, there is limited evidence in the literature to verify this statement. Here, proteomic assessment combined with pathway‐enrichment analysis is used to demonstrate that the protein cargo of exosomes reflects the epithelial/mesenchymal phenotype of secreting breast cancer cells. Given that epithelial‐mesenchymal plasticity is known to implicate various stages of cancer progression, the results suggest that breast cancer subtypes with distinct epithelial and mesenchymal phenotypes may be distinguished by directly assessing the protein content of exosomes. Additionally, the work is a substantial step toward verifying the statement that cell‐derived exosomes reflect the phenotype of the cells‐of‐origin.     

SARS‐CoV‐2 lgM/lgG antibody detection confirms the infection after three negative nucleic acid detection

An ongoing outbreak of viral pneumonia was caused by a novel coronavirus in China in 2019. By March 19, over 200 thousand confirmed cases of SARS‐CoV‐2 infection and over 9000 deaths have been reported throughout the world. For this infectious disease, nucleic acid detection is still the gold standard for pathogenic detection. However, nucleic acid detection takes a long time and has relatively high "false negative"; therefore, we need urgently a convenient and accurate detection method to make up for this deficiency. In this article, we will show such technical characteristics of lgM/lgG serum antibody detection, compared with nucleic acid detection.     

Spatio‐temporal Modelling of Weeds by Shot‐noise G Cox processes

Tractable space‐time point processes models are needed in various fields. For example in weed science for gaining biological knowledge, for prediction of weed development in order to optimize local treatments with herbicides or in epidemiology for prediction of the risk of a disease. Motivated by the spatio‐temporal point patterns for two weed species, we propose a spatio‐temporal Cox model with intensity based on gamma random fields. The model is an extension of Neyman–Scott and shot‐noise Cox processes to the space‐time domain and it allows spatial and temporal inhomogeneity. We use the weed example to give a first intuitive interpretation of the model and then show how the model is constructed more rigorously and how to estimate the parameters. The weed data are analysed using the proposed model, and both spatially and temporally the model shows a good fit to the data using classical goodness‐of‐fit tests.     

Selection of bead‐displayed,PNA‐encoded chemicals

The lack of efficient identification and isolation methods for specific molecular binders has fundamentally limited drug discovery. Here, we have developed a method to select peptide nucleic acid (PNA) encoded molecules with specific functional properties from combinatorially generated libraries. This method consists of three essential stages: (1) creation of a Lab‐on‐Bead? library, a one‐bead, one‐sequence library that, in turn, displays a library of candidate molecules, (2) fluorescence microscopy‐aided identification of single target‐bound beads and the extraction – wet or dry – of these beads and their attached candidate molecules by a micropipette manipulator, and (3) identification of the target‐binding candidate molecules via amplification and sequencing. This novel integration of techniques harnesses the sensitivity of DNA detection methods and the multiplexed and miniaturized nature of molecule screening to efficiently select and identify target‐binding molecules from large nucleic acid encoded chemical libraries. Beyond its potential to accelerate assays currently used for the discovery of new drug candidates, its simple bead‐based design allows for easy screening over a variety of prepared surfaces that can extend this technique's application to the discovery of diagnostic reagents and disease markers. Copyright © 2009 John Wiley & Sons, Ltd.     

Heterojunction‐Depleted Lead‐Free Perovskite Solar Cells with Coarse‐Grained B‐γ‐CsSnI3 Thin Films

Perovskite solar cells (PSCs) have been emerging as a breakthrough photovoltaic technology, holding unprecedented promise for low‐cost, high‐efficiency renewable electricity generation. However, potential toxicity associated with the state‐of‐the‐art lead‐containing PSCs has become a major concern. The past research in the development of lead‐free PSCs has met with mixed success. Herein, the promise of coarse‐grained B‐γ‐CsSnI₃ perovskite thin films as light absorber for efficient lead‐free PSCs is demonstrated. Thermally‐driven solid‐state coarsening of B‐γ‐CsSnI₃ perovskite grains employed here is accompanied by an increase of tin‐vacancy concentration in their crystal structure, as supported by first‐principles calculations. The optimal device architecture for the efficient photovoltaic operation of these B‐γ‐CsSnI₃ thin films is identified through exploration of several device architectures. Via modulation of the B‐γ‐CsSnI₃ grain coarsening, together with the use of the optimal PSC architecture, planar heterojunction‐depleted B‐γ‐CsSnI₃ PSCs with power conversion efficiency up to 3.31% are achieved without the use of any additives. The demonstrated strategies provide guidelines and prospects for developing future high‐performance lead‐free PVs.     

Disordered nucleiome: Abundance of intrinsic disorder in the DNA‐ and RNA‐binding proteins in 1121 species from Eukaryota,Bacteria and Archaea

Intrinsically disordered proteins (IDPs) are abundant in various proteomes, where they play numerous important roles and complement biological activities of ordered proteins. Among functions assigned to IDPs are interactions with nucleic acids. However, often, such assignments are made based on the guilty‐by‐association principle. The validity of the extension of these correlations to all nucleic acid binding proteins has never been analyzed on a large scale across all domains of life. To fill this gap, we perform a comprehensive computational analysis of the abundance of intrinsic disorder and intrinsically disordered domains in nucleiomes (～548 000 nucleic acid binding proteins) of 1121 species from Archaea, Bacteria and Eukaryota. Nucleiome is a whole complement of proteins involved in interactions with nucleic acids. We show that relative to other proteins in the corresponding proteomes, the DNA‐binding proteins have significantly increased disorder content and are significantly enriched in disordered domains in Eukaryotes but not in Archaea and Bacteria. The RNA‐binding proteins are significantly enriched in the disordered domains in Bacteria, Archaea and Eukaryota, while the overall abundance of disorder in these proteins is significantly increased in Bacteria, Archaea, animals and fungi. The high abundance of disorder in nucleiomes supports the notion that the nucleic acid binding proteins often require intrinsic disorder for their functions and regulation.     

All‐atom and coarse‐grained simulations of the forced unfolding pathways of the SNARE complex

The SNARE complex, consisting of three proteins (VAMP2, syntaxin, and SNAP‐25), is thought to drive membrane fusion by assembling into a four‐helix bundle through a zippering process. In support of the above zippering model, a recent single‐molecule optical tweezers experiment by Gao et al. revealed a sequential unzipping of SNARE along VAMP2 in the order of the linker domain → the C‐terminal domain → the N‐terminal domain. To offer detailed structural insights to this unzipping process, we have performed all‐atom and coarse‐grained steered molecular dynamics (sMD) simulations of the forced unfolding pathways of SNARE using different models and force fields. Our findings are summarized as follows: First, the sMD simulations based on either an all‐atom force field (with an implicit solvent model) or a coarse‐grained Go model were unable to capture the forced unfolding pathway of SNARE as observed by Gao et al., which may be attributed to insufficient simulation time and inaccurate force fields. Second, the sMD simulations based on a reparameterized coarse‐grained model (i.e., modified elastic network model) were able to predict a sequential unzipping of SNARE in good agreement with the findings by Gao et al. The key to this success is to reparameterize the intrahelix and interhelix nonbonded force constants against the pair‐wise residue–residue distance fluctuations collected from all‐atom MD simulations of SNARE. Therefore, our finding supports the importance of accurately describing the inherent dynamics/flexibility of SNARE (in the absence of force), in order to correctly simulate its unfolding behaviors under force. This study has established a useful computational framework for future studies of the zippering function of SNARE and its perturbations by point mutations with amino‐acid level of details, and more generally the forced unfolding pathways of other helix bundle proteins. Proteins 2014; 82:1376–1386. © 2014 Wiley Periodicals, Inc.     

A hybrid clustering algorithm for multiple‐source resolving in bioluminescence tomography

Bioluminescence tomography is a preclinical imaging modality to locate and quantify internal bioluminescent sources from surface measurements, which experienced rapid growth in the last 10 years. However, multiple‐source resolving remains a challenging issue in BLT. In this study, it is treated as an unsupervised pattern recognition problem based on the reconstruction result, and a novel hybrid clustering algorithm combining the advantages of affinity propagation (AP) and K‐means is developed to identify multiple sources automatically. Moreover, we incorporate the clustering analysis into a general multiple‐source reconstruction framework, which can provide stable reconstruction and accurate resolving result without providing the number of targets. Numerical simulations and in vivo experiments on 4T1‐luc2 mouse model were conducted to assess the performance of the proposed method in multiple‐source resolving. The encouraging results demonstrate significant effectiveness and potential of our method in preclinical BLT applications.      

Binding Pockets and Permeation Channels for Dioxygen through Cofactorless 3‐Hydroxy‐2‐methylquinolin‐4‐one 2,4‐Dioxygenase in Association with Its Natural Substrate, 3‐Hydroxy‐2‐methylquinolin‐4(1H)‐one. A Perspective from Molecular Dynamics Simulations

This work describes an investigation of pathways and binging pockets (BPs) for dioxygen (O₂) through the cofactorless oxygenase 3‐hydroxy‐2‐methylquinolin‐4‐one 2,4‐dioxygenase in complex with its natural substrate, 3‐hydroxy‐2‐methylquinolin‐4(1H)‐one, in aqueous solution. The investigation tool was random‐acceleration molecular dynamics (RAMD), whereby a tiny, randomly oriented external force is applied to O₂ in order to accelerate its movements. In doing that, care was taken that the external force only continues, if O₂ moves along a direction for a given period of time, otherwise the force changed direction randomly. Gates for expulsion of O₂ from the protein, which can also be taken as gates for O₂ uptake, were found throughout almost the whole external surface of the protein, alongside a variety of BPs for O₂. The most exploited gates and BPs were not found to correspond to the single gate and BP proposed previously from the examination of the static model from X‐ray diffraction analysis of this system. Therefore, experimental investigations of this system that go beyond the static model are urgently needed.