A general method to quantify quasi-equivalence in icosahedral viruses

A quantitative, atom-based, method is described for comparing protein subunit interfaces in icosahedral virus capsids with quasi-equivalent surface lattices. An integrated, normalized value (between 0 and 1) based on equivalent residue contacts (Q-score) is computed for every pair of subunit interactions and scores that are significantly above zero readily identify interfaces that are quasi-equivalent to each other. The method was applied to all quasi-equivalent capsid structures (T=3, 4, 7 and 13) in the Protein Data Bank and the Q-scores were interpreted in terms of their structural underpinnings. The analysis allowed classification of T=3 structures into three groups with architectures that resemble different polyhedra with icosahedral symmetry. The preference of subunits to form dimers in the T=4 human Hepatitis B virus capsid (HBV) was clearly reflected in high Q-scores of quasi-equivalent dimers. Interesting differences between the classical T=7 capsid and polyoma-like capsids were also identified. Application of the method to the outer-shell of the T=13 Blue tongue virus core (BTVC) highlighted the modest distortion between the interfaces of the general trimers and the strict trimers of VP7 subunits. Furthermore, the method identified the quasi 2-fold symmetry in the inner capsids of the BTV and reovirus cores. The results show that the Q-scores of various quasi-symmetries represent a "fingerprint" for a particular virus capsid architecture allowing particle classification into groups based on their underlying structural and geometric features.

Binding preferences of hydroxamate inhibitors of the matrix metalloproteinase human fibroblast collagenase

In this paper we report molecular dynamics (MD) and free energy perturbation (FEP) studies carried out on enzyme-inhibitor (two hydroxamates that only differ by a carbon-carbon double bond) complexes of human fibroblast collagenase to obtain insights into the structural and energetic preferences of these inhibitors. We have developed a bonded model for the catalytic and structural zinc centers (Hoops, S. C.; et al. J. Am. Chem. Soc. 1991, 113, 8262-8270) where the electrostatic representation for this model was derived using a novel quantum-mechanical/molecular-mechanical (QM/MM) minimization procedure followed by electrostatic potential fitting. The resulting bonded model for the zinc ions was then used to generate MD trajectories for structural analysis and FEP studies. This model has satisfactorily reproduced the structural features of the active site, and furthermore, the FEP simulations gave relative free energies of binding in good agreement with experimental results. MD simulations in conjunction with the FEP are able to provide a structural explanation regarding why one hydroxamate inhibitor is favored over the other, and we are also able to make predictions about changes in the inhibitor that would enhance protein-inhibitor interactions.

Interaction of small peptides with lipid bilayers

Molecular dynamics simulations of the tripeptide Ala-Phe-Ala-O-tert-butyl interacting with dimyristoylphosphatidylcholine lipid bilayers have been carried out. The lipid and aqueous environments of the peptide, the alkyl chain order, and the lipid and peptide dynamics have been investigated with use of density profiles, radial distribution functions, alkyl chain order parameter profiles, and time correlation functions. It appears that the alkyl chain region accommodates the peptides in the bilayer with minimal perturbation to this region. The peptide dynamics in the bilayer bound form has been compared with that of the free peptide in water. The peptide structure does not vary on the simulation time scale (of the order of hundreds of picoseconds) compared with the solution structure in which a random structure is observed.

A comparison of DMPC- and DLPE-based lipid bilayers

A 250 ps molecular dynamics simulation of the dimyristoylphosphatidylcholine (DMPC)-based lipid bilayer, including explicit water molecules, is reported. The solvent environment of the head groups and other structural properties of the bilayer have been analyzed and compared with experimental results as well as our previous simulation of the dilauroylphosphatidylethanolamine (DLPE)-based bilayer. From this comparison we find that the solvent structure around the DMPC head group (clathrate shell) is significantly different than that around the DLPE head group (typical hydrogen bonding interactions). We have modeled the probable relationship between the different solvent environments around the R-N(CH3)3+ (DMPC) and R-NH3+ (DLPE) head groups and the different interlammelar distances in these systems by performing potential of mean force (PMF) simulations on two N(CH3)4+ and NH4+ ions in water. From the PMF simulations it appears that the differences in the hydration of the DMPC and DLPE head groups is not responsible for the differences in the hydration force observed for these systems. We also find that the orientational polarization of DLPE and DMPC is similar, which suggests that solvent polarization is not responsible for the differences in the hydration repulsion behavior observed in these systems. We also examined the order parameters for DMPC and found them to be in reasonable agreement with experiment. Given the different characteristics of the DLPE and DMPC head groups, we suggest an explanation of the differences in the interlammellar spacings of bilayers composed of these like-charged lipids. From our DLPE simulations we find that the R-NH3+ head groups can interact with the nonesterified oxygens of the phosphate group in an intraleaflet or an interleaflet manner. For the latter a "cross link" between two leaflets can be formed, which causes a stabilization of the interlamellar spacings at fairly short distances. Moreover, due to the strong intraleaflet interaction we find that the DLPE interface is relatively "flat" (as opposed to DMPC-based bilayers), which results in a surface that has regions of positive and negative charge that reside in the same plane along the bilayer normal. Based on this we propose that the DLPE bilayer interface can correlate itself with another DLPE interface by alignment of the regions of positive (or negative) charge on one leaflet with the opposite charges on the opposing leaflet.

Structure and dynamics of the dilauroylphosphatidylethanolamine lipid bilayer

A 200-ps molecular dynamics (MD) simulation trajectory of a model dilauroylphosphatidylethanolamine (DLPE) bilayer in water at 315 K has been generated. Segmental order parameters, electron density profiles, and water pair distribution functions have been calculated. Comparison to experiment is made where possible. The dynamics of the system has been studied by analyzing the velocity autocorrelation functions (VAF) of both water and lipid atoms. Furthermore, the diffusive properties of water have been analyzed by computing the mean square displacement (MSD) and orientational correlation function (OCF) of water in two regions around the bilayer. The calculated order parameters show a behavior similar to the liquid crystalline phase of other bilayers, but the region around C1-C3 does not show the expected behavior. The electron density profile shows features that are characteristic of the liquid crystalline phase. The radial distribution functions suggest ordering of water near the charged head groups, which results in about 15 water molecules solvating each lipid molecule. We find from the VAF, MSD, and OCF calculation that the water molecules near the head groups of the lipid bilayer move more slowly than those further away. The VAF of the hydrocarbon chains have features of low-frequency motions that are probably cooperative nature in addition to the high-frequency motions associated with bond angle and torsional motions.