Solvent Interactions and Protein Dynamics in Spin-labeled T4 Lysozyme

Abstract

Aspects of T4 lysozyme dynamics and solvent interaction are investigated using atomically detailed Molecular Dynamics (MD) simulations. Two spin-labeled mutants of T4 lysozyme are analyzed (T4L-N40C and T4L-K48C), which have been found from electronic paramagnetic resonance (EPR) experiments to exhibit different mobilities at the site of spin probe attachment (N- and C-terminus of helix B, respectively). Similarities and differences in solvent distribution and diffusion around the spin label, as well as around exposed and buried residues within the protein, are discussed. The purpose is to capture possible strong interactions between the spin label (ring) and solvent molecules, which may affect EPR lineshapes. The effect of backbone motions on the water density profiles is also investigated. The focus is on the domain closure associated with the T4 lysozyme hinge-bending motion, which is analyzed by Essential Dynamics (ED). The N-terminus of helix B is found to be a “hinge” residue, which explains the high degree of flexibility and motional freedom at this site.     

NMR Determines Transient Structure and Dynamics in the Disordered C-Terminal Domain of WASp Interacting Protein

WASp-interacting protein (WIP) is a 503-residue proline-rich polypeptide expressed in human T cells. The WIP C-terminal domain binds to Wiskott-Aldrich syndrome protein (WASp) and regulates its activation and degradation, and the WIP-WASp interaction has been shown to be critical for actin polymerization and implicated in the onset of WAS and X-linked thrombocytopenia. WIP is predicted to be an intrinsically disordered protein, a class of polypeptides that are of great interest because they violate the traditional structure-function paradigm. In this first (to our knowledge) study of WIP in its unbound state, we used NMR to investigate the biophysical behavior of WIP^C, a C-terminal domain fragment of WIP that includes residues 407–503 and contains the WASp-binding site. In light of the poor spectral dispersion exhibited by WIP^C and the high occurrence (25%) of proline residues, we employed 5D-NMR¹³C-detected NMR experiments with nonuniform sampling to accomplish full resonance assignment. Secondary chemical-shift analysis, ¹⁵N relaxation rates, and protection from solvent exchange all concurred in detecting transient structure located in motifs that span the WASp-binding site. Residues 446–456 exhibited a propensity for helical conformation, and an extended conformation followed by a short, capped helix was observed for residues 468–478. The ¹³C-detected approach allows chemical-shift assignment in the WIP^C polyproline stretches and thus sheds light on their conformation and dynamics. The effects of temperature on chemical shifts referenced to a denatured sample of the polypeptide demonstrate that heating reduces the structural character of WIP^C. Thus, we conclude that the disordered WIP^C fragment is comprised of regions with latent structure connected by flexible loops, an architecture with implications for binding affinity and function.     

Sampling of Protein Dynamics in Nanosecond Time Scale by 15N NMR Relaxation and Self-Diffusion Measurements

Abstract

This paper presents a procedure for detection of intermediate nanosecond internal dynamics in globular proteins. The procedure uses ¹H-¹⁵N relaxation measurements at several spectrometer frequencies and hydrodynamic calculations based on experimental self-diffusion coefficients. New heteronuclear experiments, using pulse field gradients, are introduced for the measurement of translation diffusion coefficients of ¹⁵N labeled proteins. An advanced interpretation of recently published (Luginbühl et al., Biochemistry, 36, 7305–7312 (1997)) backbone amide ¹⁵N relaxation data, measured at two spectrometers (400 and 750 MHz for ¹H) for N-terminal DNA-binding domain (1–63) of 434 repressor, is presented. Non-applicability of commonly used fast (picosecond) dynamics model (FD) was justified by (i) poor fit of relaxation data by the FD model-free spectral density function both for isotropic and anisotropic models of the overall molecular tumbling; (ii) specific dependence of the overall rotation correlation times calculated from T₁/T₂ ratio on the spectrometer frequency; (iii) mismatch of the ratio of longitudinal ¹⁵N relaxation times T₁, measured at different spectrometer frequencies, in comparison with that anticipated for the IT) model; (iv) significantly underestimated overall rotation correlation time provided by the FD model (5.50±0.15 and 5.80±0.15 ns for 750 and 400 MHz spectrometer frequency respectively) in comparison with correlation time obtained from hydrodynamics. On the other hand, all relaxation and hydrodynamics data are in good correspondence with the model of intermediate (nanoseconds) dynamics. Overall rotation correlation time of 7.5±0.7 ns was calculated from experimental translation self-diffusion rate using hydrodynamics formalism (Garcia de la Torre, J. and Bloomfield, V.A. Quart. Rev. Biophys., 14, 81–139 (1981)). The statistical analysis of ¹⁵N relaxation data along with the hydrodynamic consideration clearly revealed that most of the residues in 434(1–63) repressor are involved in the nanosecond internal dynamics characterized by the the mean order parameters of 0.59±0.06 and the correlation times of ca. 5 ns.     

NMR Structure of the SARS-CoV Nonstructural Protein 7 in Solution at pH 6.5

The NMR structure of the severe acute respiratory syndrome coronavirus nonstructural protein (nsp) 7 in aqueous solution at pH 6.5 was determined and compared with the results of previous structure determinations of nsp7 in solution at pH 7.5 and in the crystals of a hexadecameric nsp7/nsp8 complex obtained from a solution at pH 7.5. All three structures contain four helices as the only regular secondary structures, but there are differences in the lengths and sequence locations of the four helices, as well as between the tertiary folds. The present study includes data on conformational equilibria and intramolecular rate processes in nsp7 in solution at pH 6.5, which provide further insights into the polymorphisms implicated by a comparison of the three presently available nsp7 structures.     

Integrating NMR,SAXS, and Atomistic Simulations: Structure and Dynamics of a Two-Domain Protein

Multidomain proteins with two or more independently folded functional domains are prevalent in nature. Whereas most multidomain proteins are linked linearly in sequence, roughly one-tenth possess domain insertions where a guest domain is implanted into a loop of a host domain, such that the two domains are connected by a pair of interdomain linkers. Here, we characterized the influence of the interdomain linkers on the structure and dynamics of a domain-insertion protein in which the guest LysM domain is inserted into a central loop of the host CVNH domain. Expanding upon our previous crystallographic and NMR studies, we applied SAXS in combination with NMR paramagnetic relaxation enhancement to construct a structural model of the overall two-domain system. Although the two domains have no fixed relative orientation, certain orientations were found to be preferred over others. We also assessed the accuracies of molecular mechanics force fields in modeling the structure and dynamics of tethered multidomain proteins by integrating our experimental results with microsecond-scale atomistic molecular dynamics simulations. In particular, our evaluation of two different combinations of the latest force fields and water models revealed that both combinations accurately reproduce certain structural and dynamical properties, but are inaccurate for others. Overall, our study illustrates the value of integrating experimental NMR and SAXS studies with long timescale atomistic simulations for characterizing structural ensembles of flexibly linked multidomain systems.     

The Role of Solvent in Protein Folding and in Aggregation

We discuss features of the effect of solvent on protein folding andaggregation, highlighting the physics related to the particulate nature and the peculiar structure of the aqueous solvent, and the biological significance of interactions between solvent and proteins. To this purpose we use a generalized energy landscape of extended dimensionality. A closer look at the properties of solvent induced interactions and forces proves useful for understanding the physical grounds of `ad hoc' interactions and for devising realistic ways of accounting for solvent effects. The solvent has long been known to be a crucially important part of biological systems, and times appear mature for it to be adequately accounted for in the protein folding problem. Use of the extended dimensionality energy landscape helpseliciting the possibility of coupling among conformational changes and aggregation, such as proved by experimental data in the literature.     

The Role of Solvent Interactions in Protein Conformatio

Abstract

For many years, the ease with which the pH of enzyme reactions could be varied has led to innumerable reports of the pH-dependence of the kinetic parameters of such reactions. These studies have been provided with a veneer of respectability by suggestions that they allow deductions to be made about the nature of the ionizing groups upon which the catalytic activity of the enzyme depends. It is the purpose of this article to summarize some of the reasons why a large proportion of pH-dependence studies yields little information about the real pK_a-values of ionizing groups in enzymes, and still less about the identity of these groups.     

Role of Megafauna and Frozen Soil in the Atmospheric CH4 Dynamics

Modern wetlands are the world’s strongest methane source. But what was the role of this source in the past? An analysis of global ¹⁴C data for basal peat combined with modelling of wetland succession allowed us to reconstruct the dynamics of global wetland methane emission through time. These data show that the rise of atmospheric methane concentrations during the Pleistocene-Holocene transition was not connected with wetland expansion, but rather started substantially later, only 9 thousand years ago. Additionally, wetland expansion took place against the background of a decline in atmospheric methane concentration. The isotopic composition of methane varies according to source. Owing to ice sheet drilling programs past dynamics of atmospheric methane isotopic composition is now known. For example over the course of Pleistocene-Holocene transition atmospheric methane became depleted in the deuterium isotope, which indicated that the rise in methane concentrations was not connected with activation of the deuterium-rich gas clathrates. Modelling of the budget of the atmospheric methane and its isotopic composition allowed us to reconstruct the dynamics of all main methane sources. For the late Pleistocene, the largest methane source was megaherbivores, whose total biomass is estimated to have exceeded that of present-day humans and domestic animals. This corresponds with our independent estimates of herbivore density on the pastures of the late Pleistocene based on herbivore skeleton density in the permafrost. During deglaciation, the largest methane emissions originated from degrading frozen soils of the mammoth steppe biome. Methane from this source is unique, as it is depleted of all isotopes. We estimated that over the entire course of deglaciation (15,000 to 6,000 year before present), soils of the mammoth steppe released 300–550 Pg (10¹⁵ g) of methane. From current study we conclude that the Late Quaternary Extinction significantly affected the global methane cycle.     

Protein dynamics from NMR

This review surveys recent investigations of conformational fluctuations of proteins in solution using NMR techniques. Advances in experimental methods have provided more accurate means of characterizing fast and slow internal motions as well as overall diffusion. The information obtained from NMR dynamics experiments provides insights into specific structural changes or configurational energetics associated with function. A variety of applications illustrate that studies of protein dynamics provide insights into protein-protein interactions, target recognition, ligand binding, and enzyme function.     

Thermodynamic Framework of the Interaction between Protein and Solvent Drives Protein Folding

Abstract

We use internal coordinate molecular mechanics calculations to study the impact of abasic sites on the conformation and the mechanics of the DNA double helix. Abasic sites, which are common mutagenic lesions, are shown to locally modify both the groove geometry and the curvature of DNA in a sequence dependent manner. By controlled twisting and bending, it is also shown that these lesions modify the deformability of the duplex, generally increasing its flexibility, but again to an extent which depends on the nature of the abasic site and on the surrounding base sequence. Both the conformational and mechanical influence of this type of DNA damage may be significant for recognition and repair mechanisms.     

Structure and Dynamics of Ribosomal Protein L12: An Ensemble Model Based on SAXS and NMR Relaxation

Ribosomal protein L12 is a two-domain protein that forms dimers mediated by its N-terminal domains. A 20-residue linker separates the N- and C-terminal domains. This linker results in a three-lobe topology with significant flexibility, known to be critical for efficient translation. Here we present an ensemble model of spatial distributions and correlation times for the domain reorientations of L12 that reconciles experimental data from small-angle x-ray scattering and nuclear magnetic resonance. We generated an ensemble of L12 conformations in which the structure of each domain is fixed but the domain orientations are variable. The ensemble reproduces the small-angle x-ray scattering data and the optimized correlation times of its reorientational eigenmodes fit the ¹⁵N relaxation data. The ensemble model reveals intrinsic conformational properties of L12 that help explain its function on the ribosome. The two C-terminal domains sample a large volume and extend further away from the ribosome anchor than expected for a random-chain linker, indicating that the flexible linker has residual order. Furthermore, the distances between each C-terminal domain and the anchor are anticorrelated, indicating that one of them is more retracted on average. We speculate that these properties promote the function of L12 to recruit translation factors and control their activity on the ribosome.     

Solution NMR conformation of glycosaminoglycans

Nuclear magnetic resonance (NMR) spectroscopy has been giving a pivotal contribution to the progress of glycomics, mostly by elucidating the structural, dynamical, conformational and intermolecular binding aspects of carbohydrates. Particularly in the field of conformation, NOE resonances, scalar couplings, residual dipolar couplings, and chemical shift anisotropy offsets have been the principal NMR parameters utilized. Molecular dynamics calculations restrained by NMR-data input are usually employed in conjunction to generate glycosidic bond dihedral angles. Glycosaminoglycans (GAGs) are a special class of sulfated polysaccharides extensively studied worldwide. Besides regulating innumerous physiological processes, these glycans are also widely explored in the global market as either clinical or nutraceutical agents. The conformational aspects of GAGs are key regulators to the quality of interactions with the functional proteins involved in biological events. This report discusses the solution conformation of each GAG type analyzed by one or more of the above-mentioned methods.