Search

Molecular dynamics (MD) force fields for lipids and ions are typically developed independently of one another. In simulations consisting of both lipids and ions, lipid-ion interaction energies are estimated using a predefined set of mixing rules for Lennard-Jones (LJ) interactions. This, however, does not guarantee their reliability. In fact, compared to the quantum mechanical reference data, Lorentz-Berthelot mixing rules substantially underestimate binding energies of Na(+) ions with small molecule analogues of lipid headgroups, yielding errors on the order of 80 and 130 kJ/mol, respectively for methyl acetate and diethyl phosphate. Previously, errors associated with mixing force fields have been reduced using approaches like ‘NB-fix’ in which LJ interactions are computed using explicit cross terms rather than those from mixing rules. Building on this idea, we derive explicit lipid-ion cross terms that also may implicitly include many-body cooperativity effects. Additionally, to account for interdependency between cross terms, we optimize all cross terms simultaneously by performing high-dimensional searches using our ParOpt software. The cross terms we obtain reduce the errors due to mixing rules to below 10 kJ/mol. MD simulation of lipid bilayer conducted using these optimized cross terms resolve the structural discrepancies between our previous simulations and small-angle X-ray and neutron scattering experiments. These results demonstrate that simulations of lipid bilayers with ions that are accurate up to structural data from scattering experiments can be performed without explicit polarization terms. However, it is worth noting that such NB-fix cross terms are not based on any physical principle; a polarizable lipid model would be more realistic, and is still desired. Our approach is generic and can be applied to improve accuracies of simulations employing mixed force fields.

Bloom syndrome helicase (BLM) is a RecQ-family helicase implicated in a variety of cellular processes, including DNA replication, DNA repair, and telomere maintenance. Mutations in humanBLMcause Bloom syndrome (BS), an autosomal recessive disorder that leads to myriad negative health impacts including a predisposition to cancer. BS-causing mutations inBLMoften negatively impact BLM ATPase and helicase activity. WhileBLMmutations that cause BS have been well characterized bothin vitroandin vivo, there are other less studiedBLMmutations that exist in the human population that do not lead to BS. Two of these non-BS mutations, encoding BLM P868L and BLM G1120R, when homozygous, increase sister chromatid exchanges in human cells. To characterize these naturally occurring BLM mutant proteinsin vitro, we purified the BLM catalytic core (BLMcore, residues 636-1298) with either the P868L or G1120R substitution. We also purified a BLMcoreK869A K870A mutant protein, which alters a lysine-rich loop proximal to the P868 residue. We found that BLMcoreP868L and G1120R proteins were both able to hydrolyze ATP, bind diverse DNA substrates, and unwind G-quadruplex and duplex DNA structures. Molecular dynamics simulations suggest that the P868L substitution weakens the DNA interaction with the winged-helix domain of BLM and alters the orientation of one lobe of the ATPase domain. Because BLMcoreP868L and G1120R retain helicase functionin vitro, it is likely that the increased genome instability is caused by specific impacts of the mutant proteinsin vivo. Interestingly, we found that BLMcoreK869A K870A has diminished ATPase activity, weakened binding to duplex DNA structures, and less robust helicase activity compared to wild-type BLMcore. Thus, the lysine-rich loop may have an important role in ATPase activity and specific binding and DNA unwinding functions in BLM.

Lysine methylation can modify noncovalent interactions by altering lysine's hydrophobicity as well as its electronic structure. Although the ramifications of the former are documented, the effects of the latter remain largely unknown. Understanding the electronic structure is important for determining how biological methylation modulates protein-protein binding, and the impact of artificial methylation experiments in which methylated lysines are used as spectroscopic probes and protein crystallization facilitators. The benchmarked first-principles calculations undertaken here reveal that methyl-induced polarization weakens the electrostatic attraction of amines with protein functional groups - salt bridges, hydrogen bonds and cation-π interactions weaken by as much as 10.3, 7.9 and 3.5 kT, respectively. Multipole analysis shows that weakened electrostatics is due to the altered inductive effects, which overcome increased attraction from methyl-enhanced polarizability and dispersion. Due to their fundamental nature, these effects are expected to be present in many cases. A survey of methylated lysines in protein structures reveals several cases in which methyl-induced polarization is the primary driver of altered noncovalent interactions; in these cases, destabilizations are found to be in the 0.6-4.7 kT range. The clearest case of where methyl-induced polarization plays a dominant role in regulating biological function is that of the PHD1-PHD2 domain, which recognizes lysine-methylated states on histones. These results broaden our understanding of how methylation modulates noncovalent interactions.

Collagen is heavily hydroxylated. Experiments show that proline hydroxylation is important to triple helix (monomer) stability, fibril assembly, and interaction of fibrils with other molecules. Nevertheless, experiments also show that even without hydroxylation, type I collagen does assemble into its native D-banded fibrillar structure. This raises two questions. Firstly, even though hydroxylation removal marginally affects macroscopic structure, how does such an extensive chemical change, which is expected to substantially reduce hydrogen bonding capacity, affect local structure? Secondly, how does such a chemical perturbation, which is expected to substantially decrease electrostatic attraction between monomers, affect collagen’s mechanical properties? To address these issues, we conduct a benchmarked molecular dynamics study of rat type I fibrils in the presence and absence of hydroxylation. Our simulations reproduce the experimental observation that hydroxylation removal has a minimal effect on collagen’s D-band length. We also find that the gap-overlap ratio, monomer width and monomer length are minimally affected. Surprisingly, we find that de-hydroxylation also has a minor effect on the fibril’s Young’s modulus, and elastic stress build up is also accompanied by tightening of triple-helix windings. In terms of local structure, de-hydroxylation does result in a substantial drop (23%) in inter-monomer hydrogen bonding. However, at the same time, the local structures and inter-monomer hydrogen bonding networks of non-hydroxylated amino acids are also affected. It seems that it is this intrinsic plasticity in inter-monomer interactions that preclude fibrils from undergoing any large changes in macroscopic properties. Nevertheless, changes in local structure can be expected to directly impact collagen’s interaction with extra-cellular matrix proteins. In general, this study highlights a key challenge in tissue engineering and medicine related to mapping collagen chemistry to macroscopic properties but suggests a path forward to address it using molecular dynamics simulations.

Lysine methylation is a key regulator of protein-protein binding. The amine group of lysine can accept up to three methyl groups, and experiments show that protein-protein binding free energies are sensitive to the extent of methylation. These sensitivities have been rationalized in terms of chemical and structural features present in the binding pockets of methyllysine binding domains. However, understanding their specific roles requires an energetic analysis. Here we propose a theoretical framework to combine quantum and molecular mechanics methods, and compute the effect of methylation on protein-protein binding free energies. The advantages of this approach are that it derives contributions from all local non-trivial effects of methylation on induction, polarizability and dispersion directly from self-consistent electron densities, and at the same time determines contributions from well-characterized hydration effects using a computationally efficient classical mean field method. Limitations of the approach are discussed, and we note that predicted free energies of fourteen out of the sixteen cases agree with experiment. Critical assessment of these cases leads to the following overarching principles that drive methylation-state recognition by protein domains. Methylation typically reduces the pairwise interaction between proteins. This biases binding toward lower methylated states. Simultaneously, however, methylation also makes it easier to partially dehydrate proteins and place them in protein-protein complexes. This latter effect biases binding in favor of higher methylated states. The overall effect of methylation on protein-protein binding depends ultimately on the balance between these two effects, which is observed to be tuned via several combinations of local features.

Ion descriptors in molecular mechanics models are calibrated against reference data on ion-water interactions. It is then typically assumed that these descriptors will also satisfactorily describe interactions of ions with other functional groups, such as those present in biomolecules. However, several studies now demonstrate that this transfer-ability assumption produces, in many different cases, large errors. Here we address this issue in a representative polarizable model and focus on transferability of cationic interactions from water to a series of alcohols. Both water and alcohols use hydroxyls for ion-coordination, and, therefore, this set of molecules constitutes the simplest possible case of transferability. We obtain gas phase reference data systematically from “gold-standard” quantum Monte Carlo and CCSD(T) methods, followed by bench-marked vdW-corrected DFT. We learn that the original polarizable model yields large gas phase water→alcohol transferability errors – the RMS and maximum errors are 2.3 and 5.1 kcal/mol, respectively. These errors are, nevertheless, systematic in that ionalcohol interactions are over-stabilized, and systematic errors typically imply that some essential physics is either missing or misrepresented. A comprehensive analysis shows that when both low and high-field responses of ligand dipole polarization are described accurately, then transferability improves significantly – the RMS and maximum errors in the gas phase reduce, respectively, to 0.9 and 2.5 kcal/mol. Additionally, predictions of condensed phase transfer free energies also improve. Nevertheless, within the limits of the extra-thermodynamic assumptions necessary to separate experimental estimates of salt dissolution into constituent cationic and anionic contributions, we note that the error in the condensed phase is systematic, which we attribute, at least, partially to the parameterization in long range electrostatics. Overall, this work demonstrates a rational approach to boosting transferability of ionic interactions that will be applicable broadly to improving other polarizable and non-polarizable models.

A distinguishing feature of Archaeal plasma membranes is that their phospho-lipids contain ether-links, as opposed to bacterial and eukaryotic plasma mem-branes where phospholipids primarily contain ester-links. Experiments show that this chemical difference in headgroup-tail linkage does produce distinct differ-ences in model bilayer properties. Here we examine the effects of salt on bilayer structure in the case of an ether-linked lipid bilayer. We use molecular dynamics simulations and compare equilibrium properties of two model lipid bilayers in NaCl salt solution - POPC and its ether-linked analog that we refer to as HOPC. We make the following key observations. The headgroup region of HOPC “ad-sorbs” fewer ions compared to the headgroup region of POPC. Consistent with this, we note that the Debye screening length in the HOPC system is ~10% shorter than that in the POPC system. Herein, we introduce a protocol to identify the lipid-water interfacial boundary that reproduces the bulk salt distribution consis-tent with Gouy-Chapman theory. We also note that the HOPC bilayer has ex-cess solvent in the headgroup region when compared to POPC, coinciding with a trough in the electrostatic potential. Waters in this region have longer autocor-relation times and smaller lateral diffusion rates compared to the corresponding region in the POPC bilayer, suggesting that the waters in HOPC are more strongly coordinated to the lipid headgroups. Furthermore, we note that it is this region of tightly coordinated waters in the HOPC system that has a lower density of Na(+) ions. Based on these observations we conclude that an ether-linked lipid bilayer has a lower binding affinity for Na(+) compared to an ester-linked lipid bilayer.

Phosphohistidine phosphatase 1 (PHPT1), the only known phosphohistidine phosphatase in mammals, regulates phosphohistidine levels of several proteins including those involved in signaling, lipid metabolism, and potassium ion transport. While the high-resolution structure of human PHPT1 (hPHPT1) is available and residues important for substrate binding and catalytic activity have been reported, little is known about post-translational modifications that modulate hPHPT1 activity. Here we characterize the structural and functional impact of hPHPT1 oxidation upon exposure to a reactive oxygen species, hydrogen peroxide (H2O2). Specifically, liquid chromatography-tandem mass spectrometry was used to quantify site-specific oxidation of redox-sensitive residues of hPHPT1. Results from this study revealed that H2O2 exposure induces selective oxidation of hPHPT1 at Met95, a residue within the substrate binding region. Explicit solvent molecular dynamics simulations, however, predict only a minor effect of Met95 oxidation in the structure and dynamics of the apo-state of the hPHPT1 catalytic site, suggesting that if Met95 oxidation alters hPHPT1 activity, then it will do so by altering the stability of an intermediate state. Employing a novel mass spectrometry-based assay, we determined that H2O2–induced oxidation does not impact hPHPT1 function negatively; a result contrary to the common conception that protein oxidation is typically a loss-of-function modification.

The physical properties of lipid bilayers can be remodeled by a variety of environmental factors. Here we investigate using molecular dynamics simulations the specific effects of nanoscopic substrates or external contact points on lipid membranes. We expose palmitoyl-oleoyl phosphatidylcholine bilayers unilaterally and separately to various model nano-sized substrates differing in surface hydroxyl densities. We find that a surface hydroxyl density as low as 10% is sufficient to keep the bilayer juxtaposed to the substrate. The bilayer interacts with the substrate indirectly through multiple layers of water molecules, however, despite such buffered interaction, the bilayers exhibit certain properties different from unsupported bilayers. The substrates modify transverse lipid fluctuations, charge density profiles and lipid diffusion rates, although differently in the two leaflets, which creates an asymmetry between bilayer leaflets. Other properties that include lipid cross-sectional areas, component volumes and order parameters are minimally affected. The extent of asymmetry that we observe between bilayer leaflets is well beyond what has been reported for bilayers adsorbed on infinite solid supports. This is perhaps because the bilayers are much closer to our nano-sized finite supports than to infinite solid supports, resulting in a stronger support-bilayer electrostatic coupling. The exposure of membranes to nanoscopic contact points, therefore, cannot be considered as a simple linear interpolation between unsupported membranes and membranes supported on infinite supports. In the biological context, this suggests that the exposure of membranes to non-intercalating proteins, such as those belonging to the cytoskeleton, should not always be considered as passive non-consequential interactions.

During Ca(2+) activation, the Ca(2+)-binding sites of C2 domains typically bind multiple Ca(2+) ions in close proximity. These binding events exhibit positive cooperativity, despite the strong charge repulsion between the adjacent divalent cations. Using both experimental and computational approaches, the present study probes the detailed mechanisms of Ca(2+) activation and positive cooperativity for the C2 domain of cytosolic phospholipase A(2), which binds two Ca(2+) ions in sites I and II, separated by only 4.1 A. First, each of the five coordinating side chains in the Ca(2+)-binding cleft is individually mutated and the effect on Ca(2+)-binding affinity and cooperativity is measured. The results identify Asp 43 as the major contributor to Ca(2+) affinity, while the two coordinating side chains that provide bridging coordination to both Ca(2+) ions, Asp 43 and Asp 40, are observed to make the largest contributions to positive cooperativity. Electrostatic calculations reveal that Asp 43 possesses the highest pseudo-pK(a) of the coordinating acidic residues, as well as the highest general cation affinity, due to its relatively buried location within 3.5 A of seven protein oxygens with full or partial negative charges. These calculations therefore explain the greater importance of Asp 43 in defining the Ca(2+) affinity. Overall, the experimental and computational results support an activation model in which the first Ca(2+) ion binds usually to site I, thereby preordering both bridging side chains Asp 40 and 43, and partially or fully deprotonating the three coordinating Asp residues. This initial binding event prepares the conformation and protonation state of the remaining site for Ca(2+) binding, enabling the second Ca(2+) ion to bind with higher affinity than the first as required for positive cooperativity.

Probabilities of numbers of ligands proximal to an ion lead to simple, general formulae for the free energy of ion selectivity between different media. That free energy does not depend on the definition of an inner shell for ligand-counting, but other quantities of mechanistic interest do. If analysis is restricted to a specific coordination number, then two distinct probabilities are required to obtain the free energy in addition. The normalizations of those distributions produce partition function formulae for the free energy. Quasi-chemical theory introduces concepts of chemical equilibrium, then seeks the probability that is simplest to estimate, that of the most probable coordination number. Quasi-chemical theory establishes the utility of distributions of ligand-number, and sharpens our understanding of quasi-chemical calculations based on electronic structure methods. This development identifies contributions with clear physical interpretations, and shows that evaluation of those contributions can establish a mechanistic understanding of the selectivity in ion channels.

The hydration of K(+) is studied using a hierarchy of theoretical approaches, including ab initio Born-Oppenheimer molecular dynamics and Car-Parrinello molecular dynamics, a polarizable force field model based on classical Drude oscillators, and a nonpolarizable fixed-charge potential based on the TIP3P water model. While models based more directly on quantum mechanics offer the possibility to account for complex electronic effects, polarizable and fixed-charges force fields allow for simulations of large systems and the calculation of thermodynamic observables with relatively modest computational costs. A particular emphasis is placed on investigating the sensitivity of the polarizable model to reproduce key aspects of aqueous K(+), such as the coordination structure, the bulk hydration free energy, and the self diffusion of K(+). It is generally found that, while the simple functional form of the polarizable Drude model imposes some restrictions on the range of properties that can simultaneously be fitted, the resulting hydration structure for aqueous K(+) agrees well with experiment and with more sophisticated computational models. A counterintuitive result, seen in Car-Parrinello molecular dynamics and in simulations with the Drude polarizable force field, is that the average induced molecular dipole of the water molecules within the first hydration shell around K(+) is slightly smaller than the corresponding value in the bulk. In final analysis, the perspective of K(+) hydration emerging from the various computational models is broadly consistent with experimental data, though at a finer level there remain a number of issues that should be resolved to further our ability in modeling ion hydration accurately.

Transferring Na(+) and K(+) ions from their preferred coordination states in water to states having different coordination numbers incurs a free energy cost. In several examples in nature, however, these ions readily partition from aqueous-phase coordination states into spatial regions having much higher coordination numbers. Here we utilize statistical theory of solutions, quantum chemical simulations, classical mechanics simulations, and structural informatics to understand this aspect of ion partitioning. Our studies lead to the identification of a specific role of the solvation environment in driving transitions in ion coordination structures. Although ion solvation in liquid media is an exergonic reaction overall, we find it is also associated with considerable free energy penalties for extracting ligands from their solvation environments to form coordinated ion complexes. Reducing these penalties increases the stabilities of higher-order coordinations and brings down the energetic cost to partition ions from water into overcoordinated binding sites in biomolecules. These penalties can be lowered via a reduction in direct favorable interactions of the coordinating ligands with all atoms other than the ions themselves. A significant reduction in these penalties can, in fact, also drive up ion coordination preferences. Similarly, an increase in these penalties can lower ion coordination preferences, akin to a Hofmeister effect. Since such structural transitions are effected by the properties of the solvation phase, we anticipate that they will also occur for other ions. The influence of other factors, including ligand density, ligand chemistry, and temperature, on the stabilities of ion coordination structures are also explored.