Heat-induced structural and chemical changes to a computationally designed miniprotein.

The de novo design of miniprotein inhibitors has recently emerged as a new technology to create proteins that bind with high affinity to specific therapeutic targets. Their size, ease of expression, and apparent high stability makes them excellent candidates for a new class of protein drugs. However, beyond circular dichroism melts and hydrogen/deuterium exchange experiments, little is known about their dynamics, especially at the elevated temperatures they seemingly tolerate quite well. To address that and gain insight for future designs, we have focused on identifying unintended and previously overlooked heat-induced structural and chemical changes in a particularly stable model miniprotein, EHEE_rd2_0005. Nuclear magnetic resonance (NMR) studies suggest the presence of dynamics on multiple time and temperature scales. Transiently elevating the temperature results in spontaneous chemical deamidation visible in the NMR spectra, which we validate using both capillary electrophoresis and mass spectrometry (MS) experiments. High temperatures also result in greatly accelerated intrinsic rates of hydrogen exchange and signal loss in NMR heteronuclear single quantum coherence spectra from local unfolding. These losses are in excellent agreement with both room temperature hydrogen exchange experiments and hydrogen bond disruption in replica exchange molecular dynamics simulations. Our analysis reveals important principles for future miniprotein designs and the potential for high stability to result in long-lived alternate conformational states.

Dynamics in natural and designed elastins and their relation to elastic fiber structure and recoil.

Elastin fibers assemble in the extracellular matrix from the precursor protein tropoelastin and provide the flexibility and spontaneous recoil required for arterial function. Unlike many proteins, a structure-function mechanism for elastin has been elusive. We have performed detailed NMR relaxation studies of the dynamics of the minielastins 24x' and 20x' using solution NMR, and of purified bovine elastin fibers in the presence and absence of mechanical stress using solid state NMR. The low sequence complexity of the minielastins enables us to determine average dynamical timescales and degrees of local ordering in the cross-link and hydrophobic modules separately using NMR relaxation by taking advantage of their residue-specific resolution. We find an extremely high degree of disorder, with order parameters for the entirety of the hydrophobic domains near zero, resembling that of simple chemical polymers and less than the order parameters that have been observed in other intrinsically disordered proteins. We find that average backbone order parameters in natural, purified elastin fibers are comparable to those found in 24x' and 20x' in solution. The difference in dynamics, compared with the minielastins, is that backbone correlation times are significantly slowed in purified elastin. Moreover, when elastin is mechanically stretched, the high chain disorder in purified elastin is retained, showing that any change in local ordering is below that detectable in our experiment. Combined with our previous finding of a 10-fold increase in the ordering of water when fully hydrated elastin fibers are stretched by 50%, these results support the hypothesis that stretch induced solvent ordering, i.e., the hydrophobic effect, is a key player in the elastic recoil of elastin as opposed to configurational entropy loss.

The hTERT core promoter forms three parallel G-quadruplexes.

The structure of the 68 nt sequence with G-quadruplex forming potential within the hTERT promoter is disputed. One model features a structure with three stacked parallel G-quadruplex units, while another features an unusual duplex hairpin structure adjoined to two stacked parallel and antiparallel quadruplexes. We report here the results of an integrated structural biology study designed to distinguish between these possibilities. As part of our study, we designed a sequence with an optimized hairpin structure and show that its biophysical and biochemical properties are inconsistent with the structure formed by the hTERT wild-type sequence. By using circular dichroism, thermal denaturation, nuclear magnetic resonance spectroscopy, analytical ultracentrifugation, small-angle X-ray scattering, molecular dynamics simulations and a DNase I cleavage assay we found that the wild type hTERT core promoter folds into a stacked, three-parallel G-quadruplex structure. The hairpin structure is inconsistent with all of our experimental data obtained with the wild-type sequence. All-atom models for both structures were constructed using molecular dynamics simulations. These models accurately predicted the experimental hydrodynamic properties measured for each structure. We found with certainty that the wild-type hTERT promoter sequence does not form a hairpin structure in solution, but rather folds into a compact stacked three-G-quadruplex conformation.

Recombinant expression and purification of AF1q and its interaction with T-cell Factor 7.

The protein ALL1 fused from chromosome 1q (AF1q) is overexpressed in a variety of cancers and acts to activate several signaling pathways that lead to oncogenesis. For example, AF1q has been shown to interact with T-cell Factor 7 (TCF7; also known as TCF1) from the Wnt/ß-catenin pathway resulting in the transcriptional activation of the CD44 and the enhancement of breast cancer metastasis. Despite the importance of AF1q in facilitating oncogenesis and metastasis, the structural and biophysical properties of AF1q remain largely unexplored due to the absence of a viable method for producing recombinant protein. Here, we report the overexpression of AF1q in E. coli as a fusion to a N-terminal His6-tag, which forms inclusion bodies (IBs) during expression. The AF1q protein was purified from IBs under denaturing conditions by immobilized metal affinity chromatography followed by a successful one-step dialysis refolding. Refolded AF1q was further purified to homogeneity by gel filtration chromatography resulting in an overall yield of 35 mg/L culture. Our nuclear magnetic resonance (NMR) and analytical ultracentrifugation (AUC) measurements reveal AF1q interacts with TCF7, specifically with TCF7's high-mobility group (HMG) domain (residues 154-237), which is, to our knowledge, the first biophysical characterization of the AF1q and TCF7 interaction.

Solution structure and functional investigation of human guanylate kinase reveals allosteric networking and a crucial role for the enzyme in cancer.

Human guanylate kinase (hGMPK) is the only known enzyme responsible for cellular GDP production, making it essential for cellular viability and proliferation. Moreover, hGMPK has been assigned a critical role in metabolic activation of antiviral and antineoplastic nucleoside-analog prodrugs. Given that hGMPK is indispensable for producing the nucleotide building blocks of DNA, RNA, and cGMP and that cancer cells possess elevated GTP levels, it is surprising that a detailed structural and functional characterization of hGMPK is lacking. Here, we present the first high-resolution structure of hGMPK in the apo form, determined with NMR spectroscopy. The structure revealed that hGMPK consists of three distinct regions designated as the LID, GMP-binding (GMP-BD), and CORE domains and is in an open configuration that is nucleotide binding-competent. We also demonstrate that nonsynonymous single-nucleotide variants (nsSNVs) of the hGMPK CORE domain distant from the nucleotide-binding site of this domain modulate enzymatic activity without significantly affecting hGMPK's structure. Finally, we show that knocking down the hGMPK gene in lung adenocarcinoma cell lines decreases cellular viability, proliferation, and clonogenic potential while not altering the proliferation of immortalized, noncancerous human peripheral airway cells. Taken together, our results provide an important step toward establishing hGMPK as a potential biomolecular target, from both an orthosteric (ligand-binding sites) and allosteric (location of CORE domain-located nsSNVs) standpoint.

Thrombin Exosite Maturation and Ligand Binding at ABE II Help Stabilize PAR-Binding Competent Conformation at ABE I.

Thrombin, derived from zymogen prothrombin (ProT), is a serine protease involved in procoagulation, anticoagulation, and platelet activation. Thrombin's actions are regulated through anion-binding exosites I and II (ABE I and ABE II) that undergo maturation during activation. Mature ABEs can utilize exosite-based communication to fulfill thrombin functions. However, the conformational basis behind such long-range communication and the resultant ligand binding affinities are not well understood. Protease activated receptors (PARs), involved in platelet activation and aggregation, are known to target thrombin ABE I. Unexpectedly, PAR3 (44-56) can already bind to pro-ABE I of ProT. Nuclear magnetic resonance (NMR) ligand-enzyme titrations were used to characterize how individual PAR1 (49-62) residues interact with pro-ABE I and mature ABE I. 1D proton line broadening studies demonstrated that binding affinities for native PAR1P (49-62, P54) and for the weak binding variant PAR1G (49-62, P54G) increased as ProT was converted to mature thrombin. 1H,15N-HSQC titrations revealed that PAR1G residues K51, E53, F55, D58, and E60 exhibited less affinity to pro-ABE I than comparable residues in PAR3G (44-56, P51G). Individual PAR1G residues then displayed tighter binding upon exosite maturation. Long-range communication between thrombin exosites was examined by saturating ABE II with phosphorylated GpIbα (269-282, 3Yp) and monitoring the binding of PAR1 and PAR3 peptides to ABE I. Individual PAR residues exhibited increased affinities in this dual-ligand environment supporting the presence of interexosite allostery. Exosite maturation and beneficial long-range allostery are proposed to help stabilize an ABE I conformation that can effectively bind PAR ligands.

Utilizing dipole-dipole cross-correlated relaxation for the measurement of angles between pairs of opposing CαHα-CαHα bonds in anti-parallel ß-sheets.

Dipole-dipole cross-correlated relaxation (CCR) between two spin pairs is rich with macromolecular structural and dynamic information on inter-nuclear bond vectors. Measurement of short range dipolar CCR rates has been demonstrated for a variety of inter-nuclear vector spin pairs in proteins and nucleic acids, where the multiple quantum coherence necessary for observing the CCR rate is created by through-bond scalar coupling. In principle, CCR rates can be measured for any pair of inter-nuclear vectors where coherence can be generated between one spin of each spin pair, regardless of both the distance between the two spin pairs and the distance of the two spins forming the multiple quantum coherence. In practice, however, long range CCR (lrCCR) rates are challenging to measure due to difficulties in linking spatially distant spin pairs. By utilizing through-space relaxation allowed coherence transfer (RACT), we have developed a new method for the measurement of lrCCR rates involving CαHα bonds on opposing anti-parallel ß-strands. The resulting lrCCR rates are straightforward to interpret since only the angle between the two vectors modulates the strength of the interference effect. We applied our lrCCR measurement to the third immunoglobulin-binding domain of the streptococcal protein G (GB3) and utilize published NMR ensembles and static NMR/X-ray structures to highlight the relationship between the lrCCR rates and the CαHα-CαHα inter-bond angle and bond mobility. Furthermore, we employ the lrCCR rates to guide the selection of sub-ensembles from the published NMR ensembles for enhancing the structural and dynamic interpretation of the data. We foresee this methodology for measuring lrCCR rates as improving the generation of structural ensembles by providing highly accurate details concerning the orientation of CαHα bonds on opposing anti-parallel ß-strands.

1H, 13C and 15N resonance assignment of human guanylate kinase.

Human guanylate kinase (hGMPK) is a critical enzyme that, in addition to phosphorylating its physiological substrate (d)GMP, catalyzes the second phosphorylation step in the conversion of anti-viral and anti-cancer nucleoside analogs to their corresponding active nucleoside analog triphosphates. Until now, a high-resolution structure of hGMPK is unavailable and thus, we studied free hGMPK by NMR and assigned the chemical shift resonances of backbone and side chain 1H, 13C, and 15N nuclei as a first step towards the enzyme's structural and mechanistic analysis with atomic resolution.

Deciphering Conformational Changes Associated with the Maturation of Thrombin Anion Binding Exosite I.

Thrombin participates in procoagulation, anticoagulation, and platelet activation. This enzyme contains anion binding exosites, ABE I and ABE II, which attract regulatory biomolecules. As prothrombin is activated to thrombin, pro-ABE I is converted into mature ABE I. Unexpectedly, certain ligands can bind to pro-ABE I specifically. Moreover, knowledge of changes in conformation and affinity that occur at the individual residue level as pro-ABE I is converted to ABE I is lacking. Such changes are transient and were not captured by crystallography. Therefore, we employed nuclear magnetic resonance (NMR) titrations to monitor development of ABE I using peptides based on protease-activated receptor 3 (PAR3). Proton line broadening NMR revealed that PAR3 (44-56) and more weakly binding PAR3G (44-56) could already interact with pro-ABE I on prothrombin. 1H-15N heteronuclear single-quantum coherence NMR titrations were then used to probe binding of individual 15N-labeled PAR3G residues (F47, E48, L52, and D54). PAR3G E48 and D54 could interact electrostatically with prothrombin and tightened upon thrombin maturation. The higher affinity for PAR3G D54 suggests the region surrounding thrombin R77a is better oriented to bind D54 than the interaction between PAR3G E48 and thrombin R75. Aromatic PAR3G F47 and aliphatic L52 both reported on significant changes in the chemical environment upon conversion of prothrombin to thrombin. The ABE I region surrounding the 30s loop was more affected than the hydrophobic pocket (F34, L65, and I82). Our NMR titrations demonstrate that PAR3 residues document structural rearrangements occurring during exosite maturation that are missed by reported X-ray crystal structures.

Kinetics of the Antibody Recognition Site in the Third IgG-Binding Domain of Protein G.

Protein dynamics occurring on a wide range of timescales play a crucial role in governing protein function. Particularly, motions between the globular rotational correlation time (τc ) and 40 µs (supra-τc window), strongly influence molecular recognition. This supra-τc window was previously hidden, owing to a lack of experimental methods. Recently, we have developed a high-power relaxation dispersion (RD) experiment for measuring kinetics as fast as 4 µs. For the first time, this method, performed under super-cooled conditions, enabled us to detect a global motion in the first ß-turn of the third IgG-binding domain of protein G (GB3), which was extrapolated to 371±115 ns at 310 K. Furthermore, the same residues show the plasticity in the model-free residual dipolar coupling (RDC) order parameters and in an ensemble encoding the supra-τc dynamics. This ß-turn is involved in antibody binding, exhibiting the potential link of the observed supra-τc motion with molecular recognition.

Population shuffling between ground and high energy excited states.

Stochastic processes powered by thermal energy lead to protein motions traversing time-scales from picoseconds to seconds. Fundamental to protein functionality is the utilization of these dynamics for tasks such as catalysis, folding, and allostery. A hierarchy of motion is hypothesized to connect and synergize fast and slow dynamics toward performing these essential activities. Population shuffling predicts a "top-down" temporal hierarchy, where slow time-scale conformational interconversion leads to a shuffling of the free energy landscape for fast time-scale events. Until now, population shuffling was only applied to interconverting ground states. Here, we extend the framework of population shuffling to be applicable for a system interconverting between low energy ground and high energy excited states, such as the SH3 domain mutants G48M and A39V/N53P/V55L from the Fyn tyrosine kinase, providing another tool for accessing the structural dynamics of high energy excited states. Our results indicate that the higher energy gauche - rotameric state for the leucine χ2 dihedral angle contributes significantly to the distribution of rotameric states in both the major and minor forms of the SH3 domain. These findings are corroborated with unrestrained molecular dynamics (MD) simulations on both the major and minor states of the SH3 domain demonstrating high correlations between experimental and back-calculated leucine χ2 rotameric populations. Taken together, we demonstrate how fast time-scale rotameric side-chain population distributions can be extracted from slow time-scale conformational exchange data further extending the scope and the applicability of the population shuffling model.

Sampling of Glycan-Bound Conformers by the Anti-HIV Lectin Oscillatoria agardhii agglutinin in the Absence of Sugar.

Lectins from different sources have been shown to interfere with HIV infection by binding to the sugars of viral-envelope glycoproteins. Three-dimensional atomic structures of a number of HIV-inactivating lectins have been determined, both as free proteins and in glycan-bound forms. However, details on the mechanism of recognition and binding to sugars are elusive. Herein we focus on the anti-HIV lectin OAA from Oscillatoria agardhii: We show that in the absence of sugars in solution, both the sugar-free and sugar-bound protein conformations that were observed in the X-ray crystal structures exist as conformational substates. Our results suggest that glycan recognition occurs by conformational selection within the ground state; this model differs from the popular "excited-state" model. Our findings provide further insight into molecular recognition of the major receptor on the HIV virus by OAA. These details can potentially be used for the optimization and/or development of preventive anti-HIV therapeutics.

A pre-structured helix in the intrinsically disordered 4EBP1.

The eIF4E-binding protein 1 (4EBP1) has long been known to be completely unstructured without any secondary structures, which contributed significantly to the proposal of the induced fit mechanism for target binding of intrinsically disordered proteins. We show here that 4EBP1 is not completely unstructured, but contains a pre-structured helix.

Measuring membrane protein bond orientations in nanodiscs via residual dipolar couplings.

Membrane proteins are involved in numerous vital biological processes. To understand membrane protein functionality, accurate structural information is required. Usually, structure determination and dynamics of membrane proteins are studied in micelles using either solution state NMR or X-ray crystallography. Even though invaluable information has been obtained by this approach, micelles are known to be far from ideal mimics of biological membranes often causing the loss or decrease of membrane protein activity. Recently, nanodiscs, which are composed of a lipid bilayer surrounded by apolipoproteins, have been introduced as a more physiological alternative than micelles for NMR investigations on membrane proteins. Here, we show that membrane protein bond orientations in nanodiscs can be obtained by measuring residual dipolar couplings (RDCs) with the outer membrane protein OmpX embedded in nanodiscs using Pf1 phage as an alignment medium. The presented collection of membrane protein RDCs in nanodiscs represents an important step toward more comprehensive structural and dynamical NMR-based investigations of membrane proteins in a natural bilayer environment.

ORIUM: optimized RDC-based Iterative and Unified Model-free analysis.

Residual dipolar couplings (RDCs) are NMR parameters that provide both structural and dynamic information concerning inter-nuclear vectors, such as N-H(N) and Cα-Hα bonds within the protein backbone. Two approaches for extracting this information from RDCs are the model free analysis (MFA) (Meiler et al. in J Am Chem Soc 123:6098-6107, 2001; Peti et al. in J Am Chem Soc 124:5822-5833, 2002) and the direct interpretation of dipolar couplings (DIDCs) (Tolman in J Am Chem Soc 124:12020-12030, 2002). Both methods have been incorporated into iterative schemes, namely the self-consistent RDC based MFA (SCRM) (Lakomek et al. in J Biomol NMR 41:139-155, 2008) and iterative DIDC (Yao et al. in J Phys Chem B 112:6045-6056, 2008), with the goal of removing the influence of structural noise in the MFA and DIDC formulations. Here, we report a new iterative procedure entitled Optimized RDC-based Iterative and Unified Model-free analysis (ORIUM). ORIUM unifies theoretical concepts developed in the MFA, SCRM, and DIDC methods to construct a computationally less demanding approach to determine these structural and dynamic parameters. In all schemes, dynamic averaging reduces the actual magnitude of the alignment tensors complicating the determination of the absolute values for the generalized order parameters. To readdress this scaling issue that has been previously investigated (Lakomek et al. in J Biomol NMR 41:139-155, 2008; Salmon et al. in Angew Chem Int Edit 48:4154-4157, 2009), a new method is presented using only RDC data to establish a lower bound on protein motion, bypassing the requirement of Lipari-Szabo order parameters. ORIUM and the new scaling procedure are applied to the proteins ubiquitin and the third immunoglobulin domain of protein G (GB3). Our results indicate good agreement with the SCRM and iterative DIDC approaches and signify the general applicability of ORIUM and the proposed scaling for the extraction of inter-nuclear vector structural and dynamic content.

Measuring dynamic and kinetic information in the previously inaccessible supra-τ(c) window of nanoseconds to microseconds by solution NMR spectroscopy.

Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful tool that has enabled experimentalists to characterize molecular dynamics and kinetics spanning a wide range of time-scales from picoseconds to days. This review focuses on addressing the previously inaccessible supra-tc window (defined as τ(c) < supra-τ(c) < 40 µs; in which tc is the overall tumbling time of a molecule) from the perspective of local inter-nuclear vector dynamics extracted from residual dipolar couplings (RDCs) and from the perspective of conformational exchange captured by relaxation dispersion measurements (RD). The goal of the first section is to present a detailed analysis of how to extract protein dynamics encoded in RDCs and how to relate this information to protein functionality within the previously inaccessible supra-τ(c) window. In the second section, the current state of the art for RD is analyzed, as well as the considerable progress toward pushing the sensitivity of RD further into the supra-τ(c) scale by up to a factor of two (motion up to 25 µs). From the data obtained with these techniques and methodology, the importance of the supra-τ(c) scale for protein function and molecular recognition is becoming increasingly clearer as the connection between motion on the supra-τ(c) scale and protein functionality from the experimental side is further strengthened with results from molecular dynamics simulations.

Enhanced accuracy of kinetic information from CT-CPMG experiments by transverse rotating-frame spectroscopy.

Micro-to-millisecond motions of proteins transmit pivotal signals for protein function. A powerful technique for the measurement of these motions is nuclear magnetic resonance spectroscopy. One of the most widely used methodologies for this purpose is the constant-time Carr-Purcell-Meiboom-Gill (CT-CPMG) relaxation dispersion experiment where kinetic and structural information can be obtained at atomic resolution. Extraction of accurate kinetics determined from CT-CPMG data requires refocusing frequencies that are much larger than the nuclei's exchange rate between states. We investigated the effect when fast processes are probed by CT-CPMG experiments via simulation and show that if the intrinsic relaxation rate (R(CT-CPMG)(2,0)) is not known a priori the extraction of accurate kinetics is hindered. Errors on the order of 50 % in the exchange rate are attained when processes become fast, but are minimized to 5 % with a priori (CT-CPMG)(2,0)) information. To alleviate this shortcoming, we developed an experimental scheme probing (CT-CPMG)(2,0)) with large amplitude spin-lock fields, which specifically contains the intrinsic proton longitudinal Eigenrelaxation rate. Our approach was validated with ubiquitin and the Oscillatoria agardhii agglutinin (OAA). For OAA, an underestimation of 66 % in the kinetic rates was observed if (CT-CPMG)(2,0)) is not included during the analysis of CT-CPMG data and result in incorrect kinetics and imprecise amplitude information. This was overcome by combining CT-CPMG with (CT-CPMG)(2,0)) measured with a high power R1ρ experiment. In addition, the measurement of (CT-CPMG)(2,0)) removes the ambiguities in choosing between different models that describe CT-CPMG data.

Ligand binding to anion-binding exosites regulates conformational properties of thrombin.

Thrombin participates in coagulation, anticoagulation, and initiation of platelet activation. To fulfill its diverse roles and maintain hemostasis, this serine protease is regulated via the extended active site region and anion-binding exosites (ABEs) I and II. For the current project, amide proton hydrogen-deuterium exchange coupled with MALDI-TOF mass spectrometry was used to characterize ligand binding to individual exosites and to investigate the presence of exosite-active site and exosite-exosite interactions. PAR3(44-56) and PAR1(49-62) were observed to bind to thrombin ABE I and then to exhibit long range effects over to ABE II. By contrast, Hirudin(54-65) focused more on ABE I and did not transmit influences over to ABE II. Although these three ligands were each directed to ABE I, they did not promote the same conformational consequences. D-Phe-Pro-Arg-chloromethyl ketone inhibition at the thrombin active site led to further local and long range consequences to thrombin-ABE I ligand complexes with the autolysis loop often most affected. When Hirudin(54-65) was bound to ABE I, it was still possible to bind GpIbα(269-286) or fibrinogen γ'(410-427) to ABE II. Each ligand exerted its predominant influences on thrombin and also allowed interexosite communication. The results obtained support the proposal that thrombin is a highly dynamic protein. The transmission of ligand-specific local and long range conformational events is proposed to help regulate this multifunctional enzyme.

Thermal coefficients of the methyl groups within ubiquitin.

Physiological processes such as protein folding and molecular recognition are intricately linked to their dynamic signature, which is reflected in their thermal coefficient. In addition, the local conformational entropy is directly related to the degrees of freedom, which each residue possesses within its conformational space. Therefore, the temperature dependence of the local conformational entropy may provide insight into understanding how local dynamics may affect the stability of proteins. Here, we analyze the temperature dependence of internal methyl group dynamics derived from the cross-correlated relaxation between dipolar couplings of two CH bonds within ubiquitin. Spanning a temperature range from 275 to 308 K, internal methyl group dynamics tend to increase with increasing temperature, which translates to a general increase in local conformational entropy. With this data measured over multiple temperatures, the thermal coefficient of the methyl group order parameter, the characteristic thermal coefficient, and the local heat capacity were obtained. By analyzing the distribution of methyl group thermal coefficients within ubiquitin, we found that the N-terminal region has relatively high thermostability. These results indicate that methyl groups contribute quite appreciably to the total heat capacity of ubiquitin through the regulation of local conformational entropy.