PubMed Health⌕ Search

Biomedical subjects

Yong Duan

Publications and source records attributed to Yong Duan.

At least 19 recordsLinked to original sources

Folding transition-state and denatured-state ensembles of FSD-1 from folding and unfolding simulations.

Characterization of the folding transition-state ensemble and the denatured-state ensemble is an important step toward a full elucidation of protein folding mechanisms. We report herein an investigation of the free-energy landscape of FSD-1 protein by a total of four sets of folding and unfolding molecular dynamics simulations with explicit solvent. The transition-state ensemble was initially identified from unfolding simulations at 500 K and was verified by simulations at 300 K starting from the ensemble structures. The denatured-state ensemble and the early-stage folding were studied by a combination of unfolding simulations at 500 K and folding simulations at 300 K starting from the extended conformation. A common feature of the transition-state ensemble was the substantial formation of the native secondary structures, including both the alpha-helix and beta-sheet, with partial exposure of the hydrophobic core in the solvent. Both the native and non-native secondary structures were observed in the denatured-state ensemble and early-stage folding, consistent with the smooth experimental melting curve. Interestingly, the contact orders of the transition-state ensemble structures were similar to that of the native structure and were notably lower than those of the compact structures found in early-stage folding, implying that chain and topological entropy might play significant roles in protein folding. Implications for FSD-1 folding mechanisms and the rate-limiting step are discussed. Analyses further revealed interesting non-native interactions in the denatured-state ensemble and early-stage folding and the possibility that destabilization of these interactions could help to enhance the stability and folding rate of the protein.

DNA-Binding Proteins↗

Phenol red interacts with the protofibril-like oligomers of an amyloidogenic hexapeptide NFGAIL through both hydrophobic and aromatic contacts.

Amyloid-associated diseases affect millions of people worldwide. Phenol red exhibits modest inhibition toward fibril formation of human Islet amyloid polypeptide (hIAPP) and its toxicity, which is associated with type II diabetes mellitus. However, the molecular level mechanisms of interactions remain elusive. The binding of phenol red molecules to the protofibrils of an amyloidogenic fragment (NFGAIL) of hIAPP has been investigated by molecular dynamics simulations with explicit solvent. The phenol red molecules were observed to bind primarily along either beta-sheet stacking or beta-strand directions. Through its three aromatic rings, the phenol red molecule preferentially interacted with the hydrophobic side chains of Phe, Leu, and Ile; and the polar sulfone and hydroxyl groups were mainly exposed in solvent. Thus, phenol red improves the solubility of the early protofibrils and represses further growth. Interestingly, there was no obvious preference toward the aromatic Phe residue in comparison to the hydrophobic Leu or Ile residues. The lack of binding along the hydrogen bond direction indicates that phenol red does not directly block the beta-sheet extension. Further free energy analysis suggested that a phenol red analog may potentially improve the binding affinity.

Amyloid↗

New-generation amber united-atom force field.

We have developed a new-generation Amber united-atom force field for simulations involving highly demanding conformational sampling such as protein folding and protein-protein binding. In the new united-atom force field, all hydrogens on aliphatic carbons in all amino acids are united with carbons except those on Calpha. Our choice of explicit representation of all protein backbone atoms aims at minimizing perturbation to protein backbone conformational distributions and to simplify development of backbone torsion terms. Tests with dipeptides and solvated proteins show that our goal is achieved quite successfully. The new united-atom force field uses the same new RESP charging scheme based on B3LYP/cc-pVTZ//HF/6-31g** quantum mechanical calculations in the PCM continuum solvent as that in the Duan et al. force field. van der Waals parameters are empirically refitted starting from published values with respect to experimental solvation free energies of amino acid side-chain analogues. The suitability of mixing new point charges and van der Waals parameters with existing Amber covalent terms is tested on alanine dipeptide and is found to be reasonable. Parameters for all new torsion terms are refitted based on the new point charges and the van der Waals parameters. Molecular dynamics simulations of three small globular proteins in the explicit TIP3P solvent are performed to test the overall stability and accuracy of the new united-atom force field. Good agreements between the united-atom force field and the Duan et al. all-atom force field for both backbone and side-chain conformations are observed. In addition, the per-step efficiency of the new united-atom force field is demonstrated for simulations in the implicit generalized Born solvent. A speedup around two is observed over the Duan et al. all-atom force field for the three tested small proteins. Finally, the efficiency gain of the new united-atom force field in conformational sampling is further demonstrated with a well-known toy protein folding system, an 18 residue polyalanine in distance-dependent dielectric. The new united-atom force field is at least a factor of 200 more efficient than the Duan et al. all-atom force field for ab initio folding of the tested peptide.

Models, Molecular↗

Strike a balance: optimization of backbone torsion parameters of AMBER polarizable force field for simulations of proteins and peptides.

Based on the AMBER polarizable model (ff02), we have re-optimized the parameters related to the main-chain (Phi, Psi) torsion angles by fitting to the Boltzmann-weighted average quantum mechanical (QM) energies of the important regions (i.e., beta, P(II), alpha(R), and alpha(L) regions). Following the naming convention of the AMBER force field series, this release will be called ff02pol.rl The force field has been assessed both by energetic comparison against the QM data and by the replica exchange molecular dynamics simulations of short alanine peptides in water. For Ace-Ala-Nme, the simulated populations in the beta, P(II) and alpha(R) regions were approximately 30, 43, and 26%, respectively. For Ace-(Ala)(7)-Nme, the populations in these three regions were approximately 24, 49, and 26%. Both were in qualitative agreement with the NMR and CD experimental conclusions. In comparison with the previous force field, ff02pol.rl demonstrated good balance among these three important regions. The optimized torsion parameters, together with those in ff02, allow us to carry out simulations on proteins and peptides with the consideration of polarization.

Computer Simulation↗

Grow to Fit Molecular Dynamics (G2FMD): an ab initio method for protein side-chain assignment and refinement.

The rough energy landscapes and tight packing of protein interiors are two of the critical factors that have prevented the wide application of physics-based models in protein side-chain assignment and protein structure prediction in general. Complementing the rotamer-based methods, we propose an ab initio method that utilizes molecular mechanics simulations for protein side-chain assignment and refinement. By reducing the side-chain size, a smooth energy landscape was obtained owing to the increased distances between the side chains. The side chains then gradually grow back during molecular dynamics simulations while adjusting to their surrounding driven by the interaction energies. The method overcomes the barriers due to tight packing that limit conformational sampling of physics-based models. A key feature of this approach is that the resulting structures are free from steric collisions and allow the application of all-atom models in the subsequent refinement. Tests on a small set of proteins showed nearly 100% accuracy on both chi1 and chi2 of buried residues and 94% of them were within 20 degrees from the native conformation, 79% were within 10 degrees and 42% were within 5 degrees . However, the accuracy decreased when exposed side chains were involved. Further improvement and application of the method and the possible reasons that affect the accuracy on the exposed side chains are discussed.

Chemical Phenomena↗

[Study on manpower allocation criteria of center of disease prevention and control in the context of province level, city level and county level].

OBJECTIVE: To set the manpower allocation criteria of center of disease prevention and control. METHODS: Expected allocation manpower criteria was obtained through adjusting the current manpower allocation of disease prevention and control centers. The principle was to fulfill public function and promote professional efficiency. RESULTS: Based on function requirement, in 3 - 5 years, the manpower allocation criteria of center of disease prevention and control at provincial-level is 336 persons, at city-level is 102 persons, and at county-level is 33 persons, that means in whole country 140016 persons should be needed. In 10 years, the manpower allocation criteria of center of disease prevention and control at provincial-level is 386 persons, at city-level is 112 persons, and at county-level is 38 persons, that means in whole country 159086 persons should be needed. CONCLUSION: The manpower allocation criteria advanced in the study indicated that current manpower quantity should be greatly reduced. It is an inevitable trend that disease prevention and control centers reduce the staff quantity and promote their quality.

China↗

[Inhibited proliferation of B-lymphoma Raji cells and down-regulated expression of VEGF by arsenic trioxide].

This study was aimed to investigate the effects of As2O3 on proliferation of B lymphoma Raji cell and to study the expression changes of VEGF mRNA. The modified MTT was adopted to evaluate the effect of As2O3 on proliferation of Raji cells. Semi-quantitative RT-PCR was used to detect the expression of VEGF121 and VEGF165 mRNA in Raji cells exposed to As2O3. The results showed that the As2O3 significantly inhibited the proliferation of Raji cells, the relationship between the inhibition rate and the concentrations of As2O3 was dose- and time-dependent. The VEGF121 and VEGF165 mRNA expressions were found in Raji cells, and both were well-matched in expression levels. After treatment of As2O3 at 1 micromol/L for 48 hours, the expression levels of VEGF121 and VEGF165 mRNA were significantly down-regulated and demonstrated negative correlation with the time of exposure to As2O3. Otherwise some difference were observed in effects of As2O3 on VEGF121 and VEGF165, the expression of VEGF165 mRNA was down-regulated earlier and longer than that of VEGF121, while no similar correlation was found in selected concentration groups. It is concluded that As2O3 can significantly inhibit the growth of Raji cells and may exerted its anti-angiogenesis effects by down-regulating the VEGF mRNA expression even in a low concentration for a long term.

Antineoplastic Agents↗

Molecular dynamics simulations and free energy analyses on the dimer formation of an amyloidogenic heptapeptide from human beta2-microglobulin: implication for the protofibril structure.

Amyloid formation is associated with many neurodegenerative diseases. Recent findings suggest that early oligomeric aggregates could be major sources of toxicity. We present a computational investigation of the first step of amyloid initiation-dimer formation of a seven residue peptide (NHVTLSQ) from human beta2-microglobulin at pH 2.0, which renders +2.0 units charges to each peptide. A total of over 1.2 micros of simulations with explicit solvent and 1.0 micros of simulations with implicit solvent were conducted. Main-chain conformational restraint was applied to facilitate the formation of ordered dimers. An antiparallel beta-sheet with six main-chain hydrogen bonds was dominant in the implicit solvent simulations. In contrast, no stable dimers were observed in the two negative controls, the mouse heptapeptide (KHDSMAE, +3.0 units charges) and the scrambled human heptapeptide (QVLHTSN). Explicit solvent simulations presented a more complex scenario. The wild-type human heptapeptide formed predominantly antiparallel beta-sheets ( approximately 38%) although parallel ones ( approximately 12%) were also observed. Hydrophobic contacts preceded hydrogen bond saturation in the majority of the association events in the explicit solvent simulations, highlighting the important role of hydrophobic interaction in amyloid initiation. The fact that the mouse dimer dissociated immediately after the removal of conformational restraint suggests that the higher conformational entropy barrier, along with the stronger charge repulsion and weaker hydrophobic interaction, contributed to its inability to form amyloid fibril. The closeness of positive charge pairs in the dimers of the scrambled human heptapeptide may prohibit further beta-sheet extension and fibril growth. Combining the results from simulations and free energy analyses, we propose that the building block for this amyloid fibril is an antiparallel dimer with a two-residue register shift and six main-chain hydrogen bonds. A double-layer protofibril structure is also proposed in which two antiparallel beta-sheets face each other and are held together by hydrophobic staples and hydrogen bonds of the polar side-chains.

Amino Acid Sequence↗

Multiple roles of a conserved GAF domain tyrosine residue in cyanobacterial and plant phytochromes.

The phytochrome family of red/far-red photoreceptors has been optimized to support photochemical isomerization of a bound bilin chromophore, a process that triggers a conformational change and modulates biochemical output from the surrounding protein scaffold. Recent studies have established that the efficiency of this photochemical process is profoundly altered by mutation of a conserved tyrosine residue (Tyr176) within the bilin-binding GAF domain of the cyanobacterial phytochrome Cph1 [Fischer, A. J., and Lagarias, J. C. (2004) Harnessing phytochrome's glowing potential, Proc. Natl. Acad. Sci. U.S.A. 101, 17334-17339]. Here, we show that the equivalent mutation in plant phytochromes behaves similarly, indicating that the function of this tyrosine in the primary photochemical mechanism is conserved. Saturation mutagenesis of Tyr176 in Cph1 establishes that no other residue can support comparably efficient photoisomerization. The spectroscopic consequences of Tyr176 mutations also reveal that Tyr176 regulates the conversion of the porphyrin-like conformation of the bilin precursor to a more extended conformation. The porphyrin-binding ability of the Tyr176Arg mutant protein indicates that Tyr176 also regulates the ligand-binding specificity of apophytochrome. On the basis of the hydrogen-bonding ability of Tyr176 substitutions that support the nonphotochemical C15-Z,syn to C15-Z,anti interconversion, we propose that Tyr176 orients the carboxyl side chain of a conserved acidic residue to stabilize protonation of the bilin chromophore. A homology model of the GAF domain of Cph1 predicts a C5-Z,syn, C10-Z,syn, C15-Z,anti configuration for the chromophore and implicates Glu189 as the proposed acidic residue stabilizing the extended conformation, an interpretation consistent with site-directed mutagenesis of this conserved acidic residue.

Amino Acid Substitution↗

Convergence of replica exchange molecular dynamics.

Replica exchange molecular dynamics (REMD) method is one of the generalized-ensemble algorithms which performs random walk in energy space and helps a system to escape from local energy traps. In this work, we studied the accuracy and efficiency of REMD by examining its ability to reproduce the results of multiple extended conventional molecular dynamics (MD) simulations and to enhance conformational sampling. Two sets of REMD simulations with different initial configurations, one from the fully extended and the other from fully helical conformations, were conducted on a fast-folding 21-amino-acid peptide with a continuum solvent model. Remarkably, the two REMD simulation sets started to converge even within 1.0 ns, despite their dramatically different starting conformations. In contrast, the conventional MD within the same time and with identical starting conformations did not show obvious signs of convergence. Excellent convergence between the REMD sets for T>300 K was observed after 14.0 ns REMD simulations as measured by the average helicity and free-energy profiles. We also conducted a set of 45 MD simulations at nine different temperatures with each trajectory simulated to 100.0 and 200.0 ns. An excellent agreement between the REMD and the extended MD simulation results was observed for T>300 K, showing that REMD can accurately reproduce long-time MD results with high efficiency. The autocorrelation times of the calculated helicity demonstrate that REMD can significantly enhance the sampling efficiency by 14.3+/-6.4, 35.1+/-0.2, and 71.5+/-20.4 times at, respectively, approximately 360, approximately 300, and approximately 275 K in comparison to the regular MD. Convergence was less satisfactory at low temperatures (T<300 K) and a slow oscillatory behavior suggests that longer simulation time was needed to reach equilibrium. Other technical issues, including choice of exchange frequency, were also examined.

Algorithms↗

Elongation of ordered peptide aggregate of an amyloidogenic hexapeptide NFGAIL observed in molecular dynamics simulations with explicit solvent.

The mechanisms by which amyloidogenic peptides and proteins form soluble toxic oligomers remain elusive. We have studied the formation of partially ordered tetramers and well-ordered octamers of an amyloidogenic hexapeptide NFGAIL (residues 22-27 of the human islet amyloid polypeptide) in our previous work. Continuing the effort, we here probe the beta-sheet elongation process by a combined total of 2.0 micros molecular dynamics simulations with explicit solvent. In a set of 10 simulations with the peptides restrained to the extended conformation, we observed that the main growth mode was elongation along the beta-sheet hydrogen bonds through primarily a two-stage process. Driven by hydrophobic forces, the peptides initially attached to the surface of the ordered oligomer, moved quickly to the beta-sheet edges, and formed stable beta-sheet hydrogen bonds. Addition of peptides to the existing oligomer notably improved the order of the peptide aggregate in which labile outer layer beta-sheets were stabilized, which provides good templates for further elongation. These simulations suggested that elongation along the beta-sheet hydrogen bonds occurs at the intermediate stage when low-weight oligomers start to form. We did not observe significant preference toward either parallel or antiparallel beta-sheets at the elongation stage for this peptide. In another set of 10 unrestrained simulations, the dominant growth mode was disordered aggregation. Taken together, these results offered a glimpse at the molecular events leading to the formation of ordered and disordered low-weight oligomers.

Amyloid↗

Statistical characterization of salt bridges in proteins.

The structure and folding mechanism of a given protein are determined by many factors, including the electrostatic interactions between charged residues of protein molecules known in general as salt bridges. In this study, analyses were conducted on 10,370 salt bridges in 2017 proteins and the results compared to previous statistical surveys of 36 protein structures. Although many of the general trends remained consistent with other studies, more detailed information was illuminated by the larger dataset. In particular, it was shown that there is a strong correlation between secondary structure and salt bridge formation, and that salt bridges display preferential formation in an environment of about 30% solvent accessible surface area.

Amino Acids↗

Comparison between Generalized-Born and Poisson-Boltzmann methods in physics-based scoring functions for protein structure prediction.

Continuum solvent models such as Generalized-Born and Poisson-Boltzmann methods hold the promise to treat solvation effect efficiently and to enable rapid scoring of protein structures when they are combined with physics-based energy functions. Yet, direct comparison of these two approaches on large protein data set is lacking. Building on our previous work with a scoring function based on a Generalized-Born (GB) solvation model, and short molecular-dynamics simulations, we further extended the scoring function to compare with the MM-PBSA method to treat the solvent effect. We benchmarked this scoring function against seven publicly available decoy sets. We found that, somewhat surprisingly, the results of MM-PBSA approach are comparable to the previous GB-based scoring function. We also discussed the effect to the scoring function accuracy due to presence of large ligands and ions in some native structures of the decoy sets.

Chemical Phenomena↗

Denatured-state ensemble and the early-stage folding of the G29A mutant of the B-domain of protein A.

The folding mechanism of the G29A mutant of the B-domain of protein A (BdpA) has been studied by all-atom molecular dynamics simulation using AMBER force field (ff03) and generalized Born continuum solvent model. Started from the extended chain conformation, a total of 16 simulations (400 ns each) at 300 K captured some early folding events of the G29A mutant of BdpA. In one of the 16 trajectories, the G29A mutant folded within 2.8 A (root mean square) of the wild-type NMR structure. We observed that the fast burial of hydrophobic residues was the driving force to bring the distant residues into close proximity. The initiation of the helix I and III occurred during the stage of hydrophobic collapse. The initiation and growth of the helix II was slow. Both the secondary structure formation and the development of the native tertiary contacts suggested a multistage folding process. Clustering analysis indicated that two helix species (helices I and III) could be intermediates. Further analysis revealed that the hydrophobic residues of partially folded helix II formed nativelike hydrophobic contacts with helices I and III that stabilized a nativelike state and delayed the completion of folding of the entire protein. The details of the early folding process were compared with other theoretical and experimental studies. It was found that a nativelike hydrophobic cluster was formed by residues including F(30), I(31), L(34), L(44), L(45), and A(48) that prevented further development of the native structures, and breaking the hydrophobic cluster like this one contributed to the rate-limiting step. This was in complete agreement with the recent kinetic measurements in which mutations of these residues to Gly and Ala substantially increased the folding rates by as much as 60 times. Apparently, destabilization of nonnative states dramatically enhanced the folding rates.

Cluster Analysis↗

Loop conformation and dynamics of the Escherichia coli HPPK apo-enzyme and its binary complex with MgATP.

Comparison of the crystallographic and NMR structures of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) suggests that the enzyme may undergo significant conformational change upon binding to its first substrate, ATP. Two of the three surface loops (loop 2 and loop 3) accounting for most of the conformational differences appear to be confined by crystal contacts, raising questions about the putative large-scale induced-fit conformational change of HPPK and the functional roles of the conserved side-chain residues on the loops. To investigate the loop dynamics in crystal-free environment, we carried out molecular dynamics and locally enhanced sampling simulations of the apo-enzyme and the HPPK.MgATP complex. Our simulations showed that the crystallographic B-factors underestimated the loop dynamics considerably. We found that the open-conformation of loop 3 in the binary complex is accessible to the apo-enzyme and is the favored conformation in solution phase. These results revise our previous view of HPPK-substrate interactions and the associated functional mechanism of conformational change. The lessons learned here offer valuable structural insights into the workings of HPPK and should be useful for structure-based drug design.

Adenosine Triphosphate↗

Large-scale conformational dynamics of the HIV-1 integrase core domain and its catalytic loop mutants.

HIV-1 integrase is one of the three essential enzymes required for viral replication and has great potential as a novel target for anti-HIV drugs. Although tremendous efforts have been devoted to understanding this protein, the conformation of the catalytic core domain around the active site, particularly the catalytic loop overhanging the active site, is still not well characterized by experimental methods due to its high degree of flexibility. Recent studies have suggested that this conformational dynamics is directly correlated with enzymatic activity, but the details of this dynamics is not known. In this study, we conducted a series of extended-time molecular dynamics simulations and locally enhanced sampling simulations of the wild-type and three loop hinge mutants to investigate the conformational dynamics of the core domain. A combined total of >480 ns of simulation data was collected which allowed us to study the conformational changes that were not possible to observe in the previously reported short-time molecular dynamics simulations. Among the main findings are a major conformational change (>20 A) in the catalytic loop, which revealed a gatinglike dynamics, and a transient intraloop structure, which provided a rationale for the mutational effects of several residues on the loop including Q(148), P(145), and Y(143). Further, clustering analyses have identified seven major conformational states of the wild-type catalytic loop. Their implications for catalytic function and ligand interaction are discussed. The findings reported here provide a detailed view of the active site conformational dynamics and should be useful for structure-based inhibitor design for integrase.

Binding Sites↗

Incorporating intermolecular distance into protein-protein docking.

In this work, intermolecular distance was integrated into the docking of protein-protein complexes. To develop an efficient docking procedure, 22 enzyme-inhibitor targets and 15 antibody-antigen targets were taken from a benchmark set. A three-step approach was adopted, which included global sampling by FTDOCK, filtering by intermolecular distance and ranking by a composite scoring function. For the enzyme-inhibitor targets, the composite scoring function consists of geometry and energy terms. In the set composed of the approximately 100 highest ranked candidates for each target, correct complexes were identified for all of the 22 enzyme-inhibitor targets. This docking strategy also succeeded on the four test targets, of which three are CAPRI targets with the same receptor but different binding modes. Interestingly, all three binding modes were correctly predicted. For the antibody-antigen targets, CDR and physical energy were also used in the filtering process and informatics terms were added to the scoring function. The composite score had successful prediction for 13 of the 15 antibody-antigen targets.

Algorithms↗