Predicting molecular properties of α-synuclein using force fields for intrinsically disordered proteins

Development of a visual Adhesion/Invasion Inhibition Assay to assess the functionality of Shigella-specific antibodies

Introduction

Shigella is the etiologic agent of a bacillary dysentery known as shigellosis, which causes millions of infections and thousands of deaths worldwide each year due to Shigella's unique lifestyle within intestinal epithelial cells. Cell adhesion/invasion assays have been extensively used not only to identify targets mediating host-pathogen interaction, but also to evaluate the ability of Shigella-specific antibodies to reduce virulence. However, these assays are time-consuming and labor-intensive and fail to assess differences at the single-cell level.

Objectives and methods

Here, we developed a simple, fast and high-content method named visual Adhesion/Invasion Inhibition Assay (vAIA) to measure the ability of anti-Shigellaantibodies to inhibit bacterial adhesion to and invasion of epithelial cells by using the confocal microscope Opera Phenix.

Results

We showed that vAIA performed well with a pooled human serum from subjects challenged with S. sonnei and that a specific anti-IpaD monoclonal antibody effectively reduced bacterial virulence in a dose-dependent manner.

Discussion

vAIA can therefore inform on the functionality of polyclonal and monoclonal responses thereby supporting the discovery of pathogenicity mechanisms and the development of candidate vaccines and immunotherapies. Lastly, this assay is very versatile and may be easily applied to other Shigella species or serotypes and to different pathogens.

Using NMR diffusion data to validate MD models of disordered proteins: Test case of N-terminal tail of histone H4

MD simulations can provide uniquely detailed models of intrinsically disordered proteins (IDPs). However, these models need careful experimental validation. The coefficient of translational diffusion Dtr, measurable by pulsed field gradient NMR, offers a potentially useful piece of experimental information related to the compactness of the IDP's conformational ensemble. Here, we investigate, both experimentally and via the MD modeling, the translational diffusion of a 25-residue N-terminal fragment from histone H4 (N-H4). We found that the predicted values of Dtr, as obtained from mean-square displacement of the peptide in the MD simulations, are largely determined by the viscosity of the MD water (which has been reinvestigated as a part of our study). Beyond that, our analysis of the diffusion data indicates that MD simulations of N-H4 in the TIP4P-Ew water give rise to an overly compact conformational ensemble for this peptide. In contrast, TIP4P-D and OPC simulations produce the ensembles that are consistent with the experimental Dtr result. These observations are supported by the analyses of the 15N spin relaxation rates. We also tested a number of empirical methods to predict Dtr based on IDP's coordinates extracted from the MD snapshots. In particular, we show that the popular approach involving the program HYDROPRO can produce misleading results. This happens because HYDROPRO is not intended to predict the diffusion properties of highly flexible biopolymers such as IDPs. Likewise, recent empirical schemes that exploit the relationship between the small-angle x-ray scattering-informed conformational ensembles of IDPs and the respective experimental Dtr values also prove to be problematic. In this sense, the first-principle calculations of Dtr from the MD simulations, such as demonstrated in this work, should provide a useful benchmark for future efforts in this area.

Benchmarking of force fields to characterize the intrinsically disordered R2-FUS-LC region

Intrinsically Disordered Proteins (IDPs) play crucial roles in numerous diseases like Alzheimer's and ALS by forming irreversible amyloid fibrils. The effectiveness of force fields (FFs) developed for globular proteins and their modified versions for IDPs varies depending on the specific protein. This study assesses 13 FFs, including AMBER and CHARMM, by simulating the R2 region of the FUS-LC domain (R2-FUS-LC region), an IDP implicated in ALS. Due to the flexibility of the region, we show that utilizing multiple measures, which evaluate the local and global conformations, and combining them together into a final score are important for a comprehensive evaluation of force fields. The results suggest c36m2021s3p with mTIP3p water model is the most balanced FF, capable of generating various conformations compatible with known ones. In addition, the mTIP3P water model is computationally more efficient than those of top-ranked AMBER FFs with four-site water models. The evaluation also reveals that AMBER FFs tend to generate more compact conformations compared to CHARMM FFs but also more non-native contacts. The top-ranking AMBER and CHARMM FFs can reproduce intra-peptide contacts but underperform for inter-peptide contacts, indicating there is room for improvement.

How Cu(II) binding affects structure and dynamics of α-synuclein revealed by molecular dynamics simulations

We report accelerated molecular dynamics simulations of α-Synuclein and its complex with two Cu(II) ions bound to experimentally determined binding sites. Adding two Cu(II) ions, one bound to the N-terminal region and one to the C-terminus, decreases size and flexibility of the peptide while introducing significant new contacts within and between N-terminus and non-Aβ component (NAC). Cu(II) ions also alter the pattern of secondary structure within the peptide, inducing more and longer-lasting elements of secondary structure such as β-strands and hairpins. Free energy surfaces, obtained from reweighting the accelerated molecular dynamics boost potential, further demonstrate the restriction on size and flexibility that results from binding of copper ions.