RedMDStream: Parameterization and Simulation Toolbox for Coarse-Grained Molecular Dynamics Models

AbstractCoarse-grained (CG) models in molecular dynamics (MD) are powerful tools to simulate the dynamics of large biomolecular systems on micro- to millisecond timescales. However, the CG model, potential energy terms, and parameters are typically not transferable between different molecules and problems. So parameterizing CG force fields, which is both tedious and time-consuming, is often necessary. We present RedMDStream, a software for developing, testing, and simulating biomolecules with CG MD models. Development includes an automatic procedure for the optimization of potential energy parameters based on metaheuristic methods. As an example we describe the parameterization of a simple CG MD model of an RNA hairpin

Peptidomimetic inhibitors targeting the membrane-binding site of the neutrophil proteinase 3

Proteinase 3 (PR3), together with other serine proteases, such as neutrophil elastase (NE) and cathepsin G (CG), regulates inflammatory and immune responses. However, in comparison with NE and CG, there is increasing evidence that PR3 functions significantly differ. In particular, PR3 can bind to cell membranes and such membrane-bound PR3 (mbPR3) might be differently involved in the activation of cytokines, growth factors, cellular receptors, and in the regulation of cell apoptosis. For instance, PR3 membrane binding can block some “eat me” signals, notably, phosphatidylserine membrane lipid, and facilitate non-resolving inflammation. Based on the clear evidence that PR3 membrane binding affects the biological functions of PR3, we designed peptidomimetic inhibitors that can remove mbPR3 from the membrane surface in vitro without influencing PR3 catalytic activity. Such inhibitors, which specifically target PR3 binding to membranes, are still lacking. In particular, we found peptidomimetics that inhibit binding of PR3 to POPC:PS liposomes, which mimic the biological environment of PR3.publishedVersio

Dependency Map of Proteins in the Small Ribosomal Subunit

The assembly of the ribosome has recently become an interesting target for antibiotics in several bacteria. In this work, we extended an analytical procedure to determine native state fluctuations and contact breaking to investigate the protein stability dependence in the 30S small ribosomal subunit of Thermus thermophilus. We determined the causal influence of the presence and absence of proteins in the 30S complex on the binding free energies of other proteins. The predicted dependencies are in overall agreement with the experimentally determined assembly map for another organism, Escherichia coli. We found that the causal influences result from two distinct mechanisms: one is pure internal energy change, the other originates from the entropy change. We discuss the implications on how to target the ribosomal assembly most effectively by suggesting six proteins as targets for mutations or other hindering of their binding. Our results show that by blocking one out of this set of proteins, the association of other proteins is eventually reduced, thus reducing the translation efficiency even more. We could additionally determine the binding dependency of THX—a peptide not present in the ribosome of E. coli—and suggest its assembly path

Antibacterial Peptide Nucleic Acids—Facts and Perspectives

Antibiotic resistance is an escalating, worldwide problem. Due to excessive use of antibiotics, multidrug‐resistant bacteria have become a serious threat and a major global healthcare problem of the 21st century. This fact creates an urgent need for new and effective antimicrobials. The common strategies for antibiotic discovery are based on either modifying existing antibiotics or screening compound libraries, but these strategies have not been successful in recent decades. An alternative approach could be to use gene‐specific oligonucleotides, such as peptide nucleic acid (PNA) oligomers, that can specifically target any single pathogen. This approach broadens the range of potential targets to any gene with a known sequence in any bacterium, and could significantly reduce the time required to discover new antimicrobials or their redesign, if resistance arises. We review the potential of PNA as an antibacterial molecule. First, we describe the physicochemical properties of PNA and modifications of the PNA backbone and nucleobases. Second, we review the carriers used to transport PNA to bacterial cells. Furthermore, we discuss the PNA targets in antibacterial studies focusing on antisense PNA targeting bacterial mRNA and rRNA

Irish Area Section (Protein Interactions in Biology) Interactions of aminoglycoside antibiotics with rRNA

Abstract Aminoglycoside antibiotics are protein synthesis inhibitors applied to treat infections caused mainly by aerobic Gram-negative bacteria. Due to their adverse side effects they are last resort antibiotics typically used to combat pathogens resistant to other drugs. Aminoglycosides target ribosomes. We describe the interactions of aminoglycoside antibiotics containing a 2-deoxystreptamine (2-DOS) ring with 16S rRNA. We review the computational studies, with a focus on molecular dynamics (MD) simulations performed on RNA models mimicking the 2-DOS aminoglycoside binding site in the small ribosomal subunit. We also briefly discuss thermodynamics of interactions of these aminoglycosides with their 16S RNA target

Editorial: Combining Simulations, Theory, and Experiments into Multiscale Models of Biological Events

10.3389/fmolb.2021.797754Frontiers in Molecular Biosciences879775

Collective Dynamics of the Ribosomal Tunnel Revealed by Elastic Network Modeling

The collective dynamics of the nascent polypeptide exit tunnel are investigated with the computationally efficient elastic network model using normal mode analysis. The calculated normal modes are considered individually and in linear combinations with different coefficients mimicking the phase angles between modes, in order to follow the mechanistic motions of tunnel wall residues. The low frequency fluctuations indicate three distinct regions along the tunnel - the entrance, the neck and the exit – each having distinctly different domain motions. Generally the lining of the entrance region moves in the exit direction, with the exit region having significantly larger motions, but in a perpendicular direction, whereas the confined neck region generally has rotational motions. Especially the universally conserved extensions of ribosomal proteins L4 and L22 located at the narrowest and mechanistically strategic region of tunnel undergo generally anti- or non-correlated motions, which may have an important role in nascent polypeptide gating mechanism. These motions appear to be sufficiently robust so as to be unaffected by the presence of a peptide chain in the tunnel

Dynamics of the Acetylcholinesterase Tetramer

Acetylcholinesterase rapidly hydrolyzes the neurotransmitter acetylcholine in cholinergic synapses, including the neuromuscular junction. The tetramer is the most important functional form of the enzyme. Two low-resolution crystal structures have been solved. One is compact with two of its four peripheral anionic sites (PAS) sterically blocked by complementary subunits. The other is a loose tetramer with all four subunits accessible to solvent. These structures lacked the C-terminal amphipathic t-peptide (WAT domain) that interacts with the proline-rich attachment domain (PRAD). A complete tetramer model (AChEt) was built based on the structure of the PRAD/WAT complex and the compact tetramer. Normal mode analysis suggested that AChEt could exist in several conformations with subunits fluctuating relative to one another. Here, a multiscale simulation involving all-atom molecular dynamics and Cα-based coarse-grained Brownian dynamics simulations was carried out to investigate the large-scale intersubunit dynamics in AChEt. We sampled the ns-μs timescale motions and found that the tetramer indeed constitutes a dynamic assembly of monomers. The intersubunit fluctuation is correlated with the occlusion of the PAS. Such motions of the subunits “gate” ligand-protein association. The gates are open more than 80% of the time on average, which suggests a small reduction in ligand-protein binding. Despite the limitations in the starting model and approximations inherent in coarse graining, these results are consistent with experiments which suggest that binding of a substrate to the PAS is only somewhat hindered by the association of the subunits

A-Site Residues Move Independently from P-Site Residues in all-Atom Molecular Dynamics Simulations of the 70S Bacterial Ribosome

The ribosome is a large macromolecular machine, and correlated motion between residues is necessary for coordinating function across multiple protein and RNA chains. We ran two all-atom, explicit solvent molecular dynamics simulations of the bacterial ribosome and calculated correlated motion between residue pairs by using mutual information. Because of the short timescales of our simulation (ns), we expect that dynamics are largely local fluctuations around the crystal structure. We hypothesize that residues that show coupled dynamics are functionally related, even on longer timescales. We validate our model by showing that crystallographic B-factors correlate well with the entropy calculated as part of our mutual information calculations. We reveal that A-site residues move relatively independently from P-site residues, effectively insulating A-site functions from P-site functions during translation

Molecular dynamics of ribosomal elongation factors G and Tu

Translation on the ribosome is controlled by external factors. During polypeptide lengthening, elongation factors EF-Tu and EF-G consecutively interact with the bacterial ribosome. EF-Tu binds and delivers an aminoacyl-tRNA to the ribosomal A site and EF-G helps translocate the tRNAs between their binding sites after the peptide bond is formed. These processes occur at the expense of GTP. EF-Tu:tRNA and EF-G are of similar shape, share a common binding site, and undergo large conformational changes on interaction with the ribosome. To characterize the internal motion of these two elongation factors, we used 25 ns long all-atom molecular dynamics simulations. We observed enhanced mobility of EF-G domains III, IV, and V and of tRNA in the EF-Tu:tRNA complex. EF-Tu:GDP complex acquired a configuration different from that found in the crystal structure of EF-Tu with a GTP analogue, showing conformational changes in the switch I and II regions. The calculated electrostatic properties of elongation factors showed no global similarity even though matching electrostatic surface patches were found around the domain I that contacts the ribosome, and in the GDP/GTP binding region