RESUMO
The emergence of new viruses can lead to the outbreak of pandemics as occurred at the end of 2019 with the coronavirus disease (or COVID-19). The fastest way to effectively control viral infections is to develop broad-spectrum antivirals that can fight at least an entire class of viruses. Profurin, the furin precursor propeptide, is responsible for the autoactivation step which is crucial for the maturation of several viral substrates. This role makes the study of furin and profurin interactions interesting for the development of new potential broad-spectrum antivirals for the treatment against several human viral diseases. Since there is no 3D model of profurin published in the literature or deposited in a database, this work reports the development, validation and analysis of a profurin 3D model using comparative modeling and molecular dynamics. The model is available in ModelArchive at https://www.modelarchive.org/doi/10.5452/ma-ct8l7. The usage of this model will make possible further studies of molecular docking and MD simulations of the profurin-furin system, in the design of new potential broad-spectrum antivirals for the treatment against several human viral diseases.Communicated by Ramaswamy H. Sarma.
RESUMO
In this research, a combined study on structures and vibrational spectra of antiviral rimantadine have been performed using hybrid B3LYP/6-311++G∗∗ calculations and the scaled quantum force field (SQMFF) procedure. Harmonic force fields and scaled force constants of Free Base (FB), Cationic (CA) and Hydrochloride (HCl) species derived from the antiviral rimantadine have been calculated in gas phase and in aqueous solution using normal internal coordinates and scaling factors. Good correlations were acquired comparing the theoretical IR, Raman, 1H- 13C-NMR and UV spectra of three species with the analogous experimental ones, suggesting probably, the presence of all them in both phases. The main force constants of three species have evidenced lower values than the corresponding to antiviral amantadine. The ionic character of N1-H33â¯Cl36 bond of HCl species in aqueous solution evidence positive Mulliken charge on N1 atom indicating that this species is as CA one. Rimantadine presents higher solvation energies in water than other antiviral species, such as chloroquin, niclosamide, cidofovir and brincidofovir. The FB and HCl species of rimantadine are slightly less reactive than the corresponding to amantadine while the opposite is observed for the CA species. The predicted ECD spectra for the FB and CA species show positive Cotton effect different from the negative observed for the HCl one. These different behaviours of three species of rimantadine could probably explain the differences observed in the intensities of bands predicted in the electronic spectra of these species.
RESUMO
We analyze the influence of different groups on the intermolecular energy of aromatic homodimers and on the interaction between a single aromatic molecule and a graphene surface. The analysis is performed for benzene, phenol, catechol, and dopamine. For calculating the energies, we employ density functional theory within the local density approximation (LDA-DFT). Our results show that the lowest intermolecular energies between the aromatic molecules are related to the T-shaped configurations. This lower energy results from the quadrupole interaction. In the case of the interaction between the graphene sheet and the aromatic molecules, the lowest energy configuration is the face to face. The adsorption energy of a molecule on a graphene surface involves π - π interactions that explain the face to face arrangement. These results provide insight into the manner by which substituents can be utilized in crystal engineering, supramolecular chemistry, bioinspired materials, formation of various molecular clusters, parameterization of force fields suitable for classical simulations, and design of novel sensing, drug delivery, and filters based on graphene.
RESUMO
X-ray diffraction crystallography is the primary technique to determine the three-dimensional structures of biomolecules. Although a robust method, X-ray crystallography is not able to access the dynamical behavior of macromolecules. To do so, we have to carry out molecular dynamics simulations taking as an initial system the three-dimensional structure obtained from experimental techniques or generated using homology modeling. In this chapter, we describe in detail a tutorial to carry out molecular dynamics simulations using the program NAMD2. We chose as a molecular system to simulate the structure of human cyclin-dependent kinase 2.
Assuntos
Simulação de Dinâmica Molecular , Software , Trifosfato de Adenosina/química , Algoritmos , Cristalografia por Raios X , Quinase 2 Dependente de Ciclina/química , Humanos , Conformação Proteica , Eletricidade Estática , Interface Usuário-ComputadorRESUMO
Implicit-solvent models are widely used to study the electrostatics in dissolved biomolecules, which are parameterized using force fields. Standard force fields treat the charge distribution with point charges; however, other force fields have emerged which offer a more realistic description by considering polarizability. In this work, we present the implementation of the polarizable and multipolar force field atomic multipole optimized energetics for biomolecular applications (AMOEBA), in the boundary integral Poisson-Boltzmann solver PyGBe. Previous work from other researchers coupled AMOEBA with the finite-difference solver APBS, and found difficulties to effectively transfer the multipolar charge description to the mesh. A boundary integral formulation treats the charge distribution analytically, overlooking such limitations. This becomes particularly important in simulations that need high accuracy, for example, when the quantity of interest is the difference between solvation energies obtained from separate calculations, like happens for binding energy. We present verification and validation results of our software, compare it with the implementation on APBS, and assess the efficiency of AMOEBA and classical point-charge force fields in a Poisson-Boltzmann solver. We found that a boundary integral approach performs similarly to a volumetric method on CPU. Also, we present a GPU implementation of our solver. Moreover, with a boundary element method, the mesh density to correctly resolve the electrostatic potential is the same for standard point-charge and multipolar force fields. Finally, we saw that for binding energy calculations, a boundary integral approach presents more consistent results than a finite difference approximation for multipolar force fields. © 2019 Wiley Periodicals, Inc.
RESUMO
Essential to understanding life, the biomolecular phenomena have been an important subject in science, therefore a necessary path to be covered to make progress in human knowledge. To fully comprehend these processes, the non-covalent interactions are the key. In this review, we discuss how specific protein-ligand interactions can be efficiently described by low computational cost methods, such as Molecular Mechanics (MM). We have taken as example the case of the halogen bonds (XB). Albeit generally weaker than the hydrogen bonds (HB), the XBs play a key role to drug design, enhancing the affinity and selectivity toward the biological target. Along with the attraction between two electronegative atoms in XBs explained by the σ-hole model, important orbital interactions, as well as relief of Pauli repulsion take place. Nonetheless, such electronic effects can be only well-described by accurate quantum chemical methods that have strong limitations dealing with supramolecular systems due to their high computational cost. To go beyond the poor description of XBs by MM methods, reparametrizing the force-fields equations can be a way to keep the balance between accuracy and computational cost. Thus, we have shown the steps to be considered when parametrizing force-fields to achieve reliable results of complex non-covalent interactions at MM level for In Silico drug design methods.
RESUMO
Prions are proteins that cause a group of invariably fatal neurodegenerative diseases, one of the most known being bovine spongiform encephalopathy. The three-dimensional structure of PrPSc , the altered isoform of the prion protein, has not been fully elucidated yet, and studies on prion conversion mechanisms must rely on hypothetical ß-rich structures. Experimental and computational studies indicate that the use of low pH is capable to produce a gain of ß-structure content in the otherwise unstructured N-terminal region. These in silico studies have used different PrP fragments from distinct organisms, and with different lengths and simulation protocols, making it difficult to identify the influence of the force fields on the formation of such structures. Here, we performed a systematic study of the influence of six well-established force fields (GROMOS96 53a6, GROMOS96 43a1, AMBER99SB, AMBER99SB-ILDN, CHARMM27, and OPLS-AA/L) on the process of structural conversion of the Syrian hamster cellular prion protein simulated at acidic and neutral pH. From our analysis, we observe a strong dependence of the results with the different force fields employed. Additionally, only GROMOS96 53A6 and AMBER99SB force fields are capable to capture a high ß-sheet formation at acidic pH and adequately reproduce the neutral pH. In both cases, the ß-sheet elongation seems to be guided by the movement of the N-terminal tail toward the N-terminal of α-helix HB under acidic condition. These results comprise the most wide-ranging study to date correlating force fields to structural changes in the cellular prion protein. © 2018 Wiley Periodicals, Inc.