Conserved allosteric perturbation of the GTPase domains by region 1 of Ras hypervariable regions.

Ras proteins are important intracellular signaling hubs that can interact with numerous downstream effectors and upstream regulators through their GTPase domains (G-domains) anchored to plasma membranes by the C-terminal hypervariable regions (HVRs). The biological functions of Ras were proposed to be regulated at multiple levels including the intramolecular G-domain-HVR interactions, of which the exact mechanism and specificity are still controversial. Here, we demonstrate that the HVRs, instead of having direct contacts, can weakly perturb the G-domains via an allosteric interaction that is restricted to a ∼20 Å range and highly conserved in the tested Ras isoforms (HRas and KRas4B) and nucleotide-bound states. The origin of this allosteric perturbation has been localized to a short segment (residues 167-171) coinciding with region 1 of HVRs, which exhibits moderate to weak α-helical propensities. A charge-reversal mutation (E168K) of KRas4B in region 1, previously described in the Catalog of Somatic Mutations in Cancer database, was found to induce similar chemical shift perturbations as truncation of the HVR does. Further membrane paramagnetic relaxation enhancement (mPRE) data show that this region 1 mutation alters the membrane orientations of KRas4B and moderately increases the relative population of the signaling-compatible state.

Quantitative Ensemble Interpretation of Membrane Paramagnetic Relaxation Enhancement (mPRE) for Studying Membrane-Associated Intrinsically Disordered Proteins.

An understanding of the functional role played by a membrane-associated intrinsically disordered protein (IDP) requires characterization of its heterogeneous conformations as well as its poses relative to the membranes, which is of great interest but technically challenging. Here, we explore the membrane paramagnetic relaxation enhancement (mPRE) for constructing ensembles of IDPs that dynamically associate with membrane mimetics incorporating spin-labeled lipids. To accurately interpret the mPRE Γ2 rates, both the dynamics of IDPs and spin probe molecules are taken into account, with the latter described by a weighted three-dimensional (3D) grid model built based on all-atom simulations. The IDP internal conformations, orientations, and immersion depths in lipid bilayers are comprehensively optimized in the Γ2-based ensemble modeling. Our approach is tested and validated on the example of POPG bicelle-bound disordered cytoplasmic domain of CD3ε (CD3εCD), a component of the T-cell receptor (TCR) complex. The mPRE-derived CD3εCD ensemble provides new insights into the IDP-membrane fuzzy association, in particular for the tyrosine-based signaling motif that plays a critical role in TCR signaling. The comparative analysis of the ensembles for wild-type CD3εCD and mutants that mimic the mono- and dual-phosphorylation effects suggests a delicate membrane regulatory mechanism for activation and inhibition of the TCR activity.

Unveiling the "invisible" druggable conformations of GDP-bound inactive Ras.

The prevalent view on whether Ras is druggable has gradually changed in the recent decade with the discovery of effective inhibitors binding to cryptic sites unseen in the native structures. Despite the promising advances, therapeutics development toward higher potency and specificity is challenged by the elusive nature of these binding pockets. Here we derive a conformational ensemble of guanosine diphosphate (GDP)-bound inactive Ras by integrating spin relaxation-validated atomistic simulation with NMR chemical shifts and residual dipolar couplings, which provides a quantitative delineation of the intrinsic dynamics up to the microsecond timescale. The experimentally informed ensemble unequivocally demonstrates the preformation of both surface-exposed and buried cryptic sites in Ras•GDP, advocating design of inhibition by targeting the transient druggable conformers that are invisible to conventional experimental methods. The viability of the ensemble-based rational design has been established by retrospective testing of the ability of the Ras•GDP ensemble to identify known ligands from decoys in virtual screening.

Quantitative Paramagnetic NMR-Based Analysis of Protein Orientational Dynamics on Membranes: Dissecting the KRas4B-Membrane Interactions.

Peripheral membrane proteins can adopt distinct orientations on the surfaces of lipid bilayers that are often short-lived and challenging to characterize by conventional experimental methods. Here we describe a robust approach for mapping protein orientational landscapes through quantitative interpretation of paramagnetic relaxation enhancement (PRE) data arising from membrane mimetics with spin-labeled lipids. Theoretical analysis, followed by experimental verification, reveals insights into the distinct properties of the PRE observables that are generally distorted in the case of stably membrane-anchored proteins. To suppress the artifacts, we demonstrate that undistorted Γ2 values can be obtained via transient membrane anchoring, based on which a computational framework is established for deriving accurate orientational ensembles obeying Boltzmann statistics. Application of the approach to KRas4B, a classical peripheral membrane protein whose orientations are critical for its functions and drug design, reveals four distinct orientational states that are close but not identical to those reported previously. Similar orientations are also found for a truncated KRas4B without the hypervariable region (HVR) that can sample a broader range of orientations, suggesting a confinement role of the HVR geometrically prohibiting severe tilting. Comparison of the KRas4B Γ2 rates measured using nanodiscs containing different types of anionic lipids reveals identical Γ2 patterns for the G-domain but different ones for the HVR, indicating only the latter is able to selectively interact with anionic lipids.

Exploring the state- and allele-specific conformational landscapes of Ras: understanding their respective druggabilities.

Recent advances in direct inhibition of Ras benefit from the protein's intrinsic dynamic nature that derives therapeutically vulnerable conformers bearing transiently formed cryptic pockets. Hotspot mutants of Ras are major tumor drivers and are hyperactivated in cells at variable levels, which may require allele-specific strategies for effective targeting. However, it remains unclear how the prevalent oncogenic mutations and activation states perturb the free energy landscape governing the protein dynamics and druggability. Here we characterized the nucleotide state- and allele-dependent alterations of Ras conformational dynamics using a combined NMR experimental and computational approach and constructed quantitative ensembles revealing the conservation of the cryptic SI/II-P and SII-P pockets in different states and alleles. Highly local but critical conformational reorganizations that undermine the SII-P accessibility to residue 12 have been identified as a common mechanism resulting in the low reactivities of Ras·GTP as well as Ras(G12D)·GDP with covalent SII-P inhibitors. Our results strongly support the conformational selection scenario for interactions between Ras and the previously reported binders and offer insights for the future development of state- and allele-specific, as well as pan-Ras, inhibitors.

Millisecond Allosteric Dynamics of Activated Ras Reproduced with a Slowly Hydrolyzable GTP Analogue.

The millisecond timescale dynamics of activated Ras transiently sample a low-populated conformational state that has distinct surface property from the major state and represents a promising target for binding of small-molecule compounds. To avoid the complications of hydrolysis, dynamics and other properties of active Ras have so far been routinely investigated by using non-hydrolyzable GTP analogues, which, however, were previously reported to alter both the kinetics and distribution of the conformational exchange. In this study, we quantitatively measured and validated the internal dynamics of Ras complexed with a slowly hydrolyzable GTP analogue, GTPγS, which increases the lifetime of active Ras by 23 times relative to that of native GTP. It was found that GTPγS, in addition to its better mimicking of the exchange kinetics than the commonly used non-hydrolyzable analogues GppNHp and GppCH2 p, can rigorously reproduce the natural dynamics network in active Ras, thus indicating its fitness for use in the development of allosteric inhibitors.

Atomistic Insights into the Functional Instability of the Second Helix of Fatty Acid Binding Protein.

Structural dynamics of fatty acid binding proteins (FABPs), which accommodate poorly soluble ligands in the internalized binding cavities, are intimately related to their function. Recently, local unfolding of the α-helical cap in a variant of human intestinal FABP (IFABP) has been shown to correlate with the kinetics of ligand association, shedding light on the nature of the critical conformational reorganization. Yet, the physical origin and mechanism of the functionally relevant transient unfolding remain elusive. Here, we investigate the intrinsic structural instability of the second helix (αII) of IFABP in comparison with other segments of the protein using hydrogen-exchange NMR spectroscopy, microsecond molecular dynamics simulations, and enhanced sampling techniques. Although tertiary interactions positively contribute to the stability of helices in IFABP, the intrinsic unfolding tendency of αII is encoded in its primary sequence and can be described by the Lifson-Roig theory in the absence of tertiary interactions. The unfolding pathway of αII in intact proteins involves an on-pathway intermediate state that is characterized with the fraying of the last helical turn, captured by independent enhanced sampling methods. The simulations in this work, combined with hydrogen-exchange NMR data, provide new, to our knowledge, atomistic insights into the functional local unfolding of FABPs.

Effects of rumen cannulation on dissolved gases and methanogen community in dairy cows.

Rumen cannulation is a widely employed technique in ruminant nutrition research. However, the gap between skin and rumen cannula can cause leakage of fermentation gases and influx of atmospheric air, which may adversely affect the anaerobic environment in the rumen. The present study was designed to investigate the effects of rumen cannulation on headspace gases, dissolved gases, fermentation end products, and methanogen community in the rumen of dairy cows. Eight Holstein cows were used in the experiment. Four cows were surgically fitted with rumen cannulas, whereas the other 4 intact cows were used as control. Rumen cannulation decreased gaseous hydrogen and methane concentrations, dissolved carbon dioxide concentration, and relative abundances of Methanosphaera, and increased the saturation factor of dissolved hydrogen and dissolved methane, dissolved methane concentration, volatile fatty acid concentration, 16S ribosomal RNA gene copies of methanogens, and Simpson index of methanogen community. In summary, rumen cannulation causes a reduction in headspace gaseous hydrogen and gaseous methane, which may not decrease dissolved gas concentrations due to an increase in saturation factors. Furthermore, rumen cannulation alters methanogen community with increased methanogen population and decreased relative abundances of Methanosphaera.

Extending the Lifetime of Native GTP-Bound Ras for Site-Resolved NMR Measurements: Quantifying the Allosteric Dynamics.

Characterization of native GTP-bound Ras is important for an appreciation of its cellular signaling and for the design of inhibitors, which however has been depressed by its intrinsic instability. Herein, an effective approach for extending the lifetime of Ras⋅GTP samples by exploiting the active role of Son of Sevenless (Sos) is demonstrated that sustains the activated state of Ras. This approach, combined with a postprocessing method that suppresses residual Ras⋅GDP signals, is applied to the site-resolved NMR measurement of the allosteric dynamics of Ras⋅GTP. The observed network of concerted motions well covers the recently identified allosteric inhibitor-binding pockets, but the motions are more confined than those of Ras⋅GppNHp, advocating the use of native GTP for development of allosteric inhibitors. The Sos-based approach is anticipated to generally facilitate experiments on active Ras when native GTP is preferred.

Characterization of the interaction interface and conformational dynamics of human TGIF1 homeodomain upon the binding of consensus DNA.

The TG interacting factor-1 homeodomain (TGIF1-HD) binds with the consensus DNA motif 5'-TGTCA-3' in gene promoters through its three-amino acid loop extension (TALE) type homeodomain, and then recruits co-regulators to regulate gene expression. Although the solution NMR structure of human TGIF1-HD has been reported previously, little is known about its DNA binding mechanism. NMR titrations have been extensively used to study mechanisms of ligand binding to target proteins; however, an intermediate exchange occurred predominantly between TGIF1-HD in the free and bound states when titrated with the consensus DNA, which resulted in poor-quality NMR spectra and precluded further exploration of its interaction interface and conformational dynamics. Here, the helix α3 of TGIF1-HD was speculated as the specific DNA binding interface by hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments, and subsequently confirmed by chemical exchange saturation transfer (CEST) spectroscopy. In addition, simultaneous conformational changes in other regions, including α1 and α2, were induced by DNA binding, explaining the observation of chemical shift perturbations from extensive residues besides those located in α3. Further, low-populated DNA-bound TGIF1-HD undergoing a slow exchange at a rate of 130.2 ±â€¯3.6 s-1 was derived from the analysis of the CEST data, and two residues, R220 and R221, located in the middle of α3 were identified to be crucial for DNA binding. Our study provides structural and dynamic insights into the mechanisms of TGIF1-HD recognition of extensive promoter DNA.

Nitrate improves ammonia incorporation into rumen microbial protein in lactating dairy cows fed a low-protein diet.

Generation of ammonia from nitrate reduction is slower compared with urea hydrolysis and may be more efficiently incorporated into ruminal microbial protein. We hypothesized that nitrate supplementation could increase ammonia incorporation into microbial protein in the rumen compared with urea supplementation of a low-protein diet fed to lactating dairy cows. Eight multiparous Chinese Holstein dairy cows were used in a crossover design to investigate the effect of nitrate or an isonitrogenous urea inclusion in the basal low-protein diet on rumen fermentation, milk yield, and ruminal microbial community in dairy cows fed a low-protein diet in comparison with an isonitrogenous urea control. Eight lactating cows were blocked in 4 pairs according to days in milk, parity, and milk yield and allocated to urea (7.0 g urea/kg of dry matter of basal diet) or nitrate (14.6 g of NO3-/kg of dry matter of basal diet, supplemented as sodium nitrate) treatments, which were formulated on 75% of metabolizable protein requirements. Nitrate supplementation decreased ammonia concentration in the rumen liquids (-33.1%) and plasma (-30.6%) as well as methane emissions (-15.0%) and increased dissolved hydrogen concentration (102%), microbial N (22.8%), propionate molar percentage, milk yield, and 16S rRNA gene copies of Selenomonas ruminantium. Ruminal dissolved hydrogen was positively correlated with the molar proportion of propionate (r = 0.57), and negatively correlated with acetate-to-propionate ratio (r = -0.57) and estimated net metabolic hydrogen production relative to total VFA produced (r = -0.58). Nitrate reduction to ammonia redirected metabolic hydrogen away from methanogenesis, enhanced ammonia incorporation into rumen microbial protein, and shifted fermentation from acetate to propionate, along with increasing S. ruminantium 16S rRNA gene copies, likely leading to the increased milk yield.

Measuring hydrogen exchange rates in invisible protein excited states.

Hydrogen exchange rates have become a valuable probe for studying the relationship between dynamics and structure and for dissecting the mechanism by which proteins fold to their native conformation. Typically measured rates correspond to averages over all protein states from which hydrogen exchange can occur. Here we describe a new NMR experiment based on chemical exchange saturation transfer that provides an avenue for obtaining uncontaminated, per-residue amide hydrogen exchange rates for interconverting native and invisible states so long as they can be separated on the basis of distinct (15)N chemical shifts. The approach is applied to the folding reaction of the Fyn SH3 domain that exchanges between a highly populated, NMR-visible native state and a conformationally excited, NMR-invisible state, corresponding to the unfolded ensemble. Excellent agreement between experimentally derived hydrogen exchange rates of the excited state at a pair of pHs is obtained, taking into account the expected dependence of exchange on pH. Extracted rates for the unfolded ensemble have been used to test hydrogen exchange predictions based on the primary protein sequence that are used in many analyses of solvent exchange rates, with a Pearson correlation coefficient of 0.84 obtained.

Late-Onset Nail Changes Associated with Hand, Foot, and Mouth Disease: A Clinical Analysis of 56 Cases.

OBJECTIVES: To investigate the development and clinical characteristics of nail changes in hand, foot, and mouth disease (HFMD). METHODS: A telephone survey was conducted with the parents of patients diagnosed with HFMD in the Fourth General Hospital of Nanhai from June to August 2013 to document nail changes within 3 months of diagnosis of HFMD. RESULTS: Valid survey results were obtained from 273 cases. Definitive nail changes were identified in 56 patients (20.5%). More boys (25.8%) than girls (10.6%) (p < 0.01) showed changes. The age distribution ranged from 1 to 5 years, and nail changes were rare in children younger than 1 year of age (p < 0.01). Nail changes were usually seen 1 to 2 months after the onset of HFMD and lasted for 1 to 8 weeks, most for approximately 4 weeks. Toenails or fingernails could be affected and the changes were more likely to occur synchronously. Fingernails were more commonly involved than toenails. When both fingernails and toenails were involved, this typically occurred synchronously. Although there were cases with all toenails and fingernails involved (16.1%), we did not encounter any instances involving 13 to 19 nails. The nail changes mainly presented as onychomadesis. Spontaneous recovery without special treatment was the course for all patients. No relapse or new nail involvement was identified. CONCLUSIONS: Nail change associated with HFMD usually occurs within 1 to 2 months after onset, mainly presents as onychomadesis, and is a self-limited process. Possible mechanisms are discussed.

Microsecond Timescale Dynamics of GDP-Bound Ras Underlies the Formation of Novel Inhibitor-Binding Pockets.

The recent discovery of inhibitory compounds binding to distinct pockets on GDP-bound Ras has renewed the view on the druggability of this crucial cancer driver. However, the origin of these pockets, which are not readily formed in the crystal structure in the absence of the compounds, is yet unclear. Herein, we explored the intrinsic flexibility of Ras⋅GDP on microsecond to millisecond timescales using relaxation-based NMR experiments, and identified substantial slow dynamics with τex of 34 µs at 5 °C, which maps to the regions showing a high level of correlation with those displaying conformational differences between the inhibitor-bound and free states. These findings, which have been demonstrated in both wild-type Ras and the oncogenic mutant (G12V), support the mechanism of extended conformational selection for Ras-inhibitor interactions where the small molecules redistribute the protein conformational ensemble favoring the final bound states.

Probing Invisible, Excited Protein States by Non-Uniformly Sampled Pseudo-4D CEST Spectroscopy.

Chemical exchange saturation transfer (CEST) NMR spectroscopy is a powerful tool for studies of slow timescale protein dynamics. Typical experiments are based on recording a large number of 2D data sets and quantifying peak intensities in each of the resulting planes. A weakness of the method is that peaks must be resolved in 2D spectra, limiting applications to relatively small proteins. Resolution is significantly improved in 3D spectra but recording uniformly sampled data is time-prohibitive. Here we describe non-uniformly sampled HNCO-based pseudo-4D CEST that provides excellent resolution in reasonable measurement times. Data analysis is done through fitting in the time domain, without the need of reconstructing the frequency dimensions, exploiting previously measured accurate peak positions in reference spectra. The methodology is demonstrated on several protein systems, including a nascent form of superoxide dismutase that is implicated in neurodegenerative disease.

Triple resonance-based ¹³C(α) and ¹³C(ß) CEST experiments for studies of ms timescale dynamics in proteins.

A pair of triple resonance based CEST pulse schemes are presented for measuring ¹³C(α) and ¹³C(ß) chemical shifts of sparsely populated and transiently formed conformers that are invisible to traditional NMR experiments. CEST profiles containing dips at resonance positions of ¹³C(α) or ¹³C(ß) spins of major (ground) and minor (excited) conformers are obtained in a pseudo 3rd dimension that is generated by quantifying modulations of cross peaks in ¹5N, ¹H(N) correlation spectra. An application to the folding reaction of a G48A mutant of the Fyn SH3 domain is presented, illustrating and validating the methodology.

In silico elucidation of the recognition dynamics of ubiquitin.

Elucidation of the mechanism of biomacromolecular recognition events has been a topic of intense interest over the past century. The inherent dynamic nature of both protein and ligand molecules along with the continuous reshaping of the energy landscape during the binding process renders it difficult to characterize this process at atomic detail. Here, we investigate the recognition dynamics of ubiquitin via microsecond all-atom molecular dynamics simulation providing both thermodynamic and kinetic information. The high-level of consistency found with respect to experimental NMR data lends support to the accuracy of the in silico representation of the conformational substates and their interconversions of free ubiquitin. Using an energy-based reweighting approach, the statistical distribution of conformational states of ubiquitin is monitored as a function of the distance between ubiquitin and its binding partner Hrs-UIM. It is found that extensive and dense sampling of conformational space afforded by the µs MD trajectory is essential for the elucidation of the binding mechanism as is Boltzmann sampling, overcoming inherent limitations of sparsely sampled empirical ensembles. The results reveal a population redistribution mechanism that takes effect when the ligand is at intermediate range of 1-2 nm from ubiquitin. This mechanism, which may be depicted as a superposition of the conformational selection and induced fit mechanisms, also applies to other binding partners of ubiquitin, such as the GGA3 GAT domain.

Toward a predictive understanding of slow methyl group dynamics in proteins.

The development of the most recent generation of molecular mechanics force fields promises an increasingly predictive understanding of the protein dynamics-function relationship. Based on extensive validation against various types of experimental data, the AMBER force field ff99SB was benchmarked in recent years as a favorable force field for protein simulations. Recent improvements of the side chain and backbone potentials, made by different groups, led to the ff99SB-ILDN and ff99SBnmr1 force fields, respectively. The combination of these potentials into a unified force field, termed ff99SBnmr1-ILDN, was used in this study to perform a microsecond time scale molecular dynamics simulation of free ubiquitin in explicit solvent for validation against an extensive set of experimental NMR methyl group residual dipolar couplings. Our results show a high level of consistency between the experimental data and the values predicted from the molecular dynamics trajectory reflecting a systematically improved performance of ff99SBnmr1-ILDN over the original ff99SB force field. Moreover, the unconstrained ff99SBnmr1-ILDN MD ensemble achieves a similar level of agreement as the recently introduced EROS ensemble, which was constructed based on a large body of NMR data as constraints, including the methyl residual dipolar couplings. This suggests that ff99SBnmr1-ILDN provides a high-quality representation of the motions of methyl-bearing protein side chains, which are sensitive probes of protein-protein and protein-ligand interactions.

Atomistic kinetic model for population shift and allostery in biomolecules.

Allosteric signaling in biomolecules is a key mechanism for a myriad of cellular processes. We present a general yet compact model for protein allostery at atomic detail to quantitatively explain and predict structural-dynamics properties of allosteric signal propagation. The master equation-based approach for allostery by population shift (MAPS) is introduced that derives the time scales, amplitudes, and pathways of signal transmission in peptides and proteins from dihedral angle dynamics observed in extended molecular dynamics simulations. The MAPS approach is first applied to the alanine-pentapeptide, and the results are tested against an explicit simulation in the presence of local conformational constraints, confirming the validity and accuracy of the model. We then apply the approach to a larger Markovian system based on a millisecond all-atom protein molecular dynamics trajectory of BPTI (Shaw et al. Science 2010, 330, 341-346). We use MAPS to illustrate in silico the propagation of a local perturbation over medium- to long-range distances across a disulfide bridge linking loops L1 and L2, which constitute the binding interface of BPTI.

Conversion of hexagonal Sb2Te3 nanoplates into nanorings driven by growth temperature.

We describe a novel route for the conversion of hexagonal Sb(2)Te(3) nanoplates into nanorings driven by growth temperature in a simple solvothermal process. The transmission electron microscopy was employed to investigate systemically the morphology, size, crystallinity, and microstructure of the as-prepared products. The experiments indicated that the growth temperature had a great effect on the morphology of antimony telluride nanostructures. When the experiments were conducted at 200 °C, the hexagonal antimony telluride nanoplates were obtained. However, if the experiments were carried out at higher temperature of 230 °C, the hexagonal antimony telluride nanorings were achieved by dissolution of the inner part with a higher density of defects of the hexagonal nanoplates for the first time. A possible formation mechanism was proposed on the basis of experimental results and analysis. This work may open a new rational route for the synthesis of the hexagonal antimony telluride nanorings, which may have scientific and technological applications in various functional devices.