Computational design and characterization of nanobody-derived peptides that stabilize the active conformation of the β2-adrenergic receptor (β2-AR)

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 in apo. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganizations in key salt bridge and hydrogen bond networks contributing to the CB1 activation/inactivation. For instance, D213-Y224 hydrogen bond and D184-K192 salt bridge showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

In Silico Exploration of CD200 as a Therapeutic Target for COVID-19

SARS-CoV-2, the pathogen causing COVID-19, continues to pose a significant threat to public health and has had major economic implications. Developing safe and effective vaccines and therapies offers a path forward for overcoming the COVID-19 pandemic. The presented study, performed by using the informational spectrum method (ISM), representing an electronic biology-based tool for analysis of protein-protein interactions, identified the highly conserved region of spike protein (SP) from SARS-CoV-2 virus, which is essential for recognition and targeting between the virus and its protein interactors on the target cells. This domain is suggested as a promising target for the drug therapy and vaccines, which could be effective against all currently circulating variants of SARS-CoV-2 viruses. The analysis of the virus/host interaction, performed by the ISM, also revealed OX-2 membrane glycoprotein (CD200) as a possible interactor of SP, which could serve as a novel therapeutic target for COVID-19 disease.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

The cannabinoid receptor CB1 is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 states in apo conditions. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganization of key salt bridge and hydrogen bond networks known to control CB1 activation between states. For instance, a conserved D213-Y224 hydrogen bond and D184-K192 salt bridge interactions showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

Computational Modeling and Characterization of Peptides Derived from Nanobody Complementary-Determining Region 2 (CDR2) Targeting Active-State Conformation of the β2-Adrenergic Receptor (β2AR)

This study assessed the suitability of the complementarity-determining region 2 (CDR2) of the nanobody (Nb) as a template for the derivation of nanobody-derived peptides (NDPs) targeting active-state β2-adrenergic receptor (β2AR) conformation. Sequences of conformationally selective Nbs favoring the agonist-occupied β2AR were initially analyzed by the informational spectrum method (ISM). The derived NDPs in complex with β2AR were subjected to protein-peptide docking, molecular dynamics (MD) simulations, and metadynamics-based free-energy binding calculations. Computational analyses identified a 25-amino-acid-long CDR2-NDP of Nb71, designated P4, which exhibited the following binding free-energy for the formation of the β2AR:P4 complex (ΔG = -6.8 ± 0.8 kcal/mol or a Ki = 16.5 μM at 310 K) and mapped the β2AR:P4 amino acid interaction network. In vitro characterization showed that P4 (i) can cross the plasma membrane, (ii) reduces the maximum isoproterenol-induced cAMP level by approximately 40% and the isoproterenol potency by up to 20-fold at micromolar concentration, (iii) has a very low affinity to interact with unstimulated β2AR in the cAMP assay, and (iv) cannot reduce the efficacy and potency of the isoproterenol-mediated β2AR/β-arrestin-2 interaction in the BRET2-based recruitment assay. In summary, the CDR2-NDP, P4, binds preferentially to agonist-activated β2AR and disrupts Gαs-mediated signaling.

Intracellular VHHs to monitor and modulate GPCR signaling

Single-domain antibody fragments, also known as VHHs or nanobodies, have opened promising avenues in therapeutics and in exploration of intracellular processes. Because of their unique structural properties, they can reach cryptic regions in their cognate antigen. Intracellular VHHs/antibodies primarily directed against cytosolic proteins or transcription factors have been described. In contrast, few of them target membrane proteins and even less recognize G protein-coupled receptors. These receptors are major therapeutic targets, which reflects their involvement in a plethora of physiological responses. Hence, they elicit a tremendous interest in the scientific community and in the industry. Comprehension of their pharmacology has been obscured by their conformational complexity, that has precluded deciphering their structural properties until the early 2010's. To that respect, intracellular VHHs have been instrumental in stabilizing G protein-coupled receptors in active conformations in order to solve their structure, possibly bound to their primary transducers, G proteins or β-arrestins. In contrast, the modulatory properties of VHHs recognizing the intracellular regions of G protein-coupled receptors on the induced signaling network have been poorly studied. In this review, we will present the advances that the intracellular VHHs have permitted in the field of GPCR signaling and trafficking. We will also discuss the methodological hurdles that linger the discovery of modulatory intracellular VHHs directed against GPCRs, as well as the opportunities they open in drug discovery.

Microscale thermophoresis as a powerful growing analytical technique for the investigation of biomolecular interaction and the determination of binding parameters

The in vitro panel of technologies to address biomolecular interactions are in play, however microscale thermophoresis is continuously increasing in use to represent a key player in this arena. This review highlights the usefulness of microscale thermophoresis in the determination of molecular and biomolecular affinity interactions. This work reviews the literature from January 2016 to January 2022 about microscale thermophoresis. It gives a summarized overview about both the state-of the art and the development in the field of microscale thermophoresis. The principle of microscale thermophoresis is also described supported with self-created illustrations. Moreover, some recent advances are mentioned that showing application of the technique in investigating biomolecular interactions in different fields. Finally, advantages as well as drawbacks of the technique in comparison with other competing techniques are summarized.

Downsizing antibodies: Towards complementarity-determining region (CDR)-based peptide mimetics

Monoclonal antibodies emerged as an important therapeutic drug class with remarkable specificity and binding affinity. Nonetheless, these heterotetrameric immunoglobulin proteins come with high manufacturing and therapeutic costs which can take extraordinary proportions, besides other limitations such as their limited in cellulo access imposed by their molecular size (ca. 150 kDa). These drawbacks stimulated the development of downsized functional antibody fragments (ca. 15-50 kDa), together with smaller synthetic peptides (ca. 1-3 kDa) derived from the antibodies' crucial complementarity-determining regions (CDR). Despite the general lack of success in the literal translation of CDR loops in peptide mimetics, rational structure-based and computational approaches have shown their potential for obtaining functional CDR-based peptide mimetics. In this review, we describe the efforts made in the development of antibody and nanobody paratope-derived peptide mimetics with particular focus on the used design strategies, in addition to highlighting the challenges associated with their development.

Antibody VH domain sequence analysis by a bioinformatics approach based on electronic amino acid properties may help to predict paratop location

Gene as the basic functional unit of DNA encodes information about the product such as protein. The majority of proteins realize function through protein-protein interactions involving short protein motifs. However, some proteins such as antibodies are established by the rearrangement of several (V-D-J) gene segments with the potential addition of nontemplated nucleotides that may change information encoded by the respective gene segment used. Antibody VH domain sequence analysis by ISM bioinformatics approach that is based on amino acids physicochemical features, enable to distinguish the contribution of the information encoded by VH gene or generated during VDJ gene recombination for antibody-antigen interaction. The data presented in this report revealed the significance of CDRH3 for the interaction of antibody specific for immunogenic molecules while CDRH3 contribution is minor for antibody interaction with nonimmunogenic molecules such as haptens and native mammalian dsDNA. Thus, paratopes might be located in the CDRH3 or VH regions.

A new class of sensing elements for sensors: Clamp peptides for Zika virus

The design of a new class of selective and high affinity antibody mimetics termed clamp peptide (CP) that incorporate three short peptides structurally and mechanically mimicking a clamp is proposed as sensing elements for a reliable detection sensor platform. The CPs consist of two short peptides functioning as arms that recognize two different epitopes in the target protein and are connected by a third short peptide that acts as a hinge between the peptide arms. For the construction of CPs, we employed a rational design combined with computational methods. To illustrate our approach, we designed a CP that binds selectively to the envelope protein of the Zika virus (ZIKV). The virtual docking cycles were run maximizing the discrimination between ZIKV and Dengue virus (DENV) envelope proteins. DENV was chosen among the flavivirus family because it has high structural similarity with ZIKV. When employed in a colorimetric binding assay or in label-free electrochemical impedance sensor format, the CP was selective for ZIKV vs DENV particles showing detection limit under 10⁴ copies/mL, comparable to anti-ZIKV antibodies. Apparent dissociation binding constants (Kd) confirmed a better performance of CPs than mono-arm peptides (Kd of best CP = 162 nM ± 23 nM; Kd of best mono-arm peptide = 11.15 ± 2.76 μM). The performance of the assays based on CPs was also verified in serum and urine (diluted 1:10 and 1:1 respectively). The detection limits of CPs decreased about one order of magnitude for ZIKV detection in serum or urine, with a distinct analytical signal starting from 10⁵ copies/mL of ZIKV.

Identification of ligand-specific G protein-coupled receptor states and prediction of downstream efficacy via data-driven modeling

Ligand binding stabilizes different G protein-coupled receptor states via a complex allosteric process that is not completely understood. Here, we have derived free energy landscapes describing activation of the β2 adrenergic receptor bound to ligands with different efficacy profiles using enhanced sampling molecular dynamics simulations. These reveal shifts toward active-like states at the Gprotein-binding site for receptors bound to partial and full agonists, and that the ligands modulate the conformational ensemble of the receptor by tuning protein microswitches. We indeed find an excellent correlation between the conformation of the microswitches close to the ligand binding site and in the transmembrane region and experimentally reported cyclic adenosine monophosphate signaling responses. Dimensionality reduction further reveals the similarity between the unique conformational states induced by different ligands, and examining the output of classifiers highlights two distant hotspots governing agonism on transmembrane helices 5 and 7.

Dimeric Transmembrane Orientations of APP/C99 Regulate γ-Secretase Processing Line Impacting Signaling and Oligomerization

Amyloid precursor protein (APP) cleavage by the β-secretase produces the C99 transmembrane (TM) protein, which contains three dimerization-inducing Gly-x-x-x-Gly motifs. We demonstrate that dimeric C99 TM orientations regulate the precise cleavage lines by γ-secretase. Of all possible dimeric orientations imposed by a coiled-coil to the C99 TM domain, the dimer containing the ³³Gly-x-x-x-Gly³⁷ motif in the interface promoted the Aβ₄₂ processing line and APP intracellular domain-dependent gene transcription, including the induction of BACE1 mRNA, enhancing amyloidogenic processing and signaling. Another orientation exhibiting the ²⁵Gly-x-x-x-Gly²⁹ motif in the interface favored processing to Aβ_43/40. It induced significantly less gene transcription, while promoting formation of SDS-resistant "Aβ-like" oligomers, reminiscent of Aβ peptide oligomers. These required both Val24 of a pro-β motif and the ²⁵Gly-x-x-x-Gly²⁹ interface. Thus, crossing angles imposed by precise dimeric orientations control γ-secretase initial cleavage at Aβ₄₈ or Aβ_49, linking the former to enhanced signaling and Aβ₄₂ production.

Critical Impact of Different Conserved Endoplasmic Retention Motifs and Dopamine Receptor Interacting Proteins (DRIPs) on Intracellular Localization and Trafficking of the D2 Dopamine Receptor (D2-R) Isoforms

The type 2 dopamine receptor D₂ (D₂-R), member of the G protein-coupled receptor (GPCR) superfamily, exists in two isoforms, short (D_2S-R) and long (D_2L-R). They differ by an additional 29 amino acids (AA) in the third cytoplasmic loop (ICL3) of the D_2L-R. These isoforms differ in their intracellular localization and trafficking functionality, as D_2L-R possesses a larger intracellular pool, mostly in the endoplasmic reticulum (ER). This review focuses on the evolutionarily conserved motifs in the ICL3 of the D₂-R and proteins interacting with the ICL3 of both isoforms, specifically with the 29 AA insert. These motifs might be involved in D₂-R exit from the ER and have an impact on cell-surface and intracellular localization and, therefore, also play a role in the function of dopamine receptor signaling, ligand binding and possible homo/heterodimerization. Our recent bioinformatic data on potential new interaction partners for the ICL3 of D₂-Rs are also presented. Both are highly relevant, and have clinical impacts on the pathophysiology of several diseases such as Parkinson’s disease, schizophrenia, Tourette’s syndrome, Huntington’s disease, manic depression, and others, as they are connected to a variety of essential motifs and differences in communication with interaction partners.

Docking of peptides to GPCRs using a combination of CABS-dock with FlexPepDock refinement

The structural description of peptide ligands bound to G protein-coupled receptors (GPCRs) is important for the discovery of new drugs and deeper understanding of the molecular mechanisms of life. Here we describe a three-stage protocol for the molecular docking of peptides to GPCRs using a set of different programs: (1) CABS-dock for docking fully flexible peptides; (2) PD2 method for the reconstruction of atomistic structures from C-alpha traces provided by CABS-dock and (3) Rosetta FlexPepDock for the refinement of protein–peptide complex structures and model scoring. We evaluated the proposed protocol on the set of seven different GPCR–peptide complexes (including one containing a cyclic peptide), for which crystallographic structures are available. We show that CABS-dock produces high resolution models in the sets of top-scored models. These sets of models, after reconstruction to all-atom representation, can be further improved by Rosetta high-resolution refinement and/or minimization, leading in most of the cases to sub-Angstrom accuracy in terms of interface root-mean-square-deviation measure.