Pyroglutamylation modulates electronic properties and the conformational ensemble of the amyloid ß-peptide.

Alzheimer's disease (AD) is a neurodegenerative disorder that is characterized by the formation of extracellular amyloid-ß (Aß) plaques. The underlying cause of AD is unknown, however, post-translational modifications (PTMs) of Aß have been found in AD patients and are thought to play a role in protein aggregation. One such PTM is pyroglutamylation, which can occur at two sites in Aß, Glu3 and Glu11. This modification of Aß involves the truncation and charge-neutralization of N-terminal glutamate, causing Aß to become more hydrophobic and prone to aggregation. The molecular mechanism by which the introduction of pyroglutamate (pE) promotes aggregation has not been determined. To gain a greater understanding of the role that charge neutralization and truncation of the N-terminus plays on Aß conformational sampling, we used the Drude polarizable force field (FF) to perform molecular dynamics simulations on AßpE3-42 and AßpE11-42 and comparing their properties to previous simulations of Aß1-42. The Drude polarizable FF allows for a more accurate representation of electrostatic interactions, therefore providing novel insights into the role that charge plays in protein dynamics. Here, we report the parametrization of pE in the Drude polarizable FF and the effect of pyroglutamylation on Aß. We found that AßpE3-42 and AßpE11-42 alter the permanent and induced dipoles of the peptide. Specifically, we found that AßpE3-42 and AßpE11-42 have modification-specific backbone and sidechain polarization response and perturbed solvation properties that shift the Aß conformational ensemble.

Monomeric and dimeric states of human ZO1-PDZ2 are functional partners of the SARS-CoV-2 E protein.

The Envelope (E) protein of SARS-CoV-2 plays a key role in virus maturation, assembly, and virulence mechanisms. The E protein is characterized by the presence of a PDZ-binding motif (PBM) at its C-terminus that allows it to interact with several PDZ-containing proteins in the intracellular environment. One of the main binding partners of the SARS-CoV-2 E protein is the PDZ2 domain of ZO1, a protein with a crucial role in the formation of epithelial and endothelial tight junctions (TJs). In this work, through a combination of analytical ultracentrifugation analysis and equilibrium and kinetic folding experiments, we show that ZO1-PDZ2 domain is able to fold in a monomeric state, an alternative form to the dimeric conformation that is reported to be functional in the cell for TJs assembly. Importantly, surface plasmon resonance (SPR) data indicate that the PDZ2 monomer is fully functional and capable of binding the C-terminal portion of the E protein of SARS-CoV-2, with a measured affinity in the micromolar range. Moreover, we present a detailed computational analysis of the complex between the C-terminal portion of E protein with ZO1-PDZ2, both in its monomeric conformation (computed as a high confidence AlphaFold2 model) and dimeric conformation (obtained from the Protein Data Bank), by using both polarizable and nonpolarizable simulations. Together, our results indicate both the monomeric and dimeric states of PDZ2 to be functional partners of the E protein, with similar binding mechanisms, and provide mechanistic and structural information about a fundamental interaction required for the replication of SARS-CoV-2.

Effects of Familial Alzheimer's Disease Mutations on the Folding Free Energy and Dipole-Dipole Interactions of the Amyloid ß-Peptide.

Familial Alzheimer's disease (FAD) mutations of the amyloid ß-peptide (Aß) are known to lead to early onset and more aggressive Alzheimer's disease. FAD mutations such as "Iowa" (D23N), "Arctic" (E22G), "Italian" (E22K), and "Dutch" (E22Q) have been shown to accelerate Aß aggregation relative to the wild-type (WT). The mechanism by which these mutations facilitate increased aggregation is unknown, but each mutation results in a change in the net charge of the peptide. Previous studies have used nonpolarizable force fields to study Aß, providing some insight into how this protein unfolds. However, nonpolarizable force fields have fixed charges that lack the ability to redistribute in response to changes in local electric fields. Here, we performed polarizable molecular dynamics simulations on the full-length Aß42 of WT and FAD mutations and calculated folding free energies of the Aß15-27 fragment via umbrella sampling. By studying both the full-length Aß42 and a fragment containing mutations and the central hydrophobic cluster (residues 17-21), we were able to systematically study how these FAD mutations impact secondary and tertiary structure and the thermodynamics of folding. Electrostatic interactions, including those between permanent and induced dipoles, affected side-chain properties, salt bridges, and solvent interactions. The FAD mutations resulted in shifts in the electronic structure and solvent accessibility at the central hydrophobic cluster and the hydrophobic C-terminal region. Using umbrella sampling, we found that the folding of the WT and E22 mutants is enthalpically driven, whereas the D23N mutant is entropically driven, arising from a different unfolding pathway and peptide-bond dipole response. Together, the unbiased, full-length, and umbrella sampling simulations of fragments reveal that the FAD mutations perturb nearby residues and others in hydrophobic regions to potentially alter solubility. These results highlight the role electronic polarizability plays in amyloid misfolding and the role of heterogeneous microenvironments that arise as conformational change takes place.

Impact of Electronic Polarization on Preformed, ß-Strand Rich Homogenous and Heterogenous Amyloid Oligomers.

Amyloids are a subset of intrinsically disordered proteins (IDPs) that self-assemble into cross-ß oligomers and fibrils. The structural plasticity of amyloids leads to sampling of metastable, low-molecular-weight oligomers that contribute to cytotoxicity. Of interest are amyloid-ß (Aß) and islet amyloid polypeptide (IAPP), which are involved in the pathology of Alzheimer's disease and Type 2 Diabetes Mellitus, respectively. In addition to forming homogenous oligomers and fibrils, these species have been found to cross-aggregate in heterogeneous structures. Biophysical properties, including electronic effects, that are unique or conserved between homogenous and heterogenous amyloids oligomers are thus far unexplored. Here, we simulated homogenous and heterogenous amyloid oligomers of Aß16-22 and IAPP20-29 fragments using the Drude oscillator model to investigate the impact of electronic polarization on the structural morphology and stability of preformed hexamers. Upon simulation of preformed, ß-strand rich oligomers with Drude, structural rearrangement occurred causing some loss of ß-strand structure in favor of random coil content for all oligomers. Homogenous Aß16-22 was the most stable system, deriving stability from low polarization in hydrophobic residues and through salt bridge formation. Changes in polarization were observed primarily for Aß16-22 residues in heterogenous cross-amyloid systems, displaying a decrease in charged residue dipole moments and an increase in hydrophobic sidechain dipole moments. This work is the first study utilizing the Drude-2019 force field with amyloid oligomers, providing insight into the impact of electronic effects on oligomer structure and highlighting the importance of different microenvironments on amyloid oligomer stability.

Insights into Stabilizing Forces in Amyloid Fibrils of Differing Sizes from Polarizable Molecular Dynamics Simulations.

Pathological aggregation of amyloid-forming proteins is a hallmark of a number of human diseases, including Alzheimer's, type 2 diabetes, Parkinson's, and more. Despite having very different primary amino acid sequences, these amyloid proteins form similar supramolecular, fibril structures that are highly resilient to physical and chemical denaturation. To better understand the structural stability of disease-related amyloids and to gain a greater understanding of factors that stabilize functional amyloid assemblies, insights into tertiary and quaternary interactions are needed. We performed molecular dynamics simulations on human tau, amyloid-ß, and islet amyloid polypeptide fibrils to determine key physicochemical properties that give rise to their unique characteristics and fibril structures. These simulations are the first of their kind in employing a polarizable force field to explore properties of local electric fields on dipole properties and other electrostatic forces that contribute to amyloid stability. Across these different amyloid fibrils, we focused on how the underlying forces stabilize fibrils to elucidate the driving forces behind the protein aggregation. The polarizable model allows for an investigation of how side-chain dipole moments, properties of structured water molecules in the fibril core, and the local environment around salt bridges contribute to the formation of interfaces essential for fibril stability. By systematically studying three amyloidogenic proteins of various fibril sizes for key structural properties and stabilizing forces, we shed light on properties of amyloid structures related to both diseased and functional states at the atomistic level.