Elastic network normal mode dynamics reveal the GPCR activation mechanism

Computational and experimental approaches to probe GPCR activation and signaling

G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales.

Nucleic Acid Aptamers Emerging as Modulators of G-Protein-Coupled Receptors: Challenge to Difficult Cell Surface Proteins

G-protein-coupled receptors (GPCRs), among various cell surface proteins, are essential targets in the fields of basic science and drug discovery. The discovery and development of modulators for the receptors have provided deep insights into the mechanism of action of receptors and have led to a new therapeutic option for human diseases. Although various modulators against GPCRs have been developed to date, the identification of new modulators for GPCRs remains a challenge due to several technical problems and limitations. To overcome this situation, a variety of strategies have been developed by several modalities, including nucleic acid aptamers, which are emerging as unique molecules isolated by a repetitive selection process against various types of targets from an enormous combinatorial library. This review summarized the achievements in the development of aptamers targeting GPCRs, and discussed their isolation methods and the diverse functional features of aptamers against GPCRs.

Using normal mode analysis on protein structural models. How far can we go on our predictions?

Normal mode analysis (NMA) is a fast and inexpensive approach that is largely used to gain insight into functional protein motions, and more recently to create conformations for further computational studies. However, when the protein structure is unknown, the use of computational models is necessary. Here, we analyze the capacity of NMA in internal coordinate space to predict protein motion, its intrinsic flexibility, and atomic displacements, using protein models instead of native structures, and the possibility to use it for model refinement. Our results show that NMA is quite insensitive to modeling errors, but that calculations are strictly reliable only for very accurate models. Our study also suggests that internal NMA is a more suitable tool for the improvement of structural models, and for integrating them with experimental data or in other computational techniques, such as protein docking or more refined molecular dynamics simulations.

Virtual screening for potential inhibitors of Mcl-1 conformations sampled by normal modes, molecular dynamics, and nuclear magnetic resonance

Myeloid cell leukemia-1 (Mcl-1) is often overexpressed in human cancer and is an important target for developing antineoplastic drugs. In this study, a data set containing 2.3 million lead-like molecules and a data set of all the US Food and Drug Administration (FDA)-approved drugs are virtually screened for potential Mcl-1 ligands using Protein Data Bank (PDB) ID 2MHS. The potential Mcl-1 ligands are evaluated and computationally docked on to three conformation ensembles generated by normal mode analysis (NMA), molecular dynamics (MD), and nuclear magnetic resonance (NMR), respectively. The evaluated potential Mcl-1 ligands are then compared with their clinical use. Remarkably, half of the top 30 potential drugs are used clinically to treat cancer, thus partially validating our virtual screen. The partial validation also favors the idea that the other half of the top 30 potential drugs could be used in the treatment of cancer. The normal mode-, MD-, and NMR-based conformation greatly expand the conformational sampling used herein for in silico identification of potential Mcl-1 inhibitors.

Recent Progress in Understanding the Conformational Mechanism of Heterotrimeric G Protein Activation

Heterotrimeric G proteins are key intracellular coordinators that receive signals from cells through activation of cognate G protein-coupled receptors (GPCRs). The details of their atomic interactions and structural mechanisms have been described by many biochemical and biophysical studies. Specifically, a framework for understanding conformational changes in the receptor upon ligand binding and associated G protein activation was provided by description of the crystal structure of the β2-adrenoceptor-Gs complex in 2011. This review focused on recent findings in the conformational dynamics of G proteins and GPCRs during activation processes.

Similarity in Shape Dictates Signature Intrinsic Dynamics Despite No Functional Conservation in TIM Barrel Enzymes

The conservation of the intrinsic dynamics of proteins emerges as we attempt to understand the relationship between sequence, structure and functional conservation. We characterise the conservation of such dynamics in a case where the structure is conserved but function differs greatly. The triosephosphate isomerase barrel fold (TBF), renowned for its 8 β-strand-α-helix repeats that close to form a barrel, is one of the most diverse and abundant folds found in known protein structures. Proteins with this fold have diverse enzymatic functions spanning five of six Enzyme Commission classes, and we have picked five different superfamily candidates for our analysis using elastic network models. We find that the overall shape is a large determinant in the similarity of the intrinsic dynamics, regardless of function. In particular, the β-barrel core is highly rigid, while the α-helices that flank the β-strands have greater relative mobility, allowing for the many possibilities for placement of catalytic residues. We find that these elements correlate with each other via the loops that link them, as opposed to being directly correlated. We are also able to analyse the types of motions encoded by the normal mode vectors of the α-helices. We suggest that the global conservation of the intrinsic dynamics in the TBF contributes greatly to its success as an enzymatic scaffold both through evolution and enzyme design.

Computational approaches to detect allosteric pathways in transmembrane molecular machines

Many of the functions of transmembrane proteins involved in signal processing and transduction across the cell membrane are determined by allosteric couplings that propagate the functional effects well beyond the original site of activation. Data gathered from breakthroughs in biochemistry, crystallography, and single molecule fluorescence have established a rich basis of information for the study of molecular mechanisms in the allosteric couplings of such transmembrane proteins. The mechanistic details of these couplings, many of which have therapeutic implications, however, have only become accessible in synergy with molecular modeling and simulations. Here, we review some recent computational approaches that analyze allosteric coupling networks (ACNs) in transmembrane proteins, and in particular the recently developed Protein Interaction Analyzer (PIA) designed to study ACNs in the structural ensembles sampled by molecular dynamics simulations. The power of these computational approaches in interrogating the functional mechanisms of transmembrane proteins is illustrated with selected examples of recent experimental and computational studies pursued synergistically in the investigation of secondary active transporters and GPCRs. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Use of network model to explore dynamic and allosteric properties of three GPCR homodimers

We used an elastic network model and protein structure network to study three class A GPCR homodimers.

Adaptability of protein structures to enable functional interactions and evolutionary implications

Several studies in recent years have drawn attention to the ability of proteins to adapt to intermolecular interactions by conformational changes along structure-encoded collective modes of motions. These so-called soft modes, primarily driven by entropic effects, facilitate, if not enable, functional interactions. They represent excursions on the conformational space along principal low-ascent directions/paths away from the original free energy minimum, and they are accessible to the protein even before protein-protein/ligand interactions. An emerging concept from these studies is the evolution of structures or modular domains to favor such modes of motion that will be recruited or integrated for enabling functional interactions. Structural dynamics, including the allosteric switches in conformation that are often stabilized upon formation of complexes and multimeric assemblies, emerge as key properties that are evolutionarily maintained to accomplish biological activities, consistent with the paradigm sequence→structure→dynamics→function where 'dynamics' bridges structure and function.

WEBnm@ v2.0: Web server and services for comparing protein flexibility

BACKGROUND

Normal mode analysis (NMA) using elastic network models is a reliable and cost-effective computational method to characterise protein flexibility and by extension, their dynamics. Further insight into the dynamics-function relationship can be gained by comparing protein motions between protein homologs and functional classifications. This can be achieved by comparing normal modes obtained from sets of evolutionary related proteins.

RESULTS

We have developed an automated tool for comparative NMA of a set of pre-aligned protein structures. The user can submit a sequence alignment in the FASTA format and the corresponding coordinate files in the Protein Data Bank (PDB) format. The computed normalised squared atomic fluctuations and atomic deformation energies of the submitted structures can be easily compared on graphs provided by the web user interface. The web server provides pairwise comparison of the dynamics of all proteins included in the submitted set using two measures: the Root Mean Squared Inner Product and the Bhattacharyya Coefficient. The Comparative Analysis has been implemented on our web server for NMA, WEBnm@, which also provides recently upgraded functionality for NMA of single protein structures. This includes new visualisations of protein motion, visualisation of inter-residue correlations and the analysis of conformational change using the overlap analysis. In addition, programmatic access to WEBnm@ is now available through a SOAP-based web service. Webnm@ is available at http://apps.cbu.uib.no/webnma .

CONCLUSION

WEBnm@ v2.0 is an online tool offering unique capability for comparative NMA on multiple protein structures. Along with a convenient web interface, powerful computing resources, and several methods for mode analyses, WEBnm@ facilitates the assessment of protein flexibility within protein families and superfamilies. These analyses can give a good view of how the structures move and how the flexibility is conserved over the different structures.

Microsecond Molecular Dynamics Simulations Provide Insight into the Allosteric Mechanism of the Gs Protein Uncoupling from the β2 Adrenergic Receptor

Experiments have revealed that in the β(2) adrenergic receptor (β(2)AR)-Gs protein complex the α subunit (Gαs) of the Gs protein can adopt either an "open" conformation or a "closed" conformation. In the "open" conformation the Gs protein prefers to bind to the β(2)AR, while in the "closed" conformation an uncoupling of the Gs protein from the β(2)AR occurs. However, the mechanism that leads to such different behaviors of the Gs protein remains unclear. Here, we report results from microsecond molecular dynamics simulations and community network analysis of the β(2)AR-Gs complex with Gαs in the "open" and "closed" conformations. We observed that the complex is stabilized differently in the "open" and "closed" conformations. The community network analysis reveals that in the "closed" conformation there exists strong allosteric communication between the β(2)AR and Gβγ, mediated by Gαs. We suggest that such high information flows are necessary for the Gs protein uncoupling from the β(2)AR.

Comparing the intrinsic dynamics of multiple protein structures using elastic network models

BACKGROUND

Elastic network models (ENMs) are based on the simple idea that a protein can be described as a set of particles connected by springs, which can then be used to describe its intrinsic flexibility using, for example, normal mode analysis. Since the introduction of the first ENM by Monique Tirion in 1996, several variants using coarser protein models have been proposed and their reliability for the description of protein intrinsic dynamics has been widely demonstrated. Lately an increasing number of studies have focused on the meaning of slow dynamics for protein function and its potential conservation through evolution. This leads naturally to comparisons of the intrinsic dynamics of multiple protein structures with varying levels of similarity.

SCOPE OF REVIEW

We describe computational strategies for calculating and comparing intrinsic dynamics of multiple proteins using elastic network models, as well as a selection of examples from the recent literature.

MAJOR CONCLUSIONS

The increasing interest for comparing dynamics across protein structures with various levels of similarity, has led to the establishment and validation of reliable computational strategies using ENMs. Comparing dynamics has been shown to be a viable way for gaining greater understanding for the mechanisms employed by proteins for their function. Choices of ENM parameters, structure alignment or similarity measures will likely influence the interpretation of the comparative analysis of protein motion.

GENERAL SIGNIFICANCE

Understanding the relation between protein function and dynamics is relevant to the fundamental understanding of protein structure-dynamics-function relationship. This article is part of a Special Issue entitled Recent developments of molecular dynamics.