Dynamics of RNA with A6 bulge loop

In order to understand the behavior of RNAs with large bulges In solution, molecular dynamics was performed on the RNA molecule in water with A6 bulge. The result of simulation showed that nonstacked conformation Is the main conformation in large bulges, and the backbone of large bulge is of great conformational flexibility, but bulges-induced bends are relatively rigid. The fluctuation in bulge has little influence on the bend angle of RNAs.

Structural and functional characterization of the human CCR5 receptor in complex with HIV gp120 envelope glycoprotein and CD4 receptor by molecular modeling studies

The entry of human immunodeficiency virus (HIV) into cells depends on a sequential interaction of the gp120 envelope glycoprotein with the cellular receptors CD4 and members of the chemokine receptor family. The CC chemokine receptor CCR5 is such a receptor for several chemokines and a major coreceptor for the entry of R5 HIV type-1 (HIV-1) into cells. Although many studies focus on the interaction of CCR5 with HIV-1, the corresponding interaction sites in CCR5 and gp120 have not been matched. Here we used an approach combining protein structure modeling, docking and molecular dynamics simulation to build a series of structural models of the CCR5 in complexes with gp120 and CD4. Interactions such as hydrogen bonds, salt bridges and van der Waals contacts between CCR5 and gp120 were investigated. Three snapshots of CCR5-gp120-CD4 models revealed that the initial interactions of CCR5 with gp120 are involved in the negatively charged N-terminus (Nt) region of CCR5 and positively charged bridging sheet region of gp120. Further interactions occurred between extracellular loop2 (ECL2) of CCR5 and the base of V3 loop regions of gp120. These interactions may induce the conformational changes in gp120 and lead to the final entry of HIV into the cell. These results not only strongly support the two-step gp120-CCR5 binding mechanism, but also rationalize extensive biological data about the role of CCR5 in HIV-1 gp120 binding and entry, and may guide efforts to design novel inhibitors.

Fusion dynamics of symmetric systems near barrier energies

Within the framework of the improved isospin-dependent quantum molecular dynamics (ImIQMD) model,he fusion dynamics of symmetric reaction systems are investigated systematically. Calculations show that the number of nucleon transfer in the neck region is appreciably dependent on the incident energies, but strongly on he reaction systems. A comparison of the neck dynamics is performed for the symmetric reactions 58Ni+58Niand 64Ni+64Ni at energies in the vicinity of the Coulomb barrier. An increase of the ratios of the neutron to proton in the neck region at initial collision stage is observed and obvious for the latter system, which reduces the fusion barrier of two colliding nuclei. The distribution of the dynamical fusion barriers and the fusion excitation functions are calculated and compared with the available experimental data.

Nuclear reaction dynamics induced by halo-nuclei at intermediate energy heavy ion collisions

We studied systematically the reaction dynamics induced by neutron-halo nuclei and proton-halo nuclei within the isospin dependent quantum molecular dynamics, such as the effects of loose bound halo-nuclei on the fragmentation reaction and momentum dissipation for different colliding systems with different beam energies and different impact parameters. In order to emphasize the roles of neutron-halo nucleus B-19 and proton-halo nucleus Al-23 on the reaction dynamics we also calculated the the reaction dynamics induced by the stable nuclei F-19 and Na-23 with equal mass under identical incident channel conditions. Based on the comparison of results of reaction dynamics induced by halo-nucleus colliding systems and stable nucleus collidinmg systems we found that the roles of loose bound halo-nucleus structure on the fragmentation multiplicity and nuclear stopping (momentum dissipation) are important for all of colliding systems with different beam energies and minor impact parameters, such as, the loose bound halo-nuclei structure increases the fragmentation multiplicity, but reduces the nuclear stopping.

Dynamics of strangeness production in heavy-ion collisions near threshold energies

Within the framework of the improved isospin-dependent quantum molecular dynamics (ImIQMD) model, the dynamics of strangeness (K-0,K-+, Lambda, and Sigma(-,0,+)) production in heavy-ion collisions near threshold energies is investigated systematically, with the strange particles considered to be produced mainly by inelastic collisions of baryon-baryon and pion-baryon. Collisions in the region of suprasaturation densities of the dense baryonic matter formed in heavy-ion collisions dominate the yields of strangeness production. Total multiplicities as functions of incident energies and collision centralities are calculated with the Skyrme parameter SLy6. The excitation function of strangeness production is analyzed and also compared with the KaoS data for K+ production in the reactions C-12 + C-12 and Au-197 + Au-197.

Dynamics of pion production in heavy-ion collisions around 1A GeV energies

Within the framework of the improved isospin-dependent quantum molecular dynamics model, the dynamics of pion emission in heavy-ion collisions in the region of 1A GeV energies as a probe of nuclear symmetry energy at suprasaturation densities is investigated systematically. The total pion multiplicities and the pi(-)/pi(+) yields are calculated for selected Skyrme parameters SkP, SLy6, Ska, and SIII and also for the cases of different stiffness of symmetry energy with the parameter SLy6. The influence of Coulomb potential, symmetry energy, and in-medium pion potential on the pion production is investigated and compared to each other by analyzing the distributions of transverse momentum and longitudinal rapidity and also the excitation functions of the total pion and the pi(-)/pi(+) ratio. The directed flow, elliptic flow, and polar-angle distributions are calculated for the cases of different collision centralities and also the various stiffnesses of the symmetry energies. A comparison of the calculations with the available experimental data is performed.

Influence of molecular topology on the static and dynamic properties of single polymer chain in solution

The influence of molecular topology on the structural and dynamic properties of polymer chain in solution with ring structure, three-arm branched structure, and linear structure are studied by molecular dynamics simulation. At the same degree of polymerization (N), the ring-shaped chain possesses the smallest size and largest diffusion coefficient. With increasing N, the difference of the radii of gyration between the three types of polymer chains increases, whereas the difference of the diffusion coefficients among them decreases. However, the influence of the molecular topology on the static and the dynamic scaling exponents is small. The static scaling exponents decrease slightly, and the dynamic scaling exponents increase slightly, when the topology of the polymer chain is changed from linear to ring-shaped or three-arm branched architecture. The dynamics of these three types of polymer chain in solution is Zimm-like according to the dynamic scaling exponents and the dynamic structure factors.

Exploring the mechanism of flexible biomolecular recognition with single molecule dynamics

Combining a single-molecule study of protein binding with a coarse grained molecular dynamics model including solvent (water molecules) effects, we find that biomolecular recognition is determined by flexibilities in addition to structures. Our single-molecule study shows that binding of CBD (a fragment of Wiskott-Aldrich syndrome protein) to Cdc42 involves bound and loosely bound states, which can be quantitatively explained in our model as a result of binding with large conformational changes. Our model identified certain key residues for binding consistent with mutational experiments. Our study reveals the role of flexibility and a new scenario of dimeric binding between the monomers: first bind and then fold.

Energetic basis of molecular recognition in a DNA aptamer

The thermal stability and ligand binding properties of the L-argininamide-binding DNA aptamer (5'-GATCGAAACGTAGCGCCTTCGATC3') were studied by spectroscopic and calorimetric methods. Differential calorimetric studies showed that the uncomplexed aptamer melted in a two-state reaction with a melting temperature T-m = 50.2 +/- 0.2 degrees C and a folding enthalpy Delta H degrees(fold) = -49.0 +/- 2.1 kcal mol(-1). These values agree with values of T-m = 49.6 degrees C and Delta H degrees(fold) = -51.2 kcal mol(-1) predicted for a simple hairpin structure. Melting of the uncomplexed aptamer was dependent upon salt concentration, but independent of strand concentration. The T of aptamer melting was found to increase as L-argininamide concentrations increased. Analysis of circular dichroism titration data using a single-site binding model resulted in the determination of a binding free energy Delta G degrees(bind) = -5.1 kcal mol(-1). Isothermal titration calorimetry studies revealed an exothermic binding reaction with Delta H degrees(bind) = -8.7 kcal mol(-1). Combination of enthalpy and free energy produce ail unfavorable entropy of -T Delta S degrees = +3.6 kcal mol(-1). A molar heat capacity change of -116 cal mol(-1) K-1 was determined from calorimetric measurements at four temperatures over the range of 15-40 degrees C. Molecular dynamics simulations were used to explore the structures of the unligated and ligated aptamer structures.

Single-molecule dynamics reveals cooperative binding-folding in protein recognition

The study of associations between two biomolecules is the key to understanding molecular function and recognition. Molecular function is often thought to be determined by underlying structures. Here, combining a single-molecule study of protein binding with an energy-landscape-inspired microscopic model, we found strong evidence that biomolecular recognition is determined by flexibilities in addition to structures. Our model is based on coarse-grained molecular dynamics on the residue level with the energy function biased toward the native binding structure ( the Go model). With our model, the underlying free-energy landscape of the binding can be explored. There are two distinct conformational states at the free-energy minimum, one with partial folding of CBD itself and significant interface binding of CBD to Cdc42, and the other with native folding of CBD itself and native interface binding of CBD to Cdc42. This shows that the binding process proceeds with a significant interface binding of CBD with Cdc42 first, without a complete folding of CBD itself, and that binding and folding are then coupled to reach the native binding state.

Dynamics of Delta/Notch signaling on endomesoderm segregation in the sea urchin embryo.

Endomesoderm is the common progenitor of endoderm and mesoderm early in the development of many animals. In the sea urchin embryo, the Delta/Notch pathway is necessary for the diversification of this tissue, as are two early transcription factors, Gcm and FoxA, which are expressed in mesoderm and endoderm, respectively. Here, we provide a detailed lineage analysis of the cleavages leading to endomesoderm segregation, and examine the expression patterns and the regulatory relationships of three known regulators of this cell fate dichotomy in the context of the lineages. We observed that endomesoderm segregation first occurs at hatched blastula stage. Prior to this stage, Gcm and FoxA are co-expressed in the same cells, whereas at hatching these genes are detected in two distinct cell populations. Gcm remains expressed in the most vegetal endomesoderm descendant cells, while FoxA is downregulated in those cells and activated in the above neighboring cells. Initially, Delta is expressed exclusively in the micromeres, where it is necessary for the most vegetal endomesoderm cell descendants to express Gcm and become mesoderm. Our experiments show a requirement for a continuous Delta input for more than two cleavages (or about 2.5 hours) before Gcm expression continues in those cells independently of further Delta input. Thus, this study provides new insights into the timing mechanisms and the molecular dynamics of endomesoderm segregation during sea urchin embryogenesis and into the mode of action of the Delta/Notch pathway in mediating mesoderm fate.

On the structure and dynamics of ionic liquids

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.