Biomedical subjects

Xiaolin Cheng

Publications and source records attributed to Xiaolin Cheng.

7 recordsLinked to original sources

PubMed2006

Order N algorithm for computation of electrostatic interactions in biomolecular systems.

Poisson-Boltzmann electrostatics is a well established model in biophysics; however, its application to large-scale biomolecular processes such as protein-protein encounter is still limited by the efficiency and memory constraints of existing numerical techniques. In this article, we present an efficient and accurate scheme that incorporates recently developed numerical techniques to enhance our computational ability. In particular, a boundary integral equation approach is applied to discretize the linearized Poisson-Boltzmann equation; the resulting integral formulas are well conditioned and are extended to systems with arbitrary numbers of biomolecules. The solution process is accelerated by Krylov subspace methods and a new version of the fast multipole method. In addition to the electrostatic energy, fast calculations of the forces and torques are made possible by using an interpolation procedure. Numerical experiments show that the implemented algorithm is asymptotically optimal O(N) in both CPU time and required memory, and application to the acetylcholinesterase-fasciculin complex is illustrated.

Acetylcholinesterase↗

PubMed2006

Targeted molecular dynamics study of C-loop closure and channel gating in nicotinic receptors.

The initial coupling between ligand binding and channel gating in the human alpha7 nicotinic acetylcholine receptor (nAChR) has been investigated with targeted molecular dynamics (TMD) simulation. During the simulation, eight residues at the tip of the C-loop in two alternating subunits were forced to move toward a ligand-bound conformation as captured in the crystallographic structure of acetylcholine binding protein (AChBP) in complex with carbamoylcholine. Comparison of apo- and ligand-bound AChBP structures shows only minor rearrangements distal from the ligand-binding site. In contrast, comparison of apo and TMD simulation structures of the nAChR reveals significant changes toward the bottom of the ligand-binding domain. These structural rearrangements are subsequently translated to the pore domain, leading to a partly open channel within 4 ns of TMD simulation. Furthermore, we confirmed that two highly conserved residue pairs, one located near the ligand-binding pocket (Lys145 and Tyr188), and the other located toward the bottom of the ligand-binding domain (Arg206 and Glu45), are likely to play important roles in coupling agonist binding to channel gating. Overall, our simulations suggest that gating movements of the alpha7 receptor may involve relatively small structural changes within the ligand-binding domain, implying that the gating transition is energy-efficient and can be easily modulated by agonist binding/unbinding.

Amino Acid Sequence↗

PubMed2006

Insight through molecular mechanics Poisson-Boltzmann surface area calculations into the binding affinity of triclosan and three analogues for FabI, the E. coli enoyl reductase.

Keeping pace with emerging drug resistance in clinically important pathogens will be greatly aided by inexpensive yet reliable computational methods that predict the binding affinities of ligands for drug targets. We present results using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method to calculate the affinity of a series of triclosan analogues for the E. coli enoyl reductase FabI, spanning a 450000-fold range of binding affinities. Significantly, a high correlation is observed between the calculated binding energies and those determined experimentally. Further examination indicates that the van der Waals energies are the most correlated component of the total affinity (r2 = 0.74), indicating that the shape of the inhibitor is very important in defining the binding energies for this system. The validation of MM-PBSA for the E coli FabI system serves as a platform for inhibitor design efforts focused on the homologous enzyme in Staphylococcus aureus and Mycobacterium tuberculosis.

Binding Sites↗

PubMed2005

Channel opening motion of alpha7 nicotinic acetylcholine receptor as suggested by normal mode analysis.

The gating motion of the human nicotinic acetylcholine receptor (nAChR) alpha7 was investigated with normal mode analysis (NMA) of two homology models. The first model, referred to as model I, was built from both the Lymnaea stagnalis acetylcholine binding protein (AChBP) and the transmembrane (TM) domain of the Torpedo marmorata nAChR. The second model, referred to as model C, was based solely on the recent electron microscopy structure of the T. marmorata nAChR. Despite structural differences, both models exhibit nearly identical patterns of flexibility and correlated motions. In addition, both models show a similar global twisting motion that may represent channel gating. The similar results obtained for the two models indicate that NMA is most sensitive to the contact topology of the structure rather than its finer detail. The major difference between the low-frequency motions sampled for the two models is that a symmetrical pore-breathing motion, favoring channel opening, is present as the second most dominant motion in model I, whilst largely absent from model C. The absence of this mode in model C can be attributed to its less symmetrical architecture. Finally, as a further goal of the present study, an approximate open channel model, consistent with many experimental findings, has been produced.

Humans↗

PubMed2005

Dynamic behavior of DNA base pairs containing 8-oxoguanine.

The process by which DNA repair enzymes recognize and selectively excise damaged bases in duplex DNA is fundamental to our mechanistic understanding of these critical biological reactions. 8-Oxoguanine (8-oxoG) is the most common form of oxidative DNA damage; unrepaired, this lesion generates a G:C-->T:A mutation. Central to the recognition and repair of DNA damage is base extrusion, a process in which the damaged base lesion or, in some cases, its partner disengages from the helix and is bound to the enzyme's active site where base excision takes place. The conformation adopted by 8-oxoG in duplex DNA is affected by the base positioned opposite this lesion; conformational changes may also take place when the damaged base binds to its cognate repair enzyme. We performed unrestrained molecular dynamics simulations for several 13-mer DNA duplexes. Oligomers containing G:C and 8oxoG:C pairs adopted Watson-Crick geometries in stable B-form duplexes; 8oxoG showed increased local and global flexibility and a reduced barrier to base extrusion. Duplexes containing the G:A mismatch showed much larger structural fluctuations and failed to adopt a well-defined structure. For the 8oxoG:A mismatch that is recognized by the DNA glycosylase MutY, the damaged nucleoside underwent spontaneous and reproducible anti-->syn transitions. The syn conformation is thermodynamically preferred. Steric hindrance and unfavorable electrostatics associated with the 8oxoG O8 atom in the anti conformation were the major driving forces for this transition. Transition events follow two qualitatively different pathways. The overall anti-->syn transition rate and relative probability of the two transition paths were dependent on local sequence context. These simulations indicate that both the dynamic and equilibrium behavior of the duplex change as a result of oxidation; these differences may provide valuable new insight into the selective action of enzymes on damaged DNA.

Base Pairing↗

PubMed2005

Calculation of the Maxwell stress tensor and the Poisson-Boltzmann force on a solvated molecular surface using hypersingular boundary integrals.

The electrostatic interaction among molecules solvated in ionic solution is governed by the Poisson-Boltzmann equation (PBE). Here the hypersingular integral technique is used in a boundary element method (BEM) for the three-dimensional (3D) linear PBE to calculate the Maxwell stress tensor on the solvated molecular surface, and then the PB forces and torques can be obtained from the stress tensor. Compared with the variational method (also in a BEM frame) that we proposed recently, this method provides an even more efficient way to calculate the full intermolecular electrostatic interaction force, especially for macromolecular systems. Thus, it may be more suitable for the application of Brownian dynamics methods to study the dynamics of protein/protein docking as well as the assembly of large 3D architectures involving many diffusing subunits. The method has been tested on two simple cases to demonstrate its reliability and efficiency, and also compared with our previous variational method used in BEM.

Algorithms↗

PubMed2005

Modified replica exchange simulation methods for local structure refinement.

Parallel tempering, also known as replica exchange molecular dynamics (REMD), has recently been successfully used to study the structure and thermodynamic properties of biomolecules such as peptides and small proteins. For large systems, however, applying REMD can be costly since the number of replicas needed increases as the square root of the number of degrees of freedom in the system. Often, enhanced sampling is only needed for a subset of atoms, such as a loop region of a large protein or a small ligand binding to a receptor. In such applications, it is often reasonable to assume a weak dependence of the structure of the larger region on the instantaneous conformation of the smaller region of interest. For these cases, we derived two variant replica exchange methods, partial replica exchange molecular dynamics (PREMD) and local replica exchange molecular dynamics (LREMD). The Hamiltonian for the system is separated, with replica exchange carried out only for terms involving the subsystem of interest while the remainder of the system is maintained at a single temperature. The number of replicas required for efficient exchange thus depends on the number of degrees of freedom in the fragment needing refinement rather than on the size of the full system. The method can be applied to much larger systems than was previously practical. This also provides a means to preserve the integrity of the structure outside the refinement region without introduction of restraints. LREMD takes this weak coupling approximation a step further, employing only a single representation of the large fragment that simultaneously interacts with all of the replicas of the subsystem of interest. This is obtained by combining replica exchange with the locally enhanced sampling approximation (LES), reducing the computational expense of replica exchange simulations to near that of a single standard molecular dynamics (MD) simulation. Use of LREMD also permits the use of LES without requiring the specification of a single temperature, a known difficulty for standard LES simulations. We tested these two methods on the loop region of an RNA hairpin model system and find significant advantages over standard MD and REMD simulations.

Algorithms↗