Isotopic Depletion Increases the Spatial Resolution of FPOP Top-Down Mass Spectrometry Analysis

Protein radical labeling, like fast photochemical oxidation of proteins (FPOP), coupled to a top-down mass spectrometry (MS) analysis offers an alternative analytical method for probing protein structure or protein interaction with other biomolecules, for instance, proteins and DNA. However, with the increasing mass of studied analytes, the MS/MS spectra become complex and exhibit a low signal-to-noise ratio. Nevertheless, these difficulties may be overcome by protein isotope depletion. Thus, we aimed to use protein isotope depletion to analyze FPOP-oxidized samples by top-down MS analysis. For this purpose, we prepared isotopically natural (IN) and depleted (ID) forms of the FOXO4 DNA binding domain (FOXO4-DBD) and studied the protein-DNA interaction interface with double-stranded DNA, the insulin response element (IRE), after exposing the complex to hydroxyl radicals. As shown by comparing tandem mass spectra of natural and depleted proteins, the ID form increased the signal-to-noise ratio of useful fragment ions, thereby enhancing the sequence coverage by more than 19%. This improvement in the detection of fragment ions enabled us to detect 22 more oxidized residues in the ID samples than in the IN sample. Moreover, less common modifications were detected in the ID sample, including the formation of ketones and lysine carbonylation. Given the higher quality of ID top-down MSMS data set, these results provide more detailed information on the complex formation between transcription factors and DNA-response elements. Therefore, our study highlights the benefits of isotopic depletion for quantitative top-down proteomics. Data are available via ProteomeXchange with the identifier PXD044447.

One bead per residue can describe all-atom protein structures

Atomistic resolution is the standard for high-resolution biomolecular structures, but experimental structural data are often at lower resolution. Coarse-grained models are also used extensively in computational studies to reach biologically relevant spatial and temporal scales. This study explores the use of advanced machine learning networks for reconstructing atomistic models from reduced representations. The main finding is that a single bead per amino acid residue allows construction of accurate and stereochemically realistic all-atom structures with minimal loss of information. This suggests that lower resolution representations of proteins may be sufficient for many applications when combined with a machine learning framework that encodes knowledge from known structures. Practical applications include the rapid addition of atomistic detail to low-resolution structures from experiment or computational coarse-grained models. The application of rapid, deterministic all-atom reconstruction within multi-scale frameworks is further demonstrated with a rapid protocol for the generation of accurate models from cryo-EM densities close to experimental structures.

RNAStat: An Integrated Tool for Statistical Analysis of RNA 3D Structures

The 3D architectures of RNAs are essential for understanding their cellular functions. While an accurate scoring function based on the statistics of known RNA structures is a key component for successful RNA structure prediction or evaluation, there are few tools or web servers that can be directly used to make comprehensive statistical analysis for RNA 3D structures. In this work, we developed RNAStat, an integrated tool for making statistics on RNA 3D structures. For given RNA structures, RNAStat automatically calculates RNA structural properties such as size and shape, and shows their distributions. Based on the RNA structure annotation from DSSR, RNAStat provides statistical information of RNA secondary structure motifs including canonical/non-canonical base pairs, stems, and various loops. In particular, the geometry of base-pairing/stacking can be calculated in RNAStat by constructing a local coordinate system for each base. In addition, RNAStat also supplies the distribution of distance between any atoms to the users to help build distance-based RNA statistical potentials. To test the usability of the tool, we established a non-redundant RNA 3D structure dataset, and based on the dataset, we made a comprehensive statistical analysis on RNA structures, which could have the guiding significance for RNA structure modeling. The python code of RNAStat, the dataset used in this work, and corresponding statistical data files are freely available at GitHub (https://github.com/RNA-folding-lab/RNAStat).

Structural dynamics and allostery of Rab proteins: strategies for drug discovery and design

Rab proteins represent the largest family of the Rab superfamily guanosine triphosphatase (GTPase). Aberrant human Rab proteins are associated with multiple diseases, including cancers and neurological disorders. Rab subfamily members display subtle conformational variations that render specificity in their physiological functions and can be targeted for subfamily-specific drug design. However, drug discovery efforts have not focused much on targeting Rab allosteric non-nucleotide binding sites which are subjected to less evolutionary pressures to be conserved, hence are likely to offer subfamily specificity and may be less prone to undesirable off-target interactions and side effects. To discover druggable allosteric binding sites, Rab structural dynamics need to be first incorporated using multiple experimentally and computationally obtained structures. The high-dimensional structural data may necessitate feature extraction methods to identify manageable representative structures for subsequent analyses. We have detailed state-of-the-art computational methods to (i) identify binding sites using data on sequence, shape, energy, etc., (ii) determine the allosteric nature of these binding sites based on structural ensembles, residue networks and correlated motions and (iii) identify small molecule binders through structure- and ligand-based virtual screening. To benefit future studies for targeting Rab allosteric sites, we herein detail a refined workflow comprising multiple available computational methods, which have been successfully used alone or in combinations. This workflow is also applicable for drug discovery efforts targeting other medically important proteins. Depending on the structural dynamics of proteins of interest, researchers can select suitable strategies for allosteric drug discovery and design, from the resources of computational methods and tools enlisted in the workflow.

Molecular Dynamics Simulations Suggest a Non-Doublet Decoding Model of -1 Frameshifting by tRNASer3

In-frame decoding in the ribosome occurs through canonical or wobble Watson-Crick pairing of three mRNA codon bases (a triplet) with a triplet of anticodon bases in tRNA. Departures from the triplet-triplet interaction can result in frameshifting, meaning downstream mRNA codons are then read in a different register. There are many mechanisms to induce frameshifting, and most are insufficiently understood. One previously proposed mechanism is doublet decoding, in which only codon bases 1 and 2 are read by anticodon bases 34 and 35, which would lead to -1 frameshifting. In E. coli, tRNA^Ser3_GCU can induce -1 frameshifting at alanine (GCA) codons. The logic of the doublet decoding model is that the Ala codon's GC could pair with the tRNA^Ser3's GC, leaving the third anticodon residue U36 making no interactions with mRNA. Under that model, a U36C mutation would still induce -1 frameshifting, but experiments refute this. We perform all-atom simulations of wild-type tRNA^Ser3, as well as a U36C mutant. Our simulations revealed a hydrogen bond between U36 of the anticodon and G1 of the codon. The U36C mutant cannot make this interaction, as it lacks the hydrogen-bond-donating H3. The simulation thus suggests a novel, non-doublet decoding mechanism for -1 frameshifting by tRNA^Ser3 at Ala codons.

Elastic continuum stiffness of contractile tail sheaths from molecular dynamics simulations

Contractile tails are key components of the biological nanomachinery involved in cell membrane puncturing, where they provide a means to deliver molecules and ions inside cells. Two intriguing examples of contractile tails are those from bacteriophage T4 and R2-pyocin. Although the two systems are different in terms of biological activity, they share a fascinatingly similar injection mechanism, during which the tail sheaths of both systems contract from a so-called extended state to around half of their length (the contracted state), accompanied by release of elastic energy originally stored in the sheath. Despite the great prevalence and biomedical importance of contractile delivery systems, many fundamental details of their injection machinery and dynamics are still unknown. In this work, we calculate the bending and torsional stiffness constants of a helical tail sheath strand of bacteriophage T4 and R2-pyocin, in both extended and contracted states, using molecular dynamics simulations of about one-sixth of the entire sheath. Differences in stiffness constants between the two systems are rationalized by comparing their all-atom monomer structures, changes in sheath architecture on contraction, and differences in interstrand interactions. The calculated coefficients indicate that the T4 strand is stiffer for both bending and torsion than the corresponding R2-pyocin strands in both extended and contracted conformations. The sheath strands also have greater stiffness in the contracted state for both systems. As the main application of this study, we describe how the stiffness constants can be incorporated in a model to simulate the dynamics of contractile nanoinjection machineries.

Bioinformatics Tools and Benchmarks for Computational Docking and 3D Structure Prediction of RNA-Protein Complexes

RNA-protein (RNP) interactions play essential roles in many biological processes, such as regulation of co-transcriptional and post-transcriptional gene expression, RNA splicing, transport, storage and stabilization, as well as protein synthesis. An increasing number of RNP structures would aid in a better understanding of these processes. However, due to the technical difficulties associated with experimental determination of macromolecular structures by high-resolution methods, studies on RNP recognition and complex formation present significant challenges. As an alternative, computational prediction of RNP interactions can be carried out. Structural models obtained by theoretical predictive methods are, in general, less reliable compared to models based on experimental measurements but they can be sufficiently accurate to be used as a basis for to formulating functional hypotheses. In this article, we present an overview of computational methods for 3D structure prediction of RNP complexes. We discuss currently available methods for macromolecular docking and for scoring 3D structural models of RNP complexes in particular. Additionally, we also review benchmarks that have been developed to assess the accuracy of these methods.

Investigating Substitutions in Antibody-Antigen Complexes Using Molecular Dynamics: A Case Study with Broad-spectrum, Influenza A Antibodies

In studying the binding of host antibodies to the surface antigens of pathogens, the structural and functional characterization of antibody–antigen complexes by X-ray crystallography and binding assay is important. However, the characterization requires experiments that are typically time consuming and expensive: thus, many antibody–antigen complexes are under-characterized. For vaccine development and disease surveillance, it is often vital to assess the impact of amino acid substitutions on antibody binding. For example, are there antibody substitutions capable of improving binding without a loss of breadth, or antigen substitutions that lead to antigenic escape? The questions cannot be answered reliably from sequence variation alone, exhaustive substitution assays are usually impractical, and alanine scans provide at best an incomplete identification of the critical residue–residue interactions. Here, we show that, given an initial structure of an antibody bound to an antigen, molecular dynamics simulations using the energy method molecular mechanics with Generalized Born surface area (MM/GBSA) can model the impact of single amino acid substitutions on antibody–antigen binding energy. We apply the technique to three broad-spectrum antibodies to influenza A hemagglutinin and examine both previously characterized and novel variant strains observed in the human population that may give rise to antigenic escape. We find that in some cases the impact of a substitution is local, while in others it causes a reorientation of the antibody with wide-ranging impact on residue–residue interactions: this explains, in part, why the change in chemical properties of a residue can be, on its own, a poor predictor of overall change in binding energy. Our estimates are in good agreement with experimental results—indeed, they approximate the degree of agreement between different experimental techniques. Simulations were performed on commodity computer hardware; hence, this approach has the potential to be widely adopted by those undertaking infectious disease research. Novel aspects of this research include the application of MM/GBSA to investigate binding between broadly binding antibodies and a viral glycoprotein; the development of an approach for visualizing substrate–ligand interactions; and the use of experimental assay data to rescale our predictions, allowing us to make inferences about absolute, as well as relative, changes in binding energy.

Advances in Imaging Techniques and Genetically Encoded Probes for Photoacoustic Imaging

Photoacoustic (PA) imaging is a rapidly emerging biomedical imaging modality that is capable of visualizing cellular and molecular functions with high detection sensitivity and spatial resolution in deep tissue. Great efforts and progress have been made on the development of various PA imaging technologies with improved resolution and sensitivity over the past two decades. Various PA probes with high contrast have also been extensively developed, with many important biomedical applications. In comparison with chemical dyes and nanoparticles, genetically encoded probes offer easier labeling of defined cells within tissues or proteins of interest within a cell, have higher stability in vivo, and eliminate the need for delivery of exogenous substances. Genetically encoded probes have thus attracted increasing attention from researchers in engineering and biomedicine. In this review, we aim to provide an overview of the existing PA imaging technologies and genetically encoded PA probes, and describe further improvements in PA imaging techniques and the near-infrared photochromic protein BphP1, the most sensitive genetically encoded probe thus far, as well as the potential biomedical applications of BphP1-based PA imaging in vivo.

MMB-GUI: a fast morphing method demonstrates a possible ribosomal tRNA translocation trajectory

Easy-to-use macromolecular viewers, such as UCSF Chimera, are a standard tool in structural biology. They allow rendering and performing geometric operations on large complexes, such as viruses and ribosomes. Dynamical simulation codes enable modeling of conformational changes, but may require considerable time and many CPUs. There is an unmet demand from structural and molecular biologists for software in the middle ground, which would allow visualization combined with quick and interactive modeling of conformational changes, even of large complexes. This motivates MMB-GUI. MMB uses an internal-coordinate, multiscale approach, yielding as much as a 2000-fold speedup over conventional simulation methods. We use Chimera as an interactive graphical interface to control MMB. We show how this can be used for morphing of macromolecules that can be heterogeneous in biopolymer type, sequence, and chain count, accurately recapitulating structural intermediates. We use MMB-GUI to create a possible trajectory of EF-G mediated gate-passing translocation in the ribosome, with all-atom structures. This shows that the GUI makes modeling of large macromolecules accessible to a wide audience. The morph highlights similarities in tRNA conformational changes as tRNA translocates from A to P and from P to E sites and suggests that tRNA flexibility is critical for translocation completion.

GARN: Sampling RNA 3D Structure Space with Game Theory and Knowledge-Based Scoring Strategies

Cellular processes involve large numbers of RNA molecules. The functions of these RNA molecules and their binding to molecular machines are highly dependent on their 3D structures. One of the key challenges in RNA structure prediction and modeling is predicting the spatial arrangement of the various structural elements of RNA. As RNA folding is generally hierarchical, methods involving coarse-grained models hold great promise for this purpose. We present here a novel coarse-grained method for sampling, based on game theory and knowledge-based potentials. This strategy, GARN (Game Algorithm for RNa sampling), is often much faster than previously described techniques and generates large sets of solutions closely resembling the native structure. GARN is thus a suitable starting point for the molecular modeling of large RNAs, particularly those with experimental constraints. GARN is available from: http://garn.lri.fr/.

Inferring the microscopic surface energy of protein-protein interfaces from mutation data

Mutations at protein-protein recognition sites alter binding strength by altering the chemical nature of the interacting surfaces. We present a simple surface energy model, parameterized with empirical ΔΔG values, yielding mean energies of -48 cal mol(-1) Å(-2) for interactions between hydrophobic surfaces, -51 to -80 cal mol(-1) Å(-2) for surfaces of complementary charge, and 66-83 cal mol(-1) Å(-2) for electrostatically repelling surfaces, relative to the aqueous phase. This places the mean energy of hydrophobic surface burial at -24 cal mol(-1) Å(-2) . Despite neglecting configurational entropy and intramolecular changes, the model correlates with empirical binding free energies of a functionally diverse set of rigid-body interactions (r = 0.66). When used to rerank docking poses, it can place near-native solutions in the top 10 for 37% of the complexes evaluated, and 82% in the top 100. The method shows that hydrophobic burial is the driving force for protein association, accounting for 50-95% of the cohesive energy. The model is available open-source from http://life.bsc.es/pid/web/surface_energy/ and via the CCharpPPI web server http://life.bsc.es/pid/ccharppi/.

Protein-RNA complexes and efficient automatic docking: expanding RosettaDock possibilities

Protein-RNA complexes provide a wide range of essential functions in the cell. Their atomic experimental structure solving, despite essential to the understanding of these functions, is often difficult and expensive. Docking approaches that have been developed for proteins are often challenging to adapt for RNA because of its inherent flexibility and the structural data available being relatively scarce. In this study we adapted the RosettaDock protocol for protein-RNA complexes both at the nucleotide and atomic levels. Using a genetic algorithm-based strategy, and a non-redundant protein-RNA dataset, we derived a RosettaDock scoring scheme able not only to discriminate but also score efficiently docking decoys. The approach proved to be both efficient and robust for generating and identifying suitable structures when applied to two protein-RNA docking benchmarks in both bound and unbound settings. It also compares well to existing strategies. This is the first approach that currently offers a multi-level optimized scoring approach integrated in a full docking suite, leading the way to adaptive fully flexible strategies.

LipidWrapper: an algorithm for generating large-scale membrane models of arbitrary geometry

As ever larger and more complex biological systems are modeled in silico, approximating physiological lipid bilayers with simple planar models becomes increasingly unrealistic. In order to build accurate large-scale models of subcellular environments, models of lipid membranes with carefully considered, biologically relevant curvature will be essential. In the current work, we present a multi-scale utility called LipidWrapper capable of creating curved membrane models with geometries derived from various sources, both experimental and theoretical. To demonstrate its utility, we use LipidWrapper to examine an important mechanism of influenza virulence. A copy of the program can be downloaded free of charge under the terms of the open-source FreeBSD License from http://nbcr.ucsd.edu/lipidwrapper. LipidWrapper has been tested on all major computer operating systems.

A multiscale approach to predicting affinity changes in protein-protein interfaces

Substitution mutations in protein-protein interfaces can have a substantial effect on binding, which has consequences in basic and applied biomedical research. Experimental expression, purification, and affinity determination of protein complexes is an expensive and time-consuming means of evaluating the effect of mutations, making a fast and accurate in silico method highly desirable. When the structure of the wild-type complex is known, it is possible to economically evaluate the effect of point mutations with knowledge based potentials, which do not model backbone flexibility, but these have been validated only for single mutants. Substitution mutations tend to induce local conformational rearrangements only. Accordingly, ZEMu (Zone Equilibration of Mutants) flexibilizes only a small region around the site of mutation, then computes its dynamics under a physics-based force field. We validate with 1254 experimental mutants (with 1-15 simultaneous substitutions) in a wide variety of different protein environments (65 protein complexes), and obtain a significant improvement in the accuracy of predicted ΔΔG.

Towards structural systems pharmacology to study complex diseases and personalized medicine

Genome-Wide Association Studies (GWAS), whole genome sequencing, and high-throughput omics techniques have generated vast amounts of genotypic and molecular phenotypic data. However, these data have not yet been fully explored to improve the effectiveness and efficiency of drug discovery, which continues along a one-drug-one-target-one-disease paradigm. As a partial consequence, both the cost to launch a new drug and the attrition rate are increasing. Systems pharmacology and pharmacogenomics are emerging to exploit the available data and potentially reverse this trend, but, as we argue here, more is needed. To understand the impact of genetic, epigenetic, and environmental factors on drug action, we must study the structural energetics and dynamics of molecular interactions in the context of the whole human genome and interactome. Such an approach requires an integrative modeling framework for drug action that leverages advances in data-driven statistical modeling and mechanism-based multiscale modeling and transforms heterogeneous data from GWAS, high-throughput sequencing, structural genomics, functional genomics, and chemical genomics into unified knowledge. This is not a small task, but, as reviewed here, progress is being made towards the final goal of personalized medicines for the treatment of complex diseases.

Coarse-grain modelling of protein-protein interactions

Here, we review recent advances towards the modelling of protein-protein interactions (PPI) at the coarse-grained (CG) level, a technique that is now widely used to understand protein affinity, aggregation and self-assembly behaviour. PPI models of soluble proteins and membrane proteins are separately described, but we note the parallel development that is present in both research fields with three important themes: firstly, combining CG modelling with knowledge-based approaches to predict and refine protein-protein complexes; secondly, using physics-based CG models for de novo prediction of protein-protein complexes; and thirdly modelling of large scale protein aggregates.

Computational modeling of protein-RNA complex structures

Protein-RNA interactions play fundamental roles in many biological processes, such as regulation of gene expression, RNA splicing, and protein synthesis. The understanding of these processes improves as new structures of protein-RNA complexes are solved and the molecular details of interactions analyzed. However, experimental determination of protein-RNA complex structures by high-resolution methods is tedious and difficult. Therefore, studies on protein-RNA recognition and complex formation present major technical challenges for macromolecular structural biology. Alternatively, protein-RNA interactions can be predicted by computational methods. Although less accurate than experimental measurements, theoretical models of macromolecular structures can be sufficiently accurate to prompt functional hypotheses and guide e.g. identification of important amino acid or nucleotide residues. In this article we present an overview of strategies and methods for computational modeling of protein-RNA complexes, including software developed in our laboratory, and illustrate it with practical examples of structural predictions.

Fast fitting to low resolution density maps: elucidating large-scale motions of the ribosome

Determining the conformational rearrangements of large macromolecules is challenging experimentally and computationally. Case in point is the ribosome; it has been observed by high-resolution crystallography in several states, but many others are known only from low-resolution methods including cryo-electron microscopy. Combining these data into dynamical trajectories that may aid understanding of its largest-scale conformational changes has so far remained out of reach of computational methods. Most existing methods either model all atoms explicitly, resulting in often prohibitive cost, or use approximations that lose interesting structural and dynamical detail. In this work, I introduce Internal Coordinate Flexible Fitting, which uses full atomic forces and flexibility in limited regions of a model, capturing extensive conformational rearrangements at low cost. I use it to turn multiple low-resolution density maps, crystallographic structures and biochemical information into unified all-atoms trajectories of ribosomal translocation. Internal Coordinate Flexible Fitting is three orders of magnitude faster than the most comparable existing method.

Coarse-grained simulations of protein-protein association: an energy landscape perspective

Understanding protein-protein association is crucial in revealing the molecular basis of many biological processes. Here, we describe a theoretical simulation pipeline to study protein-protein association from an energy landscape perspective. First, a coarse-grained model is implemented and its applications are demonstrated via molecular dynamics simulations for several protein complexes. Second, an enhanced search method is used to efficiently sample a broad range of protein conformations. Third, multiple conformations are identified and clustered from simulation data and further projected on a three-dimensional globe specifying protein orientations and interacting energies. Results from several complexes indicate that the crystal-like conformation is favorable on the energy landscape even if the landscape is relatively rugged with metastable conformations. A closer examination on molecular forces shows that the formation of associated protein complexes can be primarily electrostatics-driven, hydrophobics-driven, or a combination of both in stabilizing specific binding interfaces. Taken together, these results suggest that the coarse-grained simulations and analyses provide an alternative toolset to study protein-protein association occurring in functional biomolecular complexes.