The double cubic lattice method: Efficient approaches to numerical integration of surface area and volume and to dot surface contouring of molecular assemblies

Molecular dynamics simulation of interactions between a sodium dodecyl sulfate micelle and a poly(ethylene oxide) polymer

We have performed atomistic molecular dynamics simulations of an anionic sodium dodecyl sulfate (SDS) micelle and a nonionic poly(ethylene oxide) (PEO) polymer in aqueous solution. The micelle consisted of 60 surfactant molecules, and the polymer chain lengths varied from 20 to 40 monomers. The force field parameters for PEO were adjusted by using 1,2-dimethoxymethane (DME) as a model compound and matching its hydration enthalpy and conformational behavior to experiment. Excellent agreement with previous experimental and simulation work was obtained through these modifications. The simulated scaling behavior of the PEO radius of gyration was also in close agreement with experimental results. The SDS-PEO simulations show that the polymer resides on the micelle surface and at the hydrocarbon-water interface, leading to a selective reduction in the hydrophobic contribution to the solvent-accessible surface area of the micelle. The association is mainly driven by hydrophobic interactions between the polymer and surfactant tails, while the interaction between the polymer and sulfate headgroups on the micelle surface is weak. The 40-monomer chain is mostly wrapped around the micelle, and nearly 90% of the monomers are adsorbed at low PEO concentration. Simulations were also performed with multiple 20-monomer chains, and gradual addition of polymer indicates that about 120 monomers are required to saturate the micelle surface. The stoichiometry of the resulting complex is in close agreement with experimental results, and the commonly accepted "beaded necklace" structure of the SDS-PEO complex is recovered by our simulations.

The conformation of the extracellular binding domain of Death Receptor 5 in the presence and absence of the activating ligand TRAIL: a molecular dynamics study

The Death Receptor 5 (DR5), a member of tumor necrosis factor receptor (TNFR) superfamily of receptors, triggers apoptosis (programmed cell death) when stimulated by its tridentate ligand TRAIL. Until recently it was generally assumed that the activation of DR5 resulted from the recruitment of three independent receptor units, leading to the trimerization of intracellular domains. However, there is mounting evidence to suggest that, in the absence of ligand, such cytokine receptors primarily reside as preformed complexes. In this work, molecular dynamics simulations of the TRAIL-DR5 complex, the unbound receptor trimer and individual receptor monomers are compared to gain insight in the mechanism of activation. The results suggest that, in the absence of TRAIL, DR5 has a strong propensity to self-associate and that this is primarily mediated through interactions of the membrane proximal domains. The association of the free receptors leads to a loss of the threefold symmetry found within the receptor-ligand complex. The simulations suggest that the primary role of TRAIL is to induce threefold-symmetry within the DR5 complex and to constrain the receptor to a specific conformation. The implications of this in terms of the mechanism by which the receptor switches from an inactive to an active state are discussed.

Epitope mapping and structural analysis of an anti‐ErbB2 antibody A21: Molecular basis for tumor inhibitory mechanism

Anti-ErbB2 antibodies targeting distinct epitopes can have different biological functions on cancer cells. A21 prepared by surface epitope masking (SEM) method is a tumor-inhibitory anti-ErbB2 monoclonal antibody. Previously we engineered a single chain chimeric antibody chA21 with potential for therapy of ErbB2-overexpressing tumors. Here, we mapped the A21 epitope on ErbB2 extracellular domain (ECD) by screening a combinatorial phage display peptide library, serial subdomain deletion, and mutagenesis scanning. X-ray crystal structure of the A21 scFv fragment at 2.1 A resolution was also determined. A molecular model of Ag-Ab complex was then constructed based on the crystal structures of the A21 scFv and ErbB2 ECD. Some of biological functions of the A21 mAb and its derivative antibodies including their tumor cell growth inhibition and effects on the expression, internalization, and phosphorylation of ErbB2 receptor were also investigated. The results showed that A21 recognized a conformational epitope comprising a large region mostly from ErbB2 extracellular subdomain I with several surface-exposed residues important for the binding affinity. These data provide unique functional properties of A21 that are quite different from two broadly used anti-ErbB2 mAbs, Herceptin and 2C4. It suggested that the A21 epitope may be another valuable target for designing new anti-ErbB2 therapeutics.

coliSNP database server mapping nsSNPs on protein structures

We have developed coliSNP, a database server (http://yayoi.kansai.jaea.go.jp/colisnp) that maps non-synonymous single nucleotide polymorphisms (nsSNPs) on the three-dimensional (3D) structure of proteins. Once a week, the SNP data from the dbSNP database and the protein structure data from the Protein Data Bank (PDB) are downloaded, and the correspondence of the two data sets is automatically tabulated in the coliSNP database. Given an amino acid sequence, protein name or PDB ID, the server will immediately provide known nsSNP information, including the amino acid mutation caused by the nsSNP, the solvent accessibility, the secondary structure and the flanking residues of the mutated residue in a single page. The position of the nsSNP within the amino acid sequence and on the 3D structure of the protein can also be observed. The database provides key information with which to judge whether an observed nsSNP critically affects protein function and/or stability. As far as we know, this is the only web-based nsSNP database that automatically compiles SNP and protein information in a concise manner.

A hydrophobic loop in acyl-CoA binding protein is functionally important for binding to palmitoyl-coenzyme A: a molecular dynamics study

Acyl-CoA binding protein (ACBP) plays a key role in lipid metabolism, interacting via a partly unknown mechanism with high affinity with long chain fatty acyl-CoAs (LCFA-CoAs). At present there is no study of the microscopic way ligand binding is accomplished. We analyzed this process by molecular dynamics (MDs) simulations. We proposed a computational model of ligand, able to reproduce some evidence from nuclear magnetic resonance (NMR) data, quantitative time resolved fluorometry and X-ray crystallography. We found that a hydrophobic loop, not in the active site, is important for function. Besides, multiple sequence alignment shows hydrophobicity (and not the residues itselves) conservation.

Molecular dynamics simulations from putative transition states of alpha-spectrin SH3 domain

A series of molecular dynamics simulations in explicit solvent were started from nine structural models of the transition state of the SH3 domain of alpha-spectrin, which were generated by Lindorff-Larsen et al. (Nat Struct Mol Biol 2004;11:443-449) using molecular dynamics simulations in which experimental Phi - values were incorporated as restraints. Two of the nine models were simulated 10 times for 200 ns and the remaining models simulated two times for 200 ns. Complete folding was observed in one case, while in the other simulations partial folding or unfolding events were observed, which were characterized by a regularization of elements of secondary structure. These results are consistent with recent experimental evidence that the folding of SH3 domains involves low populated intermediate states.

Characterization of the protein unfolding processes induced by urea and temperature

Correct folding is critical for the biological activities of proteins. As a contribution to a better understanding of the protein (un)folding problem, we studied the effect of temperature and of urea on peptostreptococcal Protein L destructuration. We performed standard molecular dynamics simulations at 300 K, 350 K, 400 K, and 480 K, both in 10 M urea and in water. Protein L followed at least two alternative unfolding pathways. Urea caused the loss of secondary structure acting preferentially on the beta-sheets, while leaving the alpha-helices almost intact; on the contrary, high temperature preserved the beta-sheets and led to a complete loss of the alpha-helices. These data suggest that urea and high temperature act through different unfolding mechanisms, and protein secondary motives reveal a differential sensitivity to various denaturant treatments. As further validation of our results, replica-exchange molecular dynamics simulations of the temperature-induced unfolding process in the presence of urea were performed. This set of simulations allowed us to compute the thermodynamical parameters of the process and confirmed that, in the configurational space of Protein L unfolding, both of the above pathways are accessible, although to a different relative extent.

Joint neighbors approximation of macromolecular solvent accessible surface area

A new method for approximate analytical calculations of solvent accessible surface area (SASA) for arbitrary molecules and their gradients with respect to their atomic coordinates was developed. This method is based on the recursive procedure of pairwise joining of neighboring atoms. Unlike other available methods of approximate SASA calculations, the method has no empirical parameters, and therefore can be used with comparable accuracy in calculations of SASA in folded and unfolded conformations of macromolecules of any chemical nature. As shown by tests with globular proteins in folded conformations, average errors in absolute atomic surface area is around 1 A2, while for unfolded protein conformations it varies from 1.65 to 1.87 A2. Computational times of the method are comparable with those by GETAREA, one of the fastest exact analytical methods available today.

Improved evaluation of liquid densities using van der Waals molecular models

A new approach is presented to estimate molar volumes and densities of liquids in ambient conditions from van der Waals models, taking advantage of the correlation between the intermolecular volume and the atomic contributions to the molecular surface area. Using this approach, the role of hydrogen bonds can be quantified. The densities obtained prove remarkably close to the values derived from the ACD group contribution method. However, the present approach requires much less empirical parameters and may be applied to arbitrary H-C-N-O-F-S-Cl-Br compounds. Moreover, it is more reliable than earlier empirical procedures, including quantitative structure property relationships (QSPR). While it does not correct the deficiencies of group contribution methods, it provides a natural approach to introduce temperature effects.

Aqueous urea solutions: structure, energetics, and urea aggregation

Urea is ubiquitously used as a protein denaturant. To study the structure and energetics of aqueous urea solutions, we have carried out molecular dynamics simulations for a wide range of urea concentrations and temperatures. The hydrogen bonds between urea and water were found to be significantly weaker than those between water molecules, which drives urea self-aggregation due to the hydrophobic effect. From the reduction of the water exposed urea surface area, urea was found to exhibit an aggregation degree of ca. 20% at concentrations commonly used for protein denaturation. Structurally, three distinct urea pair conformations were identified and their populations were analyzed by translational and orientational pair distribution functions. Furthermore, urea was found to strengthen water structure in terms of hydrogen bond energies and population of solvation shells. Our findings are consistent with a direct interaction between urea and the protein as the main driving force for protein denaturation. As an additional, more indirect effect, urea was found to enhance water structure, which would suggest a weakening of the hydrophobic effect.

How sensitive are nanosecond molecular dynamics simulations of proteins to changes in the force field?

The sensitivity of molecular dynamics simulations to variations in the force field has been examined in relation to a set of 36 structures corresponding to 31 proteins simulated by using different versions of the GROMOS force field. The three parameter sets used (43a1, 53a5, and 53a6) differ significantly in regard to the nonbonded parameters for polar functional groups and their ability to reproduce the correct solvation and partitioning behavior of small molecular analogues of the amino acid side chains. Despite the differences in the force field parameters no major differences could be detected in a wide range of structural properties such as the root-mean-square deviation from the experimental structure, radii of gyration, solvent accessible surface, secondary structure, or hydrogen bond propensities on a 5 to 10 ns time scale. The small differences that were observed correlated primarily with the presence of charged residues as opposed to residues that differed most between the parameter sets. The work highlights the variation that can be observed in nanosecond simulations of protein systems and implications of this for force field validation, as well as for the analysis of protein simulations in general.

Tailoring cutinase activity towards polyethylene terephthalate and polyamide 6,6 fibers

Cutinase from Fusarium solani pisi was genetically modified near the active site, by site-directed mutagenesis, to enhance its activity towards polyethylene terephthalate (PET) and polyamide 6,6 (PA 6,6) fibers. The mutations L81A, N84A, L182A, V184A and L189A were done to enlarge the active site in order to better fit a larger polymer chain. Modeling studies have shown enhanced free energy stabilization of model substrate tetrahedral intermediate (TI) bound at the enzyme active site for all mutants, for both model polymers. L81A and L182A showed an activity increase of four- and five-fold, respectively, when compared with the wild type, for PET fibers. L182A showed the one- and two-fold higher ability to biodegrade aliphatic polyamide substrates. Further studies in aliphatic polyesters seem to indicate that cutinase has higher ability to recognize aliphatic substrates.

Protein Structures under Electrospray Conditions

During electrospray ionization (ESI), proteins are transferred from solution into vacuum, a process that influences the conformation of the protein. Exactly how much the conformation changes due to the dehydration process, and in what way, is difficult to determine experimentally. The aim of this study is therefore to monitor what happens to protein structures as the surrounding waters gradually evaporate, using computer simulations of the transition of proteins from water to vacuum. Five different proteins have been simulated with water shells of varying thickness, enabling us to mimic the entire dehydration process. We find that all protein structures are affected, at least to some extent, by the transfer but that the major features are preserved. A water shell with a thickness of roughly two molecules is enough to emulate bulk water and to largely maintain the solution phase structure. The conformations obtained in vacuum are quite similar and make up an ensemble which differs from the structure obtained by experimental means, and from the solution phase structure as found in simulations. Dehydration forces the protein to make more intramolecular hydrogen bonds, at the expense of exposing more hydrophobic area (to vacuum). Native hydrogen bonds usually persist in vacuum, yielding an easy route to refolding upon rehydration. The findings presented here are promising for future bio-imaging experiments with X-ray free electron lasers, and they strongly support the validity of mass spectrometry experiments for studies of intra- and intermolecular interactions.

Influence of the Hofmeister Anions on Protein Stability As Studied by Thermal Denaturation and Chemical Shift Perturbation†

The influence of external cosolutes on the thermal stability of the B1 domain of protein L (ProtL) has been studied by circular dichroism, fluorescence spectroscopy, and differential scanning calorimetry. The thermal denaturation midpoint is effectively modulated by the addition of a suite of anions and follows the Hofmeister series. The maximum increase in thermostability (corresponding to 14 degrees C) was observed in the presence of 1 M sodium sulfate. After conversion of the experimental data into the change in the virial coefficient, a mechanistic model was used to estimate the relative contributions from excluded volume and preferential anion solvation for each anion. As expected, the excluded volume term stabilizes the native conformation of ProtL for all the cosolutes, but opposite effects on protein stability arise from the anion's solvation depending on their tendency to interact with or to become excluded from the protein surface. This behavior is in agreement with the results of independent NMR experiments: the anions that strongly interact with the protein surface produce significant perturbations in the amide protein chemical shift (delta d23(HN)). A correlation obtained between delta d23(HN) and the temperature coefficients for the different amide protons provides qualitative information about the structural determinants for the interaction between the protein surface and the cosolute.

Temperature and structural changes of water clusters in vacuum due to evaporation

This paper presents a study on evaporation of pure water clusters. Molecular dynamics simulations between 20 ns and 3 micros of clusters ranging from 125 to 4096 molecules in vacuum were performed. Three different models (SPC, TIP4P, and TIP5P) were used to simulate water, starting at temperatures of 250, 275, and 300 K. We monitored the temperature, the number of hydrogen bonds, the tetrahedral order, the evaporation, the radial distribution functions, and the diffusion coefficients. The three models behave very similarly as far as temperature and evaporation are concerned. Clusters starting at a higher temperature show a higher initial evaporation rate and therefore reach the point where evaporation stop (around 240 K) sooner. The radius of the clusters is decreased by 0.16-0.22 nm after 0.5 micros (larger clusters tend to decrease their radius slightly more), which corresponds to around one evaporated molecule per nm(2). The cluster temperature seems to converge towards 215 K independent of cluster size, when starting at 275 K. We observe only small structural changes, but the clusters modeled by TIP5P show a larger percentage of molecules with low diffusion coefficient as t-->infinity, than those using the two other water models. TIP4P seems to be more structured and more hydrogen bonds are formed than in the other models as the temperature falls. The cooling rates are in good agreement with experimental results, and evaporation rates agree well with a phenomenological expression based on experimental observations.

Modeling H3 histone N-terminal tail and linker DNA interactions

Molecular dynamics computer simulations were performed for the 25-residue N-terminal tail of the H3 histone protein in the proximity of a DNA segment of 10 base pairs (bp), representing a model for the linker DNA in chromatin. Several least biased configurations were used as initial configurations. The secondary structure content of the protein was increased by the presence of DNA close to it, but the locations of the secondary motifs were different for different initial orientations of the DNA grooves with respect to the protein. As a common feature to all simulations, the electrostatic attraction between negatively charged DNA and positively charged protein was screened by the water solvent and counterbalanced by the intrinsic compaction of the protein due to hydrophobic effects. The protein secondary structure limited the covering of DNA by the protein to 4-5 bp. The degree of compaction and charge density of the bound protein suggests a possible role of H3 tail in a nonspecific bending and plasticity of the linker DNA when the protein is located in the crowded dense chromatin.

Molecular dynamics analysis of HIV-1 matrix protein: clarifying differences between crystallographic and solution structures

One of the main structural features of the mature HIV-1 virion is the matrix protein (p17). This partially globular protein presents four helixes centrally organized and a fifth one, H5, projecting away from the packed bundle of helixes. Comparison between solution and crystallographic data of p17 indicates a 6 A displacement of a short 3(10) helix and a partial unfolding of H5 in solution related to crystal. While the behavior of the 3(10) helix has been previously addressed to virion assembly, the relevance and origin of H5 partial unfolding is possibly related to the contacts between p17 and other viral elements, such as p24. In this context, we present a 40 ns conformational sampling of monomeric p17 using molecular dynamics simulations. The performed simulations presented a progressive conversion of the p17 crystallographic structure to the NMR conformation, suggesting that the biological form of this protein may have its C-terminal portion partially unfolded.

Conformational polymorphism of the PrP106-126 peptide in different environments: a molecular dynamics study

Extensive molecular dynamic simulations (approximately 240 ns) have been used to investigate the conformational behavior of PrP106-126 prion peptide in four different environments (water, dimethyl sulfoxide, hexane, and trifluoroethanol) and under both neutral and acidic conditions. The conformational polymorphism of PrP106-126 in solution observed in the simulations supports the role of this fragment in the structural transition of the native to the abnormal form of prion protein in response to changes in the local environmental conditions. The peptide in solution is primarily unstructured. The simulations show an increased presence of helical structure in an apolar solvent, in agreement with the results from circular dichroism spectroscopy. In water solution, beta-sheet elements were observed between residues 108-112 and either residues 115-121 or 121-126. An alpha-beta transition was observed under neutral conditions. In DMSO, the peptide adopted an extended conformation, in agreement with nuclear magnetic resonance experiments.

Molecular dynamics study of a calmodulin-like protein with an IQ peptide: spontaneous refolding of the protein around the peptide

The Calmodulin (CaM) is a small (16.7 kDa), highly acidic protein that is crucial to all eukaryotes by serving as a prototypical calcium sensor. In the present study, we investigated, through molecular dynamics simulations, the dynamics of a complex between the Mlc1p protein, which is a CaM-like protein, and the IQ4 peptide. This protein-peptide interaction is of high importance because IQ motifs are widely distributed among different kinds of CaM-binding proteins. The Mlc1p-IQ4 complex, which had been resolved by crystallography to 2.1 A, confers to a Ca(+2)-independent stable structure. During the simulations, the complex undergoes a complicated modulation process, which involves bending of the angles between the alpha-helices of the protein, breaking of the alpha-helical structure of the IQ4 peptide into two sections, and formation of new contact points between the protein and the peptide. The dynamics of the process consist of fast sub picosecond events and much slower ones that take a few nanoseconds to completion. Our study expands the information embedded in the crystal structure of the Mlc1p-IQ4 complex by describing its dynamic behavior as it evolves from the crystal structure to a form stable in solution. The article shows that careful application of molecular dynamics simulations can be used for extending the structural information presented by the crystal structure, thereby revealing the dynamic configuration of the protein in its physiological environment.

A molecular dynamics study and free energy analysis of complexes between the Mlc1p protein and two IQ motif peptides

The Mlc1p protein from the budding yeast Saccharomyces cerevisiae is a Calmodulin-like protein, which interacts with IQ-motif peptides located at the yeast's myosin neck. In this study, we report a molecular dynamics study of the Mlc1p-IQ2 protein-peptide complex, starting with its crystal structure, and investigate its dynamics in an aqueous solution. The results are compared with those obtained by a previous study, where we followed the solution structure of the Mlc1p-IQ4 protein-peptide complex by molecular dynamics simulations. After the simulations, we performed an interaction free-energy analysis using the molecular mechanics Poisson-Boltzmann surface area approach. Based on the dynamics of the Mlc1p-IQ protein-peptide complexes, the structure of the light-chain-binding domain of myosin V from the yeast S. cerevisiae is discussed.

Molecular simulation of the binding of nerve growth factor peptide mimics to the receptor tyrosine kinase A

Nerve growth factor (NGF) mimics play an important role for therapies that target the receptor tyrosine kinase A (trkA). The N-terminal fragment of the NGF (N-term@NGF) was previously demonstrated to be an important determinant for affinity and specificity in the binding to trkA. Here we use a variety of computational tools (contact surface analysis and free energy predictions) to identify residues playing a key role for the binding to the receptor. Molecular dynamics simulations are then used to investigate the stability of complexes between trkA and peptides mimicking N-term@NGF. Steered molecular dynamics calculations are finally performed to investigate the process of detaching the peptide from the receptor. Three disruptive events are observed, the first involving the breaking of all intermolecular interactions except two salt bridges, which break subsequently.

Lifting a wet glass from a table: a microscopic picture

Why is it so hard to lift a wet glass from a table? Is it easier when there is whiskey between the glass and the table? Macroscopically, the picture is quite simple: two surfaces have to be disrupted that are connected indirectly through hydrogen bonds and/or van der Waals forces. In the beginning, a surface has to be created leading to surface tension, and after that a liquid bridge has to be broken. Here we study the phenomenon at the microscopic level using molecular dynamics simulations. The effective force between two quartz plates is measured at different distances and with different alcohol/water mixtures between them. This allows us to compute the total work necessary to "lift the glass from the table". Different aspects of the process, such as clustering and liquid ordering are discussed. We compare the structure of the liquid/glass interface to that of a liquid/vapor interface, for which we present simulation results, like surface tension, as well. On the basis of the simulations, we are able to provide a detailed description of the energetics during the separation process as a function of alcohol concentration. It is shown that there is a net entropy loss upon separating two plates with water or a 10% MeOH solution between them, whereas for higher alcohol concentrations, there is net entropy gain. These findings increase our understanding of the properties of colloid suspensions which is important for process technology.

Structure and assembly of P-pili: a protruding hinge region used for assembly of a bacterial adhesion filament

High-resolution structures of macromolecular complexes offer unparalleled insight into the workings of biological systems and hence the interplay of these systems in health and disease. We have adopted a multifaceted approach to understanding the pathogenically important structure of P-pili, the class I adhesion pili from pyelonephritic Escherichia coli. Our approach combines electron cryomicroscopy, site-directed mutagenesis, homology modeling, and energy calculations, resulting in a high-resolution model of PapA, the major structural element of these pili. Fitting of the modeled PapA subunit into the electron cryomicroscopy data provides a detailed view of these pilins within the supramolecular architecture of the pilus filament. A structural hinge in the N-terminal region of the subunit is located at the site of a newly resolved electron density that protrudes from the P-pilus surface. The structural flexibility provided by this hinge is necessary for assembly of P-pili, illustrating one solution to construction of large macromolecular complexes from small repeating units. These data support our hypothesis that domain-swapped pilin subunits transit the outer cell membrane vertically and rotate about the hinge for final positioning into the pilus filament. Our data confirm and supply a structural basis for much previous genetic, biochemical, and structural data. This model of the P-pilus filament provides an insight into the mechanism of assembly of a macromolecular complex essential for initiation of kidney infection by these bacteria.

Optimized parallel tempering simulations of proteins

We apply a recently developed adaptive algorithm that systematically improves the efficiency of parallel tempering or replica exchange methods in the numerical simulation of small proteins. Feedback iterations allow us to identify an optimal set of temperatures/replicas which are found to concentrate at the bottlenecks of the simulations. A measure of convergence for the equilibration of the parallel tempering algorithm is discussed. We test our algorithm by simulating the 36-residue villin headpiece subdomain HP-36 where we find a lowest-energy configuration with a root-mean-square deviation of less than 4 A to the experimentally determined structure.

PSI-DOCK: towards highly efficient and accurate flexible ligand docking

We have developed a new docking method, Pose-Sensitive Inclined (PSI)-DOCK, for flexible ligand docking. An improved SCORE function has been developed and used in PSI-DOCK for binding free energy evaluation. The improved SCORE function was able to reproduce the absolute binding free energies of a training set of 200 protein-ligand complexes with a correlation coefficient of 0.788 and a standard error of 8.13 kJ/mol. For ligand binding pose exploration, a unique searching strategy was designed in PSI-DOCK. In the first step, a tabu-enhanced genetic algorithm with a rapid shape-complementary scoring function is used to roughly explore and store potential binding poses of the ligand. Then, these predicted binding poses are optimized and compete against each other by using a genetic algorithm with the accurate SCORE function to determine the binding pose with the lowest docking energy. The PSI-DOCK 1.0 program is highly efficient in identifying the experimental binding pose. For a test dataset of 194 complexes, PSI-DOCK 1.0 achieved a 67% success rate (RMSD < 2.0 A) for only one run and a 74% success rate for 10 runs. PSI-DOCK can also predict the docking binding free energy with high accuracy. For a test set of 64 complexes, the correlation between the experimentally observed binding free energies and the docking binding free energies for 64 complexes is r = 0.777 with a standard deviation of 7.96 kJ/mol. Moreover, compared with other docking methods, PSI-DOCK 1.0 is extremely easy to use and requires minimum docking preparations. There is no requirement for the users to add hydrogen atoms to proteins because all protein hydrogen atoms and the flexibility of the terminal protein atoms are intrinsically taken into account in PSI-DOCK. There is also no requirement for the users to calculate partial atomic charges because PSI-DOCK does not calculate an electrostatic energy term. These features are not only convenient for the users but also help to avoid the influence of different preparation methods.

A molecular dynamics study of the effect of Ca2+ removal on calmodulin structure

Calmodulin is a small (148 residues), ubiquitous, highly-conserved Ca(2+) binding protein serving as a modulator of many calcium-dependent processes. In this study, we followed, by means of molecular dynamics, the structural stability of the protein when one of its four bound Ca(2+) ions is removed, and compared it to a simulation of the fully Ca(2+) bound protein. We found that the removal of a single Ca(2+) ion from the N-lobe of the protein, which has a lower affinity for the ion, is sufficient to initiate a considerable structural rearrangement. Although the overall structure of the fully 4 Ca(2+) bound protein remained intact in the extended conformation, the Ca(2+)-removed protein changed its conformation into a compact state. The observation that the 3 Ca(2+) loaded protein assumes a compacted solution state is in accord with experimental observation that the NSCP protein, which binds only three Ca(2+) ions, is natively in a compact state. Examination of the folding dynamics reveals a cooperation between the C-lobe, N-lobe, and the interdomain helix that enable the conformation change. The forces driving this conformational change are discussed.

cAMP Modulation of the cytoplasmic domain in the HCN2 channel investigated by molecular simulations

The hyperpolarization-activated cyclic nucleotide-modulated (HCN) cation channels are opened by membrane hyperpolarization, while their activation is modulated by the binding of cyclic adenosine monophosphate (cAMP) in the cytoplasm. Here we investigate the molecular basis of cAMP channel modulation by performing molecular dynamics simulations of a segment comprising the C-linker and the cyclic nucleotide binding domain (CNBD) in the presence and absence of cAMP, based on the available crystal structure of HCN2 from mouse. In presence of cAMP, the protein undergoes an oscillation of the quaternary structure on the order of 10 ns, not observed in the apoprotein. In contrast, the absence of ligand causes conformational rearrangements within the CNBDs, driving these domains to a more flexible state, similar to that described in CNBDs of other proteins. This increased flexibility causes a rather disordered movement of the CNBDs, resulting in an inhibitory effect on the channel. We propose that the cAMP-triggered large-scale oscillation plays an important role for the channel's function, being coupled to a motion of the C-linker which, in turn, modulates the gating of the channel.

Effect of pressure on the conformation of proteins. A molecular dynamics simulation of lysozyme

The effect of pressure on the structure and mobility of lysozyme was studied by molecular dynamics computer simulation at 1 and 3 kbar (1 atm = 1.01325 bar = 101.325 kPa). The results have good agreement with the available experimental data, allowing the analysis of other features of the effect of pressure on the protein solution. The studies of mobility show that although the general mobility is restricted under pressure this is not true for some particular residues. From the analysis of secondary structure along the trajectories it is observed that the conformation under pressure is more stable, suggesting that pressure acts as a 'conformer selector' on the protein. The difference in solvent-accessed surface (SAS) with pressure shows a clear inversion of the hydrophilic/hydrophobic SAS ratio, which consequently shows that the hydrophobic interaction is considerably weaker under high hydrostatic pressure conditions.

The solvation interface is a determining factor in peptide conformational preferences

The 21 residue polyalanine-based F(s) peptide was studied using thousands of long, explicit solvent, atomistic molecular dynamics simulations that reached equilibrium at the ensemble level. Peptide conformational preference as a function of hydrophobicity was examined using a spectrum of explicit solvent models, and the peptide length-dependence of the hydrophilic and hydrophobic components of solvent-accessible surface area for several ideal conformational types was considered. Our results demonstrate how the character of the solvation interface induces several conformational preferences, including a decrease in mean helical content with increased hydrophilicity, which occurs predominantly through reduced nucleation tendency and, to a lesser extent, destabilization of helical propagation. Interestingly, an opposing effect occurs through increased propensity for 3(10)-helix conformations, as well as increased polyproline structure. Our observations provide a framework for understanding previous reports of conformational preferences in polyalanine-based peptides including (i) terminal 3(10)-helix prominence, (ii) low pi-helix propensity, (iii) increased polyproline conformations in short and unfolded peptides, and (iv) membrane helix stability in the presence and absence of water. These observations provide physical insight into the role of water in peptide conformational equilibria at the atomic level, and expand our view of the complexity of even the most "simple" of biopolymers. Whereas previous studies have focused predominantly on hydrophobic effects with respect to tertiary structure, this work highlights the need for consideration of such effects at the secondary structural level.

A salt-bridge motif involved in ligand binding and large-scale domain motions of the maltose-binding protein

The uptake of nutrients is essential for the survival of bacterial cells. Many specialized systems have evolved, such as the maltose-dependent ABC transport system that transfers oligosaccharides through the cytoplasmic membrane. The maltose/maltodextrin-binding protein (MBP) serves as an initial high-affinity binding component in the periplasm that delivers the bound sugar into the cognate ABC transporter MalFGK(2). We have investigated the domain motions induced by the binding of the ligand maltotriose into the binding cleft using molecular dynamics simulations. We find that MBP is predominantly in the open state without ligand and in the closed state with ligand bound. Oligosaccharide binding induces a closure motion (30.0 degrees rotation), whereas ligand removal leads to domain opening (32.6 degrees rotation) around a well-defined hinge affecting key areas relevant for chemotaxis and transport. Our simulations suggest that a "hook-and-eye" motif is involved in the binding. A salt bridge between Glu-111 and Lys-15 forms that effectively locks the protein-ligand complex in a semiclosed conformation inhibiting any further opening and promoting complete closure. This previously unrecognized feature seems to secure the ligand in the binding site and keeps MBP in the closed conformation and suggests a role in the initial steps of substrate transport.

Reorganization and conformational changes in the reduction of tetraheme cytochromes

Molecular dynamics simulation (MD) constitutes an alternative to time-consuming experiments for studying conformational changes. We apply MD on a redox system where experimental information exists for the fully oxidized and fully reduced states: tetraheme cytochrome c3. Instead of doing one simulation for each state, we apply 10 4-ns replicas for both states, which provides robust statistics to characterize the redox changes. Besides these long simulations, we perform 120 short ones (50 ps), where an equilibrated oxidized state is perturbed to a reduced state. This allows the application of a nonequilibrium method, the subtraction technique, which makes it possible to characterize the different timescales of conformational changes. Reduction induces conformational changes in the N-terminus and on the loops spanning residues 36-42 and 88-93, which correlate very well with experiments, demonstrating the applicability of this methodology. We also analyze the effect of reduction on hydrogen bonds, solvent accessible surface and bound water, the changes being found to involve the hemes and propionate groups. Redox-induced protonation is also investigated, by protonating the propionates D from hemes I and IV. Although this change in the former does not have major conformational consequences, it induces in the latter conformational changes beyond the ones obtained with reduction.

A molecular dynamics study of the formation, stability, and oligomerization state of two designed coiled coils: possibilities and limitations

The formation, relative stability, and possible stoichiometries of two (self-)complementary peptide sequences (B and E) designed to form either a parallel homodimeric (B + B) or an antiparallel heterodimeric (B + E) coiled coil have been investigated. Peptide B shows a characteristic coiled coil pattern in circular dichroism spectra at pH 7.4, whereas peptide E is apparently random coiled under these conditions. The peptides are complementary to each other, with peptide E forming a coiled coil when mixed with peptide B. Molecular dynamics simulations show that combinations of B + B and B + E readily form a dimeric coiled coil, whereas E + E does not fall in line with the experimental data. However, the simulations strongly suggest the preferred orientation of the helices in the homodimeric coiled coil is antiparallel, with interactions at the interface quite different to that of the idealized model. In addition, molecular dynamics simulations suggest equilibrium between dimers, trimers, and tetramers of alpha-helices for peptide B.

Evidence of a double-lid movement in Pseudomonas aeruginosa lipase: insights from molecular dynamics simulations

Pseudomonas aeruginosa lipase is a 29-kDa protein that, following the determination of its crystal structure, was postulated to have a lid that stretched between residues 125 and 148. In this paper, using molecular dynamics simulations, we propose that there exists, in addition to the above-mentioned lid, a novel second lid in this lipase. We further show that the second lid, covering residues 210–222, acts as a triggering lid for the movement of the first. We also investigate the role of hydrophobicity in the movement of the lids and show that two residues, Phe214 and Ala217, play important roles in lid movement. To our knowledge, this is the first time that a double-lid movement of the type described in our manuscript has been presented to the scientific community. This work also elucidates the interplay of hydrophobic interactions in the dynamics, and hence the function, of an enzyme.

Lipases hydrolyse long-chain fatty acid esters at water-oil interfaces through the mechanism of interfacial activation mediated by the movement of a lid subdomain that covers the active site. Studying lid movement is an area of active research in the field of protein dynamics. The lipase from Pseudomonas aeruginosa is a 29-kDa protein that was previously crystallized in the open conformation, and as expected, an approximately 20-residue lid subdomain was identified. In the present study, the authors report extensive molecular dynamics simulations of the P. aeruginosa lipase. They show that this protein has two lids covering the substrate-binding pocket. The first lid is the one proposed from the known crystal structure. The second lid, a much shorter one, lies over the binding pocket facing the first lid. Furthermore, using position-restrained simulations, these authors show that movement of the second lid may actually be a trigger for the movement of the first, and that this triggering action is driven by hydrophobic contacts between the two lids. This computational study paves a way for experimentalists to study the structure and dynamics of this protein in greater detail in order to understand coupled subdomain movements in a comprehensive fashion.

On the Use of Different Dielectric Constants for Computing Individual and Pairwise Terms in Poisson−Boltzmann Studies of Protein Ionization Equilibrium

Poisson-Boltzmann (PB) models are a fast and common tool for studying electrostatic processes in proteins, particularly their ionization equilibrium (protonation and/or reduction), often yielding quite good results when compared with more detailed models. Yet, they are conceptually very simple and necessarily approximate, their empirical character being most evident when it comes to the choice of the dielectric constant assigned to the protein region. The present study analyzes several factors affecting the ability of PB-based methods to model protein ionization equilibrium. We give particular attention to a suggestion made by Warshel and co-workers (e.g., Sham et al. J. Phys. Chem. B 1997, 101, 4458) of using different protein dielectric constants for computing the individual (site) and the pairwise (site-site) terms of the ionization free energies. Our prediction of pK(a) values for several proteins indicates that no advantage is obtained by such a procedure, even for sites that are buried and/or display large pK(a) shifts relative to the solution values. In particular, the present methodology gives the best predictions using a dielectric constant around 20, for shifted/buried and nonshifted/exposed sites alike. The similarities and differences between the PB model and Warshel's PDLD/S model are discussed, as well as the reasons behind their apparently discrepant results. The present PB model is shown to predict also good reduction potentials in redox proteins.

Generalized-ensemble simulations of the human parathyroid hormone fragment PTH(1-34)

A generalized-ensemble technique, multicanonical sampling, is used to study the folding of a 34-residue human parathyroid hormone fragment. An all-atom model of the peptide is employed and the protein-solvent interactions are approximated by an implicit solvent. Our results demonstrate that generalized-ensemble simulations are well suited to sample low-energy structures of such large polypeptides. Configurations with a root-mean-square deviation to the crystal structure of less than 1 A are found. Finally, we discuss limitations of our implicit solvent model.

The exchangeability of amino acids in proteins

The comparative analysis of protein sequences depends crucially on measures of amino acid similarity or distance. Many such measures exist, yet it is not known how well these measures reflect the operational exchangeability of amino acids in proteins, since most are derived by methods that confound a variety of effects, including effects of mutation. In pursuit of a pure measure of exchangeability, we present (1) a compilation of data on the effects of 9671 amino acid exchanges engineered and assayed in a set of 12 proteins; (2) a statistical procedure to combine results from diverse assays of exchange effects; (3) a matrix of "experimental exchangeability" values EX(ij) derived from applying this procedure to the compiled data; and (4) a set of three tests designed to evaluate the power of an exchangeability measure to (i) predict the effects of amino acid exchanges in the laboratory, (ii) account for the disease-causing potential of missense mutations in the human population, and (iii) model the probability of fixation of missense mutations in evolution. EX not only captures useful information on exchangeability while remaining free of other effects, but also outperforms all measures tested except for the best-performing alignment scoring matrix, which is comparable in performance.

Binding behaviour and conformational properties of globular proteins in the presence of immobilised non-polar ligands used in reversed-phase liquid chromatography

The thermodynamic and extra-thermodynamic dependencies of five types of cytochrome c in water-acetonitrile mixtures of different composition in the presence of immobilised n-octyl ligands as a function of temperature from 278 K to 338 K have been investigated. The corresponding enthalpic, entropic and heat capacity parameters, deltaHdegrees assoc, deltaS degrees assoc and delta C degrees p, have been evaluated from the observed non-linear Van't Hoff plots of these globular proteins in these heterogeneous systems. The relationships between the free energy dependencies, various molecular parameters and extra-thermodynamic dependencies (empirical correlations) of these protein-non-polar ligand interactions have also been examined. Thus, the involvement of enthalpy-entropy compensation effects has been documented for the binding of these cytochrome cs to solvated n-octyl ligands. Moreover, the results confirm that this experimental approach permits changes in molecular surface area due to the unfolding of these proteins on association with non-polar ligands as a function of temperature to be correlated with other biophysical properties. This study thus provides a general procedure whereby the corresponding free energy dependencies of globular proteins on association with solvated non-polar ligands in heterogeneous two-phase systems can be quantitatively evaluated in terms of fundamental molecular parameters.

Distribution-Based Descriptors of the Molecular Shape

A rational design of economically cost-effective chemical libraries as well as successful data mining during a process of drug discovery employs a vast array of the molecular descriptors. Despite the huge importance of this area of the research there is still a need for the further development of the simple, intuitive, easily calculable, specific and size-invariant parameters of the molecular shape. Here we present ab initio calculation of the molecular volumes and expectations for the molecular areas of projection. These molecular size parameters were used as a basis to define a group of novel descriptors of the molecular shape. A set of molecular descriptors was developed: ovality, roughness, size-corrected parameters, and the parameters derived from the higher central momenta of the distributions of the size descriptors-skewness and kurtosis. The rationale for the construction of the descriptors was first to calculate the descriptors of the molecular size along many directions in the space and second to use the statistical parameters of the distribution of those descriptors as the shape descriptors. The size descriptors well suited for the above purpose and discussed in this paper are generalized molecular radii and the molecular areas of projection. Molecular volume and projection area-derived descriptors were calculated and their applicability as shape descriptors was illustrated using exploratory methods (factor analysis and hierarchical cluster analysis). The shape descriptors appear to be promising in their ability to discriminate and classify the molecular shapes e.g. spheroids, disklike, rodlike, starlike, crosslike, anglelike, etc.

Molecular Dynamics Simulations of Helix-Forming, Amine-Functionalized m-Poly(phenyleneethynylene)s

We present here the results of all-atom and united-atom molecular dynamics (MD) simulations that were used to examine the folding behavior of an amine-functionalized m-poly(phenyleneethynylene) (m-PPE) oligomer in aqueous environment. The parallelized GROMACS MD simulation code and OPLS force field were used for multiple MD simulations of m-PPE oligomers containing 24 phenyl rings in extended, coiled and helix conformations separately in water to determine the minimum energy conformation of the oligomer in aqueous solvent and what interactions are most important in determining this structure. Simulation results showed that the helix is the preferred minimum energy conformation of a single oligomer in water and that Lennard-Jones interactions are the dominant forces for the stabilization of the helix. In addition, these solvophobic interactions are strong enough to maintain the helix conformation at temperatures up to 523 K.

Modeling the backbone dynamics of reduced and oxidized solvated rat microsomal cytochrome b5

In this article, a description of the statistics and dynamics of cytochrome b(5) in both reduced and oxidized forms is given. Results of molecular dynamics computer simulations in the explicit solvent have been combined with mode-coupling diffusion models including and neglecting the molecule-solvent correlations. R(1) and R(1 rho) nuclear magnetic relaxation parameters of (15)N in the protein backbone have been calculated and compared with experiments. Slight changes in charge density in the heme upon oxidation produces a cascade of changes in charge distributions from heme propionates up to charged residues approximately 1.5 nm from Fe. These changes in charge distributions modify the molecular surface and the water shell surrounding the protein. The statistical changes upon oxidation can be included in diffusive models that physically explain the upper and lower limits of R(1 rho) relaxation parameters at high off-resonance fields.

Prediction of charge-induced molecular alignment of biomolecules dissolved in dilute liquid-crystalline phases

Alignment of macromolecules in nearly neutral aqueous lyotropic liquid-crystalline media such as bicelles, commonly used in macromolecular NMR studies, can be predicted accurately by a steric obstruction model (Zweckstetter and Bax, 2000). A simple extension of this model is described that results in improved predictions for both the alignment orientation and magnitude of protein and DNA solutes in charged nematic media, such as the widely used medium of filamentous phage Pf1. The extended model approximates the electrostatic interaction between a solute and an ordered phage particle as that between the solute's surface charges and the electric field of the phage. The model is evaluated for four different proteins and a DNA oligomer. Results indicate that alignment in charged nematic media is a function not only of the solute's shape, but also of its electric multipole moments of net charge, dipole, and quadrupole. The relative importance of these terms varies greatly from one macromolecule to another, and evaluation of the experimental data indicates that these terms scale differently with ionic strength. For several of the proteins, the calculated alignment is sensitive to the precise position of the charged groups on the protein surface. This suggests that NMR alignment measurements can potentially be used to probe protein electrostatics. Inclusion of electrostatic interactions in addition to steric effects makes the extended model applicable to all liquid crystals used in biological NMR to date.

A new analytical method for computing solvent-accessible surface area of macromolecules and its gradients

In the calculation of thermodynamic properties and three-dimensional structures of macromolecules, such as proteins, it is important to have an efficient algorithm for computing the solvent-accessible surface area of macromolecules. Here, we propose a new analytical method for this purpose. In the proposed algorithm we consider the transformation that maps the spherical circles formed by intersection of the atomic surfaces in three-dimensional space onto the circles on a two-dimensional plane, and the problem of computing the solvent-accessible surface area is reduced to the problem of computing the corresponding curve integrals on the plane. This allows to consider only the integrals along the circular trajectories on the plane. The algorithm is suitable for parallelization. Testings on many proteins as well as the comparison to the other analogous algorithms have shown that our method is accurate and efficient.