Rui P. P. Neves

Protein-Ligand Docking in the New Millennium – A Retrospective of 10 Years in the Field

Protein-ligand docking is currently an important tool in drug discovery efforts and an active area of research that has been the subject of important developments over the last decade. These are well portrayed in the rising number of available protein-ligand docking software programs, increasing level of sophistication of its most recent applications, and growing number of users. While starting by summarizing the key concepts in protein-ligand docking, this article presents an analysis of the evolution of this important field of research over the past decade. Particular attention is given to the massive range of alternatives, in terms of protein-ligand docking software programs currently available. The emerging trends in this field are the subject of special attention, while old established docking alternatives are critically revisited. Current challenges in the field of protein-ligand docking such as the treatment of protein flexibility, the presence of structural water molecules and its effect in docking, and the entropy of binding are dissected and discussed, trying to anticipate the next years in the field.

Application of quantum mechanics/molecular mechanics methods in the study of enzymatic reaction mechanisms

Quantum mechanics/molecular mechanics (QM/MM) methods offer a very appealing option for the computational study of enzymatic reaction mechanisms, by separating the problem into two parts that can be treated with different computational methods. Hence, in a QM/MM formalism, the part of the system in which catalysis actually occurs and that involves the active site, substrates and directly participating amino acid residues is treated at an adequate quantum mechanical level to describe the chemistry taking place. For the remaining of the enzyme, which does not participate directly in the reaction, but that typically involves a much larger number of atoms, molecular mechanics is employed, traditionally through the application of a biomolecular force field. When applied with care, QM/MM methods can be used with great advantage in comparing, at a structural and energetic level, different mechanistic proposals, discarding mechanistic alternatives and proposing new mechanistic pathways that are consistent with the available experimental data. With time, diverse flavors within the QM/MM methods have emerged, differing in a variety of technical and conceptual aspects. Hence present alternatives differ between additive and subtractive QM/MM schemes, the type of boundary schemes, and the way in which the electrostatic interactions between the two regions are accounted for. Also, single‐conformation QM/MM, multi‐PES approaches, and QM/MM Molecular Dynamics coexist today, each type with its own advantages and limitations. This review focuses on the application of QM/MM methods in the study of enzymatic reaction mechanisms, briefly presenting also the most important technical aspects involved in these calculations. Particular attention is dedicated to the application of the single‐conformation QM/MM, multi‐PES QM/MM studies, and QM/MM‐FEP methods and to the advantages and disadvantages of the different types of QM/MM. Recent breakthroughs are also introduced. A selection of hand‐picked examples is used to illustrate such features. WIREs Comput Mol Sci 2017, 7:e1281. doi: 10.1002/wcms.1281

Cytological and Expression Studies and Quantitative Analysis of the Temporal and Stage-Specific Effects of Follicle-Stimulating Hormone and Testosterone During Cocultures of the Normal Human Seminiferous Epithelium

Abstract In vitro culturing of normal human seminiferous epithelium remains largely unexplored. To study normal human spermatogenesis in vitro, we used a micromethod for the purification and culture of Sertoli cells, spermatogonia A, spermatocytes, and early round spermatids. Cytological quantitative data for Sertoli and premeiotic germ cell cocultures isolated from normal testicular biopsies demonstrated that cells were able to proliferate (4%), complete meiosis (6.7%), and differentiate into late round (54%), elongating (49%), and elongated (17%) spermatids at similar in vivo time delays (up to 16 days) in response to FSH + testosterone stimulation. Cells maintained normal meiotic segregation, chromosome complements, and specific gene expression profiles. Follicle-stimulating hormone + testosterone stimulated spermatogonia proliferation and Sertoli cell survival. Follicle-stimulating hormone and especially FSH + testosterone increased diploid germ cell survival during the first week, whereas only FSH + testosterone was able to inhibit cell death during the second week of culture. Follicle-stimulating hormone and especially FSH + testosterone also stimulated meiosis resumption, although this was restricted to late pachytene and secondary spermatocytes. In contrast, spermiogenesis was only stimulated by FSH + testosterone. Expression studies showed that apoptosis was induced in the nucleus of diploid cells, and in nuclear and cytoplasmic compartments of spermatids, mainly triggered by the Fas pathway. Although junctional complexes between Sertoli and premeiotic germ cells were partially reacquired, the same did not apply to spermatids, suggesting that FSH potentiated by testosterone was unable to render Sertoli cells competent to bind round spermatids.

Parameters for Molecular Dynamics Simulations of Manganese-Containing Metalloproteins.

A set of geometrical parameters has been determined for single manganese metalloproteins for the AMBER force field, and ultimately to other force fields with a similar philosophy. Twelve (12) models from 9 different single-cluster manganese proteins were optimized and parametrized, using a bonded model approach. Mn-ligand bonds, Mn-ligand angles, and Restrained Electrostatic Potential charges for all the 74 residues in the first metal coordination sphere of each Mn metalloprotein were parametrized. The determined parameters were validated with molecular dynamics simulations and several statistics strategies were used to analyze the results. In addition, to validate the parametrized models, frequency and normal mode calculations were performed and comparisons were obtained for the overall structures both with quantum mechanics and molecular mechanics calculations. Linear and polynomial fittings were performed to estimate Mn-ligand bond force constants for generic manganese centers. Furthermore, averages are proposed for the main Mn-ligand angle interactions of typical manganese coordination centers: axial, square and triangular equatorial planes, and tetrahedral positions, for the different combinations of donor atoms from waters and hard ligands.

Benchmarking of Density Functionals for the Accurate Description of Thiol-Disulfide Exchange.

A set of 92 density functionals was employed to accurately characterize thiol-disulfide exchange. The properties we have benchmarked throughout the study include the geometry of a 15 atoms model system, the potential energy surface, the activation barrier, and the energy of reaction for thiol-disulfide exchange. Reference energies were determined at the CCSD(T)/CBS//MP2/aug-cc-pVDZ level of theory, and reference geometries were calculated at the MP2/aug-cc-pVTZ level. M11-L, M06-2X, M06-HF, N12-SX, PBE1PBE, PBEh1PBE, and OHSE2PBE described better the geometry of the model system, with average deviations of 0.06 Å in bond lengths (0.06 Å in bond-breaking lengths) and 1.9° in bond angles. On the other hand, the potential energy surface and its gradient were more accurately described by the hybrid density functional BHandH, closely followed by mPW1N, mPW1K, and mPWB1K. The barrier height and energy of reaction were better reproduced by the BMK and M06-2X functionals (deviations of 0.17 and 0.07 kcal·mol(-1), respectively) for a set of 10 Poples basis sets. MN12-SX and M11-L showed very good results for the widely used 6-311++G(2d,2p) basis set, with deviations of 0.02 and 0.05 kcal·mol(-1), respectively. We studied the effect of the split-valence, diffuse, and polarized functions in the activation barrier of thiol-disulfide exchange, for a set of 10 Poples basis sets. While increasing the splitting and polarization may increase the activation barrier in approximately 1 kcal·mol(-1), diffuse functions generally contribute to decreasing it no more than 0.10 kcal·mol(-1). In general, 13 functionals provided energies within 1 kcal·mol(-1) of the reference value. The BB1K density functional is one of the best density functionals to characterize thiol-disulfide exchange reactions; however, several density functionals with modified Perdew-Wang exchange and about 40% Hartree-Fock exchange, such as mPW1K, mPW1N, and mPWB1K, show a good performance, too.

Improving the Biodesulfurization of Crude Oil and Derivatives: A QM/MM Investigation of the Catalytic Mechanism of NADH-FMN Oxidoreductase (DszD)

The development of biocatalytic desulfurization strategies of petroleum and its derivatives could result in more economic alternatives than the widely used chemical desulfurization. The organism Rhodococcus erythropolis IGTS8 has been shown to metabolize organic sulfur compounds through a mechanism known as 4S pathway, which involves four enzymes (DszA, DszB, DszC, and DszD) and has been explored in biodesulfurization. Here we have applied QM/MM methods to study the catalytic mechanism of the enzyme DszD, a NADH-FMN oxidoreductase that occupies a central place on the 4S pathway by catalyzing the formation of the FMNH2 that is used by the two monooxynases in the cycle: DszA and DszC. In addition, to clarify the catalytic mechanism of this enzyme, this study analyzed in detail the role played by the active site Thr residue and of Asn and Ala enzyme mutants. The results help to explain previous experimental evidence and suggest new strategies for improving biodesulfurization through an increase in the activity of DszD.

Mechanistic insights on the reduction of glutathione disulfide by protein disulfide isomerase

Significance Protein disulfide isomerase (PDI) is a ubiquitous enzyme involved in disulfide bond formation during protein folding. It has been related to neurological diseases (Parkinson or Alzheimer’s) because of unfolded protein response phenomena. It also participates in the regulation of the glutathione redox buffer [glutathione/glutathione disulfide (GSH/GSSG)] in the endoplasmic reticulum, also important in protein folding. PDI catalyzes the nucleophilic attack of thiolates to disulfide bonds (thiol-disulfide exchange), enhancing the formation of correct disulfide links that drive protein folding and ensure protein function. This reaction is ubiquitous to disulfide-oxidoreductases across several organisms, and it shows a distinctive chemistry. We present the enzymatic mechanism of PDI in the GSH/GSSG buffer and discuss the chemistry of thiol-disulfide exchange. We explore the enzymatic mechanism of the reduction of glutathione disulfide (GSSG) by the reduced a domain of human protein disulfide isomerase (hPDI) with atomistic resolution. We use classical molecular dynamics and hybrid quantum mechanics/molecular mechanics calculations at the mPW1N/6–311+G(2d,2p):FF99SB//mPW1N/6–31G(d):FF99SB level. The reaction proceeds in two stages: (i) a thiol-disulfide exchange through nucleophilic attack of the Cys53-thiolate to the GSSG-disulfide followed by the deprotonation of Cys56-thiol by Glu47-carboxylate and (ii) a second thiol-disulfide exchange between the Cys56-thiolate and the mixed disulfide intermediate formed in the first step. The Gibbs activation energy for the first stage was 18.7 kcal·mol−1, and for the second stage, it was 7.2 kcal·mol−1, in excellent agreement with the experimental barrier (17.6 kcal·mol−1). Our results also suggest that the catalysis by protein disulfide isomerase (PDI) and thiol-disulfide exchange is mostly enthalpy-driven (entropy changes below 2 kcal·mol−1 at all stages of the reaction). Hydrogen bonds formed between the backbone of His55 and Cys56 and the Cys56-thiol result in an increase in the Gibbs energy barrier of the first thiol-disulfide exchange. The solvent plays a key role in stabilizing the leaving glutathione thiolate formed. This role is not exclusively electrostatic, because an explicit inclusion of several water molecules at the density-functional theory level is a requisite to form the mixed disulfide intermediate. In the intramolecular oxidation of PDI, a transition state is only observed if hydrogen bond donors are nearby the mixed disulfide intermediate, which emphasizes that the thermochemistry of thiol-disulfide exchange in PDI is influenced by the presence of hydrogen bond donors.

Modelling of Human Fatty Acid Synthase and in Silico Docking of Acyl Carrier Protein Domain and Its Partner Catalytic Domains

Human fatty acid synthase (hFAS) is a megasynthase whose main function is de novo biosynthesis of saturated fatty acids. Interest has been drawn to this enzyme beyond its physiological role due to the association between high levels of hFAS and clinical conditions such as obesity, diabetes, and cancer. Thus, it has become an undeniably attractive pharmacological target. Until now, no crystal structure of the complete hFAS is available, hindering attempts to fully understand this protein. Using homology modeling, we built a model of the entire megasynthase, encompassing all of its domains, including the acyl carrier protein (ACP) and thioesterase (TE) mobile domains absent in the crystal structure of mammalian fatty acid synthase (FAS). On a second stage, we used data-driven protein-protein docking between the substrate shuttling domain ACP and every catalytic domain in the protein. We also propose sets of amino acids at the interface of each domain that we believe are important to favor the interaction between ACP and each domain of hFAS. After inspection, we validated each complex between ACP and MAT/KS/KR/DH/ER domains through classical molecular dynamics simulations and RMSd analysis. Additionally, we mapped the interactions between the residues at the active site of each catalytic domain and its intermediaries. In every docking, we ensured that the distance between catalytic residues and the intermediaries was maintained. Until now, there was not a complete 3D model of this megasynthase. This study is the first to present a homology model for the whole hFAS, including its two mobile domains and possible poses of ACP throughout the cycle of fatty acid biosynthesis, thus mapping obligatory checkpoints in its trajectory. Hence, we believe that these structural insights will allow for future studies of the catalytic mechanism of the overall hFAS.