Computational Study of the pH-Dependent Ionic Environment around β-Lactoglobulin.

Ions are involved in multiple biological processes and may exist bound to biomolecules or may be associated with their surface. Although the presence of ions in nucleic acids has traditionally gained more interest, ion-protein interactions, often with a marked dependency on pH, are beginning to gather attention. Here we present a detailed analysis on the binding and distribution of ions around β-lactoglobulin using a constant-pH MD (CpHMD) method, at a pH range 3-8, and compare it with the more traditional Poisson-Boltzmann (PB) model and the existing experimental data.

FixBox: A General Algorithm to Fix Molecular Systems in Periodic Boxes.

Periodic boundary conditions (PBCs) are a standard feature of molecular simulations, and their mathematical and computational aspects are well-understood and relatively straightforward. However, they can in practice be a nuisance when simulating heterogeneous systems, especially when different types of molecules change their relative positions during the simulation. Although the translation required to fix a broken molecular complex of interest can in most cases be easily inferred by visual inspection, it typically depends on the type of system, its configuration, and the box geometry, making automated procedures problematic.

Approach to Study pH-Dependent Protein Association Using Constant-pH Molecular Dynamics: Application to the Dimerization of β-Lactoglobulin.

Protein-protein association is often mediated by electrostatic interactions and modulated by pH. However, experimental and computational studies have often overlooked the effect of association on the protonation state of the protein. In this work, we present a methodological approach based on constant-pH molecular dynamics (MD), which aims to provide a detailed description of a pH-dependent protein-protein association, and apply it to the dimerization of β-lactoglobulin (BLG).

Effect of pH on the influenza fusion peptide properties unveiled by constant-pH molecular dynamics simulations combined with experiment.

The influenza virus fusion process, whereby the virus fuses its envelope with the host endosome membrane to release the genetic material, takes place in the acidic late endosome environment. Acidification triggers a large conformational change in the fusion protein, hemagglutinin (HA), which enables the insertion of the N-terminal region of the HA2 subunit, known as the fusion peptide, into the membrane of the host endosome. However, the mechanism by which pH modulates the molecular properties of the fusion peptide remains unclear.

Molecular Anatomy of Plant Photoprotective Switches: The Sensitivity of PsbS to the Environment, Residue by Residue.

Under strong sunlight, plants avoid photooxidation by quenching the excess absorbed energy. Quenching is triggered by PsbS, a membrane protein that is activated and deactivated by the light-dependent pH changes in the thylakoid lumen. The mechanism of action of this protein is unknown, but it was suggested that several glutamates act as pH sensors.

Role of Counterions in Constant-pH Molecular Dynamics Simulations of PAMAM Dendrimers.

Electrostatic interactions play a pivotal role in the structure and mechanism of action of most biomolecules. There are several conceptually different methods to deal with electrostatics in molecular dynamics simulations. Ionic strength effects are usually introduced using such methodologies and can have a significant impact on the quality of the final conformation space obtained.

Effect of a pH Gradient on the Protonation States of Cytochrome c Oxidase: A Continuum Electrostatics Study.

Cytochrome c oxidase (CcO) couples the reduction of dioxygen to water with transmembrane proton pumping, which leads to the generation of an electrochemical gradient. In this study we analyze how one of the components of the electrochemical gradient, the difference in pH across the membrane, or ΔpH, influences the protonation states of residues in CcO. We modified our continuum electrostatics/Monte Carlo (CE/MC) method in order to include the ΔpH and applied it to the study of CcO, in what is, to our best knowledge, the first CE/MC study of CcO in the presence of a pH gradient.

Predicting protein partition coefficients in aqueous two phase system.

The present work aims to achieve an additional insight into the protein partitioning behavior in aqueous two phase systems (ATPSs), together with a study on the viability of a semi-empirical model based on continuum electrostatics to predict the protein partition characteristics. The partitioning behaviors of 14 globular proteins, with different properties, were explored in three polymer/polymer ATPSs. By the Collander equation, a linear correlation between protein partitioning coefficients in all systems was observed.

Structuring Peptide Dendrimers through pH Modulation and Substrate Binding.

Dendrimers are a family of ramified synthetic molecules. pH effects and electrostatic interactions are known to be crucial players to explain the conformational and functional behaviors observed in these systems. Nonetheless, to date, no computational study involving these systems has explicitly addressed the protonation equilibrium taking place at different pH values for dendrimers containing multiple ionizable sites.

Coupling between protonation and conformation in cytochrome c oxidase: Insights from constant-pH MD simulations.

Cytochrome c oxidases (CcOs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes reduce dioxygen (O(2)) to water and, simultaneously, generate a transmembrane electrochemical proton gradient. Despite their importance in the aerobic metabolism and the large amount of structural and biochemical data available for the A1-type CcO family, there is still no consensually accepted description of the molecular mechanisms operating in this protein.

Constant-pH MD Simulations Portray the Protonation and Structural Behavior of Four Decapeptides Designed to Coordinate Cu(2+).

The cyclic decapeptide C-Asp, containing one Asp residue and three His residues, was designed by Fragoso et al. (Chem. Eur.

The Effect of Membrane Environment on Surfactant Protein C Stability Studied by Constant-pH Molecular Dynamics.

Pulmonary surfactant protein C (SP-C) is a small peptide with two covalently linked fatty acyl chains that plays a crucial role in the formation and stabilization of the pulmonary surfactant reservoirs during the compression and expansion steps of the respiratory cycle. Although its function is known to be tightly related to its highly hydrophobic character and key interactions maintained with specific lipid components, much is left to understand about its molecular mechanism of action. Also, although it adopts a mainly helical structure while associated with the membrane, factors as pH variation and deacylation have been shown to affect its stability and function.

The anti-inflammatory action of the analgesic kyotorphin neuropeptide derivatives: insights of a lipid-mediated mechanism.

Recently, a designed class of efficient analgesic drugs derived from an endogenous neuropeptide, kyotorphin (KTP, Tyr-Arg) combining C-terminal amidation (KTP-NH2) and N-terminal conjugation to ibuprofen (Ib), IbKTP-NH2, was developed. The Ib moiety is an enhancer of KTP-NH2 analgesic action. In the present study, we have tested the hypothesis that KTP-NH2 is an enhancer of the Ib anti-inflammatory action.

Photonic activation of plasminogen induced by low dose UVB.

Activation of plasminogen to its active form plasmin is essential for several key mechanisms, including the dissolution of blood clots. Activation occurs naturally via enzymatic proteolysis. We report that activation can be achieved with 280 nm light.

Modeling the partitioning of amino acids in aqueous two phase systems.

A new model to obtain fast prediction of partition coefficients in polymer/polymer aqueous two phase systems (ATPSs) is presented, using amino acids as test systems. In particular, the partitioning behavior of eleven amino acids (glycine, alanine, leucine, phenylalanine, lysine, arginine, histidine, aspartic acid, glutamic acid, glutamine and serine) has been studied in 6 polymer/polymer ATPSs, formed by different pairs of nonionic polymers, including polyethylene glycol (PEG), Dextran, Ucon and Ficoll at 0.15M NaCl in 0.

Structural effects of pH and deacylation on surfactant protein C in an organic solvent mixture: a constant-pH MD study.

The pulmonary surfactant protein C (SP-C) is a small highly hydrophobic protein that adopts a mainly helical structure while associated with the membrane but misfolds into a β-rich metastable structure upon deacylation, membrane dissociation, and exposure to the neutral pH of the aqueous alveolar subphase, eventually leading to the formation of amyloid aggregates associated with pulmonary alveolar proteinosis. The present constant-pH MD study of the acylated and deacylated isoforms of SP-C in a chloroform/methanol/water mixture, often used to mimic the membrane environment, shows that the loss of the acyl groups has a structural destabilizing effect and that the increase of pH promotes intraprotein contacts which contribute to the loss of helical structure in solution. These contacts result from the poor solvation of charged groups by the solvent mixture, which exhibits a limited membrane-mimetic character.

Water encapsulation in a polyoxapolyaza macrobicyclic compound.

A new heteroditopic macrobicyclic compound (t(2)pN(5)O(3)) containing two separate polyoxa and polyaza compartments was synthesized in good yield through a [1 + 1] "tripod-tripod coupling" strategy. The X-ray crystal structure of H(3)t(2)pN(5)O(3)(3+) revealed the presence of one encapsulated water molecule accepting two hydrogen bonds from two protonated secondary amines and donating a hydrogen bond to one amino group. The acid-base behavior of the compound was studied by potentiometry at 298.

Reversibility of prion misfolding: insights from constant-pH molecular dynamics simulations.

The prion protein (PrP) is the cause of a group of diseases known as transmissible spongiform encephalopathies (TSEs). Creutzfeldt-Jakob disease and bovine spongiform encephalopathy are examples of TSEs. Although the normal form of PrP (PrP(C)) is monomeric and soluble, it can misfold into a pathogenic form (PrP(Sc)) that has a high content of β-structure and can aggregate forming amyloid fibrils.

Constant-pH molecular dynamics simulations reveal a β-rich form of the human prion protein.

The misfolding of the prion protein (PrP) into a pathogenic β-rich form (PrP(Sc)) has been suggested to occur in the endocytic pathway, triggered by low pH. In this work we performed several constant-pH molecular dynamics simulations of human PrP 90-231 in the pH range 2-7, totaling more than 2 μs. We observed a strong conformational pH dependence where on average the helix content decreased and the β content increased toward acidic pH.

Membrane-induced conformational changes of kyotorphin revealed by molecular dynamics simulations.

The analgesic dipeptide kyotorphin (l-Tyr-l-Arg) was studied in the two most relevant protonation states at physiological pH, both in water and in a membrane model, using molecular dynamics simulations. Kyotorphin is found to exhibit a remarkable conformational freedom even when strongly interacting with the bilayer. Nevertheless, we observe a strong decrease in the population of the tyrosine's chi(1) torsion angle around 60 degrees that could be correlated with the dipeptide biological function.

Conformational analysis in a multidimensional energy landscape: study of an arginylglutamate repeat.

The identification of the distinct conformation classes of a molecule is a common and often crucial step in establishing structure-function relationships. Many different methods have been suggested for that purpose which differ in their choice of a (dis)similarity measure and clustering algorithm. The present study discusses and analyzes these issues, proposing a method based on principal component analysis (PCA), which is applied to conformations obtained from molecular dynamics (MD) simulations of an arginylglutamate repeat.