Search

Understanding the plasmonic coupling between a set of metallic nanoparticles (NPs) in a 2D array, and how a substrate affects such coupling, is fundamental for the development of optimized optoelectronic structures. Here, a simple semi-analytical procedure based on discrete dipole approximation (DDA) is reported to simulate the far-field and near-field properties of arrays of NPs, considering the coupling between particles, and the effect of the presence of a semiconductor substrate based on the image dipole approach. The method is validated for Ag NP dimers and single Ag NPs on a gallium nitride (GaN) substrate, a semiconductor widely used in optical devices, by comparison with the results obtained by the finite element method (FEM), indicating a good agreement in the weak coupling regime. Next, the method is applied to square and random arrays of Ag NPs on a GaN substrate. The increase in the surface density of NPs on a GaN substrate mainly results in a redshift of the dipolar resonance frequency and an increase in the near-field enhancement. This model, based on a single dipole approach, grants very low computational times, representing an advantage to predict the optical properties of large NP arrays on a semiconductor substrate for different applications.

Cyclic GMP-AMP Synthase (cGAS) is activated upon DNA binding and catalyzes the synthesis of 2',3'-cGAMP from GTP and ATP. This cyclic dinucleotide is a messenger that triggers the autoimmune system of eukaryotic cells. In this study, we propose a Molecular Dynamics (MD) investigation of cGAS activation. We notably provide insights into the motion of the activation loop, both from a mechanical point of view and considering its role in the catalysis of cGAMP production. We finally shed light on the reaction resulting in cGAMP synthesis. Two possible catalytic routes (referred to as GTP-ATP and ATP-GTP) are proposed based on the active site occupancy, paving the way toward further exploration of the reaction mechanism.

Human aldehyde oxidase (hAOX1) is a molybdenum dependent enzyme that plays an important role in the metabolism of various compounds either endogenous or xenobiotics. Due to its promiscuity, hAOX1 plays a major role in the pharmacokinetics of many drugs and therefore has gathered a lot of attention from the scientific community and, particularly, from the pharmaceutical industry. In this work, homology modelling, molecular docking and molecular dynamics simulations were used to study the structure of the monomer and dimer of human AOX. The results with the monomer of hAOX1 allowed to shed some light on the role played by thioridazine and two malonate ions that are co-crystalized in the recent X-ray structure of hAOX1. The results show that these molecules endorse several conformational rearrangements in the binding pocket of the enzyme and these changes have an impact in the active site topology as well as in the stability of the substrate (phthalazine). The results show that the presence of both molecules open two gates located at the entrance of the binding pocket, from which results the flooding of the active site. They also endorse several modifications in the shape of the binding pocket (namely the position of Lys893) that, together with the presence of the solvent molecules, favour the release of the substrate to the solvent. Further insights were also obtained with the assembled homodimer of hAOX1. The allosteric inhibitor (THI) binds closely to the region where the dimerization of both monomers occur. These findings suggest that THI can interfere with protein dimerization.

The development of docking algorithms to predict near-native structures of protein:protein complexes from the structure of the isolated monomers is of paramount importance for molecular biology and drug discovery. In this study, we assessed the capacity of the interfacial area of protein:protein complexes and of Molecular Mechanics-Poisson Boltzmann Surface Area (MM-PBSA)-derived properties, to rank docking poses. We used a set of 48 protein:protein complexes, and a total of 67 docking experiments distributed among bound:bound, bound:unbound, and unbound:unbound test cases. The MM-PBSA binding free energy of protein monomers has been shown to be very convenient to predict high-quality structures with a high success rate. In fact, considering solely the top-ranked pose of more than 200 docking solutions of each of 39 protein:protein complexes, the success rate was 77% in the prediction of high-quality poses, or 90% if considering high- or medium-quality poses. If considering high- or medium-quality poses as the top-one prediction, a success rate of 87% was obtained for a scoring scheme based on computational alanine scanning mutagenesis data. Such ranking accuracy highlights the ability of these properties to predict near-native poses in protein:protein docking.

In this paper, we report a theoretical investigation of the catalytic mechanism of peptide amidases that involve a Ser-(cis)Ser-Lys catalytic triad. Previous suggestions propose that these enzymes should follow a distinct catalytic mechanism from the one that is present in the classic Ser-His-Asp catalytic triad. The theoretical and computational results obtained in this work indicate the opposite idea, showing that both mechanisms are very similar and only few differences are observed between both reactions. The results reveal that the different alignment of the Ser-(cis)Ser-Lys catalytic triad in relation to the classical Ser-His-Asp triad may provide a better stabilisation of the reaction intermediates, and therefore make these enzymes catalytically more efficient. The catalytic mechanism has been determined at the M06-2X/6-311++G**//B3LYP/6-31G* level of theory and requires five sequential steps instead of the two that are generally proposed: (i) nucleophilic attack of serine on the carbonyl group of the substrate, forming the first tetrahedral intermediate, (ii) formation of an acyl-enzyme complex, (ii) release of an ammonia product, (iv) nucleophilic attack of a water molecule forming the second tetrahedral intermediate, and (iv) the release of the product of the reaction, the carboxylic acid. The computational results suggest that the rate-limiting step is the first one that requires an activation free energy of 15.93 kcal mol -1 . This result agrees very well with the available experimental data that predict a reaction rate of 2200 s -1 , which corresponds to a free energy barrier of 14 kcal mol -1 .

Human fatty acid synthase (hFAS) is a multifunctional enzyme involved in a wide diversity of biological functions. For instance, it is a precursor of phospholipids and other complex processes such as the de novo synthesis of long chain fatty acid. Human FAS is also a component of biological membranes and it is implicated in the overexpression of several types of cancers. In this work, we describe the catalytic mechanism of β-ketoreductase (KR), which is a catalytic domain of the hFAS enzyme that catalyzes the reduction of β-ketoacyl to β-hydroxyacyl with the concomitant oxidation of the NADPH cofactor. The catalysis by KR is an intermediate step in the cycle of reactions that elongate the substrate's carbon chain until the final product is obtained. We study and propose the catalytic mechanism of the KR domain determined using the hybrid QM/MM methodology, at the ONIOM(B3LYP/6-311+G(2d,2p):AMBER) level of theory. The results indicate that the reaction mechanism occurs in two stages: (i) nucleophilic attack by a NADPH hydride to the β-carbon of the substrate, together with an asynchronous deprotonation of the Tyr2034 by the oxygen of the β-alkoxide to hold the final alcohol product; and (ii) an asynchronous deprotonation of the hydroxyl in the NADP + 's ribose by Tyr2034, and of the Lys1995 by the resulting alkoxide in the former ribose to restore the protonation state of Tyr2034. The reduction step occurs with a Gibbs energy barrier of 11.7 kcal mol -1 and a Gibbs reaction energy of -10.6 kcal mol -1 . These results have provided an understanding of the catalytic mechanism of the KR hFAS domain, a piece of the heavy hFAS biosynthetic machinery.

3-O-Sulfotransferase (3-OST) is one of the enzymes involved in heparan sulfate (HS) biosynthesis. HSs are polysaccharides with variable patterns of sulfation and acetylation that serve as entry receptors for herpes simplex virus type 1 (HSV-1). 3-OST is responsible for the transfer of a sulfate group from 3'-phosphoadenosine-5'-phosphosulfate (PAPS) to glucosamine units of HS. In this work, the catalytic mechanism of 3-OST was studied with atomic detail, using computational methods. We investigated the protonation state of key residues using the H++ web-based pKa prediction tool and molecular dynamics (MD) simulations and estimated the most relevant protonation state of the catalytic residues during catalysis. Catalytic histidine (His186) is predominantly protonated, while catalytic aspartate and glutamate (Asp189 and Glu184) are predominantly deprotonated. Subsequently, to study the catalytic mechanism, we applied a QM/MM method at the ONIOM(B3LYP/6-31G(d):ff94) level, starting from three geometries extracted from the 3, 6 and 8 ns point on the MD simulation. The results show that the reaction mechanism of 3-OST occurs by a single elementary step, consisting of an associative SN2 transfer of the sulfate group from PAPS to the HS glucosamine units, with the transfer of a proton from glucosamine to the catalytic Glu184. The activation free energies for this reaction were determined at the ONIOM(M06-2X-D3/6-311++G(2d,2p):ff94//B3LYP/6-31G(d):ff94) level of theory. Despite the free energy differences among the three conformations (10.2, 20.9 and 16.1 kcal mol(-1)), our results are consistent with the upper limit determined experimentally for the full cycle (20.4 kcal mol(-1)). The data obtained in this study will be useful for further studies on the inhibition of this enzyme, which is a useful target for drugs that block HSV-1 viral infections.

Protein-protein (P-P) 3D structures are fundamental to structural biology and drug discovery. However, most of them have never been determined. Many docking algorithms were developed for that purpose, but they have a very limited accuracy in generating native-like structures and identifying the most correct one, in particular when a single answer is asked for. With such a low success rate it is difficult to point out one docked structure as being native-like. Here we present a new, high accuracy, scoring method to identify the 3D structure of P-P complexes among a set of trial poses. It incorporates alanine scanning mutagenesis experimental data that need to be obtained a priori. The scoring scheme works by matching the computational and the experimental alanine scanning mutagenesis results. The size of the trial P-P interface area is also taken into account. We show that the method ranks the trial structures and identifies the native-like structures with unprecedented accuracy (∼94%), providing the correct P-P 3D structures that biochemists and molecular biologists need to pursue their studies. With such a success rate, the bottleneck of protein-protein docking moves from the scoring to searching algorithms.

We report on the viability of breaking selected bonds in biological systems using tailored electromagnetic radiation. We first demonstrate, by performing large-scale simulations, that pulsed electric fields cannot produce selective bond breaking. Then, we present a theoretical framework for describing selective energy concentration on particular bonds of biomolecules upon application of tailored electromagnetic radiation. The theory is based on the mapping of biomolecules to a set of coupled harmonic oscillators and on optimal control schemes to describe optimization of temporal shape, the phase and polarization of the external radiation. We have applied this theory to demonstrate the possibility of selective bond breaking in the active site of bacterial DNA topoisomerase. For this purpose, we have focused on a model that was built based on a case study. Results are given as a proof of concept.

Enzymes play a biologically essential role in performing and controlling an important share of the chemical processes occurring in life. However, despite their critical role in nature, attaining a clear understanding of the way an enzyme acts, i.e. its catalytic mechanism, is a cumbersome task that requires the cooperative efforts of a large number of different scientific techniques. Computational methods offer a particularly insightful way to study such mechanisms, always beautifully complementing the information arising from experimental techniques and working as an excellent alternative for assessing the viability of different mechanistic proposals. This review highlights two important computational strategies to study enzymatic catalysis - the cluster modeling approach and the hybrid quantum mechanical/molecular mechanical (QM/MM) method - complemented with a selection of hand-picked examples of our own work.

Hypertension is a chronic condition that affects nearly 25% of adults worldwide. As the Renin-Angiotensin-Aldosterone System is implicated in the control of blood pressure and body fluid homeostasis, its combined blockage is an attractive therapeutic strategy currently in use for the treatment of several cardiovascular conditions. We have performed QM/MM calculations to study the mouse renin catalytic mechanism in atomistic detail, using the N-terminal His6-Asn14 segment of angiotensinogen as substrate. The enzymatic reaction (hydrolysis of the peptidic bond between residues in the 10th and 11th positions) occurs through a general acid/base mechanism and, surprisingly, it is characterized by three mechanistic steps: it begins with the creation of a first very stable tetrahedral gem-diol intermediate, followed by protonation of the peptidic bond nitrogen, giving rise to a second intermediate. In a final step the peptidic bond is completely cleaved and both gem-diol hydroxyl protons are transferred to the catalytic dyad (Asp32 and Asp215). The final reaction products are two separate peptides with carboxylic acid and amine extremities. The activation energy for the formation of the gem-diol intermediate was calculated as 23.68 kcal mol(-1), whereas for the other steps the values were 15.51 kcal mol(-1) and 14.40 kcal mol(-1), respectively. The rate limiting states were the reactants and the first transition state. The associated barrier (23.68 kcal mol(-1)) is close to the experimental values for the angiotensinogen substrate (19.6 kcal mol(-1)). We have also tested the influence of the density functional on the activation and reaction energies. All eight density functionals tested (B3LYP, B3LYP-D3, X3LYP, M06, B1B95, BMK, mPWB1K and B2PLYP) gave very similar results.