PubMed Health⌕ Search

Biomedical subjects

Guanglei Cui

Publications and source records attributed to Guanglei Cui.

14 recordsLinked to original sources

Exploring the role of the active site cysteine in human muscle creatine kinase.

All known guanidino kinases contain a conserved cysteine residue that interacts with the non-nucleophilic eta1-nitrogen of the guanidino substrate. Site-directed mutagenesis studies have shown that this cysteine is important, but not essential for activity. In human muscle creatine kinase (HMCK) this residue, Cys283, forms part of a conserved cysteine-proline-serine (CPS) motif and has a pKa about 3 pH units below that of a regular cysteine residue. Here we employ a computational approach to predict the contribution of residues in this motif to the unusually low cysteine pKa. We calculate that hydrogen bonds to the hydroxyl and to the backbone amide of Ser285 would both contribute approximately 1 pH unit, while the presence of Pro284 in the motif lowers the pKa of Cys283 by a further 1.2 pH units. Using UV difference spectroscopy the pKa of the active site cysteine in WT HMCK and in the P284A, S285A, and C283S/S285C mutants was determined experimentally. The pKa values, although consistently about 0.5 pH unit lower, were in broad agreement with those predicted. The effect of each of these mutations on the pH-rate profile was also examined. The results show conclusively that, contrary to a previous report (Wang et al. (2001) Biochemistry 40, 11698-11705), Cys283 is not responsible for the pKa of 5.4 observed in the WT V/K(creatine) pH profile. Finally we use molecular dynamics simulations to demonstrate that, in order to maintain the linear alignment necessary for associative inline transfer of a phosphoryl group, Cys283 needs to be ionized.

Adenosine Diphosphate↗

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↗

X-ray diffraction "fingerprinting" of DNA structure in solution for quantitative evaluation of molecular dynamics simulation.

Solution state x-ray diffraction fingerprinting is demonstrated as a method for experimentally assessing the accuracy of molecular dynamics (MD) simulations. Fourier transforms of coordinate data from MD simulations are used to produce reciprocal space "fingerprints" of atomic pair distance correlations that are characteristic of the ensemble and are the direct numerical analogues of experimental solution x-ray diffraction (SXD). SXD experiments and MD simulations were carried out to test the ability of experiment and simulation to resolve sequence-dependent modifications in helix conformation for B-form DNA. SXD experiments demonstrated that solution-state poly(AT) and poly(A)-poly(T) duplex DNA sequences exist in ensembles close to canonical B-form and B'-form structures, respectively. In contrast, MD simulations analyzed in terms of SXD fingerprints are shown to deviate from experiment, most significantly for poly(A)-poly(T) duplex DNA. Compared with experiment, MD simulation shortcomings were found to include both mismatches in simulated conformer structures and number population within the ensembles. This work demonstrates an experimental approach for quantitatively evaluating MD simulations and other coordinate models to simulate biopolymer structure in solution and suggests opportunities to use solution diffraction data as experimental benchmarks for developing supramolecular force fields optimized for a range of in situ applications.

DNA↗

Critical assessment of quantum mechanics based energy restraints in protein crystal structure refinement.

A critical evaluation of the performance of X-ray refinement protocols using various energy functions is presented using the bovine pancreatic trypsin inhibitor (BPTI) protein. The four potential energy functions we explored include: (1) fully quantum mechanical calculations; (2) one based on an incomplete molecular mechanics (MM) energy function employed in the Crystallography and NMR System (CNS) with empirical parameters developed by Engh and Huber (EH), which lacks electrostatic and attractive van der Waals terms; (3) one based on a complete MM energy function (AMBER ff99 parameter set); and (4) the same as 3, with the addition of a Generalized Born (GB) implicit solvation term. The R, R (free), real space R values of the refined structures and deviations from the original experimental structure were used to assess the relative performance. It was found that at 1 Angstrom resolution the physically based energy functions 1, 3, and 4 performed better than energy function 2, which we attribute to the better representation of key interactions, particularly electrostatics. The observed departures from the experimental structure were similar for the refinements with physically based energy functions and were smaller than the structure refined with EH. A test refinement was also performed with the reflections truncated at a high-resolution cutoff of 2.5 Angstrom and with random perturbations introduced into the initial coordinates, which showed that low-resolution refinements with physically based energy functions held the structure closer to the experimental structure solved at 1 Angstrom resolution than the EH-based refinements.

Animals↗

Computational studies of the farnesyltransferase ternary complex part I: substrate binding.

Farnesyltransferase (FTase) catalyzes the transfer of farnesyl from farnesyl diphosphate (FPP) to cysteine residues at or near the C-terminus of protein acceptors with a CaaX motif (a, aliphatic; X, Met). Farnesylation is a critical modification to many switch proteins involved in the extracellular signal transduction pathway, which facilitates their fixation on the cell membrane where the extracellular signal is processed. The malfunction caused by mutations in these proteins often causes uncontrolled cell reproduction and leads to tumor formation. FTase inhibitors have been extensively studied as potential anticancer agents in recent years, several of which have advanced to different phases of clinical trials. However, the mechanism of this biologically important enzyme has not been firmly established. Understanding how FTase recruits the FPP substrate is the first and foremost step toward further mechanistic investigations and the design of effective FTase inhibitors. Molecular dynamic simulations were carried out on the ternary structure of FTase complexed with the FPP substrate and an acetyl-capped tetrapeptide (acetyl-CVIM), which revealed that the FPP substrate maintains an inactive conformation and the binding of the diphosphate group can be largely attributed to residues R291beta, K164alpha, K294beta, and H248beta. The FPP substrate assumes an extended conformation in the binding site with restricted rotation of the backbone dihedral angles; however, it does not have a well-defined conformation when unbound in solution. This is evident from multinanosecond MD simulations of the FPP substrate in a vacuum and solution. Our conclusion is further supported by theoretical J coupling calculations. Our results on the FPP binding are in good agreement with previous experimental kinetic studies on FTase mutants. The hypothetical conformational activation of the FPP substrate is currently under investigation.

Farnesyltranstransferase↗

A highly pi-stacked organic semiconductor for field-effect transistors based on linearly condensed pentathienoacene.

We present the synthesis and characterization of a fused-ring compound, dithieno[2,3-d:2',3'-d']thieno[3,2-b:4,5-b']dithiophene (pentathienoacene, PTA). In contrast to pentacene, PTA has a larger band gap than most semiconductors used in organic field-effect transistors (OFETs) and therefore is expected to be stable in air. The large pi-conjugated and planar molecular structure of PTA would also form higher molecular orders that are conductive for carrier transport. X-ray diffraction and atomic force microscopy experiments on its films show that the molecules stack in layers with their long axis upright from the surface. X-ray photoelectron spectroscopy suggests that there are no chemical bonds at the PTA/Au interface. OFETs based on the PTA have been constructed, and their performances as p-type semiconductors are also presented. A high mobility of 0.045 cm(2)/V s and an on/off ratio of 10(3) for a PTA OFET have been achieved, demonstrating the potential of PTA for application in future organic electronics.

Electric Conductivity↗

High-performance field-effect transistors based on Langmuir-Blodgett films of cyclo[8]pyrrole.

We demonstrate the field-effect transistors (FETs) made of cyclo[8]pyrrole thin films prepared by the Langmuir-Blodgett (LB) method. The cyclo[8]pyrrole molecule possesses a 30-pi-electron system and narrower highest-occupied molecular orbital-lowest-unoccupied molecular orbital energy gap (0.63 eV), forms a stable, reproducible monolayer at the air-water interface, and transfers onto a substrate with a nearly unity transfer ratio and face-to-face configuration due to its strong pi-pi interaction. The LB films are uniform characterized by atomic force microscopy and in ordered form confirmed by X-ray diffraction. The FET exhibited high performances with one of the highest hole mobilities (0.68 cm2 V(-1) s(-1)) for thin-film transistors and a high on/off ratio, implying a promising material in the FET family.

Journal Article↗

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↗

Third-order nonlinear optical properties of an ultrathin film containing a porphyrin derivative.

An ultrathin nanoscopic multilayer film has been fabricated through the electrostatic layer-by-layer self-assembly of negatively charged 5,10,15,20-tetrakis(3,4,5-trihydroxyphenyl)porphyrin (DHP) and oppositely charged polyethylenimine (PEI). UV-vis spectra showed a continuous and uniform deposition process of PEI and DHP. The film structure was characterized by small-angle X-ray diffraction measurement and AFM images. The nonlinear optical properties of ultrathin film were studied by Z-scan technique with laser duration of 8 ns at a wavelength of 532 nm. The film sample exhibited strong nonlinear saturated absorption and a self-defocusing effect. The nonlinear absorption coefficient and refractive index of the self-assembly ultrathin film are -9.7 x 10(-5) m/W and -7.56 x 10(-12) m(2)/W, respectively.

Journal Article↗

Fabrication of ultrathin films with large third-order nonlinear optical properties.

An ultrathin composite film containing anionic Ag-His complexes (His: L-Histidine) and oppositely charged BH-PPV was fabricated by means of the electrostatic layer-by-layer self-assembly technique. UV/Vis spectra showed a continuous deposition process of Ag-His complexes and BH-PPV. The film structure was characterized by using small-angle X-ray diffraction, AFM, and SEM. The nonlinear optical properties of the ultrathin film were studied by using the Z-scan technique with a laser duration of 8 ns at a wavelength of 532 nm. The film sample exhibited a strong nonlinear saturated absorption, with an alpha2 value of -3.9 x 10(-5) mW(-1) and a self-defocusing effect with an n2 value of -4.78 x 10 (-12) m2W(-1).

Journal Article↗

Inhibition of the bacterial enoyl reductase FabI by triclosan: a structure-reactivity analysis of FabI inhibition by triclosan analogues.

To explore the molecular basis for the picomolar affinity of triclosan for FabI, the enoyl reductase enzyme from the type II fatty acid biosynthesis pathway in Escherichia coli, an SAR study has been conducted using a series of triclosan analogues. Triclosan (1) is a slow, tight-binding inhibitor of FabI, interacting specifically with the E.NAD(+) form of the enzyme with a K(1) value of 7 pM. In contrast, 2-phenoxyphenol (2) binds with equal affinity to the E.NAD(+) (K(1) = 0.5 microM) and E.NADH (K(2) = 0.4 microM) forms of the enzyme and lacks the slow-binding step observed for triclosan. Thus, removal of the three triclosan chlorine atoms reduces the affinity of the inhibitor for FabI by 70,000-fold and removes the preference for the E.NAD(+) FabI complex. 5-Chloro-2-phenoxyphenol (3) is a slow, tight-binding inhibitor of FabI and binds to the E.NAD(+) form of the enzyme (K(1) = 1.1 pM) 7-fold more tightly than triclosan. Thus, while the two ring B chlorine atoms are not required for FabI inhibition, replacement of the ring A chlorine increases binding affinity by 450,000-fold. Given this remarkable observation, the SAR study was extended to the 5-fluoro-2-phenoxyphenol (4) and 5-methyl-2-phenoxyphenol (5) analogues to further explore the role of the ring A substituent. While both 4 and 5 are slow, tight-binding inhibitors, they bind substantially less tightly to FabI than triclosan. Compound 4 binds to both E.NAD(+) and E.NADH forms of the enzyme with K(1) and K(2) values of 3.2 and 240 nM, respectively, whereas compound 5 binds exclusively to the E.NADH enzyme complex with a K(2) value of 7.2 nM. Thus, the ring A substituent is absolutely required for slow, tight-binding inhibition. In addition, pK(a) measurements coupled with simple electrostatic calculations suggest that the interaction of the ring A substituent with F203 is a major factor in governing the affinity of analogues 3-5 for the FabI complex containing the oxidized form of the cofactor.

Anti-Bacterial Agents↗

Photocurrent generation in multilayer self-assembly films fabricated from water-soluble poly(phenylene vinylene).

A novel, water-soluble, cationic PPV derivative poly[(2,5-bis(3-bromotrimethylammoniopropoxy)-phenylene-1,4-divinylene)-alt-1,4-(2,5-bis(2-(2-hydroxyethoxy)ethoxy))phenylene vinylene] (BH-PPV) has been synthesized by a Heck coupling reaction. Multilayered assemblies of the BH-PPV and the sodium salt of hexa(sulfobutyl)fullerenes (C(60)-HS) were fabricated successfully by an alternate deposition technique. The multilayer structures were studied by UV/Vis spectroscopy, small angle X-ray diffraction, and atomic force microscopy. The photoinduced charge transfer property of the self-assembled multilayer film was also measured by a three-electrode cell technique. A steady and rapid cathodic 5.5 microA cm(-2) photocurrent response was measured as the irradiation of the multilayer film was switched on and off. Importantly, the response of on/off cycling is prompt and reproducible. A possible mechanism for the electron-transfer process is proposed.

Journal Article↗

Conformational heterogeneity observed in simulations of a pyrene-substituted DNA.

NMR studies previously carried out for a DNA system with a pyrene-substituted base did not observe NOEs involving the adenine located 5' to the pyrene, and thus the conformation of the adenine was poorly defined in the resulting family of refined structures. However, chemical shift data suggested that an AT base pair may be present. We have carried out fully unrestrained molecular dynamics simulations starting from several members of the family of structures, and these simulations support the existence of an AT base pair for this region. Simulations in both explicit and implicit solvent were carried out, with each converging to either anti or syn conformation for adenine and base pairing in all cases. During these simulations, large and dramatic conformational changes are observed that suggest pathways for complex conformational changes in the highly packed DNA interior. Our analysis reveals little difference in the energies of these syn and anti conformations, in contrast to control calculations carried out for standard DNA (in the absence of a neighboring pyrene). While no interconversion between the conformations was observed in standard simulations, reversible anti/syn exchange was directly simulated using the locally enhanced sampling approach. No exchange was seen in the non-pyrene control sequence. Together, these results suggest that an increased flexibility is introduced as a consequence of the pyrene substitution, offering an explanation that is consistent with the available experimental data. These results increase our optimism that simulations in atomic detail may provide accurate models for experimental observations in complex systems.

Computer Simulation↗