Unveiling functional motions based on point mutations in biased signaling systems: A normal mode study on nerve growth factor bound to TrkA

Unveiling mutation effects on the structural dynamics of the main protease from SARS-CoV-2 with hybrid simulation methods

The main protease of SARS-CoV-2 (called M^pro or 3CL^pro) is essential for processing polyproteins encoded by viral RNA. Several M^pro mutations were found in SARS-CoV-2 variants, which are related to higher transmissibility, pathogenicity, and resistance to neutralization antibodies. Macromolecules adopt several favored conformations in solution depending on their structure and shape, determining their dynamics and function. In this study, we used a hybrid simulation method to generate intermediate structures along the six lowest frequency normal modes and sample the conformational space and characterize the structural dynamics and global motions of WT SARS-CoV-2 M^pro and 48 mutations, including mutations found in P.1, B.1.1.7, B.1.351, B.1.525 and B.1.429+B.1.427 variants. We tried to contribute to the elucidation of the effects of mutation in the structural dynamics of SARS-CoV-2 M^pro. A machine learning analysis was performed following the investigation regarding the influence of the K90R, P99L, P108S, and N151D mutations on the dimeric interface assembling of the SARS-CoV-2 M^pro. The parameters allowed the selection of potential structurally stable dimers, which demonstrated that some single surface aa substitutions not located at the dimeric interface (K90R, P99L, P108S, and N151D) are able to induce significant quaternary changes. Furthermore, our results demonstrated, by a Quantum Mechanics method, the influence of SARS-CoV-2 M^pro mutations on the catalytic mechanism, confirming that only one of the chains of the WT and mutant SARS-CoV-2 M^pros are prone to cleave substrates. Finally, it was also possible to identify the aa residue F140 as an important factor related to the increasing enzymatic reactivity of a significant number of SARS-CoV-2 M^pro conformations generated by the normal modes-based simulations.

Effect of phosphorylation on the structural dynamics, thermal stability of human dopamine transporter: A simulation study using normal modes, molecular dynamics and Markov State Model

The Human Dopamine Transporter (hDAT) plays an essential role in modulating the Influx/Efflux of dopamine, and it is involved in the mechanism of certain neurodegenerative diseases such as Parkinson's disease. Several studies have reported important states for Dopamine transport: outward-facing open state (OFo), the outward-facing closed state (OFc), the holo-occluded state closed (holo), and the inward-facing open state (IFo). Furthermore, experimental assays have shown that different phosphorylation conditions in hDAT can affect the rate of dopamine absorption. We present a protocol using hybrid simulation methods to study the conformational dynamics and stability of states of hDAT under different phosphorylation sites. With this protocol, we explored the conformational space of hDAT, identified the states, and evaluated the free energy differences and the transition probabilities between them in each of the phosphorylation cases. We also presented the conformational changes and correlated them with those described in the literature. There is a thesis/hypothesis that the phosphorylation condition corresponding to NP-333 system (where all sites Ser/Thr from residue 2 to 62 and 254 to 613 are phosphorylated, except residue 333) would decrease the rate of dopamine transport from the extracellular medium to the intracellular medium by hDAT as previously described in the literature by Lin et al., 2003. Our results corroborated this thesis/hypothesis and the data reported. It is probably due to the affectation/changes/alteration of the conformational dynamics of this system that makes the intermediate states more likely and makes it difficult to initial states associated with the uptake of dopamine in the extracellular medium, corroborating the experimental results. Furthermore, our results showed that just single phosphorylation/dephosphorylation could alter intrinsic protein motions affecting the sampling of one or more states necessary for dopamine transport. In this sense, the modification of phosphorylation influences protein movements and conformational preferences, affecting the stability of states and the transition between them and, therefore, the transport.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

HARP: a database of structural impacts of systematic missense mutations in drug targets of Mycobacterium leprae

Computational Saturation Mutagenesis is an in-silico approach that employs systematic mutagenesis of each amino acid residue in the protein to all other amino acid types, and predicts changes in thermodynamic stability and affinity to the other subunits/protein counterparts, ligands and nucleic acid molecules. The data thus generated are useful in understanding the functional consequences of mutations in antimicrobial resistance phenotypes. In this study, we applied computational saturation mutagenesis to three important drug-targets in Mycobacterium leprae (M. leprae) for the drugs dapsone, rifampin and ofloxacin namely Dihydropteroate Synthase (DHPS), RNA Polymerase (RNAP) and DNA Gyrase (GYR), respectively. M. leprae causes leprosy and is an obligate intracellular bacillus with limited protein structural information associating mutations with phenotypic resistance outcomes in leprosy. Experimentally solved structures of DHPS, RNAP and GYR of M. leprae are not available in the Protein Data Bank, therefore, we modelled the structures of these proteins using template-based comparative modelling and introduced systematic mutations in each model generating 80,902 mutations and mutant structures for all the three proteins. Impacts of mutations on stability and protein-subunit, protein-ligand and protein-nucleic acid affinities were computed using various in-house developed and other published protein stability and affinity prediction software. A consensus impact was estimated for each mutation using qualitative scoring metrics for physicochemical properties and by a categorical grouping of stability and affinity predictions. We developed a web database named HARP (a database of Hansen's Disease Antimicrobial Resistance Profiles), which is accessible at the URL - https://harp-leprosy.org and provides the details to each of these predictions.