Cosolvent Effects on the Growth of Protein Aggregates Formed by a Single Domain Globular Protein and an Intrinsically Disordered Protein

Understanding the relationship between preferential interactions of peptides in water-acetonitrile mixtures with protein-solvent contact surface area

The influence of polar, water-miscible organic solvents (POS) on protein structure, stability, and functional activity is a subject of significant interest and complexity. This study examines the effects of acetonitrile (ACN), a semipolar, aprotic solvent, on the solvation properties of blocked Ace-Gly-X-Gly-Nme tripeptides (where Ace and Nme stands for acetyl and N-methyl amide groups respectively and X is any amino acid) through extensive molecular dynamics simulations. Individual simulations were conducted for each peptide, encompassing five different ACN concentrations within the range of χACN = 0.1-0.9. The preferential solvation parameter (Γ) calculated using the Kirkwood-Buff integral method was used for the assessment of peptide interactions with water/ACN. Additionally, weighted Voronoi tessellation was applied to obtain a three-way data set containing four time-averaged contact surface area types between peptide atoms and water/ACN atoms. A mathematical technique known as N-way Partial Least Squares (NPLS) was utilized to anticipate the preferential interactions between peptides and water/ACN from the contact surface areas. Furthermore, the temperature dependency of peptide-solvent interactions was investigated using a subset of 10 amino acids representing a range of hydrophobicities. MD simulations were conducted at five temperatures, spanning from 283 to 343 K, with subsequent analysis of data focusing on both preferential solvation and peptide-solvent contact surface areas. The results demonstrate the efficacy of utilizing contact surface areas between the peptide and solvent constituents for successfully predicting preferential interactions in water/ACN mixtures across various ACN concentrations and temperatures.

Cooperative Membrane Binding of HIV-1 Matrix Proteins

The HIV-1 assembly process begins with a newly synthesized Gag polyprotein being targeted to the inner leaflet of the plasma membrane of the infected cells to form immature viral particles. Gag-membrane interactions are mediated through the myristoylated (Myr) N-terminal matrix (MA) domain of Gag, which eventually multimerize on the membrane to form trimers and higher order oligomers. The study of the structure and dynamics of peripheral membrane proteins like MA has been challenging for both experimental and computational studies due to the complex transient dynamics of protein-membrane interactions. Although the roles of anionic phospholipids (PIP2, PS) and the Myr group in the membrane targeting and stable membrane binding of MA are now well-established, the cooperative interactions between the MA monomers and MA-membrane remain elusive in the context of viral assembly and release. Our present study focuses on the membrane binding dynamics of a higher order oligomeric structure of MA protein (a dimer of trimers), which has not been explored before. Employing time-lagged independent component analysis (tICA) to our microsecond-long trajectories, we investigate conformational changes of the matrix protein induced by membrane binding. Interestingly, the Myr switch of an MA monomer correlates with the conformational switch of adjacent monomers in the same trimer. Together, our findings suggest complex protein dynamics during the formation of the immature HIV-1 lattice; while MA trimerization facilitates Myr insertion, MA trimer-trimer interactions in the immature lattice can hinder the same.

Cooperative Membrane Binding of HIV-1 Matrix Proteins

The HIV-1 assembly process begins with a newly synthesized Gag polyprotein being targeted to the inner leaflet of the plasma membrane of the infected cells to form immature viral particles. Gag-membrane interactions are mediated through the myristoylated(Myr) N-terminal matrix (MA) domain of Gag which eventually multimerize on the membrane to form trimers and higher-order oligomers. The study of the structure and dynamics of peripheral membrane proteins like MA has been challenging for both experimental and computational studies due to the complex dynamics of protein-membrane interactions. Although the roles of anionic phospholipids (PIP2, PS) and the Myr group in the membrane targeting and stable membrane binding of MA are now well-established, the cooperative interactions between MA monomers and MA-membrane still remain elusive. Our present study focuses on the membrane binding dynamics of a higher-order oligomeric structure of MA protein (a dimer of trimers), which has not been explored before. Employing time-lagged independent component analysis (tICA) to our microsecond-long trajectories, we investigate conformational changes of the matrix protein induced by membrane binding. Interestingly, the Myr switch of a MA monomer correlates with the conformational switch of adjacent monomers in the same trimer. Together, our findings suggest that MA trimerization facilitates Myr insertion, but MA trimer-trimer interactions in the lattice of immature HIV-1 particles can hinder the same. Additionally, local lipid density patterns of different lipid species provide a signature of the initial stage of lipid-domain formation upon membrane binding of the protein complex.

TOC

pH Induced Switch in the Conformational Ensemble of Intrinsically Disordered Protein Prothymosin-α and Its Implications for Amyloid Fibril Formation

Aggregation of intrinsically disordered proteins (IDPs) can lead to neurodegenerative diseases. Although there is experimental evidence that acidic pH promotes IDP monomer compaction leading to aggregation, the general mechanism is unclear. We studied the pH effect on the conformational ensemble of prothymosin-α (proTα), which is involved in multiple essential functions, and probed its role in aggregation using computer simulations. We show that compaction in the proTα dimension at low pH is due to the protein's collapse in the intermediate region (E41-D80) rich in glutamic acid residues, enhancing its β-sheet content. We observed by performing dimer simulations that the conformations with high β-sheet content could act as aggregation-prone (N*) states and nucleate the aggregation process. The simulations initiated using N* states form dimers within a microsecond time scale, whereas the non-N* states do not form dimers within this time scale. This study contributes to understanding the general principles of pH-induced IDP aggregation.

Trimethylamine N-oxide alters structure-function integrity of β-casein: Structural disorder co-regulates the aggregation propensity and chaperone activity

Intrinsically disordered proteins (IDPs), involved in the regulation and function of various cellular processes like transcription, translation, cell cycle etc., exist as ensembles of rapidly interconverting structures with functional plasticity. Among numerous cellular regulatory mechanisms involved in structural and functional regulation of IDPs, osmolytes are emerging as promising regulatory agents due to their ability to affect the structure-function integrity of IDPs. The present study investigated the effect of methylamine osmolytes on β-casein, an IDP essential for maintaining the overall stability of casein complex in milk. It was observed that trimethylamine N-oxide induces a compact structural state in β-casein with slightly decreased chaperone activity and insignificant aggregation propensity. However, the other two osmolytes from this group, i.e., sarcosine and betaine, had no significant effect on the overall structure and chaperone activity of the IDP. The present study hints towards the possible evolutionary selection of higher structural disorder in β-casein, compared to α-casein, for stability of the casein complex and prevention of amyloidosis in the mammary gland.

Double Domain Swapping in Human γC and γD Crystallin Drives Early Stages of Aggregation

Human γD (HγD) and γC (HγC) are two-domain crystallin (Crys) proteins expressed in the nucleus of the eye lens. Structural perturbations in the protein often trigger aggregation, which eventually leads to cataract. To decipher the underlying molecular mechanism, it is important to characterize the partially unfolded conformations, which are aggregation-prone. Using a coarse grained protein model and molecular dynamics simulations, we studied the role of on-pathway folding intermediates in the early stages of aggregation. The multidimensional free energy surface revealed at least three different folding pathways with the population of partially structured intermediates. The two dominant pathways confirm sequential folding of the N-terminal [Ntd] and the C-terminal domains [Ctd], while the third, least favored, pathway involves intermediates where both the domains are partially folded. A native-like intermediate (I*), featuring the folded domains and disrupted interdomain contacts, gets populated in all three pathways. I* forms domain swapped dimers by swapping the entire Ntds and Ctds with other monomers. Population of such oligomers can explain the increased resistance to unfolding resulting in hysteresis observed in the folding experiments of HγD Crys. An ensemble of double domain swapped dimers are also formed during refolding, where intermediates consisting of partially folded Ntds and Ctds swap secondary structures with other monomers. The double domain swapping model presented in our study provides structural insights into the early events of aggregation in Crys proteins and identifies the key secondary structural swapping elements, where introducing mutations will aid in regulating the overall aggregation propensity.

MEPSAnd: minimum energy path surface analysis over n-dimensional surfaces

SUMMARY

n-dimensional energy surfaces are becoming computationally accessible, yet interpreting their information is not straightforward. We present minimum energy path surface analysis over n-dimensional surfaces (MEPSAnd), an open source GUI-based program that natively calculates minimum energy paths across energy surfaces of any number of dimensions. Among other features, MEPSAnd can compute the path through lowest barriers and automatically provide a set of alternative paths. MEPSAnd offers distinct plotting solutions as well as direct python scripting.

AVAILABILITY AND IMPLEMENTATION

MEPSAnd is freely available (under GPLv3 license) at: http://bioweb.cbm.uam.es/software/MEPSAnd/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Recent Advances in Computational Protocols Addressing Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDP) are abundant in the human genome and have recently emerged as major therapeutic targets for various diseases. Unlike traditional proteins that adopt a definitive structure, IDPs in free solution are disordered and exist as an ensemble of conformations. This enables the IDPs to signal through multiple signaling pathways and serve as scaffolds for multi-protein complexes. The challenge in studying IDPs experimentally stems from their disordered nature. Nuclear magnetic resonance (NMR), circular dichroism, small angle X-ray scattering, and single molecule Förster resonance energy transfer (FRET) can give the local structural information and overall dimension of IDPs, but seldom provide a unified picture of the whole protein. To understand the conformational dynamics of IDPs and how their structural ensembles recognize multiple binding partners and small molecule inhibitors, knowledge-based and physics-based sampling techniques are utilized in-silico, guided by experimental structural data. However, efficient sampling of the IDP conformational ensemble requires traversing the numerous degrees of freedom in the IDP energy landscape, as well as force-fields that accurately model the protein and solvent interactions. In this review, we have provided an overview of the current state of computational methods for studying IDP structure and dynamics and discussed the major challenges faced in this field.