DFT Studies on Schiff Base Formation of Vitamin B6 Analogues

Gaseous inhibition of the transsulfuration pathway by cystathionine β-synthase

The transsulfuration pathway plays a key role in mammals for maintaining the balance between cysteine and homocysteine, whose concentrations are critical in several biochemical processes. Human cystathionine β-synthase is a heme-containing, pyridoxal 5'-phosphate (PLP)-dependent enzyme found in this pathway. The heme group does not participate directly in catalysis, but has a regulatory function, whereby CO or NO binding inhibits the PLP-dependent reactions. In this study, we explore the detailed structural changes responsible for inhibition using quantum chemical calculations to validate the experimentally observed bonding patterns associated with heme CO and NO binding and molecular dynamics simulations to explore the medium-range structural changes triggered by gas binding and propagating to the PLP active site, which is more than 20 Å distant from the heme group. Our results support a previously proposed mechanical signaling model, whereby the cysteine decoordination associated with gas ligand binding leads to breaking of a hydrogen bond with an arginine residue on a neighbouring helix. In turn, this leads to a shift in position of the helix, and hence also of the PLP cofactor, ultimately disrupting a key hydrogen bond that stabilizes the PLP in its catalytically active form.

Physiological Associations between Vitamin B Deficiency and Diabetic Kidney Disease

The number of people living with chronic kidney disease (CKD) is growing as our global population continues to expand. With aging, diabetes, and cardiovascular disease being major harbingers of kidney disease, the number of people diagnosed with diabetic kidney disease (DKD) has grown concurrently. Poor clinical outcomes in DKD could be influenced by an array of factors-inadequate glycemic control, obesity, metabolic acidosis, anemia, cellular senescence, infection and inflammation, cognitive impairment, reduced physical exercise threshold, and, importantly, malnutrition contributing to protein-energy wasting, sarcopenia, and frailty. Amongst the various causes of malnutrition in DKD, the metabolic mechanisms of vitamin B (B1 (Thiamine), B2 (Riboflavin), B3 (Niacin/Nicotinamide), B5 (Pantothenic Acid), B6 (Pyridoxine), B8 (Biotin), B9 (Folate), and B12 (Cobalamin)) deficiency and its clinical impact has garnered greater scientific interest over the past decade. There remains extensive debate on the biochemical intricacies of vitamin B metabolic pathways and how their deficiencies may affect the development of CKD, diabetes, and subsequently DKD, and vice-versa. Our article provides a review of updated evidence on the biochemical and physiological properties of the vitamin B sub-forms in normal states, and how vitamin B deficiency and defects in their metabolic pathways may influence CKD/DKD pathophysiology, and in reverse how CKD/DKD progression may affect vitamin B metabolism. We hope our article increases awareness of vitamin B deficiency in DKD and the complex physiological associations that exist between vitamin B deficiency, diabetes, and CKD. Further research efforts are needed going forward to address the knowledge gaps on this topic.

A DFT study of the active role of the phosphate group of an internal aldimine in a transamination reaction

A transamination reaction from an internal aldimine ([PLP]) and (S)-alanine to pyridoxamine phosphate (PMP) and pyruvic acid was investigated by DFT calculations. As [PLP], a model where the lysine (-Lys) part was approximated by -CH[-NH-C(O)-CH₃]-C(O)-NH-CH₃ was adopted. (H₂O)₄ was also included to trace reaction paths involving proton transfers. 13 elementary processes were obtained. For (the external aldimine → quinoid), (quinoid → ketimine) and (ketimine → carbinol amine) processes, the water dimer was found to connect a phosphate-group oxygen with the moving proton. The connection promoted the Grotthuss-type proton transfer in transition states. It was revealed that the phosphate group is not a mere substituent but has the central role in the transfer.

The conversion of L-lysine into L-β-lysine: the role of 5'-deoxyadenosyl radical and water-a DFT study

In the current study, we deal with the crucial role of 5'-deoxyadenosyl radical and water in the mechanism of the conversion of L-lysine into L-β-lysine. The DFT (density functional theory)-B3LYP method coupled with 6-31G(d) basis set has been performed to investigate the optimized structures of transition states (TSs) and intermediates (IMs) of two processes in water: (i) the attack of 5'-deoxyadenosyl radical to complex PLP-L-lysine and (ii) hydrolysis to liberate L-β-lysine. Meanwhile, M062X/6-311++g(3df,2p) level of theory is applied to compute the relative Gibbs energy ΔG. Procedure (i) has undergone various steps but includes two main structural aziridinyl rings TS2 (ΔG = 4.1 kcal/mol) and TS3 (ΔG = 2.3 kcal/mol). In stage (ii), hydroxy group of water would help to break the bond between β-NH₂ group of L-lysine and PLP better than that of Tyr389.

An applied quantum-chemical model for genipin-crosslinked chitosan (GCS) nanocarrier

The genipin-crosslinked chitosan (GCS) nanocarrier has received a lot of attention due to its unique biological and chemical properties as an effective drug delivery system. GCS was modeled by considering two chitosan (CS) polymer sequences with six monomer units that are crosslinked by genipin. To investigate the characteristics of this model, we considered it as a nanocarrier of the anti-cancer drug cladribine (2CdA). Seven configurations of GCS and 2CdA (GCS/2CdA1-7) were optimized at M06-2X/6-31G(d,p) in aqueous solution. The average binding energy above 100 kJ mol^-1 indicates a high drug loading amount. The high adsorption of the drug on GCS is due to the hydrogen bonds that were investigated by AIM analysis. Hydrogen bonds also allow the drug to be released more slowly. These results were confirmed by experimental evidence and the comparison of this model with the simple model of one polymer chain. Also, the mechanism of GCS formation was investigated by calculating the activation parameters, which indicates that solvent (H₂O) molecules are explicitly involved in the formation of GCS.

Structural insights into the mechanism of internal aldimine formation and catalytic loop dynamics in an archaeal Group II decarboxylase

Formation of the internal aldimine (LLP) is the first regulatory step that activates pyridoxal 5'-phosphate (PLP) dependent enzymes. The process involves a nucleophilic attack on PLP by an active site Lys residue, followed by proton transfers resulting in a carbinolamine (CBA) intermediate that undergoes dehydration to form the aldimine. Despite a general understanding of the pathway, the structural basis of the mechanistic roles of specific residues in each of these steps is unclear. Here we determined the crystal structure of the LLP form (holo-form) of a Group II PLP-dependent decarboxylase from Methanocaldococcus jannaschii (MjDC) at 1.7 Å resolution. By comparing the crystal structure of MjDC in the LLP form with that of the pyridoxal-P (non-covalently bound aldehyde) form, we demonstrate structural evidence for a water-mediated mechanism of LLP formation. A conserved extended hydrogen-bonding network around PLP coupled to the pyridinyl nitrogen influences activation and catalysis by affecting the electronic configuration of PLP. Furthermore, the two cofactor bound forms revealed open and closed conformations of the catalytic loop (CL) in the absence of a ligand, supporting a hypothesis for a regulatory link between LLP formation and CL dynamics. The evidence suggests that activation of Group II decarboxylases involves a complex interplay of interactions between the electronic states of PLP, the active site micro-environment and CL dynamics.

New views on the reaction of primary amine and aldehyde from DFT study

A general theoretical investigation on the reaction of primary amine with aldehyde was carried out by density functional theory. The calculation systems involve three kinds of primary amines (methylamine, vinylamine, and phenylamine) and three kinds of aldehydes (formaldehyde, acetaldehyde, and acrylaldehyde). The steric and electronic inductive effects on the reaction mechanism were studied. Results reveal that the nucleophilic attack of primary amine on aldehyde under neutral conditions leads to carbinolamines, rather than Schiff bases. The nucleophilic attack on the protonated aldehyde produces the protonated Schiff base. The steric hindrance of the aldehyde slows down the nucleophilic attack but allows enough time to abstract a H; consequently, the formation of the protonated Schiff base is preferred. During the carbinolamine protonation, the H(+) preferably locates on the amine nitrogen and then is abstracted by the hydroxyl oxygen over an energy barrier, leaving protonated Schiff base after timely water liberation. The formation of a prereaction potential energy well obviously softens the steric and electronic inductive effects on the active barrier for different reactants.

Nonenzymatic Reactions above Phospholipid Surfaces of Biological Membranes: Reactivity of Phospholipids and Their Oxidation Derivatives

Phospholipids play multiple and essential roles in cells, as components of biological membranes. Although phospholipid bilayers provide the supporting matrix and surface for many enzymatic reactions, their inherent reactivity and possible catalytic role have not been highlighted. As other biomolecules, phospholipids are frequent targets of nonenzymatic modifications by reactive substances including oxidants and glycating agents which conduct to the formation of advanced lipoxidation end products (ALEs) and advanced glycation end products (AGEs). There are some theoretical studies about the mechanisms of reactions related to these processes on phosphatidylethanolamine surfaces, which hypothesize that cell membrane phospholipids surface environment could enhance some reactions through a catalyst effect. On the other hand, the phospholipid bilayers are susceptible to oxidative damage by oxidant agents as reactive oxygen species (ROS). Molecular dynamics simulations performed on phospholipid bilayers models, which include modified phospholipids by these reactions and subsequent reactions that conduct to formation of ALEs and AGEs, have revealed changes in the molecular interactions and biophysical properties of these bilayers as consequence of these reactions. Then, more studies are desirable which could correlate the biophysics of modified phospholipids with metabolism in processes such as aging and diseases such as diabetes, atherosclerosis, and Alzheimer's disease.

A DFT study of the carboxymethyl-phosphatidylethanolamine formation from glyoxal and phosphatidylethanolamine surface. Comparison with the formation of N(ε)-(carboxymethyl)lysine from glyoxal and l-lysine

Mechanisms of the generation of CML and CM-PE from the reactions between glyoxal and l-lysine, and glyoxal and phosphatidylethanolamine (PE) were studied using the DFT method.

DFT study of the mechanism of the reaction of aminoguanidine with methylglyoxal

We have studied the mechanism of the reaction between aminoguanidine (AG) and methylglyoxal (MG) by carrying out Dmol3/DFT calculations, obtaining intermediates, transition-state structures, and free-energy profiles for all of the elementary steps of the reaction. Designed models included explicit water solvent, which forms hydrogen-bond networks around the reactants and intermediate molecules, facilitating intramolecular proton transfer in some steps of the reaction mechanism. The reaction take place in four steps, namely: (1) formation of a guanylhydrazone-acetylcarbinol adduct by condensation of AG and MG; (2) dehydration of the adduct; (3) formation of an 1,2,4-triazine derivative by ring closure; and (4) dehydration with the formation of 5-methyl 3-amino-1,2,4-triazine as the final product. From a microkinetic point of view, the first dehydration step was found to be the rate-determining step for the reaction, with the reaction having an apparent activation energy of 12.65 kcal mol⁻¹. Additionally, some analogous structures of intermediates and transition states for the reaction between AG and 2,3-dicarbonyl-phosphatidylethanolamine, a possible intermediate in Amadori-glycated phosphatidylethanolamine (Amadori-PE) autooxidation, were obtained to evaluate the reaction above a phosphatidylethanolamine (PE) surface. Our results are in agreement with experimental results obtaining by other authors, showing that AG is efficient at trapping dicarbonyl compounds such as methylglyoxal, and by extension these compounds joined to biomolecules such as PE in environments such as surfaces and their aqueous surroundings.

Experimental and theoretical studies of vibrational spectrum and molecular structure and related properties of pyridoxine (vitamin B6)

Vibrational spectrum of pyridoxine has been investigated using experimental IR and Raman spectroscopic and density functional theory (DFT) methods. Vibrational assignments of the observed IR and Raman bands have been proposed in light of the results obtained from computations. In order to assign the observed IR and Raman frequencies the potential energy distributions (PEDs) have also been computed. Optimized geometrical parameters suggest that if the OH groups of the two methyl groups are replaced by H atoms the resulting molecule has Cs point group symmetry. To investigate molecular stability and bond strength we have used natural bond orbital analysis (NBO). Charge transfer occurs in the molecule have been shown by the calculated highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) energies. The mapping of electron density iso-surface with electrostatic potential (ESP), has been carried out to get the information about the size, shape, charge density distribution and site of chemical reactivity of the molecule.

Solvent-catalyzed ring-chain-ring tautomerization in axially chiral compounds

The mechanism of ring-chain-ring tautomerization and the prominent effect of the solvent environment have been computationally investigated in an effort to explain the enantiomeric interconversion observed in 2-oxazolidinone derivatives, heterocyclic analogues of biphenyl atropisomers, which were isolated as single stable enantiomers and have the potential to be used as axially chiral catalysts. This study has shed light on the identity of the intermediate species involved in the ring-chain-ring tautomerization process as well as the catalytic effect of polar protic solvents. These mechanistic details will prove very useful in predicting and understanding ring-chain tautomeric equilibria in similar heterocyclic systems and will further enable experimentalists to devise appropriate experimental conditions in which axially chiral catalysts remain stable as single enantiomers.

Microsolvated transition state models for improved insight into chemical properties and reaction mechanisms

Over the years, several methods have been developed to effectively represent the chemical behavior of solutes in solvents. The environmental effects arising due to solvation can generally be achieved either through inclusion of discrete solvent molecules or by inscribing into a cavity in a homogeneous and continuum dielectric medium. In both these approaches of computational origin, the perturbations on the solute induced by the surrounding solvent are at the focus of the problem. While the rigor and method of inclusion of solvent effects vary, such solvation models have found widespread applications, as evident from modern chemical literature. A hybrid method, commonly referred to as cluster-continuum model (CCM), brings together the key advantages of discrete and continuum models. In this perspective, we intend to highlight the latent potential of CCM toward obtaining accurate estimates on a number of properties as well as reactions of contemporary significance. The objective has generally been achieved by choosing illustrative examples from the literature, besides expending efforts to bring out the complementary advantages of CCM as compared to continuum or discrete solvation models. The majority of examples emanate from the prevalent applications of CCM to organic reactions, although a handful of interesting organometallic reactions have also been discussed. In addition, increasingly accurate computations of properties like pK(a) and solvation of ions obtained using the CCM protocol are also presented.

Theoretical study on HF elimination and aromatization mechanisms: a case of pyridoxal 5' phosphate-dependent enzyme

Pyridoxal 5-phosphate (PLP), the phosphorylated and the oxidized form of vitamin B6 is an organic cofactor. PLP forms a Schiff base with the ϵ-amino group of a lysine residue of PLP-dependent enzymes. γ-Aminobutyric acid (GABA) aminotransferase is a PLP-dependent enzyme that degrades GABA to succinic semialdehyde, while reduction of GABA concentration in the brain causes convolution besides several neurological diseases. The fluorine-containing substrate analogues for the inactivation of the GABA-AT are synthesized extensively in cases where the inactivation mechanisms involve HF elimination. Although two proposed mechanisms are present for the HF elimination, the details of the base-induced HF elimination are not well identified. In this density functional theory (DFT) study, fluorine-containing substrate analogue, 5-amino-2-fluorocyclohex-3-enecarboxylic acid, is particularly chosen in order to explain the details of the HF elimination reactions. On the other hand, the experimental studies revealed that aromatization competes with Michael addition mechanism in the presence of 5-amino-2-fluorocyclohex-3-enecarboxylic acid. The results allowed us to draw a conclusion for the nature of HF elimination, besides the elucidation of the mechanism preference for the inactivation mechanism. Furthermore, the solvent phase calculations carried out in this study ensure that the proton transfer steps should be assisted either by a water molecule or a base for lower activation energy barriers.

Studies on molecular structure and tautomerism of a vitamin B6 analog with density functional theory

This work presents a computational study on the molecular structure and tautomeric equilibria of a novel Schiff base L derived from pyridoxal (PL) and o-phenylenediamine by using the density functional method B3LYP with basis sets 6-31 G(d,p), 6-31++G(d,p), 6-311 G(d,p) and 6-311++G(d,p). The optimized geometrical parameters obtained by B3LYP/6-31 G(d,p) method showed the best agreement with the experimental values. Tautomeric stability study of L inferred that the enolimine form is more stable than its ketoenamine form in both gas phase and solution. However, protonation of the pyridoxal nitrogen atom (LH) have accelerated the formation of ketoenamine form, and therefore, both ketoenamine and enolimine forms could be present in acidic media.

Reactivity of a phospholipid monolayer model under periodic boundary conditions: a density functional theory study of the Schiff base formation between phosphatidylethanolamine and acetaldehyde

A mechanism for the formation of the Schiff base between an acetaldehyde and an amine-phospholipid monolayer model based on Dmol3/density functional theory calculations under periodic boundary conditions was constructed. This is the first time such a system has been modeled to examine its chemical reactivity at this computation level. Each unit cell contains two phospholipid molecules, one acetaldehyde molecule, and nine water molecules. One of the amine-phospholipid molecules in the cell possesses a neutral amino group that is used to model the nucleophilic attack on the carboxyl group of acetaldehyde, whereas the other has a charged amino group acting as a proton donor. The nine water molecules form a hydrogen bond network along the polar heads of the phospholipids that facilitates very fast proton conduction at the interface. Using periodic boundary conditions afforded proton transfer between different cells. The reaction takes place in two steps, namely, (1) formation of a carbinolamine and (2) its dehydration to the Schiff base. The carbinolamine is the primary reaction intermediate, and dehydration is the rate-determining step of the process, consistent with available experimental evidence for similar reactions. On the basis of the results, the cell membrane surface environment may boost phospholipid glycation via a neighboring catalyst effect.

Water-Assisted Transamination of Glycine and Formaldehyde

A computational study on the transamination reaction of molecular complexes that consist of NH2CH2COOH + CH2O + nH2O, where n = 0, 1, 2, is presented. This work has allowed the description of the geometries of all the intermediates and transition states of the reactions, which can be described by five steps: carbinolamine formation, dehydration, 1,3 proton transfer, hydrolysis, and carbinolamine elimination. Among the five steps of the reaction, hydrolysis and elimination occur with the existence of general acid catalysis related to the carboxylic group. The water molecules can be involved in the reaction by performing as a proton-transfer carrier and a stabilizing zwitterion. It can be predicted from our calculations that in the transamination between alpha-amino acids and alpha-keto acids, the carbinolamine is formed with small barrier or even barrierless while the dehydration occurs easily at room temperature. However, without heating the 1,3 proton transfer could not occur as the barrier is 26.7 kcal/mol relative to the reactant complex when including two water molecules. Our results are in good agreement with experimental conclusions.