Hydration enthalpy characteristics of amino acids in solutions

The enthalpies of solvation of 17 amino acids were evaluated by using the sublimation enthalpies of amino acids and the standard enthalpies of their solution in water. An equation was derived, which relates the volume-specific enthalpy of sublimation (deltaH(subl)/V(w)) to the sum of the common bond lengths in molecules (sigman(i)l(i)) of substances examined. The results obtained are interpreted in terms of the effect of hydrophobic and hydrophilic side chain on the interactions between the zwitterions of amino acids and water molecules.     

Effect of amino acids and low molecular weight peptides on the dynamics of water

The effect of amino acids and low-molecular-weight peptides on the dynamics of water was studied. Water medium together with molecules of dipeptides and amino acids dissolved in it is considered as a complex of interacting anharmonic oscillators. It was shown that the temperature behaviour of this system is determined by nonlinear resonances, which give rise to both the phenomenon of self-synchronization in the whole system or its part and to the phenomena of phase instability and coherence decay, depending on the store of oscillatory energy. Dissolved molecules are also involved in these oscillations, and if the frequencies and amplitudes of their oscillations are within the range in which nonlinear resonance occurs, they can affect the movement of the whole system.     

Probing the hirudin-thrombin interaction by incorporation of noncoded amino acids and molecular dynamics simulation

Thrombin is a primary target for the development of novel anticoagulants, since it plays two important and opposite roles in hemostasis: procoagulant and anticoagulant. All thrombin functions are influenced by Na+ binding, which triggers the transition of this enzyme from an anticoagulant (slow) form to a procoagulant (fast) form. In previous studies, we have conveniently produced by chemical synthesis analogues of the N-terminal fragment 1-47 of hirudin HM2 containing noncoded amino acids and displaying up to approximately 2700-fold more potent antithrombin activity, comparable to that of full-length hirudin. In the work presented here, we have exploited the versatility of chemical synthesis to probe the structural and energetic properties of the S3 site of thrombin through perturbations introduced in the structure of hirudin fragment 1-47. In particular, we have investigated the effects of systematic replacement of Tyr3 with noncoded amino acids retaining the aromatic nucleus of Tyr, as well as similar hydrophobic and steric properties, but possessing different electronic (e.g., p-fluoro-, p-iodo-, or p-nitro-Phe), charge (p-aminomethyl-Phe), or conformational (homo-Phe) properties. Our results indicate that the affinity of fragment 1-47 for thrombin is proportional to the desolvation free energy change upon complex formation, and is inversely related to the electric dipole moment of the amino acid side chain at position 3 of hirudin. In this study, we have also identified the key features that are responsible for the preferential binding of hirudin to the procoagulant (fast) form of thrombin. Strikingly, shaving at position 3, by Tyr --> Ala exchange, abolishes the differences in the affinity for thrombin allosteric forms, whereas a bulkier side chain (e.g., beta-naphthylalanine) improves binding preferentially to the fast form. These results provide strong, albeit indirect, evidence that the procoagulant (fast) form of thrombin is in a more open and accessible conformation with respect to the less forgiving structure it acquires in the slow form. This view is also supported by the results of molecular dynamics simulations conducted for 18 ns on free thrombin in full explicit water, showing that after approximately 5 ns thrombin undergoes a significant conformational transition, from a more open conformation (which we propose can be related to the fast form) to a more compact and closed one (which we propose can be related to the slow form). This transition mainly involves the Trp148 and Trp60D loop, the S3 site, and the fibrinogen binding site, whereas the S1 site, the Na+-binding site, and the catalytic pocket remain essentially unchanged. In particular, our data indicate that the S3 site of the enzyme is less accessible to water in the putative slow form. This structural picture provides a reasonable molecular explanation for the fact that physiological substrates related to the procoagulant activity of thrombin (fibrinogen, thrombin receptor 1, and factor XIII) orient a bulky side chain into the S3 site of the enzyme. Taken together, our results can have important implications for the design of novel thrombin inhibitors, of practical utility in the treatment of coagulative disorders.     

Self-assembly and interactions of short antimicrobial cationic lipopeptides with membrane lipids: ITC,FTIR and molecular dynamics studies

In this work, the self-organization and the behavior of the surfactant-like peptides in the presence of biological membrane models were studied. The studies were focused on synthetic palmitic acid-containing lipopeptides, C₁₆–KK–NH₂ (I), C₁₆–KGK–NH₂ (II) and C₁₆–KKKK–NH₂ (III). The self-assembly was explored by molecular dynamics simulations using a coarse-grained force field. The critical micellar concentration was estimated by the surface tension measurements. The thermodynamics of the peptides binding to the anionic and zwitterionic lipids were established using isothermal titration calorimetry (ITC). The influence of the peptides on the lipid acyl chain ordering was determined using FTIR spectroscopy. The compounds studied show surface-active properties with a distinct CMC over the millimolar range. An increase in the steric and electrostatic repulsion between polar head groups shifts the CMC toward higher values and reduces the aggregation number. An analysis of the peptide–membrane binding revealed a unique interplay between the initial electrostatic and the subsequent hydrophobic interactions enabling the lipopeptides to interact with the lipid bilayer. In the case of C₁₆–KKKK–NH₂ (III), compensation of the electrostatic and hydrophobic interactions upon binding to the anionic membrane has been suggested and consequently no overall binding effects were noticed in ITC thermograms and FTIR spectra.     

FTIR,FT-Raman spectra and quantum chemical studies of nebivolol

Fourier transform infrared and Raman spectra of nebivolol have been recorded. The structure, conformational stability, geometry optimisation, and vibrational wave numbers have been investigated. Satisfactory vibrational assignments were made for the stable conformer of the molecule using Restricted Hartree–Fock (RHF) and density functional theory (DFT) calculation (B3LYP) with the 6-31G(d,p) basis set. Comparison of the observed fundamental vibrational wave numbers of the molecule and calculated results by RHF and DFT methods indicates that B3LYP is superior for molecular vibrational problems. Comparison of the simulated spectra with the experimental spectra provides important information about the ability of the computational method to describe the vibrational modes. The RHF and DFT-based NMR calculation procedure was also done. It was used to assign the ¹³C NMR chemical shift of nebivolol.     

Position-resolved free energy of solvation for amino acids in lipid membranes from molecular dynamics simulations

Studies of insertion and interactions of amino acids in lipid membranes are pivotal to our understanding of membrane protein structure and function. Calculating the insertion cost as a function of transmembrane helix sequence is thus an important step towards improved membrane protein prediction and eventually drug design. Here, we present position-dependent free energies of solvation for all amino acid analogs along the membrane normal. The profiles cover the entire region from bulk water to hydrophobic core, and were produced from all-atom molecular dynamics simulations. Experimental differences corresponding to mutations and costs for entire segments match experimental data well, and in addition the profiles provide the spatial resolution currently not available from experiments. Polar side-chains largely maintain their hydration and assume quite ordered conformations, which indicates the solvation cost is mainly entropic. The cost of solvating charged side-chains is not only significantly lower than for implicit solvation models, but also close to experiments, meaning these could well maintain their protonation states inside the membrane. The single notable exception to the experimental agreement is proline, which is quite expensive to introduce in vivo despite its hydrophobicity--a difference possibly explained by kinks making it harder to insert helices in the translocon.     

Testing beta-helix terminal coils stability by targeted substitutions with non-proteogenic amino acids: a molecular dynamics study

The search for new building block templates useful for nanostructures design, targets protein motifs with a wide range of structures. Stabilizing these building blocks when extracted from their natural environment becomes a fundamental goal in order to successfully control their assembly. Targeted replacements of natural residues by conformationally constrained amino acids were shown to be a successful strategy to achieve such stabilization. In this work, the effect of replacing natural amino acids by non-proteogenic residues in a beta-helix building block has been evaluated using extensive molecular dynamics simulations. Here, we focus on systematic substitutions of valine residues present in beta-sheet segments of a beta-helical building block excised from Escherichia coli galactoside acetyltransferase, residues 131-165. Four different types of non-proteogenic amino acids have been considered for substitution: (i) one dehydroamino acid, (ii) two d-amino acids, (iii) one beta-amino acid and (iv) two alpha,alpha-dialkylamino acids. Our results indicate that the ability of non-proteogenic amino acids to stabilize small building block motifs is site-dependent. We conclude that if the replacement does not alter the energy balance between attractive non-covalent interactions and steric hindrance, synthetic residues are suitable candidates to nucleate beta-helix formation.     

Hydration of counterions interacting with DNA double helix: a molecular dynamics study

In the present work, molecular dynamics simulations have been carried out to study the dependence of counterion distribution around the DNA double helix on the character of ion hydration. The simulated systems consisted of DNA fragment d(CGCGAATTCGCG) in water solution with the counterions Na⁺, K⁺, Cs⁺ or Mg²⁺. The characteristic binding sites of the counterions with DNA and the changes in their hydration shell have been determined. The results show that due to the interaction with DNA at least two hydration shells of the counterions undergo changes. The first hydration shell of Na⁺, K⁺, Cs⁺, and Mg²⁺ counterions in the bulk consists of six, seven, ten, and six water molecules, respectively, while the second one has several times higher values. The Mg²⁺ and Na⁺ counterions, constraining water molecules of the first hydration shell, mostly form with DNA water-mediated contacts. In this case the coordination numbers of the first hydration shell do not change, while the coordination numbers of the second one decrease about twofold. The Cs⁺ and K⁺ counterions that do not constrain surrounding water molecules may be easily dehydrated, and when interacting with DNA their first hydration shell may be decreased by three and five water molecules, respectively. Due to the dehydration effect, these counterions can squeeze through the hydration shell of DNA to the bottom of the double helix grooves. The character of ion hydration establishes the correlation between the coordination numbers of the first and the second hydration shells.

Open image in new window

Graphical Abstract Hydration of counterions interacting with DNA double helix

     

Further studies on pyruvoyl amino acids

Pyruvoyl-α-aminoisobutyric acid was prepared, and its absorption characteristics in the ultraviolet compared with those of pyruvoylglycine, pyruvoyl-dl-alanine, and pyruvoyl-dl-phenylalanine. All four compounds in aqueous solution possess the same type of absorption, with maxima at 242–245 mμ and at 310–337 mμ. On alkalinization of the solution, the characteristic absorption of the latter three compounds disappears, and is not restored on acidification. On the other hand, the characteristic absorption of pyruvoyl-α-aminoisobutyric acid changes very little on alkalinization of the solution, and the small change is completely reversed on acidification.These differences have been interpreted as being due to the presence of at least one hydrogen atom on the α-carbon of the amino acid residue in pyruvoylglycine, pyruvoylalanine, and pyruvoylphenylalanine which permits ring closure in alkaline solution to the corresponding γ-hydroxy-pyrrolidonecarboxylic acid derivatives. Where no such hydrogen exists, as in pyruvoyl-α-aminoisobutyric acid, ring closure cannot occur, and the original starting material can be recovered.     

Anticancer and DNA binding studies of potential amino acids based quinazolinone analogs: Synthesis,SAR and molecular docking

A novel series of amino acids conjugated quinazolinone-Schiff’s bases were synthesized and screened for their in vitro anticancer activity and validated by molecular docking and DNA binding studies. In the present investigations, compounds 32, 33, 34, 41, 42 and 43 showed most potent anticancer activity against tested cancer cell lines and DNA binding study using methyl green comparing to doxorubicin and ethidium bromide as a positive control respectively. The structure-activity relationship (SAR) revealed that the tryptophan and phenylalanine derived electron donating groups (OH and OCH₃) favored DNA binding studies and anticancer activity whereas; electron withdrawing groups (Cl, NO₂, and F) showed least anticancer activity. The molecular docking study, binding interactions of the most active compounds 33, 34, 42 and 43 stacked with A–T rich regions of the DNA minor groove by surface binding interactions were confirmed.     

Huntington's disease: studies on brain free amino acids

We studied the levels of free amino acids in putamen and Brodmann's area 10 of 12 patients who died with Huntington's disease and 13 non-neurologic controls. GABA, glutamate and alpha-amino-n-butyric acid concentrations were found to be reduced in putamen of Huntington's disease patients. In Brodmann's area 10 the levels of glutamate, histidine and lysine were decreased, but the content of aspartate, GABA, glycine, serine and taurine was increased in the same group of patients.     

Hydration effect on low-frequency protein dynamics observed in simulated neutron scattering spectra

Hydration effects on protein dynamics were investigated by comparing the frequency dependence of the calculated neutron scattering spectra between full and minimal hydration states at temperatures between 100 and 300 K. The protein boson peak is observed in the frequency range 1-4 meV at 100 K in both states. The peak frequency in the minimal hydration state shifts to lower than that in the full hydration state. Protein motions with a frequency higher than 4 meV were shown to undergo almost harmonic motion in both states at all temperatures simulated, whereas those with a frequency lower than 1 meV dominate the total fluctuations above 220 K and contribute to the origin of the glass-like transition. At 300 K, the boson peak becomes buried in the quasielastic contributions in the full hydration state but is still observed in the minimal hydration state. The boson peak is observed when protein dynamics are trapped within a local minimum of its energy surface. Protein motions, which contribute to the boson peak, are distributed throughout the whole protein. The fine structure of the dynamics structure factor is expected to be detected by the experiment if a high resolution instrument (<∼20 μeV) is developed in the near future.     

Chiroptical studies of fluorescamine labeled amino acids

Absorption and circular dichroism studies of fluorescamine condensation products with α-amino acids, dipeptides and phenylethylamine in the 300–450 nm region are reported. The results make a major revision of the previously suggested rule for absolute configuration determination necessary.     

Thermal perturbation differential spectra of ribonucleic acids. I. Hydration effects.

A relatively important change in UV absorption is observed upon thermal perturbation of nucleotide solutions. Comparison of these thermal perturbation spectra of nucleic acid residues with solvent perturbation spectra of the same compounds suggests that this spectral change can most probably be attributed to temperature induced hydration change of the bases. This conclusion is confirmed by the results obtained from acid-base perturbation spectra of these nucleotides as well as thermal perturbation spectra of nucleotides containing modified bases. It is shown that this temperature dependent change in UV absorption is also present in dinucleoside monophosphates. In that case, this effect is superimposed upon the well known change in absorbance due to the unstacking of the bases during heating.     

Hydration structure and dynamics of B- and Z-DNA in the presence of counterions via molecular dynamics simulations

Following our previous attempts at understanding the structural and dynamical properties of water and counterions hydrating nucleic acids, we have performed molecular dynamics simulations for B- and Z-DNA. In these simulations, the nucleic acids were held rigid. In the case of B-DNA, one turn of B-DNA double helix was considered in the presence of 1500 water molecules and 20 counterions (K⁺). The simulations were performed for 4.0 ps after equilibrating the system. For Z-DNA, we considered one turn of the double helix in the presence of 1851 water molecules and 24 counterions (K⁺). The simulations were carried out for 3.5 ps after equilibration. The average temperature of these simulations was ～ 360 K for Z-DNA and ～ 345 K for B-DNA. In these simulations the hydrogen atoms were explicitly taken into account. For both simulations, a fifth-order predictor-corrector was used for solving the translational equations of motion. The rotational motion of the water molecules was represented in terms of quaternion algebra and the rotational equations of motion were solved with a second-order quaternion method using a sixth-order predictor-corrector method. A time step of 0.5 · 10^?15 s was used in these simulations. The structural and the dynamical properties of water solvating the counterions, and the phosphate groups of the DNA, were computed to understand the hydration structure. Diffusion coefficients and velocity correlation functions were calculated for both ions and the water molecules. The velocity correlation functions for the ions exhibit a caged behavior. The dipole correlation functions for the water molecules indicate that the water molecules close to the helix retain the memory of their initial orientations for longer periods of time than those away from the helix. During the time period of our simulation (3–4 ps) the ion probability distributions show a well-defined pattern and suggest limited mobility for the ions, being close to the helix.