Water molecule ordering on the surface of an intrinsically disordered protein.

The dynamics and local structure of the hydration water on surfaces of folded proteins have been extensively investigated. However, our knowledge of the hydration of intrinsically disordered proteins (IDPs) is more limited. Here, we compare the local structure of water molecules hydrating a globular protein, lysozyme, and the intrinsically disordered N-terminal of c-Src kinase (SH4UD) using molecular dynamics simulation. The radial distributions from the protein surface of the first and the second hydration shells are similar for the folded protein and the IDP. However, water molecules in the first hydration shell of both the folded protein and the IDP are perturbed from the bulk. This perturbation involves a loss of tetrahedrality, which is, however, significantly more marked for the folded protein than the IDP. This difference arises from an increase in the first hydration shell of the IDP of the fraction of hydration water molecules interacting with oxygen. The water ordering is independent of the compactness of the IDP. In contrast, the lifetimes of water molecules in the first hydration shell increase with IDP compactness, indicating a significant impact of IDP configuration on water surface pocket kinetics, which here is linked to differential pocket volumes and polarities.

Disordered Domain Shifts the Conformational Ensemble of the Folded Regulatory Domain of the Multidomain Oncoprotein c-Src.

c-Src kinase is a multidomain non-receptor tyrosine kinase that aberrantly phosphorylates several signaling proteins in cancers. Although the structural properties of the regulatory domains (SH3-SH2) and the catalytic kinase domain have been extensively characterized, there is less knowledge about the N-terminal disordered region (SH4UD) and its interactions with the other c-Src domains. Here, we used domain-selective isotopic labeling combined with the small-angle neutron scattering contrast matching technique to study SH4UD interactions with SH3-SH2. Our results show that in the presence of SH4UD, the radius of gyration (Rg) of SH3-SH2 increases, indicating that it has a more extended conformation. Hamiltonian replica exchange molecular dynamics simulations provide a detailed molecular description of the structural changes in SH4UD-SH3-SH2 and show that the regulatory loops of SH3 undergo significant conformational changes in the presence of SH4UD, while SH2 remains largely unchanged. Overall, this study highlights how a disordered region can drive a folded region of a multidomain protein to become flexible, which may be important for allosteric interactions with binding partners. This may help in the design of therapeutic interventions that target the regulatory domains of this important family of kinases.

Predicting the structural basis of targeted protein degradation by integrating molecular dynamics simulations with structural mass spectrometry.

Targeted protein degradation (TPD) is a promising approach in drug discovery for degrading proteins implicated in diseases. A key step in this process is the formation of a ternary complex where a heterobifunctional molecule induces proximity of an E3 ligase to a protein of interest (POI), thus facilitating ubiquitin transfer to the POI. In this work, we characterize 3 steps in the TPD process. (1) We simulate the ternary complex formation of SMARCA2 bromodomain and VHL E3 ligase by combining hydrogen-deuterium exchange mass spectrometry with weighted ensemble molecular dynamics (MD). (2) We characterize the conformational heterogeneity of the ternary complex using Hamiltonian replica exchange simulations and small-angle X-ray scattering. (3) We assess the ubiquitination of the POI in the context of the full Cullin-RING Ligase, confirming experimental ubiquitinomics results. Differences in degradation efficiency can be explained by the proximity of lysine residues on the POI relative to ubiquitin.

Experimental mapping of short-wavelength phonons in proteins.

Phonons are quasi-particles, observed as lattice vibrations in periodic materials, that often dampen in the presence of structural perturbations. Nevertheless, phonon-like collective excitations exist in highly complex systems, such as proteins, although the origin of such collective motions has remained elusive. Here we present a picture of temperature and hydration dependence of collective excitations in green fluorescent protein (GFP) obtained by inelastic neutron scattering. Our results provide evidence that such excitations can be used as a measure of flexibility/softness and are possibly associated with the protein's activity. Moreover, we show that the hydration water in GFP interferes with the phonon propagation pathway, enhancing the structural rigidity and stability of GFP.

Full structural ensembles of intrinsically disordered proteins from unbiased molecular dynamics simulations.

Molecular dynamics (MD) simulation is widely used to complement ensemble-averaged experiments of intrinsically disordered proteins (IDPs). However, MD often suffers from limitations of inaccuracy. Here, we show that enhancing the sampling using Hamiltonian replica-exchange MD (HREMD) led to unbiased and accurate ensembles, reproducing small-angle scattering and NMR chemical shift experiments, for three IDPs of varying sequence properties using two recently optimized force fields, indicating the general applicability of HREMD for IDPs. We further demonstrate that, unlike HREMD, standard MD can reproduce experimental NMR chemical shifts, but not small-angle scattering data, suggesting chemical shifts are insufficient for testing the validity of IDP ensembles. Surprisingly, we reveal that despite differences in their sequence, the inter-chain statistics of all three IDPs are similar for short contour lengths (< 10 residues). The results suggest that the major hurdle of generating an accurate unbiased ensemble for IDPs has now been largely overcome.

Mesophilic Pyrophosphatase Function at High Temperature: A Molecular Dynamics Simulation Study.

The mesophilic inorganic pyrophosphatase from Escherichia coli (EcPPase) retains function at 353 K, the physiological temperature of hyperthermophilic Thermococcus thioreducens, whereas the homolog protein (TtPPase) from this hyperthermophilic organism cannot function at room temperature. To explain this asymmetric behavior, we examined structural and dynamical properties of the two proteins using molecular dynamics simulations. The global flexibility of TtPPase is significantly higher than its mesophilic homolog at all tested temperature/pressure conditions. However, at 353 K, EcPPase reduces its solvent-exposed surface area and increases subunit compaction while maintaining flexibility in its catalytic pocket. In contrast, TtPPase lacks this adaptability and has increased rigidity and reduced protein/water interactions in its catalytic pocket at room temperature, providing a plausible explanation for its inactivity near room temperature.

Generation of the configurational ensemble of an intrinsically disordered protein from unbiased molecular dynamics simulation.

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampled molecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function.

Using Small-Angle Scattering Data and Parametric Machine Learning to Optimize Force Field Parameters for Intrinsically Disordered Proteins.

Intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) play important roles in many aspects of normal cell physiology, such as signal transduction and transcription, as well as pathological states, including Alzheimer's, Parkinson's, and Huntington's disease. Unlike their globular counterparts that are defined by a few structures and free energy minima, IDP/IDR comprise a large ensemble of rapidly interconverting structures and a corresponding free energy landscape characterized by multiple minima. This aspect has precluded the use of structural biological techniques, such as X-ray crystallography and nuclear magnetic resonance (NMR) for resolving their structures. Instead, low-resolution techniques, such as small-angle X-ray or neutron scattering (SAXS/SANS), have become a mainstay in characterizing coarse features of the ensemble of structures. These are typically complemented with NMR data if possible or computational techniques, such as atomistic molecular dynamics, to further resolve the underlying ensemble of structures. However, over the past 10-15 years, it has become evident that the classical, pairwise-additive force fields that have enjoyed a high degree of success for globular proteins have been somewhat limited in modeling IDP/IDR structures that agree with experiment. There has thus been a significant effort to rehabilitate these models to obtain better agreement with experiment, typically done by optimizing parameters in a piecewise fashion. In this work, we take a different approach by optimizing a set of force field parameters simultaneously, using machine learning to adapt force field parameters to experimental SAXS scattering profiles. We demonstrate our approach in modeling three biologically IDP ensembles based on experimental SAXS profiles and show that our optimization approach significantly improve force field parameters that generate ensembles in better agreement with experiment.

Small-Angle Neutron Scattering Reveals Energy Landscape for Rhodopsin Photoactivation.

Knowledge of the activation principles for G-protein-coupled receptors (GPCRs) is critical to development of new pharmaceuticals. Rhodopsin is the archetype for the largest GPCR family, yet the changes in protein dynamics that trigger signaling are not fully understood. Here we show that rhodopsin can be investigated by small-angle neutron scattering (SANS) in fully protiated detergent micelles under contrast matching to resolve light-induced changes in the protein structure. In SANS studies of membrane proteins, the zwitterionic detergent [(cholamidopropyl)dimethylammonio]-propanesulfonate (CHAPS) is advantageous because of the low contrast difference between the hydrophobic core and hydrophilic head groups as compared with alkyl glycoside detergents. Combining SANS results with quasielastic neutron scattering reveals how changes in volumetric protein shape are coupled (slaved) to the aqueous solvent. Upon light exposure, rhodopsin is swollen by the penetration of water into the protein core, allowing interactions with effector proteins in the visual signaling mechanism.

Collective Excitations in Protein as a Measure of Balance Between its Softness and Rigidity.

In this article, we elucidate the protein activity from the perspective of protein softness and flexibility by studying the collective phonon-like excitations in a globular protein, human serum albumin (HSA), and taking advantage of the state-of-the-art inelastic X-ray scattering (IXS) technique. Such excitations demonstrate that the protein becomes softer upon thermal denaturation due to disruption of weak noncovalent bonds. On the other hand, no significant change in the local excitations is detected in ligand- (drugs) bound HSA compared to the ligand-free HSA. Our results clearly suggest that the protein conformational flexibility and rigidity are balanced by the native protein structure for biological activity.

Quasi-elastic Neutron Scattering Reveals Ligand-Induced Protein Dynamics of a G-Protein-Coupled Receptor.

Light activation of the visual G-protein-coupled receptor (GPCR) rhodopsin leads to significant structural fluctuations of the protein embedded within the membrane yielding the activation of cognate G-protein (transducin), which initiates biological signaling. Here, we report a quasi-elastic neutron scattering study of the activation of rhodopsin as a GPCR prototype. Our results reveal a broadly distributed relaxation of hydrogen atom dynamics of rhodopsin on a picosecond-nanosecond time scale, crucial for protein function, as only observed for globular proteins previously. Interestingly, the results suggest significant differences in the intrinsic protein dynamics of the dark-state rhodopsin versus the ligand-free apoprotein, opsin. These differences can be attributed to the influence of the covalently bound retinal ligand. Furthermore, an idea of the generic free-energy landscape is used to explain the GPCR dynamics of ligand-binding and ligand-free protein conformations, which can be further applied to other GPCR systems.

Effects of pressure on the dynamics of an oligomeric protein from deep-sea hyperthermophile.

Inorganic pyrophosphatase (IPPase) from Thermococcus thioreducens is a large oligomeric protein derived from a hyperthermophilic microorganism that is found near hydrothermal vents deep under the sea, where the pressure is up to 100 MPa (1 kbar). It has attracted great interest in biophysical research because of its high activity under extreme conditions in the seabed. In this study, we use the quasielastic neutron scattering (QENS) technique to investigate the effects of pressure on the conformational flexibility and relaxation dynamics of IPPase over a wide temperature range. The ß-relaxation dynamics of proteins was studied in the time ranges from 2 to 25 ps, and from 100 ps to 2 ns, using two spectrometers. Our results indicate that, under a pressure of 100 MPa, close to that of the native environment deep under the sea, IPPase displays much faster relaxation dynamics than a mesophilic model protein, hen egg white lysozyme (HEWL), at all measured temperatures, opposite to what we observed previously under ambient pressure. This contradictory observation provides evidence that the protein energy landscape is distorted by high pressure, which is significantly different for hyperthermophilic (IPPase) and mesophilic (HEWL) proteins. We further derive from our observations a schematic denaturation phase diagram together with energy landscapes for the two very different proteins, which can be used as a general picture to understand the dynamical properties of thermophilic proteins under pressure.