Unravelling carbohydrate binding module 21 (CBM21) dynamics of interaction with amylose.

The carbohydrate binding module 21 (CBM21) from Rhizopus oryzae is a dual-site CBM proposed to disrupt polysaccharide structures. Additionally, it serves as a purification tag in industry. CBM21 crystal structure features a Glc residue in an unusual 1S3 conformation, whose relevance for the CBM mechanism of action is unclear. In this context, we seek to contribute for the understanding of CBM21 mechanism of action by: i) investigating the role of the 1S3 conformation on carbohydrate recognition, and ii) characterize the protein-carbohydrate binding dynamics using molecular dynamics and metadynamics simulations at MM and QM/MM levels. Results indicate the 1S3 Glc conformation is unlikely to occur under biological conditions, being originated from the crystallographic environment. CBM21 binding to small ligands appears transient and unstable, while protein dimerization and polysaccharide chain size influence complex stability. In interactions with amylose, CBM21 exhibits a repeated unbinding followed by re-binding, while simultaneously alternating between binding sites I and II. These results suggest that CBM21 acts through transient interactions, directing carbohydrates to the catalytic center rather than forming strong and long-lasting bonds with carbohydrates. Accordingly, we expect such atomistic depiction of CBM21 mechanism could aid in CBM design targeting biotechnological applications.

Cycloaddition reactions via "on water" protocol reactions: A density functional theory study.

In this work, the reactions of quadricyclane with dimethyl azodicarboxylate (DMAD) and of quadricyclane with diethyl azodicarboxylate (DEAD) in gas phase and in water environments were studied by a first-principles investigation within the framework of auxiliary density functional theory (ADFT). For these type of organic reactions is known that water is required to accelerate them. Since the reason of why this occur is still unknown, this work aims to gain insight into this reaction mechanism. For this investigation, the generalized gradient approximation as well as a hybrid functional were employed. The obtained optimized structures for the reactants, of the products and of the transition states are reported, together with the corresponding frequency analysis results and the reaction profiles. Along the proposed concerted reaction mechanism, a critical points search of the electron density and a charge analysis were performed. The calculated potential energy barriers of these reactions in gas phase and in water environments are compared. In agreement with experiment, the obtained results indicate that both reactions occur faster in water than in gas phase. This study shows that there is a change in the polarity of the two most important carbon atoms of the formed compounds along the reactions and that the decrease of the activation energy barrier which occurs in liquid phase in these reactions is because the structures of the main transition states are stabilized by the water environment. Therefore, the here obtained results demonstrate the important role played by the water-molecule framework into the activation energy barrier and structures of the molecules that participate in the DMAD and DEAD cycloaddition reactions.

Cartesian constraints in QM/MM optimizations.

With the rise of quantum mechanical/molecular mechanical (QM/MM) methods, the interest in the calculation of molecular assemblies has increased considerably. The structures and dynamics of such assemblies are usually governed to a large extend by intermolecular interactions. As a result, the corresponding potential energy surfaces are topological rich and possess many shallow minima. Therefore, local structure optimizations of QM/MM molecular assemblies can be challenging, in particular if optimization constraints are imposed. To overcome this problem, structure optimization in normal coordinate space is advocated. To do so, the external degrees of freedom of a molecule are separated from the internal ones by a projector matrix in the space of the Cartesian coordinates. Here we extend this approach to Cartesian constraints. To this end, we devise an algorithm that adds the Cartesian constraints directly to the projector matrix and in this way eliminates them from the reduced coordinate space in which the molecule is optimized. To analyze the performance and stability of the constrained optimization algorithm in normal coordinate space, we present constrained minimizations of small molecular systems and amino acids in gas phase as well as water employing QM/MM constrained optimizations. All calculations are performed in the framework of auxiliary density functional theory as implemented in the program deMon2k.

Binding mechanism of disulfide species to ferric hemeproteins: The case of metmyoglobin.

The interactions of the heme iron of hemeproteins with sulfide and disulfide compounds are of potential interest as physiological signaling processes. While the interaction with hydrogen sulfide has been described computationally and experimentally, the reaction with disulfide, and specifically the molecular mechanism for ligand binding has not been studied in detail. In this work, we study the association process for disulfane and its conjugate base disulfanide at different pH conditions. Additionally, by means of advanced sampling techniques based on multiple steered molecular dynamics, we provide free energy profiles for ligand migration for both acid/base species, showing a similar behavior to the previously reported for the related H2S/HS¯ pair. Finally, we studied the ligand interchange reaction (H2O/H2S, HS¯ and H2O/HSSH, HSS¯) by means of hybrid quantum mechanics-molecular mechanics calculations. We show that the anionic species are able to displace more efficiently the H2O bound to the iron, and that the H-bond network in the distal cavity can help the neutral species to perform the reaction. Altogether, we provide a molecular explanation for the experimental information and show that the global association process depends on a fine balance between the migration towards the active site and the ligand interchange reaction.

Photophysics of tz Adenine and tz Guanine fluorescent nucleobases embedded into DNA and RNA.

UV-VIS photoinduced events of tz A and tz G embedded into DNA and RNA are described by combining the Extended Multi-State Second-Order Perturbation Theory (XMS-CASPT2) and electrostatic embedding molecular mechanics methods (QM/MM). Our results point out that the S1 1 (ππ* La ) state is the bright state in both environments. After the photoexcitation to the S1 1 (ππ* La ) state, the electronic population evolves barrierless towards its minimum, from where the excess of energy can be dissipated by fluorescence. As the minimum energy crossing point structure between the ground and first bright states lies in a high-energy region, the direct internal conversion to the ground state is an unviable mechanism. Other spectroscopic properties (for instance, absorption and Stokes shifts) and comparisons with photochemical properties of canonical nucleobases are also provided.

Unveiling the tartrazine binding mode with ds-DNA by UV-visible spectroscopy, electrochemical, and QM/MM methods.

Here, we studied the interaction between the food colorant tartrazine (TZ) and double stranded DNA (dsDNA), using spectroscopic, electrochemical, and computational methods such as QM/MM combined with TD-DFT. Despite the UV-vis spectroscopy is widely used to study the interaction between molecules, for the case of TZ there are discrepancies in the analyses presented in the literature available, presenting both hyperchromic and hypochromic effects and consequently different rationalizations for their results. Herein we propose the combination of UV-vis experiments with the design of high-level computational models capable of reproducing the experimental behavior to finally define the proper binding mode at the molecular scale together with the rationalization of the experimental optical response due to the complex formation. To complement the UV-vis experiments, we propose the use of electrochemical measurements, to support the results obtained through UV-vis spectroscopy, as it has been successfully used for the determination of interaction modes between small molecules and biomolecules in any condition. Our UV-vis spectroscopy experiments showed only a hypochromic effect of the absorption spectra of TZ after interaction with DNA, indicative of TZ being deeply buried in the DNA structure. The effect of ionic strength in the experimental procedures led to the dissociation of TZ, thus indicating that the interaction mode was groove binding. On the other hand, the electrochemical studies showed an irreversible reduction peak of TZ, which after the interaction with DNA exhibited a positive shift in potential that can be attributed to groove binding. The binding constant for TZ-DNA was calculated as 4.45x104M-1 (UV-vis) and 5.75x104M-1 (electrochemistry), in line with other groove binder azo dyes. Finally, through the QM/MM calculations we found that the minor-groove binding mode interacting in zones rich in adenine and thymine was the model best suited to reproduce the experimental UV-vis response.

Computational vaccinology guided design of multi-epitope subunit vaccine against a neglected arbovirus of the Americas.

Mayaro virus (MAYV) is an arbovirus found in the Americas that can cause debilitating arthritogenic disease. Although it is an emerging virus, the only current approach is vector control, as there are no approved vaccines to prevent MAYV infection nor therapeutics to treat it. In search of an effective vaccine candidate against MAYV, we used immunoinformatics and molecular modeling to attempt to identify promiscuous T-cell epitopes of the nonstructural polyproteins (nsP1, nsP2, nsP3, and nsP4) from 127 MAYV genomes sequenced in the Americas (08 Bolivia, 72 Brazil, 04 French Guiana, 05 Haiti, 20 Peru, 04 Trinidad and Tobago, and 14 Venezuela). For this purpose, consensus sequences of 360 proteins were used to identify short protein sequences that can bind to MHC I class (MHC II). Our analysis revealed 56 potential MHC-I/TCD8+ (29 MHC-II/TCD4+) epitopes, but only 6 (16) TCD8+ (TCD4+) epitopes showed high antigenicity and conservation, non-allergenicity, non-toxicity, and excellent population coverage. Finally, classical and quantum mechanical calculations (QM:MM) were used to improve the quality of the docking calculations, with the QM part of the simulations performed using the density functional theory formalism (DFT). These results provide insights for the advancement of diagnostic platforms, vaccine development, and immunotherapeutic interventions.Communicated by Ramaswamy H. Sarma.

Unraveling the Role of the Tyrosine Tetrad from the Binding Site of the Epigenetic Writer MLL3 in the Catalytic Mechanism and Methylation Multiplicity.

MLL3, also known as KMT2C, is a lysine mono-methyltransferase in charge of the writing of an epigenetic mark on lysine 4 from histone 3. The catalytic site of MLL3 is composed of four tyrosines, namely, Y44, Y69, Y128, and Y130. Tyrosine residues are highly conserved among lysine methyltransferases' catalytic sites, although their complete function is still unclear. The exploration of how modifications on these residues from the enzymatic machinery impact the enzymatic activity of MLL3 could shed light transversally into the inner functioning of enzymes with similar characteristics. Through the use of QMMM calculations, we focus on the effect of the mutation of each tyrosine from the catalytic site on the enzymatic activity and the product specificity in the current study. While we found that the mutations of Y44 and Y128 by phenylalanine inactivated the enzyme, the mutation of Y128 by alanine reactivated the enzymatic activity of MLL3. Moreover, according to our models, the Y128A mutant was even found to be capable of di- and tri-methylate lysine 4 from histone 3, what would represent a gain of function mutation, and could be responsible for the development of diseases. Finally, we were able to establish the inactivation mechanism, which involved the use of Y130 as a water occlusion structure, whose conformation, once perturbed by its mutation or Y128 mutant, allows the access of water molecules that sequester the electron pair from lysine 4 avoiding its methylation process and, thus, increasing the barrier height.

Explaining the interaction of mangiferin with MMP-9 and NF-ƙß: a computational study.

Mangiferin is a glycosylated xanthone widely distributed in nature, which exhibits wide pharmacological activities, highlighting its anti-cancer properties. Mangiferin interferes with inflammation, lipid, and calcium signaling, which selectively inhibits multiple NFkB target genes as interleukin-6, tumor necrosis factor, plasminogen, and matrix metalloproteinase, among others. In this work, the interactions of this polyphenol with MMP-9 and NF-κß are characterized by using computational chemistry methods. The results show MMP-9 inhibition by mangiferina is characterized for the interact with the catalytic Zn atom through a penta-coordinate structure. It is also demonstrated through a strong charge transfer established between mangiferin and Zn in the QM/MM study. Concerning the mangiferin/NF-κß system, the 92.3% of interactions between p50 sub-unity and DNA are maintained with a binding energy of - 8.04 kcal/mol. These findings indicate that mangiferin blocks the p50-p65/DNA interaction resulting in the loss of the functions of this hetero-dimeric member and suggesting inhibition of the cancer progression. Experimental results concerning the anti-cancer properties of mangiferin show that this natural compound can inhibit selectively MMP-9 and NF-ƙß. Although the anti-tumor properties of mangiferin are well defined, its molecular mechanisms of actions are not described. In this work, a computational study is carried out to characterize the interactions of mangiferin with these molecular targets. The results obtained corroborate the anti-proliferative and anti-apoptotic activity of mangiferin and provide a depiction of its mechanisms of action.

A rationale for the unlike potency of avibactam and ETX2514 against OXA-24 ß-lactamase.

Diazabicyclooctanone inhibitors such as ETX2514 and avibactam have shown enhanced inhibitory performance to fight the antibiotic resistance developed by pathogens. However, avibactam is ineffective against Acinetobacter baumannii infections, unlike ETX2514. The molecular basis for this difference has not been tackled from a molecular approach, precluding the knowledge of relevant information. In this article, the mechanisms involved in the inhibition of OXA-24 by ETX2514 and avibactam are studied theoretically by hybrid QM/MM calculations. The results show that both inhibitors share the same inhibition mechanisms, comprising acylation a deacylation stages. The involved mechanisms include the same number of steps, transition states and intermediates; although they differ in the involved activation barriers. This difference accounts for the dissimilar inhibitory ability of both inhibitors. The molecular reason for this is the endocyclic double bond in the piperidine ring of ETX2514 increasing the ring strain and chemical reactivity on the N6 and C7 atoms, besides the methyl substituent, which enhance the hydrophobic character of the ring. Furthermore, Lys218 and the carboxylated Lys84 of ETX2514, play a crucial role in the mechanism by coordinating their protonation states in an on/off (protonated/deprotonated) manner, favoring the proton transference between the residues and the inhibitor.

Self-Organization Dynamics of Collagen-like Peptides Crosslinking Is Driven by Rose-Bengal-Mediated Electrostatic Bridges.

The present work focuses on the computational study of the structural micro-organization of hydrogels based on collagen-like peptides (CLPs) in complex with Rose Bengal (RB). In previous studies, these hydrogels computationally and experimentally demonstrated that when RB was activated by green light, it could generate forms of stable crosslinked structures capable of regenerating biological tissues such as the skin and cornea. Here, we focus on the structural and atomic interactions of two collagen-like peptides (collagen-like peptide I (CLPI), and collagen-like peptide II, (CLPII)) in the presence and absence of RB, highlighting the acquired three-dimensional organization and going deep into the stabilization effect caused by the dye. Our results suggest that the dye could generate a ternary ground-state complex between collagen-like peptide fibers, specifically with positively charged amino acids (Lys in CLPI and Arg in CLPII), thus stabilizing ordered three-dimensional structures. The discoveries generated in this study provide the structural and atomic bases for the subsequent rational development of new synthetic peptides with improved characteristics for applications in the regeneration of biological tissues during photochemical tissue bonding therapies.

A theoretical approach for the acylation/deacylation mechanisms of avibactam in the reversible inhibition of KPC-2.

Klebsiella pneumoniae carbapenemase (KPC-2) is the most commonly encountered class A ß-lactamase variant worldwide, which confer high-level resistance to most available antibiotics. In this article we address the issue by a combined approach involving molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics calculations. The study contributes to improve the understanding, at molecular level, of the acylation and deacylation stages of avibactam involved in the inhibition of KPC-2. The results show that both mechanisms, acylation and deacylation, the reaction occur via the formation of a tetrahedral intermediate. The formation of this intermediate corresponds to the rate limiting stage. The activation barriers are 19.5 kcal/mol and 23.0 kcal/mol for the acylation and deacylation stages, respectively. The associated rate constants calculated, using the Eyring equation, are 1.2 × 10-1 and 3.9 × 10-4 (s-1). These values allow estimating a value of 3.3 × 10-3 for the inhibition constant, in good agreement with the experimental value.

Neutron structures of Leishmania mexicana triosephosphate isomerase in complex with reaction-intermediate mimics shed light on the proton-shuttling steps.

Triosephosphate isomerase (TIM) is a key enzyme in glycolysis that catalyses the interconversion of glyceraldehyde 3-phosphate and dihydroxy-acetone phosphate. This simple reaction involves the shuttling of protons mediated by protolysable side chains. The catalytic power of TIM is thought to stem from its ability to facilitate the deprotonation of a carbon next to a carbonyl group to generate an enediolate intermediate. The enediolate intermediate is believed to be mimicked by the inhibitor 2-phosphoglycolate (PGA) and the subsequent enediol intermediate by phosphoglycolohydroxamate (PGH). Here, neutron structures of Leishmania mexicana TIM have been determined with both inhibitors, and joint neutron/X-ray refinement followed by quantum refinement has been performed. The structures show that in the PGA complex the postulated general base Glu167 is protonated, while in the PGH complex it remains deprotonated. The deuteron is clearly localized on Glu167 in the PGA-TIM structure, suggesting an asymmetric hydrogen bond instead of a low-barrier hydrogen bond. The full picture of the active-site protonation states allowed an investigation of the reaction mechanism using density-functional theory calculations.

Computational modeling of carbohydrate processing enzymes reactions.

Carbohydrate processing enzymes are of biocatalytic interest. Glycoside hydrolases and the recently discovered lytic polysaccharide monooxygenase for their use in biomass degradation to obtain biofuels or valued chemical entities. Glycosyltransferases or engineered glycosidases and phosphorylases for the synthesis of carbohydrates and glycosylated products. Quantum mechanics-molecular mechanics (QM/MM) methods are highly contributing to establish their different chemical reaction mechanisms. Other computational methods are also used to study enzyme conformational changes, ligand pathways, and processivity, e.g. for processive glycosidases like cellobiohydrolases. There is still a long road to travel to fully understand the role of conformational dynamics in enzyme activity and also to disclose the variety of reaction mechanisms these enzymes employ. Additionally, computational tools for enzyme engineering are beginning to be applied to evaluate substrate specificity or aid in the design of new biocatalysts with increased thermostability or tailored activity, a growing field where molecular modeling is finding its way.

Distal lysine (de)coordination in the algal hemoglobin THB1: A combined computer simulation and experimental study.

THB1 is a monomeric truncated hemoglobin from the green alga Chlamydomonas reinhardtii. In the absence of exogenous ligands and at neutral pH, the heme group of THB1 is coordinated by two protein residues, Lys53 and His77. THB1 is thought to function as a nitric oxide dioxygenase, and the distal binding of O2 requires the cleavage of the Fe-Lys53 bond accompanied by protonation and expulsion of the lysine from the heme cavity into the solvent. Nuclear magnetic resonance spectroscopy and crystallographic data have provided dynamic and structural insights of the process, but the details of the mechanism have not been fully elucidated. We applied a combination of computer simulations and site-directed mutagenesis experiments to shed light on this issue. Molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics restrained optimizations were performed to explore the nature of the transition between the decoordinated and lysine-bound states of the ferrous heme in THB1. Lys49 and Arg52, which form ionic interactions with the heme propionates in the X-ray structure of lysine-bound THB1, were observed to assist in maintaining Lys53 inside the protein cavity and play a key role in the transition. Lys49Ala, Arg52Ala and Lys49Ala/Arg52Ala THB1 variants were prepared, and the consequences of the replacements on the Lys (de)coordination equilibrium were characterized experimentally for comparison with computational prediction. The results reinforced the dynamic role of protein-propionate interactions and strongly suggested that cleavage of the Fe-Lys53 bond and ensuing conformational rearrangement is facilitated by protonation of the amino group inside the distal cavity.

Mechanistic study of the biosynthesis of R-phenylcarbinol by acetohydroxyacid synthase enzyme using hybrid quantum mechanics/molecular mechanics simulations.

The biosynthesis of R-phenylacetylcarbinol (R-PAC) by the acetohydroxy acid synthase, (AHAS) is addressed by molecular dynamics simulations (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and QM/MM free energy calculations. The results show the reaction starts with the nucleophilic attack of the C2α atom of the HEThDP intermediate on the Cß atom of the carbonyl group of benzaldehyde substrate via the formation of a transition state (TS1) with the HEThDP intermediate under 4'-aminopyrimidium (APH+) form. The calculated activation free energy for this step is 17.4kcal mol-1 at 27 °C. From this point, the reaction continues with the abstraction of Hß atom of the HEThDP intermediate by the Oß atom of benzaldehyde to form the intermediate I. The reaction is completed with the cleavage of the bond C2α-C2 to form the product R-PAC and to regenerate the ylide intermediate under the APH+ form, allowing in this way to reinitiate to the catalytic cycle once more. The calculated activation barrier for this last step is 15.9kcal mol-1 at 27 °C.

A QM/MM study of the conformation stability and electronic structure of the photochromic switches derivatives of DHA/VHF in acetonitrile solution.

We present a detailed theoretical study of the electronic absorption spectra and thermochemistry of molecular photoswitches composed of one and two photochromic units of dihydroazulene (DHA)/vinylheptafulvene (VHF) molecules. Six different isomers are considered depending on the ring opening/closure forms of the DHA units. The solvent effect of acetonitrile is investigated using a sequential Molecular Mechanics/Quantum Mechanics approach. The thermochemical investigations of these photochromic molecules were performed using the Free Energy Perturbation method, and the simulations were performed using Configurational Bias Monte Carlo. We show that to open the 5-member ring of the DHA, there is no significant gain in thermal release of energy for the back reaction when a unit or two DHA units are considered. Overall, we found agreement between the solvation free energy based on Monte Carlo simulations and the continuum solvent model. However, the cavitation term in the continuum model is shown to be a source of disagreement when the non-electrostatic terms are compared. The electronic absorption spectra are calculated using TDDFT CAM-B3LYP/cc-pVDZ. Agreement with experiment is obtained within 0.1 eV, considering statistically uncorrelated configurations from the simulations. Inhomogeneous broadening is also considered and found to be well described in all cases.

Exploring the Catalytic Mechanism of the RNA Cap Modification by nsp16-nsp10 Complex of SARS-CoV-2 through a QM/MM Approach.

The inhibition of key enzymes that may contain the viral replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have assumed central importance in drug discovery projects. Nonstructural proteins (nsps) are essential for RNA capping and coronavirus replication since it protects the virus from host innate immune restriction. In particular, nonstructural protein 16 (nsp16) in complex with nsp10 is a Cap-0 binding enzyme. The heterodimer formed by nsp16-nsp10 methylates the 5'-end of virally encoded mRNAs to mimic cellular mRNAs and thus it is one of the enzymes that is a potential target for antiviral therapy. In this study, we have evaluated the mechanism of the 2'-O methylation of the viral mRNA cap using hybrid quantum mechanics/molecular mechanics (QM/MM) approach. It was found that the calculated free energy barriers obtained at M062X/6-31+G(d,p) is in agreement with experimental observations. Overall, we provide a detailed molecular analysis of the catalytic mechanism involving the 2'-O methylation of the viral mRNA cap and, as expected, the results demonstrate that the TS stabilization is critical for the catalysis.

Multi-scale simulation reveals that an amino acid substitution increases photosensitizing reaction inputs in Rhodopsins.

Evaluating the availability of molecular oxygen (O2 ) and energy of excited states in the retinal binding site of rhodopsin is a crucial challenging first step to understand photosensitizing reactions in wild-type (WT) and mutant rhodopsins by absorbing visible light. In the present work, energies of the ground and excited states related to 11-cis-retinal and the O2 accessibility to the ß-ionone ring are evaluated inside WT and human M207R mutant rhodopsins. Putative O2 pathways within rhodopsins are identified by using molecular dynamics simulations, Voronoi-diagram analysis, and implicit ligand sampling while retinal energetic properties are investigated through density functional theory, and quantum mechanical/molecular mechanical methods. Here, the predictions reveal that an amino acid substitution can lead to enough energy and O2 accessibility in the core hosting retinal of mutant rhodopsins to favor the photosensitized singlet oxygen generation, which can be useful in understanding retinal degeneration mechanisms and in designing blue-lighting-absorbing proteic photosensitizers.

Reactivation of VX-Inhibited Human Acetylcholinesterase by Deprotonated Pralidoxime. A Complementary Quantum Mechanical Study.

In the present work, we performed a complementary quantum mechanical (QM) study to describe the mechanism by which deprotonated pralidoxime (2-PAM) could reactivate human (Homo sapiens sapiens) acetylcholinesterase (HssAChE) inhibited by the nerve agent VX. Such a reaction is proposed to occur in subsequent addition-elimination steps, starting with a nucleophile bimolecular substitution (SN2) mechanism through the formation of a trigonal bipyramidal transition state (TS). A near attack conformation (NAC), obtained in a former study using molecular mechanics (MM) calculations, was taken as a starting point for this project, where we described the possible formation of the TS. Together, this combined QM/MM study on AChE reactivation shows the feasibility of the reactivation occurring via attack of the deprotonated form of 2-PAM against the Ser203-VX adduct of HssAChE.