Simulation of Complex Biomolecular Systems: The Ribosome Challenge.

Large biomolecular systems are at the heart of many essential cellular processes. The dynamics and energetics of an increasing number of these systems are being studied by computer simulations. Pushing the limits of length- and timescales that can be accessed by current hard- and software has expanded the ability to describe biomolecules at different levels of detail. We focus in this review on the ribosome, which exemplifies the close interplay between experiment and various simulation approaches, as a particularly challenging and prototypic nanomachine that is pivotal to cellular biology due to its central role in translation. We sketch widely used simulation methods and demonstrate how the combination of simulations and experiments advances our understanding of the function of the translation apparatus based on fundamental physics.

Binding of the peptide deformylase on the ribosome surface modulates the exit tunnel interior.

Proteosynthesis on ribosomes is regulated at many levels. Conformational changes of the ribosome, possibly induced by external factors, may transfer over large distances and contribute to the regulation. The molecular principles of this long-distance allostery within the ribosome remain poorly understood. Here, we use structural analysis and atomistic molecular dynamics simulations to investigate peptide deformylase (PDF), an enzyme that binds to the ribosome surface near the ribosomal protein uL22 during translation and chemically modifies the emerging nascent peptide. Our simulations of the entire ribosome-PDF complex reveal that the PDF undergoes a swaying motion on the ribosome surface at the submicrosecond timescale. We show that the PDF affects the conformational dynamics of parts of the ribosome over distances of more than 5 nm. Using a supervised-learning algorithm, we demonstrate that the exit tunnel is influenced by the presence or absence of PDF. Our findings suggest a possible effect of the PDF on the nascent peptide translocation through the ribosome exit tunnel.

The effect of off-center σ -hole on the atom-centered partial charges in halogenated molecules.

Partial atomic charges belong to key concepts of computational chemistry. In some cases, however, they fail in describing the electrostatics of molecules. One such example is the σ -hole, a region of positive electrostatic potential located on halogens and other atoms. In molecular mechanics, the σ -hole is often modeled as a pseudo-atom with a positive partial charge located off the halogen nucleus. Here we address a question, to what extent the pseudo-atom affects partial charges of other atoms in the molecule. To this aim, we have thoroughly analyzed partial charges of over 2300 halogenated molecules from the ZINC database calculated by the restricted electrostatic potential (RESP) method and compared them with the charges fitted by RESP including the pseudo-atom. We show that the pseudo-atom improves charge fitting for a vast majority of molecules. The σ -hole, modeled as the off-center charge, affects the atoms within three covalent bonds from the halogen.

Folding of VemP into translation-arresting secondary structure is driven by the ribosome exit tunnel.

The ribosome is a fundamental biomolecular complex that synthesizes proteins in cells. Nascent proteins emerge from the ribosome through a tunnel, where they may interact with the tunnel walls or small molecules such as antibiotics. These interactions can cause translational arrest with notable physiological consequences. Here, we studied the arrest caused by the regulatory peptide VemP, which is known to form α-helices inside the ribosome tunnel near the peptidyl transferase center under specific conditions. We used all-atom molecular dynamics simulations of the entire ribosome and circular dichroism spectroscopy to study the driving forces of helix formation and how VemP causes the translational arrest. To that aim, we compared VemP dynamics in the ribosome tunnel with its dynamics in solution. We show that the VemP peptide has a low helical propensity in water and that the propensity is higher in mixtures of water and trifluorethanol. We propose that helix formation within the ribosome is driven by the interactions of VemP with the tunnel and that a part of VemP acts as an anchor. This anchor might slow down VemP progression through the tunnel enabling α-helix formation, which causes the elongation arrest.

Antitubercular polyhalogenated phenothiazines and phenoselenazine with reduced binding to CNS receptors.

Targeting energy metabolism in Mycobacterium tuberculosis (Mtb) is a new paradigm in the search for innovative anti-TB drugs. NADH:menaquinone oxidoreductase is a non-proton translocating type II NADH dehydrogenase (NDH-2) that is an essential enzyme in the respiratory chain of Mtb and is not found in mammalian mitochondria. Phenothiazines (PTZs) represent one of the most known class of NDH-2 inhibitors, but their use as anti-TB drugs is currently limited by the wide range of potentially serious off-target effects. In this work, we designed and synthesized a series of new PTZs by decorating the scaffold in an unconventional way, introducing various halogen atoms. By replacing the sulfur atom with selenium, a dibromophenoselenazine 20 was also synthesized. Among the synthesized poly-halogenated PTZs (HPTZs), dibromo and tetrachloro derivatives 9 and 11, along with the phenoselenazine 20, emerged with a better anti-TB profile than the therapeutic thioridazine (TZ). They targeted non-replicating Mtb, were bactericidal, and synergized with rifampin and bedaquiline. Moreover, their anti-TB activity was found to be related to the NDH-2 inhibition. Most important, they showed a markedly reduced affinity to dopaminergic and serotonergic receptors respect to the TZ. From this work emerged, for the first time, as the poly-halogenation of the PTZ core, while permitting to maintain good anti-TB profile could conceivably lead to fewer CNS side-effects risk, making more tangible the use of PTZs for this alternative therapeutic application.

Molecular simulations of the ribosome and associated translation factors.

The ribosome is a macromolecular complex which is responsible for protein synthesis in all living cells according to their transcribed genetic information. Using X-ray crystallography and, more recently, cryo-electron microscopy (cryo-EM), the structure of the ribosome was resolved at atomic resolution in many functional and conformational states. Molecular dynamics simulations have added information on dynamics and energetics to the available structural information, thereby have bridged the gap to the kinetics obtained from single-molecule and bulk experiments. Here, we review recent computational studies that brought notable insights into ribosomal structure and function.

Halogen Bonding in Nucleic Acid Complexes.

Halogen bonding (X-bonding) has attracted notable attention among noncovalent interactions. This highly directional attraction between a halogen atom and an electron donor has been exploited in knowledge-based drug design. A great deal of information has been gathered about X-bonds in protein-ligand complexes, as opposed to nucleic acid complexes. Here we provide a thorough analysis of nucleic acid complexes containing either halogenated building blocks or halogenated ligands. We analyzed close contacts between halogens and electron-rich moieties. The phosphate backbone oxygen is clearly the most common halogen acceptor. We identified 21 X-bonds within known structures of nucleic acid complexes. A vast majority of the X-bonds is formed by halogenated nucleobases, such as bromouridine, and feature excellent geometries. Noncovalent ligands have been found to form only interactions with suboptimal interaction geometries. Hence, the first X-bonded nucleic acid binder remains to be discovered.

Reaction Path Averaging: Characterizing the Structural Response of the DNA Double Helix to Electron Transfer.

A polarizable environment, prominently the solvent, responds to electronic changes in biomolecules rapidly. The knowledge of conformational relaxation of the biomolecule itself, however, may be scarce or missing. In this work, we describe in detail the structural changes in DNA undergoing electron transfer between two adjacent nucleobases. We employ an approach based on averaging of tens to hundreds of thousands of nonequilibrium trajectories generated with molecular dynamics simulation, and a reduction of dimensionality suitable for DNA. We show that the conformational response of the DNA proceeds along a single collective coordinate that represents the relative orientation of two consecutive base pairs, namely, a combination of helical parameters shift and tilt. The structure of DNA relaxes on time scales reaching nanoseconds, contributing marginally to the relaxation of energies, which is dominated by the modes of motion of the aqueous solvent. The concept of reaction path averaging (RPA), conveniently exploited in this context, makes it possible to filter out any undesirable noise from the nonequilibrium data, and is applicable to any chemical process in general.

Computer Modeling of Halogen Bonds and Other σ-Hole Interactions.

In the field of noncovalent interactions a new paradigm has recently become popular. It stems from the analysis of molecular electrostatic potentials and introduces a label, which has recently attracted enormous attention. The label is σ-hole, and it was first used in connection with halogens. It initiated a renaissance of interest in halogenated compounds, and later on, when found also on other groups of atoms (chalcogens, pnicogens, tetrels and aerogens), it resulted in a new direction of research of intermolecular interactions. In this review, we summarize advances from about the last 10 years in understanding those interactions related to σ-hole. We pay particular attention to theoretical and computational techniques, which play a crucial role in the field.

Polar flattening and the strength of halogen bonding.

The effect of polar flattening on the stability of 32 halogen-bonded complexes was investigated by utilizing CCSD(T)/CBS, DFT, and DFT-SAPT/CBS methods. It is shown that the value of polar flattening increases with the decreasing value of studied isodensity. For the complexes investigated, the polar flattening based on the isodensity of 0.001 au reaches 0.2-0.3 Å and 10-15% in absolute and relative values, respectively. These geometrical changes induce differences in the stabilization energy up to 20%.

Substrate Specificity, Inhibitor Selectivity and Structure-Function Relationships of Aldo-Keto Reductase 1B15: A Novel Human Retinaldehyde Reductase.

Human aldo-keto reductase 1B15 (AKR1B15) is a newly discovered enzyme which shares 92% amino acid sequence identity with AKR1B10. While AKR1B10 is a well characterized enzyme with high retinaldehyde reductase activity, involved in the development of several cancer types, the enzymatic activity and physiological role of AKR1B15 are still poorly known. Here, the purified recombinant enzyme has been subjected to substrate specificity characterization, kinetic analysis and inhibitor screening, combined with structural modeling. AKR1B15 is active towards a variety of carbonyl substrates, including retinoids, with lower kcat and Km values than AKR1B10. In contrast to AKR1B10, which strongly prefers all-trans-retinaldehyde, AKR1B15 exhibits superior catalytic efficiency with 9-cis-retinaldehyde, the best substrate found for this enzyme. With ketone and dicarbonyl substrates, AKR1B15 also shows higher catalytic activity than AKR1B10. Several typical AKR inhibitors do not significantly affect AKR1B15 activity. Amino acid substitutions clustered in loops A and C result in a smaller, more hydrophobic and more rigid active site in AKR1B15 compared with the AKR1B10 pocket, consistent with distinct substrate specificity and narrower inhibitor selectivity for AKR1B15.

The Effect of Halogen-to-Hydrogen Bond Substitution on Human Aldose Reductase Inhibition.

The effect of halogen-to-hydrogen bond substitution on the binding energetics and biological activity of a human aldose reductase inhibitor has been studied using X-ray crystallography, IC50 measurements, advanced binding free energy calculations, and simulations. The replacement of Br or I atoms by an amine (NH2) group has not induced changes in the original geometry of the complex, which made it possible to study the isolated features of selected noncovalent interactions in a biomolecular complex.

Characteristics of a σ-Hole and the Nature of a Halogen Bond.

The nature of halogen bonding in 128 complexes was investigated using advanced quantum mechanical calculations. First, isolated halogen donors were studied and their σ-holes were described in terms of size and magnitude. Later, both partners in the complex were considered and their interaction was described in terms of DFT-SAPT decomposition. The whole set of complexes under study was split into two categories on the basis of their stabilisation energy. The first subset with 38 complexes possesses stabilisation energies in the range 7-32 kcal/mol, while the second subset with 90 complexes has stabilisation energies smaller than 7 kcal/mol. The first subset is characterised by small intermolecular distances (less than 2.5 Å) and a significant contraction of van der Waals (vdW) distance (sum of vdW radii). Here the polarisation/electrostatic energy is dominant, mostly followed by induction and dispersion energies. The importance of induction energy reflects the charge-transfer character of the respective halogen bonds. Intermolecular distances in the second subset are large and the respective contraction of vdW distance upon the formation of a halogen bond is much smaller. Here the dispersion energy is mostly dominant, followed by polarisation and induction energies. Considering the whole set of complexes, we conclude that the characteristic features of their halogen bonds arise from the concerted action of polarisation and dispersion energies and neither of these energies can be considered as dominant. Finally, the magnitude of the σ-hole and DFT-SAPT stabilisation energy correlates only weakly within the whole set of complexes.

Statistical analysis of σ-holes: a novel complementary view on halogen bonding.

To contribute to the understanding of noncovalent binding of halogenated molecules with biological activity, electrostatic potential (ESP) maps of more than 2500 compounds were thoroughly analysed. A peculiar region of positive ESP, called the σ-hole, is a concept of central importance for halogen bonding. We aim to simplify the view on σ-holes and to provide general trends in organic drug-like molecules. The results are in fair agreement with crystallographic surveys of small molecules as well as biomolecular complexes and may help to improve the intuition of chemists when dealing with halogenated compounds.