Joanna Trylska

Binding Pathways of Ligands to HIV-1 Protease: Coarse-grained and Atomistic Simulations

Multiscale simulations (coarse‐grained Brownian dynamics simulations and all‐atom molecular dynamics simulations in implicit solvent) were applied to reveal the binding processes of ligands as they enter the binding site of the HIV‐1 protease. The initial structures used for the molecular dynamics simulations were generated based on the Brownian dynamics trajectories, and this is the first molecular dynamics simulation of modeling the association of a ligand with the protease. We found that a protease substrate successfully binds to the protein when the flaps are fully open. Surprisingly, a smaller cyclic urea inhibitor (XK263) can reach the binding site when the flaps are not fully open. However, if the flaps are nearly closed, the inhibitor must rearrange or binding can fail because the inhibitor cannot attain proper conformations to enter the binding site. Both the peptide substrate and XK263 can also affect the proteins internal motion, which may help the flaps to open. Simulations allow us to efficiently study the ligand binding processes and may help those who study drug discovery to find optimal association pathways and to design those ligands with the best binding kinetics.

Coarse‐grained force field for the nucleosome from self‐consistent multiscaling

A coarse‐grained simulation model for the nucleosome is developed, using a methodology modified from previous work on the ribosome. Protein residues and DNA nucleotides are represented as beads, interacting through harmonic (for neighboring) or Morse (for nonbonded) potentials. Force‐field parameters were estimated by Boltzmann inversion of the corresponding radial distribution functions obtained from a 5‐ns all‐atom molecular dynamics (MD) simulation, and were refined to produce agreement with the all‐atom MD simulation. This self‐consistent multiscale approach yields a coarse‐grained model that is capable of reproducing equilibrium structural properties calculated from a 50‐ns all‐atom MD simulation. This coarse‐grained model speeds up nucleosome simulations by a factor of 103 and is expected to be useful in examining biologically relevant dynamical nucleosome phenomena on the microsecond timescale and beyond.

The role of hydrogen bonding in the enzymatic reaction catalyzed by HIV-1 protease.

The hydrogen‐bond network in various stages of the enzymatic reaction catalyzed by HIV‐1 protease was studied through quantum‐classical molecular dynamics simulations. The approximate valence bond method was applied to the active site atoms participating directly in the rearrangement of chemical bonds. The rest of the protein with explicit solvent was treated with a classical molecular mechanics model. Two possible mechanisms were studied, general‐acid/general‐base (GA/GB) with Asp 25 protonated at the inner oxygen, and a direct nucleophilic attack by Asp 25. Strong hydrogen bonds leading to spontaneous proton transfers were observed in both reaction paths. A single‐well hydrogen bond was formed between the peptide nitrogen and outer oxygen of Asp 125. The proton was diffusely distributed with an average central position and transferred back and forth on a picosecond scale. In both mechanisms, this interaction helped change the peptide‐bond hybridization, increased the partial charge on peptidyl carbon, and in the GA/GB mechanism, helped deprotonate the water molecule. The inner oxygens of the aspartic dyad formed a low‐barrier, but asymmetric hydrogen bond; the proton was not positioned midway and made a slightly elongated covalent bond, transferring from one to the other aspartate. In the GA/GB mechanism both aspartates may help deprotonate the water molecule. We observed the breakage of the peptide bond and found that the protonation of the peptidyl amine group was essential for the peptide‐bond cleavage. In studies of the direct nucleophilic mechanism, the peptide carbon of the substrate and oxygen of Asp 25 approached as close as 2.3 Å.

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.

Molecular dynamics study of the ribosomal A-site.

Many aminoglycosidic antibiotics target the A-site of 16S RNA in the small ribosomal subunit and affect the fidelity of protein translation in bacteria. Upon binding, aminoglycosides displace two adenines (A1492 and A1493 for E. coli numbering) that are involved in tRNA anticodon loop recognition. The major difference in the aminoglycosidic binding site between the prokaryota and eukaryota is an adenine into guanine substitution in the position 1408. This mutation likely affects the dynamics of near A1492 and A1493 and hinders the binding of aminoglycosides to eukaryotic ribosomes. With multiple 20 ns long all-atom molecular dynamics simulations, we study the flexibility of a 22 nucleotide RNA fragment which mimics the aminoglycosidic binding site. Simulations are carried out for both native and A1408G mutated RNA as well as for their complexes with aminoglycosidic representative paromomycin. We observe intra- and extrahelical configurations of A1492 and A1493, which differ between the prokaryotic and the mutated structure. We obtain configurations of the A-site that are also observed in the NMR and crystal structures. Our studies show the differences in the internal mobility of the A-site, as well as that in ion and water density distributions inside of the binding cleft, between the prokaryotic and mutated RNA. We also compare the performance of two force field parameters for RNA, Amber and Charmm.

Coarse-grained models to study dynamics of nanoscale biomolecules and their applications to the ribosome

Biopolymers are of dynamic nature and undergo functional motions spanning a large spectrum of timescales. To study the internal dynamics of nano-sized molecular complexes that exceed hundred thousands of atoms with atomic detail is computationally inefficient. Therefore, to achieve both the spatial and temporal scales of biological interest coarse-grained models of macromolecules are often used. By uniting groups of atoms into single interacting centers one decreases the resolution of the system and gets rid of the irrelevant degrees of freedom. This simplification, even though it requires parameterization, makes the studies of biomolecular dynamics computationally tractable and allows us to reach beyond the microsecond time frame. Here, I review the coarse-grained models of macromolecules composed of proteins and nucleic acids. I give examples of one-bead models that were developed to investigate the internal dynamics and focus on their applications to the ribosome--the nanoscale protein synthesis machine.

Brownian dynamics simulations on CPU and GPU with BD_BOX

There has been growing interest in simulating biological processes under in vivo conditions due to recent advances in experimental techniques dedicated to study single particle behavior in crowded environments. We have developed a software package, BD_BOX, for multiscale Brownian dynamics simulations. BD_BOX can simulate either single molecules or multicomponent systems of diverse, interacting molecular species using flexible, coarse‐grained bead models. BD_BOX is written in C and employs modern computer architectures and technologies; these include MPI for distributed‐memory architectures, OpenMP for shared‐memory platforms, NVIDIA CUDA framework for GPGPU, and SSE vectorization for CPU.

Secondary structures of native and pathogenic huntingtin N-terminal fragments.

Huntingtons disease is a neurodegenerative disorder caused by a polyglutamine (polyQ) expansion in the N-terminal fragment of the Huntingtin (Htt) protein. Structural properties of Htt N-terminal regions and the molecular mechanism leading to protein aggregation have not been fully explained yet. We performed all-atom replica exchange molecular dynamics to investigate the structures of Htt N-terminal parts with polyQ tracts of nonpathogenic and pathogenic lengths. The monomers were composed of the headpiece (17 N-terminal residues), a polyQ tract (polyQ(17) for native and polyQ(55) for pathogenic sequence), and a polyP(11) region, followed by 17 amino acids of mixed sequence. We found that corresponding regions in both fragments fold to similar secondary structures; the headpiece and polyQ stretch adopt mainly α-helical conformations, and polyP(11) forms the PP II-type helix. The native N-terminal fragment is more compact and stabilized by hydrophobic interactions between the surface of polyP(11) and the amphipathic helix of the headpiece. In the pathogenic fragment the headpiece is solvent exposed and does not interact with polyP(11). The predicted structure of the native N-terminal fragment agrees with the X-ray structure of the Htt first exon containing polyQ(17). The structure of the pathogenic fragment adheres to an aggregation model that is mediated by the Htt headpiece.

Understanding the origins of bacterial resistance to aminoglycosides through molecular dynamics mutational study of the ribosomal A-site.

Paromomycin is an aminoglycosidic antibiotic that targets the RNA of the bacterial small ribosomal subunit. It binds in the A-site, which is one of the three tRNA binding sites, and affects translational fidelity by stabilizing two adenines (A1492 and A1493) in the flipped-out state. Experiments have shown that various mutations in the A-site result in bacterial resistance to aminoglycosides. In this study, we performed multiple molecular dynamics simulations of the mutated A-site RNA fragment in explicit solvent to analyze changes in the physicochemical features of the A-site that were introduced by substitutions of specific bases. The simulations were conducted for free RNA and in complex with paromomycin. We found that the specific mutations affect the shape and dynamics of the binding cleft as well as significantly alter its electrostatic properties. The most pronounced changes were observed in the U1406C∶U1495A mutant, where important hydrogen bonds between the RNA and paromomycin were disrupted. The present study aims to clarify the underlying physicochemical mechanisms of bacterial resistance to aminoglycosides due to target mutations.

RedMD—Reduced molecular dynamics package

We developed a software package (RedMD) to perform molecular dynamics simulations and normal mode analysis of reduced models of proteins, nucleic acids, and their complexes. With RedMD one can perform molecular dynamics simulations in a microcanonical ensemble, with Berendsen and Langevin thermostats, and with Brownian dynamics. We provide force field and topology generators which are based on the one‐bead per residue/nucleotide elastic network model and its extensions. The user can change the force field parameters with the command line options that are passed to generators. Also, the generators can be modified, for example, to add new potential energy functions. Normal mode analysis tool is available for elastic or anisotropic network models. The program is written in C and C++ languages and the structure/topology of a molecule is based on an XML format. OpenMP technology for shared‐memory architectures was used for code parallelization. The code is distributed under GNU public licence and available at http://bionano.icm.edu.pl/software/.