Predicting allosteric pockets in protein biological assemblages.

MOTIVATION: Allostery enables changes to the dynamic behavior of a protein at distant positions induced by binding. Here, we present APOP, a new allosteric pocket prediction method, which perturbs the pockets formed in the structure by stiffening pairwise interactions in the elastic network across the pocket, to emulate ligand binding. Ranking the pockets based on the shifts in the global mode frequencies, as well as their mean local hydrophobicities, leads to high prediction success when tested on a dataset of allosteric proteins, composed of both monomers and multimeric assemblages. RESULTS: Out of the 104 test cases, APOP predicts known allosteric pockets for 92 within the top 3 rank out of multiple pockets available in the protein. In addition, we demonstrate that APOP can also find new alternative allosteric pockets in proteins. Particularly interesting findings are the discovery of previously overlooked large pockets located in the centers of many protein biological assemblages; binding of ligands at these sites would likely be particularly effective in changing the protein's global dynamics. AVAILABILITY AND IMPLEMENTATION: APOP is freely available as an open-source code (https://github.com/Ambuj-UF/APOP) and as a web server at https://apop.bb.iastate.edu/.

Precise druggability of the PTH type 1 receptor.

Class B G protein-coupled receptors (GPCRs) are notoriously difficult to target by small molecules because their large orthosteric peptide-binding pocket embedded deep within the transmembrane domain limits the identification and development of nonpeptide small molecule ligands. Using the parathyroid hormone type 1 receptor (PTHR) as a prototypic class B GPCR target, and a combination of molecular dynamics simulations and elastic network model-based methods, we demonstrate that PTHR druggability can be effectively addressed. Here we found a key mechanical site that modulates the collective dynamics of the receptor and used this ensemble of PTHR conformers to identify selective small molecules with strong negative allosteric and biased properties for PTHR signaling in cell and PTH actions in vivo. This study provides a computational pipeline to detect precise druggable sites and identify allosteric modulators of PTHR signaling that could be extended to GPCRs to expedite discoveries of small molecules as novel therapeutic candidates.

ClustENMD: efficient sampling of biomolecular conformational space at atomic resolution.

SUMMARY: Efficient sampling of conformational space is essential for elucidating functional/allosteric mechanisms of proteins and generating ensembles of conformers for docking applications. However, unbiased sampling is still a challenge especially for highly flexible and/or large systems. To address this challenge, we describe a new implementation of our computationally efficient algorithm ClustENMD that is integrated with ProDy and OpenMM softwares. This hybrid method performs iterative cycles of conformer generation using elastic network model for deformations along global modes, followed by clustering and short molecular dynamics simulations. ProDy framework enables full automation and analysis of generated conformers and visualization of their distributions in the essential subspace. AVAILABILITY AND IMPLEMENTATION: ClustENMD is open-source and freely available under MIT License from https://github.com/prody/ProDy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ProDy 2.0: increased scale and scope after 10 years of protein dynamics modelling with Python.

SUMMARY: ProDy, an integrated application programming interface developed for modelling and analysing protein dynamics, has significantly evolved in recent years in response to the growing data and needs of the computational biology community. We present major developments that led to ProDy 2.0: (i) improved interfacing with databases and parsing new file formats, (ii) SignDy for signature dynamics of protein families, (iii) CryoDy for collective dynamics of supramolecular systems using cryo-EM density maps and (iv) essential site scanning analysis for identifying sites essential to modulating global dynamics. AVAILABILITY AND IMPLEMENTATION: ProDy is open-source and freely available under MIT License from https://github.com/prody/ProDy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Allosteric interactions in the parathyroid hormone GPCR-arrestin complex formation.

Peptide ligands of class B G-protein-coupled receptors act via a two-step binding process, but the essential mechanisms that link their extracellular binding to intracellular receptor-arrestin interactions are not fully understood. Using NMR, crosslinking coupled to mass spectrometry, signaling experiments and computational approaches on the parathyroid hormone (PTH) type 1 receptor (PTHR), we show that initial binding of the PTH C-terminal part constrains the conformation of the flexible PTH N-terminal signaling epitope before a second binding event occurs. A 'hot-spot' PTH residue, His9, that inserts into the PTHR transmembrane domain at this second step allosterically engages receptor-arrestin coupling. A conformational change in PTHR intracellular loop 3 permits favorable interactions with ß-arrestin's finger loop. These results unveil structural determinants for PTHR-arrestin complex formation and reveal that the two-step binding mechanism proceeds via cooperative fluctuations between ligand and receptor, which extend to other class B G-protein-coupled receptors.

Modulation of Toroidal Proteins Dynamics in Favor of Functional Mechanisms upon Ligand Binding.

Toroidal proteins serve as molecular machines and play crucial roles in biological processes such as DNA replication and RNA transcription. Despite progress in the structural characterization of several toroidal proteins, we still lack a mechanistic understanding of the significance of their architecture, oligomerization states, and intermolecular interactions in defining their biological function. In this work, we analyze the collective dynamics of toroidal proteins with different oligomerization states, namely, dimeric and trimeric DNA sliding clamps, nucleocapsid proteins (4-, 5-, and 6-mers) and Trp RNA-binding attenuation proteins (11- and 12-mers). We observe common global modes, among which cooperative rolling stands out as a mechanism enabling DNA processivity, and clamshell motions as those underlying the opening/closure of the sliding clamps. Alterations in global dynamics due to complexation with DNA or the clamp loader are shown to assist in enhancing motions to enable robust function. The analysis provides new insights into the differentiation and enhancement of functional motions upon intersubunit and intermolecular interactions.

Protein-Ligand Complexes as Constrained Dynamical Systems.

This study focuses on how the low-frequency end of the vibrational spectrum related to the functional motions changes as a protein binds to a small ligand(s). Our recently proposed residue-specific (RESPEC) elastic network model provides a natural laboratory for this aim due to its systematic mixed coarse-graining approach and parametrization. Current analysis on a large data set of protein-ligand complexes reveals a universal curve enclosing the frequency distributions, which bears the features of previous computational and experimental studies. We mostly observe positive frequency shifts in the collective modes of the protein upon ligand binding. This observation, conforming to the Rayleigh-Courant-Fisher theorem, points to a constraining effect imposed by ligands on protein dynamics, which may be accompanied by a negative vibrational entropy difference. Positive frequency shifts in the global modes can thus be linked to the harmonic well getting steeper, because of interactions with the ligand(s).

Investigation of allosteric coupling in human ß2-adrenergic receptor in the presence of intracellular loop 3.

BACKGROUND: This study investigates the allosteric coupling that exists between the intra- and extracellular parts of human ß2-adrenergic receptor (ß2-AR), in the presence of the intracellular loop 3 (ICL3), which is missing in all crystallographic experiments and most of the simulation studies reported so far. Our recent 1 µs long MD run has revealed a transition to the so-called very inactive state of the receptor, in which ICL3 packed under the G protein's binding cavity and completely blocked its accessibility to G protein. Simultaneously, an outward tilt of transmembrane helix 5 (TM5) caused an expansion of the extracellular ligand-binding site. In the current study, we performed independent runs with a total duration of 4 µs to further investigate the very inactive state with packed ICL3 and the allosteric coupling event (three unrestrained runs and five runs with bond restraints at the ligand-binding site). RESULTS: In all three independent unrestrained runs (each 500 ns long), ICL3 preserved its initially packed/closed conformation within the studied time frame, suggesting an inhibition of the receptor's activity. Specific bond restraints were later imposed between some key residues at the ligand-binding site, which have been experimentally determined to interact with the ligand. Restraining the binding site region to an open state facilitated ICL3 closure, whereas a relatively constrained/closed binding site hindered ICL3 packing. However, the reverse operation, i.e. opening of the packed ICL3, could not be realized by restraining the binding site region to a closed state. Thus, any attempt failed to free the ICL3 from its locked state due to the presence of persistent hydrogen bonds. CONCLUSIONS: Overall, our simulations indicated that starting with very inactive states, the receptor stayed almost irreversibly inhibited, which in turn decreased the overall mobility of the receptor. Bond restraints which represented the geometric restrictions caused by ligands of various sizes when bound at the ligand-binding site, induced the expected conformational changes in TM5, TM6 and consequently, ICL3. Still, once ICL3 was packed, the allosteric coupling became ineffective due to strong hydrogen bonds connecting ICL3 to the core of the receptor.

How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics.

The tunnel region at triosephosphate isomerase (TIM)'s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM's catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand's effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop's dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.

Features of large hinge-bending conformational transitions. Prediction of closed structure from open state.

We performed a detailed analysis of conformational transition pathways for a set of 10 proteins, which undergo large hinge-bending-type motions with 4-12 Å RMSD (root mean-square distance) between open and closed crystal structures. Anisotropic network model-Monte Carlo (ANM-MC) algorithm generates a targeted pathway between two conformations, where the collective modes from the ANM are used for deformation at each iteration and the conformational energy of the deformed structure is minimized via an MC algorithm. The target structure was approached successfully with an RMSD of 0.9-4.1 Å when a relatively low cutoff radius of 10 Å was used in ANM. Even though one predominant mode (first or second) directed the open-to-closed conformational transition, changes in the dominant mode character were observed for most cases along the transition. By imposing radius of gyration constraint during mode selection, it was possible to predict the closed structure for eight out of 10 proteins (with initial 4.1-7.1 Å and final 1.7-2.9 Å RMSD to target). Deforming along a single mode leads to most successful predictions. Based on the previously reported free energy surface of adenylate kinase, deformations along the first mode produced an energetically favorable path, which was interestingly facilitated by a change in mode shape (resembling second and third modes) at key points. Pathway intermediates are provided in our database of conformational transitions (http://safir.prc.boun.edu.tr/anmmc/method/1).

Influence of point mutations on PR65 conformational adaptability: Insights from molecular simulations and nanoaperture optical tweezers.

PR65 is the HEAT repeat scaffold subunit of the heterotrimeric protein phosphatase 2A (PP2A) and an archetypal tandem repeat protein. Its conformational mechanics plays a crucial role in PP2A function by opening/closing substrate binding/catalysis interface. Using in silico saturation mutagenesis, we identified PR65 "hinge" residues whose substitutions could alter its conformational adaptability and thereby PP2A function, and selected six mutations that were verified to be expressed and soluble. Molecular simulations and nanoaperture optical tweezers revealed consistent results on the specific effects of the mutations on the structure and dynamics of PR65. Two mutants observed in simulations to stabilize extended/open conformations exhibited higher corner frequencies and lower translational scattering in experiments, indicating a shift toward extended conformations, whereas another displayed the opposite features, confirmed by both simulations and experiments. The study highlights the power of single-molecule nanoaperture-based tweezers integrated with in silico approaches for exploring the effect of mutations on protein structure and dynamics.

Effect of intracellular loop 3 on intrinsic dynamics of human ß2-adrenergic receptor.

BACKGROUND: To understand the effect of the long intracellular loop 3 (ICL3) on the intrinsic dynamics of human ß2-adrenergic receptor, molecular dynamics (MD) simulations were performed on two different models, both of which were based on the inactive crystal structure in complex with carazolol (after removal of carazolol and T4-lysozyme). In the so-called loop model, the ICL3 region that is missing in available crystal structures was modeled as an unstructured loop of 32-residues length, whereas in the clipped model, the two open ends were covalently bonded to each other. The latter model without ICL3 was taken as a reference, which has also been commonly used in recent computational studies. Each model was embedded into POPC bilayer membrane with explicit water and subjected to a 1 µs molecular dynamics (MD) simulation at 310 K. RESULTS: After around 600 ns, the loop model started a transition to a "very inactive" conformation, which is characterized by a further movement of the intracellular half of transmembrane helix 6 (TM6) towards the receptor core, and a close packing of ICL3 underneath the membrane completely blocking the G-protein's binding site. Concurrently, the binding site at the extracellular part of the receptor expanded slightly with the Ser207-Asp113 distance increasing to 18 Å from 11 Å, which was further elaborated by docking studies. CONCLUSIONS: The essential dynamics analysis indicated a strong coupling between the extracellular and intracellular parts of the intact receptor, implicating a functional relevance for allosteric regulation. In contrast, no such transition to the "very inactive" state, nor any structural correlation, was observed in the clipped model without ICL3. Furthermore, elastic network analysis using different conformers for the loop model indicated a consistent picture on the specific ICL3 conformational change being driven by global modes.

Coupling between catalytic loop motions and enzyme global dynamics.

Catalytic loop motions facilitate substrate recognition and binding in many enzymes. While these motions appear to be highly flexible, their functional significance suggests that structure-encoded preferences may play a role in selecting particular mechanisms of motions. We performed an extensive study on a set of enzymes to assess whether the collective/global dynamics, as predicted by elastic network models (ENMs), facilitates or even defines the local motions undergone by functional loops. Our dataset includes a total of 117 crystal structures for ten enzymes of different sizes and oligomerization states. Each enzyme contains a specific functional/catalytic loop (10-21 residues long) that closes over the active site during catalysis. Principal component analysis (PCA) of the available crystal structures (including apo and ligand-bound forms) for each enzyme revealed the dominant conformational changes taking place in these loops upon substrate binding. These experimentally observed loop reconfigurations are shown to be predominantly driven by energetically favored modes of motion intrinsically accessible to the enzyme in the absence of its substrate. The analysis suggests that robust global modes cooperatively defined by the overall enzyme architecture also entail local components that assist in suitable opening/closure of the catalytic loop over the active site.

Cooperative mechanics of PR65 scaffold underlies the allosteric regulation of the phosphatase PP2A.

PR65, a horseshoe-shaped scaffold composed of 15 HEAT (observed in Huntingtin, elongation factor 3, protein phosphatase 2A, and the yeast kinase TOR1) repeats, forms, together with catalytic and regulatory subunits, the heterotrimeric protein phosphatase PP2A. We examined the role of PR65 in enabling PP2A enzymatic activity with computations at various levels of complexity, including hybrid approaches that combine full-atomic and elastic network models. Our study points to the high flexibility of this scaffold allowing for end-to-end distance fluctuations of 40-50 Å between compact and extended conformations. Notably, the intrinsic dynamics of PR65 facilitates complexation with the catalytic subunit and is retained in the PP2A complex enabling PR65 to engage the two domains of the catalytic subunit and provide the mechanical framework for enzymatic activity, with support from the regulatory subunit. In particular, the intra-repeat coils at the C-terminal arm play an important role in allosterically mediating the collective dynamics of PP2A, pointing to target sites for modulating PR65 function.

Influence of Point Mutations on PR65 Conformational Adaptability: Insights from Nanoaperture Optical Tweezer Experiments and Molecular Simulations.

PR65 is the HEAT-repeat scaffold subunit of the heterotrimeric protein phosphatase 2A (PP2A) and an archetypal tandem-repeat protein, forming a spring-like architecture. PR65 conformational mechanics play a crucial role in PP2A function by opening/closing the substrate-binding/catalysis interface. Using in-silico saturation mutagenesis we identified "hinge" residues of PR65, whose substitutions are predicted to restrict its conformational adaptability and thereby disrupt PP2A function. Molecular simulations revealed that a subset of hinge mutations stabilized the extended/open conformation, whereas another had the opposite effect. By trapping in nanoaperture optical tweezer, we characterized PR65 motion and showed that the former mutants exhibited higher corner frequencies and lower translational scattering, indicating a shift towards extended conformations, whereas the latter showed the opposite behavior. Thus, experiments confirm the conformations predicted computationally. The study highlights the utility of nanoaperture-based tweezers for exploring structure and dynamics, and the power of integrating this single-molecule method with in silico approaches.

Activation and Speciation Mechanisms in Class A GPCRs.

Accurate development of allosteric modulators of GPCRs require a thorough assessment of their sequence, structure, and dynamics, toward gaining insights into their mechanisms of actions shared by family members, as well as dynamic features that distinguish subfamilies. Building on recent progress in the characterization of the signature dynamics of proteins, we analyzed here a dataset of 160 Class A GPCRs to determine their sequence similarities, structural landscape, and dynamic features across different species (human, bovine, mouse, squid, and rat), different activation states (active/inactive), and different subfamilies. The two dominant directions of variability across experimentally resolved structures, identified by principal component analysis of the dataset, shed light to cooperative mechanisms of activation, subfamily differentiation, and speciation of Class A GPCRs. The analysis reveals the functional significance of the conformational flexibilities of specific structural elements, including: the dominant role of the intracellular loop 3 (ICL3) together with the cytoplasmic ends of the adjoining helices TM5 and TM6 in enabling allosteric activation; the role of particular structural motifs at the extracellular loop 2 (ECL2) connecting TM4 and TM5 in binding ligands specific to different subfamilies; or even the differentiation of the N-terminal conformation across different species. Detailed analyses of the modes of motions accessible to the members of the dataset and their variations across members demonstrate how the active and inactive states of GPCRs obey distinct conformational dynamics. The collective fluctuations of the GPCRs are robustly defined in the active state, while the inactive conformers exhibit broad variance among members.

GPCR kinases generate an APH1A phosphorylation barcode to regulate amyloid-ß generation.

Emerging evidence suggests that G protein-coupled receptor (GPCR) kinases (GRKs) are associated with the pathophysiology of Alzheimer's disease (AD). However, GRKs have not been directly implicated in regulation of the amyloid-ß (Aß) pathogenic cascade in AD. Here, we determine that GRKs phosphorylate a non-canonical substrate, anterior pharynx-defective 1A (APH1A), an integral component of the γ-secretase complex. Significantly, we show that GRKs generate distinct phosphorylation barcodes in intracellular loop 2 (ICL2) and the C terminus of APH1A, which differentially regulate recruitment of the scaffolding protein ß-arrestin 2 (ßarr2) to APH1A and γ-secretase-mediated Aß generation. Further molecular dynamics simulation studies reveal an interaction between the ßarr2 finger loop domain and ICL2 and ICL3 of APH1A, similar to a GPCR-ß-arrestin complex, which regulates γ-secretase activity. Collectively, these studies provide insight into the molecular and structural determinants of the APH1A-ßarr2 interaction that critically regulate Aß generation.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods.

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

Effect of ligand binding on the intraminimum dynamics of proteins.

Effects of ligand binding on protein dynamics are studied via molecular dynamics (MD) simulations on two different enzymes, dihydrofolate reductase (DHFR) and triosephosphate isomerase (TIM), in their unliganded (free) and liganded states. Domain motions in MD trajectories are analyzed by collectivities and rotation angles along the principal components (PCs). DHFR in the free state has well-defined domain rotations, whereas rotations are slightly damped in the binary complex with nicotinamide adenine dinucleotide phosphate (NADPH), and remarkably distorted in the presence of NADP(+) , showing that NADP(+) is solely responsible for the loss of correlation of the domains in DHFR. Although mean square fluctuations of MD simulations in the same PC subspaces are similar for different ligation states, linear stochastic time series models show that backbone flexibility along the first five PCs is decreased upon NADPH and NADP(+) binding in subpicosecond scale. This shows that mobility of the protein along the PCs is closely related with intraminimum dynamics, and alterations in ligation states may change the intraminimum dynamics significantly. Low vibrational frequencies of the alpha-carbon atoms of DHFR are determined from the time series models of a larger number of low indexed PCs, and it is found that number of modes in the lowest frequencies is reduced upon ligand binding. A similar result is obtained for TIM in the unliganded and dihydroxyacetone phosphate bound states. We suggest that stochastic time series modeling is a promising method to be used in determining subtle perturbations in protein dynamics.

Essential site scanning analysis: A new approach for detecting sites that modulate the dispersion of protein global motions.

Despite the wealth of methods developed for exploring the molecular basis of allostery in biomolecular systems, there is still a need for structure-based predictive tools that can efficiently detect susceptible sites for triggering allosteric responses. Toward this goal, we introduce here an elastic network model (ENM)-based method, Essential Site Scanning Analysis (ESSA). Essential sites are here defined as residues that would significantly alter the protein's global dynamics if bound to a ligand. To mimic the crowding induced upon substrate binding, the heavy atoms of each residue are incorporated as additional network nodes into the α-carbon-based ENM, and the resulting shifts in soft mode frequencies are used as a metric for evaluating the essentiality of each residue. Results on a dataset of monomeric proteins indicate the enrichment of allosteric and orthosteric binding sites, as well as global hinge regions among essential residues, highlighting the significant role of these sites in controlling the overall structural dynamics. Further integration of ESSA with information on predicted pockets and their local hydrophobicity density enables successful predictions of allosteric pockets for both ligand-bound and -unbound structures. ESSA can be efficiently applied to large multimeric systems. Three case studies, namely (i) G-protein binding to a GPCR, (ii) heterotrimeric assembly of the Ser/Thr protein phosphatase PP2A, and (iii) allo-targeting of AMPA receptor, demonstrate the utility of ESSA for identifying essential sites and narrowing down target allosteric sites identified by druggability simulations.