Modeling of proteins and their assemblies with the integrative modeling platform

Detection of autism spectrum disorder-related pathogenic trio variants by a novel structure-based approach

BACKGROUND

Glutamatergic synapse dysfunction is believed to underlie the development of Autism Spectrum Disorder (ASD) and Intellectual Disability (ID) in many individuals. However, identification of genetic markers that contribute to synaptic dysfunction in these individuals is notoriously difficult. Based on genomic analysis, structural modeling, and functional data, we recently established the involvement of the TRIO-RAC1 pathway in ASD and ID. Furthermore, we identified a pathological de novo missense mutation hotspot in TRIO's GEF1 domain. ASD/ID-related missense mutations within this domain compromise glutamatergic synapse function and likely contribute to the development of ASD/ID. The number of ASD/ID cases with mutations identified within TRIO's GEF1 domain is increasing. However, tools for accurately predicting whether such mutations are detrimental to protein function are lacking.

METHODS

Here we deployed advanced protein structural modeling techniques to predict potential de novo pathogenic and benign mutations within TRIO's GEF1 domain. Mutant TRIO-9 constructs were generated and expressed in CA1 pyramidal neurons of organotypic cultured hippocampal slices. AMPA receptor-mediated postsynaptic currents were examined in these neurons using dual whole-cell patch clamp electrophysiology. We also validated these findings using orthogonal co-immunoprecipitation and fluorescence lifetime imaging (FLIM-FRET) experiments to assay TRIO mutant overexpression effects on TRIO-RAC1 binding and on RAC1 activity in HEK293/T cells.

RESULTS

Missense mutations in TRIO's GEF1 domain that were predicted to disrupt TRIO-RAC1 binding or stability were tested experimentally and found to greatly impair TRIO-9's influence on glutamatergic synapse function. In contrast, missense mutations in TRIO's GEF1 domain that were predicted to have minimal effect on TRIO-RAC1 binding or stability did not impair TRIO-9's influence on glutamatergic synapse function in our experimental assays. In orthogonal assays, we find most of the mutations predicted to disrupt binding display loss of function but mutants predicted to disrupt stability do not reflect our results from neuronal electrophysiological data.

LIMITATIONS

We present a method to predict missense mutations in TRIO's GEF1 domain that may compromise TRIO function and test for effects in a limited number of assays. Possible limitations arising from the model systems employed here can be addressed in future studies. Our method does not provide evidence for whether these mutations confer ASD/ID risk or the likelihood that such mutations will result in the development of ASD/ID.

CONCLUSIONS

Here we show that a combination of structure-based computational predictions and experimental validation can be employed to reliably predict whether missense mutations in the human TRIO gene impede TRIO protein function and compromise TRIO's role in glutamatergic synapse regulation. With the growing accessibility of genome sequencing, the use of such tools in the accurate identification of pathological mutations will be instrumental in diagnostics of ASD/ID.

Using Integrative Modeling Platform to compute, validate, and archive a model of a protein complex structure

Biology is advanced by producing structural models of biological systems, such as protein complexes. Some systems are recalcitrant to traditional structure determination methods. In such cases, it may still be possible to produce useful models by integrative structure determination that depends on simultaneous use of multiple types of data. An ensemble of models that are sufficiently consistent with the data is produced by a structural sampling method guided by a data-dependent scoring function. The variation in the ensemble of models quantified the uncertainty of the structure, generally resulting from the uncertainty in the input information and actual structural heterogeneity in the samples used to produce the data. Here, we describe how to generate, assess, and interpret ensembles of integrative structural models using our open source Integrative Modeling Platform program (https://integrativemodeling.org).

Exposing the Elusive Exocyst Structure

A major challenge for a molecular understanding of membrane trafficking has been the elucidation of high-resolution structures of large, multisubunit tethering complexes that spatially and temporally control intracellular membrane fusion. Exocyst is a large hetero-octameric protein complex proposed to tether secretory vesicles at the plasma membrane to provide quality control of soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE)-mediated membrane fusion. Breakthroughs in methodologies, including sample preparation, biochemical characterization, fluorescence microscopy, and single-particle cryoelectron microscopy, are providing critical insights into the structure and function of the exocyst. These studies now pose more questions than answers for understanding fundamental functional mechanisms, and they open wide the door for future studies to elucidate interactions with protein and membrane partners, potential conformational changes, and molecular insights into tethering reactions.

Iterative Molecular Dynamics-Rosetta Membrane Protein Structure Refinement Guided by Cryo-EM Densities

Knowing atomistic details of proteins is essential not only for the understanding of protein function but also for the development of drugs. Experimental methods such as X-ray crystallography, NMR, and cryo-electron microscopy (cryo-EM) are the preferred forms of protein structure determination and have achieved great success over the most recent decades. Computational methods may be an alternative when experimental techniques fail. However, computational methods are severely limited when it comes to predicting larger macromolecule structures with little sequence similarity to known structures. The incorporation of experimental restraints in computational methods is becoming increasingly important to more reliably predict protein structure. One such experimental input used in structure prediction and refinement is cryo-EM densities. Recent advances in cryo-EM have arguably revolutionized the field of structural biology. Our previously developed cryo-EM-guided Rosetta-MD protocol has shown great promise in the refinement of soluble protein structures. In this study, we extended cryo-EM density-guided iterative Rosetta-MD to membrane proteins. We also improved the methodology in general by picking models based on a combination of their score and fit-to-density during the Rosetta model selection. By doing so, we have been able to pick models superior to those with the previous selection based on Rosetta score only and we have been able to further improve our previously refined models of soluble proteins. The method was tested with five membrane spanning protein structures. By applying density-guided Rosetta-MD iteratively we were able to refine the predicted structures of these membrane proteins to atomic resolutions. We also showed that the resolution of the density maps determines the improvement and quality of the refined models. By incorporating high-resolution density maps (∼4 Å), we were able to more significantly improve the quality of the models than when medium-resolution maps (6.9 Å) were used. Beginning from an average starting structure root mean square deviation (RMSD) to native of 4.66 Å, our protocol was able to refine the structures to bring the average refined structure RMSD to 1.66 Å when 4 Å density maps were used. The protocol also successfully refined the HIV-1 CTD guided by an experimental 5 Å density map.

Applications of contact predictions to structural biology

Evolutionary pressure on residue interactions, intramolecular or intermolecular, that are important for protein structure or function can lead to covariance between the two positions. Recent methodological advances allow much more accurate contact predictions to be derived from this evolutionary covariance signal. The practical application of contact predictions has largely been confined to structural bioinformatics, yet, as this work seeks to demonstrate, the data can be of enormous value to the structural biologist working in X-ray crystallo-graphy, cryo-EM or NMR. Integrative structural bioinformatics packages such as Rosetta can already exploit contact predictions in a variety of ways. The contribution of contact predictions begins at construct design, where structural domains may need to be expressed separately and contact predictions can help to predict domain limits. Structure solution by molecular replacement (MR) benefits from contact predictions in diverse ways: in difficult cases, more accurate search models can be constructed using ab initio modelling when predictions are available, while intermolecular contact predictions can allow the construction of larger, oligomeric search models. Furthermore, MR using supersecondary motifs or large-scale screens against the PDB can exploit information, such as the parallel or antiparallel nature of any β-strand pairing in the target, that can be inferred from contact predictions. Contact information will be particularly valuable in the determination of lower resolution structures by helping to assign sequence register. In large complexes, contact information may allow the identity of a protein responsible for a certain region of density to be determined and then assist in the orientation of an available model within that density. In NMR, predicted contacts can provide long-range information to extend the upper size limit of the technique in a manner analogous but complementary to experimental methods. Finally, predicted contacts can distinguish between biologically relevant interfaces and mere lattice contacts in a final crystal structure, and have potential in the identification of functionally important regions and in foreseeing the consequences of mutations.

Docking Screens for Novel Ligands Conferring New Biology

It is now plausible to dock libraries of 10 million molecules against targets over several days or weeks. When the molecules screened are commercially available, they may be rapidly tested to find new leads. Although docking retains important liabilities (it cannot calculate affinities accurately nor even reliably rank order high-scoring molecules), it can often can distinguish likely from unlikely ligands, often with hit rates above 10%. Here we summarize the improvements in libraries, target quality, and methods that have supported these advances, and the open access resources that make docking accessible. Recent docking screens for new ligands are sketched, as are the binding, crystallographic, and in vivo assays that support them. Like any technique, controls are crucial, and key experimental ones are reviewed. With such controls, docking campaigns can find ligands with new chemotypes, often revealing the new biology that may be docking's greatest impact over the next few years.

Assembly and Molecular Architecture of the Phosphoinositide 3-Kinase p85α Homodimer

Phosphoinositide 3-kinases (PI3Ks) are a family of lipid kinases that are activated by growth factor and G-protein-coupled receptors and propagate intracellular signals for growth, survival, proliferation, and metabolism. p85α, a modular protein consisting of five domains, binds and inhibits the enzymatic activity of class IA PI3K catalytic subunits. Here, we describe the structural states of the p85α dimer, based on data from in vivo and in vitro solution characterization. Our in vitro assembly and structural analyses have been enabled by the creation of cysteine-free p85α that is functionally equivalent to native p85α. Analytical ultracentrifugation studies showed that p85α undergoes rapidly reversible monomer-dimer assembly that is highly exothermic in nature. In addition to the documented SH3-PR1 dimerization interaction, we identified a second intermolecular interaction mediated by cSH2 domains at the C-terminal end of the polypeptide. We have demonstrated in vivo concentration-dependent dimerization of p85α using fluorescence fluctuation spectroscopy. Finally, we have defined solution conditions under which the protein is predominantly monomeric or dimeric, providing the basis for small angle x-ray scattering and chemical cross-linking structural analysis of the discrete dimer. These experimental data have been used for the integrative structure determination of the p85α dimer. Our study provides new insight into the structure and assembly of the p85α homodimer and suggests that this protein is a highly dynamic molecule whose conformational flexibility allows it to transiently associate with multiple binding proteins.

Function and dynamics of macromolecular complexes explored by integrative structural and computational biology

Three vignettes exemplify the potential of combining EM and X-ray crystallographic data with molecular dynamics (MD) simulation to explore the architecture, dynamics and functional properties of multicomponent, macromolecular complexes. The first two describe how EM and X-ray crystallography were used to solve structures of the ribosome and the Arp2/3-actin complex, which enabled MD simulations that elucidated functional dynamics. The third describes how EM, X-ray crystallography, and microsecond MD simulations of a GPCR:G protein complex were used to explore transmembrane signaling by the β-adrenergic receptor. Recent technical advancements in EM, X-ray crystallography and computational simulation create unprecedented synergies for integrative structural biology to reveal new insights into heretofore intractable biological systems.

Accurate SAXS profile computation and its assessment by contrast variation experiments

A major challenge in structural biology is to characterize structures of proteins and their assemblies in solution. At low resolution, such a characterization may be achieved by small angle x-ray scattering (SAXS). Because SAXS analyses often require comparing profiles calculated from many atomic models against those determined by experiment, rapid and accurate profile computation from molecular structures is needed. We developed fast open-source x-ray scattering (FoXS) for profile computation. To match the experimental profile within the experimental noise, FoXS explicitly computes all interatomic distances and implicitly models the first hydration layer of the molecule. For assessing the accuracy of the modeled hydration layer, we performed contrast variation experiments for glucose isomerase and lysozyme, and found that FoXS can accurately represent density changes of this layer. The hydration layer model was also compared with a SAXS profile calculated for the explicit water molecules in the high-resolution structures of glucose isomerase and lysozyme. We tested FoXS on eleven protein, one DNA, and two RNA structures, revealing superior accuracy and speed versus CRYSOL, AquaSAXS, the Zernike polynomials-based method, and Fast-SAXS-pro. In addition, we demonstrated a significant correlation of the SAXS score with the accuracy of a structural model. Moreover, FoXS utility for analyzing heterogeneous samples was demonstrated for intrinsically flexible XLF-XRCC4 filaments and Ligase III-DNA complex. FoXS is extensively used as a standalone web server as a component of integrative structure determination by programs IMP, Chimera, and BILBOMD, as well as in other applications that require rapidly and accurately calculated SAXS profiles.

Protein modeling: what happened to the "protein structure gap"?

Computational modeling of three-dimensional macromolecular structures and complexes from their sequence has been a long-standing vision in structural biology. Over the last 2 decades, a paradigm shift has occurred: starting from a large "structure knowledge gap" between the huge number of protein sequences and small number of known structures, today, some form of structural information, either experimental or template-based models, is available for the majority of amino acids encoded by common model organism genomes. With the scientific focus of interest moving toward larger macromolecular complexes and dynamic networks of interactions, the integration of computational modeling methods with low-resolution experimental techniques allows the study of large and complex molecular machines. One of the open challenges for computational modeling and prediction techniques is to convey the underlying assumptions, as well as the expected accuracy and structural variability of a specific model, which is crucial to understanding its limitations.

Modeling of proteins and their assemblies with the Integrative Modeling Platform

To understand the workings of the living cell, we need to characterize protein assemblies that constitute the cell (for example, the ribosome, 26S proteasome, and the nuclear pore complex). A reliable high-resolution structural characterization of these assemblies is frequently beyond the reach of current experimental methods, such as X-ray crystallography, NMR spectroscopy, electron microscopy, footprinting, chemical cross-linking, FRET spectroscopy, small angle X-ray scattering, and proteomics. However, the information garnered from different methods can be combined and used to build models of the assembly structures that are consistent with all of the available datasets, and therefore more accurate, precise, and complete. Here, we describe a protocol for this integration, whereby the information is converted to a set of spatial restraints and a variety of optimization procedures can be used to generate models that satisfy the restraints as well as possible. These generated models can then potentially inform about the precision and accuracy of structure determination, the accuracy of the input datasets, and further data generation. We also demonstrate the Integrative Modeling Platform (IMP) software, which provides the necessary computational framework to implement this protocol, and several applications for specific use cases.

Probabilistic determination of probe locations from distance data

Distance constraints, in principle, can be employed to determine information about the location of probes within a three-dimensional volume. Traditional methods for locating probes from distance constraints involve optimization of scoring functions that measure how well the probe location fits the distance data, exploring only a small subset of the scoring function landscape in the process. These methods are not guaranteed to find the global optimum and provide no means to relate the identified optimum to all other optima in scoring space. Here, we introduce a method for the location of probes from distance information that is based on probability calculus. This method allows exploration of the entire scoring space by directly combining probability functions representing the distance data and information about attachment sites. The approach is guaranteed to identify the global optimum and enables the derivation of confidence intervals for the probe location as well as statistical quantification of ambiguities. We apply the method to determine the location of a fluorescence probe using distances derived by FRET and show that the resulting location matches that independently derived by electron microscopy.

Integrative structural modelling of the cardiac thin filament: energetics at the interface and conservation patterns reveal a spotlight on period 2 of tropomyosin

Cardiomyopathies are a major health problem, with inherited cardiomyopathies, many of which are caused by mutations in genes encoding sarcomeric proteins, constituting an ever-increasing fraction of cases. To begin to study the mechanisms by which these mutations cause disease, we have employed an integrative modelling approach to study the interactions between tropomyosin and actin. Starting from the existing blocked state model, we identified a specific zone on the actin surface which is highly favourable to support tropomyosin sliding from the blocked/closed states to the open state. We then analysed the predicted actin-tropomyosin interface regions for the three states. Each quasi-repeat of tropomyosin was studied for its interaction strength and evolutionary conservation to focus on smaller surface zones. Finally, we show that the distribution of the known cardiomyopathy mutations of α-tropomyosin is consistent with our model. This analysis provides structural insights into the possible mode of interactions between tropomyosin and actin in the open state for the first time.

Integrative structural modeling with small angle X-ray scattering profiles

Recent technological advances enabled high-throughput collection of Small Angle X-ray Scattering (SAXS) profiles of biological macromolecules. Thus, computational methods for integrating SAXS profiles into structural modeling are needed more than ever. Here, we review specifically the use of SAXS profiles for the structural modeling of proteins, nucleic acids, and their complexes. First, the approaches for computing theoretical SAXS profiles from structures are presented. Second, computational methods for predicting protein structures, dynamics of proteins in solution, and assembly structures are covered. Third, we discuss the use of SAXS profiles in integrative structure modeling approaches that depend simultaneously on several data types.

Networks of ProteinProtein Interactions: From Uncertainty to Molecular Details

Proteins are the bricks and mortar of cells. The work of proteins is structural and functional, as they are the principal element of the organization of the cell architecture, but they also play a relevant role in its metabolism and regulation. To perform all these functions, proteins need to interact with each other and with other bio-molecules, either to form complexes or to recognize precise targets of their action. For instance, a particular transcription factor may activate one gene or another depending on its interactions with other proteins and not only with DNA. Hence, the ability of a protein to interact with other bio-molecules, and the partners they have at each particular time and location can be crucial to characterize the role of a protein. Proteins rarely act alone; they rather constitute a mingled network of physical interactions or other types of relationships (such as metabolic and regulatory) or signaling cascades. In this context, understanding the function of a protein implies to recognize the members of its neighborhood and to grasp how they associate, both at the systemic and atomic level. The network of physical interactions between the proteins of a system, cell or organism, is defined as the interactome. The purpose of this review is to deepen the description of interactomes at different levels of detail: from the molecular structure of complexes to the global topology of the network of interactions. The approaches and techniques applied experimentally and computationally to attain each level are depicted. The limits of each technique and its integration into a model network, the challenges and actual problems of completeness of an interactome, and the reliability of the interactions are reviewed and summarized. Finally, the application of the current knowledge of protein-protein interactions on modern network medicine and protein function annotation is also explored.

Genome structure determination via 3C-based data integration by the Integrative Modeling Platform

The three-dimensional (3D) architecture of a genome determines the spatial localization of regulatory elements and the genes they regulate. Thus, elucidating the 3D structure of a genome may result in significant insights about how genes are regulated. The current state-of-the art in experimental methods, including light microscopy and cell/molecular biology, are now able to provide detailed information on the position of genes and their interacting partners. However, such methods by themselves are not able to determine the high-resolution 3D structure of genomes or genomic domains. Here we describe a computational module of the Integrative Modeling Platform (IMP, http://www.integrativemodeling.org) that uses chromosome conformation capture data to determine the 3D architecture of genomic domains and entire genomes at unprecedented resolutions. This approach, through the visualization of looping interactions between distal regulatory elements, allows characterizing global chromatin features and their relation to gene expression. We illustrate our work by outlining the determination of the 3D architecture of the α-globin domain in the human genome.

Excited states of ribosome translocation revealed through integrative molecular modeling

The dynamic nature of biomolecules leads to significant challenges when characterizing the structural properties associated with function. While X-ray crystallography and imaging techniques (such as cryo-electron microscopy) can reveal the structural details of stable molecular complexes, strategies must be developed to characterize configurations that exhibit only marginal stability (such as intermediates) or configurations that do not correspond to minima on the energy landscape (such as transition-state ensembles). Here, we present a methodology (MDfit) that utilizes molecular dynamics simulations to generate configurations of excited states that are consistent with available biophysical and biochemical measurements. To demonstrate the approach, we present a sequence of configurations that are suggested to be associated with transfer RNA (tRNA) movement through the ribosome (translocation). The models were constructed by combining information from X-ray crystallography, cryo-electron microscopy, and biochemical data. These models provide a structural framework for translocation that may be further investigated experimentally and theoretically to determine the precise energetic character of each configuration and the transition dynamics between them.