Biomedical subjects

Yuguang Mu

Publications and source records attributed to Yuguang Mu.

8 recordsLinked to original sources

PubMed2006

Charge structure and counterion distribution in hexagonal DNA liquid crystal.

A hexagonal liquid crystal of DNA fragments (double-stranded, 150 basepairs) with tetramethylammonium (TMA) counterions was investigated with small angle neutron scattering (SANS). We obtained the structure factors pertaining to the DNA and counterion density correlations with contrast matching in the water. Molecular dynamics (MD) computer simulation of a hexagonal assembly of nine DNA molecules showed that the inter-DNA distance fluctuates with a correlation time around 2 ns and a standard deviation of 8.5% of the interaxial spacing. The MD simulation also showed a minimal effect of the fluctuations in inter-DNA distance on the radial counterion density profile and significant penetration of the grooves by TMA. The radial density profile of the counterions was also obtained from a Monte Carlo (MC) computer simulation of a hexagonal array of charged rods with fixed interaxial spacing. Strong ordering of the counterions between the DNA molecules and the absence of charge fluctuations at longer wavelengths was shown by the SANS number and charge structure factors. The DNA-counterion and counterion structure factors are interpreted with the correlation functions derived from the Poisson-Boltzmann equation, MD, and MC simulation. Best agreement is observed between the experimental structure factors and the prediction based on the Poisson-Boltzmann equation and/or MC simulation. The SANS results show that TMA is too large to penetrate the grooves to a significant extent, in contrast to what is shown by MD simulation.

Computer Simulation↗

PubMed2006

Folding, misfolding, and amyloid protofibril formation of WW domain FBP28.

We study the folding mechanism of a triple beta-strand WW domain from the Formin binding protein 28 (FBP28) at atomic resolution with explicit water model using replica exchange molecular dynamics computer simulations. Extended sampling over a wide range of temperatures to obtain the free energy, enthalpy, and entropy surfaces as a function of structural coordinates has been performed. Simulations were started from different configurations covering the folded and unfolded states. In the free energy landscape a transition state is identified and its structures and -values are compared with experimental data from a homologous protein, the prolyl-isomerase Pin1 WW domain. A stable intermediate state is found to accumulate during the simulation characterized by the carboxyl-terminal beta-strand 3 having misregistered hydrogen bonds and where the structural heterogeneity is due to nonnative turn II formation. Furthermore, the aggregation behavior of the FBP28 WW domain may be related to one such misfolded structure, which has a much lower free energy of dimer formation than that of the native dimer. Based on the misfolded dimer, aggregation to form protofibril structure is discussed.

Amyloid↗

PubMed2005

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study.

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.

Binding Sites↗

PubMed2005

Structure and energy landscape of a photoswitchable peptide: a replica exchange molecular dynamics study.

A replica exchange molecular dynamics (REMD) simulation of a bicyclic azobenzene peptide in explicit dimethyl sulfoxide solution is presented in order to characterize the conformational structures and energy landscape of a photoswitchable peptide. It is shown that an enhanced-sampling technique such as the REMD method is essential to obtain a converged conformational sampling of the peptide at room temperature. This is because conventional MD simulations of less than approximately 100-ns length are either trapped in local minima (at 295 K) or-if run at high temperature-do not resemble the room-temperature REMD results. Calculating various nuclear Overhauser effects (NOEs) and (3)J-couplings, a good overall agreement between the REMD simulations and the NMR experiments of Renner et al. (Biopolymers 2000;54:501-514) is found. In particular, the REMD study confirms the general picture drawn by Renner et al. that the trans-isomer of the azobenzene peptide exhibits a well-defined structure, while the cis-isomer is a conformational heterogeneous system; that is, the trans-isomer occurs in 2 well-defined conformers, while the cis-isomer represents an energetically frustrated system that leads to an ensemble of conformational structures. Employing a principal component analysis of the REMD data, the free energy landscape of the systems is studied at various temperatures. The implications for the folding and unfolding pathways of the system are discussed.

Algorithms↗

PubMed2005

Energy landscape of a small peptide revealed by dihedral angle principal component analysis.

A 100 ns molecular dynamics simulation of penta-alanine in explicit water is performed to study the reversible folding and unfolding of the peptide. Employing a standard principal component analysis (PCA) using Cartesian coordinates, the resulting free-energy landscape is found to have a single minimum, thus suggesting a simple, relatively smooth free-energy landscape. Introducing a novel PCA based on a transformation of the peptide dihedral angles, it is found, however, that there are numerous free energy minima of comparable energy (less than or approximately 1 kcal/mol), which correspond to well-defined structures with characteristic hydrogen-bonding patterns. That is, the true free-energy landscape is actually quite rugged and its smooth appearance in the Cartesian PCA represents an artifact of the mixing of internal and overall motion. Well-separated minima corresponding to specific conformational structures are also found in the unfolded part of the free energy landscape, revealing that the unfolded state of penta-alanine is structured rather than random. Performing a connectivity analysis, it is shown that neighboring states are connected by low barriers of similar height and that each state typically makes transitions to three or four neighbor states. Several principal pathways for helix nucleation are identified and discussed in some detail.

Energy Transfer↗

PubMed2005

Structure and dynamics of an RNA tetraloop: a joint molecular dynamics and NMR study.

Molecular dynamics simulations of the RNA tetraloop 5'-CGCUUUUGCG-3' with high melting temperature and significant conformational heterogeneity in explicit water solvent are presented and compared to NMR studies. The NMR data allow for a detailed test of the theoretical model, including the quality of the force field and the conformational sampling. Due to the conformational heterogeneity of the tetraloop, high temperature (350 K) and locally enhanced sampling simulations need to be invoked. The Amber98 force field leads to a good overall agreement with experimental data. Based on NMR data and a principal component analysis of the 350 K trajectory, the dynamic structure of the tetraloop is revealed. The principal component free energy surface exhibits four minima, which correspond to well-defined conformational structures that differ mainly by their base stacking in the loop region. No correlation between the motion of the sugar rings and the stacking dynamics of the loop bases is found.

Nuclear Magnetic Resonance, Biomolecular↗

PubMed2005

Synthesis of spin-labeled RNAs for long range distance measurements by peldor.

Site directed spin labeled RNA duplexes with different interspin distances were synthesized. The radical 2,2,5,5-tetramethyl-pyrrolin- 1-yloxyl-3-acetylene (TPA) was introduced during the solid-phase synthesis through a Sonogashira cross-coupling with 5-iodo-uridine. Tm and CD studies showed that the spin label does not to disturb significantly the A-form of these duplexes. 4-Pulse Electron Double Resonance (PELDOR) was then used to measure intramolecular spin-spin distances of 19.3, 33.0 and 40.9 A, which are in very good agreement with the calculated values of 17.6, 32.1 and 39.1 A, obtained from Molecular Dynamics (MD) simulations.

Biochemistry↗

PubMed2004

A PELDOR-based nanometer distance ruler for oligonucleotides.

A pulsed electron paramagnetic resonance (EPR) spectroscopic ruler for oligonucleotides was developed using a series of duplex DNAs. The spin-labeling is accomplished during solid-phase synthesis of the oligonucleotides utilizing a palladium-catalyzed cross-coupling reaction between 5-iodo-2'-deoxyuridine and the rigid spin-label 2,2,5,5-tetramethyl-pyrrolin-1-yloxyl-3-acetylene (TPA). 4-Pulse electron double resonance (PELDOR) was then used to measure the intramolecular spin-spin distances via the dipolar coupling, yielding spin-spin distances of 19.2, 23.3, 34.7, 44.8, and 52.5 A. Employing a full-atom force field with explicit water, molecular dynamic (MD) simulations on the same spin-labeled oligonucleotides in their duplex B-form gave spin-spin distances of 19.6, 21.4, 33.0, 43.3, and 52.5 A, respectively, in very good agreement with the measured distances. This shows that the oligonucleotides adopt a B-form duplex structure also in frozen aqueous buffer solution. It also demonstrates that the combined use of site-directed spin-labeling, PELDOR experiments, and MD simulations can yield a microscopic picture about the overall structure of oligonucleotides. The technique is also applicable to more complex systems, like ribozymes or DNA/RNA-protein complexes, which are difficult to access by NMR or X-ray crystallography.

Base Sequence↗