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

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

10.3389/fmolb.2021.797754Frontiers in Molecular Biosciences879775

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

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

Non-Bulk-Like Solvent Behavior in the Ribosome Exit Tunnel

As nascent proteins are synthesized by the ribosome, they depart via an exit tunnel running through the center of the large subunit. The exit tunnel likely plays an important part in various aspects of translation. Although water plays a key role in many bio-molecular processes, the nature of water confined to the exit tunnel has remained unknown. Furthermore, solvent in biological cavities has traditionally been characterized as either a continuous dielectric fluid, or a discrete tightly bound molecule. Using atomistic molecular dynamics simulations, we predict that the thermodynamic and kinetic properties of water confined within the ribosome exit tunnel are quite different from this simple two-state model. We find that the tunnel creates a complex microenvironment for the solvent resulting in perturbed rotational dynamics and heterogenous dielectric behavior. This gives rise to a very rugged solvation landscape and significantly retarded solvent diffusion. We discuss how this non-bulk-like solvent is likely to affect important biophysical processes such as sequence dependent stalling, co-translational folding, and antibiotic binding. We conclude with a discussion of the general applicability of these results to other biological cavities

Quantum mechanical studies of lincosamides

Lincosamides are a class of antibiotics used both in clinical and veterinary practice for a wide range of pathogens. This group of drugs inhibits the activity of the bacterial ribosome by binding to the 23S RNA of the large ribosomal subunit and blocking protein synthesis. Currently, three X-ray structures of the ribosome in complex with clindamycin are available in the Protein Data Bank, which reveal that there are two distinct conformations of the pyrrolidinyl propyl group of the bound clindamycin. In this work, we used quantum mechanical methods to investigate the probable conformations of clindamycin in order to explain the two binding modes in the ribosomal 23S RNA. We studied three lincosamide antibiotics: clindamycin, lincomycin, and pirlimycin at the B3LYP level with the 6-31G** basis set. The focus of our work was to connect the conformational landscape and electron densities of the two clindamycin conformers found experimentally with their physicochemical properties. For both functional conformers, we applied natural bond orbital (NBO) analysis and the atoms in molecules (AIM) theory, and calculated the NMR parameters. Based on the results obtained, we were able to show that the structure with the intramolecular hydrogen bond C=O…H–O is the most stable conformer of clindamycin. The charge transfer between the pyrrolidine-derivative ring and the six-atom sugar (methylthiolincosamide), which are linked via an amide bond, was found to be the dominant factor influencing the high stability of this conformer

Inhibitory Potential of Prodomain of Plasmodium falciparum Protease Serine Repeat Antigen 5 for Asexual Blood Stages of Parasite

Plasmodium falciparum serine repeat antigen 5 (SERA5) is a target for both drug and vaccine intervention against malaria. SERA5 is secreted in the parasitophorous vacuole where it is proteolytically processed before schizont rupture. Among the processed products is a 50.8-kDa central domain of the protease, which possesses chymotrypsin-like activity and consists of a 28.9-kDa catalytic domain with a 21.9-kDa N-terminal prodomain, which remain attached together. Because SERA5 has been implicated in merozoite egress from host erythrocytes, the effect of the prodomain and a heptapeptide derived from its C-terminus spanning from D560 to F566 (DNSDNMF) on parasite growth was studied. When E. coli-expressed prodomain was incubated with parasite culture, a significant delay in transition from schizont to ring stages was observed up to nanomolar concentrations. The peptide, DNSDNMF also showed similar effects but at nearly 1000-fold higher concentrations. The peptide was also found to interact with the catalytic domain. These data demonstrate the crucial role of SERA5 prodomain for the egress process. Given the inhibitory potential of the prodomain for the parasite, we suggest that peptidomimetic inhibitors based on SERA5 prodomain sequences can be developed as future therapeutics against malaria

A Computational Investigation on the Connection between Dynamics Properties of Ribosomal Proteins and Ribosome Assembly

Assembly of the ribosome from its protein and RNA constituents has been studied extensively over the past 50 years, and experimental evidence suggests that prokaryotic ribosomal proteins undergo conformational changes during assembly. However, to date, no studies have attempted to elucidate these conformational changes. The present work utilizes computational methods to analyze protein dynamics and to investigate the linkage between dynamics and binding of these proteins during the assembly of the ribosome. Ribosomal proteins are known to be positively charged and we find the percentage of positive residues in r-proteins to be about twice that of the average protein: Lys+Arg is 18.7% for E. coli and 21.2% for T. thermophilus. Also, positive residues constitute a large proportion of RNA contacting residues: 39% for E. coli and 46% for T. thermophilus. This affirms the known importance of charge-charge interactions in the assembly of the ribosome. We studied the dynamics of three primary proteins from E. coli and T. thermophilus 30S subunits that bind early in the assembly (S15, S17, and S20) with atomic molecular dynamic simulations, followed by a study of all r-proteins using elastic network models. Molecular dynamics simulations show that solvent-exposed proteins (S15 and S17) tend to adopt more stable solution conformations than an RNA-embedded protein (S20). We also find protein residues that contact the 16S rRNA are generally more mobile in comparison with the other residues. This is because there is a larger proportion of contacting residues located in flexible loop regions. By the use of elastic network models, which are computationally more efficient, we show that this trend holds for most of the 30S r-proteins

Visualization of Early Events in Acetic Acid Denaturation of HIV-1 Protease: A Molecular Dynamics Study

Protein denaturation plays a crucial role in cellular processes. In this study, denaturation of HIV-1 Protease (PR) was investigated by all-atom MD simulations in explicit solvent. The PR dimer and monomer were simulated separately in 9 M acetic acid (9 M AcOH) solution and water to study the denaturation process of PR in acetic acid environment. Direct visualization of the denaturation dynamics that is readily available from such simulations has been presented. Our simulations in 9 M AcOH reveal that the PR denaturation begins by separation of dimer into intact monomers and it is only after this separation that the monomer units start denaturing. The denaturation of the monomers is flagged off by the loss of crucial interactions between the α-helix at C-terminal and surrounding β-strands. This causes the structure to transit from the equilibrium dynamics to random non-equilibrating dynamics. Residence time calculations indicate that denaturation occurs via direct interaction of the acetic acid molecules with certain regions of the protein in 9 M AcOH. All these observations have helped to decipher a picture of the early events in acetic acid denaturation of PR and have illustrated that the α-helix and the β-sheet at the C-terminus of a native and functional PR dimer should maintain both the stability and the function of the enzyme and thus present newer targets for blocking PR function