Use of the parmbsc0 force field and trajectory analysis to study the binding of netropsin to the DNA fragment (5'CCAATTGG)(2) in the presence of excess NaCl salt in aqueous solution

The parmbsc0 force field was applied to study in detail the binding of netropsin, at a salt concentration of 0.28M Na(+), to the minor groove of an 8-mer (5'CCAATTGG)(2) DNA duplex forming a netropsin·DNA complex which previously has been characterized by X-ray crystallography, albeit with the use of closely related DNA duplexes. The X-ray structure revealed that the terminal guanidinium and amidinium groups of netropsin interact with the extreme ends of the palindromic AATT sequence of the receptor DNA. The parmbsc0 parameters of B-DNA and AMBER v9 parameters of netropsin generated a stable 6ns molecular dynamics (MD) trajectory for a 1:1 class I binding motif of this complex. Trajectory analysis for the salt and hydration effects on the binding of netropsin to the 8-mer DNA duplex revealed that 18 water molecules and 2 Na(+) are displaced from the DNA upon netropsin binding. A hydration density map of the complex parallels the X-ray data showing that two structured water molecules are localized near the netropsin guanidinium and amidinium groups forming H-bond bridges between the receptor and the ligand.     

The ABCs of molecular dynamics simulations on B-DNA, circa 2012

This article provides a retrospective on the ABC initiative in the area of all-atom molecular dynamics (MD) simulations including explicit solvent on all tetranucleotide steps of duplex B-form DNA duplex, ca. 2012. The ABC consortium has completed two phases of simulations, the most current being a set of 50-100 trajectories based on the AMBER ff99 force field together with the parmbsc0 modification. Some general perspectives on the field of MD on DNA and sequence effects on DNA structure are provided, followed by an overview our MD results, including a detailed comparison of the ff99/parmbsc0 results with crystal and NMR structures available for d(CGCGAATTCGCG). Some projects inspired by or related to the ABC initiative and database are also reviewed, including methods for the trajectory analyses, informatics of dealing with the large database of results, compressions of trajectories for efficacy of distribution, DNA solvation by water and ions, parameterization of coarse-grained models with applications and gene finding and genome annotation.     

Refinement of the AMBER force field for nucleic acids: improving the description of alpha/gamma conformers

We present here the parmbsc0 force field, a refinement of the AMBER parm99 force field, where emphasis has been made on the correct representation of the alpha/gamma concerted rotation in nucleic acids (NAs). The modified force field corrects overpopulations of the alpha/gamma = (g+,t) backbone that were seen in long (more than 10 ns) simulations with previous AMBER parameter sets (parm94-99). The force field has been derived by fitting to high-level quantum mechanical data and verified by comparison with very high-level quantum mechanical calculations and by a very extensive comparison between simulations and experimental data. The set of validation simulations includes two of the longest trajectories published to date for the DNA duplex (200 ns each) and the largest variety of NA structures studied to date (15 different NA families and 97 individual structures). The total simulation time used to validate the force field includes near 1 mus of state-of-the-art molecular dynamics simulations in aqueous solution.     

Towards a molecular dynamics consensus view of B-DNA flexibility

We present a systematic study of B-DNA flexibility in aqueous solution using long-scale molecular dynamics simulations with the two more recent versions of nucleic acids force fields (CHARMM27 and parmbsc0) using four long duplexes designed to contain several copies of each individual base pair step. Our study highlights some differences between pambsc0 and CHARMM27 families of simulations, but also extensive agreement in the representation of DNA flexibility. We also performed additional simulations with the older AMBER force fields parm94 and parm99, corrected for non-canonical backbone flips. Taken together, the results allow us to draw for the first time a consensus molecular dynamics picture of B-DNA flexibility.     

A modeled structure for amidase-03 from Bacillus anthracis

Homology models of amidase-03 from Bacillus anthracis were constructed using Modeller (9v2). Modeller constructs protein models using an automated approach for comparative protein structure modeling by the satisfaction of spatial restraints. A template structure of Listeria monocytogenes bacteriophage PSA endolysin PlyPSA (PDB ID: 1XOV) was selected from protein databank (PDB) using BLASTp with BLOSUM62 sequence alignment scoring matrix. We generated five models using the Modeller default routine in which initial coordinates are randomized and evaluated by pseudo-energy parameters. The protein models were validated using PROCHECK and energy minimized using the steepest descent method in GROMACS 3.2 (flexible SPC water model in cubic box of size 1 Å instead of rigid SPC model). We used G43a1 force field in GROMACS for energy calculations and the generated structure was subsequently analyzed using the VMD software for stereo-chemistry, atomic clash and misfolding. A detailed analysis of the amidase-03 model structure from Bacillus anthracis will provide insight to the molecular design of suitable inhibitors as drug candidates.     

Quantifying Molecular Forces with Serially Connected Force Sensors

To understand the mechanical forces involved in cell adhesion, molecular force sensors have been developed to study tension through adhesion proteins. Recently, a class of molecular force sensors called tension gauge tethers (TGTs) have been developed that rely on irreversible force-dependent dissociation of a DNA duplex to study cell adhesion forces. Although the TGT offers a high signal-to-noise ratio and is ideal for studying fast/single-molecular adhesion processes, quantitative interpretation of experimental results has been challenging. Here, we use a computational approach to investigate how TGT fluorescence readout can be quantitatively interpreted. In particular, we studied force sensors made of a single TGT, multiplexed single TGTs, and two TGTs connected in series. Our results showed that fluorescence readout using a single TGT can result from drastically different combinations of force history and adhesion event density that span orders of magnitude. In addition, the apparent behavior of the TGT is influenced by the tethered receptor-ligand, making it necessary to calibrate the TGT with every new receptor-ligand. To solve this problem, we proposed a system of two serially connected TGTs. Our result shows that not only is the ratiometric readout of serial TGT independent of the choice of receptor-ligand, it is able to reconstruct force history with sub-pN force resolution. This is also not possible by simply multiplexing different types of TGTs together. Last, we systematically investigated how the sequence composition of the two serially connected TGTs can be tuned to achieve different dynamic range. This computational study demonstrated how serially connected irreversible molecular dissociation processes can accurately quantify molecular force and laid the foundation for subsequent experimental studies.     

Local and global effects of strong DNA bending induced during molecular dynamics simulations

DNA bending plays an important role in many biological processes, but its molecular and energetic details as a function of base sequence remain to be fully understood. Using a recently developed restraint, we have studied the controlled bending of four different B-DNA oligomers using molecular dynamics simulations. Umbrella sampling with the AMBER program and the recent parmbsc0 force field yield free energy curves for bending. Bending 15-base pair oligomers by 90° requires roughly 5kcalmol⁻¹, while reaching 150° requires of the order of 12kcalmol⁻¹. Moderate bending occurs mainly through coupled base pair step rolls. Strong bending generally leads to local kinks. The kinks we observe all involve two consecutive base pair steps, with disruption of the central base pair (termed Type II kinks in earlier work). A detailed analysis of each oligomer shows that the free energy of bending only varies quadratically with the bending angle for moderate bending. Beyond this point, in agreement with recent experiments, the variation becomes linear. An harmonic analysis of each base step yields force constants that not only vary with sequence, but also with the degree of bending. Both these observations suggest that DNA is mechanically more complex than simple elastic rod models would imply.     

Cysteine 265 Is in the Active Site of, But Is Not Essential for Catalysis by tRNA-Guanine Transglycosylase (TGT) from Escherichia coli

Site-directed mutagenesis and X-ray absorption spectroscopy studies have previously shown that the tRNA-guanine transglycosylase (TGT) from Escherichia coli is a zinc metalloprotein and identified the enzymic ligands to the zinc [Chong et al. (1995), Biochemistry 34, 3694–3701; Garcia et al. (1966), Biochemistry 35, 3133–3139]. During these studies one mutant, TGT (C265A), was found to exhibit a significantly lower specific activity, but was not found to be involved in the zinc site. The present report demonstrates that TGT is inactivated by treatment with thiol reagents (e.g., DTNB, MMTS, and N-ethylmaleimide). Further, this inactivation is shown to be due to modification of cysteine 265. The kinetic parameters for the mutants TGT (C265A) and TGT (C265S), however, suggest that this residue is not performing a critical role in the TGT reaction. We conclude that cysteine 265 is in the active site of TGT, but is not performing a critical catalytic function. This conclusion is supported by the recent determination of the X-ray crystal structure of the TGT from Zymomonas mobilis [Romier et al. (1966), EMBO J. 15, 2850–2857], which reveals that the residue corresponding to cysteine 265 is distant from the putative catalytic site, but is in the middle of a region of the enzyme surface proposed to bind tRNA.     

Cysteine 265 Is in the Active Site of, But Is Not Essential for Catalysis by tRNA-Guanine Transglycosylase (TGT) from Escherichia coli

Site-directed mutagenesis and X-ray absorption spectroscopy studies have previously shown that the tRNA-guanine transglycosylase (TGT) from Escherichia coli is a zinc metalloprotein and identified the enzymic ligands to the zinc [Chong et al. (1995), Biochemistry 34, 3694–3701; Garcia et al. (1966), Biochemistry 35, 3133–3139]. During these studies one mutant, TGT (C265A), was found to exhibit a significantly lower specific activity, but was not found to be involved in the zinc site. The present report demonstrates that TGT is inactivated by treatment with thiol reagents (e.g., DTNB, MMTS, and N-ethylmaleimide). Further, this inactivation is shown to be due to modification of cysteine 265. The kinetic parameters for the mutants TGT (C265A) and TGT (C265S), however, suggest that this residue is not performing a critical role in the TGT reaction. We conclude that cysteine 265 is in the active site of TGT, but is not performing a critical catalytic function. This conclusion is supported by the recent determination of the X-ray crystal structure of the TGT from Zymomonas mobilis [Romier et al. (1966), EMBO J. 15, 2850–2857], which reveals that the residue corresponding to cysteine 265 is distant from the putative catalytic site, but is in the middle of a region of the enzyme surface proposed to bind tRNA.     

Experimental validation of molecular dynamics simulations of lipid bilayers: a new approach

A novel protocol has been developed for comparing the structural properties of lipid bilayers determined by simulation with those determined by diffraction experiments, which makes it possible to test critically the ability of molecular dynamics simulations to reproduce experimental data. This model-independent method consists of analyzing data from molecular dynamics bilayer simulations in the same way as experimental data by determining the structure factors of the system and, via Fourier reconstruction, the overall transbilayer scattering-density profiles. Multi-nanosecond molecular dynamics simulations of a dioleoylphosphatidylcholine bilayer at 66% RH (5.4 waters/lipid) were performed in the constant pressure and temperature ensemble using the united-atom GROMACS and the all-atom CHARMM22/27 force fields with the GROMACS and NAMD software packages, respectively. The quality of the simulated bilayer structures was evaluated by comparing simulation with experimental results for bilayer thickness, area/lipid, individual molecular-component distributions, continuous and discrete structure factors, and overall scattering-density profiles. Neither the GROMACS nor the CHARMM22/27 simulations reproduced experimental data within experimental error. The widths of the simulated terminal methyl distributions showed a particularly strong disagreement with the experimentally observed distributions. A comparison of the older CHARMM22 with the newer CHARMM27 force fields shows that significant progress is being made in the development of atomic force fields for describing lipid bilayer systems empirically.     

The RNA-binding PUA domain of archaeal tRNA-guanine transglycosylase is not required for archaeosine formation

Bacterial tRNA-guanine transglycosylase (TGT) replaces the G in position 34 of tRNA with preQ(1), the precursor to the modified nucleoside queuosine. Archaeal TGT, in contrast, substitutes preQ(0) for the G in position 15 of tRNA as the first step in archaeosine formation. The archaeal enzyme is about 60% larger than the bacterial protein; a carboxyl-terminal extension of 230 amino acids contains the PUA domain known to contact the four 3'-terminal nucleotides of tRNA. Here we show that the C-terminal extension of the enzyme is not required for the selection of G15 as the site of base exchange; truncated forms of Pyrococcus furiosus TGT retain their specificity for guanine exchange at position 15. Deletion of the PUA domain causes a 4-fold drop in the observed k(cat) (2.8 x 10(-3) s(-1)) and results in a 75-fold increased K(m) for tRNA(Asp)(1.2 x 10(-5) m) compared with full-length TGT. Mutations in tRNA(Asp) altering or abolishing interactions with the PUA domain can compete with wild-type tRNA(Asp) for binding to full-length and truncated TGT enzymes. Whereas the C-terminal domains do not appear to play a role in selection of the modification site, their relevance for enzyme function and their role in vivo remains to be discovered.     

Development of platelet contractile force as a research and clinical measure of platelet function

This article reviews work performed at the Medical College of Virginia of Virginia Commonwealth University during the development of a whole-blood assay of platelet function. The new assay is capable of assessing platelet function during clotting and thus allows measurement of the contribution of platelets to thrombin generation. Because platelets are monitored in the presence of thrombin, the test gages platelets under conditions of maximal activation. Three parameters are simultaneously assessed on one 700-μL sample of citrated whole blood. Platelet contractile force (PCF), the force produced by platelets during clot retraction, is directly measured as a function of time. This parameter is sensitive to platelet number, platelet metabolic status, glycoprotein IIb/IIIa status, and the presence of antithrombin activities. Clot elastic modulus (CEM), also measured as a function of time, is sensitive to fibrinogen concentration, platelet concentration, the rate of thrombin generation, the flexibility of red cells, and the production of force by platelets. The third parameter, the thrombin generation time (TGT) is determined from the PCF kinetics curve. Because PCF is absolutely thrombin dependent (no thrombin—no force), the initial upswing in PCF occurs at the moment of thrombin production. TGT is sensitive to clotting factor deficiencies, clotting factor inhibitors, and the presence of antithrombins, all of which prolong the TGT and are known to be hemophilic states. Treatment of hemophilic states with hemostatic agents shortens, the TGT toward normal. TGT has been demonstrated to be shorter and PCF to be increased in coronary artery disease, diabetes mellitus, and several other thrombophilic states. Treatment of thrombophilic states with a variety of heparin and nonheparin anticoagulants prolongs the TGT toward normal. The combination of PCF, CEM, and TGT measured on the same sample may allow rapid assessment of global hemostasis and the response to a variety of procoagulant and anticoagulant medications.     

Characterization of cDNA encoding the human tRNA-guanine transglycosylase (TGT) catalytic subunit

Queuosine (Q) is a 7-deazaguanosine found in the first position of the anticodon of tRNAs that recognize NAU and NAC codons (Tyr, Asn, Asp and His). Eukaryotes synthesize Q by the base-for-base exchange of queuine (Q base) for guanine in the unmodified tRNA, a reaction catalyzed by TGT. A search of the human EST database for sequences with significant homology to the well studied TGT from Escherichia coli identified several candidates for full-length (1.3-1.4 kb) cDNA clones. Three candidate cDNA clones, available from IMAGE Consortium, LLNL, (Lennon et al., 1996, Genomics 33, 151-152) were obtained: IMAGE Clone Id Nos. 611146, 1422928, and 72154. Here we report the complete sequences of these clones. IMAGE:72154 contains an ORF encoding a 44 kDa polypeptide with high homology to bacterial TGTs and was subcloned into the mammalian expression vector pMAMneo-Cat. When this construct was transfected into the TGT-negative cell line, GC(3)/c1 (Gündüz et al., 1992, Biochim. Biophys. Acta 1139, 229-238), it restored the ability of the cells to form Q-containing tRNA. This TGT cDNA sequence is encoded in human chromosome 19 clone CTC-539A10 (GenBank accession no. AC011475), enabling determination of the exon-intron boundaries for the TGT gene. The sequence of IMAGE:611146 is 5'-truncated by 76 bp compared to that from IMAGE:72154 and, except for two differences in the 3'-non-coding region, the remainder of the sequence is identical to that of IMAGE:72154. IMAGE:1422928 is a 1390 bp chimera: the 5'-portion, bp 1-708, is identical to a genomic DNA sequence from chromosome 15 (GenBank accession no. AC067805, bp 148976-149683); the 3'-end, bp 726-1390, is identical to the 3'-end of the TGT cDNA sequence from IMAGE:611146.     

Molecular dynamics simulation of a decasaccharide fragment of heparin in aqueous solution

Molecular dynamics (MD) simulations on heparin-water-sodium systems were carried out in order to establish a simulation protocol able to represent heparin solution conformation under physiological conditions. Atomic charges suitable for heparin oligosaccharides were obtained from ab initio quantum-mechanical computations, at the 6-31G(**) level. The GROMACS forcefield, the SPC, and SPC/E water models were employed. Also heparin was simulated with IdoA residues in 1C(4) or 2S(0) conformational states. The results of the performed MD simulations are in agreement with the available experimental data, suggesting that this approach can be applied for the study of heparin interactions with its target proteins and thus play a role in the development of new antithrombotic agents.     

Exploring electron transfer between heme proteins of cytochrome C super family in biosensors: a molecular mechanics study

Electron transfer between heme proteins with mediators plays an important role in the fabrication of sensitive bio-nano sensors. Heme protein Cytochrome c (pdb code - 1HRC) was chosen as the mediator with Cytochrome c' (pdb code - 1A7V) as the probe protein for our investigation on the electron transfer process. We used the software GRAMM, HEX, and MACRODOX to build the protein complex with further evaluation by GROMACS potential. After molecular mechanics refinement by GROMACS the protein complexes were evaluated in terms of the following criteria: Hydrophobic packing, proximity of the hemes, hydrogen bonds, enthalpy and entropy of binding. The free energy was calculated for each complex to derive the feasible stable models. The combined electron transport of the chosen geometric models was evaluated to choose the possible models. Electrostatic potential was calculated using the program APBS around the heme in the presence and absence of other proteins. From our studies, we derived multiple feasible models and possible electronic path. These studies helped us to understand the relay mechanism between the two proteins and to design mutant proteins by rational site directed mutagenesis to enhance the redox potential and thereby improving the signal to noise ratio in amperometric bionano sensors.     

Vinculin transmits high-level integrin tensions that are dispensable for focal adhesion formation

Focal adhesions (FAs) transmit force and mediate mechanotransduction between cells and the matrix. Previous studies revealed that integrin-transmitted force is critical to regulate FA formation. As vinculin is a prominent FA protein implicated in integrin tension transmission, this work studies the relation among integrin tensions (force), vinculin (protein), and FA formation (structure) by integrin tension manipulation, force visualization and vinculin knockout (KO). Two DNA-based integrin tension tools are adopted: tension gauge tether (TGT) and integrative tension sensor (ITS), with TGT restricting integrin tensions under a designed T_tol (tension tolerance) value and ITS visualizing integrin tensions above the T_tol value by fluorescence. Results show that large FAs (area >1 μm²) were formed on the TGT surface with T_tol of 54 pN but not on those with lower T_tol values. Time-series analysis of FA formation shows that focal complexes (area <0.5 μm²) appeared on all TGT surfaces 20 min after cell plating, but only matured to large FAs on TGT with T_tol of 54 pN. Next, we tested FA formation in vinculin KO cells on TGT surfaces. Surprisingly, the T_tol value of TGT required for large FA formation is drastically decreased to 23 pN. To explore the cause, we visualized integrin tensions in both wild-type and vinculin KO cells using ITS. The results showed that integrin tensions in FAs of wild-type cells frequently activate ITS with T_tol of 54 pN. With vinculin KO, however, integrin tensions in FAs became lower and unable to activate 54 pN ITS. Force signal intensities of integrin tensions reported by 33 and 43 pN ITS were also significantly reduced with vinculin KO, suggesting that vinculin is essential to transmit high-level integrin tensions and involved in transmitting intermediate-level integrin tensions in FAs. However, the high-level integrin tensions transmitted by vinculin are not required by FA formation.