Testing the inter-operability of the CHARMM and SPC/Fw force fields for conformational sampling

Effects of the structure of lipid-based agents in their complexation with a single stranded mRNA fragment: a computational study

In this work we employed fully atomistic molecular dynamics simulations, aiming towards a better understanding of the mechanisms associated with the formation and the stability of lipid-based RNA nanoassemblies, in an aqueous environment. We examined two groups of lipid-based complexation agents, differing in the degree of hydrophobicity and in the overall charge. The first group was comprised of cationic ionizable agents while the second included electrically neutral amphoteric phosphatidylcholine lipids. It was found that the overall charge of the complexation agents played the most decisive role in the energetics of the lipid/RNA association, while their degree of hydrophobicity affected their self-assembly and their complexation kinetics. The latter also affected the structural stability of the formed complexes since the water entrapped within the clusters of the less hydrophobic agents appeared to reduce the coherence of the lipid-RNA nanoassemblies. The combined effects of the aforementioned attributes dictated also the RNA conformation after complexation. The results from the present study provide thus new insight towards controlling the morphology, the energetic stability and the structural integrity of the formed complexes.

Impact of Glutamate Carboxylation in the Adsorption of the α-1 Domain of Osteocalcin to Hydroxyapatite and Titania

One proposed mechanism of implant fouling is attributed to the nonspecific adsorption of non-collagenous bone matrix proteins (NCPs) onto a newly implanted interface. With the goal of capturing the fundamental mechanistic and thermodynamic forces that govern changes in these NCP recognition domains as a function of γ-carboxyglutamic acid (Gla) post-translational modification and surface chemistry, we probe the adsorption process of the most commonly occurring NCP, osteocalcin, onto a mineral and metal oxide surface. Here, we apply two enhanced sampling methods to independently probe the effects of post-translational modification and peptide structure on adsorption. First, well-tempered metadynamics was used to capture the binding of acetyl and N-methylamide capped glutamic acid and Gla single amino acids onto crystalline hydroxyapatite and titania model surfaces at physiological pH. Following this, parallel tempering metadynamics in the well-tempered ensemble (PTMetaD-WTE) was used to study adsorption of the α-1 domain of osteocalcin onto hydroxyapatite and titania. Simulations were performed for the α-1 domain of osteocalcin in both its fully decarboxylated (dOC) and fully carboxylated (OC) form. Our simulations find that increased charge density due to carboxylation results in increased interactions at the interface, and stronger adsorption of the single amino acids to both surfaces. Interestingly, the role of Gla in promoting compact and helical structure in the α-1 domain resulted in disparate binding modes at the two surfaces, which is attributed to differences in interfacial water behavior. Overall, this work provides a benchmark for understanding the mechanisms that drive adsorption of Gla-containing mineralizing proteins onto different surface chemistries.

Caffeine destabilizes preformed Aβ protofilaments: insights from all atom molecular dynamics simulations

The aggregation and deposition of neurotoxic Aβ fibrils are key in the etiology of Alzheimer's disease (AD). It has been clinically recognized as a major form of dementia across the globe. Finding and testing various natural compounds to target Aβ fibrils to disrupt their stable structures seems to be a promising and attractive therapeutic strategy. The destabilization effects of caffeine on Aβ fibrils are investigated via in silico studies, where a series of molecular dynamics (MD) simulations, each of 100 ns, was conducted. The simulation outcomes obtained henceforth clearly indicated the drift of the terminal chains from the protofibrils, leading to disorganization of the characteristically organized cross-β structures of Aβ fibrils. The structural instability of Aβ17-42 protofibrils is explained through enhanced fluctuations in the RMSD, radius of gyration and RMSF values in the presence of caffeine. The key interactions providing stability, comprising D23-K28 salt bridges, intra- and inter-chain hydrogen bonding and hydrophobic interactions involving interchain A21-V36 and F19-G38 and intrachain L34-V36, were found to be disrupted due to increases in the distances between the participating components. The loss of β-sheet structure with the introduction of turns and α-helices in terminal chains may further inhibit the formation of higher order aggregates, which is necessary to stop the progression of the disease. The atomistic details obtained via MD studies relating to the mechanism behind the underlying destabilization of Aβ17-42 protofibrils by caffeine encourage further investigations exploring the potency of natural compounds to treat AD via disrupting preformed neurotoxic Aβ protofibrils.

Pathways to Structure-Property Relationships of Peptide-Materials Interfaces: Challenges in Predicting Molecular Structures

An in-depth appreciation of how to manipulate the molecular-level recognition between peptides and aqueous materials interfaces, including nanoparticles, will advance technologies based on self-organized metamaterials for photonics and plasmonics, biosensing, catalysis, energy generation and harvesting, and nanomedicine. Exploitation of the materials-selective binding of biomolecules is pivotal to success in these areas and may be particularly key to producing new hierarchically structured biobased materials. These applications could be accomplished by realizing preferential adsorption of a given biomolecule onto one materials composition over another, one surface facet over another, or one crystalline polymorph over another. Deeper knowledge of the aqueous abiotic-biotic interface, to establish clear structure-property relationships in these systems, is needed to meet this goal. In particular, a thorough structural characterization of the surface-adsorbed peptides is essential for establishing these relationships but can often be challenging to accomplish via experimental approaches alone. In addition to myriad existing challenges associated with determining the detailed molecular structure of any molecule adsorbed at an aqueous interface, experimental characterization of materials-binding peptides brings new, complex challenges because many materials-binding peptides are thought to be intrinsically disordered. This means that these peptides are not amenable to experimental techniques that rely on the presence of well-defined secondary structure in the peptide when in the adsorbed state. To address this challenge, and in partnership with experiment, molecular simulations at the atomistic level can bring complementary and critical insights into the origins of this abiotic/biotic recognition and suggest routes for manipulating this phenomenon to realize new types of hybrid materials. For the reasons outlined above, molecular simulation approaches also face challenges in their successful application to model the biotic-abiotic interface, related to several factors. For instance, simulations require a plausible description of the chemistry and the physics of the interface, which comprises two very different states of matter, in the presence of liquid water. Also, it is essential that the conformational ensemble be comprehensively characterized under these conditions; this is especially challenging because intrinsically disordered peptides do not typically admit one single structure or set of structures. Moreover, a plausible structural model of the substrate is required, which may require a high level of detail, even for single-element materials such as Au surfaces or graphene. Developing and applying strategies to make credible predictions of the conformational ensemble of adsorbed peptides and using these to construct structure-property relationships of these interfaces have been the goals of our efforts. We have made substantial progress in developing interatomic potentials for these interfaces and adapting advanced conformational sampling approaches for these purposes. This Account summarizes our progress in the development and deployment of interfacial force fields and molecular simulation techniques for the purpose of elucidating these insights at biomolecule-materials interfaces, using examples from our laboratories ranging from noble-metal interfaces to graphitic substrates (including carbon nanotubes and graphene) and oxide materials (such as titania). In addition to the well-established application areas of plasmonic materials, biosensing, and the production of medical implant materials, we outline new directions for this field that have the potential to bring new advances in areas such as energy materials and regenerative medicine.

Facet-Specific Adsorption of Tripeptides at Aqueous Au Interfaces: Open Questions in Reconciling Experiment and Simulation

The adsorption of three homo-tripeptides, HHH, YYY, and SSS, at the aqueous Au interface is investigated, using molecular dynamics simulations. We find that consideration of surface facet effects, relevant to experimental conditions, opens up new questions regarding interpretations of current experimental findings. Our well-tempered metadynamics simulations predict the rank ordering of the tripeptide binding affinities at aqueous Au(111) to be YYY > HHH > SSS. This ranking differs with that obtained from existing experimental data which used surface-immobilized Au nanoparticles as the target substrate. The influence of Au facet on these experimental findings is then considered, via our binding strength predictions of the relevant amino acids at aqueous Au(111) and Au(100)(1 × 1). The Au(111) interface supports an amino acid ranking of Tyr > HisA ≃ HisH > Ser, matching that of the tripeptides on Au(111), while the ranking on Au(100) is HisA > Ser ≃ Tyr ≃ HisH, with only HisA showing non-negligible binding. The substantial reduction in Tyr amino acid affinity for Au(100) vs Au(111) offers one possible explanation for the experimentally observed weaker adsorption of YYY on the nanoparticle-immobilized substrate compared with HHH. In a separate set of simulations, we predict the structures of the adsorbed tripeptides at the two aqueous Au facets, revealing facet-dependent differences in the adsorbed conformations. Our findings suggest that Au facet effects, where relevant, may influence the adsorption structures and energetics of biomolecules, highlighting the possible influence of the structural model used to interpret experimental binding data.