Understanding the Differences of Danusertib's Residence Time in Aurora Kinases A/B: Dissociation Paths and Key Residues Identified using Conventional and Enhanced Molecular Dynamics Simulations.

The accurate experimental estimation of protein-ligand systems' residence time (τ) has become very relevant in drug design projects due to its importance in the last stages of refinement of the drug's pharmacodynamics and pharmacokinetics. It is now well-known that it is not sufficient to estimate the affinity of a protein-drug complex in the thermodynamic equilibrium process in in vitro experiments (closed systems), where the concentrations of the drug and protein remain constant. On the contrary, it is mandatory to consider the conformational dynamics of the system in terms of the binding and unbinding processes between protein and drugs in in vivo experiments (open systems), where their concentrations are in constant flux. This last model has been proven to dictate much of several drugs' pharmacological activities in vivo. At the atomistic level, molecular dynamics simulations can explain why some drugs are more effective than others or unveil the molecular aspects that make some drugs work better in one molecular target. Here, the protein kinases Aurora A/B, complexed with its inhibitor Danusertib, were studied using conventional and enhanced molecular dynamics (MD) simulations to estimate the dissociation paths and, therefore, the computational τ values and their comparison with experimental ones. Using classical molecular dynamics (cMD), three differential residues within the Aurora A/B active site, which seems to play an essential role in the observed experimental Danusertib's residence time against these kinases, were characterized. Then, using WT-MetaD, the relative Danusertib's residence times against Aurora A/B kinases were measured in a nanosecond time scale and were compared to those τ values observed experimentally. In addition, the potential dissociation paths of Danusertib in Aurora A and B were characterized, and differences that might be explained by the differential residues in the enzyme's active sites were found. In perspective, it is expected that this computational protocol can be applied to other protein-ligand complexes to understand, at the molecular level, the differences in residence times and amino acids that may contribute to it.

Understanding the Interactions of Ubiquitin-Specific Protease 7 with Its Substrates through Molecular Dynamics Simulations: Insights into the Role of Its C-Terminal Domains in Substrate Recognition.

Ubiquitin-specific protease 7 (USP7) is a deubiquitinase enzyme that plays a critical role in regulating various cellular processes by cleaving ubiquitin molecules from target proteins. The C-terminal loop (CTL) motif is a specific region at the C-terminal end of the USP7 enzyme. Recent experiments suggest that the CTL motif plays a role in USP7's catalytic activity by contributing to the enzyme's structural stability, substrate recognition, and catalytic efficiency. The objective of this work is to elucidate these roles through the utilization of computational methods for molecular simulations. For this, we conducted extensive molecular dynamics (MD) simulations to investigate the conformational dynamics and protein-protein interactions within the USP7 enzyme-substrate complex with the substrate consisting of the ubiquitin tagged with the fluorescent label rhodamine 110-gly (Ub-Rho). Our results demonstrate that the CTL motif plays a crucial role in stabilizing the Ubl domains' conformation and augmenting the stability of active conformations within the enzyme-substrate complex. Conversely, the absence of the CTL motif results in increased flexibility and variability in Ubl domains' motion, leading to a reduced percentage of active conformations. Furthermore, our analysis of protein-protein interactions highlights the significance of the CTL motif in anchoring the Ubl45 domains to the catalytic domain (CD), thereby facilitating stable interactions with the substrate. Overall, our findings provide valuable insights into the conformational dynamics and protein-protein interactions inherent in the USP7 enzyme-substrate complex. These insights shed light on some mechanistic details of USP7 concerning the substrate's recognition before its catalytic action.

On the mechanical response of graphene-capped copper nanoparticles.

In this study, we investigated the mechanical behavior of pristine copper (Cu) nanoparticles (NPs) and Cu@graphene (Cu@G) hybrid NPs using molecular dynamics simulations. The longitudinal engineering strain was calculated as a measure of compression until reaching 25% of the initial size of the NPs. The stress-strain curves revealed the elastic-to-plastic transition in the Cu NPs at a longitudinal strain of 3.57% with a yield strength of 6.15 GPa. On the other hand, the Cu@G NPs exhibited a maximum average load point at a longitudinal strain of 6.81% with a yield strength of 8.26 GPa. The hybrid Cu@G NPs showed increased strength and resistance to plastic deformation compared to the pure Cu NPs, while the calculation of the elastic modulus indicated a higher load resistance provided by the graphene coverage for the Cu@G NPs. Furthermore, the analysis of atomic configurations, dislocations, and stress distribution demonstrated that the graphene flakes play a crucial role in preventing dislocation events and faceting in the Cu@G NPs by acting as a shock absorber, distributing the applied force on themselves, and producing a more homogeneous stress distribution on the Cu NPs; additionally, they prevent the movement of Cu atoms, reducing the occurrence of dislocations and surface faceting, thanks to their supportive effect. Overall, our findings highlight the potential of hybrid nanomaterials, such as Cu@G, for enhancing the mechanical properties of metallic NPs, which could have significant implications for the development of advanced nanomaterials with improved performance in a variety of applications.

Effect of strand register in the stability and reactivity of crystals from peptides forming amyloid fibrils.

Amyloids are cytotoxic protein aggregates that deposit in human tissues, leading to several health disorders. Their aggregates can also exhibit catalytic properties, and they have been used as candidates for the development of functional biomaterials. Despite being polymorphic, amyloids often assemble as cross-ß fibrils formed by in-register ß sheet layers. Recent studies of some amyloidogenic protein segments revealed that they crystallize as antiparallel out-of-register ß sheets. Such arrangement has been proposed to be responsible for the cytotoxicity in amyloid diseases, however, there is still no consensus on the molecular mechanism. Interestingly, two amyloidogenic peptide segments, NFGAILS and FGAILSS, arrange into out-of-register and in-register ß sheets, respectively, even though they solely differ by one aminoacid residue at both termini. In this work, we used density functional theory (DFT) to address how the strand register contributes into the packing and molecular properties of the NFGAILS and FGAILSS crystals. Our results show that the out-of-register structure is substantially more stable, at 0 K, than the in-register one due to stronger inter-strand contacts. Based on an analysis of the electrostatic potential of the crystal slabs, it is suggested that the out-of-register may potentially interact with negatively charged groups, like those found in cell membranes. Moreover, calculated reactivity descriptors indicate a similar outcome, where only the out-of-register peptide exhibits intrinsic reactive surface sites at the exposed amine and carboxylic groups. It is therefore suggested that the out-of-register arrangement may indeed be crucial for amyloid cytotoxicity. The findings presented here could help to further our understanding of amyloid aggregation, function, and toxicity.

Can graphene improve the thermal conductivity of copper nanofluids?

Copper (Cu) nanofluids (NFs) have attracted attention due to their high thermal conductivity, which has conferred a wide variety of applications. However, their high reactivity favors oxidation, corrosion and aggregation, leading them to lose their properties of interest. Copper capped by graphene (Cu@G) core@shell nanoparticles (NPs) have also attracted interest from the medical and industrial sectors because graphene can shield the Cu NPs from undesired phenomena. Additionally, they share some properties that expand the range of applications of Cu NFs. In this work, new Morse potentials are reported to reproduce the behavior of Cu@G NPs through molecular dynamics. Coordination-dependent Morse parameters were fitted for C, H, and Cu based on density functional theory calculations. Then, these parameters were implemented to evaluate the thermal conductivity of Cu@G NFs employing the Green-Kubo formalism, with NPs from 1.5 to 6.1 nm at 100 to 800 K, varying the size, the number of layers and the orientation of the graphene flakes. It was found that Cu@G NFs are stable and have an improved thermal conductivity compared to the Cu NFs, being 3.7 to 18.2 times higher at 300 K with only one graphene layer and above 26.2 times higher for the graphene-trilayered NPs. These values can be higher for temperatures below 300 K. Oppositely, the size, homogeneity and orientations of the graphene flakes did not affect the thermal conductivity of the Cu@G NFs.

Structural Insights into the Inhibition Site in the Phosphorylcholine Phosphatase Enzyme of Pseudomonas aeruginosa.

Pseudomonas aeruginosa is a highly pathogenic Gram-negative microorganism associated with high mortality levels in burned or immunosuppressed patients or individuals affected by cystic fibrosis. Studies support a colonization mechanism whereby P. aeruginosa can breakdown the host cell membrane phospholipids through the sequential action of two enzymes: (I) hemolytic phospholipase C acting upon phosphatidylcholine or sphingomyelin to produce phosphorylcholine (Pcho) and (II) phosphorylcholine phosphatase (PchP) that hydrolyzes Pcho to generate choline and inorganic phosphate. This coordinated action provides the bacteria with carbon, nitrogen, and inorganic phosphate to support growth. Furthermore, PchP exhibits a distinctive inhibition mechanism by high substrate concentration. Here, we combine kinetic assays and computational approaches such as molecular docking, molecular dynamics, and free-energy calculations to describe the inhibitory site of PchP, which shares specific residues with the enzyme's active site. Our study provides insights into a coupled inhibition mechanism by the substrate, allowing us to postulate that the integrity of the inhibition site is needed to the correct functioning of the active site. Our results allow us to gain a better understanding of PchP function and provide the basis for a rational drug design that might contribute to the treatment of infections caused by this important opportunistic pathogen.

Crystal structure and molecular dynamics simulations of a promiscuous ancestor reveal residues and an epistatic interaction involved in substrate binding and catalysis in the ATP-dependent vitamin kinase family members.

Enzymes with hydroxymethylpyrimidine/phosphomethylpyrimidine kinase activity (HMPPK) are essential in the vitamin B1 (thiamine pyrophosphate) biosynthesis and recycling pathways. In contrast, enzymes with pyridoxal kinase activity (PLK) produce pyridoxal phosphate (vitamin B6), an essential cofactor for various biochemical reactions. In the ATP-dependent vitamin kinases family, the members of PLK/HMPPK-like subfamily have both enzymatic activities. It has been proposed that the promiscuous PLK activity of ancestral HMPPK enzymes could have been the starting point for this activity. In earlier work, we reconstructed the ancestral sequences of this family and characterized the substrate specificity of the common ancestor between PLK/HMPPK-like and HMPPK enzymes (AncC). From these studies, the Gln45Met mutation was proposed as a critical event for the PLK activity emergence. Here, we crystallize and determine the AncC structure by X-ray crystallography and assess the role of the Gln45Met mutation by site-directed mutagenesis. Kinetic characterization of this mutant shows a significant increase in the PL affinity. Through molecular dynamics simulation and MM/PBSA calculations some residues, important for substrate interactions and catalysis, were identified in the wild type and in the mutated ancestor. Interestingly, a strong epistatic interaction responsible for the evolutionary pathway of the PLK activity in PLK/HMPPK-like enzymes was revealed. Also, other putative mutations relevant to PLK activity in modern PLK/HMPPK-like enzymes were identified.

Coarse-Grained Parameters for Divalent Cations within the SIRAH Force Field.

Although molecular dynamics simulations allow for the study of interactions among virtually all biomolecular entities, metal ions still pose significant challenges in achieving an accurate structural and dynamical description of many biological assemblies, particularly to coarse-grained (CG) models. Although the reduced computational cost of CG methods often makes them the technique of choice for the study of large biomolecular systems, the parameterization of metal ions is still very crude or not available for the vast majority of CG force fields. Here, we show that incorporating statistical data retrieved from the Protein Data Bank (PDB) to set specific Lennard-Jones interactions can produce structurally accurate CG molecular dynamics simulations using the SIRAH force field. We provide a set of interaction parameters for calcium, magnesium, and zinc ions, which cover more than 80% of the metal-bound structures reported in the PDB. Simulations performed on several proteins and DNA systems show that it is possible to preclude the use of topological constraints by modifying specific Lennard-Jones interactions.

Characterization of hydroxymethylpyrimidine phosphate kinase from mesophilic and thermophilic bacteria and structural insights into their differential thermal stability.

The hydroxymethylpyrimidine phosphate kinases (HMPPK) encoded by the thiD gene are involved in the thiamine biosynthesis pathway, can perform two consecutive phosphorylations of 4-amino-5-hydroxymethyl-2-methyl pyrimidine (HMP) and are found in thermophilic and mesophilic bacteria, but only a few characterizations of mesophilic enzymes are available. The presence of another homolog enzyme (pyridoxal kinase) that can only catalyze the first phosphorylation of HMP and encoded by pdxK gene, has hampered a precise annotation in this enzyme family. Here we report the kinetic characterization of two HMPPK with structure available, the mesophilic and thermophilic enzyme from Salmonella typhimurium (StHMPPK) and Thermus thermophilus (TtHMPPK), respectively. Also, given their high structural similarity, we have analyzed the structural determinants of protein thermal stability in these enzymes by molecular dynamics simulation. The results show that pyridoxal kinases (PLK) from gram-positive bacteria (PLK/HMPPK-like enzymes) constitute a phylogenetically separate group from the canonical PLK, but closely related to the HMPPK, so the PLK/HMPPK-like and canonical PLK, both encoded by pdxK genes, are different and must be annotated distinctly. The kinetic characterization of StHMPPK and TtHMPPK, shows that they perform double phosphorylation on HMP, both enzymes are specific for HMP, not using pyridoxal-like molecules as substrates and their kinetic mechanism involves the formation of a ternary complex. Molecular dynamics simulation shows that StHMPPK and TtHMPPK have striking differences in their conformational flexibility, which can be correlated with the hydrophobic packing and electrostatic interaction network given mainly by salt bridge bonds, but interestingly not by the number of hydrogen bond interactions as reported for other thermophilic enzymes. ENZYMES: EC 2.7.1.49, EC 2.7.4.7, EC 2.7.1.35, EC 2.7.1.50.

Molecular Insights into the Trapping Effect of Ca2+ in Protein Kinase A: A Molecular Dynamics Study.

Protein kinase A has become a model system for the study of kinases, and therefore, a comprehensive understanding of the underlying molecular mechanisms in its catalytic cycle is of crucial importance. One of the aspects that has received recent attention is the role that metal cofactors play in the catalytic cycle. Although Mg2+ is the well-known physiological ion used by protein kinases, Ca2+ ions can also assist the phosphoryl transfer reaction but with lower catalytic activities. This inhibitory effect has been attributed to the ability of Ca2+ to trap the reaction products at the active site, and it has been proposed as a possible regulatory mechanism of the enzyme. Thus, in order to get a clearer understanding of these molecular events, computational simulations in the product state of PKA, in the presence of Mg2+ and Ca2+ ions, were performed through molecular dynamics (MD). Different protonation states of the active site were considered in order to model the different mechanistic pathways that have been proposed. Our results show that different protonation states of the phosphorylated serine residue at the peptide substrate (pSer21), as well as the protonation state of residue Asp166, can have a marked influence on the flexibility of regions surrounding the active site. This is the case of the glycine-rich loop, a structural motif that is directly involved in the release of the products from the PKA active site. MD simulations were capable to reproduce the crystallographic conformations but also showed other conformations not previously reported in the crystal structures that may be involved in enhancing the affinity of pSP20 to PKA in the presence of Ca2+. Hydrogen bonding interactions at the PKA-pSP20 interface were influenced whether by the protonation state of the active site or by the metal cofactor used by the enzyme. Altogether, our results provide molecular aspects into the inhibitory mechanism of Ca2+ in PKA and suggest which is the most probable protonation state of the phosphorylated product at the active site.

Tetrahydroquinoline-Isoxazole/Isoxazoline Hybrid Compounds as Potential Cholinesterases Inhibitors: Synthesis, Enzyme Inhibition Assays, and Molecular Modeling Studies.

A series of 44 hybrid compounds that included in their structure tetrahydroquinoline (THQ) and isoxazole/isoxazoline moieties were synthesized through the 1,3-dipolar cycloaddition reaction (1,3-DC) from the corresponding N-allyl/propargyl THQs, previously obtained via cationic Povarov reaction. In vitro cholinergic enzymes inhibition potential of all compounds was tested. Enzyme inhibition assays showed that some hybrids exhibited significant potency to inhibit acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). Especially, the hybrid compound 5n presented the more effective inhibition against AChE (4.24 µM) with an acceptable selectivity index versus BChE (SI: 5.19), while compound 6aa exhibited the greatest inhibition activity on BChE (3.97 µM) and a significant selectivity index against AChE (SI: 0.04). Kinetic studies were carried out for compounds with greater inhibitory activity of cholinesterases. Structure-activity relationships of the molecular hybrids were analyzed, through computational models using a molecular cross-docking algorithm and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) binding free energy approach, which indicated a good correlation between the experimental inhibition values and the predicted free binding energy.

Mechanistic insights into the phosphoryl transfer reaction in cyclin-dependent kinase 2: A QM/MM study.

Cyclin-dependent kinase 2 (CDK2) is an important member of the CDK family exerting its most important function in the regulation of the cell cycle. It catalyzes the transfer of the gamma phosphate group from an ATP (adenosine triphosphate) molecule to a Serine/Threonine residue of a peptide substrate. Due to the importance of this enzyme, and protein kinases in general, a detailed understanding of the reaction mechanism is desired. Thus, in this work the phosphoryl transfer reaction catalyzed by CDK2 was revisited and studied by means of hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. Our results suggest that the base-assisted mechanism is preferred over the substrate-assisted pathway when one Mg2+ is present in the active site, in agreement with a previous theoretical study. The base-assisted mechanism resulted to be dissociative, with a potential energy barrier of 14.3 kcal/mol, very close to the experimental derived value. An interesting feature of the mechanism is the proton transfer from Lys129 to the phosphoryl group at the second transition state, event that could be helping in neutralizing the charge on the phosphoryl group upon the absence of a second Mg2+ ion. Furthermore, important insights into the mechanisms in terms of bond order and charge analysis were provided. These descriptors helped to characterize the synchronicity of bond forming and breaking events, and to characterize charge transfer effects. Local interactions at the active site are key to modulate the charge distribution on the phosphoryl group and therefore alter its reactivity.

Binding and Catalytic Mechanisms of Veratryl Alcohol Oxidation by Lignin Peroxidase: A Theoretical and Experimental Study.

Lignin peroxidase (LiP) and its natural substrate veratryl alcohol (VA) play a crucial role in lignin degradation by white-rot fungi. Understanding the molecular determinants for the interaction of this enzyme with its substrates is essential in the rational design of engineered peroxidases for biotechnological application. Here, we combine computational and experimental approaches to analyze the interaction of Phanerochaete chrysosporium LiP (isoenzyme H8) with VA and its radical cation (VA•+, resulting from substrate oxidation by the enzyme). Interaction energy calculations at semiempirical quantum mechanical level (SQM) between LiP and VA/VA•+ enabled to identify those residues at the acidic environment of catalytic Trp171 involved in the main interactions. Then, a battery of variants, with single and multiple mutations at these residues (Glu168, Asp165, Glu250, Asp264, and Phe267), was generated by directed mutagenesis, and their kinetics parameters were estimated on VA and two additional substrates. The experimental results show that Glu168 and Glu250 are crucial for the binding of VA, with Glu250 also contributing to the turnover of the enzyme. The experimental results were further rationalized through new calculations of interaction energies between VA/VA•+ and LiP with each of the single mutations. Finally, the delocalization of spin density was determined with quantum mechanics/molecular mechanics calculations (QM/MM), further supporting the contribution of Glu250 to VA oxidation at Trp171.

Experimental and Theoretical Approaches in the Study of Phenanthroline-Tetrahydroquinolines for Alzheimer's Disease.

The imino-Diels-Alder reaction is one of the most common strategies in organic chemistry and is an important tool for providing a broad spectrum of biologically active heterocyclic systems. A combined theoretical and experimental study of the imino-Diels-Alder reaction is described. The new phenanthroline-tetrahydroquinolines were evaluated as cholinesterase inhibitors. Their cytotoxicity in human neuroblastoma SH-SY5Y cells was also evaluated. The theoretical results suggest that compounds formation in stages can be explained by endo cycloadducts under the established reaction conditions, thereby confirming experimental results obtained for percentage yield. These results allowed us to establish that pyridine substituent remarkably influences activation energy and reaction yield, as well as in acetylcholinesterase (AChE) activity. Among these derivatives, compounds with 4-pyridyl and 4-nitrophenyl showed favorable AChE activity and proved to be non-cytotoxic.

Studying the phosphoryl transfer mechanism of the E. coli phosphofructokinase-2: from X-ray structure to quantum mechanics/molecular mechanics simulations.

Phosphofructokinases (Pfks) catalyze the ATP-dependent phosphorylation of fructose-6-phosphate (F6P) and they are regulated in a wide variety of organisms. Although numerous aspects of the kinetics and regulation have been characterized for Pfks, the knowledge about the mechanism of the phosphoryl transfer reaction and the transition state lags behind. In this work, we describe the X-ray crystal structure of the homodimeric Pfk-2 from E. coli, which contains products in one site and reactants in the other, as well as an additional ATP molecule in the inhibitory allosteric site adjacent to the reactants. This complex was previously predicted when studying the kinetic mechanism of ATP inhibition. After removing the allosteric ATP, molecular dynamic (MD) simulations revealed conformational changes related to domain packing, as well as stable interactions of Lys27 and Asp256 with donor (ATP) and acceptor (fructose-6-) groups, and of Asp166 with Mg2+. The phosphoryl transfer reaction mechanism catalyzed by Pfk-2 was investigated through Quantum Mechanics/Molecular Mechanics (QM/MM) simulations using a combination of the string method and a path-collective variable for the exploration of its free energy surface. The calculated activation free energies showed that a dissociative mechanism, occurring with a metaphosphate intermediate formation followed by a proton transfer to Asp256, is more favorable than an associative one. The structural analysis reveals the role of Asp256 acting as a catalytic base and Lys27 stabilizing the transition state of the dissociative mechanism.

Classical and Quantum Mechanical Calculations of the Stacking Interaction of NdIII Complexes with Regular and Mismatched DNA Sequences.

The design of organometallic complexes used as selective intercalators to bind and react at DNA mismatch sites has concentrated efforts in the last few years. In this context, lanthanides have received attention to be employed as active optical centers due to their spectroscopic properties. Despite the fact that there are several experimental data about synthesis and DNA binding of these compounds, theoretical analyses describing their interaction with DNA are scarce. To understand the binding to regular and mismatched DNA sequences as well as to determine the effect of the intercalation on the spectroscopic properties of the complexes, a complete theoretical study going from classical to relativistic quantum mechanics calculations has been performed on some lanthanide complexes with phenanthroline derivatives synthesized and characterized herein, viz. [Nd(NO3)3(H2O)(dppz-R)] with R = H, NO2-, CN- and their [Nd(NO3)3(H2O)(dpq)] analogue, which was computationally modeled. The results were in correct agreement with the available experimental data showing that dppz complexes have higher binding affinities to DNA than dpq one and supporting the idea that these complexes are not selective to mismatch sites in the sampled time scale. Finally, the spectroscopic analysis evidence an intercalative binding mode and made possible the elucidation of the emission mechanism of these systems. This approach is proposed as a benchmark study to extend this methodology on similar systems and constitutes the first theoretical insight in the interaction between DNA and lanthanide complexes.

Inhibition of key enzymes in the inflammatory pathway by hybrid molecules of terpenes and synthetic drugs: In vitro and in silico studies.

The aim of this work was to compare the anti-inflammatory activity of compounds prepared from terpenes and the synthetic drugs ibuprofen and naproxen. The anti-inflammatory activity of the hybrid compounds was compared with the activity of the parent compounds. This was accomplished using in vitro inhibition of lipoxygenases (LOX) and COX-2, and in silico docking studies in 15-LOX and COX-2. The synthesized hybrids showed an inhibition of COX-2 and LOX between 9.8%-57.4% and 0.0%-97.7%, respectively. None of the hybrids showed an improvement in the inhibitory effect toward these pro-inflammatory enzymes, compared to the parent terpenes and non-steroidal anti-inflammatory drugs. The docking studies allowed us to predict the potential binding modes of hybrids 6-15 within COX-2 and 15-LOX active sites. The relative affinity of the compounds inside the binding sites could be explained by forming non-covalent interactions with most important and known amino acids reported for those enzymes. A good correlation (r2  = 0.745) between docking energies and inhibition percentages against COX-2 was found. The high inhibition obtained for compound 10 against COX-2 was explained by hydrogen bond interactions at the enzyme binding site. New synthetic possibilities could be obtained from our in silico models, improving the potency of these hybrid compounds.

Folding and Lipid Composition Determine Membrane Interaction of the Disordered Protein COR15A.

Plants from temperate climates, such as the model plant Arabidopsis thaliana, are challenged with seasonal low temperatures that lead to increased freezing tolerance in fall in a process termed cold acclimation. Among other adaptations, this involves the accumulation of cold-regulated (COR) proteins, such as the intrinsically disordered chloroplast-localized protein COR15A. Together with its close homolog COR15B, it stabilizes chloroplast membranes during freezing. COR15A folds into amphipathic α-helices in the presence of high concentrations of low-molecular-mass crowders or upon dehydration. Under these conditions, the (partially) folded protein binds peripherally to membranes. In our study, we have used coarse-grained molecular dynamics simulations to elucidate the details of COR15A-membrane binding and its effects on membrane structure and dynamics. Simulation results indicate that at least partial folding of COR15A and the presence of highly unsaturated galactolipids in the membranes are necessary for efficient membrane binding. The bound protein is stabilized on the membrane by interactions of charged and polar amino acids with galactolipid headgroups and by interactions of hydrophobic amino acids with the upper part of the fatty acyl chains. Experimentally, the presence of liposomes made from a mixture of lipids mimicking chloroplast membranes induces additional folding in COR15A under conditions of partial dehydration, in agreement with the simulation results.

Multiscale Simulations of Clavulanate Inhibition Identify the Reactive Complex in Class A ß-Lactamases and Predict the Efficiency of Inhibition.

Clavulanate is used as an effective drug in combination with ß-lactam antibiotics to treat infections of some antibiotic resistant bacteria. Here, we perform combined quantum mechanics/molecular mechanics simulations of several covalent complexes of clavulanate with class A ß-lactamases KPC-2 and TEM-1. Simulations of the deacylation reactions identify the decarboxylated trans-enamine complex as being responsible for inhibition. Further, the obtained free energy barriers discriminate clinically relevant inhibition (TEM-1) from less effective inhibition (KPC-2).

Design, facile synthesis, and evaluation of novel spiro- and pyrazolo[1,5-c]quinazolines as cholinesterase inhibitors: Molecular docking and MM/GBSA studies.

Given the wide spectrum of biological uses of pyrazolo[1,5-c]quinazoline and spiro-quinazoline derivatives as anticancer, anti-inflammatory analgesic agents, and their therapeutic applications in neurodegenerative disorders, it is compulsory to find easy, efficient, and simple methods to obtain and chemically diversify these families of compounds, thereby improving their biological applications. In this paper, we report the design and eco-friendly two-step synthesis of novel, fused spiro-pyrazolo[1,5-c]quinazoline derivatives as cholinesterase inhibitors. In addition, we studied their protein-ligand interactions via molecular docking and MM/GBSA calculations for a further rational design of more potent inhibitors. In first step, 2-(1H-pyrazol-5-yl)anilines were obtained through microwave (MW) assisted solvent-free/catalyst-free conditions and the second step involved the synthesis of the spiro-pyrazolo[1,5-c]quinazolines by a cyclocondensation reaction between 2-(1H-pyrazol-5-yl)anilines and cyclic ketones, or acetophenones, using stirring at room temperature. The compounds were obtained in high purity, good yields (50-97%), and at varying reaction times. The spiro-compounds were evaluated as acetylcholinesterase and butyrylcholinesterase inhibitors (AChEIs/BuChEIs) respectively, and the most potent compound exhibited a moderate AChE inhibitory activity (5f: IC50 = 84 µM). Molecular docking studies indicated that the binding mode of the compound 5f share common characteristics with the galantamine/donepezil-AChE complexes. Moreover, free binding energy (ΔG) calculations showed a good agreement with the experimental biological activity values. Our theoretical results indicated that halogen bond interactions could be involved with differential potency of these compounds and provide a new starting point to design novel pyrazolo[1,5-c]quinazolines as new anti-Alzheimer agents.