Charged Amino Acid Residues Research Articles

Effect of Surface-Immobilized States of Antimicrobial Peptides on Their Ability to Disrupt Bacterial Cell Membrane Structure

Antimicrobial peptide (AMP) surfaces are widely used to inhibit biofilm formation and bacterial infection. However, endpoint-immobilized AMPs on surfaces are totally different from free-state AMPs due to the constraints of the surface. In this work, the interactions between AMPs and bacterial cell membranes were analyzed through coarse-grained molecular dynamics and all-atom molecular dynamics simulations. This AMP disrupted membrane structure by altering the thickness and curvature of the membrane. Furthermore, the effect of surface-immobilized states of AMPs on their ability to disrupt membrane structure was revealed. The immobilized AMPs in the freeze-N system could bind to the membrane and disrupt the membrane structure through electrostatic forces between positively charged N-terminal amino acid residues and the negatively charged membrane, while the immobilized AMPs in the freeze-C system were repelled. The results will aid in the rational design of new AMP surfaces with enhanced efficacy and stability.

Effect of peptide hydrophilicity on membrane curvature and permeation.

Using a well-developed reaction coordinate in umbrella sampling, we studied the single peptide permeation through a model cancerous cell membrane, varying the hydrophilicity and the charge of the peptides. Two peptides, melittin and pHD108, were studied. The permeation mechanism differs from a barrel-stave-like mechanism to toroidal pore and vesicle formation based on the number and the placement of the hydrophilic amino acids in the peptide. Membrane curvature changes dynamically as the permeation process occurs. In the case of vesicles, the peptide traverses along a smooth, homogenous pathway, whereas a rugged, steep pathway was found when lipid molecules did not line up along the wall of the membrane (barrel-stave-like mechanism). A mechanism similar to a toroidal pore consists of multiple minima. Higher free energy was found for the permeating terminal containing charged amino acid residues. Vesicle formation was found for pHD108 peptide N-terminal with a maximum membrane thinning effect of 54.4% with free energy cost of 8.20 ± 0.10kcal mol-1 and pore radius of 2.33 ± 0.07nm. Insights gained from this study can help to build synthetic peptides for drug delivery.

Deciphering the acidophilia and acid resistance in Acetilactobacillus jinshanensis dominating baijiu fermentation through multi-omics analysis

Lactic acid bacteria (LAB) are pivotal in constructing the intricate bio-catalytic networks underlying traditional fermented foods such as Baijiu. However, LAB and their metabolic mechanisms are partially understood in Moutai flavor Baijiu fermentation. Here, we found that Acetilactobacillus jinshanensis became the· dominant species with relative abundance reaching 92%, where the acid accumulated rapidly and peaked at almost 30 g/kg in Moutai flavor Baijiu. After separation, purification, and cultivation, A. jinshanensis exhibited pronounced acidophilia and higher acid resistance compared to other LAB. Further integrated multi-omics analysis revealed that fatty acid synthesis, cell membrane integrity, pHi and redox homeostasis maintenance, protein and amide syntheses were possibly crucial acid-resistant mechanisms in A. jinshanensis. Structural proteomics indicated that the surfaces of A. jinshanensis proteases contained more positively charged amino acid residues to maintain protein stability in acidic environments. The genes HSP20 and acpP were identified as acid-resistant genes for A. jinshanensis by heterologous expression analysis. These findings not only enhance our understanding of LAB in Baijiu, providing a scientific basis for acid regulation for production process, but also offer valuable insights for studying core species in other fermentation systems.

Statistical analysis of the unique characteristics of secondary structures in proteins

Protein folding is a complex process influenced by the primary sequence of amino acids. Early studies focused on understanding whether the specificity or the conservation of properties of amino acids was crucial for folding into secondary structures such as α-helices, β-sheets, turns, and coils. However, with the advent of artificial intelligence (AI) and machine learning (ML), the emphasis has shifted towards the precise nature and occurrence of specific amino acids. In our study, we analyzed a large set of proteins from diverse organisms to identify unique features of secondary structures, particularly in terms of the distribution of polar, non-polar, and charged amino acid residues. We found that α-helices tend to have a higher proportion of charged and non-polar groups compared to other secondary structures and that the presence of oppositely charged amino acid residues in helices stabilizes them, facilitating the formation of longer helices. These characteristics are distinct to α-helices. This study offers valuable insights for researchers in the field of protein design, enabling the de-novo creation of short helical peptides for a range of applications. We have also developed a web server for extensive analysis of proteins from different databases. The web server is housed at https://proseqanalyser.iitgn.ac.in/

Interactions of Potential Endocrine-Disrupting Chemicals with Whole Human Proteome Predicted by AlphaFold2 Using an In Silico Approach.

Binding with proteins is a critical molecular initiating event through which environmental pollutants exert toxic effects in humans. Previous studies have been limited by the availability of three-dimensional (3D) protein structures and have focused on only a small set of environmental contaminants. Using the highly accurate 3D protein structure predicted by AlphaFold2, this study explored over 60 million interactions obtained through molecular docking between 20,503 human proteins and 1251 potential endocrine-disrupting chemicals. A total of 66,613,773 docking results were obtained, 1.2% of which were considered to be high binding, as their docking scores were lower than -7. Monocyte to macrophage differentiation factor 2 (MMD2) was predicted to interact with the highest number of environmental pollutants (526), with polychlorinated biphenyls and polychlorinated dibenzofurans accounting for a significant proportion. Dimension reduction and clustering analysis revealed distinct protein profiles characterized by high binding affinities for perfluoroalkyl and polyfluoroalkyl substances (PFAS), phthalate-like chemicals, and other pollutants, consistent with their uniquely enriched pathways. Further structural analysis indicated that binding pockets with a high proportion of charged amino acid residues, relatively low α-helix content, and high β-sheet content were more likely to bind to PFAS than others. This study provides insights into the toxicity pathways of various pollutants impacting human health and offers novel perspectives for the establishment and expansion of adverse outcome pathway-based models.

Enhancing the Intrinsic Antiplasmodial Activity and Improving the Stability and Selectivity of a Tunable Peptide Scaffold Derived from Human Platelet Factor 4.

The control of malaria, a disease caused by Plasmodium parasites that kills over half a million people every year, is threatened by the continual emergence and spread of drug resistance. Therefore, new molecules with different mechanisms of action are needed in the antimalarial drug development pipeline. Peptides developed from host defense molecules are gaining traction as anti-infectives due to theood of inducing drug resistance. Human platelet factor 4 (PF4) has intrinsic activity against P. falciparum, and a macrocyclic helix-loop-helix peptide derived from its active domain recapitulates this activity. In this study, we used a stepwise approach to optimize first-generation PF4-derived internalization peptides (PDIPs) by producing analogues with substitutions to charged and hydrophobic amino acid residues or with modifications to terminal residues including backbone cyclization. We evaluated the in vitro activity of PDIP analogues against P. falciparum compared to their overall helical structure, resistance to breakdown by serum proteases, selective binding to negatively charged membranes, and hemolytic activity. Next, we combined antiplasmodial potency-enhancing substitutions that retained favorable membrane and cell-selective properties onto the most stable scaffold to produce a backbone cyclic PDIP analogue with four-fold improved activity against P. falciparum compared to first-generation peptides. These studies demonstrate the ability to modify PDIP to select for and combine desirable properties and further validate the suitability of this unique peptide scaffold for developing a new molecule class that is distinct from existing antimalarial drugs.

Crystal structure and characterization of monascin from the extracts of Monascus purpureus-fermented rice.

We present a novel solid form of monascin, an azaphilonoid derivative extracted from Monascus purpureus-fermented rice. The crystal structure, C21H26O5, was characterized by single-crystal X-ray diffraction and belongs to the orthorhombic space group P212121. To gain insight into the electronic properties of the short contacts in the crystalline state of monascin, we utilized the Experimental Library of Multipolar Atom Model 2 (ELMAM2) database to transfer the electron density of monascin in its crystalline state. Hirshfeld surface analysis, fingerprint analysis, electronic properties and energetic characterization reveal that intermolecular C-H...O hydrogen bonds play a crucial role in the noncovalent bonding interactions by connecting molecules into two- and three-dimensional networks. The molecular electrostatic potential (MEP) map of the monascin molecule demonstrates that negatively charged regions located at four O atoms are favoured binding sites for more positively charged amino acid residues during molecular recognition. In addition, powder X-ray diffraction confirms that no transformation occurs during the crystallization of monascin.

Combined clinical, structural, and cellular studies discriminate pathogenic and benign TRPV4 variants.

Dominant mutations in the calcium-permeable ion channel TRPV4 (transient receptor potential vanilloid 4) cause diverse and largely distinct channelopathies, including inherited forms of neuromuscular disease, skeletal dysplasias, and arthropathy. Pathogenic TRPV4 mutations cause gain of ion channel function and toxicity that can be rescued by small molecule TRPV4 antagonists in cellular and animal models, suggesting that TRPV4 antagonism could be therapeutic for patients. Numerous variants in TRPV4 have been detected with targeted and whole exome/genome sequencing, but for the vast majority, their pathogenicity remains unclear. Here, we used a combination of clinical information and experimental structure-function analyses to evaluate 30 TRPV4 variants across various functional protein domains. We report clinical features of seven patients with TRPV4 variants of unknown significance and provide extensive functional characterization of these and an additional 17 variants, including structural position, ion channel function, subcellular localization, expression level, cytotoxicity, and protein-protein interactions. We find that gain-of-function mutations within the TRPV4 intracellular ankyrin repeat domain target charged amino acid residues important for RhoA interaction, whereas ankyrin repeat domain residues outside of the RhoA interface have normal or reduced ion channel activity. We further identify a cluster of gain-of-function variants within the intracellular intrinsically disordered region that may cause toxicity via altered interactions with membrane lipids. In contrast, assessed variants in the transmembrane domain and other regions of the intrinsically disordered region do not cause gain of function and are likely benign. Clinical features associated with gain of function and cytotoxicity include congenital onset of disease, vocal cord weakness, and motor predominant disease, whereas patients with likely benign variants often demonstrated late-onset and sensory-predominant disease. These results provide a framework for assessing additional TRPV4 variants with respect to likely pathogenicity, which will yield critical information to inform patient selection for future clinical trials for TRPV4 channelopathies.

Miniprotein engineering for inhibition of PD-1/PD-L1 interaction.

Miniproteins constitute an excellent basis for the development of structurally demanding functional molecules. The engrailed homeodomain, a three-helix-containing miniprotein, was applied as a scaffold for constructing programmed cell death protein 1/programmed death-ligand 1 (PD-1/PD-L1) interaction inhibitors. PD-L1 binders were initially designed using the computer-aided approach and subsequently optimized iteratively. The conformational stability was assessed for each obtained miniprotein using circular dichroism spectroscopy, indicating that numerous mutations could be introduced. The formation of a sizable hydrophobic surface at the inhibitor that fits the molecular target imposed the necessity for the incorporation of additional charged amino acid residues to retain its appropriate solubility. Finally, the miniprotein effectively binding to PD-L1 (KD = 51.4 nM) that inhibits PD-1/PD-L1 interaction in cell-based studies with EC50 = 3.9 μM, was discovered.

Insights into the role of glycerophospholipids on the iron export function of SLC40A1 and the molecular mechanisms of ferroportin disease.

SLC40A1 is the sole iron export protein reported in mammals. In humans, its dysfunction is responsible for ferroportin disease, an inborn error of iron metabolism transmitted as an autosomal dominant trait and observed in different ethnic groups. As a member of the major facilitator superfamily, SLC40A1 requires a series of conformational changes to enable iron translocation across the plasma membrane. The influence of lipids on protein stability and its conformational changes has been little investigated to date. Here, we combine molecular dynamics simulations of SLC40A1 embedded in membrane bilayers with experimental alanine scanning mutagenesis to analyze the specific role of glycerophospholipids. We identify four basic residues (Lys90, Arg365, Lys366, and Arg371) that are located at the membrane-cytosol interface and consistently interact with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) molecules. These residues surround a network of salt bridges and hydrogens bonds that play a critical role in stabilizing SLC40A1 in its basal outward-facing conformation. More deeply embedded in the plasma membrane, we identify Arg179 as a charged amino acid residue also tightly interacting with lipid polar heads. This results in a local deformation of the lipid bilayer. Interestingly, Arg179 is adjacent to Arg178, which forms a functionally important salt-bridge with Asp473 and is a recurrently associated with ferroportin disease when mutated to glutamine. We demonstrate that the two p.Arg178Gln and p.Arg179Thr missense variants have similar functional behaviors. These observations provide insights into the role of phospholipids in the formation/disruption of the SLC40A1 inner gate, and give a better understanding of the diversity of molecular mechanisms of ferroportin disease.

The road to evolution of ProTx2: how to be a subtype-specific inhibition of human Nav1.7.

The human voltage-gated sodium channel Nav1.7 is a widely proven target for analgesic drug studies. ProTx2, a 30-residue polypeptide from Peruvian green tarantula venom, shows high specificity to activity against human Nav1.7, suggesting its potential to become a non-addictive analgesic. However, its high sensitivity to human Nav1.4 raises concerns about muscle side effects. Here, we engineered three mutants (R13A, R13D, and K27Y) of ProTx2 to evaluate their pharmacological activities toward Nav1.7 and Nav1.4. It is demonstrated that the mutant R13D maintained the analgesic effect in mice while dramatically reducing its muscle toxicity compared with ProTx2. The main reason is the formation of a strong electrostatic interaction between R13D and the negatively charged amino acid residues in DII/S3-S4 of Nav1.7, which is absent in Nav1.4. This study advances our understanding and insights on peptide toxins, paving the way for safer, effective non-addictive analgesic development.

Effect of caffeine on the aggregation of amyloid-β-A 3D RISM study.

Alzheimer's disease is a detrimental neurological disorder caused by the formation of amyloid fibrils due to the aggregation of amyloid-β peptide. The primary therapeutic approaches for treating Alzheimer's disease are targeted to prevent this amyloid fibril formation using potential inhibitor molecules. The discovery of such inhibitor molecules poses a formidable challenge to the design of anti-amyloid drugs. This study investigates the effect of caffeine on dimer formation of the full-length amyloid-β using a combined approach of all-atom, explicit water molecular dynamics simulations and the three-dimensional reference interaction site model theory. The change in the hydration free energy of amyloid-β dimer, with and without the inhibitor molecules, is calculated with respect to the monomeric amyloid-β, where the hydration free energy is decomposed into energetic and entropic components, respectively. Dimerization is accompanied by a positive change in the partial molar volume. Dimer formation is spontaneous, which implies a decrease in the hydration free energy. However, a reverse trend is observed for the dimer with inhibitor molecules. It is observed that the negatively charged residues primarily contribute for the formation of the amyloid-β dimer. A residue-wise decomposition reveals that hydration/dehydration of the side-chain atoms of the charged amino acid residues primarily contribute to dimerization.

AGRAMP: machine learning models for predicting antimicrobial peptides against phytopathogenic bacteria.

Antimicrobial peptides (AMPs) are promising alternatives to traditional antibiotics for combating plant pathogenic bacteria in agriculture and the environment. However, identifying potent AMPs through laborious experimental assays is resource-intensive and time-consuming. To address these limitations, this study presents a bioinformatics approach utilizing machine learning models for predicting and selecting AMPs active against plant pathogenic bacteria. N-gram representations of peptide sequences with 3-letter and 9-letter reduced amino acid alphabets were used to capture the sequence patterns and motifs that contribute to the antimicrobial activity of AMPs. A 5-fold cross-validation technique was used to train the machine learning models and to evaluate their predictive accuracy and robustness. The models were applied to predict putative AMPs encoded by intergenic regions and small open reading frames (ORFs) of the citrus genome. Approximately 7% of the 10,000-peptide dataset from the intergenic region and 7% of the 685,924-peptide dataset from the whole genome were predicted as probable AMPs. The prediction accuracy of the reported models range from 0.72 to 0.91. A subset of the predicted AMPs was selected for experimental test against Spiroplasma citri, the causative agent of citrus stubborn disease. The experimental results confirm the antimicrobial activity of the selected AMPs against the target bacterium, demonstrating the predictive capability of the machine learning models. Hydrophobic amino acid residues and positively charged amino acid residues are among the key features in predicting AMPs by the Random Forest Algorithm. Aggregation propensity appears to be correlated with the effectiveness of the AMPs. The described models would contribute to the development of effective AMP-based strategies for plant disease management in agricultural and environmental settings. To facilitate broader accessibility, our model is publicly available on the AGRAMP (Agricultural Ngrams Antimicrobial Peptides) server.

An extensive ion-pair/hydrogen-bond network contributes to the thermostability of the MutL ATPase domain from Aquifex aeolicus.

Proteins from hyperthermophiles often contain a large number of ionic interactions. Close examination of the previously determined crystal structure of the ATPase domain of MutL from a hyperthermophile, Aquifex aeolicus, revealed that the domain contains a continuous ion-pair/hydrogen-bond network consisting of 11 charged amino acid residues on a β-sheet. Mutations were introduced to disrupt the network, showing that the more extensively the network was disrupted, the greater the thermostability of the protein was decreased. Based on urea denaturation analysis, a thermodynamic parameter, energy for the conformational stability, was evaluated, which indicated that amino acid residues in the network contributed additively to the protein stability. A continuous network rather than a cluster of isolated interactions would pay less entropic penalty upon fixing the side chains to make the same number of ion pairs/hydrogen bonds, which might contribute more favorably to the structural formation of thermostable proteins.

Changes in total charge on spike protein of SARS-CoV-2 in emerging lineages.

Charged amino acid residues on the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been shown to influence its binding to different cell surface receptors, its non-specific electrostatic interactions with the environment, and its structural stability and conformation. It is therefore important to obtain a good understanding of amino acid mutations that affect the total charge on the spike protein which have arisen across different SARS-CoV-2 lineages during the course of the virus' evolution. We analyse the change in the number of ionizable amino acids and the corresponding total charge on the spike proteins of almost 2200 SARS-CoV-2 lineages that have emerged over the span of the pandemic. Our results show that the previously observed trend toward an increase in the positive charge on the spike protein of SARS-CoV-2 variants of concern has essentially stopped with the emergence of the early omicron variants. Furthermore, recently emerged lineages show a greater diversity in terms of their composition of ionizable amino acids. We also demonstrate that the patterns of change in the number of ionizable amino acids on the spike protein are characteristic of related lineages within the broader clade division of the SARS-CoV-2 phylogenetic tree. Due to the ubiquity of electrostatic interactions in the biological environment, our findings are relevant for a broad range of studies dealing with the structural stability of SARS-CoV-2 and its interactions with the environment. The data underlying the article are available in the Supplementary material.

Conformational control of antimicrobial peptide amphiphilicity: consequences for boosting membrane interactions and antimicrobial effects of photocatalytic TiO2 nanoparticles.

This study reports on the effects of conformationally controlled amphiphilicity of antimicrobial peptides (AMPs) on their ability to coat TiO2 nanoparticles (NPs) and boost the photocatalytic antimicrobial effects of such NPs. For this, TiO2 NPs were combined with AMP EFK17 (EFKRIVQRIKDFLRNLV), displaying a disordered conformation in aqueous solution but helix formation on interaction with bacterial membranes. The membrane-bound helix is amphiphilic, with all polar and charged amino acid residues located at one side and all non-polar and hydrophobic residues on the other. In contrast, the d-enantiomer variant EFK17-d (E(dF)KR(dI)VQR(dI)KD(dF)LRNLV) is unable to form the amphiphilic helix on bacterial membrane interaction, whereas the W-residues in EFK17-W (EWKRWVQRWKDFLRNLV) boost hydrophobic interactions of the amphiphilic helix. Circular dichroism results showed the effects displayed for the free peptide, to also be present for peptide-coated TiO2 NPs, causing peptide binding to decrease in the order EFK17-W > EFK17 > EFK17-d. Notably, the formation of reactive oxygen species (ROS) by the TiO2 NPs was essentially unaffected by the presence of peptide coating, for all the peptides investigated, and the coatings stabilized over hours of UV exposure. Photocatalytic membrane degradation from TiO2 NPs coated with EFK17-W and EFK17 was promoted for bacteria-like model bilayers containing anionic phosphatidylglycerol but suppressed in mammalian-like bilayers formed by zwitterionic phosphatidylcholine and cholesterol. Structural aspects of these effects were further investigated by neutron reflectometry with clear variations observed between the bacteria- and mammalian-like model bilayers for the three peptides. Mirroring these results in bacteria-like model membranes, combining TiO2 NPs with EFK17-W and EFK17, but not with non-adsorbing EFK17-d, resulted in boosted antimicrobial effects of the resulting cationic composite NPs already in darkness, effects enhanced further on UV illumination.

Elucidating the microscopic origin of observed anti-correlation between the protein-protein and protein-water interaction energies

Water as a solvent plays an important role in dictating the structure and dynamics of the large biomolecules like proteins. It has been proposed that the protein dynamics and energetics are closely coupled to the water molecules in the hydration shell. In particular, it has been shown earlier that the fluctuations in intra-protein (protein-protein) interaction energy is strongly anti-correlated with the protein-water interaction energy (solvation energy). Moreover, the dynamical signatures of water gets reflected into the protein dynamics in terms of a dominant power law behavior. In this work, we have carried out systematic molecular dynamics (MD) simulations of three different protein systems to understand the key physical factors that contribute to such anti-correlations. First, we compare the dynamics in globular folded proteins like Lysozyme with intrinsically disordered peptides (IDPs) like Amyloid beta to show that enhanced conformational fluctuations lead to significantly stronger anti-correlation. Trp-cage miniprotein shows a similar behavior distributed over two sub-ensembles of folded and unfolded states. We also investigate what type of interactions and molecular processes are responsible for the observed anti-correlation. In particular, we observe that the extent of correlation decreases significantly in the order of charged, polar and non-polar amino acid residues. Charged residues that form salt bridge interactions show a dominant anti-correlation. Surprisingly, a few non-polar (hydrophobic and buried) residues demonstrate significant anti-correlation as well. We attribute this to the structural fluctuations that lead to an equilibrium between buried and solvent exposed states of these residues, where the protein-protein and protein-water interactions compensate each other. Our studies indicate that the molecular origin of the observed anti-correlation can be vary depending on the type of the interaction. The overall coupling/correlation is dominated by the electrostatic interactions and charged residues due to a general dielectric screening effect, whereas a few buried hydrophobic residues may also show such coupling due to structural dynamics between buried and exposed states.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases.

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Effects of Hydrogen Peroxide on the Hydrogen Bonding Structure and Dynamics of Water and Its Influence on the Aqueous Solvation of the Insulin Monomer.

The hydrogen bond structure and dynamics of water and hydrogen peroxide (H2O2) in their binary mixtures have been studied at 298 K by classical molecular dynamics simulations. Twelve different concentrations of aqueous-H2O2 solutions are considered for this study. We have analyzed the interactions between water and H2O2 by site-site pair correlation functions and observed that the probability of formation of OW···HP hydrogen bonds are higher compared to OP···HW. The second solvation shell of water is strongly affected by increasing H2O2 concentrations (XP > 0.50), which signifies the destruction of the tetrahedral network structure of water. The translational and rotational dynamics of water and H2O2 do not significantly change up to 25% of H2O2 in aqueous mixtures. The hydrogen bond lifetime of water-water, water-H2O2, and H2O2-H2O2 in the aqueous-H2O2 solutions shows a very minimal change with increasing H2O2 concentrations. In addition to this, we also investigated the effect of H2O2 on the insulin monomer and observed that higher concentrations of H2O2 (XP = 0.10) change the secondary structure. The influence of H2O2 is more on chain-B than that on chain-A in the insulin monomer. The H2O2 occupancy at the protein surface is higher for negatively charged (GLU) and polar (ASN and THR) amino acid residues compared with that for positively charged and neutral residues in the solutions.

Interaction of Glycated Albumin with Receptor for Glycation End Products According to Molecular Modeling Data

In diabetes mellitus (DM) patients, the accumulation of advanced glycation end products (AGE) leads to inflammation and oxidative stress through the activation of specific receptors for AGE (RAGE). Glycated albumin (gHSA) makes a significant contribution to the overall level of AGE in human body and, as a result, to the pathogenesis of DM and concomitant diseases. The mechanism of interaction of gHSA with RAGE is practically not studied. The purpose of the present paper is to study the binding of gHSA to RAGE using molecular modeling methods, to find the main sites of interaction and structural features of glycation sites that determine the efficiency of complex formation with RAGE. Ten gHSA models were constructed using molecular docking and molecular dynamics (MD) methods; each model corresponded to one modified lysine residue (carboxymethyl-lysine): Lys64, Lys73, Lys137, Lys233, Lys262, Lys317, Lys378, Lys525, Lys573, Lys574. Complexes of gHSA with the V-domain of RAGE were constructed using the macromolecular docking method, and their stability was studied using MD simulation. In the constructed gHSA models, the carboxyl groups of glycated Lys317 and Lys525 form intramolecular salt bridges with surrounding amino acids; in other cases, the carboxyl groups of the modified lysines are free to interact with positively charged amino acid residues on the RAGE surface. According to the macromolecular docking data and subsequent MD simulation, the complex of RAGE with gHSA glycated at Lys233 is most effective in terms of strength and specificity. Specific RAGE complexes with gHSA glycated at Lys317 and Lys574 are not formed. The obtained data on the interaction of gHSA with RAGE will help to understand the role of albumin in the pathophysiology of DM and advance towards the prevention and development of effective therapy for this disease.