Biomedical subjects

Lijiang Yang

Publications and source records attributed to Lijiang Yang.

3 recordsLinked to original sources

PubMed2006

How well does Poisson-Boltzmann implicit solvent agree with explicit solvent? A quantitative analysis.

We have quantitatively studied the performance of a finite-difference Poisson-Boltzmann implicit solvent with respect to the TIP3P explicit solvent in a range of systems of biochemical interest. An overall agreement was found between the tested implicit and explicit solvents for hydrogen-bonding/salt-bridging dimers and peptide monomers and dimers of different conformations and different lengths. These comparative analyses also indicate a good transferability of empirically optimized parameters for the implicit solvent from small training molecules to large testing peptides. However, deviations between the two tested solvents are also apparent. Specifically, a consistent deviation was observed when hydrogen-bonding or salt-bridging dimers are within 4-6 A. The deviation reaches a maximum at about 5.5 A, the so-called water-bridging distance. The tested implicit solvent, even with optimized parameters, cannot capture the subtle fluctuation in the distance-dependent reaction field energy profiles, although smoothed profiles can still be obtained and are in overall agreement with those in the explicit solvent. Interestingly, the same mechanism underlining the above discrepancy is also responsible for the larger deviations of certain peptide conformations, such as parallel beta-strand dimers. It is likely that the observed discrepancy may cause improper conformational distributions in simulations with the implicit solvent when hydrogen-bonding or salt-bridging interactions are crucial, such as secondary structure populations in proteins. Validation of the implicit solvent with optimized parameters in dynamics simulations will be the next step to study the influences of the observed discrepancy at biological conditions.

Alanine↗

PubMed2006

On the enhanced sampling over energy barriers in molecular dynamics simulations.

We present here calculations of free energies of multidimensional systems using an efficient sampling method. The method uses a transformed potential energy surface, which allows an efficient sampling of both low and high energy spaces and accelerates transitions over barriers. It allows efficient sampling of the configuration space over and only over the desired energy range(s). It does not require predetermined or selected reaction coordinate(s). We apply this method to study the dynamics of slow barrier crossing processes in a disaccharide and a dipeptide system.

Journal Article↗

PubMed2006

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↗