Improved GROMOS 54A7 Charge Sets for Phosphorylated Tyr, Ser, and Thr to Deal with pH-Dependent Binding Phenomena

Histidine Protonation and Conformational Switching in Diphtheria Toxin Translocation Domain

Protonation of key histidine residues has been long implicated in the acid-mediated cellular action of the diphtheria toxin translocation (T-) domain, responsible for the delivery of the catalytic domain into the cell. Here, we use a combination of computational (constant-pH Molecular Dynamics simulations) and experimental (NMR, circular dichroism, and fluorescence spectroscopy along with the X-ray crystallography) approaches to characterize the initial stages of conformational change happening in solution in the wild-type T-domain and in the H223Q/H257Q double mutant. This replacement suppresses the acid-induced transition, resulting in the retention of a more stable protein structure in solutions at pH 5.5 and, consequently, in reduced membrane-disrupting activity. Here, for the first time, we report the pKa values of the histidine residues of the T-domain, measured by NMR-monitored pH titrations. Most peaks in the histidine side chain spectral region are titrated with pKas ranging from 6.2 to 6.8. However, the two most up-field peaks display little change down to pH 6, which is a limiting pH for this protein in solution at concentrations required for NMR. These peaks are absent in the double mutant, suggesting they belong to H223 and H257. The constant-pH simulations indicate that for the T-domain in solution, the pKa values for histidine residues range from 3.0 to 6.5, with those most difficult to protonate being H251 and H257. Taken together, our experimental and computational data demonstrate that previously suggested cooperative protonation of all six histidines in the T-domain does not occur.

Phosphorylation Modification Force Field FB18CMAP Improving Conformation Sampling of Phosphoproteins

Phosphorylation of proteins plays an important regulatory role at almost all levels of cellular organization. Molecular dynamics (MD) simulation is a promising tool to reveal the mechanism of how phosphorylation regulates many key biological processes at the atomistic level. MD simulation accuracy depends on force field precision, while the current force fields for phospho-amino acids have resulted in notable inconsistency with experimental data. Here, a new force field parameter (named FB18CMAP) is generated by fitting against quantum mechanics (QM) energy in aqueous solution with φ/ψ dihedral potential-energy surfaces optimized using CMAP parameters. MD simulations of phosphorylated dipeptides, intrinsically disordered proteins (IDPs), and ordered (folded) proteins show that FB18CMAP can mimic NMR observables and structural characteristics of phosphorylated dipeptides and proteins more accurately than the FB18 force field. These findings suggest that FB18CMAP performs well in both the simulation of ordered and disordered states of phosphorylated proteins.

Constant pH molecular dynamics simulations: Current status and recent applications

Many important protein functions are carried out through proton-coupled conformational dynamics. Thus, the ability to accurately model protonation states dynamically has wide-ranging implications. Over the past two decades, two main types of constant pH methods (discrete and continuous) have been developed to enable proton-coupled molecular dynamics (MD) simulations. In this short review, we discuss the current status of the development and highlight recent applications that have advanced our understanding of protein structure-function relationships. We conclude the review by outlining the remaining challenges in the method development and projecting important areas for future applications.

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 (pK^ins ∼ 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 pK_a estimations and was a pivotal step in bridging the in silico data and the in vivo cell experiments.

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 (pK_a 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.

Novel US-CpHMD Protocol to Study the Protonation-Dependent Mechanism of the ATP/ADP Carrier

We have designed a protocol combining constant-pH molecular dynamics (CpHMD) simulations with an umbrella sampling (US) scheme (US-CpHMD) to study the mechanism of ADP/ATP transport (import and export) by their inner mitochondrial membrane carrier protein [ADP/ATP carrier (AAC)]. The US scheme helped overcome the limitations of sampling the slow kinetics involved in these substrates' transport, while CpHMD simulations provided an unprecedented realism by correctly capturing the associated protonation changes. The import of anionic substrates along the mitochondrial membrane has a strong energetic disadvantage due to a smaller substrate concentration and an unfavorable membrane potential. These limitations may have created an evolutionary pressure on AAC to develop specific features benefiting the import of ADP. In our work, the potential of mean force profiles showed a clear selectivity in the import of ADP compared to ATP, while in the export, no selectivity was observed. We also observed that AAC sequestered both substrates at longer distances in the import compared to the export process. Furthermore, only in the import process do we observe transient protonation of both substrates when going through the AAC cavity, which is an important advantage to counteract the unfavorable mitochondrial membrane potential. Finally, we observed a substrate-induced disruption of the matrix salt-bridge network, which can promote the conformational transition (from the C- to M-state) required to complete the import process. This work unraveled several important structural features where the complex electrostatic interactions were pivotal to interpreting the protein function and illustrated the potential of applying the US-CpHMD protocol to other transport processes involving membrane proteins.

Benzo[a]pyrene and heavy metal ion adsorption on nanoplastics regulated by humic acid: Cooperation/competition mechanisms revealed by molecular dynamics simulations

Nanoplastics adsorb pollutants and organic matter to aggravate or alleviate impact to the eco-environment and human health. However, the interaction mechanisms remain unclear and difficult to study using current experimental techniques. By means of molecular dynamics simulation, here we investigate adsorption of benzo[a]pyrene (BaP) and heavy metal ions (Cu²⁺) on nanoplastics of different materials and surface charges regulated by humic acid (HA). Among considered materials, polystyrene shows the highest capacity of adsorbing BaPs via forming sandwiched π-stacking structures with benzene rings. Driven by hydrophobic, electrostatic and hydrogen bonding interactions, HAs spontaneously aggregate into micelle-like structures with hydrophobic core and charged exterior accessible to BaPs and Cu²⁺, respectively. Cationic and neutral nanoplastics adsorb more HAs to form eco-coronas, which modulate BaP and Cu²⁺ adsorption via following cooperation/competition mechanisms. On one hand, the direct binding of BaPs to nanoplastics is hindered by HAs through BaP encapsulation plus competitive adsorption. On the other hand, adsorbed HAs expose carboxyl groups to offer rich binding sites to promote Cu²⁺ adsorption on neutral and cationic nanoplastics, while unbound HAs compete with anionic nanoplastics to inhibit Cu²⁺ adsorption. These results provide molecular level insights into transport, transformation and accessibility of nanoplastics with coexisting contaminants in the aqueous environment.

Development and Validation of AMBER-FB15-Compatible Force Field Parameters for Phosphorylated Amino Acids

Phosphorylation of select amino acid residues is one of the most common biological mechanisms for regulating protein structures and functions. While computational modeling can be used to explore the detailed structural changes associated with phosphorylation, most molecular mechanics force fields developed for the simulation of phosphoproteins have been noted to be inconsistent with experimental data. In this work, we parameterize force fields for the phosphorylated forms of the amino acids serine, threonine, and tyrosine using the ForceBalance software package with the goal of improving agreement with experiments for these residues. Our optimized force field, denoted as FB18, is parameterized using high-quality ab initio potential energy scans and is designed to be fully compatible with the AMBER-FB15 protein force field. When utilized in MD simulations together with the TIP3P-FB water model, we find that FB18 consistently enhances the prediction of experimental quantities such as ³J NMR couplings and intramolecular hydrogen-bonding propensities in comparison to previously published models. As was reported with AMBER-FB15, we also see improved agreement with the reference QM calculations in regions at and away from local minima. We thus believe that the FB18 parameter set provides a promising route for the further investigation of the varied effects of protein phosphorylation.

Designing new antitubercular isoniazid derivatives with improved reactivity and membrane trafficking abilities

Isoniazid (INH) is one of the two most effective first-line antitubercular drugs and is still used at the present time as a scaffold for developing new compounds to fight TB. In a previous study, we have observed that an INH derivative, an hydrazide N'-substituted with a C₁₀acyl chain, was able to counterbalance its smaller reactivity with a higher membrane permeability. This resulted in an improved performance against the most prevalent Mycobacterium tuberculosis (Mtb) resistant strain (S315T), compared to INH. In this work, we have designed two new series of INH derivatives (alkyl hydrazides and hydrazones) with promising in silico properties, namely membrane permeabilities and spontaneous IN* radical formation. The kinetics, cytotoxicity, and biological activity evaluations confirmed the in silico predictions regarding the very high reactivity of the alkyl hydrazides. The hydrazones, on the other hand, showed very similar behavior compared to INH, particularly in biological tests that take longer to complete, indicating that these compounds are being hydrolyzed back to INH. Despite their improved membrane permeabilities, the reactivities of these two series are too high, impairing their overall performance. Nevertheless, the systematic data gathered about these compounds have showed us the need to find a balance between lipophilicity and reactivity, which is paramount to devise better INH-based derivatives aimed at circumventing Mtb resistance.