Bridging from molecular simulation to biochemical networks

Structural Systems Biology Toolkit (SSBtoolkit): From Molecular Structure to Subcellular Signaling Pathways

Here, we introduce the Structural Systems Biology (SSB) toolkit, a Python library that integrates structural macromolecular data with systems biology simulations to model signal-transduction pathways of G-protein-coupled receptors (GPCRs). Our framework streamlines simulation and analysis of the mathematical models of GPCRs cellular pathways, facilitating the exploration of the signal-transduction kinetics induced by ligand-GPCR interactions: the dose-response of the ligand can be modeled, along with the corresponding change in the concentration of other signaling molecular species over time, like for instance [Ca2+] or [cAMP]. SSB toolkit brings to light the possibility of easily investigating the subcellular effects of ligand binding on receptor activation, even in the presence of genetic mutations, thereby enhancing our understanding of the intricate relationship between ligand-target interactions at the molecular level and the higher-level cellular and (patho)physiological response mechanisms.

Combining hypothesis- and data-driven neuroscience modeling in FAIR workflows

Modeling in neuroscience occurs at the intersection of different points of view and approaches. Typically, hypothesis-driven modeling brings a question into focus so that a model is constructed to investigate a specific hypothesis about how the system works or why certain phenomena are observed. Data-driven modeling, on the other hand, follows a more unbiased approach, with model construction informed by the computationally intensive use of data. At the same time, researchers employ models at different biological scales and at different levels of abstraction. Combining these models while validating them against experimental data increases understanding of the multiscale brain. However, a lack of interoperability, transparency, and reusability of both models and the workflows used to construct them creates barriers for the integration of models representing different biological scales and built using different modeling philosophies. We argue that the same imperatives that drive resources and policy for data - such as the FAIR (Findable, Accessible, Interoperable, Reusable) principles - also support the integration of different modeling approaches. The FAIR principles require that data be shared in formats that are Findable, Accessible, Interoperable, and Reusable. Applying these principles to models and modeling workflows, as well as the data used to constrain and validate them, would allow researchers to find, reuse, question, validate, and extend published models, regardless of whether they are implemented phenomenologically or mechanistically, as a few equations or as a multiscale, hierarchical system. To illustrate these ideas, we use a classical synaptic plasticity model, the Bienenstock-Cooper-Munro rule, as an example due to its long history, different levels of abstraction, and implementation at many scales.

Multiscale kinetic analysis of proteins

The stochasticity of molecular motion results in the existence of multiple kinetically relevant pathways in many biomolecular mechanisms. Because it is highly demanding to characterize them for complex systems, mechanisms are often described with a single-pathway perspective. However, kinetic network analysis and sub-ensemble experimental insight are increasingly demonstrating not only the existence of competing pathways but also the importance of kinetic selection in biology. This review focuses on advances in multiscale kinetic analysis of proteins, which connects molecular level information from simulations to macroscopic data to characterize mechanistic reaction networks and the reactive flux through them. We describe a range of methods used and highlight several examples where kinetic modeling has revealed functional importance of pathway heterogeneity.

Regulation of adenylyl cyclase 5 in striatal neurons confers the ability to detect coincident neuromodulatory signals

Long-term potentiation and depression of synaptic activity in response to stimuli is a key factor in reinforcement learning. Strengthening of the corticostriatal synapses depends on the second messenger cAMP, whose synthesis is catalysed by the enzyme adenylyl cyclase 5 (AC5), which is itself regulated by the stimulatory Gαolf and inhibitory Gαi proteins. AC isoforms have been suggested to act as coincidence detectors, promoting cellular responses only when convergent regulatory signals occur close in time. However, the mechanism for this is currently unclear, and seems to lie in their diverse regulation patterns. Despite attempts to isolate the ternary complex, it is not known if Gαolf and Gαi can bind to AC5 simultaneously, nor what activity the complex would have. Using protein structure-based molecular dynamics simulations, we show that this complex is stable and inactive. These simulations, along with Brownian dynamics simulations to estimate protein association rates constants, constrain a kinetic model that shows that the presence of this ternary inactive complex is crucial for AC5's ability to detect coincident signals, producing a synergistic increase in cAMP. These results reveal some of the prerequisites for corticostriatal synaptic plasticity, and explain recent experimental data on cAMP concentrations following receptor activation. Moreover, they provide insights into the regulatory mechanisms that control signal processing by different AC isoforms.

Computational algorithms for in silico profiling of activating mutations in cancer

Methods to catalog and computationally assess the mutational landscape of proteins in human cancers are desirable. One approach is to adapt evolutionary or data-driven methods developed for predicting whether a single-nucleotide polymorphism (SNP) is deleterious to protein structure and function. In cases where understanding the mechanism of protein activation and regulation is desired, an alternative approach is to employ structure-based computational approaches to predict the effects of point mutations. Through a case study of mutations in kinase domains of three proteins, namely, the anaplastic lymphoma kinase (ALK) in pediatric neuroblastoma patients, serine/threonine-protein kinase B-Raf (BRAF) in melanoma patients, and erythroblastic oncogene B 2 (ErbB2 or HER2) in breast cancer patients, we compare the two approaches above. We find that the structure-based method is most appropriate for developing a binary classification of several different mutations, especially infrequently occurring ones, concerning the activation status of the given target protein. This approach is especially useful if the effects of mutations on the interactions of inhibitors with the target proteins are being sought. However, many patients will present with mutations spread across different target proteins, making structure-based models computationally demanding to implement and execute. In this situation, data-driven methods-including those based on machine learning techniques and evolutionary methods-are most appropriate for recognizing and illuminate mutational patterns. We show, however, that, in the present status of the field, the two methods have very different accuracies and confidence values, and hence, the optimal choice of their deployment is context-dependent.

Integrating Structural Information to Study the Dynamics of Protein-Protein Interactions in Cells

The information of how two proteins interact is embedded in the atomic details of their binding interfaces. These interactions, spatial-temporally coordinating each other as a network in a variable cytoplasmic environment, dominate almost all biological functions. A feasible and reliable computational model is highly demanded to realistically simulate these cellular processes and unravel the complexities beneath them. We therefore present a multiscale framework that integrates simulations on two different scales. The higher-resolution model incorporates structural information of proteins and energetics of their binding, while the lower-resolution model uses a highly simplified representation of proteins to capture the long-time-scale dynamics of a system with multiple proteins. Through a systematic benchmark test and two practical applications of biomolecular systems with specific cellular functions, we demonstrated that this method could be a powerful approach to understand molecular mechanisms of dynamic interactions between biomolecules and their functional impacts with high computational efficiency.

The Role of Multiscale Protein Dynamics in Antigen Presentation and T Lymphocyte Recognition

T lymphocytes are stimulated when they recognize short peptides bound to class I proteins of the major histocompatibility complex (MHC) protein, as peptide-MHC complexes. Due to the diversity in T-cell receptor (TCR) molecules together with both the peptides and MHC proteins they bind to, it has been difficult to design vaccines and treatments based on these interactions. Machine learning has made some progress in trying to predict the immunogenicity of peptide sequences in the context of specific MHC class I alleles but, as such approaches cannot integrate temporal information and lack explanatory power, their scope will always be limited. Here, we advocate a mechanistic description of antigen presentation and TCR activation which is explanatory, predictive, and quantitative, drawing on modeling approaches that collectively span several length and time scales, being capable of furnishing reliable biological descriptions that are difficult for experimentalists to provide. It is a form of multiscale systems biology. We propose the use of chemical rate equations to describe the time evolution of the foreign and host proteins to explain how the original proteins end up being presented on the cell surface as peptide fragments, while we invoke molecular dynamics to describe the key binding processes on the molecular level, including those of peptide-MHC complexes with TCRs which lie at the heart of the immune response. On each level, complementary methods based on machine learning are available, and we discuss the relationship between these divergent approaches. The pursuit of predictive mechanistic modeling approaches requires experimentalists to adapt their work so as to acquire, store, and expose data that can be used to verify and validate such models.

Dynamic force matching: A method for constructing dynamical coarse-grained models with realistic time dependence

Coarse-grained (CG) models of molecular systems, with fewer mechanical degrees of freedom than an all-atom model, are used extensively in chemical physics. It is generally accepted that a coarse-grained model that accurately describes equilibrium structural properties (as a result of having a well constructed CG potential energy function) does not necessarily exhibit appropriate dynamical behavior when simulated using conservative Hamiltonian dynamics for the CG degrees of freedom on the CG potential energy surface. Attempts to develop accurate CG dynamic models usually focus on replacing Hamiltonian motion by stochastic but Markovian dynamics on that surface, such as Langevin or Brownian dynamics. However, depending on the nature of the system and the extent of the coarse-graining, a Markovian dynamics for the CG degrees of freedom may not be appropriate. In this paper, we consider the problem of constructing dynamic CG models within the context of the Multi-Scale Coarse-graining (MS-CG) method of Voth and coworkers. We propose a method of converting a MS-CG model into a dynamic CG model by adding degrees of freedom to it in the form of a small number of fictitious particles that interact with the CG degrees of freedom in simple ways and that are subject to Langevin forces. The dynamic models are members of a class of nonlinear systems interacting with special heat baths that were studied by Zwanzig [J. Stat. Phys. 9, 215 (1973)]. The properties of the fictitious particles can be inferred from analysis of the dynamics of all-atom simulations of the system of interest. This is analogous to the fact that the MS-CG method generates the CG potential from analysis of equilibrium structures observed in all-atom simulation data. The dynamic models generate a non-Markovian dynamics for the CG degrees of freedom, but they can be easily simulated using standard molecular dynamics programs. We present tests of this method on a series of simple examples that demonstrate that the method provides realistic dynamical CG models that have non-Markovian or close to Markovian behavior that is consistent with the actual dynamical behavior of the all-atom system used to construct the CG model. Both the construction and the simulation of such a dynamic CG model have computational requirements that are similar to those of the corresponding MS-CG model and are good candidates for CG modeling of very large systems.

Achievements and challenges in structural bioinformatics and computational biophysics

Motivation: The field of structural bioinformatics and computational biophysics has undergone a revolution in the last 10 years. Developments that are captured annually through the 3DSIG meeting, upon which this article reflects.

Results: An increase in the accessible data, computational resources and methodology has resulted in an increase in the size and resolution of studied systems and the complexity of the questions amenable to research. Concomitantly, the parameterization and efficiency of the methods have markedly improved along with their cross-validation with other computational and experimental results.

Conclusion: The field exhibits an ever-increasing integration with biochemistry, biophysics and other disciplines. In this article, we discuss recent achievements along with current challenges within the field.

Contact:Rafael.Najmanovich@USherbrooke.ca

The Theory of Ultra-Coarse-Graining. 2. Numerical Implementation

The increasing interest in the modeling of complex macromolecular systems in recent years has spurred the development of numerous coarse-graining (CG) techniques. However, many of the CG models are constructed assuming that all details beneath the resolution of CG degrees of freedom are fast and average out, which sets limits on the resolution of feasible coarse-grainings and on the range of applications of the CG models. Ultra-coarse-graining (UCG) makes it possible to construct models at any desired resolution while accounting for discrete conformational or chemical changes within the CG sites that can modulate the interactions between them. Here, we discuss the UCG methodology and its numerical implementation. We pay particular attention to the numerical mechanism for including state transitions between different conformations within CG sites because this has not been discussed previously. Using a simple example of 1,2-dichloroethane, we demonstrate the ability of the UCG model to reproduce the multiconfigurational behavior of this molecular liquid, even when each molecule is modeled with only one CG site. The methodology can also be applied to other molecular liquids and macromolecular systems with time scale separation between conformational transitions and other intramolecular motions and rotations.

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.

The scoring of poses in protein-protein docking: current capabilities and future directions

BACKGROUND

Protein-protein docking, which aims to predict the structure of a protein-protein complex from its unbound components, remains an unresolved challenge in structural bioinformatics. An important step is the ranking of docked poses using a scoring function, for which many methods have been developed. There is a need to explore the differences and commonalities of these methods with each other, as well as with functions developed in the fields of molecular dynamics and homology modelling.

RESULTS

We present an evaluation of 115 scoring functions on an unbound docking decoy benchmark covering 118 complexes for which a near-native solution can be found, yielding top 10 success rates of up to 58%. Hierarchical clustering is performed, so as to group together functions which identify near-natives in similar subsets of complexes. Three set theoretic approaches are used to identify pairs of scoring functions capable of correctly scoring different complexes. This shows that functions in different clusters capture different aspects of binding and are likely to work together synergistically.

CONCLUSIONS

All functions designed specifically for docking perform well, indicating that functions are transferable between sampling methods. We also identify promising methods from the field of homology modelling. Further, differential success rates by docking difficulty and solution quality suggest a need for flexibility-dependent scoring. Investigating pairs of scoring functions, the set theoretic measures identify known scoring strategies as well as a number of novel approaches, indicating promising augmentations of traditional scoring methods. Such augmentation and parameter combination strategies are discussed in the context of the learning-to-rank paradigm.

Characterizing changes in the rate of protein-protein dissociation upon interface mutation using hotspot energy and organization

Predicting the effects of mutations on the kinetic rate constants of protein-protein interactions is central to both the modeling of complex diseases and the design of effective peptide drug inhibitors. However, while most studies have concentrated on the determination of association rate constants, dissociation rates have received less attention. In this work we take a novel approach by relating the changes in dissociation rates upon mutation to the energetics and architecture of hotspots and hotregions, by performing alanine scans pre- and post-mutation. From these scans, we design a set of descriptors that capture the change in hotspot energy and distribution. The method is benchmarked on 713 kinetically characterized mutations from the SKEMPI database. Our investigations show that, with the use of hotspot descriptors, energies from single-point alanine mutations may be used for the estimation of off-rate mutations to any residue type and also multi-point mutations. A number of machine learning models are built from a combination of molecular and hotspot descriptors, with the best models achieving a Pearson's Correlation Coefficient of 0.79 with experimental off-rates and a Matthew's Correlation Coefficient of 0.6 in the detection of rare stabilizing mutations. Using specialized feature selection models we identify descriptors that are highly specific and, conversely, broadly important to predicting the effects of different classes of mutations, interface regions and complexes. Our results also indicate that the distribution of the critical stability regions across protein-protein interfaces is a function of complex size more strongly than interface area. In addition, mutations at the rim are critical for the stability of small complexes, but consistently harder to characterize. The relationship between hotregion size and the dissociation rate is also investigated and, using hotspot descriptors which model cooperative effects within hotregions, we show how the contribution of hotregions of different sizes, changes under different cooperative effects.

Reaching new levels of realism in modeling biological macromolecules in cellular environments

An increasing number of studies are aimed at modeling cellular environments in a comprehensive and realistic fashion. A major challenge in these efforts is how to bridge spatial and temporal scales over many orders of magnitude. Furthermore, there are additional challenges in integrating different aspects ranging from questions about biomolecular stability in crowded environments to the description of reactive processes on cellular scales. In this review, recent studies with models of biomolecules in cellular environments at different levels of detail are discussed in terms of their strengths and weaknesses. In particular, atomistic models, implicit representations of cellular environments, coarse-grained and spheroidal models of biomolecules, as well as the inclusion of reactive processes via reaction-diffusion models are described. Furthermore, strategies for integrating the different models into a comprehensive description of cellular environments are discussed.

A structural systems biology approach for quantifying the systemic consequences of missense mutations in proteins

Gauging the systemic effects of non-synonymous single nucleotide polymorphisms (nsSNPs) is an important topic in the pursuit of personalized medicine. However, it is a non-trivial task to understand how a change at the protein structure level eventually affects a cell's behavior. This is because complex information at both the protein and pathway level has to be integrated. Given that the idea of integrating both protein and pathway dynamics to estimate the systemic impact of missense mutations in proteins remains predominantly unexplored, we investigate the practicality of such an approach by formulating mathematical models and comparing them with experimental data to study missense mutations. We present two case studies: (1) interpreting systemic perturbation for mutations within the cell cycle control mechanisms (G2 to mitosis transition) for yeast; (2) phenotypic classification of neuron-related human diseases associated with mutations within the mitogen-activated protein kinase (MAPK) pathway. We show that the application of simplified mathematical models is feasible for understanding the effects of small sequence changes on cellular behavior. Furthermore, we show that the systemic impact of missense mutations can be effectively quantified as a combination of protein stability change and pathway perturbation.

Kinetic rate constant prediction supports the conformational selection mechanism of protein binding

The prediction of protein-protein kinetic rate constants provides a fundamental test of our understanding of molecular recognition, and will play an important role in the modeling of complex biological systems. In this paper, a feature selection and regression algorithm is applied to mine a large set of molecular descriptors and construct simple models for association and dissociation rate constants using empirical data. Using separate test data for validation, the predicted rate constants can be combined to calculate binding affinity with accuracy matching that of state of the art empirical free energy functions. The models show that the rate of association is linearly related to the proportion of unbound proteins in the bound conformational ensemble relative to the unbound conformational ensemble, indicating that the binding partners must adopt a geometry near to that of the bound prior to binding. Mirroring the conformational selection and population shift mechanism of protein binding, the models provide a strong separate line of evidence for the preponderance of this mechanism in protein-protein binding, complementing structural and theoretical studies.

Coarse-grained (multiscale) simulations in studies of biophysical and chemical systems

Recent years have witnessed an explosion in computational power, leading to attempts to model ever more complex systems. Nevertheless, there remain cases for which the use of brute-force computer simulations is clearly not the solution. In such cases, great benefit can be obtained from the use of physically sound simplifications. The introduction of such coarse graining can be traced back to the early usage of a simplified model in studies of proteins. Since then, the field has progressed tremendously. In this review, we cover both key developments in the field and potential future directions. Additionally, particular emphasis is given to two general approaches, namely the renormalization and reference potential approaches, which allow one to move back and forth between the coarse-grained (CG) and full models, as these approaches provide the foundation for CG modeling of complex systems.

Multi-scale modeling in biology: how to bridge the gaps between scales?

Human physiological functions are regulated across many orders of magnitude in space and time. Integrating the information and dynamics from one scale to another is critical for the understanding of human physiology and the treatment of diseases. Multi-scale modeling, as a computational approach, has been widely adopted by researchers in computational and systems biology. A key unsolved issue is how to represent appropriately the dynamical behaviors of a high-dimensional model of a lower scale by a low-dimensional model of a higher scale, so that it can be used to investigate complex dynamical behaviors at even higher scales of integration. In the article, we first review the widely-used different modeling methodologies and their applications at different scales. We then discuss the gaps between different modeling methodologies and between scales, and discuss potential methods for bridging the gaps between scales.

Predicting kinetic constants of protein-protein interactions based on structural properties

Elucidating kinetic processes of protein-protein interactions (PPI) helps to understand how basic building blocks affect overall behavior of living systems. In this study, we used structure-based properties to build predictive models for kinetic constants of PPI. A highly diverse PPI dataset, protein-protein kinetic interaction data and structures (PPKIDS), was built. PPKIDS contains 62 PPI with complex structures and kinetic constants measured experimentally. The influence of structural properties on kinetics of PPI was studied using 35 structure-based features, describing different aspects of complex structures. Linear models for the prediction of kinetic constants were built by fitting with selected subsets of structure-based features. The models gave correlation coefficients of 0.801, 0.732, and 0.770 for k(off), k(on), and K(d), respectively, in leave-one-out cross validations. The predictive models reported here use only protein complex structures as input and can be generally applied in PPI studies as well as systems biology modeling. Our study confirmed that different properties play different roles in the kinetic process of PPI. For example, k(on) was affected by overall structural features of complexes, such as the composition of secondary structures, the change of translational and rotational entropy, and the electrostatic interaction; while k(off) was determined by interfacial properties, such as number of contacted atom pairs per 100 Ų. This information provides useful hints for PPI design.

Bridging the gap: linking molecular simulations and systemic descriptions of cellular compartments

Metabolic processes in biological cells are commonly either characterized at the level of individual enzymes and metabolites or at the network level. Often these two paradigms are considered as mutually exclusive because concepts from neither side are suited to describe the complete range of scales. Additionally, when modeling metabolic or regulatory cellular systems, often a large fraction of the required kinetic parameters are unknown. This even applies to such simple and extensively studied systems like the photosynthetic apparatus of purple bacteria. Using the chromatophore vesicles of Rhodobacter sphaeroides as a model system, we show that a consistent kinetic model emerges when fitting the dynamics of a molecular stochastic simulation to a set of time dependent experiments even though about two thirds of the kinetic parameters in this system are not known from experiment. Those kinetic parameters that were previously known all came out in the expected range. The simulation model was built from independent protein units composed of elementary reactions processing single metabolites. This pools-and-proteins approach naturally compiles the wealth of available molecular biological data into a systemic model and can easily be extended to describe other systems by adding new protein or nucleic acid types. The automated parameter optimization, performed with an evolutionary algorithm, reveals the sensitivity of the model to the value of each parameter and the relative importances of the experiments used. Such an analysis identifies the crucial system parameters and guides the setup of new experiments that would add most knowledge for a systemic understanding of cellular compartments. The successful combination of the molecular model and the systemic parametrization presented here on the example of the simple machinery for bacterial photosynthesis shows that it is actually possible to combine molecular and systemic modeling. This framework can now straightforwardly be applied to other currently less well characterized but biologically more relevant systems.

Computational approaches for protein function prediction: a combined strategy from multiple sequence alignment to molecular docking-based virtual screening

The functional characterization of proteins represents a daily challenge for biochemical, medical and computational sciences. Although finally proved on the bench, the function of a protein can be successfully predicted by computational approaches that drive the further experimental assays. Current methods for comparative modeling allow the construction of accurate 3D models for proteins of unknown structure, provided that a crystal structure of a homologous protein is available. Binding regions can be proposed by using binding site predictors, data inferred from homologous crystal structures, and data provided from a careful interpretation of the multiple sequence alignment of the investigated protein and its homologs. Once the location of a binding site has been proposed, chemical ligands that have a high likelihood of binding can be identified by using ligand docking and structure-based virtual screening of chemical libraries. Most docking algorithms allow building a list sorted by energy of the lowest energy docking configuration for each ligand of the library. In this review the state-of-the-art of computational approaches in 3D protein comparative modeling and in the study of protein-ligand interactions is provided. Furthermore a possible combined/concerted multistep strategy for protein function prediction, based on multiple sequence alignment, comparative modeling, binding region prediction, and structure-based virtual screening of chemical libraries, is described by using suitable examples. As practical examples, Abl-kinase molecular modeling studies, HPV-E6 protein multiple sequence alignment analysis, and some other model docking-based characterization reports are briefly described to highlight the importance of computational approaches in protein function prediction.

Modelling the molecular mechanisms of synaptic plasticity using systems biology approaches

Synaptic plasticity is thought to underlie learning and memory, but the complexity of the interactions between the ion channels, enzymes and genes that are involved in synaptic plasticity impedes a deep understanding of this phenomenon. Computer modelling has been used to investigate the information processing that is performed by the signalling pathways involved in synaptic plasticity in principal neurons of the hippocampus, striatum and cerebellum. In the past few years, new software developments that combine computational neuroscience techniques with systems biology techniques have allowed large-scale, kinetic models of the molecular mechanisms underlying long-term potentiation and long-term depression. We highlight important advancements produced by these quantitative modelling efforts and introduce promising approaches that use advancements in live-cell imaging.

A force field for virtual atom molecular mechanics of proteins

Activities of many biological macromolecules involve large conformational transitions for which crystallography can specify atomic details of alternative end states, but the course of transitions is often beyond the reach of computations based on full-atomic potential functions. We have developed a coarse-grained force field for molecular mechanics calculations based on the virtual interactions of C alpha atoms in protein molecules. This force field is parameterized based on the statistical distribution of the energy terms extracted from crystallographic data, and it is formulated to capture features dependent on secondary structure and on residue-specific contact information. The resulting force field is applied to energy minimization and normal mode analysis of several proteins. We find robust convergence in minimizations to low energies and energy gradients with low degrees of structural distortion, and atomic fluctuations calculated from the normal mode analyses correlate well with the experimental B-factors obtained from high-resolution crystal structures. These findings suggest that the virtual atom force field is a suitable tool for various molecular mechanics applications on large macromolecular systems undergoing large conformational changes.

Sending signals dynamically

Proteins mediate transmission of signals along intercellular and intracellular pathways and between the exterior and the interior of a cell. The dynamic properties of signaling proteins are crucial to their functions. We discuss emerging paradigms for the role of protein dynamics in signaling. A central tenet is that proteins fluctuate among many states on evolutionarily selected energy landscapes. Upstream signals remodel this landscape, causing signaling proteins to transmit information to downstream partners. New methods provide insight into the dynamic properties of signaling proteins at the atomic scale. The next stages in the signaling hierarchy-how multiple signals are integrated and how cellular signaling pathways are organized in space and time-present exciting challenges for the future, requiring bold multidisciplinary approaches.

The role of predictive modelling in rationally re-engineering biological systems

Technologies to synthesize and transplant a complete genome into a cell have opened limitless potential to redesign organisms for complex, specialized tasks. However, large-scale re-engineering of a biological circuit will require systems-level optimization that will come from a deep understanding of operational relationships among all the constituent parts of a cell. The integrated framework necessary for conducting such complex bioengineering requires the convergence of systems and synthetic biology. Here, we review the status of these rapidly developing interdisciplinary fields of biology and provide a perspective on plausible venues for their merger.

A framework for evolutionary systems biology

BACKGROUND

Many difficult problems in evolutionary genomics are related to mutations that have weak effects on fitness, as the consequences of mutations with large effects are often simple to predict. Current systems biology has accumulated much data on mutations with large effects and can predict the properties of knockout mutants in some systems. However experimental methods are too insensitive to observe small effects.

RESULTS

Here I propose a novel framework that brings together evolutionary theory and current systems biology approaches in order to quantify small effects of mutations and their epistatic interactions in silico. Central to this approach is the definition of fitness correlates that can be computed in some current systems biology models employing the rigorous algorithms that are at the core of much work in computational systems biology. The framework exploits synergies between the realism of such models and the need to understand real systems in evolutionary theory. This framework can address many longstanding topics in evolutionary biology by defining various 'levels' of the adaptive landscape. Addressed topics include the distribution of mutational effects on fitness, as well as the nature of advantageous mutations, epistasis and robustness. Combining corresponding parameter estimates with population genetics models raises the possibility of testing evolutionary hypotheses at a new level of realism.

CONCLUSION

EvoSysBio is expected to lead to a more detailed understanding of the fundamental principles of life by combining knowledge about well-known biological systems from several disciplines. This will benefit both evolutionary theory and current systems biology. Understanding robustness by analysing distributions of mutational effects and epistasis is pivotal for drug design, cancer research, responsible genetic engineering in synthetic biology and many other practical applications.

Fast predictions of thermodynamics and kinetics of protein-protein recognition from structures: from molecular design to systems biology

The increasing call for an overall picture of the interactions between the components of a biological system that give rise to the observed function is often summarized by the expression systems biology. Both the interpretative and predictive capabilities of holistic models of biochemical systems, however, depend to a large extent on the level of physico-chemical knowledge of the individual molecular interactions making up the network. This review is focused on the structure-based quantitative characterization of protein-protein interactions, ubiquitous in any biochemical pathway. Recently developed, fast and effective computational methods are reviewed, which allow the assessment of kinetic and thermodynamic features of the association-dissociation processes of protein complexes, both in water soluble and membrane environments. The performance and the accuracy of fast and semi-empirical structure-based methods have reached comparable levels with respect to the classical and more elegant molecular simulations. Nevertheless, the broad accessibility and lower computational cost provide the former methods with the advantageous possibility to perform systems-level analyses including extensive in silico mutagenesis screenings and large-scale structural predictions of multiprotein complexes.

Calculating enzyme kinetic parameters from protein structures

Enzyme kinetic parameters can differ between different species and isoenzymes for the same catalysed reaction. Computational approaches to calculate enzymatic kinetic parameters from the three-dimensional structures of proteins will be reviewed briefly here. Enzyme kinetic parameters may be derived by modelling and simulating the rate-determining process. An alternative, approximate, but more computationally efficient approach is the comparison of molecular interaction fields for experimentally characterized enzymes and those for which parameters should be determined. A correlation between differences in interaction fields and experimentally determined kinetic parameters can be used to determine parameters for orthologous enzymes from other species. The estimation of enzymatic kinetic parameters is an important step in setting up mathematical models of biochemical pathways in systems biology.

Molecular simulations of protein dynamics: new windows on mechanisms in biology

Recent advances in computer hardware and software have led to the development of increasingly successful molecular simulations of protein structural dynamics that are intrinsic to biological processes. These simulations have resulted in models that increasingly agree with experimental observations, suggest new experiments and provide insights into biological mechanisms. Used in combination with data obtained with sophisticated experimental techniques, simulations are helping us to understand biological complexity at the atomic and molecular levels and are giving promising insights into the genetic, thermodynamic and functional/mechanistic behaviour of biological processes. Here, we highlight some examples of such approaches that illustrate the current state and potential of the field of molecular simulation.