DIRECT‐ID: An automated method to identify and quantify conformational variations—application to β2‐adrenergic GPCR

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT‐ID, was used to identify relevant ligand‐modulated structural variations in the β₂‐adrenergic (β₂AR) G‐protein coupled receptor. Micro‐second MD simulations of the β₂AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI‐167107 (agonist), epinephrine (agonist), salbutamol (long‐acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β₂AR differently and DIRECT‐ID analysis of the inverse‐agonist vs. agonist‐modulated β₂AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand‐specific microswitches across the GPCR. This study demonstrates the utility of DIRECT‐ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems. © 2015 Wiley Periodicals, Inc.     

Identification of a Conformational Equilibrium That Determines the Efficacy and Functional Selectivity of the μ‐Opioid Receptor

G‐protein‐coupled receptor (GPCR) ligands impart differing degrees of signaling in the G‐protein and arrestin pathways, in phenomena called “biased signaling”. However, the mechanism underlying the biased signaling of GPCRs is still unclear, although crystal structures of GPCRs bound to the G protein or arrestin are available. In this study, we observed the NMR signals from methionine residues of the μ‐opioid receptor (μOR) in the balanced‐ and biased‐ligand‐bound states. We found that the intracellular cavity of μOR exists in an equilibrium between closed and multiple open conformations with coupled conformational changes on the transmembrane helices 3, 5, 6, and 7, and that the population of each open conformation determines the G‐protein‐ and arrestin‐mediated signaling levels in each ligand‐bound state. These findings provide insight into the biased signaling of GPCRs and will be helpful for development of analgesics that stimulate μOR with reduced tolerance and dependence.     

Isotopic splitting patterns in the 13C NMR spectra of some partially deuterated 1‐aryl‐2‐(phenyldiazenyl)butane‐1,3‐dione and 4‐hydroxy‐3‐(phenyldiazenyl)‐2H‐chromen‐2‐one: evidence for elucidation of tautomeric forms

Nuclear magnetic resonance spectra of synthesized azo dyes derived from aniline derivatives in reaction with benzoylacetone and 4‐hydroxycoumarin were studied in both CDCl₃ and (CD₃)₂SO (two drops of D₂O were added into solutions of dyes). All dyes showed intramolecular hydrogen bonding. Dyes derived from o‐nitro aniline in the reaction with benzoylacetone, and 4‐hydroxycoumarin showed bifurcated intramolecular hydrogen bonds. The solvent‐substrate proton exchange of dyes derived from benzoylacetone and 4‐hydroxycoumarin was examined in the presence of two drops of D₂O. Among ten dye samples, two dyes derived from benzoylacetone did not show deuteration, three dyes showed partial deuteration and five dyes showed full deuteration under similar conditions. For the partially deuterated dyes the β‐isotope effect in ¹³C splitting was investigated and was used for the determination of the predominant tautomeric form. Copyright © 2016 John Wiley & Sons, Ltd.     

Nanobody‐Enabled Reverse Pharmacology on G‐Protein‐Coupled Receptors

The conformational complexity of transmembrane signaling of G‐protein‐coupled receptors (GPCRs) is a central hurdle for the design of screens for receptor agonists. In their basal states, GPCRs have lower affinities for agonists compared to their G‐protein‐bound active state conformations. Moreover, different agonists can stabilize distinct active receptor conformations and do not uniformly activate all cellular signaling pathways linked to a given receptor (agonist bias). Comparative fragment screens were performed on a β₂‐adrenoreceptor–nanobody fusion locked in its active‐state conformation by a G‐protein‐mimicking nanobody, and the same receptor in its basal‐state conformation. This simple biophysical assay allowed the identification and ranking of multiple novel agonists and permitted classification of the efficacy of each hit in agonist, antagonist, or inverse agonist categories, thereby opening doors to nanobody‐enabled reverse pharmacology.     

In‐Membrane Chemical Modification (IMCM) for Site‐Specific Chromophore Labeling of GPCRs

We present in‐membrane chemical modification (IMCM) for obtaining selective chromophore labeling of intracellular surface cysteines in G‐protein‐coupled receptors (GPCRs) with minimal mutagenesis. This method takes advantage of the natural protection of most cysteines by the membrane environment. Practical use of IMCM is illustrated with the site‐specific introduction of chromophores for NMR and fluorescence spectroscopy in the human κ‐opioid receptor (KOR) and the human A_2A adenosine receptor (A_2AAR). IMCM is applicable to a wide range of in vitro studies of GPCRs, including single‐molecule spectroscopy, and is a promising platform for in‐cell spectroscopy experiments.     

Using open‐tubular capillary electrochromatography with part‐coating column for binding constants determination of β2‐adrenergic receptor with seven drugs

An open‐tubular capillary electrochromatography method has been developed for the determination of binding constants between β₂‐adrenergic receptor (β₂‐AR) and seven drugs. β₂‐AR was oriented immobilized onto one part of inner surface of capillary via microwave‐assisted technical synthesis. According to the linear relationship between coating length and the apparent mobility of analyte, the binding constant (K_b) can be obtained by related theories and equations. The order of K_b values between drugs such as adrenaline hydrochloride, norepinephrine bitartrate, and propranolol hydrochloride with β₂‐AR is well consistent with that reported in the literature. By the method, K_b values between four extracts of Radix Paeoniae Rubra and β₂‐AR were also successfully obtained. Subsequently, computer models were applied to interpret the CEC experiments. And the results proved to be in good agreement with the method. The work, herein, demonstrates the potential of the method in drug‐receptor affinity interactions evaluation and screening of lead compounds from natural sources.     

High‐expression β1 adrenergic receptor/cell membrane chromatography method based on a target receptor to screen active ingredients from traditional Chinese medicines

β‐Adrenergic receptors are important targets for drug discovery. We have developed a new β₁‐adrenergic receptor cell membrane chromatography (β₁AR‐CMC) with offline ultra‐performance LC (UPLC) and MS method for screening active ingredients from traditional Chinese medicines. In this study, Chinese hamster ovary‐S cells with high β₁AR expression levels were established and used to prepare a cell membrane stationary phase in a β₁AR‐CMC model. The retention fractions were separated and identified by the UPLC–MS system. The screening results found that isoimperatorin from Rhizoma et Radix Notopterygii was the targeted component that could act on β₁AR in similar manner of metoprolol as a control drug. In addition, the biological effects of active component were also investigated in order to search for a new type of β₁AR antagonist. It will be a useful method for drug discovery as a leading compound resource.     

G-Protein-Coupled Receptor–Membrane Interactions Depend on the Receptor Activation State

G-protein-coupled receptors (GPCRs) are the largest family of human membrane proteins and serve as primary targets of approximately one-third of currently marketed drugs. In particular, adenosine A₁ receptor (A₁AR) is an important therapeutic target for treating cardiac ischemia–reperfusion injuries, neuropathic pain, and renal diseases. As a prototypical GPCR, the A₁AR is located within a phospholipid membrane bilayer and transmits cellular signals by changing between different conformational states. It is important to elucidate the lipid–protein interactions in order to understand the functional mechanism of GPCRs. Here, all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method were performed on both the inactive (antagonist bound) and active (agonist and G-protein bound) A₁AR, which was embedded in a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) lipid bilayer. In the GaMD simulations, the membrane lipids played a key role in stabilizing different conformational states of the A₁AR. Our simulations further identified important regions of the receptor that interacted distinctly with the lipids in highly correlated manner. Activation of the A₁AR led to differential dynamics in the upper and lower leaflets of the lipid bilayer. In summary, GaMD enhanced simulations have revealed strongly coupled dynamics of the GPCR and lipids that depend on the receptor activation state. © 2019 Wiley Periodicals, Inc.     

Reduction of conformational mobility and aggregation in W60G β2‐microglobulin: assessment by 15N NMR relaxation

The amyloid pathology associated with long‐term haemodialysis is due to the deposition of β₂‐microglobulin, the non‐polymorphic light chain of class I major histocompatibility complex, that accumulates at bone joints into amyloid fibrils. Several lines of evidence show the relevance of the tryptophan residue at position 60 for the fibrillogenic transition of the protein. A comparative ¹⁵N NMR relaxation analysis is presented for wild‐type human β₂‐microglobulin and W60G β₂‐microglobulin, i.e. the mutant with a glycyne replacing the natural tryptophan residue at position 60. The experimental data, collected at 11.4 T and 310 K, were analyzed by means of the reduced spectral density approach. Molecular dynamics (MD) simulations and corresponding thermodynamic integration, together with hydrodynamic calculations were performed to support data interpretation. The analysis results for the mutant protein are consistent with a reduced aggregation with respect to the wild‐type counterpart, as a consequence of an increased conformational rigidity probed by either NMR relaxation and MD simulations. Although dynamics in solution is other than fibrillar competence, the assessed properties of the mutant protein can be related with its reduced ability of forming fibrils when seeded in 20% trifluoroethanol. Copyright © 2013 John Wiley & Sons, Ltd.     

Conformational Preferences of a β‐Octapeptide as Function of Solvent and Force‐Field Parameters

The ability to design properly folded β‐peptides with specific biological activities requires detailed insight into the relationship between the amino acid sequence and the secondary and/or tertiary structure of the peptide. One of the most frequently used spectroscopic techniques for resolving the structure of a biomolecule is NMR spectroscopy. Because only signal intensities and frequencies are recorded in the experiment, a conformational interpretation of the recorded data is not straightforward, especially for flexible molecules. The occurrence of conformational and/or time averaging, and the limited amount and accuracy of experimental data hamper the precise conformational determination of a biomolecule. In addition, the relation between experimental observables with the underlying conformational ensemble is often only approximately known, thereby aggravating the difficulty of structure determination of biomolecules. The problematic aspects of structure refinement based on NMR nuclear Overhauser effect (NOE) intensities and ³J‐coupling data are illustrated by simulating a β‐octapeptide in explicit MeOH and H₂O as solvents using three different force fields. NMR Data indicated that this peptide would fold into a 3₁₄‐helix in MeOH and into a hairpin in H₂O. Our analysis focused on the conformational space visited by the peptide, on structural properties of the peptide, and on agreement of the MD trajectories with available NMR data. We conclude that 1) although the 3₁₄‐helical structure is present when the peptide is solvated in MeOH, it is not the only relevant conformation, and that 2) the NMR data set available for the peptide, when solvated in H₂O, does not provide sufficient information to derive a single secondary structure, but rather a multitude of folds that fulfill the NOE data set.     

On the Use of a Supramolecular Coarse‐Grained Model for the Solvent in Simulations of the Folding Equilibrium of an Octa‐β‐peptide in MeOH and H2O

Molecular dynamics (MD) simulation can give a detailed picture of conformational equilibria of biomolecules, but it is only reliable if the force field used in the simulation is accurate, and the sampling of the conformational space accessible to the biomolecule shows many (un)folding transitions to allow for precise averages of observable quantities. Here, the use of coarse‐grained (CG) solvent MeOH and H₂O models to speed up the sampling of the conformational equilibria of an octa‐β‐peptide is investigated. This peptide is thought to predominantly adopt a 3₁₄‐helical fold when solvated in MeOH, and a hairpin fold when solvated in H₂O on the basis of the NMR data. Various factors such as the chirality of a residue, a force‐field modification for the solute, coarse‐graining of the solvent model, and an extension of the nonbonded interaction cut‐off radius are shown to influence the simulated conformational equilibria and the agreement with the experimental NMR data for the octa‐β‐peptide.     

An Example of Head-to-Head Dimerization of β4.4-Helices

¹H-NMR and vapor-pressure osmometry results are presented, which indicate the occurrence of a rapid equilibrium involving the head-to-head dimerization of β⁴⁴-helices in chloroform solutions of HCO-_L-Ile-(_D-AIle-_L-Ile)₄-OMe. This equilibrium typifies the one that, in Urry's view, would be responsible for the formation and breaking down of the ion-conducting channels formed by gramicidin A in lipid bilayers.     

A novel azocompound, 2‐(4‐phenylazoaniline)‐4‐phenylphenol: Spectroscopic and quantum‐chemical approach

A novel azocompound with two nonequivalents azo groups, 2‐(4‐phenylazoaniline)‐4‐phenylphenol, was synthesized and characterized by spectroscopic and computational analysis. An intramolecular hydrogen bonding (HB), ? O₁? H₁ ··· N₁? , involving the ? N₁?N₂? group and the proton in a neighbor hydroxyl moiety, was identified. It was found responsible for a characteristic π‐conjugated H₁? O₁? C₁₈?C₁₃? N₂?N₁? six‐membered cyclic fragment. It is worth noting that this azo group is involved in an azo‐hydrazo equilibrium, being the azo form the most stable one. This resonance‐assisted HB was characterized using the OH‐related infrared bands and the corresponding signals in ¹H NMR. In addition, conformational studies and geometrical and electronic parameter calculations were performed using the density functional theory, at B3LYP/6‐311++G** level. Bond and ring critical points were identified using the atoms in molecules theory, which allowed confirming the intramolecular HB. The second azo‐group cannot be involved in HB, but it also presents two stereoisomerics forms corresponding to cis (Z) and trans (E) configurations, with the later being the one with the lowest energy. © 2013 Wiley Periodicals, Inc.     

Unwinding of the C‐Terminal Residues of Neuropeptide Y is critical for Y2 Receptor Binding and Activation

Despite recent breakthroughs in the structural characterization of G‐protein‐coupled receptors (GPCRs), there is only sparse data on how GPCRs recognize larger peptide ligands. NMR spectroscopy, molecular modeling, and double‐cycle mutagenesis studies were integrated to obtain a structural model of the peptide hormone neuropeptide Y (NPY) bound to its human G‐protein‐coupled Y₂ receptor (Y₂R). Solid‐state NMR measurements of specific isotope‐labeled NPY in complex with in vitro folded Y₂R reconstituted into phospholipid bicelles provided the bioactive structure of the peptide. Guided by solution NMR experiments, it could be shown that the ligand is tethered to the second extracellular loop by hydrophobic contacts. The C‐terminal α‐helix of NPY, which is formed in a membrane environment in the absence of the receptor, is unwound starting at T³² to provide optimal contacts in a deep binding pocket within the transmembrane bundle of the Y₂R.     

NMR spectroscopic conformational analysis of 4‐methylene‐cyclohexyl pivalate—The effect of sp2 hybridization

The conformational equilibrium of the axial/equatorial conformers of 4‐methylene‐cyclohexyl pivalate is studied by dynamic NMR spectroscopy in a methylene chloride/freon mixture. At 153 K, the ring interconversion gets slow on the nuclear magnetic resonance timescale, the conformational equilibrium (?ΔG°) can be examined, and the barrier to ring interconversion (ΔG^#) can be determined. The structural influence of sp² hybridization on both ΔG° and ΔG^# of the cyclohexyl moiety can be quantified.     

New Open‐Chain and Cyclic Tetrapeptides,Consisting of α‐, β2‐, and β3‐Amino‐Acid Residues,as Somatostatin Mimics – A Survey

Cyclo‐β‐tetrapeptides are known to adopt a conformation with an intramolecular transannular hydrogen bond in solution. Analysis of this structure reveals that incorporation of a β²‐amino‐acid residue should lead to mimics of ‘α‐peptidic β‐turns’ (cf. A, B, C ). It is also known that short‐chain mixed β/α‐peptides with appropriate side chains can be used to mimic interactions between α‐peptidic hairpin turns and G protein‐coupled receptors. Based on these facts, we have now prepared a number of cyclic and open‐chain tetrapeptides, 7 – 20 , consisting of α‐, β²‐, and β³‐amino‐acid residues, which bear the side chains of Trp and Lys, and possess backbone configurations such that they should be capable of mimicking somatostatin in its affinity for the human SRIF receptors (hsst_1–5). All peptides were prepared by solid‐phase coupling by the Fmoc strategy. For the cyclic peptides, the three‐dimensional orthogonal methodology (Scheme 3) was employed with best success. The new compounds were characterized by high‐resolution mass spectrometry, NMR and CD spectroscopy, and, in five cases, by a full NMR‐solution‐structure determination (in MeOH or H₂O; Fig. 4). The affinities of the new compounds for the receptors hsst_1–5 were determined by competition with [¹²⁵I]LTT‐SRIF₂₈ or [¹²⁵I] [Tyr¹⁰]‐CST₁₄. In Table 1, the data are listed, together with corresponding values of all β‐ and γ‐peptidic somatostatin/Sandostatin^® mimics measured previously by our groups. Submicromolar affinities have been achieved for most of the human SRIF receptors hsst_1–5. Especially high, specific binding affinities for receptor hsst₄ (which is highly expressed in lung and brain tissue, although still of unknown function!) was observed with some of the β‐peptidic mimics. In view of the fact that numerous peptide‐activated G protein‐coupled receptors (GPCRs) recognize ligands with turn structure (Table 2), the results reported herein are relevant far beyond the realm of somatostatin: many other peptide GPCRs should be ‘reached’ with β‐ and γ‐peptidic mimics as well, and these compounds are proteolytically and metabolically stable, and do not need to be cell‐penetrating for this purpose (Fig. 5).     

Chalcogenide Derivatives of the seco‐Cubane [Sn3(μ2‐NHtBu)2(μ2‐NtBu)(μ3‐NtBu)]

A series of monochalcogenide derivatives of the seco‐cubane [Sn₃(μ₂‐NHtBu)₂(μ₂‐NtBu)(μ₃‐NtBu)] has been prepared and characterized by NMR and X‐ray crystallographic studies. These complexes exhibit different tin‐chalcogen bonding modes. In the case of the monotelluride, a terminal Sn=Te bond was observed in solution and in the solid state, whereas for the monosulfide, a μ₂ bridging mode was adopted by the sulfur atoms. The monoselenide was found to employ both bonding modes in solution, although only the terminal Sn=Se bonding mode was structurally characterized. The complexes undergo chalcogen exchange between tin atoms in solution, and this process was studied by variable temperature NMR.     

Effect of pH on the rotational conformations of 1,3‐diamino‐2‐hydroxypropane

The effect of pH on the rotational conformations of 1,3‐diamino‐2‐hydroxypropane in aqueous solution was investigated by proton NMR. Both the observed chemical shifts and coupling constants were used to calculate experimental pK_a values. The observed couplings were correlated with the expected couplings for the various possible staggered conformations to try to determine the pattern of conformations for the diamine and its conjugate acids. The best fits suggested a modest preference for the gauche–gauche conformation, especially at low pH, where the diprotonated hydroxydiamine predominates. In methanol, dimethyl sulfoxide and trichloromethane solutions, it was only possible to evaluate the conformational equilibria of the diamine. Slow proton exchange, which caused uncertainties in both chemical shifts and couplings for the monoprotonated and unprotonated diamine, nullified efforts to determine whether or not hydrogen bonding was important for these species in less polar solvents. Copyright © 2002 John Wiley & Sons, Ltd.     

Diastereoselective Alkyl Grignard 1,4‐Additions to para‐Substituted (2R)‐N‐Cinnamoylbornane‐10,2‐sultam Derivatives: Influence of N‐Atom Pyramidalization

Several typical ¹³C‐NMR displacements (of C?O, C(α), C(β), and C_ipso), as well as conformational or energy properties (S? N? C?O dihedral angle, ΔE syn/anti; HOMO/LUMO) could be correlated with the electronic parameters of p‐substituted N‐cinnamoylbornane‐10,2‐sultams 2 . Even under nonchelating conditions, the pyramidalization of the sultam N‐atom decreases for electron‐attracting p‐substituents, inducing a modification of the sultam‐ring puckering. Detailed comparison of the X‐ray structure analyses of 2b, 2d , and 2m showed that the orientation of the sterically directing pseudo‐axial S?O(2) and H? C(2) is modified and precludes any conclusion about the π‐facial stereoelectronic influence of the N lone pair on the alkyl Grignard 1,4‐addition. We also showed that the aggregating alkyl Grignard reagent may be used in equimolar fashion, demonstrating that the sultam moiety is chelated with a Lewis acid such as MgBr₂. The Schlenk equilibrium may also be used to generate the appropriate conditions of effective 1,4‐diastereoselectivity. Although the anti‐s‐cis/syn‐s‐cis difference of conformational energies for N‐cinnamoyl derivatives 2 is higher than for the simple N‐crotonoyl analogue, an X‐ray structure analysis of the SO₂/C?O syn derivative 10 confirms the predictive validity of our conformational calculations for ΔE≤1.8 kcal/mol.     

A zwitterion produced by a strong intramolecular N→B interaction in 1‐hydroxy‐2‐(pyridin‐2‐ylcarbonyl)benzo[d][1,2,3]diazaborinine

2‐Acylated 2,3,1‐benzodiazaborines can display unusual structures and reactivities. The crystal structure analysis of the boron heterocycle obtained by condensing 2‐formylphenylboronic acid and picolinohydrazide reveals it to be an N→B‐chelated zwitterionic tetracycle (systematic name: 1‐hydroxy‐11‐oxo‐9,10,17λ⁵‐triaza‐1λ⁴‐boratetracyclo[8.7.0.0^2,7.0^12,17]heptadeca‐3,5,7,12,14,16‐hexaen‐17‐ylium‐1‐uide), C₁₃H₁₀BN₃O₂, produced by the intramolecular addition of the Lewis basic picolinoyl N atom of 1‐hydroxy‐2‐(pyridin‐2‐ylcarbonyl)benzo[d][1,2,3]diazaborinine to the boron heterocycle B atom acting as a Lewis acid. Neither of the other two pyridinylcarbonyl isomers (viz. nicotinoyl and isonicotinoyl) are able to adopt such a structure for geometric reasons. A favored yet reversible chelation equilibrium provides an explanation for the slow D₂O exchange observed for the OH resonance in the ¹H NMR spectrum, as well as for its unusual upfield chemical shift. Deuterium exchange may take place solely in the minor open (unchelated) species present in solution.