Hydration and ion-binding of glycine betaine: How they may be involved into protection of proteins under abiotic stresses

Enhancing osmotic stress tolerance of cell mimetics by modulating lipid bilayer

Seeking effective ways to maintain cellular homeostasis is crucial to the survival of organisms when they encounter osmotic stress. Glycine betaine (GB) is a widely generated natural osmolyte, but its endogenous production and action are limited. Herein, a kind of nonionic surfactant dodecyl-β-d-glucopyranoside (DG) and a common polymer polyethylene glycol (PEG) are proven to have the ability to enhance the osmotic stress (induced by sugar concentration changes) tolerance of cell and organism models, those are giant unilamellar vesicles (GUVs) and gram-negative Escherichia coli. DG or PEG only induces small size decrease and certain shape change of GUVs. Importantly, DG or PEG at the concentration 100 times lower than that of GB effectively increases the survival rate of bacteria under both hypoosmotic and hyperosmotic conditions. This intriguing result is attributed to the insertion of DG or adsorption of PEG in the lipid bilayer membrane, leading to enhanced membrane permeability. These exogenous substances can replace GB to facilely and highly efficiently augment adaptation of organisms to osmotic stress.

The Three-Dimensional Reference Interaction Site Model Approach as a Promising Tool for Studying Hydrated Viruses and Their Complexes with Ligands

Viruses are the most numerous biological form living in any ecosystem. Viral diseases affect not only people but also representatives of fauna and flora. The latest pandemic has shown how important it is for the scientific community to respond quickly to the challenge, including critically assessing the viral threat and developing appropriate measures to counter this threat. Scientists around the world are making enormous efforts to solve these problems. In silico methods, which allow quite rapid obtention of, in many cases, accurate information in this field, are effective tools for the description of various aspects of virus activity, including virus-host cell interactions, and, thus, can provide a molecular insight into the mechanism of virus functioning. The three-dimensional reference interaction site model (3D-RISM) seems to be one of the most effective and inexpensive methods to compute hydrated viruses, since the method allows us to provide efficient calculations of hydrated viruses, remaining all molecular details of the liquid environment and virus structure. The pandemic challenge has resulted in a fast increase in the number of 3D-RISM calculations devoted to hydrated viruses. To provide readers with a summary of this literature, we present a systematic overview of the 3D-RISM calculations, covering the period since 2010. We discuss various biophysical aspects of the 3D-RISM results and demonstrate capabilities, limitations, achievements, and prospects of the method using examples of viruses such as influenza, hepatitis, and SARS-CoV-2 viruses.

Specific Proton-Donor Properties of Glycine Betaine. Metric Parameters and Enthalpy of Noncovalent Interactions in its Dimer, Water Complexes and Crystalline Hydrate

Trimethylglycine (glycine betaine, GB) is an important organic osmolyte that accumulates in various plant species in response to environmental stresses and has significant potential as a bioactive agent with low environmental impact. It is assumed that the hydration of GB is playing an important role in the protective mechanism. The hydration and aggregation properties of GB have not yet been studied in detail at the atomistic level. In this work, noncovalent interactions in the GB dimer and its complexes with water and crystalline monohydrate are studied. Depending on the object, periodic and non-periodic DFT calculations are used. Particular attention is paid to the metric parameters and enthalpies of intermolecular hydrogen bonds. The identification of noncovalent interactions is carried out by means of the Bader analysis of periodic or non-periodic electron density. The enthalpy of hydrogen bonds is estimated using the Rosenberg formula (PCCP 2 (2000) 2699). The specific proton donor properties of glycine betaine are due to its ability to form intermolecular C-H∙∙∙O bonds with the oxygen atom of a water molecule or the carboxylate group of a neighboring GB. The enthalpy of these bonds can be significantly greater than 10 kJ/mol. The water molecule that forms a hydrogen bond with the carboxylate group of GB also interacts with its CH groups through lone pairs of electrons. The C-H∙∙∙O bonds contribute up to 40% of the total entropy of the GB-water interaction, which is about 45 kJ/mol. The possibility of identifying C-H∙∙∙O bonds by the proton nuclear magnetic resonance method is discussed.

Protein 3D Hydration: A Case of Bovine Pancreatic Trypsin Inhibitor

Characterization of the hydrated state of a protein is crucial for understanding its structural stability and function. In the present study, we have investigated the 3D hydration structure of the protein BPTI (bovine pancreatic trypsin inhibitor) by molecular dynamics (MD) and the integral equation method in the three-dimensional reference interaction site model (3D-RISM) approach. Both methods have found a well-defined hydration layer around the protein and revealed the localization of BPTI buried water molecules corresponding to the X-ray crystallography data. Moreover, under 3D-RISM calculations, the obtained positions of waters bound firmly to the BPTI sites are in reasonable agreement with the experimental results mentioned above for the BPTI crystal form. The analysis of the 3D hydration structure (thickness of hydration shell and hydration numbers) was performed for the entire protein and its polar and non-polar parts using various cut-off distances taken from the literature as well as by a straightforward procedure proposed here for determining the thickness of the hydration layer. Using the thickness of the hydration shell from this procedure allows for calculating the total hydration number of biomolecules properly under both methods. Following this approach, we have obtained the thickness of the BPTI hydration layer of 3.6 Å with 369 water molecules in the case of MD simulation and 3.9 Å with 333 water molecules in the case of the 3D-RISM approach. The above procedure was also applied for a more detailed description of the BPTI hydration structure near the polar charged and uncharged radicals as well as non-polar radicals. The results presented for the BPTI as an example bring new knowledge to the understanding of protein hydration.

Metabolic changes in antennal glands of Caribbean spiny lobsters Panulirus argus infected by Panulirus argus virus 1 (PaV1)

Panulirus argus virus 1 (PaV1) (Family Mininucleoviridae) causes chronic and systemic infection in wild juvenile spiny lobsters Panulirus argus (Latreille, 1804), ending in death by starvation and metabolic wasting. In marine decapods, the antennal gland is involved in osmoregulation and excretion. In this compact organ, fluid is filtered from the hemolymph, and ions are reabsorbed to produce a hypotonic urine. Although PaV1 is released with the urine in infected individuals, little is known regarding the metabolic effect of PaV1 in the antennal gland. The objective of this study was to perform a comparative evaluation of the metabolic profile of the antennal gland of clinically PaV1-infected lobsters versus those with no clinical signs of infection, using proton nuclear magnetic resonance analysis. Overall, 48 compounds were identified, and the most represented metabolites were those involved in carbohydrate, amino acid, energy, and nucleotide metabolism. Most of the metabolites that were down-regulated in the infected group were essential and non-essential amino acids. Some metabolites involved in the urea cycle and carbohydrate metabolism were also altered. This study represents a first approach to the metabolic evaluation of the antennal gland. We broadly discuss alterations in the content of several proteinogenic and non-proteinogenic amino acids and other key metabolites involved in energetic and nucleotide metabolism.

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

Ion Pairing of the Neurotransmitters Acetylcholine and Glutamate in Aqueous Solutions

Neurotransmitters (NTs) play an important role in neural communication, regulating a variety of functions such as motivation, learning, memory, and muscle contraction. Their intermolecular interactions in biological media are an important factor affecting their biological activity. However, the available information on the features of these interactions is scarce and contradictory, especially, in an estimation of possible ion binding. In this paper, we present the results of a study for two well-known NTs, acetylcholine (ACh) and glutamate (Glu), with relation to the NT-inorganic ion and the NT-NT binding in a water environment. The features of NT pairing are investigated in aqueous AChCl and NaGlu solutions over a wide concentration range using the integral equation method in 1D- and 3D- reference interaction site model (RISM) approaches. The data for ACh are given for its two bioactive TG (trans, gauche) and TT (trans, trans) conformers. As was found, for both NTs, the results indicate the NT-inorganic counterion contact pair to be the predominant associate type in the concentrated solutions. In this case, the counterions occupy the vacated "water" space in the hydration shell of the onium moiety (ACh) or carboxylate groups (Glu). For ACh, the "unfolded" TT conformer demonstrates a slightly greater possibility for counterion pairing in comparison with the "folded" TG conformer. For Glu, the probability of its binding with a counterion is slightly stronger for the "side-chain" carboxylate group than for the "backbone" group. The obtained results also revealed an insignificant probability of Glu^--Glu^- pairing. Namely, the RISM data indicate Glu^--Glu^- binding by NH₃⁺-COO^- interactions. A link between the ion binding of NTs and their biological activity is discussed. This contribution adds new knowledge to our understanding of the interactions between the NTs and their molecular environment, providing further insights into the behavior of these compounds in biological media.

Hydration of the Zwitterionic and Protonated Forms of Glycine Betaine Probed by Soft X-ray Emission Spectroscopy Coupled with Chemometrics

Soft X-ray absorption and emission spectra of glycine betaine (GB) have been measured at the O K-edge in neutral and strongly acidic solutions. The absorption spectra of the neutral solutions have a resonance peak at 532.6 eV, assigned to the transition to the π* orbital, whereas in the acidic solutions, the peak is shifted by -0.3 eV. The emission spectra taken as a function of the GB concentration have been analyzed by means of a modified classical least-squares regression method to obtain the hydration number of the solute. The analysis is successful when the emission spectra have been acquired at the energy of a slightly detuned resonance, giving 28 and 24 as the minimum values for the zwitterionic and protonated GB, respectively. The number of 28 accords with the reported values for the number of water molecules in the first hydration layer of the zwitterion and is greater than that obtained by other experimental techniques. The obtained numbers are used to discuss the hydration structure of GB with the aid of ab initio molecular orbital calculations. The hydration structure of the protonated form of GB is explored for the first time.