scholarly journals Conformational stability and dynamics in crystals recapitulate protein behaviour in solution

2020 ◽  
Author(s):  
BM Sala ◽  
T Le Marchand ◽  
G Pintacuda ◽  
C Camilloni ◽  
A Natalello ◽  
...  

AbstractA growing body of evidences has established that in many cases proteins may preserve most of their function and flexibility in a crystalline environment, and several techniques are today capable to detect transiently-populated states of macromolecules in tightly packed lattices. Intriguingly, in the case of amyloidogenic precursors, the presence of these conformations (hidden to conventional crystallographic studies) can be correlated to the pathological fate of the native fold.It remains unclear, however, to which extent these minor conformations reflect the protein behaviour that is more commonly studied in solution. Here, we address this question by investigating some biophysical properties of a prototypical amyloidogenic system, β2-microglobulin (β2m) in solution and in microcrystalline state.By combining NMR chemical shifts with Molecular Dynamics (MD) simulations, we confirmed that conformational dynamics of β2m native state in the crystal lattice is in keeping with what observed in solution.A comparative study of protein stability in solution and in crystallo is then carried out, monitoring the change in protein secondary structure at increasing temperature by Fourier transform infrared (FTIR) spectroscopy. The increased structural order of the crystalline state contributes to provide better resolved spectral components compared to those collected in solution and crucially, the crystalline samples display thermal stabilities in good agreement with the trend observed in solution.Overall, this work shows that protein stability and occurrence of pathological hidden states in crystals parallel their solution counterpart, confirming the interest of crystals as a platform for the biophysical characterisation of processes such as unfolding and aggregation.

2020 ◽  
Vol 14 (3) ◽  
pp. 216-226
Author(s):  
Priyanka Borah ◽  
Venkata S.K. Mattaparthi

Background: Aggregation of misfolded proteins under stress conditions in the cell might lead to several neurodegenerative disorders. Amyloid-beta (Aβ1-42) peptide, the causative agent of Alzheimer’s disease, has the propensity to fold into β-sheets under stress, forming aggregated amyloid plaques. This is influenced by factors such as pH, temperature, metal ions, mutation of residues, and ionic strength of the solution. There are several studies that have highlighted the importance of ionic strength in affecting the folding and aggregation propensity of Aβ1-42 peptide. Objective: To understand the effect of ionic strength of the solution on the aggregation propensity of Aβ1-42 peptide, using computational approaches. Materials and Methods: In this study, Molecular Dynamics (MD) simulations were performed on Aβ1-42 peptide monomer placed in (i) 0 M, (ii) 0.15 M, and (iii) 0.30 M concentration of NaCl solution. To prepare the input files for the MD simulations, we have used the Amberff99SB force field. The conformational dynamics of Aβ1-42 peptide monomer in different ionic strengths of the solutions were illustrated from the analysis of the corresponding MD trajectory using the CPPtraj tool. Results: From the MD trajectory analysis, we observe that with an increase in the ionic strength of the solution, Aβ1-42 peptide monomer shows a lesser tendency to undergo aggregation. From RMSD and SASA analysis, we noticed that Aβ1-42 peptide monomer undergoes a rapid change in conformation with an increase in the ionic strength of the solution. In addition, from the radius of gyration (Rg) analysis, we observed Aβ1-42 peptide monomer to be more compact at moderate ionic strength of the solution. Aβ1-42 peptide was also found to hold its helical secondary structure at moderate and higher ionic strengths of the solution. The diffusion coefficient of Aβ1-42 peptide monomer was also found to vary with the ionic strength of the solution. We observed a relatively higher diffusion coefficient value for Aβ1-42 peptide at moderate ionic strength of the solution. Conclusion: Our findings from this computational study highlight the marked effect of ionic strength of the solution on the conformational dynamics and aggregation propensity of Aβ1-42 peptide monomer.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Raghavendar Reddy Sanganna Gari ◽  
Joel José Montalvo‐Acosta ◽  
George R. Heath ◽  
Yining Jiang ◽  
Xiaolong Gao ◽  
...  

AbstractConformational changes in ion channels lead to gating of an ion-conductive pore. Ion flux has been measured with high temporal resolution by single-channel electrophysiology for decades. However, correlation between functional and conformational dynamics remained difficult, lacking experimental techniques to monitor sub-millisecond conformational changes. Here, we use the outer membrane protein G (OmpG) as a model system where loop-6 opens and closes the β-barrel pore like a lid in a pH-dependent manner. Functionally, single-channel electrophysiology shows that while closed states are favored at acidic pH and open states are favored at physiological pH, both states coexist and rapidly interchange in all conditions. Using HS-AFM height spectroscopy (HS-AFM-HS), we monitor sub-millisecond loop-6 conformational dynamics, and compare them to the functional dynamics from single-channel recordings, while MD simulations provide atomistic details and energy landscapes of the pH-dependent loop-6 fluctuations. HS-AFM-HS offers new opportunities to analyze conformational dynamics at timescales of domain and loop fluctuations.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yogeeshwar Ajjugal ◽  
Kripi Tomar ◽  
D. Krishna Rao ◽  
Thenmalarchelvi Rathinavelan

AbstractBase pair mismatches in DNA can erroneously be incorporated during replication, recombination, etc. Here, the influence of A…A mismatch in the context of 5′CAA·5′TAG sequence is explored using molecular dynamics (MD) simulation, umbrella sampling MD, circular dichroism (CD), microscale thermophoresis (MST) and NMR techniques. MD simulations reveal that the A…A mismatch experiences several transient events such as base flipping, base extrusion, etc. facilitating B–Z junction formation. A…A mismatch may assume such conformational transitions to circumvent the effect of nonisostericity with the flanking canonical base pairs so as to get accommodated in the DNA. CD and 1D proton NMR experiments further reveal that the extent of B–Z junction increases when the number of A…A mismatch in d(CAA)·d(T(A/T)G) increases (1–5). CD titration studies of d(CAA)·d(TAG)n=5 with the hZαADAR1 show the passive binding between the two, wherein, the binding of protein commences with B–Z junction recognition. Umbrella sampling simulation indicates that the mismatch samples anti…+ syn/+ syn…anti, anti…anti & + syn…+ syn glycosyl conformations. The concomitant spontaneous transitions are: a variety of hydrogen bonding patterns, stacking and minor or major groove extrahelical movements (with and without the engagement of hydrogen bonds) involving the mismatch adenines. These transitions frequently happen in anti…anti conformational region compared with the other three regions as revealed from the lifetime of these states. Further, 2D-NOESY experiments indicate that the number of cross-peaks diminishes with the increasing number of A…A mismatches implicating its dynamic nature. The spontaneous extrahelical movement seen in A…A mismatch may be a key pre-trapping event in the mismatch repair due to the accessibility of the base(s) to the sophisticated mismatch repair machinery.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Mei Dang ◽  
Yifan Li ◽  
Jianxing Song

AbstractTDP-43 and hnRNPA1 contain tandemly-tethered RNA-recognition-motif (RRM) domains, which not only functionally bind an array of nucleic acids, but also participate in aggregation/fibrillation, a pathological hallmark of various human diseases including amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD), alzheimer's disease (AD) and Multisystem proteinopathy (MSP). Here, by DSF, NMR and MD simulations we systematically characterized stability, ATP-binding and conformational dynamics of TDP-43 and hnRNPA1 RRM domains in both tethered and isolated forms. The results reveal three key findings: (1) upon tethering TDP-43 RRM domains become dramatically coupled and destabilized with Tm reduced to only 49 °C. (2) ATP specifically binds TDP-43 and hnRNPA1 RRM domains, in which ATP occupies the similar pockets within the conserved nucleic-acid-binding surfaces, with the affinity slightly higher to the tethered than isolated forms. (3) MD simulations indicate that the tethered RRM domains of TDP-43 and hnRNPA1 have higher conformational dynamics than the isolated forms. Two RRM domains become coupled as shown by NMR characterization and analysis of inter-domain correlation motions. The study explains the long-standing puzzle that the tethered TDP-43 RRM1–RRM2 is particularly prone to aggregation/fibrillation, and underscores the general role of ATP in inhibiting aggregation/fibrillation of RRM-containing proteins. The results also rationalize the observation that the risk of aggregation-causing diseases increases with aging.


RSC Advances ◽  
2021 ◽  
Vol 11 (3) ◽  
pp. 1320-1331
Author(s):  
Shahzaib Ahamad ◽  
Hema Kanipakam ◽  
Vijay Kumar ◽  
Dinesh Gupta

MD simulations of TTBK2 mutants to study its impact on stability of the protein.


1998 ◽  
Vol 543 ◽  
Author(s):  
T. Çağin ◽  
Y. Zhou ◽  
E. S. Yamaguchi ◽  
R. Frazier ◽  
A. Ho ◽  
...  

AbstractTo understand antiwear phenomena in motor engines at the atomic level and provide evidence inselecting future ashless wear inhibitors, we studied the thermal stability of the self-assembled monolayer(SAM) model for dithiophosphate (DTP) and dithiocarbamate (DTC) molecules on the iron oxidesurface using molecular dynamics. The interactions for DTP, DTC and Fe2O3 are evaluated based on aforce field derived from fitting to ab initio quantum chemical calculations of dimethyl DTP (and DTC)and Fe(OH)2(H2O)2-DTP (DTC) clusters. MD simulations at constant-NPT are conducted to assesrelative thermal stabilities of the DTP and DTC with different pendant groups (n-propyl, i-propyl, npentyl.and i-pentyl). To investigate frictional process, we employ a steady state MD method, in whichone of the Fe2O3 slabs maintained at a constant linear velocity. We obtain the time averaged normaland frictional forces from the interatomic forces. Then, we calculated the friction coefficient at theinterface between SAMs of DTP and the confined lubricant, hexadecane, to assess the shear stability ofDTPs with different pendant groups.


1969 ◽  
Vol 47 (4) ◽  
pp. 480-484 ◽  
Author(s):  
Ian C. P. Smith ◽  
Tetsuo Yamane ◽  
R. G. Shulman

Proton magnetic resonance spectra at 220 MHz of alanine transfer RNA do not permit assignments of individual peaks due to each of the common bases; only a peak attributable to protons at position eight in adenine can be assigned with certainty. Measurements of the relative areas of proton magnetic resonance peaks due to the base and ribose-1′ protons indicate that the ribose moieties of tRNA are not involved in bonds stronger than those experienced by the bases. Proton magnetic resonance peaks attributable to the methyl and dihydro protons of the rare bases can be distinguished in the 220 MHz spectra; the variation of their line widths and chemical shifts with increasing temperature indicates that the rare bases are located in regions of the alanine transfer RNA molecule which are more highly organized than indicated by an open cloverleaf model.


2021 ◽  
Author(s):  
Zachary Smith ◽  
Pratyush Tiwary

Molecular dynamics (MD) simulations provide a wealth of high-dimensional data at all-atom and femtosecond resolution but deciphering mechanistic information from this data is an ongoing challenge in physical chemistry and biophysics. Theoretically speaking, joint probabilities of the equilibrium distribution contain all thermodynamic information, but they prove increasingly difficult to compute and interpret as the dimensionality increases. Here, inspired by tools in probabilistic graphical modeling, we develop a factor graph trained through belief propagation that helps factorize the joint probability into an approximate tractable form that can be easily visualized and used. We validate the study through the analysis of the conformational dynamics of two small peptides with 5 and 9 residues. Our validations include testing the conditional dependency predictions through an intervention scheme inspired by Judea Pearl. Secondly we directly use the belief propagation based approximate probability distribution as a high-dimensional static bias for enhanced sampling, where we achieve spontaneous back-and-forth motion between metastable states that is up to 350 times faster than unbiased MD. We believe this work opens up useful ways to thinking about and dealing with high-dimensional molecular simulations.


Sign in / Sign up

Export Citation Format

Share Document