Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: Calculation of structural and dynamic properties

Abstract

Long dynamics simulations were carried out on the B1 immunoglobulin-binding domain of streptococcal protein G (ProtG) and bovine pancreatic trypsin inhibitor (BPTI) using atomistic descriptions of the proteins and a continuum representation of solvent effects. To mimic frictional and random collision effects, Langevin dynamics (LD) were used. The main goal of the calculations was to explore the stability of tens-of-nanosecond trajectories as generated by this molecular mechanics approximation and to analyze in detail structural and dynamical properties. Conformational fluctuations, order parameters, cross correlation matrices, residue solvent accessibilities, pKa values of titratable groups, and hydrogen-bonding (HB) patterns were calculated from all of the trajectories and compared with available experimental data. The simulations comprised over 40 ns per trajectory for ProtG and over 30 ns per trajectory for BPTI. For comparison, explicit water molecular dynamics simulations (EW/MD) of 3 ns and 4 ns, respectively, were also carried out. Two continuum simulations were performed on each protein using the CHARMM program, one with the all-atom PAR22 representation of the protein force field (here referred to as PAR22/LD simulations) and the other with the modifications introduced by the recently developed CMAP potential (CMAP/LD simulations). The explicit solvent simulations were performed with PAR22 only. Solvent effects are described by a continuum model based on screened Coulomb potentials (SCP) reported earlier, i.e., the SCP-based implicit solvent model (SCP-ISM). For ProtG, both the PAR22/LD and the CMAP/LD 40-ns trajectories were stable, yielding C α root mean square deviations (RMSD) of about 1.0 and 0.8 Å respectively along the entire simulation time, compared to 0.8 Å for the EW/MD simulation. For BPTI, only the CMAP/LD trajectory was stable for the entire 30-ns simulation, with a C α RMSD of ≈ 1.4 Å, while the PAR22/LD trajectory became unstable early in the simulation, reaching a C α RMSD of about 2.7 Å and remaining at this value until the end of the simulation; the C α RMSD of the EW/MD simulation was about 1.5 Å. The source of the instabilities of the BPTI trajectories in the PAR22/LD simulations was explored by an analysis of the backbone torsion angles. To further validate the findings from this analysis of BPTI, a 35-ns SCP-ISM simulation of Ubiquitin (Ubq) was carried out. For this protein, the CMAP/LD simulation was stable for the entire simulation time (C α RMSD of ≈1.0 Å), while the PAR22/LD trajectory showed a trend similar to that in BPTI, reaching a C α RMSD of ≈1.5 Å at 7 ns. All the calculated properties were found to be in agreement with the corresponding experimental values, although local deviations were also observed. HB patterns were also well reproduced by all the continuum solvent simulations with the exception of solvent-exposed side chain-side chain (sc-sc) HB in ProtG, where several of the HB interactions observed in the crystal structure and in the EW/MD simulation wer...

Intramolecular Disulfide Bonds of the Prolactin Receptor Short Form Are Required for Its Inhibitory Action on the Function of the Long Form of the Receptor

“…The screened Coulomb potential continuum model was used to represent the electrostatic and hydrophobic effects of the solvent (14,15). The potential energy of the system was represented by the CHARMM force field (version c33b1) (6), with the CMAP correction (23,25) of the all-atom representation (26). The simulations were extended up to 50 ns, preceded by an equilibration phase of 1 ns.…”

TIGER2: An improved algorithm for temperature intervals with global exchange of replicas

“…This threshold value was chosen based on the modulation of pH dependence properties of side chains. [53][54][55] Figure 9 plots the results of SCEF calculated from the lowest potential energy REMD and TIGER2 structures and the native structure and shows that most of the SCEF values of the calculated structures closely follow the values of the native structure except at Glu5, Thr8, and Gly10 in the TIGER2 structure and Asp3 in the REMD structure. Furthermore, all SCEF values are greater than 0.3, indicating that all residues in the chignolin peptide are exposed to solvent.…”

Development and application of hybrid structure based method for efficient screening of ligands binding to G-protein coupled receptors

“…Details of the method as applied to modeling the extra and intra cellular dopamine D2 and D4 receptors have been described elsewhere [Kortagere et al (under preparation)]. This approach used the Screened Coulomb Potential-Implicit Solvent Model (SCP-ISM) [37] to represent solvent effects implemented in CHARMM (ver 32 [58]) and was validated on a number of known structures [38,39]. This multi-step loop modeling protocol utilizes a Monte-Carlo simulated annealing method to explore the various conformations available to a loop segment that is tethered at one end.…”