Biomedical subjects

Ao Ma

Publications and source records attributed to Ao Ma.

9 recordsLinked to original sources

PubMed2025

Effects of Lamina-Chromatin Attachment on Super Long-Range Chromatin Interactions.

The interactions between chromatin and lamin proteins localized on the nuclear envelope play a crucial role in the three-dimensional (3D) organization of the genome. This study investigates the influence of lamin associated domains (LADs) on genome organization at the chromosome level using 3D polymer models of mouse embryonic fibroblasts (MEFs) and embryonic stem cells (mESCs). By integrating genome-wide LAD maps from DamID assays, we simulated chromatin conformations with and without LAD attachment to the nuclear envelope. Our results show that incorporating LAD-lamin interactions yields a radial chromatin distribution consistent with experimental observations. Moreover, LAD-lamin interactions induce significant super long-range chromatin contacts across distant genomic regions. These findings suggest two distinct mechanisms driving induction of chromatin interactions by LAD-lamin attachment.

3D single cell conformations↗

PubMed2006

Implications of alternative substrate binding modes for catalysis by uracil-DNA glycosylase: an apparent discrepancy resolved.

A theoretical study showed that the base excision repair enzyme uracil-DNA glycosylase (UDG) exploits electrostatic interactions with backbone phosphate groups in the substrate for catalysis. Although experiments performed to test the calculated results confirmed the predicted importance of the -2, -1, and +1 phosphate groups, there was an apparent disagreement with regard to the +2 phosphate group. The calculations indicated that it made an important contribution, while experimentally, the effect of its deletion or neutralization was small. The +2 phosphate group interacts directly with an active site histidine (H148 in humans) in the crystal structure of UDG in complex with double-stranded (ds) DNA. We previously calculated that H148 has a strong anticatalytic effect due to its protonation, and here we use alchemical free energy simulations to estimate its site-specific pKa. We find that it is positively charged over the entire experimental pH range (4-10), so its deprotonation cannot compensate for deletion or neutralization of the +2 phosphate group. The free energy simulations are facilitated by an efficient charge-scaling procedure that allows quantitative correction for the implicit treatment of solvent far from the active site; improvements are made to that method to account carefully for differences in the truncation of electrostatic interactions in the contributing molecular-mechanical and continuum-electrostatic approaches. Additional simulations are used to demonstrate that the +2 phosphate group is fully solvent exposed in complexes with single-stranded DNA substrates like those used in the experiments. In contrast, it is well-structured and buried in the dsDNA complex used in the original simulations. Differences in solvent shielding thus account for the apparent lack of an effect observed experimentally upon neutralization or deletion of this group.

Catalysis↗

PubMed2006

Bias annealing: a method for obtaining transition paths de novo.

Computational studies of dynamics in complex systems require means for generating reactive trajectories with minimum knowledge about the processes of interest. Here, we introduce a method for generating transition paths when an existing one is not already available. Starting from biased paths obtained from steered molecular dynamics, we use a Monte Carlo procedure in the space of whole trajectories to shift gradually to sampling an ensemble of unbiased paths. Application to basin-to-basin hopping in a two-dimensional model system and nucleotide-flipping by a DNA repair protein demonstrates that the method can efficiently yield unbiased reactive trajectories even when the initial steered dynamics differ significantly. The relation of the method to others and the physical basis for its success are discussed.

Journal Article↗

PubMed2006

Dynamic coupling between coordinates in a model for biomolecular isomerization.

To understand a complex reaction, it is necessary to project the dynamics of the system onto a low-dimensional subspace of physically meaningful coordinates. We recently introduced an automatic method for identifying coordinates that relate closely to stable-state commitment probabilities and successfully applied it to a model for biomolecular isomerization, the C(7eq)-->alpha(R) transition of the alanine dipeptide [A. Ma and A. R. Dinner, J. Phys. Chem. B 109, 6769 (2005)]. Here, we explore approximate means for estimating diffusion tensors for systems subject to restraints in one and two dimensions and then use the results together with an extension of Kramers theory for unimolecular reaction rates [A. Berezhkovskii and A. Szabo, J. Chem. Phys. 122, 014503 (2005)] to show explicitly that both the potential of mean force and the diffusion tensor are essential for describing the dynamics of the alanine dipeptide quantitatively. In particular, the signficance of off-diagonal elements of the diffusion tensor suggests that the coordinates of interest are coupled by the hydrodynamic-like response of the bath of remaining degrees of freedom.

Acetamides↗

PubMed2006

Monte Carlo simulations of biomolecules: The MC module in CHARMM.

We describe the implementation of a general and flexible Monte Carlo (MC) module for the program CHARMM, which is used widely for modeling biomolecular systems with empirical energy functions. Construction and use of an almost arbitrary move set with only a few commands is made possible by providing several predefined types of moves that can be combined. Sampling can be enhanced by noncanonical acceptance criteria, automatic optimization of step sizes, and energy minimization. A systematic procedure for improving MC move sets is introduced and applied to simulations of two peptides. The resulting move sets allow MC to sample the configuration spaces of these systems much more rapidly than Langevin dynamics. The rate of convergence of the difference in free energy between ethane and methanol in explicit solvent is also examined, and comparable performances are observed for MC and the Nosé-Hoover algorithm. Its ease of use combined with its sampling efficiency make the MC module in CHARMM an attractive alternative for exploring the behavior of biomolecular systems.

Alanine↗

PubMed2005

Controlled subnanosecond isomerization of HCN to CNH in solution.

We report a study of control of the HCN-->CNH isomerization in a liquid Ar solution. We show, using molecular dynamics simulations, nearly complete conversion from HCN to CNH can be achieved in solution on the subnanosecond time scale without requiring laser pulse shaping or molecular alignment. The mechanism of the isomerization reaction involves multiphoton rovibrational excitation on the ground electronic state potential energy surface coupled with rapid rovibrational relaxation in solution. The results demonstrate the important role of rotation-vibration coupling in multiphoton excitation of small molecules and constitute the first realistic computational demonstration of fast, robust, and high-yield laser field manipulation of solution-phase molecular processes.

Journal Article↗

PubMed2005

Automatic method for identifying reaction coordinates in complex systems.

To interpret simulations of a complex system to determine the physical mechanism of a dynamical process, it is necessary to identify the small number of coordinates that distinguish the stable states from the transition states. We develop an automatic method for identifying these degrees of freedom from a database of candidate physical variables. In the method neural networks are used to determine the functional dependence of the probability of committing to a stable state (committor) on a set of coordinates, and a genetic algorithm selects the combination of inputs that yields the best fit. The method enables us to obtain the first set of coordinates that is demonstrably sufficient to specify the transition state of the C(7eq)--> alpha(R) isomerization of the alanine dipeptide in the presence of explicit water molecules. It is revealed that the solute-solvent coupling can be described by a solvent-derived electrostatic torque around one of the main-chain bonds, and the collective, long-ranged nature of this interaction accounts for previous failures to characterize this reaction.

Alanine↗

PubMed2005

Isomerization and dissociation dynamics of HCN in a picosecond infrared laser field: a full-dimensional classical study.

We report a full-dimensional study of the classical dynamics of HCN-->HNC isomerization and of HCN rovibrational dissociation driven by a strong but nonionizing picosecond infrared laser field. The dynamics of the isolated molecule and of the molecule in liquid Ar have both been studied. Our theoretical and numerical results show that when all degrees of freedom are accounted for the field induced molecular dynamics can be totally different from what was found in previous studies, where the HCN molecule is restricted to a plane containing the external field. It is shown that as HCN is driven by an infrared laser field, the rotation of the H atom around the C-N bond provides an important and highly efficient energy absorption mechanism. In the presence of a monochromatic picosecond infrared laser field with an intensity of 10(13) W/cm(2), this energy absorption mechanism generates considerable HCN-->HNC isomerization yield or high rovibrational dissociation yield without molecular preorientation or prealignment. Our study of the field induced isomerization and dissociation dynamics of the same system in liquid Ar shows that the picosecond isomerization dynamics is insignificantly affected by the surrounding atomic liquid whereas the dissociation yield may be greatly suppressed in a high density liquid. The implications of this study for full-dimensional quantum dynamics of multiphoton rovibrational excitation and dissociation of triatomics are briefly discussed.

Journal Article↗

PubMed2004

Multiphonon vibrational relaxation in liquids: should it lead to an exponential-gap law?

The profound differences between solids and liquids notwithstanding, high-frequency vibrational energy relaxation in liquids seems to be well described by assuming that the excess energy is being transferred into discrete overtones of some fundamental intermolecular vibrations-precisely the way it is in crystalline solids. In a solid-state context, this kind of analysis can be used to justify the observation that relaxation rates fall off exponentially with the energy being transferred. Liquids, however, have a substantial degree of disorder, causing their relevant intermolecular spectra to have correspondingly diffuse band edges and large bandwidths. It is therefore not at all obvious what should become of this exponential-gap-law phenomenology. We show in this paper how near exponential-gap-law behavior can still be derived for vibrational energy relaxation in liquids. To do so, we take advantage of the simple dynamics that the high-frequency relaxation has when it is launched from an individual instantaneous configuration. Interestingly, the physically relevant region turns out not to be true asymptotic limit of our formalism, but for realistic liquid parameters the behavior in the physical regime differs only slightly from an exact exponential-gap law and is strikingly independent of the details of the intermolecular spectra.

Journal Article↗