Unveiling the catalytic mechanism of GTP hydrolysis in microtubules.

Microtubules (MTs) are large cytoskeletal polymers, composed of αß-tubulin heterodimers, capable of stochastically converting from polymerizing to depolymerizing states and vice versa. Depolymerization is coupled with hydrolysis of guanosine triphosphate (GTP) within ß-tubulin. Hydrolysis is favored in the MT lattice compared to a free heterodimer with an experimentally observed rate increase of 500- to 700-fold, corresponding to an energetic barrier lowering of 3.8 to 4.0 kcal/mol. Mutagenesis studies have implicated α-tubulin residues, α:E254 and α:D251, as catalytic residues completing the ß-tubulin active site of the lower heterodimer in the MT lattice. The mechanism for GTP hydrolysis in the free heterodimer, however, is not understood. Additionally, there has been debate concerning whether the GTP-state lattice is expanded or compacted relative to the GDP state and whether a "compacted" GDP-state lattice is required for hydrolysis. In this work, extensive quantum mechanics/molecular mechanics simulations with transition-tempered metadynamics free-energy sampling of compacted and expanded interdimer complexes, as well as a free heterodimer, have been carried out to provide clear insight into the GTP hydrolysis mechanism. α:E254 was found to be the catalytic residue in a compacted lattice, while in the expanded lattice, disruption of a key salt bridge interaction renders α:E254 less effective. The simulations reveal a barrier decrease of 3.8 ± 0.5 kcal/mol for the compacted lattice compared to a free heterodimer, in good agreement with experimental kinetic measurements. Additionally, the expanded lattice barrier was found to be 6.3 ± 0.5 kcal/mol higher than compacted, demonstrating that GTP hydrolysis is variable with lattice state and slower at the MT tip.

Road-blocker HSP disease mutation disrupts pre-organization for ATP hydrolysis in kinesin through a second sphere control.

Kinesin motor proteins perform several essential cellular functions powered by the adenosine triphosphate (ATP) hydrolysis reaction. Several single-point mutations in the kinesin motor protein KIF5A have been implicated to hereditary spastic paraplegia disease (HSP), a lethal neurodegenerative disease in humans. In earlier studies, we have shown that a series of HSP-related mutations can impair the kinesin's long-distance displacement or processivity by modulating the order-disorder transition of the linker connecting the heads to the coiled coil. On the other hand, the reduction of kinesin's ATP hydrolysis reaction rate by a distal asparagine-to-serine mutation is also known to cause HSP disease. However, the molecular mechanism of the ATP hydrolysis reaction in kinesin by this distal mutation is still not fully understood. Using classical molecular dynamics simulations combined with quantum mechanics/molecular mechanics calculations, the pre-organization geometry required for optimal hydrolysis in kinesin motor bound to α/ß-tubulin is determined. This optimal geometry has only a single salt-bridge (of the possible two) between Arg203-Glu236, putting a reactive water molecule at a perfect position for hydrolysis. Such geometry is also needed to create the appropriate configuration for proton translocation during ATP hydrolysis. The distal asparagine-to-serine mutation is found to disrupt this optimal geometry. Therefore, the current study along with our previous one demonstrates how two different effects on kinesin dynamics (processivity and ATP hydrolysis), caused by a different set of genotypes, can give rise to the same phenotype leading to HSP disease.

Structures and mechanisms of actin ATP hydrolysis.

The major cytoskeleton protein actin undergoes cyclic transitions between the monomeric G-form and the filamentous F-form, which drive organelle transport and cell motility. This mechanical work is driven by the ATPase activity at the catalytic site in the F-form. For deeper understanding of the actin cellular functions, the reaction mechanism must be elucidated. Here, we show that a single actin molecule is trapped in the F-form by fragmin domain-1 binding and present their crystal structures in the ATP analog-, ADP-Pi-, and ADP-bound forms, at 1.15-Å resolutions. The G-to-F conformational transition shifts the side chains of Gln137 and His161, which relocate four water molecules including W1 (attacking water) and W2 (helping water) to facilitate the hydrolysis. By applying quantum mechanics/molecular mechanics calculations to the structures, we have revealed a consistent and comprehensive reaction path of ATP hydrolysis by the F-form actin. The reaction path consists of four steps: 1) W1 and W2 rotations; 2) PG-O3B bond cleavage; 3) four concomitant events: W1-PO3- formation, OH- and proton cleavage, nucleophilic attack by the OH- against PG, and the abstracted proton transfer; and 4) proton relocation that stabilizes the ADP-Pi-bound F-form actin. The mechanism explains the slow rate of ATP hydrolysis by actin and the irreversibility of the hydrolysis reaction. While the catalytic strategy of actin ATP hydrolysis is essentially the same as those of motor proteins like myosin, the process after the hydrolysis is distinct and discussed in terms of Pi release, F-form destabilization, and global conformational changes.

Electric fields control water-gated proton transfer in cytochrome c oxidase.

Aerobic life is powered by membrane-bound enzymes that catalyze the transfer of electrons to oxygen and protons across a biological membrane. Cytochrome c oxidase (CcO) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O2 reduction into proton pumping. Here we show that CcO creates orientated electric fields around a nonpolar cavity next to the active site, establishing a molecular switch that directs the protons along distinct pathways. By combining large-scale quantum chemical density functional theory (DFT) calculations with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations and atomistic molecular dynamics (MD) explorations, we find that reduction of the electron donor, heme a, leads to dissociation of an arginine (Arg438)-heme a3 D-propionate ion-pair. This ion-pair dissociation creates a strong electric field of up to 1 V Å-1 along a water-mediated proton array leading to a transient proton loading site (PLS) near the active site. Protonation of the PLS triggers the reduction of the active site, which in turn aligns the electric field vectors along a second, "chemical," proton pathway. We find a linear energy relationship of the proton transfer barrier with the electric field strength that explains the effectivity of the gating process. Our mechanism shows distinct similarities to principles also found in other energy-converting enzymes, suggesting that orientated electric fields generally control enzyme catalysis.

The adaptability of the ion-binding site by the Ag(I)/Cu(I) periplasmic chaperone SilF.

The periplasmic chaperone SilF has been identified as part of an Ag(I) detoxification system in Gram-negative bacteria. Sil proteins also bind Cu(I) but with reported weaker affinity, therefore leading to the designation of a specific detoxification system for Ag(I). Using isothermal titration calorimetry, we show that binding of both ions is not only tighter than previously thought but of very similar affinities. We investigated the structural origins of ion binding using molecular dynamics and QM/MM simulations underpinned by structural and biophysical experiments. The results of this analysis showed that the binding site adapts to accommodate either ion, with key interactions with the solvent in the case of Cu(I). The implications of this are that Gram-negative bacteria do not appear to have evolved a specific Ag(I) efflux system but take advantage of the existing Cu(I) detoxification system. Therefore, there are consequences for how we define a particular metal resistance mechanism and understand its evolution in the environment.

Noncovalent interactions that tune the reactivities of the flavins in bifurcating electron transferring flavoprotein.

Bifurcating electron transferring flavoproteins (Bf-ETFs) tune chemically identical flavins to two contrasting roles. To understand how, we used hybrid quantum mechanical molecular mechanical calculations to characterize noncovalent interactions applied to each flavin by the protein. Our computations replicated the differences between the reactivities of the flavins: the electron transferring flavin (ETflavin) was calculated to stabilize anionic semiquinone (ASQ) as needed to execute its single-electron transfers, whereas the Bf flavin (Bfflavin) was found to disfavor the ASQ state more than does free flavin and to be less susceptible to reduction. The stability of ETflavin ASQ was attributed in part to H-bond donation to the flavin O2 from a nearby His side chain, via comparison of models employing different tautomers of His. This H-bond between O2 and the ET site was uniquely strong in the ASQ state, whereas reduction of ETflavin to the anionic hydroquinone (AHQ) was associated with side chain reorientation, backbone displacement, and reorganization of its H-bond network including a Tyr from the other domain and subunit of the ETF. The Bf site was less responsive overall, but formation of the Bfflavin AHQ allowed a nearby Arg side chain to adopt an alternative rotamer that can H-bond to the Bfflavin O4. This would stabilize the anionic Bfflavin and rationalize effects of mutation at this position. Thus, our computations provide insights on states and conformations that have not been possible to characterize experimentally, offering explanations for observed residue conservation and raising possibilities that can now be tested.

Probing the range of applicability of structure- and energy-adjusted QM/MM link bonds III: QM/MM MD simulations of solid-state systems at the example of layered carbon structures.

The previously introduced workflow to achieve an energetically and structurally optimized description of frontier bonds in quantum mechanical/molecular mechanics (QM/MM)-type applications was extended into the regime of computational material sciences at the example of a layered carbon model systems. Optimized QM/MM link bond parameters at HSEsol/6-311G(d,p) and self-consistent density functional tight binding (SCC-DFTB) were derived for graphitic systems, enabling detailed investigation of specific structure motifs occurring in graphene-derived structures v i a quantum-chemical calculations. Exemplary molecular dynamics (MD) simulations in the isochoric-isothermic (NVT) ensemble were carried out to study the intercalation of lithium and the properties of the Stone-Thrower-Wales defect. The diffusivity of lithium as well as hydrogen and proton adsorption on a defective graphene surface served as additional example. The results of the QM/MM MD simulations provide detailed insight into the applicability of the employed link-bond strategy when studying intercalation and adsorption properties of graphitic materials.

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.

Iron-histidine bonding in bishistidyl hemoproteins-A local vibrational mode study.

We investigated the intrinsic strength of distal and proximal FeN bonds for both ferric and ferrous oxidation states of bishistidyl hemoproteins from bacteria, animals, human, and plants, including two cytoglobins, ten hemoglobins, two myoglobins, six neuroglobins, and six phytoglobins. As a qualified measure of bond strength, we used local vibrational force constants k a (FeN) based on local mode theory developed in our group. All calculations were performed with a hybrid QM/MM ansatz. Starting geometries were taken from available x-ray structures. k a (FeN) values were correlated with FeN bond lengths and covalent bond character. We also investigated the stiffness of the axial NFeN bond angle. Our results highlight that protein effects are sensitively reflected in k a (FeN), allowing one to compare trends in diverse protein groups. Moreover, k a (NFeN) is a perfect tool to monitor changes in the axial heme framework caused by different protein environments as well as different Fe oxidation states.

Modeling solvent effects and convergence of 31P-NMR shielding calculations with COBRAMM.

Solvent effects on 31P-NMR parameters for triphenylphosphine oxide and triphenylphosphine in chloroform have been extensively investigated by testing different solvation models. The solvent is described implicitly, mixed implicitly/explicitly, and using full explicit models. Polarizable continuum model (PCM), molecular dynamic (MD) simulations, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations are used to disclose the effects of solute/solvent interactions and, more generally, the role of the embedding in NMR simulations. The results show the beneficial effect of carrying out QM/MM optimizations on top of geometries directly extracted from classical MD simulations, used to ensure representative conformational sampling. The nuclear shielding convergence has been tested against a different number of snapshots and with the inclusion of solvent shells into the QM region. An automated MD//QM/MM//GIAO protocol, implemented in the COBRAMM package, is here proposed and tested on trimethyl phosphite showing that our approach boosts the convergence of nuclear shielding satisfactorily. The present work aims to be a stepping-stone to assess proper QM/MM computational strategies in simulating chemical shifts in non-homogeneous systems like supramolecular and biological systems.

Iron-Sulfur Peptides Mimicking Ferredoxin for an Efficient Electron Transfer to Hydrogenase.

In the green alga Chlamydomonas reinhardtii, hydrogenase HydA1 converts protons and electrons to H2 at the H-cluster, which includes a [4Fe-4S] cluster linked to a [2Fe] cluster. The yield of H2 is limited by the electron transfer to HydA1, mediated by the iron-sulfur unit of a photosynthetic electron transfer ferredoxin (PetF). In this study, I have investigated by molecular dynamics and the hybrid quantum mechanics/molecular mechanics method two canonical iron-sulfur peptides (PM1 and FBM) that hold potential as PetF replacements. Using a docking approach, I predict that the distance between the two iron-sulfur clusters in FBM/HydA1 is shorter than in PM1/HydA1, ensuring a greater electron transfer rate. This finding is in line with the reported higher H2 production rates for FBM/HydA1. I also show that the redox potential of these peptides, and therefore their electron transfer properties, can be changed by single-residue mutations in the secondary coordination sphere of their cluster. In particular, I have designed a PM1 variant that disrupts the hydrogen-bonding network between water and the cluster, shifting the redox potential negatively compared to PM1. These results will guide experiments aimed at replacing PetF with peptides that can unlock the biotechnological potential of the alga.

QM/MM Study Into the Mechanism of Oxidative C=C Double Bond Cleavage by Lignostilbene-α,ß-Dioxygenase.

The enzymatic biosynthesis of fragrance molecules from lignin fragments is an important reaction in biotechnology for the sustainable production of fine chemicals. In this work we investigated the biosynthesis of vanillin from lignostilbene by a nonheme iron dioxygenase using QM/MM and tested several suggested proposals via either an epoxide or dioxetane intermediate. Binding of dioxygen to the active site of the protein results in the formation of an iron(II)-superoxo species with lignostilbene cation radical. The dioxygenase mechanism starts with electrophilic attack of the terminal oxygen atom of the superoxo group on the central C=C bond of lignostilbene, and the second-coordination sphere effects in the substrate binding pocket guide the reaction towards dioxetane formation. The computed mechanism is rationalized with thermochemical cycles and valence bond schemes that explain the electron transfer processes during the reaction mechanism. Particularly, the polarity of the protein and the local electric field and dipole moments enable a facile electron transfer and an exergonic dioxetane formation pathway.

Exploring the Esterase Catalytic Activity of Minimalist Heptapeptide Amyloid Fibers.

This paper investigates the esterase activity of minimalist amyloid fibers composed of short seven-residue peptides, IHIHIHI (IH7) and IHIHIQI (IH7Q), with a particular focus on the role of the sixth residue position within the peptide sequence. Through computational simulations and analyses, we explore the molecular mechanisms underlying catalysis in these amyloid-based enzymes. Contrary to initial hypotheses, our study reveals that the twist angle of the fiber, and thus the catalytic site's environment, is not notably affected by the sixth residue. Instead, the sixth residue interacts with the p-nitrophenylacetate (pNPA) substrate, particularly through its -NO2 group, potentially enhancing catalysis. Quantum mechanics/molecular mechanics (QM/MM) simulations of the reaction mechanism suggest that the polarizing effect of glutamine enhances catalytic activity by forming a stabilizing network of hydrogen bonds with pNPA, leading to lower energy barriers and a more exergonic reaction. Our findings provide valuable insights into the intricate interplay between peptide sequence, structural arrangement, and catalytic function in amyloid-based enzymes, offering potentially valuable information for the design and optimization of biomimetic catalysts.

Bacteria at Work - Experimental and Theoretical Studies Reveal the Catalytic Mechanism of Ectoine Synthase.

Ectoine synthase (EctC) catalyses the ultimate step of ectoine biosynthesis, a kosmotropic compound produced as compatible solute by many bacteria and some archaea or eukaryotes. EctC is an Fe2+-dependent homodimeric cytoplasmic protein. Using Mössbauer spectroscopy, molecular dynamics simulations and QM/MM calculations, we determined the most likely coordination number and geometry of the Fe2+ ion and proposed a mechanism of the EctC-catalysed reaction. Most notably, we show that apart from the three amino acids binding to the iron ion (Glu57, Tyr84 and His92), one water molecule and one hydroxide ion are required as additional ligands for the reaction to occur. They fill the first coordination sphere of the Fe2+-cofactor and act as critical proton donors and acceptors during the cyclization reaction.

Resolving the Structure of a Guanine Quadruplex in TMPRSS2 Messenger RNA by Circular Dichroism and Molecular Modeling.

The presence of a guanine quadruplex in the opening reading frame of the messenger RNA coding for the transmembrane serine protease 2 (TMPRSS2) may pave the way to original anticancer and host-oriented antiviral strategy. Indeed, TMPRSS2 in addition to being overexpressed in different cancer types, is also related to the infection of respiratory viruses, including SARS-CoV-2, by promoting the cellular and viral membrane fusion through its proteolytic activity. The design of selective ligands targeting TMPRSS2 messenger RNA requires a detailed knowledge, at atomic level, of its structure. Therefore, we have used an original experimental-computational protocol to predict the first resolved structure of the parallel guanine quadruplex secondary structure in the RNA of TMPRSS2, which shows a rigid core flanked by a flexible loop. This represents the first atomic scale structure of the guanine quadruplex structure present in TMPRSS2 messenger RNA.

Electron Attachment to the Nucleobase Uracil in Diethylene Glycol: The Signature of a Doorway.

The cellular environment plays a significant role in low energy electron-mediated radiation damage to genetic materials. In this study, we have modeled the effect of the bulk medium on electron attachment to nucleobases in diethylene glycol (DEG) using uracil as a test case, in accordance with recent experimental work on the observation of dissociative quasi-free electron attachment to nucleoside via excited anion radical in solution (in DEG). Our EOM-CCSD-based quantum mechanical/molecular mechanical (QM/MM) simulations indicate that the electron scavenging by uracil in DEG is much slower than that observed in the aqueous medium due to its viscosity. This work also establishes that a doorway mechanism exists in uracil microsolvated and bulk solvated with DEG, with the dipole-bound state and solvent-bound state acting as doorway states, respectively.

Molecular Mechanism of Labelling Functional Cysteines by Heterocyclic Thiones.

Heterocyclic thiones have recently been identified as reversible covalent warheads, consistent with their mild electrophilic nature. Little is known so far about their mechanism of action in labelling nucleophilic sidechains, especially cysteines. The vast number of tractable cysteines promotes a wide range of target proteins to examine; however, our focus was put on functional cysteines. We chose the main protease of SARS-CoV-2 harboring Cys145 at the active site that is a structurally characterized and clinically validated target of covalent inhibitors. We screened an in-house, cysteine-targeting covalent inhibitor library which resulted in several covalent fragment hits with benzoxazole, benzothiazole and benzimidazole cores. Thione derivatives and Michael acceptors were selected for further investigations with the objective of exploring the mechanism of inhibition of the thiones and using the thoroughly characterized Michael acceptors for benchmarking our studies. Classical and hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations were carried out that revealed a new mechanism of covalent cysteine labelling by thione derivatives, which was supported by QM and free energy calculations and by a wide range of experimental results. Our study shows that the molecular recognition step plays a crucial role in the overall binding of both sets of molecules.

In-silico predicted mouse melanopsins with blue spectral shifts deliver efficient subcellular signaling.

Melanopsin is a photopigment belonging to the G Protein-Coupled Receptor (GPCR) family expressed in a subset of intrinsically photosensitive retinal ganglion cells (ipRGCs) and responsible for a variety of processes. The bistability and, thus, the possibility to function under low retinal availability would make melanopsin a powerful optogenetic tool. Here, we aim to utilize mouse melanopsin to trigger macrophage migration by its subcellular optical activation with localized blue light, while simultaneously imaging the migration with red light. To reduce melanopsin's red light sensitivity, we employ a combination of in silico structure prediction and automated quantum mechanics/molecular mechanics modeling to predict minimally invasive mutations to shift its absorption spectrum towards the shorter wavelength region of the visible spectrum without compromising the signaling efficiency. The results demonstrate that it is possible to achieve melanopsin mutants that resist red light-induced activation but are activated by blue light and display properties indicating preserved bistability. Using the A333T mutant, we show that the blue light-induced subcellular melanopsin activation triggers localized PIP3 generation and macrophage migration, which we imaged using red light, demonstrating the optogenetic utility of minimally engineered melanopsins.

Merocyanines form bacteriorhodopsins with strongly bathochromic absorption maxima.

There is a need to shift the absorbance of biomolecules to the optical transparency window of tissue for applications in optogenetics and photo-pharmacology. There are a few strategies to achieve the so-called red shift of the absorption maxima. Herein, a series of 11 merocyanine dyes were synthesized and employed as chromophores in place of retinal in bacteriorhodopsin (bR) to achieve a bathochromic shift of the absorption maxima relative to bR's [Formula: see text] of 568 nm. Assembly with the apoprotein bacterioopsin (bO) led to stable, covalently bound chromoproteins with strongly bathochromic absorbance bands, except for three compounds. Maximal red shifts were observed for molecules 9, 2, and 8 in bR where the [Formula: see text] was 766, 755, and 736 nm, respectively. While these three merocyanines have different end groups, they share a similar structural feature, namely, a methyl group which is located at the retinal equivalent position 13 of the polyene chain. The absorption and fluorescence data are also presented for the retinal derivatives in their aldehyde, Schiff base (SB), and protonated SB (PSB) forms in solution. According to their hemicyanine character, the PSBs and their analogue bRs exhibited fluorescence quantum yields (Φf) several orders of magnitude greater than native bR (Φf 0.02 to 0.18 versus 1.5 × 10-5 in bR) while also exhibiting much smaller Stokes shifts than bR (400 to 1000 cm-1 versus 4030 cm-1 in bR). The experimental results are complemented by quantum chemical calculations where excellent agreement between the experimental [Formula: see text] and the calculated [Formula: see text] was achieved with the second-order algebraic-diagrammatic construction [ADC(2)] method. In addition, quantum mechanics/molecular mechanics (QM/MM) calculations were employed to shed light on the origin of the bathochromic shift of merocyanine 2 in bR compared with native bR.

Rational prioritization strategy allows the design of macrolide derivatives that overcome antibiotic resistance.

Antibiotic resistance is a major threat to global health; this problem can be addressed by the development of new antibacterial agents to keep pace with the evolutionary adaptation of pathogens. Computational approaches are essential tools to this end since their application enables fast and early strategical decisions in the drug development process. We present a rational design approach, in which acylide antibiotics were screened based on computational predictions of solubility, membrane permeability, and binding affinity toward the ribosome. To assess our design strategy, we tested all candidates for in vitro inhibitory activity and then evaluated them in vivo with several antibiotic-resistant strains to determine minimal inhibitory concentrations. The predicted best candidate is synthetically more accessible, exhibits higher solubility and binding affinity to the ribosome, and is up to 56 times more active against resistant pathogens than telithromycin. Notably, the best compounds designed by us show activity, especially when combined with the membrane-weakening drug colistin, against Acinetobacter baumanii, Pseudomonas aeruginosa, and Escherichia coli, which are the three most critical targets from the priority list of pathogens of the World Health Organization.