Conformation of ring single-stranded DNA measured by DNA origami structures

Measuring the mechanical properties of single-stranded DNA (ssDNA) is a complex challenge that has been addressed lately by different methods. We measured the persistence length of ring ssDNA using a combination of a special DNA origami structure, a self-avoiding ring polymer simulation model, and nonparametric estimation statistics. The method overcomes the complexities set forth by previously used methods. We designed the DNA origami nano structures and measured the ring ssDNA polymer conformations using atomic force microscopy. We then calculated their radius of gyration, which was used as a fitting parameter for finding the persistence length. As there is no simple formulation for the radius of gyration distribution, we developed a simulation program consisting of a self-avoiding ring polymer to fit the persistence length to the experimental data. ssDNA naturally forms stem-loops, which should be taken into account in fitting a model to the experimental measurement. To overcome that hurdle, we found the possible loops using minimal energy considerations and used them in our fitting procedure of the persistence length. Due to the statistical nature of the loops formation, we calculated the persistence length for different percentages of loops that are formed. In the range of 25–75% loop formation, we found the persistence length to be 1.9–4.4 nm, and for 50% loop formation we get a persistence length of 2.83 ± 0.63 nm. This estimation narrows the previously known persistence length and provides tools for finding the conformations of ssDNA.     

Scaling of the equilibrium sedimentation distribution in dense DNA solutions

DNA molecules, several persistence lengths long in sedimentation equilibrium at speeds high enough to maintain fairly close packing, show a dense, sharply-bounded turbid phase and an isotropic phase (as with shorter fragments) and also an intermediate, somewhat turbid region. The concentration distribution in the isotropic phase is in satisfactory agreement with a simple extension of scaled particle theory in which semiflexible chains are equivalent to straight rods of the same length. The net intermolecular interactions, as inferred from the Zimm cluster integral, are purely repulsive. As in our previous study with short fragments, the results are compatible with a hard-core electrostatic radius, decreasing with increasing salt concentration. However, for the longer fragments it is necessary to infer either a slightly greater mass per unit length or a slightly smaller electrostatic radius for closest agreement with scaled particle theory. The properties of the solution at the boundary with the turbid, presumably strongly ordered phase are consistent with those found for shorter fragments and with theoretical scaling expectation for a hard, asymmetric particle.     

Shapes of semiflexible polymers in confined spaces

Biological macromolecules, living in the confines of a cell, often adopt conformations that are unlikely to occur in free space. In this paper, we investigate the effects of confinement on the shape of a semiflexible chain. Results of Monte Carlo simulations show the existence of a shape transition when the persistence length of the polymer becomes comparable to the dimensions of the box. An order parameter is introduced to quantify this behavior. A simple model is constructed to study the effect of the shape transition on the effective persistence length of the polymer.     

Computer simulation of hydrodynamic properties of semiflexible macromolecules: randomly broken chains, wormlike chains, and analysis of properties of DNA

The translational and rotational diffusion coefficients and the intrinsic viscosity of semiflexible, randomly broken, and wormlike chains have been obtained by Monte Carlo simulation in the context of the rigid-body treatment. Both approximate and rigorous rigid-body hydrodynamics are used, so that the error introduced by the approximate methods can be evaluated. A randomly broken chain and a wormlike chain having the same contour length and persistence length have the same radius of gyration but different values for any of the hydrodynamic properties. The two types of chains are compared in this regard. Considering that the cross section of the chain is represented by a cylinder better than by a string of spheres, we devise a cylindrical correction to be applied to the results simulated for chains of beads. Application is made to the analysis of experimental data for the translational and rotational coefficients of DNA fragments with up to 10(3) base pairs, obtaining the persistence length for each model. The values for the wormlike chain agree well with model-independent values obtained from radii of gyration and with other literature data at varying ionic strength. The randomly broken chain is equally able to reproduce the experimental length dependence of the properties, but the resulting persistence length may be too high.     

Transmission of poly(dT60) single-stranded DNA through polycarbonate track-etched ultrafiltration membranes

The objective of this study was to examine membrane filtration of a single stranded DNA (ssDNA) with 60 thymine nucleotides, and to elucidate the variables controlling its transmission across track-etched porous membranes. Dead end filtration measurements were performed using different pore size membranes (10, 15, and 30 nm) at different transmembrane pressures in solutions with ionic strength ranging from 0 to 1000 mM NaCl. The diffusivity of the ssDNA was determined using fluorescence recovery after photobleaching, yielding hydrodynamic radii ranging from 1.6 to 2.8 nm, with values decreasing with increasing solution ionic strength. Despite the small ssDNA/membrane pore size, nearly 100% rejection was observed for measurements performed with the 10 and 15 nm pore size membranes under low-ionic strength conditions. These high rejections can be attributed to strong repulsive electrostatic ssDNA-membrane interactions. With increasing ionic strength, electrostatic interactions as well as the effective size of the ssDNA decreases and the flexibility of the ssDNA increases, leading to a reduction in ssDNA rejection. A design of experiments approach was used to plan filtration experiments that adequately covered the variable space with a manageable number of experiments. The results yielded an empirical expression relating ssDNA rejection to pore size, solution ionic strength and transmembrane pressure. There was evidence of flow induced elongation at high-transmembrane pressures in the 30 nm pore size membranes, but not in the smaller pore size membranes. These results are consistent with critical flux estimates developed using a free draining model for the ssDNA.     

Affinity enhancement of nanobody binding to EGFR: in silico site-directed mutagenesis and molecular dynamics simulation approaches

Epidermal growth factor receptor (EGFR), a transmembrane glycoprotein, is overexpressed in many cancers such as head-neck, breast, prostate, and skin cancers for this reason it is a good target in cancer therapy and diagnosis. In nanobody-based cancer diagnosis and treatment, nanobodies with high affinity toward receptor (e.g. EGFR) results in effective treatment or diagnosis of cancer. In this regard, the main aim of this study is to develop a method based on molecular dynamic (MD) simulations for designing of 7D12 based nanobody with high affinity compared with wild-type nanobody. By surveying electrostatic and desolvation interactions between different residues of 7D12 and EGFR, the critical residues of 7D12 that play the main role in the binding of 7D12 to EGFR were elucidated and based on these residues, five logical variants were designed. Following the 50 ns MD simulations, pull and umbrella sampling simulation were performed for 7D12 and all its variants in complex with EGFR. Binding free energy of 7D12 (and all its variants) with EGFR was obtained by weighted histogram analysis method. According to binding free energy results, GLY101 to GLU mutation showed the highest binding affinity but this variant is unstable after 50 ns MD simulations. ALA100 to GLU mutation shows suitable binding enhancement with acceptable structural stability. Suitable position and orientation of GLU in residue 100 of 7D12 against related amino acids of EGFR formed some extra hydrogen and electrostatic interactions which resulted in binding enhancement.     

Solution studies of the quaternary structure and assembly of human von Willebrand factor

The reversible association of protomers of von Willebrand protein (vWF) was studied in order to analyze the forces and mechanism of vWF polymer assembly. At concentrations of vWF found in plasma (approximately 16 micrograms/mL), disulfide bond reduction with 50 mM 2-mercaptoethanol (2-ME) markedly reduced both vWF activity, as measured by ristocetin-dependent platelet agglutination, and average polymer size (Rh, the mean hydrodynamic radius) in solution, as determined by quasi-elastic light scattering (QLS) and by gel filtration chromatography. With increasing vWF concentration, activity and Rh increased despite reduction of interprotomer disulfide bonds. Changes in temperature after 2-ME treatment produced reversible changes in activity and Rh. Varying the total vWF concentration at any given temperature after 2-ME treatment changed Rh in a consistent and predictable fashion, so that estimates of the dissociation constant for vWF protomer-polymer equilibrium were obtained: Kd5 degrees C = 0.77 micrograms/mL, Kd25 degrees C = 2.4 micrograms/mL, and Kd37 degrees C = 7.7 micrograms/mL, where under the conditions of reduction presented here, the basic protomer of vWF is a dimer. Increasing ionic strength after 2-ME treatment with 1 M KCl did not change Rh, while approximately 100 microM sodium dodecyl sulfate (SDS) or approximately 300 microM sodium deoxycholate (DOC) reduced both Rh and activity compared with those of unreduced polymer. These data show that disulfide bonds are necessary to maintain vWF polymer size and activity at plasma concentrations but that noncovalent forces of association can maintain vWF polymer size and activity at higher concentrations. These forces of association may be important for polymer assembly during intracellular synthesis of vWF.     

Electrostatic contribution to the bending of DNA

A model is derived that accounts for the short-range electrostatic contribution to the bending of DNA molecule in solution and in complexes with proteins in terms of the non-linear Poisson-Boltzmann equation. We defined that the short-range electrostatic interactions depend on the changes of the polyion surface charge density under deformation, while the long-range interactions depend on the bending-induced changes in distances between each two points along the polyion axis. After an appropriate simplification of the Poisson-Boltzmann equation, the short-range term is calculated separately giving the lower limit for the electrostatic contribution to the DNA persistence length. The result is compared with the theoretical approaches developed earlier [M. Fixman, J. Chem. Phys. 76 (1982) 6346; M. Le Bret, J. Chem. Phys. 76 (1982) 6243] and with the experimental data. The conclusion is made that the results of Fixman-Le Bret, which took into account both types of the electrostatic interactions for a uniformly bent polyion, give the upper limit for the electrostatic persistence length at low ionic strength, and the actual behavior of the DNA persistence length lies between two theoretical limits. Only the short-range term is significant at moderate-to-high ionic strength where our results coincide with the predictions of Fixman-Le Bret. The bending of DNA on the protein surface that is accompanied by an asymmetric neutralization of the DNA charge is also analyzed. In this case, the electrostatic bending energy gives a significant favorite contribution to the total bending energy of DNA. Important implications to the mechanisms of DNA-protein interactions, particularly in the nucleosome particle, are discussed.     

Effect of pH on collagen flexibility determined from dilute solution viscoelastic measurements

The frequency dependences of the storage and loss shear moduli, G′ and G″ of dilute solutions of collagen at pH 7.4, ionic strength approximately 0.2, were measured at 15.0°C by the Birnboim—Schrag multiple-lumped resonator apparatus. By use of two solvents, water and 50% glycerol, the effective reduced frequency range was extended to cover 2.5 logarithmic decades. The intrinsic viscosity and longest (rotational) relaxation time were considerably smaller than those determined at pH 4.0 in an earlier study. At pH 4.0, the behaviour could be modelled by a rodlike molecule with partial flexibility along its entire length and a persistence length of 161 nm with no loose joints. The behaviour at pH 7.4 corresponds approximately to the expectation for a semiflexible rod with two loose joints near the ends and a similar persistence length (169nm) for the centre segment.     

A Mesoscale Model of DNA and Its Renaturation

A mesoscale model of DNA is presented (3SPN.1), extending the scheme previously developed by our group. Each nucleotide is mapped onto three interaction sites. Solvent is accounted for implicitly through a medium-effective dielectric constant and electrostatic interactions are treated at the level of Debye-Hückel theory. The force field includes a weak, solvent-induced attraction, which helps mediate the renaturation of DNA. Model parameterization is accomplished through replica exchange molecular dynamics simulations of short oligonucleotide sequences over a range of composition and chain length. The model describes the melting temperature of DNA as a function of composition as well as ionic strength, and is consistent with heat capacity profiles from experiments. The dependence of persistence length on ionic strength is also captured by the force field. The proposed model is used to examine the renaturation of DNA. It is found that a typical renaturation event occurs through a nucleation step, whereby an interplay between repulsive electrostatic interactions and colloidal-like attractions allows the system to undergo a series of rearrangements before complete molecular reassociation occurs.     

Flow birefringence of T7 phage DNA: Dependence on salt concentration

We have determined extinction angles and flow birefringence of T7 bacteriophage DNA over a wide range of shear, polymer concentration, and solvent ionic strength. From these data, information on the simple salt dependence of coil permeability to solvent and on short-range intrachain interactions (persistence length) was obtained. At all ionic strengths, our results are consistent with a partially draining coil in the Gaussian subchain dynamical theory of Rouse-Zimm-Tschoegl-Bloomfield. Salt dependence of persistence length is comparable to, although somewhat less than, that obtained previously using similar methods with a fivefold higher-molecular-weight DNA (T2 bacteriophage DNA). Possible reasons for observed discrepancies are analyzed, and the results of this work are compared in detail to other current studies of solvent ionic strength dependence in persistence length and hydrodynamic properties of DNA.     

Pulling-speed-dependent force-extension profiles for semiflexible chains

We present theory and simulations to describe nonequilibrium stretching of semiflexible chains that serve as models of DNA molecules. Using a self-consistent dynamical variational approach, we calculate the force-extension curves for worm-like chains as a function of the pulling speed, v(0). Due to nonequilibrium effects the stretching force, which increases with v(0), shows nonmonotonic variations as the persistence length increases. To complement the theoretical calculations we also present Langevin simulation results for extensible worm-like chain models for the dynamics of stretching. The theoretical force-extension predictions compare well with the simulation results. The simulations show that, at high enough pulling speeds, the propagation of tension along the chain conformations transverse to the applied force occurs by the Brochard-Wyart's stem-flower mechanism. The predicted nonequilibrium effects can only be observed in double-stranded DNA at large ( approximately 100 microm/s) pulling speeds.     

Molecular dynamics simulations of phenylimidazole inhibitor complexes of cytochrome P450 cam

Molecular dynamics simulations have been performed on three phenylimidazole inhibitor complexes ofP450 _cam, utilizing the X-ray structures and the AMBER suite of programs. Compared to their corresponding optimized X-ray structures, very similar features were observed for the 1-phenylimidazole (1-PI) and 2-phenylimidazole (2-PI) complexes during a 100 ps MD simulation. The 1-PI inhibitor binds as a Type II complex with the imidazole nitrogen as a ligand of the heme iron. Analysis of the inhibitor-enzyme interctions during the MD simulations reveals that electrostatic interactions of the imidazole with the heme and van der Waals interactions of the phenyl ring with nearby hydrophobic residues are dominant. By contrast, 2-PI binds as a Type I inhibitor in the substrate binding pocket, but not as a ligand of the iron. The interactions of this inhibitor are qualitatively different from that of the Type II 1-PI, being mainly electrostatic/H-bonding interactions with a bound water and polar residues. Although the third compound, 4-PI, in common with 1-PI, also binds as a Type II inhibitor, with one nitrogen of the imidazole as a ligand to the iron, the MD average binding orientation deviates significantly from the X-ray structure. The most important changes observed include: (1) the rotation of the imidazole ring of this inhibitor by about 90° to enhance electrostatic interactions of the imidazole NH group with the carbonyl group of LEU244, and (2) the rotation of the carbonyl group of ASP251 to form a H-bond with VAL254. An analysis of the H-bonding network surrounding this substrate in the optimized crystal structure revealed that there is no H-bonding partner either for the free polar NH group in the imidazole ring of 4-phenylimidazole or for the polar carbonyl group of the nearby ASP251 residue. The deviation of the dynamically averaged inhibitor-enzyme structure of the 4-PI complex from the optimized crystal structure can therefore be rationalized as a consequence of the optimization of the electrostatic interactions among the polar groups.     

Characterization of suramin binding sites on the human group IIA secreted phospholipase A2 by site-directed mutagenesis and molecular dynamics simulation

Suramin is a polysulphonated naphthylurea with inhibitory activity against the human secreted group IIA phospholipase A(2) (hsPLA2GIIA), and we have investigated suramin binding to recombinant hsPLA2GIIA using site-directed mutagenesis and molecular dynamics (MD) simulations. The changes in suramin binding affinity of 13 cationic residue mutants of the hsPLA2GIIA was strongly correlated with alterations in the inhibition of membrane damaging activity of the protein. Suramin binding to hsPLA2GIIA was also studied by MD simulations, which demonstrated that altered intermolecular potential energy of the suramin/mutant complexes was a reliable indicator of affinity change. Although residues in the C-terminal region play a major role in the stabilization of the hsPLA2GIIA/suramin complex, attractive and repulsive hydrophobic and electrostatic interactions with residues throughout the protein together with the adoption of a bent suramin conformation, all contribute to the stability of the complex. Analysis of the hsPLA2GIIA/suramin interactions allows the prediction of the properties of suramin analogues with improved binding and higher affinities which may be candidates for novel phospholipase A(2) inhibitors.     

The stability and dynamics of the human calcitonin amyloid peptide DFNKF

The stability and dynamics of the human calcitonin-derived peptide DFNKF (hCT(15-19)) are studied using molecular dynamics (MD) simulations. Experimentally, this peptide is highly amyloidogenic and forms fibrils similar to the full length calcitonin. Previous comparative MD studies have found that the parallel beta-stranded sheet is a stable organization of the DFNKF protofibril. Here, we probe the stability and dynamics of the small parallel DFNKF oligomers. The results show that even small DFNKF oligomers, such as trimers and tetramers, are stable for a sufficient time in the MD simulations, indicating that the crucial nucleus seed size for amyloid formation can be quite small. The simulations also show that the stability of DFNKF oligomers increases with their sizes. The small but stable seed may reflect the experimental rapid formation of the DFNKF fibrils. Further, a noncooperative process of parallel beta-sheet formation from the out-of-register trimer is observed in the simulations. In general, the residues of DFNKF peptides near the N-/C-termini are more flexible, whereas the interior residues are more stable. Simulations of mutants and capped peptides show that both interstrand hydrophobic and electrostatic interactions play important roles in stabilizing the DFNKF parallel oligomers. This study provides insights into amyloid formation.     

Comparison of hard-cylinder and screened coulomb interactions in the modeling of supercoiled DNAs

A 1000 base pair (bp) model supercoiled DNA is simulated using spherical screened Coulomb interactions between subunits on one hand and equivalent hard-cylinder interactions on the other. The amplitudes, or effective charges, of the spherical screened Coulomb electrostatic potentials are chosen so that the electrostatic potential surrounding the middle of a linear array of 2001 subunits (31.8 Å diameter) closely matches the solution of the nonlinear Poisson-Boltzmann equation for a cylinder with 12 Å radius and the full linear charge density of DNA at all distances beyond the 24 Å hard-core diameter. This superposition of spherical screened Coulomb potentials is practically identical to the particular solution of the cylindrical linearized Poisson-Boltzmann equation that matches the solution of the nonlinear Poisson-Boltzmann equation at large distances. The interaction energy between subunits is reckoned from the effective charges according to the standard DLVO expression. The equivalent hard-cylinder diameter is chosen following Stigter's protocol for matching second virial coefficients, but for the full linear charge density of DNA. The electrostatic persistence length of the model with screened Coulomb interactions is extremely sensitive to the (arbitrarily) chosen subunit length at the higher salt concentrations. The persistence length of the hard-cylinder model is adjusted to match that of the screened Coulomb model for each ionic condition. Simulations for a superhelix density σ = -0.05 using a spherical screened Coulomb interaction plus a 24 Å hard-cylinder core (SCPHC) potential indicate that the radius of gyration of this 1000 bp DNA actually undergoes a slight increase as the NaCl concentration is raised from 0.01 to 1.0M. Thus, merely softening the potential from hard-cylinder to screened Coulomb form does not produce a large decrease in radius of gyration with increasing NaCl concentration for DNAs of this size. Radii of gyration, static structure factors, and diffusion coefficients obtained using the equivalent hard-cylinder (EHC) potential agree well with those obtained using the SCPHC potential in 1.0M NaCl, but in 0.1M NaCl the agreement is not as good, and in 0.01M NaCl the agreement is definitely unsatisfactory. These conclusions differ in significant respects from those obtained in previous studies. © 1997 John Wiley & Sons, Inc. Biopoly 42: 455–470, 1997     

Biomechanics of actin filaments: a computational multi-level study

The actin microfilament (F-actin) is a structural and functional component of the cell cytoskeleton. Notwithstanding the primary role it plays for the mechanics of the cell, the mechanical behaviour of F-actin is still not totally explored. In particular, the relationship between the mechanics of F-actin and its molecular architecture is not completely understood. In this study, the mechanical properties of F-actin were related to the molecular topology of its building monomers (G-actin) by employing a computational multi-level approach. F-actins with lengths up to 500 nm were modelled and characterized, using a combination of equilibrium molecular dynamics (MD) simulations and normal mode analysis (NMA). MD simulations were performed to analyze the molecular rearrangements of G-actin in physiological conditions; NMA was applied to compute the macroscopic properties of F-actin from its vibrational modes of motion. Results from this multi-level approach showed that bending stiffness, bending modulus and persistence length are independent from the length of F-actin. On the contrary, the orientations and motions of selected groups of residues of G-actin play a primary role in determining the filament flexibility. In conclusion, this study (i) demonstrated that a combined computational approach of MD and NMA allows to investigate the biomechanics of F-actin taking into account the molecular topology of the filament (i.e., the molecular conformations of G-actin) and (ii) that this can be done using only crystallographic G-actin, without the need of introducing experimental parameters nor of reducing the number of residues.     

Escherichia coli replicative helicase PriA protein-single-stranded DNA complex. Stoichiometries, free energy of binding, and cooperativities

Analyses of interactions of the Escherichia coli replicative helicase, PriA protein, with a single-stranded (ss) DNA have been performed, using the quantitative fluorescence titration technique. The stoichiometry of the PriA helicase.ssDNA complex has been examined in binding experiments with a series of ssDNA oligomers. The total site-size of the PriA.ssDNA complex, i.e. the maximum number of nucleotide residues occluded by the PriA helicase in the complex, is 20 +/- 3 residues per protein monomer. However, the protein can efficiently form a complex with a minimum of 8 nucleotides. Thus, the enzyme has a strong ssDNA-binding site that engages in direct interactions with a significantly smaller number of nucleotides than the total site-size. The ssDNA-binding site is located in the center of the enzyme molecule, with the protein matrix protruding over a distance of approximately 6 nucleotides on both sides of the binding site. The analysis of the binding of two PriA molecules to long oligomers was performed using statistical thermodynamic models that take into account the overlap of potential binding sites, cooperative interactions, and the protein.ssDNA complexes with different stoichiometries. The intrinsic affinity depends little upon the length of the ssDNA. Moreover, the binding is accompanied by weak cooperative interactions.