Arginine Residues Modulate the Membrane Interactions of pHLIP Peptides.

Most processes at the water-membrane interface often involve protonation events in proteins or peptides that trigger important biological functions and events. This is the working principle behind the pHLIP peptide technology. A key titrating aspartate (Asp14 in wt) is required to protonate to induce the insertion process, increase its thermodynamic stability when membrane-embedded, and trigger the peptide's overall clinical functionality. At the core of pHLIP properties, the aspartate pKa and protonation are a consequence of the residue side chain sensing the changing surrounding environment. In this work, we characterized how the microenvironment of the key aspartate residue (Asp13 in the investigated pHLIP variants) can be modulated by a simple point mutation of a cationic residue (ArgX) at distinct sequence positions (R10, R14, R15, and R17). We carried out a multidisciplinary study using pHRE simulations and experimental measurements. Fluorescence and circular dichroism measurements were carried out to establish the stability of pHLIP variants in state III and establish the kinetics of the insertion and exit of the peptide from the membrane. We estimated the contribution of the arginine to the local electrostatic microenvironment, which promotes or hinders other electrostatic players from coexisting in the Asp interaction shell. Our data indicate that the stability and kinetics of the peptide insertion and exit from the membrane are altered when Arg is topologically available for a direct salt-bridge formation with Asp13. Hence, the position of arginine contributes to fine-tuning the pH responses of pHLIP peptides, which finds wide applications in clinics.

Corrigendum to "Unveiling molecular details behind improved activity at neutral to alkaline pH of an engineered DyP-type peroxidase" Comput Struct Biotechnol J, 20 (2022), 3899-3910.

[This corrects the article DOI: 10.1016/j.csbj.2022.07.032.].

Increasing the Realism of in Silico pHLIP Peptide Models with a Novel pH Gradient CpHMD Method.

The pH-low insertion peptides (pHLIP) are pH-dependent membrane inserting peptides, whose function depends on the cell microenvironment acidity. Several peptide variants have been designed to improve upon the wt-sequence, particularly the state transition kinetics and the selectivity for tumor pH. The variant 3 (Var3) peptide is a 27 residue long peptide, with a key titrating residue (Asp-13) that, despite showing a modest performance in liposomes (pKins ∼ 5.0), excelled in tumor cell experiments. To help rationalize these results, we focused on the pH gradient in the cell membrane, which is one of the crucial properties that are not present in liposomes. We extended our CpHMD-L method and its pH replica-exchange (pHRE) implementation to include a pH gradient and mimic the pHLIP-membrane microenvironment in a cell where the internal pH is fixed (pH 7.2) and the external pH is allowed to change. We showed that, by properly modeling the pH-gradient, we can correctly predict the experimentally observed loss and gain of performance in tumor cells experiments by the wt and Var3 sequences, respectively. In sum, the pH gradient implementation allowed for more accurate and realistic pKa estimations and was a pivotal step in bridging the in silico data and the in vivo cell experiments.

Unveiling molecular details behind improved activity at neutral to alkaline pH of an engineered DyP-type peroxidase.

DyP-type peroxidases (DyPs) are microbial enzymes that catalyze the oxidation of a wide range of substrates, including synthetic dyes, lignin-derived compounds, and metals, such as Mn2+ and Fe2+, and have enormous biotechnological potential in biorefineries. However, many questions on the molecular basis of enzyme function and stability remain unanswered. In this work, high-resolution structures of PpDyP wild-type and two engineered variants (6E10 and 29E4) generated by directed evolution were obtained. The X-ray crystal structures revealed the typical ferredoxin-like folds, with three heme access pathways, two tunnels, and one cavity, limited by three long loops including catalytic residues. Variant 6E10 displays significantly increased loops' flexibility that favors function over stability: despite the considerably higher catalytic efficiency, this variant shows poorer protein stability compared to wild-type and 29E4 variants. Constant-pH MD simulations revealed a more positively charged microenvironment near the heme pocket of variant 6E10, particularly in the neutral to alkaline pH range. This microenvironment affects enzyme activity by modulating the pK a of essential residues in the heme vicinity and should account for variant 6E10 improved activity at pH 7-8 compared to the wild-type and 29E4 that show optimal enzymatic activity close to pH 4. Our findings shed light on the structure-function relationships of DyPs at the molecular level, including their pH-dependent conformational plasticity. These are essential for understanding and engineering the catalytic properties of DyPs for future biotechnological applications.

Improved Protocol to Tackle the pH Effects on Membrane-Inserting Peptides.

Many important biological pathways rely on membrane-interacting peptides or proteins, which can alter the biophysical properties of the cell membrane by simply adsorbing to its surface to undergo a full insertion process. To study these phenomena with atomistic detail, model peptides have been used to refine the current computational methodologies. Improvements have been made with force-field parameters, enhanced sampling techniques to obtain faster sampling, and the addition of chemical-physical properties, such as pH, whose influence dramatically increases at the water/membrane interface. The pH (low) insertion peptide (pHLIP) is a peptide that inserts across a membrane bilayer depending on the pH due to the presence of a key residue (Asp14) whose acidity-induced protonation triggers the whole process. The complex nature of these peptide/membrane interactions resulted in sampling limitations of the protonation and configurational space albeit using state-of-the-art methods such as the constant-pH molecular dynamics. To address this issue and circumvent those limitations, new simulations were performed with our newly developed pH-replica exchange method using wild-type (wt)-pHLIP in different 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine membrane sizes. This technique provided enhanced sampling and allowed for the calculation of more complete Asp14 pKa profiles. The conformational heterogeneity derived from strong electrostatic interactions between Asp14 and the lipid phosphate groups was identified as the source of most pKa variability. In spite of these persistent and harder-to-equilibrate phosphate interactions, the pKa values at deeper regions (6.0-6.2) still predicted the experimental pK of insertion (6.0) since the electrostatic perturbation decays as the residue inserts further into the membrane. We also observed that reducing the system size leads to membrane deformations where it increasingly loses the ability to accommodate the pHLIP-induced perturbations. This indicates that large membrane patches, such as 256 or even 352 lipids, are needed to obtain stable and more realistic pHLIP/membrane systems. These results strengthen our method pKa predictive and analytical capabilities to study the intricate play of electrostatic effects of the peptide/membrane interface, granting confidence for future applications in similar systems.

pK a Calculations in Membrane Proteins from Molecular Dynamics Simulations.

The conformational changes of membrane proteins are crucial to their function and usually lead to fluctuations in the electrostatic environment of the protein surface. A very effective way to quantify these changes is by calculating the pK a values of the protein's titratable residues, which can be regarded as electrostatic probes. To achieve this, we need to take advantage of the fast and reliable pK a calculators developed for globular proteins and adapt them to include the explicit effects of membranes. Here, we provide a detailed linear response approximation protocol that uses our own software (PypKa) to calculate reliable pK a values from short MD simulations of membrane proteins.

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.

The Lysosomotropic Activity of Hydrophobic Weak Base Drugs is Mediated via Their Intercalation into the Lysosomal Membrane.

Lipophilic weak base therapeutic agents, termed lysosomotropic drugs (LDs), undergo marked sequestration and concentration within lysosomes, hence altering lysosomal functions. This lysosomal drug entrapment has been described as luminal drug compartmentalization. Consistent with our recent finding that LDs inflict a pH-dependent membrane fluidization, we herein demonstrate that LDs undergo intercalation and concentration within lysosomal membranes. The latter was revealed experimentally and computationally by (a) confocal microscopy of fluorescent compounds and drugs within lysosomal membranes, and (b) molecular dynamics modeling of the pH-dependent membrane insertion and accumulation of an assortment of LDs, including anticancer drugs. Based on the multiple functions of the lysosome as a central nutrient sensory hub and a degradation center, we discuss the molecular mechanisms underlying the alteration of morphology and impairment of lysosomal functions as consequences of LDs' intercalation into lysosomes. Our findings bear important implications for drug design, drug induced lysosomal damage, diseases and pertaining therapeutics.

Membrane-Induced p Ka Shifts in wt-pHLIP and Its L16H Variant.

The pH (low) insertion peptides (pHLIPs) is a family of peptides that are able to insert into a lipid bilayer at acidic pH. The molecular mechanism of pHLIPs insertion, folding, and stability in the membrane at low pH is based on multiple protonation events, which are challenging to study at the molecular level. More specifically, the relation between the experimental p K of insertion (p Kexp) of pHLIPs and the p Ka of the key residues is yet to be clarified. We carried out a computational study, complemented with new experimental data, and established the influence of (de)protonation of titrable residues on the stability of the peptide membrane-inserted state. Constant-pH molecular dynamics simulations were employed to calculate the p Ka values of these residues along the membrane normal. In the wt-pHLIP, we identified Asp14 as the key residue for the stability of the membrane-inserted state, and its p Ka value is strongly correlated with the experimental p Kexp measured in thermodynamics studies. Also, in order to narrow down the pH range at which pHLIP is stable in the membrane, we designed a new pHLIP variant, L16H, where Leu in the 16th position was replaced by a titrable His residue. Our results showed that the L16H variant undergoes two transitions. The calculated p Ka and experimentally observed p Kexp values are in good agreement. Two distinct p Kexp values delimit a pH range where the L16H peptide is stably inserted in the membrane, while, outside this range, the membrane-inserted state is destabilized and the peptide exits from the bilayer. pHLIP peptides have been successfully used to target cancer cells for the delivery of diagnostics and therapeutic agents to acidic tumors. The fine-tuning of the stability of the pHLIP inserted state and its restriction to a narrow well-defined pH range might allow the design of new peptides, able to discriminate between tissues with different extracellular pH values.

The Impact of Using Single Atomistic Long-Range Cutoff Schemes with the GROMOS 54A7 Force Field.

With the recent increase in computing power, the molecular modeling community is now more focused on improving the accuracy and overall quality of biomolecular simulations. For the available simulation packages, force fields, and all other associated methods used, this relates to how well they describe the conformational space and thermodynamic properties of a biomolecular system. The parameter sets of GROMOS force fields have been parametrized and validated with the reaction field (RF) method using charge groups and a twin-range cutoff scheme (0.8/1.4 nm). However, the most recent versions of GROMACS (since v.2016) discontinued the support for charge groups. To take full advantage of the newer and faster versions of this software package with GROMOS 54A7 and RF, we need to evaluate the impact of using a single cutoff scheme (vs twin-range) and of using the Verlet list update method (which is atomistic) compared to the group-based cutoff scheme. Our results show that the GROMOS 54A7 force field seems consistent with a single cutoff, since the resulting conformation and protonation ensembles were indistinguishable. The GROMOS parametrization procedure was also reproduced using an atomistic cutoff scheme, and we have observed that the hydration free energy values of small amino acid side-chain analogues were similar to the ones obtained with the group-based protocol. We do observe a small impact of the atomistic cutoff scheme in the conformational space of the model systems studied (G1-PAMAM and DMPC). However, since the structural properties of these systems are well converged for the cutoff range used (1.4-2.0 nm), unlike with the group-based cutoff schemes, we are confident that the atomistic cutoff can be adopted with RF for MD and constant-pH MD biomolecular simulations using the GROMOS 54A7 force field.