Protein Conformational Fluctuations Research Articles

Molecular Insights into the Conformational and Binding Behaviors of Human Serum Albumin Induced by Surface-Active Ionic Liquids.

Extensive research has been carried out to investigate the stability and function of human serum albumin (HSA) when exposed to surface-active ionic liquids (SAILs) with different head groups (imidazolium, morpholinium, and pyridinium) and alkyl chain lengths (ranging from decyl to tetradecyl). Analysis of the protein fluorescence spectra indicates noticeable changes in the secondary structure of HSA with varying concentrations of all SAILs tested. Helicity calculations based on the Fourier transform infrared (FTIR) data show that HSA becomes more organized at the micellar concentration of SAILs, leading to an increased protein activity at this level. Small-angle neutron scattering (SANS) data confirm the formation of a bead-necklace structure between the SAILs and HSA. Atomistic molecular dynamics (MD) simulation results identify several hotspots on the protein surface for interaction with SAIL, which results in the modulation of protein conformational fluctuation and stability. Furthermore, fluorescence resonance energy transfer (FRET) experiments with the intramolecular charge transfer (ICT) probe trans-ethyl p-(dimethylamino) cinnamate (EDAC) demonstrate that higher alkyl chain lengths and SAIL concentrations result in a significantly increased energy transfer efficiency. The findings of this study provide a detailed molecular-level understanding of how the protein structure and function are affected by the presence of SAILs, with potential implications for a wide range of applications involving protein-SAIL composite systems.

Regulation Mechanism of Electron Transfer Reaction from Cytochrome C to Cytochrome C Oxidase

The respiratory chain in mitochondria is quite essential for the energy generation in cells, effectively producing “the energy currency in living organisms”, ATP. The driving force for this system is a series of electron transfers form Complex I (NADH:quinone oxidoreductase) to Complex IV (Cytochrome c Oxidase; CcO) in inner mitochondria membrane, and these electron transfer reactions promote the proton pumping from the matrix to the intermembrane space in mitochondria, and, by use of the resulting proton concentration gradient across the membrane, Complex V (ATP synthase) produces ATP. These series of electron transfer reactions terminate at CcO to reduce molecular oxygen to water molecules. The electrons required for this four-electron reduction of molecular oxygen in CcO are donated from a small soluble electron transfer hemoprotein, Cytochrome c (Cyt c), and this electron transfer is also associated with the proton pumping for the ATP generation in the respiratory chain. To examine the molecular regulation mechanism of the electron transfer reaction to terminate the electron transport chain for four-electron reduction of molecular oxygen and proton pumping for the ATP generation, we kinetically followed the electron transfer reaction and analyzed the electron transfer reaction by the Michaelis-Menten equation. By utilizing the Cyt c mutants, we identified the amino acid residues to regulate the electron transfer reactions (Biochem. J., 2020, 477,1565). The kinetic analysis of the electron transfer reaction also revealed that the apparent rate for the electron transfer reaction determined by the steady-state kinetics is much slower than the intramolecular electron transfer rate in the Cyt c – CcO complex. On the other hand, we successfully determined the interaction site between Cyt c and CcO by using NMR (Proc. Natl. Acad. Sci., 2011, 108, 12271) and protein-protein docking simulation (J. Biol. Chem., 2016, 291, 15320). The Cyt c – CcO complex structure estimated from the NMR and protein-protein docking simulation clearly showed that the two redox centers, heme iron in Cyt c and CuA in subunit II of CcO, were located more than 20 Å apart, and the complexes exhibiting shorter redox center distances were energetically unstable. Combined with the results from kinetic and structural studies, we proposed a “conformationally gated electron transfer mechanism”, where the thermal fluctuation of the protein structure controls the electron transfer rate. However, direct experimental evidence that protein conformational fluctuations regulate electron transfer reactions has not been obtained. In this study, we focused on the viscosity of the protein solution as a factor to regulate the conformational fluctuation of the proteins, based on the assumption that high viscosity of the solution suppresses the conformational fluctuation. The electron transfer reaction from Cyt c to CcO was examined in the high viscosity solution and NMR relaxation measurements under high viscosity were also performed to identify the local structural fluctuation responsible for the regulation of the electron transfer rate. Together with the transient complex structures between Cyt c and CcO estimated by the MD simulation, the regulation mechanism of the electron transfer from Cyt c to CcO based on the “conformational gating” will be discussed.

NMR Studies of Aromatic Ring Flips to Probe Conformational Fluctuations in Proteins.

Aromatic residues form a significant part of the protein core, where they make tight interactions with multiple surrounding side chains. Despite the dense packing of internal side chains, the aromatic rings of phenylalanine and tyrosine residues undergo 180° rotations, or flips, which are mediated by transient and large-scale "breathing" motions that generate sufficient void volume around the aromatic ring. Forty years after the seminal work by Wagner and Wüthrich, NMR studies of aromatic ring flips are now undergoing a renaissance as a powerful means of probing fundamental dynamic properties of proteins. Recent developments of improved NMR methods and isotope labeling schemes have enabled a number of advances in addressing the mechanisms and energetics of aromatic ring flips. The nature of the transition states associated with ring flips can be described by thermodynamic activation parameters, including the activation enthalpy, activation entropy, activation volume, and also the isothermal volume compressibility of activation. Consequently, it is of great interest to study how ring flip rate constants and activation parameters might vary with protein structure and external conditions like temperature and pressure. The field is beginning to gather such data for aromatic residues in a variety of environments, ranging from surface exposed to buried. In the future, the combination of solution and solid-state NMR spectroscopy together with molecular dynamics simulations and other computational approaches is likely to provide detailed information about the coupled dynamics of aromatic rings and neighboring residues. In this Perspective, we highlight recent developments and provide an outlook toward the future.

Ultra-sensitive lock-in amplifier coupled oscillatory magnetic tweezers for piconewton force manipulation applications

We have developed lock-in amplifier coupled oscillatory magnetic tweezers, aiming to synchronize the oscillatory magnetic force application and single-molecule response detection at the applied oscillatory frequency by incorporating an optical lock-in amplifier detection. The designed home-built lock-in amplifier detection instrument enables an exact reference signal input and hence reveals an accurate and sensitive magnetic response synchronization. We have further demonstrated the approach with the rhodamine 6G stained super-paramagnetic beads by monitoring the below-the-noise-background weak fluorescence signal changes due to the magnetic response of the super-paramagnetic beads under the oscillatory force manipulation. The integration of the lock-in amplifier and the oscillating magnetic tweezers can significantly expand the application of the magnetic tweezers for signal detection below the noise background, such as adapting to the important applications in the detailed exploration of mechanical properties of biomolecules and studies of the protein conformational fluctuation dynamics.

Transition-State Compressibility and Activation Volume of Transient Protein Conformational Fluctuations.

Proteins are dynamic entities that intermittently depart from their ground-state structures and undergo conformational transitions as a critical part of their functions. Central to understanding such transitions are the structural rearrangements along the connecting pathway, where the transition state plays a special role. Using NMR relaxation at variable temperature and pressure to measure aromatic ring flips inside a protein core, we obtain information on the structure and thermodynamics of the transition state. We show that the isothermal compressibility coefficient of the transition state is similar to that of short-chain hydrocarbon liquids, implying extensive local unfolding of the protein. Our results further indicate that the required local volume expansions of the protein can occur not only with a net positive activation volume of the protein, as expected from previous studies, but also with zero activation volume by compaction of remote void volume, when averaged over the ensemble of states.

Universal Relation between Instantaneous Diffusivity and Radius of Gyration of Proteins in Aqueous Solution.

Protein conformational fluctuations are highly complex and exhibit long-term correlations. Here, molecular dynamics simulations of small proteins demonstrate that these conformational fluctuations directly affect the protein's instantaneous diffusivity D_{I}. We find that the radius of gyration R_{g} of the proteins exhibits 1/f fluctuations that are synchronous with the fluctuations of D_{I}. Our analysis demonstrates the validity of the local Stokes-Einstein-type relation D_{I}∝1/(R_{g}+R_{0}), where R_{0}∼0.3 nm is assumed to be a hydration layer around the protein. From the analysis of different protein types with both strong and weak conformational fluctuations, the validity of the Stokes-Einstein-type relation appears to be a general property.

Single-Protein Spectroscopy Reveals Multiple Protein Conformational Fluctuations Regulating Photosynthetic Light Harvesting

Single-Protein Spectroscopy Reveals Multiple Protein Conformational Fluctuations Regulating Photosynthetic Light Harvesting

Nonaffine Displacements Encode Collective Conformational Fluctuations in Proteins.

Identifying subtle conformational fluctuations underlying the dynamics of biomacromolecules is crucial for resolving their free energy landscape. We show that a collective variable, originally proposed for crystalline solids, is able to filter out essential macromolecular motions more efficiently than other approaches. While homogeneous or "affine" deformations of the biopolymer are trivial, biopolymer conformations are complicated by the occurrence of inhomogeneous or "nonaffine" displacements of atoms relative to their positions in the native structure. We show that these displacements encode functionally relevant conformations of macromolecules, and in combination with a formalism based upon time-structured independent component analysis, they quantitatively resolve the free energy landscape of a number of macromolecules of hierarchical complexity. The kinetics of conformational transitions among the basins can now be mapped within the framework of a Markov state model. The nonaffine modes, obtained by projecting out homogeneous fluctuations from the local displacements, are found to be responsible for local structural changes required for transitioning between pairs of macrostates.

Paramagnetic relaxation enhancement-assisted structural characterization of a partially disordered conformation of ubiquitin.

Nuclear magnetic resonance (NMR) is a powerful tool to study three-dimensional structures as well as protein conformational fluctuations in solution, but it is compromised by increases in peak widths and missing signals. We previously reported that ubiquitin has two folded conformations, N1 and N2 and plus another folded conformation, I, in which some amide group signals of residues 33-41 almost disappeared above 3 kbar at pH 4.5 and 273 K. Thus, well-converged structural models could not be obtained for this region owing to the absence of distance restraints. Here, we reexamine the problem using the ubiquitin Q41N variant as a model for this locally disordered conformation, I. We demonstrate that the variant shows pressure-induced loss of backbone amide group signals at residues 28, 33, 36, and 39-41 like the wild-type, with a similar but smaller effect on CαH and CβH signals. In order to characterize this I structure, we measured paramagnetic relaxation enhancement (PRE) under high pressure to obtain distance restraints, and calculated the structure assisted by Bayesian inference. We conclude that the more disordered I conformation observed at pH 4.0, 278 K, and 2.5 kbar largely retained the N2 conformation, although the amide groups at residues 33-41 have more heterogeneous conformations and more contact with water, which differ from the N1 and N2 states. The PRE-assisted strategy has the potential to improve structural characterization of proteins that lack NMR signals, especially for relatively more open and hydrated protein conformations.

Aggrescan3D standalone package for structure-based prediction of protein aggregation properties.

Aggrescan3D (A3D) standalone is a multiplatform Python package for structure-based prediction of protein aggregation properties and rational design of protein solubility. A3D allows the re-design of protein solubility by combining structural aggregation propensity and stability predictions, as demonstrated by a recent experimental study. It also enables predicting the impact of protein conformational fluctuations on the aggregation properties. The standalone A3D version is an upgrade of the original web server implementation-it introduces a number of customizable options, automated analysis of multiple mutations and offers a flexible computational framework for merging it with other computational tools. A3D standalone is distributed under the MIT license, which is free for academic and non-profit users. It is implemented in Python. The A3D standalone source code, wiki with documentation and examples of use, and installation instructions for Linux, macOS and Windows are available in the A3D standalone repository at https://bitbucket.org/lcbio/aggrescan3d.

Improving Coarse-Grained Protein Force Fields with Small-Angle X-ray Scattering Data.

Small-angle X-ray scattering (SAXS) experiments provide valuable structural data for biomolecules in solution. We develop a highly efficient maximum entropy approach to fit SAXS data by introducing minimal biases to a coarse-grained protein force field, the associative memory, water mediated, structure, and energy model (AWSEM). We demonstrate that the resulting force field, AWSEM-SAXS, succeeds in reproducing scattering profiles and models protein structures with shapes that are in much better agreement with experimental results. Quantitative metrics further reveal a modest, but consistent, improvement in the accuracy of modeled structures when SAXS data are incorporated into the force field. Additionally, when applied to a multiconformational protein, we find that AWSEM-SAXS is able to recover the population of different protein conformations from SAXS data alone. We, therefore, conclude that the maximum entropy approach is effective in fine-tuning the force field to better characterize both protein structure and conformational fluctuation.

Determination of Structural Ensembles of Proteins: Restraining vs Reweighting.

The conformational fluctuations of proteins can be described using structural ensembles. To address the challenge of determining these ensembles accurately, a wide range of strategies have recently been proposed to combine molecular dynamics simulations with experimental data. Quite generally, there are two ways of implementing this type of approach, either by applying structural restraints during a simulation, or by reweighting a posteriori the conformations from an a priori ensemble. It is not yet clear, however, whether these two approaches can offer ensembles of equivalent quality. The advantages of the reweighting method are that it can involve any type of starting simulation and that it enables the integration of experimental data after the simulations are run. A disadvantage, however, is that this procedure may be inaccurate when the a priori ensemble is of poor quality. Here, our goal is to systematically compare the restraining and reweighting approaches and to explore the conditions required for the reweighted ensembles to be accurate. Our results indicate that the reweighting approach is computationally efficient and can perform as well as the restraining approach when the a priori sampling is already relatively accurate. More generally, to enable an effective use of the reweighting approach by avoiding the pitfalls of poor sampling, we suggest metrics for the quality control of the reweighted ensembles.

Water-Protein Interactions Coupled with Protein Conformational Transition

Conformational fluctuations of proteins are crucially important for their functions. However, changes in the location and dynamics of hydrated water in many proteins accompanied by the conformational transition have not been fully understood. Here, we used phase-modulated clean chemical exchange NMR approach to investigate pressure-induced changes in water-to-amide proton exchange occurring at sub-second time scale. With the transition of ubiquitin from its native conformation (N1) to an alternative conformation (N2) at 250 MPa, proton exchange rates of residues 32–35, 40–41, and 71, which are located at the C-terminal side of the protein, were significantly increased. These observations can be explained by the destabilization of the hydrogen bonds in the backbone and partial exposure of those amide groups to solvent in N2. We conclude that phase-modulated clean chemical exchange NMR approach coupled with pressure perturbation will be a useful tool for investigations of more open and hydrated protein structures.

AGGRESCAN3D: Toward the Prediction of the Aggregation Propensities of Protein Structures.

Protein aggregation is responsible for the onset and spread of many human diseases, ranging from neurodegenerative disorders to cancer and diabetes. Moreover, it is one of the major bottlenecks for the production of protein-based therapeutics such as antibodies or enzymes. AGGRESCAN3D (A3D) is a web server aimed to identify and evaluate structural aggregation prone regions, overcoming the limitations of sequence-based algorithms in the prediction of the aggregation propensity of globular proteins. A3D allows the redesign of protein solubility by predicting in silico the impact of mutations and protein conformational fluctuations on the aggregation of native polypeptides.

Coupled Valence-Bond State Molecular Dynamics Description of an Enzyme-Catalyzed Reaction in a Non-Aqueous Organic Solvent.

Enzymes are widely used in nonaqueous solvents to catalyze non-natural reactions. While experimental measurements showed that the solvent nature has a strong effect on the reaction kinetics, the molecular details of the catalytic mechanism in nonaqueous solvents have remained largely elusive. Here we study the transesterification reaction catalyzed by the paradigm subtilisin Carlsberg serine protease in an organic apolar solvent. The rate-limiting acylation step involves a proton transfer between active-site residues and the nucleophilic attack of the substrate to form a tetrahedral intermediate. We design the first coupled valence-bond state model that simultaneously describes both reactions in the enzymatic active site. We develop a new systematic procedure to parametrize this model on high-level ab initio QM/MM free energy calculations that account for the molecular details of the active site and for both substrate and protein conformational fluctuations. Our calculations show that the reaction energy barrier changes dramatically with the solvent and protein conformational fluctuations. We find that the mechanism of the tetrahedral intermediate formation during the acylation step is similar to that determined under aqueous conditions, and that the proton transfer and nucleophilic attack reactions occur concertedly. We identify the reaction coordinate to be mostly due to the rearrangement of some residual water molecules close to the active site.

Molecular impact of covalent modifications on nonribosomal peptide synthetase carrier protein communication

Nonribosomal peptide synthesis involves the interplay between covalent protein modifications, conformational fluctuations, catalysis, and transient protein-protein interactions. Delineating the mechanisms involved in orchestrating these various processes will deepen our understanding of domain-domain communication in nonribosomal peptide synthetases (NRPSs) and lay the groundwork for the rational reengineering of NRPSs by swapping domains handling different substrates to generate novel natural products. Although many structural and biochemical studies of NRPSs exist, few studies have focused on the energetics and dynamics governing the interactions in these systems. Here, we present detailed binding studies of an adenylation domain and its partner carrier protein in apo-, holo-, and substrate-loaded forms. Results from fluorescence anisotropy, isothermal titration calorimetry, and NMR titrations indicated that covalent modifications to a carrier protein modulate domain communication, suggesting that chemical modifications to carrier proteins during NRPS synthesis may impart directionality to sequential NRPS domain interactions. Comparison of the structure and dynamics of an apo-aryl carrier protein with those of its modified forms revealed structural fluctuations induced by post-translational modifications and mediated by modulations of protein dynamics. The results provide a comprehensive molecular description of a carrier protein throughout its life cycle and demonstrate how a network of dynamic residues can propagate the molecular impact of chemical modifications throughout a protein and influence its affinity toward partner domains.

Decomposition of total solvation energy into core, side-chains and water contributions: Role of cross correlations and protein conformational fluctuations in dynamics of hydration layer

Dynamical coupling between water and amino acid side-chain residues in solvation dynamics is investigated by selecting residues often used as natural probes, namely tryptophan, tyrosine and histidine, located at different positions on protein surface. Such differently placed residues are found to exhibit different timescales of relaxation. The total solvation response measured by the probe is decomposed in terms of its interactions with (i) protein core, (ii) side-chain and (iii) water. Significant anti cross-correlation among these contributions are observed. When the motion of the protein side-chains is quenched, solvation either becomes faster or slower depending on the location of the probe.

High-pressure NMR techniques for the study of protein dynamics, folding and aggregation

High-pressure is a well-known perturbation method used to destabilize globular proteins and dissociate protein complexes or aggregates. The heterogeneity of the response to pressure offers a unique opportunity to dissect the thermodynamic contributions to protein stability. In addition, pressure perturbation is generally reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility of monitoring at an atomic resolution the structural transitions occurring upon unfolding and determining the kinetic properties of the process. The recent development of commercially available high-pressure sample cells greatly increased the potential applications for high-pressure NMR experiments that can now be routinely performed. This review summarizes the recent applications and future directions of high-pressure NMR techniques for the characterization of protein conformational fluctuations, protein folding and the stability of protein complexes and aggregates.

How Conformational Dynamics Influences the Protein Search for Targets on DNA

Protein search and association to specific target sequences on DNA is essential for all fundamental biological processes. The detailed qualitative and quantitative mechanism of fast target finding is still unknown in spite of numerous experimental and theoretical efforts. Particularly, the role of protein conformational fluctuations in the search dynamics remains uncovered. We extend developed earlier theoretical method to analyze how the conformational dynamic affects the process of finding the specific targets on DNA. Our approach is based on discrete-state stochastic model that takes into account main physical-chemical processes. This allows us to evaluate explicitly the protein search for the targets on DNA at different conditions. Our calculations show that conformational fluctuations might strongly affect the protein search dynamics. We discuss the contribution of the shift in the conformational equilibrium in the target search kinetics. We utilized extensive Monte Carlo computer simulations to validate our theoretical predictions.

静水圧及びガス圧NMR法を用いたタンパク質の動的性質の研究

Pressure and temperature are fundamental physical parameters to control molecular properties and equilibria of systems. However, pressure approaches to life science are still very limited compared with temperature, due probably to their technical difficulties. Here, nuclear magnetic resonance (NMR) approaches combined with hydrostatic pressure and gas pressure are reviewed. Using NMR spectroscopy with variable pressure is a useful method for characterizing conformational fluctuation of proteins. Gas-pressure NMR spectroscopy using molecular oxygen provides a general and highly sensitive method for detecting internal hydrophobic cavities, which might be sources of protein's conformational fluctuation.