Biomedical subjects

Joanna Trylska

Publications and source records attributed to Joanna Trylska.

11 recordsLinked to original sources

PubMed2006

Binding of aminoglycosidic antibiotics to the oligonucleotide A-site model and 30S ribosomal subunit: Poisson-Boltzmann model, thermal denaturation, and fluorescence studies.

The binding of paromomycin and similar antibiotics to the oligonucleotide A-site model and the small (30S) ribosomal subunit has been studied using continuum electrostatics methods. Crystallographic information from complexes of paromomycin, tobramycin, and Geneticin bound to an A-site oligonucleotide, and paromomycin and streptomycin complexed to the 30S subunit was used as a foundation to develop structures of similar antibiotics in the same ribosomal binding site. Relative binding free energies were calculated by combining the electrostatic contribution, which was obtained by solving the Poisson-Boltzmann equation, with a surface-area-dependent apolar term and contributions from conformational changes. These computed results showed good correlation with the experimental data resulting from fluorescence binding assays and thermal denaturation studies, demonstrating the ability of the Poisson-Boltzmann model to provide insight into the electrostatic mechanisms for aminoglycoside binding and direction for designing more effective antibiotics.

Aminoglycosides↗

PubMed2006

Flap opening dynamics in HIV-1 protease explored with a coarse-grained model.

We present a one-bead coarse-grained model that enables dynamical simulations of proteins on the time scale of tens of microseconds. The parameterization of the force field includes accurate conformational terms that allow for fast and reliable exploration of the configurational space. The model is applied to the dynamics of flap opening in HIV-1 protease. The experimental structure of the recently crystallized semi-open conformation of HIV-1 protease is well reproduced in the simulation, which supports the accuracy of our model. Thanks to very long simulations and extensive sampling of opening and closing events, we also investigate the thermodynamics and kinetics of the opening process. We have shown that the effect of the solvent slows down the dynamics to the experimentally observed time scales. The model is found to be reliable for application to substrate docking simulations, which are currently in progress.

HIV Protease↗

PubMed2006

Electrostatic properties of cowpea chlorotic mottle virus and cucumber mosaic virus capsids.

Electrostatic properties of cowpea chlorotic mottle virus (CCMV) and cucumber mosaic virus (CMV) were investigated using numerical solutions to the Poisson-Boltzmann equation. Experimentally, it has been shown that CCMV particles swell in the absence of divalent cations when the pH is raised from 5 to 7. CMV, although structurally homologous, does not undergo this transition. An analysis of the calculated electrostatic potential confirms that a strong electrostatic repulsion at the calcium-binding sites in the CCMV capsid is most likely the driving force for the capsid swelling process during the release of calcium. The binding interaction between the encapsulated genome material (RNA) inside of the capsid and the inner capsid shell is weakened during the swelling transition. This probably aids in the RNA release process, but it is unlikely that the RNA is released through capsid openings due to unfavorable electrostatic interaction between the RNA and capsid inner shell residues at these openings. Calculations of the calcium binding energies show that Ca(2+) can bind both to the native and swollen forms of the CCMV virion. Favorable binding to the swollen form suggests that Ca(2+) ions can induce the capsid contraction and stabilize the native form.

Binding Sites↗

PubMed2006

The influence of macromolecular crowding on HIV-1 protease internal dynamics.

High macromolecular concentrations, or crowded conditions, have been shown to affect a wide variety of molecular processes, including diffusion, association and dissociation, and protein folding and stability. Here, we model the effect of macromolecular crowding on the internal dynamics of a protein, HIV-1 protease, using Brownian dynamics simulations. HIV-1 protease possesses a pair of flaps which are postulated to open in the early stages of its catalytic mechanism. Compared to low concentrations, close-packed concentrations of repulsive crowding agents are found to significantly reduce the fraction of time that the protease flaps are open. Macromolecular crowding is likely to have a major effect on in vivo enzyme activity, and may play an important regulatory role in the viral life cycle.

Amino Acids↗

PubMed2006

Gated binding of ligands to HIV-1 protease: Brownian dynamics simulations in a coarse-grained model.

The internal motions of proteins may serve as a "gate" in some systems, which controls ligand-protein association. This study applies Brownian dynamics simulations in a coarse-grained model to study the gated association rate constants of HIV-1 proteases and drugs. The computed gated association rate constants of three protease mutants, G48V/V82A/I84V/L90M, G48V, and L90M with three drugs, amprenavir, indinavir, and saquinavir, yield good agreements with experiments. The work shows that the flap dynamics leads to "slow gating". The simulations suggest that the flap flexibility and the opening frequency of the wild-type, the G48V and L90M mutants are similar, but the flaps of the variant G48V/V82A/I84V/L90M open less frequently, resulting in a lower gated rate constant. The developed methodology is fast and provides an efficient way to predict the gated association rate constants for various protease mutants and ligands.

Carbamates↗

PubMed2006

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.

Computational Biology↗

PubMed2005

Exploring assembly energetics of the 30S ribosomal subunit using an implicit solvent approach.

To explore the relationship between the assembly of the 30S ribosomal subunit and interactions among the constituent components, 16S RNA and proteins, relative binding free energies of the T. thermophilus 30S proteins to the 16S RNA were studied based on an implicit solvent model of electrostatic, nonpolar, and entropic contributions. The late binding proteins in our assembly map were found not to bind to the naked 16S RNA. The 5' domain early kinetic class proteins, on average, carry the highest positive charge, get buried the most upon binding to 16S RNA, and show the most favorable binding. Some proteins (S10/S14, S6/S18, S13/S19) have more stabilizing interactions while binding as dimers. Our computed assembly map resembles that of E. coli; however, the central domain path is more similar to that of A. aeolicus, a hyperthermophilic bacteria.

Entropy↗

PubMed2005

Exploring global motions and correlations in the ribosome.

We studied slower global coupled motions of the ribosome with half a microsecond of coarse-grained molecular dynamics. A low-resolution anharmonic network model that allows for the evolution of tertiary structure and long-scale sampling was developed and parameterized. Most importantly, we find that functionally important movements of L7/L12 and L1 lateral stalks are anticorrelated. Other principal directions of motions include widening of the tRNA cleft and the rotation of the small subunit which occurs as one block and is in phase with the movement of L1 stalk. The effect of the dynamical correlation pattern on the elongation process is discussed. Small fluctuations of the 3' tRNA termini and anticodon nucleotides show tight alignment of substrates for the reaction. Our model provides an efficient and reliable way to study the dynamics of large biomolecular systems composed of both proteins and nucleic acids.

Anticodon↗

PubMed2004

Ribosome motions modulate electrostatic properties.

The electrostatic properties of the 70S ribosome of Thermus thermophilus were studied qualitatively by solving the Poisson-Boltzmann (PB) equation in aqueous solution and with physiological ionic strength. The electrostatic potential was calculated for conformations of the ribosome derived by recent normal mode analysis (Tama, F., et al. Proc Natl Acad Sci USA 2003 100, 9319-9323) of the ratchet-like reorganization that occurs during translocation (Frank, J.; Agrawal, R. K. Nature 2000 406, 318-322). To solve the PB equation, effective parameters (charges and radii), applicable to a highly charged backbone model of the ribosome, were developed. Regions of positive potential were found at the binding site of the elongation factors G and Tu, as well as where the release factors bind. Large positive potential areas are especially pronounced around the L11 and L6 proteins. The region around the L1 protein is also positively charged, supporting the idea that L1 may interact with the E-site tRNA during its release from the ribosome after translocation. Functional rearrangement of the ribosome leads to electrostatic changes which may help the translocation of the tRNAs during the elongation stage.

Binding Sites↗

PubMed2004

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 A.

Binding Sites↗

PubMed2002

Molecular dynamics simulations of the first steps of the reaction catalyzed by HIV-1 protease.

The mechanism of the first steps of the reaction catalyzed by HIV-1 protease was studied through molecular dynamics simulations. The potential energy surface in the active site was generated using the approximate valence bond method. The approximate valence bond (AVB) method was parameterized based on density functional calculations. The surrounding protein and explicit water environment was modeled with conventional, classical force field. The calculations were performed based on HIV-1 protease complexed with the MVT-101 inhibitor that was modified to a model substrate. The protonation state of the catalytic aspartates was determined theoretically. Possible reaction mechanisms involving the lytic water molecule are accounted for in this study. The modeled steps include the dissociation of the lytic water molecule and proton transfer onto Asp-125, the nucleophilic attack followed by a proton transfer onto peptide nitrogen. The simulations show that in the active site most preferable energetically are structures consisting of ionized or polarized molecular fragments that are not accounted for in conventional molecular dynamics. The mobility of the lytic water molecule, the dynamics of the hydrogen bond network, and the conformation of the aspartates in the active center were analyzed.

Aspartic Acid↗