High Charge Density in Peptide Dendrimers is Required to Destabilize Membranes: Insights into Endosome Evasion.

Peptide dendrimers are a type of branched, symmetric, and topologically well-defined molecule that have already been used as delivery systems for nucleic acid transfection. Several of the most promising sequences showed high efficiency in many key steps of transfection, namely, binding siRNA, entering cells, and evading the endosome. However, small changes to the peptide dendrimers, such as in the hydrophobic core, the amino acid chirality, or the total available charges, led to significantly different experimental results with unclear mechanistic insights. In this work, we built a computational model of several of those peptide dendrimers (MH18, MH13, and MH47) and some of their variants to study the molecular details of the structure and function of these molecules. We performed CpHMD simulations in the aqueous phase and in interaction with a lipid bilayer to assess how conformation and protonation are affected by pH in different environments. We found that while the different peptide dendrimer sequences lead to no substantial structural differences in the aqueous phase, the total charge and, more importantly, the total charge density are key for the capacity of the dendrimer to interact and destabilize the membrane. These dendrimers become highly charged when the pH changes from 7.5 to 4.5, and the presence of a high charge density, which is decreased for MH47 that has four fewer titratable lysines, is essential to trigger membrane destabilization. These findings are in excellent agreement with the experimental data and help us to understand the high efficiency of some dendrimers and why the dendrimer MH47 is unable to complete the transfection process. This evidence provides further understanding of the mode of action of these peptide dendrimers and will be pivotal for the future design of new sequences with improved transfection capabilities.

Interfacial Dynamics and Growth Modes of ß2-Microglobulin Dimers.

Protein aggregation is a complex process, strongly dependent on environmental conditions and highly structurally heterogeneous, both at the final level of fibril structure and intermediate level of oligomerization. Since the first step in aggregation is the formation of a dimer, it is important to clarify how certain properties of the latter (e.g., stability or interface geometry) may play a role in self-association. Here, we report a simple model that represents the dimer's interfacial region by two angles and combine it with a simple computational method to investigate how modulations of the interfacial region occurring on the ns-µs time scale change the dimer's growth mode. To illustrate the proposed methodology, we consider 15 different dimer configurations of the ß2m D76N mutant protein equilibrated with long Molecular Dynamics simulations and identify which interfaces lead to limited and unlimited growth modes, having, therefore, different aggregation profiles. We found that despite the highly dynamic nature of the starting configurations, most polymeric growth modes tend to be conserved within the studied time scale. The proposed methodology performs remarkably well taking into consideration the nonspherical morphology of the ß2m dimers, which exhibit unstructured termini detached from the protein's core, and the relatively weak binding affinities of their interfaces, which are stabilized by nonspecific apolar interactions. The proposed methodology is general and can be applied to any protein for which a dimer structure has been experimentally determined or computationally predicted.

Extending the Stochastic Titration CpHMD to CHARMM36m.

The impact of pH on proteins is significant but often neglected in molecular dynamics simulations. Constant-pH Molecular Dynamics (CpHMD) is the state-of-the-art methodology to deal with these effects. However, it still lacks widespread adoption by the scientific community. The stochastic titration CpHMD is one of such methods that, until now, only supported the GROMOS force field family. Here, we extend this method's implementation to include the CHARMM36m force field available in the GROMACS software package. We test this new implementation with a diverse group of proteins, namely, lysozyme, Staphylococcal nuclease, and human and E. coli thioredoxins. All proteins were conformationally stable in the simulations, even at extreme pH values. The RMSE values (pKa prediction vs experimental) obtained were very encouraging, in particular for lysozyme and human thioredoxin. We have also identified a few residues that challenged the CpHMD simulations, highlighting scenarios where the method still needs improvement independently of the force field. The CHARMM36m all-atom implementation was more computationally efficient when compared with the GROMOS 54A7, taking advantage of a shorter nonbonded interaction cutoff and a less frequent neighboring list update. The new extension will allow the study of pH effects in many systems for which this force field is particularly suited, i.e., proteins, membrane proteins, lipid bilayers, and nucleic acids.

In Silico End-to-End Protein-Ligand Interaction Characterization Pipeline: The Case of SARS-CoV-2.

SARS-CoV-2 triggered a worldwide pandemic disease, COVID-19, for which an effective treatment has not yet been settled. Among the most promising targets to fight this disease is SARS-CoV-2 main protease (Mpro), which has been extensively studied in the last few months. There is an urgency for developing effective computational protocols that can help us tackle these key viral proteins. Hence, we have put together a robust and thorough pipeline of in silico protein-ligand characterization methods to address one of the biggest biological problems currently plaguing our world. These methodologies were used to characterize the interaction of SARS-CoV-2 Mpro with an α-ketoamide inhibitor and include details on how to upload, visualize, and manage the three-dimensional structure of the complex and acquire high-quality figures for scientific publications using PyMOL (Protocol 1); perform homology modeling with MODELLER (Protocol 2); perform protein-ligand docking calculations using HADDOCK (Protocol 3); run a virtual screening protocol of a small compound database of SARS-CoV-2 candidate inhibitors with AutoDock 4 and AutoDock Vina (Protocol 4); and, finally, sample the conformational space at the atomic level between SARS-CoV-2 Mpro and the α-ketoamide inhibitor with Molecular Dynamics simulations using GROMACS (Protocol 5). Guidelines for careful data analysis and interpretation are also provided for each Protocol.

Predicting stable binding modes from simulated dimers of the D76N mutant of ß 2-microglobulin.

The D76N mutant of the ß 2 m protein is a biologically motivated model system to study protein aggregation. There is strong experimental evidence, supported by molecular simulations, that D76N populates a highly dynamic conformation (which we originally named I 2 ) that exposes aggregation-prone patches as a result of the detachment of the two terminal regions. Here, we use Molecular Dynamics simulations to study the stability of an ensemble of dimers of I 2 generated via protein-protein docking. MM-PBSA calculations indicate that within the ensemble of investigated dimers the major contribution to interface stabilization at physiological pH comes from hydrophobic interactions between apolar residues. Our structural analysis also reveals that the interfacial region associated with the most stable binding modes are particularly rich in residues pertaining to both the N- and C-terminus, as well residues from the BC- and DE-loops. On the other hand, the less stable interfaces are stabilized by intermolecular interactions involving residues from the CD- and EF-loops. By focusing on the most stable binding modes, we used a simple geometric rule to propagate the corresponding dimer interfaces. We found that, in the absence of any kind of structural rearrangement occurring at an early stage of the oligomerization pathway, some interfaces drive a self-limited growth process, while others can be propagated indefinitely allowing the formation of long, polymerized chains. In particular, the interfacial region of the most stable binding mode reported here falls in the class of self-limited growth.

Computational Analysis of the Interactions between the S100B Extracellular Chaperone and Its Amyloid ß Peptide Client.

S100B is an astrocytic extracellular Ca2+-binding protein implicated in Alzheimer's disease, whose role as a holdase-type chaperone delaying Aß42 aggregation and toxicity was recently uncovered. Here, we employ computational biology approaches to dissect the structural details and dynamics of the interaction between S100B and Aß42. Driven by previous structural data, we used the Aß25-35 segment, which recapitulates key aspects of S100B activity, as a starting guide for the analysis. We used Haddock to establish a preferred binding mode, which was studied with the full length Aß using long (1 µs) molecular dynamics (MD) simulations to investigate the structural dynamics and obtain representative interaction complexes. From the analysis, Aß-Lys28 emerged as a key candidate for stabilizing interactions with the S100B binding cleft, in particular involving a triad composed of Met79, Thr82 and Glu86. Binding constant calculations concluded that coulombic interactions, presumably implicating the Lys28(Aß)/Glu86(S100B) pair, are very relevant for the holdase-type chaperone activity. To confirm this experimentally, we examined the inhibitory effect of S100B over Aß aggregation at high ionic strength. In agreement with the computational predictions, we observed that electrostatic perturbation of the Aß-S100B interaction decreases anti-aggregation activity. Altogether, these findings unveil features relevant in the definition of selectivity of the S100B chaperone, with implications in Alzheimer's disease.

Cu2+-binding to S100B triggers polymerization of disulfide cross-linked tetramers with enhanced chaperone activity against amyloid-ß aggregation.

S100B is an extracellular protein implicated in Alzheimer's Disease and a suppressor of amyloid-ß aggregation. Herein we report a mechanism tying Cu2+ binding to a change in assembly state yielding disulfide cross-linked oligomers with higher anti-aggregation activity. This chemical control of chaperone function illustrates a regulatory process relevant under metal and proteostasis dysfunction as in neurodegeneration.