Alchemical Calculation of Relative Free Energies for Charge-Changing Mutations at Protein-Protein Interfaces Considering Fixed and Variable Protonation States.

The calculation of relative free energies (ΔΔG) for charge-changing mutations at protein-protein interfaces through alchemical methods remains challenging due to variations in the system's net charge during charging steps, the possibility of mutated and contacting ionizable residues occurring in various protonation states, and undersampling issues. In this study, we present a set of strategies, collectively termed TIRST/TIRST-H+, to address some of these challenges. Our approaches combine thermodynamic integration (TI) with the prediction of pKa shifts to calculate ΔΔG values. Moreover, special sets of restraints are employed to keep the alchemically transformed molecules separated. The accuracy of the devised approaches was assessed on a large and diverse data set comprising 164 point mutations of charged residues (Asp, Glu, Lys, and Arg) to Ala at the protein-protein interfaces of complexes with known three-dimensional structures. Mean absolute and root-mean-square errors ranging from 1.38 to 1.66 and 1.89 to 2.44 kcal/mol, respectively, and Pearson correlation coefficients of ∼0.6 were obtained when testing the approaches on the selected data set using the GPU-TI module of Amber18 suite and the ff14SB force field. Furthermore, the inclusion of variable protonation states for the mutated acid residues improved the accuracy of the predicted ΔΔG values. Therefore, our results validate the use of TIRST/TIRST-H+ in prospective studies aimed at evaluating the impact of charge-changing mutations to Ala on the stability of protein-protein complexes.

Kinetic and structural studies of Mycobacterium tuberculosis dihydroorotate dehydrogenase reveal new insights into class 2 DHODH inhibition.

Tuberculosis (TB) is a leading cause of death worldwide. TB represents a serious public health threat, and it is characterized by high transmission rates, prevalence in impoverished regions, and high co-infection rates with HIV. Moreover, the serious side effects of long-term treatment that decrease patient adherence, and the emergence of multi-resistant strains of Mycobacterium tuberculosis, the causing agent of TBs, pose several challenges for its eradication. The search for a new TB treatment is necessary and urgent. Dihydroorotate dehydrogenase (DHODH) is responsible for the stereospecific oxidation of (S)-dihydroorotate (DHO) to orotate during the fourth and only redox step of the de novo pyrimidine nucleotide biosynthetic pathway. DHODH has been considered an attractive target against infectious diseases. As a first step towards exploiting DHODH as a drug target against TB, we performed a full kinetic characterization of both bacterial MtDHODH and its human ortholog (HsDHDOH) using both substrates coenzyme Q0 (Q0) and vitamin K3 (K3). MtDHODH follows a ping-pong mechanism of catalysis and shares similar catalytic parameters with the human enzyme. Serendipitously, Q0 was found to inhibit MtDHODH (KI (Q0) = 138 ± 31 µM). To the best of our knowledge, Q0 is the first non-orotate like dihydroorotate-competitive inhibitor for class 2 DHODHs ever described. Molecular dynamics simulations along with in silico solvent mapping allowed us to successfully probe protein flexibility and correlate it with the druggability of binding sites. Together, our results provide the starting point for the design of a new generation of potent and selective inhibitors against MtDHODH.

Computational study on the allosteric mechanism of Leishmania major IF4E-1 by 4E-interacting protein-1: Unravelling the determinants of m7GTP cap recognition.

During their life cycle, Leishmania parasites display a fine-tuned regulation of the mRNA translation through the differential expression of isoforms of eukaryotic translation initiation factor 4E (LeishIF4Es). The interaction between allosteric modulators such as 4E-interacting proteins (4E-IPs) and LeishIF4E affects the affinity of this initiation factor for the mRNA cap. Here, several computational approaches were employed to elucidate the molecular bases of the previously-reported allosteric modulation in L. major exerted by 4E-IP1 (Lm4E-IP1) on eukaryotic translation initiation factor 4E 1 (LmIF4E-1). Molecular dynamics (MD) simulations and accurate binding free energy calculations (ΔGbind ) were combined with network-based modeling of residue-residue correlations. We also describe the differences in internal motions of LmIF4E-1 apo form, cap-bound, and Lm4E-IP1-bound systems. Through community network calculations, the differences in the allosteric pathways of allosterically-inhibited and active forms of LmIF4E-1 were revealed. The ΔGbind values show significant differences between the active and inhibited systems, which are in agreement with the available experimental data. Our study thoroughly describes the dynamical perturbations of LmIF4E-1 cap-binding site triggered by Lm4E-IP1. These findings are not only essential for the understanding of a critical process of trypanosomatids' gene expression but also for gaining insight into the allostery of eukaryotic IF4Es, which could be useful for structure-based design of drugs against this protein family.

The effect of acidic pH on the adsorption and lytic activity of the peptides Polybia-MP1 and its histidine-containing analog in anionic lipid membrane: a biophysical study by molecular dynamics and spectroscopy.

Antimicrobial peptides (AMPs) are part of the innate immune system of many species. AMPs are short sequences rich in charged and non-polar residues. They act on the lipid phase of the plasma membrane without requiring membrane receptors. Polybia-MP1 (MP1), extracted from a native wasp, is a broad-spectrum bactericide, an inhibitor of cancer cell proliferation being non-hemolytic and non-cytotoxic. MP1 mechanism of action and its adsorption mode is not yet completely known. Its adsorption to lipid bilayer and lytic activity is most likely dependent on the ionization state of its two acidic and three basic residues and consequently on the bulk pH. Here we investigated the effect of bulk acidic (pH 5.5) and neutral pH (7.4) solution on the adsorption, insertion, and lytic activity of MP1 and its analog H-MP1 to anionic (7POPC:3POPG) model membrane. H-MP1 is a synthetic analog of MP1 with lysines replaced by histidines. Bulk pH changes could modulate this peptide efficiency. The combination of different experimental techniques and molecular dynamics (MD) simulations showed that the adsorption, insertion, and lytic activity of H-MP1 are highly sensitive to bulk pH in opposition to MP1. The atomistic details, provided by MD simulations, showed peptides contact their N-termini to the bilayer before the insertion and then lay parallel to the bilayer. Their hydrophobic faces inserted into the acyl chain phase disturb the lipid-packing.

Modulatory Effects of Acidic pH and Membrane Potential on the Adsorption of pH-Sensitive Peptides to Anionic Lipid Membrane.

Anionic lipid membrane electrostatic potential and solution pH can influence cationic peptide adsorption to these bilayers, especially those containing simultaneously acid and basic residues. Here, we investigate the effects of the pH solution on MP1 (IDWKKLLDAAKQIL-NH2) adsorption to anionic (7POPC:3POPG) lipid vesicles in comparison to its analog H-MP1, with histidines substituting lysines. We used the association of adsorption isotherms and constant pH molecular dynamic simulations (CpHMD) to explore the effects of membrane potential and pH on peptides' adsorption on this lipid membrane. We analyzed the fluorescence and zeta potential adsorption isotherms using the Gouy-Chapman theory. In CpHMD simulations for the peptides in solution and adsorbed on the lipid bilayer, we used the conformations obtained by conventional MD simulations at a µs timescale. Non-equilibrium Monte Carlo simulations provided the protonation states of acidic and basic residues. CpHMD showed average pKa shifts of two to three units, resulting in a higher net charge for the analog than for MP1, strongly modulating the peptide adsorption. The fractions of the protonation of acidic and basic residues and the peptides' net charges obtained from the analysis of the adsorption isotherms were in reasonable agreement with those from CpHMD. MP1 adsorption was almost insensitive to solution pH. H-MP1 was much more sensitive to partitioning, at acidic pH, with an affinity ten times higher than in neutral ones.

Effects of N-terminus modifications on the conformation and permeation activities of the synthetic peptide L1A.

We investigate the effect of the N-terminus modification of the L1A, a synthetic octadecapeptide, on its helical content, affinity and lytic action in model membranes and on its hemolytic and antibacterial activities. L1A and its acetylated analog displayed a selective antibacterial activity to Gram-negative bacteria without being hemolytic. The covalently linked 2-aminobezoic acid to the N-terminus impaired the antibacterial efficacy and increased hemolysis. Despite their lower net charge (+2), N-terminus modifications resulted in enhanced affinity and improved lytic efficiency in anionic vesicles. The analogs also showed higher helical content and consequently higher amphipathicity in these vesicles. The conformational analysis by molecular dynamics simulations in 30 % of TFE/water showed that the hydrophobic faces of the peptides are in close contact with CF3 groups of TFE while the hydrophilic faces with water molecules. Due to the loss of the amino charge, the N-termini of the analogs are buried in TFE molecules. The analysis of the pair distribution functions, obtained for the center of mass of the charged groups, has evidenced that the state of the N-terminus has influenced the possibility of different ion-pairing. The higher complexity of the bacterial cells compared with anionic vesicles hampers to establish correlations structure-function for the analogs.

An integrated approach with experimental and computational tools outlining the cooperative binding between 2-phenylchromone and human serum albumin.

2-Phenylchromone (2PHE) is a flavone, found in cereals and herbs, indispensable in the human diet. Its chemical structure is the basis of all flavonoids present in black and green tea, soybean, red fruits and so on. Although offering such nutritional value, it still requires a molecular approach to understand its interactions with a specific target. The combination of experimental and computational techniques makes it possible to describe the interaction between 2PHE and human serum albumin (HSA). Fluorescence spectroscopy results show that the quenching mechanism is static, and thermodynamic analysis points to an entropically driven complex. The binding density function method provides information about a positive cooperative interaction, while drug displacement experiments indicate Sites 1 and 2 of HSA as the most probable binding sites. From the molecular dynamic study, it appears that the molecular docking is in agreement with experimental data and thus more realistic.

Structural modeling of high-affinity thyroid receptor-ligand complexes.

Understanding the molecular basis of the binding modes of natural and synthetic ligands to nuclear receptors is fundamental to our comprehension of the activation mechanism of this important class of hormone regulated transcription factors and to the development of new ligands. Thyroid hormone receptors (TR) are particularly important targets for pharmaceuticals development because TRs are associated with the regulation of metabolic rates, body weight, and circulating levels of cholesterol and triglycerides in humans. While several high-affinity ligands are known, structural information is only partially available. In this work we obtain structural models of several TR-ligand complexes with unknown structure by docking high affinity ligands to the receptors' ligand binding domain with subsequent relaxation by molecular dynamics simulations. The binding modes of these ligands are discussed providing novel insights into the development of TR ligands. The experimental binding free energies are reasonably well-reproduced from the proposed models using a simple linear interaction energy free-energy calculation scheme.