Alternative Conformations Research Articles

Quantitatively Accounting for Protein Reorganization in Computer-Aided Drug Design.

Structure-based drug design frequently operates under the assumption that a single holo structure is relevant. However, a large number of crystallographic examples clearly show that multiple conformations are possible. In those cases, the protein reorganization free energy must be known to accurately predict binding free energies for ligands. Only then can the energetic preference among these multiple protein conformations be utilized to design ligands with stronger binding potency and selectivity. Here, we present a computational method to quantify these protein reorganization free energies. We test it on two retrospective drug design cases, Abl kinase and HSP90, and illustrate how alternative holo conformations can be derisked and lead to large boosts in affinity. This method will allow computer-aided drug design to better support complex protein targets.

Unusual substructure conformations observed in crystal structures of a dicistrovirus RNA-dependent RNA polymerase suggest contribution of the N-terminal extension in proper folding

The Dicistroviridae is a virus family that includes many insect pathogens. These viruses contain a positive-sense RNA genome that is replicated by the virally encoded RNA-dependent RNA polymerase (RdRP) also named 3Dpol. Compared with the Picornaviridae RdRPs such as poliovirus (PV) 3Dpol, the Dicistroviridae representative Israeli acute paralysis virus (IAPV) 3Dpol has an additional N-terminal extension (NE) region that is about 40-residue in length. To date, both the structure and catalytic mechanism of the Dicistroviridae RdRP have remain elusive. Here we reported crystal structures of two truncated forms of IAPV 3Dpol, namely Δ85 and Δ40, both missing the NE region, and the 3Dpol protein in these structures exhibited three conformational states. The palm and thumb domains of these IAPV 3Dpol structures are largely consistent with those of the PV 3Dpol structures. However, in all structures, the RdRP fingers domain is partially disordered, while different conformations of RdRP substructures and interactions between them are also present. In particular, a large-scale conformational change occurred in the motif B-middle finger region in one protein chain of the Δ40 structure, while a previously documented alternative conformation of motif A was observed in all IAPV structures. These experimental data on one hand show intrinsic conformational variances of RdRP substructures, and on the other hand suggest possible contribution of the NE region in proper RdRP folding in IAPV.

PPARγ Corepression Involves Alternate Ligand Conformation and Inflation of H12 Ensembles.

Inverse agonists of peroxisome proliferator activated receptor γ (PPARγ) have emerged as safer alternatives to full agonists for their reduced side effects while still maintaining impressive insulin-sensitizing properties. To shed light on their molecular mechanism, we characterized the interaction of the PPARγ ligand binding domain with SR10221. X-ray crystallography revealed a novel binding mode of SR10221 in the presence of a transcriptionally repressing corepressor peptide, resulting in much greater destabilization of the activation helix, H12, than without corepressor peptide. Electron paramagnetic resonance provided in-solution complementary protein dynamic data, which revealed that for SR10221-bound PPARγ, H12 adopts a plethora of conformations in the presence of corepressor peptide. Together, this provides the first direct evidence for corepressor-driven ligand conformation for PPARγ and will allow the development of safer and more effective insulin sensitizers suitable for clinical use.

Translation regulation by a guanidine-II riboswitch is highly tunable in sensitivity, dynamic range, and apparent cooperativity.

Riboswitches function as important translational regulators in bacteria. Comprehensive mutational analysis of transcriptional riboswitches has been used to probe the energetic intricacies of interplay between the aptamer and expression platform, but translational riboswitches have been inaccessible to massively parallel techniques. The guanidine-II (gdm-II) riboswitch is an exclusively translational class. We have integrated RelE cleavage with next-generation sequencing to quantify ligand-dependent changes in translation initiation for all single and double mutations of the Pseudomonas aeruginosa gdm-II riboswitch, a total of more than 23,000 variants. This extensive mutational analysis is consistent with the prominent features of the bioinformatic consensus. These data indicate, unexpectedly, that direct sequestration of the Shine-Dalgarno sequence is dispensable for riboswitch function. Additionally, this comprehensive data set reveals important positions not identified in previous computational and crystallographic studies. Mutations in the variable linker region stabilize alternate conformations. The double mutant data reveal the functional importance of the previously modeled P0b helix formed by the 5' and 3' tails that serves as the basis for translational control. Additional mutations to GU wobble base pairs in both P1 and P2 reveal how the apparent cooperativity of the system involves an intricate network of communication between the two binding sites. This comprehensive examination of a translational riboswitch's expression platform illuminates how the riboswitch is precisely tuned and tunable with regard to ligand sensitivity, the amplitude of expression between ON and OFF states, and the cooperativity of ligand binding.

FLEXR: automated multi-conformer model building using electron-density map sampling.

Protein conformational dynamics that may inform biology often lie dormant in high-resolution electron-density maps. While an estimated ∼18% of side chains in high-resolution models contain alternative conformations, these are underrepresented in current PDB models due to difficulties in manually detecting, building and inspecting alternative conformers. To overcome this challenge, we developed an automated multi-conformer modeling program, FLEXR. Using Ringer-based electron-density sampling, FLEXR builds explicit multi-conformer models for refinement. Thereby, it bridges the gap of detecting hidden alternate states in electron-density maps and including them in structural models for refinement, inspection and deposition. Using a series of high-quality crystal structures (0.8-1.85 Å resolution), we show that the multi-conformer models produced by FLEXR uncover new insights that are missing in models built either manually or using current tools. Specifically, FLEXR models revealed hidden side chains and backbone conformations in ligand-binding sites that may redefine protein-ligand binding mechanisms. Ultimately, the tool facilitates crystallographers with opportunities to include explicit multi-conformer states in their high-resolution crystallographic models. One key advantage is that such models may better reflect interesting higher energy features in electron-density maps that are rarely consulted by the community at large, which can then be productively used for ligand discovery downstream. FLEXR is open source and publicly available on GitHub at https://github.com/TheFischerLab/FLEXR.

Structures of factor XI and prekallikrein bound to domain 6 of high–molecular weight kininogen reveal alternate domain 6 conformations and exosites

BackgroundHigh–molecular weight kininogen (HK) circulates in plasma as a complex with zymogen prekallikrein (PK). HK is both a substrate and a cofactor for activated plasma kallikrein, and the principal exosite interactions occur between PK N-terminal apple domains and the C-terminal D6 domain of HK. ObjectivesTo determine the structure of the complex formed between PK apple domains and an HKD6 fragment and compare this with the coagulation factor XI (FXI)-HK complex. MethodsWe produced recombinant FXI and PK heavy chains (HCs) spanning all 4 apple domains. We cocrystallized PKHC (and subsequently FXIHC) with a 31-amino acid synthetic peptide spanning HK residues Ser565-Lys595 and determined the crystal structure. We also analyzed the full-length FXI-HK complex in solution using hydrogen deuterium exchange mass spectrometry. ResultsThe 2.3Å PKHC-HK peptide crystal structure revealed that the HKD6 sequence WIPDIQ (Trp569-Gln574) binds to the apple 1 domain and HK FNPISDFPDT (Phe582-Thr591) binds to the apple 2 domain with a flexible intervening sequence resulting in a bent double conformation. A second 3.2Å FXIHC-HK peptide crystal structure revealed a similar interaction with the apple 2 domain but an alternate, straightened conformation of the HK peptide where residues LSFN (Leu579-Asn583) interacts with a unique pocket formed between the apple 2 and 3 domains. HDX-MS of full length FXI-HK complex in solution confirmed interactions with both apple 2 and apple 3. ConclusionsThe alternate conformations and exosite binding of the HKD6 peptide likely reflects the diverging relationship of HK to the functions of PK and FXI.

Water stabilizes an alternate turn conformation in horse heart myoglobin

Comparison of myoglobin structures reveals that protein isolated from horse heart consistently adopts an alternate turn conformation in comparison to its homologues. Analysis of hundreds of high-resolution structures discounts crystallization conditions or the surrounding amino acid protein environment as explaining this difference, that is also not captured by the AlphaFold prediction. Rather, a water molecule is identified as stabilizing the conformation in the horse heart structure, which immediately reverts to the whale conformation in molecular dynamics simulations excluding that structural water.

Comparative Study of High-Resolution LysB29(Nε-myristoyl) des(B30) Insulin Structures Display Novel Dynamic Causal Interrelations in Monomeric-Dimeric Motions

The treatment of insulin-dependent diabetes mellitus is characterized by artificial supplementation of pancreatic β-cell ability to regulate sugar levels in the blood. Even though various insulin analogs are crucial for reasonable glycemic control, understanding the dynamic mechanism of the insulin analogs may help to improve the best-protracted insulin analog to assist people with type 1 diabetes (T1D) to live comfortably while maintaining tight glycemic control. Here, we present the high-resolution crystal structure of NN304, known as insulin detemir, to 1.7 Å resolution at cryogenic temperature. We computationally further investigated our crystal structure’s monomeric-dimeric conformation and dynamic profile by comparing it with a previously available detemir structure (PDB ID: 1XDA). Our structure (PDB ID: 8HGZ), obtained at elevated pH, provides electrostatically triggered minor movements in the equilibrium between alternate conformational substates compared to the previous structure, suggesting it might induce an intermediate state in the dissociation pathway of the insulin detemir’s hexamer:dihexamer equilibrium. Supplemented with orientational cross-correlation analysis by a Gaussian network model (GNM), this alternate oligomeric conformation offers the distinct cooperative motions originated by loose coupling of distant conformational substates of a protracted insulin analog that has not been previously observed.

Modeling with Alternate Locations in X-ray Protein Structures.

In many molecular modeling applications, the standard procedure is still to handle proteins as single, rigid structures. While the importance of conformational flexibility is widely known, handling it remains challenging. Even the crystal structure of a protein usually contains variability exemplified in alternate side chain orientations or backbone segments. This conformational variability is encoded in PDB structure files by so-called alternate locations (AltLocs). Most modeling approaches either ignore AltLocs or resolve them with simple heuristics early on during structure import. We analyzed the occurrence and usage of AltLocs in the PDB and developed an algorithm to automatically handle AltLocs in PDB files enabling all structure-based methods using rigid structures to take the alternative protein conformations described by AltLocs into consideration. A respective software tool named AltLocEnumerator can be used as a structure preprocessor to easily exploit AltLocs. While the amount of data makes it difficult to show impact on a statistical level, handling AltLocs has a substantial impact on a case-by-case basis. We believe that the inspection and consideration of AltLocs is a very valuable approach in many modeling scenarios.

Synthesis, characterization, stereochemistry, antibacterial and antioxidant activity of N-chloroacetyl-2,6-diphenylpiperidines

A new series of N-chloroacetyl-2,6-diphenylpiperidines 13–18 have been synthesized from their corresponding piperidines 7–12. The synthesized compound 14 is unambiguously characterized using IR, mass, 1D and 2D (1H–1H, COSY & 1H–13C, HSQC) NMR spectral techniques and the remaining compounds 13–16 & 18 are characterized by IR, mass and 1D NMR spectra. The NMR spectral data indicate that the compound 17 prefers to adopt alternate chair conformation and remaining compounds prefer to adopt boat conformation. The X-ray crystal structure solved for N-chloroacetyl-3-isopropyl-2,6-diphenylpiperidine 16 supports the distorted boat conformation in the solid state. Anti-bacterial evaluation shows that the compounds 15 & 16 demonstrate superior activity against all the tested bacteria. All the synthesized compounds show promising antioxidant activities compared to the standard BHA (butylated hydroxyanisole).

The structure of AgamOBP5 in complex with the natural insect repellents Carvacrol and Thymol: Crystallographic, fluorescence and thermodynamic binding studies

Among several proteins participating in the olfactory perception process of insects, Odorant Binding Proteins (OBPs) are today considered valid targets for the discovery of compounds that interfere with their host-detection behavior. The 3D structures of Anopheles gambiae mosquito AgamOBP1 in complex with the known synthetic repellents DEET and Icaridin have provided valuable information on the structural characteristics that govern their selective binding.However, no structure of a plant-derived repellent bound to an OBP has been available until now. Herein, we present the novel three-dimensional crystal structures of AgamOBP5 in complex with two natural phenolic monoterpenoid repellents, Carvacrol and Thymol, and the MPD molecule.Structural analysis revealed that both monoterpenoids occupy a binding site (Site-1) by adopting two alternative conformations. An additional Carvacrol was also bound to a secondary site (Site-2) near the central cavity entrance. A protein-ligand hydrogen-bond network supplemented by van der Waals interactions spans the entire binding cavity, bridging α4, α6, and α3 helices and stabilizing the overall structure.Fluorescence competition and Differential Scanning Calorimetry experiments verified the presence of two binding sites and the stabilization effect on AgamOBP5. While Carvacrol and Thymol bind to Site-1 with equal affinity in the submicromolar range, they exhibit a significantly lower and distinct binding capacity for Site-2 with Kd's of ~7 μΜ and ~18 μΜ, respectively.Finally, a comparison of AgamOBP5 complexes with the AgamOBP4-Indole structure revealed that variations of ligand-interacting aminoacids such as A109T, I72M, A112L, and A105T cause two structurally similar and homologous proteins to display different binding specificities.

Structure-Based Regulatory Role for the 5′UTR of RCNMV RNA2

Red clover necrotic mosaic virus (RCNMV) is a segmented positive-strand RNA virus consisting of RNA1 and RNA2. Previous studies demonstrated that efficient translation of RCNMV RNA2 requires de novo synthesis of RNA2 during infections, suggesting that RNA2 replication is required for its translation. We explored a potential mechanism underlying the regulation of replication-associated translation of RNA2 by examining RNA elements in its 5' untranslated region (5'UTR). Structural analysis of the 5'UTR suggested that it can form two mutually exclusive configurations: a more thermodynamically stable conformation, termed the 5'-basal stem structure (5'BS), in which 5'-terminal sequences are base paired, and an alternative conformation, where the 5'-end segment is single stranded. Functional mutational analysis of the 5'UTR structure indicated that (i) 43S ribosomal subunits enter at the very 5'-end of RNA2; (ii) the alternative conformation, containing unpaired 5'-terminal nucleotides, mediates efficient translation; (iii) the 5'BS conformation, with a paired 5'-end segment, supresses translation; and (iv) the 5'BS conformation confers stability to RNA2 from 5'-to-3' exoribonuclease Xrn1. Based on our results, we suggest that during infections, newly synthesized RNA2s transiently adopt the alternative conformation to allow for efficient translation, then refold into the 5'BS conformation, which supresses translation and promotes efficient RNA2 replication. The potential advantages of this proposed 5'UTR-based regulatory mechanism for coordinating RNA2 translation and replication are discussed.

Room-temperature crystallography reveals altered binding of small-molecule fragments to PTP1B.

Much of our current understanding of how small-molecule ligands interact with proteins stems from X-ray crystal structures determined at cryogenic (cryo) temperature. For proteins alone, room-temperature (RT) crystallography can reveal previously hidden, biologically relevant alternate conformations. However, less is understood about how RT crystallography may impact the conformational landscapes of protein-ligand complexes. Previously, we showed that small-molecule fragments cluster in putative allosteric sites using a cryo crystallographic screen of the therapeutic target PTP1B (Keedy et al., 2018). Here, we have performed two RT crystallographic screens of PTP1B using many of the same fragments, representing the largest RT crystallographic screens of a diverse library of ligands to date, and enabling a direct interrogation of the effect of data collection temperature on protein-ligand interactions. We show that at RT, fewer ligands bind, and often more weakly - but with a variety of temperature-dependent differences, including unique binding poses, changes in solvation, new binding sites, and distinct protein allosteric conformational responses. Overall, this work suggests that the vast body of existing cryo-temperature protein-ligand structures may provide an incomplete picture, and highlights the potential of RT crystallography to help complete this picture by revealing distinct conformational modes of protein-ligand systems. Our results may inspire future use of RT crystallography to interrogate the roles of protein-ligand conformational ensembles in biological function.

Automatic and accurate ligand structure determination guided by cryo-electron microscopy maps

Advances in cryo-electron microscopy (cryoEM) and deep-learning guided protein structure prediction have expedited structural studies of protein complexes. However, methods for accurately determining ligand conformations are lacking. In this manuscript, we develop EMERALD, a tool for automatically determining ligand structures guided by medium-resolution cryoEM density. We show this method is robust at predicting ligands along with surrounding side chains in maps as low as 4.5 Å local resolution. Combining this with a measure of placement confidence and running on all protein/ligand structures in the EMDB, we show that 57% of ligands replicate the deposited model, 16% confidently find alternate conformations, 22% have ambiguous density where multiple conformations might be present, and 5% are incorrectly placed. For five cases where our approach finds an alternate conformation with high confidence, high-resolution crystal structures validate our placement. EMERALD and the resulting analysis should prove critical in using cryoEM to solve protein-ligand complexes.

Conformational multiplicity of bacterial ferric binding protein revealed by small angle x-ray scattering and molecular dynamics calculations.

This study combines molecular dynamics (MD) simulations with small angle x-ray scattering (SAXS) measurements to investigate the range of conformations that can be adopted by a pH/ionic strength (IS) sensitive protein and to quantify its distinct populations in solution. To explore how the conformational distribution of proteins may be modified in the environmental niches of biological media, we focus on the periplasmic ferric binding protein A (FbpA) from Haemophilus influenzae involved in the mechanism by which bacteria capture iron from higher organisms. We examine iron-binding/release mechanisms of FbpA in varying conditions simulating its biological environment. While we show that these changes fall within the detectable range for SAXS as evidenced by differences observed in the theoretical scattering patterns calculated from the crystal structure models of apo and holo forms, detection of conformational changes due to the point mutation D52A and changes in ionic strength (IS) from SAXS scattering profiles have been challenging. Here, to reach conclusions, statistical analyses with SAXS profiles and results from different techniques were combined in a complementary fashion. The SAXS data complemented by size exclusion chromatography point to multiple and/or alternative conformations at physiological IS, whereas they are well-explained by single crystallographic structures in low IS buffers. By fitting the SAXS data with unique conformations sampled by a series of MD simulations under conditions mimicking the buffers, we quantify the populations of the occupied substates. We also find that the D52A mutant that we predicted by coarse-grained computational modeling to allosterically control the iron binding site in FbpA, responds to the environmental changes in our experiments with conformational selection scenarios that differ from those of the wild type.

Elucidating the distinct contributions of miR-122 in the HCV life cycle reveals insights into virion assembly.

Efficient hepatitis C virus (HCV) RNA accumulation is dependent upon interactions with the human liver-specific microRNA, miR-122. MiR-122 has at least three roles in the HCV life cycle: it acts as an RNA chaperone, or 'riboswitch', allowing formation of the viral internal ribosomal entry site; it provides genome stability; and promotes viral translation. However, the relative contribution of each role in HCV RNA accumulation remains unclear. Herein, we used point mutations, mutant miRNAs, and HCV luciferase reporter RNAs to isolate each of the roles and evaluate their contribution to the overall impact of miR-122 in the HCV life cycle. Our results suggest that the riboswitch has a minimal contribution in isolation, while genome stability and translational promotion have similar contributions in the establishment phase of infection. However, in the maintenance phase, translational promotion becomes the dominant role. Additionally, we found that an alternative conformation of the 5' untranslated region, termed SLIIalt, is important for efficient virion assembly. Taken together, we have clarified the overall importance of each of the established roles of miR-122 in the HCV life cycle and provided insight into the regulation of the balance between viral RNAs in the translating/replicating pool and those engaged in virion assembly.

Characterisation of the complete mitochondrial genome of Taraxacum mongolicum revealed five repeat-mediated recombinations.

We reported the mitochondrial genome of Taraxacum mongolicum for the first time. Five pairs of repeats that can mediate recombination were validated, leading to multiple conformations of genome. Taraxacum mongolicum belongs to the Asteraceae family and has important pharmaceutical value. To explore the possible interaction between the organelle genomes, we assembled the complete mitochondrial genome (mitogenome) of T. mongolicum using Illumina and Oxford Nanopore sequencing data. This genome corresponded to a circular molecule 304,467bp long. It encodes 52 unique genes including 31 protein-coding, 3 ribosomal RNA (rRNA) and 18 transfer RNA (tRNA) genes. In addition to the single circular conformation, the existence of alternative conformations mediated by five repetitive sequences in the mitogenome was identified and validated. Recombination mediated by the inverted repeats resulted in two conformations. Conversely, recombination mediated by the two direct repeats broke one large circular molecule into two subgenomic circular molecules. Furthermore, we identified 12 homologous fragments by comparing the sequences of mitogenome and plastome, including eight complete tRNA genes. Lastly, we identified a total of 278 RNA-editing sites in protein-coding sequences based on RNA-seq data. Among them, cox1 and nad5 gene has the most sites (21), followed by the nad2 gene with 19 sites. We successfully validated 213 predicted RNA-editing sites using PCR amplification and Sanger sequencing. This project reported the first mitogenome of T. mongolicum and demonstrated its multiple conformations generated by repeat-mediated recombination. This genome could provide critical information for the molecular breeding of T. mongolicum, and also be used as a reference genome for other species of the genus Taraxacum.

The role of K288A in conformational selection in LeuT.

Upon publication of the LeuT crystal structure in 2005, researchers noticed a positively-charged lysine residue (K288) protruding into the hydrophobic membrane, an energetically unfavorable arrangement. Subsequently, it was found that a mutation of this residue (K288A) increased the kinetic turnover rate of LeuT. Furthermore, it has been suggested that this change in kinetics was due to the alanine mutant sampling the “rate-limiting” inward-facing conformation more frequently, possibly due to its unfavorable interactions with the membrane. This served as the basis for the generation of an “inward-facing” LeuT crystal structure in 2013. We sought to test these hypotheses in solution using EPR spectroscopy. Evidence provided in this series of experiments suggests that neither hypothesis was supported. The K288A mutant in total appears to sample an occluded-like conformation that is not found in the ensemble of wild type structures. Furthermore, as these experiments were conducted in a micelle, changes in the protein's conformation were occurring in the absence of a membrane. The result of this work cautions against the use of mutations of conserved residues to “select” for alternative conformations in structural studies.

High-throughput biochemistry in RNA sequence space: predicting structure and function.

RNAs are central to fundamental biological processes in all known organisms. The set of possible intramolecular interactions of RNA nucleotides defines the range of alternative structural conformations of a specific RNA that can coexist, and these structures enable functional catalytic properties of RNAs and/or their productive intermolecular interactions with other RNAs or proteins. However, the immense combinatorial space of potential RNA sequences has precluded predictive mapping between RNA sequence and molecular structure and function. Recent advances in high-throughput approaches in vitro have enabled quantitative thermodynamic and kinetic measurements of RNA-RNA and RNA-protein interactions, across hundreds of thousands of sequence variations. In this Review, we explore these techniques, how they can be used to understand RNA function and how they might form the foundations of an accurate model to predict the structure and function of an RNA directly from its nucleotide sequence. The experimental techniques and modelling frameworks discussed here are also highly relevant for the sampling of sequence-structure-function space of DNAs and proteins.

Histone H3 core domain in chromatin with different DNA linker lengths studied by 1H-Detected solid-state NMR spectroscopy.

Chromatin, a dynamic protein-DNA complex that regulates eukaryotic genome accessibility and essential functions, is composed of nucleosomes connected by linker DNA with each nucleosome consisting of DNA wrapped around an octamer of histones H2A, H2B, H3 and H4. Magic angle spinning solid-state nuclear magnetic resonance (NMR) spectroscopy can yield unique insights into histone structure and dynamics in condensed nucleosomes and nucleosome arrays representative of chromatin at physiological concentrations. Recently we used J-coupling-based solid-state NMR methods to investigate with residue-specific resolution the conformational dynamics of histone H3 N-terminal tails in 16-mer nucleosome arrays containing 15, 30 or 60 bp DNA linkers. Here, we probe the H3 core domain in the 16-mer arrays as a function of DNA linker length via dipolar coupling-based 1H-detected solid-state NMR techniques. Specifically, we established nearly complete assignments of backbone chemical shifts for H3 core residues in arrays with 15-60 bp DNA linkers reconstituted with 2H,13C,15N-labeled H3. Overall, these chemical shifts were similar irrespective of the DNA linker length indicating no major changes in H3 core conformation. Notably, however, multiple residues at the H3-nucleosomal DNA interface in arrays with 15 bp DNA linkers exhibited relatively pronounced differences in chemical shifts and line broadening compared to arrays with 30 and 60 bp linkers. These findings are consistent with increased heterogeneity in nucleosome packing and structural strain within arrays containing short DNA linkers that likely leads to side-chains of these interfacial residues experiencing alternate conformations or shifts in their rotamer populations relative to arrays with the longer DNA linkers.