Search

Thanks to the considerable research which has been undertaken in the last few years to improve our understanding of the biology and mechanism of action of SARS‐CoV‐2, we know how the virus uses its surface spike protein to infect host cells. The transmembrane prosthesis, serine 2 (TMPRSS2) protein, located on the surface of human cells, recognizes the cleavage site in the spike protein, leading to the release of the fusion peptide and entry of the virus into the host cells. Because of its role, TMPRSS2 has been proposed as a drug target to prevent infection by the virus. In this study, we aim to increase our understanding of TMPRSS2 using long scale microsecond atomistic molecular dynamics simulations, focusing on the conformational changes over time. The comparison between simulations conducted on the protein in the native (apo) and inhibited form (holo), has shown that in the holo form the inhibitor stabilizes the catalytic site and induces rearrangements in the extracellular domain of the protein. In turn, it leads to the formation of a new cavity in the vicinity of the ligand binding pocket that is stable in the microsecond time scale. Given the low specificity of known protease inhibitors, these findings suggest a new potential drug target site that can be used to improve TMPRSS2 specific recognition by newly designed inhibitors. [ABSTRACT FROM AUTHOR]

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain–domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains. [ABSTRACT FROM AUTHOR]

Available estimates indicate that 30–40 % of all proteins need at least one metal ion to perform their biological function. Therefore, they are called metalloproteins. The correct biosynthesis of metalloproteins requires living organisms to be able to cope with issues such as the limited bioavailability or the potential cytotoxicity of several essential metals. Thus, organisms have developed complex machineries that guarantee the proper intracellular concentration and distribution among compartments of each metal, i.e. metal homeostasis. To understand how the different proteins responsible for metal homeostasis carry out their function, it is necessary to investigate their three‐dimensional (3D) structure and mobility at the atomic level. Nuclear magnetic resonance spectroscopy is one of the main experimental techniques providing this information. Computer simulations of molecular dynamics (MD) complement experimental information by showing how the 3D structure fluctuates over time and as a function of environmental conditions, with the possibility of exploring a wider range of timescales and conditions than usually amenable to experiment. Here we review numerous applications of MD for the investigation of the structure and dynamics of metalloproteins, and we also mention some technical aspects related to the parametrization of metals in commonly used force fields. Molecular dynamics simulations complement experimental information providing an atomistic view on biological functions. Here we discuss some applications to the study of metalloproteins, from the determination of their 3D structures, to their role in metal storage, transport and metal‐dependent transcriptional regulation. Some technical aspects related to the parametrization of metal ions are also mentioned. [ABSTRACT FROM AUTHOR]

Force-field parameters are developed for a multisite model of Ni(II) ions to be used in molecular dynamics simulations combined to enhanced sampling methods. The performances of two charge-partitioning schemes are validated by taking into account structural, thermodynamic, and kinetic observables. One of the two models, featuring partial charges on the dummy atoms only, matches both Ni(II) free energy of solvation and water exchange rates. Such model is particularly suited to study complexation events at a fully dynamic description. © 2017 Wiley Periodicals, Inc. [ABSTRACT FROM AUTHOR]

In oxygenic photosynthesis, plastocyanin shuttles electrons between the membrane-bound complexes cytochrome b6f and photosystem I. The homologous complex between cytochrome f and plastocyanin, both from spinach, is the object of this study. The solution structure of the reduced spinach plastocyanin was determined using high field NMR spectroscopy, whereas the model structure of oxidized cytochrome f was obtained by homology modeling calculations and molecular dynamics. The model structure of the intermolecular complex was calculated using the program AUTODOCK, taking into account biological information obtained from mutagenesis experiments. The best electron transfer pathway from the heme group of cytochrome f to the copper ion of plastocyanin was calculated using the program HARLEM, obtaining a coupling decay value of 1.8 × 10-4. Possible mechanisms of interaction and electron transfer between plastocyanin and cytochrome f were discussed considering the possible formation of a supercomplex that associates one cytochrome b6f, one photosystem I, and one plastocyanin. [ABSTRACT FROM AUTHOR]

HP1043 is an essential orphan response regulator of Helicobacter pylori orchestrating multiple crucial cellular processes. Classified as a member of the OmpR/PhoB family of two-component systems, HP1043 exhibits a highly degenerate receiver domain and evolved to function independently of phosphorylation. Here, we investigated the HP1043 binding mode to a target sequence in the hp1227 promoter (Php1227). Scanning mutagenesis of HP1043 DNA-binding domain and consensus sequence led to the identification of residues relevant for the interaction of the protein with a target DNA. These determinants were used as restraints to guide a data-driven protein-DNA docking. Results suggested that, differently from most other response regulators of the same family, HP1043 binds in a head-to-head conformation to the Php1227 target promoter. HP1043 interacts with DNA largely through charged residues and contacts with both major and minor grooves of the DNA are required for a stable binding. Computational alanine scanning on molecular dynamics trajectory was performed to corroborate our findings. Additionally, in vitro transcription assays confirmed that HP1043 positively stimulates the activity of RNA polymerase. [ABSTRACT FROM AUTHOR]

The ubiquitin pathway required for most proteins' targeted degradation involves three classes of enzymes: E1-activating enzyme, E2-conjugating enzyme, and E3-ligases. The human Ark2C is the single known E3 ligase that adopts an alternative, Ub-dependent mechanism for the activation of Ub transfer in the pathway. Its RING domain binds both E2-Ub and free Ub with high affinity, resulting in a catalytic active UbR-RING-E2-UbD complex formation. We examined potential changes in the conformational plasticity of the Ark2C RING domain and its ligands in their complexed form within the ubiquitin pathway through molecular dynamics (MD). Three molecular mechanics force fields compared to previous NMR relaxation studies of RING domain of Arkadia were used for effective and accurate assessment of MDs. Our results suggest the Ark2C Ub-RING docking site has a substantial impact on maintaining the conformational rigidity of E2-E3 assembly, necessary for the E3's catalytic activity. In the UbR-RING-E2-UbD catalytic complex, the UbR molecule was found to have greater mobility than the other Ub, bound to E2. Furthermore, network-based bioinformatics helped us identify E3 RING ligase candidates which potentially exhibit similar structural modules as Ark2C, along with predicted substrates targeted by the Ub-binding RING Ark2C. Our findings could trigger a further exploration of related unrevealed functions of various other E3 RING ligases. [ABSTRACT FROM AUTHOR]

Membrane receptors constitute major targets for pharmaceutical intervention. Drug design efforts rely on the identification of ligand binding poses. However, the limited experimental structural information available may make this extremely challenging, especially when only low-resolution homology models are accessible. In these cases, the predictions may be improved by molecular dynamics simulation approaches. Here we review recent developments of multiscale, hybrid molecular mechanics/coarse-grained (MM/CG) methods applied to membrane proteins. In particular, we focus on our in-house MM/CG approach. It is especially tailored for G-protein coupled receptors, the largest membrane receptor family in humans. We show that our MM/CG approach is able to capture the atomistic details of the receptor/ligand binding interactions, while keeping the computational cost low by representing the protein frame and the membrane environment in a highly simplified manner. We close this review by discussing ongoing improvements and challenges of the current implementation of our MM/CG code. [ABSTRACT FROM AUTHOR]

Variants found in the respiratory complex I (CI) subunit genes encoded by mitochondrial DNA can cause severe genetic diseases. However, it is difficult to establish a priori whether a single or a combination of CI variants may impact oxidative phosphorylation. Here we propose a computational approach based on coarse-grained molecular dynamics simulations aimed at investigating new CI variants. One of the primary CI variants associated with the Leber hereditary optic neuropathy (m.14484T>C/ MT-ND6) was used as a test case and was investigated alone or in combination with two additional rare CI variants whose role remains uncertain. We found that the primary variant positioned in the E-channel region, which is fundamental for CI function, stiffens the enzyme dynamics. Moreover, a new mechanism for the transition between π- and α-conformation in the helix carrying the primary variant is proposed. This may have implications for the E-channel opening/closing mechanism. Finally, our findings show that one of the rare variants, located next to the primary one, further worsens the stiffening, while the other rare variant does not affect CI function. This approach may be extended to other variants candidate to exert a pathogenic impact on CI dynamics, or to investigate the interaction of multiple variants. • CI variants found in mtDNA encoded subunit genes can cause severe genetic diseases. • Here we propose a CG MD approach aimed at the study of new CI variants. • We successfully applied our approach on three LHON variants. • The primary m.14484 T > C/ MT-ND6 variant stiffens CI dynamics in the E-channel region. • One of the other variants further worsens the stiffening of CI dynamics. [ABSTRACT FROM AUTHOR]

The survival of several pathogenic bacteria, such as Helicobacter pylori (Hp), relies on the activity of the nickel-dependent enzyme urease. Nickel insertion into urease is mediated by a multimeric chaperone complex (Hp UreDFG) that is responsible for the transport of Ni(II) from a conserved metal binding motif located in the UreG dimer (CPH motif) to the catalytic site of the enzyme. The X-ray structure of Hp UreDFG revealed the presence of water-filled tunnels that were proposed as a route for Ni(II) translocation. Here, we probe the transport of Ni(II) through the internal tunnels of Hp UreDFG, from the CPH motif to the external surface of the complex, using microsecond-long enhanced molecular dynamics simulations. The results suggest a "bucket-brigade" mechanism whereby Ni(II) can be transported through a series of stations found along these internal pathways. [Display omitted] • Helicobacter pylori relies on the Ni(II)-dependent enzyme urease. • Ni(II) insertion into urease requires a multimeric chaperone complex (Hp UreDFG). • Hp UreDFG structure shows internal tunnels that can be involved in Ni(II) transport. • Ni(II) translocation was probed using μs-long enhanced MD simulations. • The results suggest that Ni(II) can be transported through a series of stations. [ABSTRACT FROM AUTHOR]