CHARMM: the biomolecular simulation program.

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

Molecular dynamics simulations of human rhinovirus and an antiviral compound.

The human rhinovirus 14 (HRV14) protomer, with or without the antiviral compound WIN 52084s, was simulated using molecular dynamics and rotational symmetry boundary conditions to model the effect of the entire icosahedral capsid. The protein asymmetrical unit, comprising four capsid proteins (VP1, VP2, VP3, and VP4) and two calcium ions, was solvated both on the exterior and the interior to fill the inside of the capsid. The stability of the simulations of this large system (~800 residues and 6,650 water molecules) is comparable to more conventional globular protein simulations. The influence of the antiviral compound on compressibility and positional fluctuations is reported. The compressibility, estimated from the density fluctuations in the region of the binding pocket, was found to be greater with WIN 52084s bound than without the drug, substantiating previous computations on reduced viral systems. An increase in compressibility correlates with an entropically more favorable system. In contrast to the increase in density fluctuations and compressibility, the positional fluctuations decreased dramatically for the external loops of VP1 and the N-terminus of VP3 when WIN 52084s is bound. Most of these VP1 and VP3 loops are found near the fivefold axis, a region whose mobility was not considered in reduced systems, but can be observed with this simulation of the full viral protomer. Altered loop flexibility is consistent with changes in proteolytic sensitivity observed experimentally. Moreover, decreased flexibility in these intraprotomeric loops is noteworthy since the externalization of VP4, part of VP1, and RNA during the uncoating process is thought to involve areas near the fivefold axis. Both the decrease in positional fluctuations at the fivefold axis and the increase in compressibility near the WIN pocket are discussed in relationship to the antiviral activity of stabilizing the virus against uncoating.

Accuracy of bound peptide structures determined by exchange transferred nuclear Overhauser data: a simulation study.

The exchange-transferred NOE method to determine the three-dimensional structure of peptides bound to proteins, or other macromolecular systems, is becoming increasingly important in drug design efforts and for large or multicomponent assemblies, such as membrane receptors, where structural analysis of the full system is intractable. The exchange-transferred nuclear Overhauser effect spectroscopy (etNOESY) method allows the determination of the bound-state conformation of the peptide from the intra-molecular NOE interactions between ligand protons. Because only ligand-ligand NOEs are generally observable, the etNOESY method is restricted to fewer NOEs per residue than direct protein structure determination. In addition, the averaging of relaxation rates between free and bound states affects the measured cross-peak intensities, and possibly the accuracy of distance estimates. Accordingly, the study reported here was conducted to examine the conditions required to define a reliable structure. The program CORONA was used to simulate etNOE data using a rate-matrix including magnetic relaxation and exchange rates for two peptide-protein complexes derived from the reference complex of cAMP-dependent protein kinase ligated with a 24-residue inhibitor peptide. The results indicate that reasonably accurate peptide structures can be determined with relatively few NOE interactions when the interactions occur between non-neighboring residues. The reliability of the structural result is suggested from the pattern of NOE interactions. A structure with an accuracy of approximately 1.3 A rms difference for the main-chain atoms can be obtained when etNOE interactions between non-neighboring residues occur over the length of the peptide. The global precision is higher (approximately 0.9 A rms difference) but is not correlated to global accuracy. A local definition of precision along the backbone appears to be a good indicator of the local accuracy.

Structure and self-association of the Rous sarcoma virus capsid protein.

BACKGROUND: The capsid protein (CA) of retroviruses, such as Rous sarcoma virus (RSV), consists of two independently folded domains. CA functions as part of a polyprotein during particle assembly and budding and, in addition, forms a shell encapsidating the genomic RNA in the mature, infectious virus. RESULTS: The structures of the N- and C-terminal domains of RSV CA have been determined by X-ray crystallography and solution nuclear magnetic resonance (NMR) spectroscopy, respectively. The N-terminal domain comprises seven alpha helices and a short beta hairpin at the N terminus. The N-terminal domain associates through a small, tightly packed, twofold symmetric interface within the crystal, different from those previously described for other retroviral CAs. The C-terminal domain is a compact bundle of four alpha helices, although the last few residues are disordered. In dilute solution, RSV CA is predominantly monomeric. We show, however, using electron microscopy, that intact RSV CA can assemble in vitro to form both tubular structures constructed from toroidal oligomers and planar monolayers. Both modes of assembly occur under similar solution conditions, and both sheets and tubes exhibit long-range order. CONCLUSIONS: The tertiary structure of CA is conserved across the major retroviral genera, yet sequence variations are sufficient to cause change in associative behavior. CA forms the exterior shell of the viral core in all mature retroviruses. However, the core morphology differs between viruses. Consistent with this observation, we find that the capsid proteins of RSV and human immunodeficiency virus type 1 exhibit different associative behavior in dilute solution and assemble in vitro into different structures.

Alkylation of cysteine-containing peptides to mimic palmitoylation.

Numerous proteins that are involved in cell signaling and viral replication require post-translational modification by palmitoylation to function properly. The molecular details by which this palmitoyl modification affects protein function remain poorly understood. To facilitate in vitro biochemical and structural studies of the role of palmitoylation on protein function, a method was developed for alkylating peptides with saturated C16 groups at cysteine residues and demonstrated using peptides derived from the palmitoylated region of Sindbis virus E2 glycoprotein. The synthetic approach takes advantage of disulfide chemistry to specifically modify only the cysteine residues within peptides and covalently links C16 groups via disulfide bridges using a new thioalkylating reagent, hexyldexyldithiopyridine. The chemistry presented here takes place in solution under mild conditions without the need for protection of the peptide functional groups. A method for purifying these modified peptides is also described. This protocol can be of general use to investigators studying the role of palmitoylation in biological systems.

Substrate recognition by the Lyn protein-tyrosine kinase. NMR structure of the immunoreceptor tyrosine-based activation motif signaling region of the B cell antigen receptor.

The immunoreceptor tyrosine-based activation motif (ITAM) plays a central role in transmembrane signal transduction in hematopoietic cells by mediating responses leading to proliferation and differentiation. An initial signaling event following activation of the B cell antigen receptor is phosphorylation of the CD79a (Ig-alpha) ITAM by Lyn, a Src family protein-tyrosine kinase. To elucidate the structural basis for recognition between the ITAM substrate and activated Lyn kinase, the structure of an ITAM-derived peptide bound to Lyn was determined using exchange-transferred nuclear Overhauser NMR spectroscopy. The bound substrate structure has an irregular helix-like character. Docking based on the NMR data into the active site of the closely related Lck kinase strongly favors ITAM binding in an orientation similar to binding of cyclic AMP-dependent protein kinase rather than that of insulin receptor tyrosine kinase. The model of the complex provides a rationale for conserved ITAM residues, substrate specificity, and suggests that substrate binds only the active conformation of the Src family tyrosine kinase, unlike the ATP cofactor, which can bind the inactive form.

Theoretical studies of viral capsid proteins.

Recent results in structural biology and increases in computer power have prompted initial theoretical studies on capsids of nonenveloped icosahedral viruses. The macromolecular assembly of 60 to 180 protein copies into a protein shell results in a structure of considerable size for molecular dynamics simulations. Nonetheless, progress has been made in examining these capsid assemblies from molecular dynamics calculations and kinetic models. The goals of these studies are to understand capsid function and structural properties, including quarternary structural stability, effects of antiviral compounds that bind the capsid and the self-assembly process. The insight that can be gained from the detailed information provided by simulations is demonstrated in studies of human rhinovirus; an entropic basis for the antiviral activity of hydrophobic compounds, predicted from calculated compressibility values, has been corroborated by experimental measurements on poliovirus.

Molecular dynamics investigation of the effect of an antiviral compound on human rhinovirus.

The factors that influence the enhanced stability observed experimentally of human rhinovirus 14 (HRV14) upon binding a hydrophobic antiviral drug have been investigated by molecular dynamics. Simulations centered about the HRV14 drug-binding pocket allow the reliable assessment of differences in capsid protein motions of HRV14 and drug-bound HRV14. We propose that the experimentally observed stabilization of the ligated virus arises from higher entropy, rather than enthalpy. Time-averaged interaction energies between the viral protein and molecules occupying the pocket are less favorable in the presence of the drug, consistent with the proposal that the observed stability arises from entropic effects. Interaction energies characterizing subunit-subunit contacts within one viral protomer are found to be substantially stronger than those between two protomers. Such distinction in subunit interaction would have clear implications on assembly and disassembly. Drug binding is found to affect large-scale, collective properties, while leaving local atomic properties unperturbed. Specifically, the simulations reveal a weakening of long-range correlations in atomic motions upon drug binding. On the other hand, neither the fast time scale RMS fluctuations of individual atomic positions nor the fluctuation build-up curves from the capsid beta-sandwich forming the drug-binding pocket show a consistent distinction between the drug-bound and drug-free viral simulations. Collectively, the detailed description available from the simulations provides an understanding of the experimental observations on the drug-induced changes in thermal stability and protease sensitivity reported for picornaviruses. The predicted significance of binding entropy can be explored experimentally and should be considered in the design of new antiviral compounds.

NMR structure of phospho-tyrosine signaling complexes.

A structural basis for activation and substrate specificity of src tyrosine kinases, and regulation of protein-protein association by tyrosine phosphorylation is described. Lyn, a src-family tyrosine kinase, recognizes and phosphorylates the immunoreceptor tyrosine-based activation motif, ITAM, a critical component in transmembrane signal transduction in hemopoietic cells. The structure of an ITAM peptide substrate bound to an active form of Lyn tyrosine kinase was determined by high-resolution NMR, and a model of the complex was generated using the crystallographic structure of Lck, a closely related Src-family kinase. The results provide a rationale for the conserved ITAM residues and specificity of Lyn, and suggest that substrate plays a role in stabilizing the kinase conformation optimal for catalysis. It is our hope that the Lck-ITAM peptide model complex will be useful in aiding structure-based drug design efforts that target substrate binding determinants in the design. Concerning the regulation of protein-protein association, we report on a complex between erythrocyte band 3 and two glycolytic enzymes, aldolase and glyceraldehyde-3-phosphate dehydrogenase. The formation of this complex is negatively regulated by tyrosine phosphorylation of band 3 by p72syk tyrosine kinase. In red blood cells, this association results in a decrease in glycolysis due to competitive inhibition of the glycolytic enzymes. The structure of band 3 recognized by the glycolytic enzymes was determined by solution NMR, and found to be a loop structure with tyrosine centrally positioned and excluded from intermolecular contact. This phosphorylation sensitive interaction, or PSI, loop may be the basis of a general mechanism for negative regulation through tyrosine phosphorylation.

Influence of an antiviral compound on the temperature dependence of viral protein flexibility and packing: a molecular dynamics study.

The antiviral activity of compounds that bind an internal pocket of picornaviruses is due in part to stabilization of the protein capsid and inhibition of the uncoating process required for virus replication. Information on the basis for this structural stabilization of the virus capsid is important to elucidate the mechanism of antiviral action and provide insights into the disassembly process. It has been proposed that this stabilization is entropically based, since binding the nonpolar antiviral compound increases the compressibility, and thus the conformational flexibility, of the virus. Such a proposal predicts a difference in the temperature dependence of the atomic positional fluctuations for free virus and drug-bound virus; nonpolar interactions are weaker and less directional, and would give rise to greater conformational disorder at low temperature. Further, the transition that has been observed in globular proteins to a state resembling a frozen liquid, in which the protein is considered "trapped" in potential energy wells, is predicted to occur at lower temperature when the antiviral compound is bound. Results described here from computer simulations of rhinovirus over a range in temperature show these predicted changes in conformational disorder and the temperature of the transition in mobility. In addition to providing independent support for the above proposal for antiviral activity, these results indicate that the mobility transition of a protein can be controlled by the binding of an appropriate ligand, an effect not previously reported.

Insights into tyrosine phosphorylation control of protein-protein association from the NMR structure of a band 3 peptide inhibitor bound to glyceraldehyde-3-phosphate dehydrogenase.

A protein-protein association regulated by phosphorylation of tyrosine is examined by NMR structural studies and biochemical studies. Binding of glyceraldehyde-3-phosphate dehydrogenase (G3PDH) and aldolase to the N-terminus of human erythrocyte anion transporter, band 3, inhibits enzyme activity. This inhibition is reversed upon phosphorylation of band 3 Y8, as shown by kinetic studies on purified components, as well as in vivo studies. Thus, tyrosine phosphorylation mediates against the intermolecular protein-protein association, in contrast to the positive control involving SH2 and PTB domains where phosphorylation is required for binding. To elucidate the basis of recognition and negative control by tyrosine phosphorylation, the structure of a synthetic peptide, B3P, corresponding to the first 15 residues of band 3 (MEELQDDYEDMMEEN-NH2), bound to G3PDH has been determined using the exchange-transferred nuclear Overhauser effect. The G3PDH-bound B3P structure was found to be very similar to the structure recognized by aldolase. A hydrophobic triad forms from side chains within a loop structure of residues 4 through 9 in both bound species. Another structural feature stabilizing the loop, in the case of the B3P-G3PDH complex, is a hydrogen bond between the side chains of Y8 and D10 associated with a beta-turn of residues 8-11. Based on the structure of this phosphorylation sensitive interaction (PSI) loop, it is suggested that tyrosine phosphorylation disrupts protein-protein association, in part, by intramolecular electrostatic destabilization. The inhibition by B3P is competitive with respect to the coenzyme NAD+ and noncompetitive with the substrate analog arsenate. Specific binding of B3P to G3PDH is demonstrated by reversion of the NMR spectral properties of bound B3P to those of the free peptide upon addition of coenzyme and substrate analog. The stoichiometry of binding for the B3P-G3PDH complex was determined from Sephadex G-50 displacement experiments to be 4:1. Collectively, these results are consistent with B3P binding the active site of G3PDH.

Energetics of quasiequivalence: computational analysis of protein-protein interactions in icosahedral viruses.

Quaternary structure polymorphism found in quasiequivalent virus capsids provides a static framework for studying the dynamics of protein interactions. The same protein subunits are found in different structural environments within these particles, and in some cases, the molecular switching required for the polymorphic quaternary interactions is obvious from high-resolution crystallographic studies. Employing atomic resolution structures, molecular mechanics, and continuum electrostatic methods, we have computed association energies for unique subunit interfaces of three icosahedral viruses, black beetle virus, southern bean virus, and human rhinovirus 14. To quantify the chemical determinants of quasiequivalence, the energetic contributions of individual residues forming quasiequivalent interfaces were calculated and compared. The potential significance of the differences in stabilities at quasiequivalent interfaces was then explored with the combinatorial assembly approach. The analysis shows that the unique association energies computed for each virus serve as a sensitive basis set that may determine distinct intermediates and pathways of virus capsid assembly. The pathways for the quasiequivalent viruses displayed isoenergetic oligomers at specific points, suggesting that these may determine the quaternary structure polymorphism required for the assembly of a quasiequivalent particle.

The conformation of NAD+ bound to lactate dehydrogenase determined by nuclear magnetic resonance with suppression of spin diffusion.

We have reinvestigated the conformation of NAD+ bound to dogfish lactate dehydrogenase (LDH) by using an NMR experiment that allows one to exploit nuclear Overhauser effects to determine internuclear distances between pairs of protons, without perturbation of spin-diffusion effects from other protons belonging either to the cofactor or to the binding pocket of the enzyme. The analysis indicates that the structure of bound NAD+ is in accord with the conformation determined in the solid state by x-ray diffraction for the adenosine moiety, but deviates significantly from that of the nicotinamide. The NMR data indicate conformational averaging about the glycosidic bond of the nicotinamide nucleotide. In view of the strict stereospecificity of catalysis by LDH and the conformational averaging of bound NAD+ that we infer from solution-state NMR, we suggest that LDH binds the cofactor in both syn and anti conformations, but that binding interactions in the syn conformation are not catalytically productive.

Catalysis by entropic guidance from enzymes.

Certain features underlying enzymatic catalysis, such as energetic stabilization from binding interactions or proximity and orientation of chemical groups, are evident in the equilibrium-averaged structure of an enzymatic complex determined by crystallography or NMR. Transient features are not apparent from an average structure. Here, we report on a catalytically relevant property of an enzymatic complex revealed by thermal fluctuations from a molecular dynamics study. The conformational fluctuations of the cofactor NADH are altered by binding the enzyme lactate dehydrogenase (LDH) compared to those of free NADH; thermal motions give rise to structures similar to that of the putative transition state. The alteration is stereospecific, in agreement with measured changes in vibrational spectra, and leads to an understanding of the correlation, established some time ago by crystallography and NMR, between the nicotinamide glycosidic bond torsion angle (anti/syn) and the stereospecificity of hydride transfer. These results suggest that one catalytic role of the enzyme is to funnel the population of NADH conformers to the transition state and reduce the entropic barrier to activation. The specific motions in an enzyme complex that might function to enhance transition state formation are described.

Amide hydrogen exchange determined by mass spectrometry: application to rabbit muscle aldolase.

The protein fragmentation/mass spectrometry method described by Zhang and Smith [(1993) Protein Sci. 2, 522-531] has been extended to measure amide hydrogen exchange rates in rabbit muscle aldolase, a homotetramer with M(r) = 157,000. Following a period of deuterium exchange, the partially deuterated protein was proteolytically fragmented into peptides whose deuterium contents were determined by directly coupled HPLC fast atom bombardment mass spectrometry. Hydrogen exchange rates were determined for amide hydrogens located in short segments derived from 85% of the aldolase backbone. Isotopic exchange rate constants spanning the range from 100 to 0.001 h-1 were determined for the exchange-in times used in this study (2.5 min to 44 h). The exchange rates for amide hydrogens located within short segments differed by as much as 10(4), demonstrating that local structural features dramatically affect the isotopic exchange rates in large proteins. A high level of correlation between the slowing of hydrogen exchange and intramolecular hydrogen bonding in aldolase was found. An exception to this correlation occurs at the subunit interface, where the amide hydrogens in one peptide segment with few amide hydrogen bonds have slower exchange rates than expected, suggesting that the amide hydrogens in this region are effectively shielded from the deuterated solvent. Isotope patterns observed for most peptides were binomial, indicating that hydrogen exchange proceeds through the EX2 mechanism (uncorrelated exchange). However, bimodal isotope patterns were found for peptides derived from three short segments of aldolase (including residues 58-64, 279-283, and 326-337), suggesting structural differences in these regions. A high level of correlation was found between crystallographic B-factors and amide hydrogen exchange rates, suggesting an isotopic exchange mechanism involving localized low-amplitude, high-frequency motions that do not require collective motion of many residues. From a methodology viewpoint, these results demonstrate that the combination of protein fragmentation with mass spectrometry is a useful method for determining the rates at which amide hydrogens located over major portions of large proteins undergo isotopic exchange.

Simulations of lysozyme: internal motions and the reaction mechanism.

Simulations of lysozyme by molecular dynamics have greatly increased our understanding of this enzyme. It has been shown how the internal motions are related to the structural elements (helices, sheets and loops) of the molecule. Comparisons of the motions in the free and substrate-bound form reveal that most are similar but that there are significant differences. Comparisons of the theoretical results with X-ray and nuclear magnetic resource data show good agreement. The hinge-bending motion, which opens and closes the binding site cleft between the two domains, is correlated with substrate binding. From analysis of simulations with bound substrate, an alternative mechanism for oligoglycoside hydrolysis was proposed. It involves cleavage of an endocyclic C-O bond, instead of the exocyclic cleavage proposed in the standard mechanism. Both mechanisms have been demonstrated in solution, but it is still unclear which is prevalent in lysozyme.

Solution structure of a band 3 peptide inhibitor bound to aldolase: a proposed mechanism for regulating binding by tyrosine phosphorylation.

Human erythrocyte band 3 inhibits glycolytic enzymes, including aldolase, by binding these cytoplasmic enzymes at its N-terminus. Phosphorylation of Y8 disrupts inhibition, and there is evidence that in vivo glycolysis levels in erythrocytes are regulated in part by a phosphorylation/dephosphorylation signaling pathway. The structural basis for control by phosphorylation has been investigated by NMR studies on a complex between aldolase and a synthetic peptide corresponding to the first 15 residues of band 3 (MEELQDDYEDMMEEN-NH2). The structure of this band 3 peptide (B3P) when it is bound to rabbit muscle aldolase was determined using the exchange-transferred nuclear Overhauser effect (ETNOE). Two hundred NMR structures for B3P were generated by simulated annealing molecular dynamics with NMR-derived distance restraints and excluding electrostatic terms. Twenty structures were further refined against a force field including full partial charges. The important conformational feature of B3P in the bound state is a folded loop structure involving residues 4-9 and M12 that surrounds Y8 and is stabilized by a hydrophobic cluster with the ring of Y8 sandwiched between the methyl groups of L4 and M12. Differential line broadening indicates that this loop structure binds aldolase in a relatively specific manner, while terminal regions are structurally heterogeneous. To better understand B3P inhibition of aldolase and the mechanism of phosphorylation control, a complex was modeled by docking B3P into the active site of aldolase and optimizing the fit using restrained molecular dynamics and energy minimization. The B3P loop is complementary in conformation to the beta-barrel central core containing the aldolase active site residues. Binding is electrostatic in nature with numerous ionic and hydrogen-bonding interactions involving several conserved lysine and arginine residues of aldolase. How phosphorylation of band 3 could disrupt inhibition was considered by modeling a phosphoryl moiety onto Y8 of B3P. An energetic analysis with respect to rigid phosphate rotation suggests that aldolase inhibition is reversed primarily because of electrostatic repulsion between B3P residues that destabilizes the B3P loop formed in the complex. This proposed intramolecular mechanism for blocking protein--protein association by electrostatic repulsion with the phosphoryl group may be applicable to other protein--protein signaling complexes.

Reexamination of induced fit as a determinant of substrate specificity in enzymatic reactions.

It has been argued that a substrate-induced conformational change involving the orientation of catalytic groups cannot affect the specificity for two substrates in an enzymatic system where the chemical step is rate limiting, because such an induced fit would alter the catalytic efficiency for both to an equal extent. To the contrary, the generalized induced-fit treatment described here shows that when critical substrate-specific conformational changes in the enzyme persist in the transition state, specificity is linked to conformational differences between the reactive complex for a good substrate and the related complex for a poor one. Conformational differences are a determinant of specificity when the reaction proceeds via an "induced-fit" transition state. Our treatment also shows that such conformational changes can enhance the specificity of an enzyme with suboptimal catalytic efficiency. If substrate-dependent conformational differences in a primative enzyme can enhance specificity, evolutionary pressure to increase specificity could inseparably link enzymatic specificity to induced conformational changes.

A novel basis of capsid stabilization by antiviral compounds.

Picornaviruses are inactivated by a family of hydrophobic drugs that bind at an internal site in the viral capsid and inhibit viral uncoating. A basis for the capsid stabilization previously unrecognized is revealed by molecular dynamics simulations of the antiviral drug WIN52084s bound to a hydrophobic pocket of solvated human rhinovirus 14. Isothermal compressibilities of the complex and human rhinovirus 14 without the antiviral drug calculated from density fluctuations show that the presence of WIN52084s increases the compressibility of the viral capsid near the antiviral drug. This counterintuitive result is understandable on the basis of the empirical evidence of thermal melting temperatures and protein-folding entropies of globular proteins. Based on this evidence, we propose that a larger compressibility from drug binding confers greater thermal stability to capsid proteins by increasing the conformational entropy of capsids, thereby diminishing the entropy gain with uncoating. We suggest that compressibility is fundamental to the structural integrity of viral capsids and that examination of compressibility and antiviral activity will provide insights into the disassembly process.