Search NASA⌕ Search

SEARCH · Search NASA

Results for “Protein aggregation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Strengthening of enterococcal biofilms by Esp

Multidrug-resistant (MDR) Enterococcus faecalis are major causes of hospital-acquired infections. Numerous clinical strains of E . faecalis harbor a large pathogenicity island that encodes enterococcal surface protein (Esp), which is suggested to promote biofilm production and virulence, but this remains controversial. To resolve this issue, we characterized the Esp N-terminal region, the portion implicated in biofilm production. Small angle X-ray scattering indicated that the N-terminal region had a globular head, which consisted of two DEv-Ig domains as visualized by X-ray crystallography, followed by an extended tail. The N-terminal region was not required for biofilm production but instead significantly strengthened biofilms against mechanical or degradative disruption, greatly increasing retention of Enterococcus within biofilms. Biofilm strengthening required low pH, which resulted in Esp unfolding, aggregating, and forming amyloid-like structures. The pH threshold for biofilm strengthening depended on protein stability. A truncated fragment of the first DEv-Ig domain, plausibly generated by a host protease, was the least stable and sufficient to strengthen biofilms at pH ≤ 5.0, while the entire N-terminal region and intact Esp on the enterococcal surface was more stable and required a pH ≤ 4.3. These results suggested a virulence role of Esp in strengthening enterococcal biofilms in acidic abiotic or host environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Effect of the electrostatic surface potential on the oligomerization of full-length human recombinant prion protein at single-molecule level

The electrostatic surface potential (ESP) of prion oligomers has critical influences on the aggregating processes of the prion molecules. The atomic force microscopy (AFM) and structural simulation were combined to investigate the molecular basis of the full-length human recombinant prion oligomerization on mica surfaces. The high resolution non-intrusive AFM images showed that the prion oligomers formed different patterns on mica surfaces at different buffer pH values. The basic binding units for the large oligomers were determined to be prion momoners (Ms), dimers (Ds), and trimers (Ts). The forming of the D and T units happened through the binding of hydrophobic β-sheets of the M units. In contrast, the α-helices of these M, D, and T units were the binding areas for the formation of large oligomers. At pH 4.5, the binding units M, D, and T showed clear polarized ESP distributions on the surface domains, while at pH 7.0, they showed more evenly distributed ESPs. Based on the conformations of oligomers observed from AFM images, the D and T units were more abundantly on mica surface at pH 4.5 because the ESP re-distribution of M units helped to stabilize these larger oligomers. The amino acid side chains involved in the binding interfaces were stabilized by hydrogen bonds and electrostatic interactions. The detailed analysis of the charged side chains at pH 4.5 indicated that the polarized ESPs induced the aggregations among M, D, and T to form larger oligomers. Therefore, the hydrogen bonds and electrostatic interactions worked together to form the stabilized prion oligomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Variable optical properties of light-harvesting complex II revisited

Understanding photosynthetic light harvesting requires knowledge of the molecular mechanisms that dissipate excess energy in thylakoids. However, it remains unclear how the physical environment of light-harvesting complex II (LHCII) influences the process of chlorophyll de-excitation. Here, we demonstrate that protein-protein interactions between LHCIIs affect the optical properties of LHCII and thus influence the total energy budget. Aggregation of LHCII in the dark altered its absorption properties, independent of the amount of prior light exposure. We also revisited the triplet excited state involved in light-induced fluorescence quenching and found another relaxation pathway involving emission in the green region, which might be related to triplet excited energy transfer to neighboring carotenoids and annihilation processes that result in photoluminescence. LHCII- containing liposomes with different protein densities exhibited altered fluorescence and scattering properties. Our results suggest that macromolecular reorganization affects overall optical properties, which need to be addressed to compare the level of energy dissipation.

Iwai, Masakazu↗
DOE OSTI2020

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Convective flow effects on protein crystal growth

A high-resolution microscopic interferometric setup for the monitoring of protein morphologies has been developed. Growth or dissolution of a crystal can be resolved with a long-term depth resolution of 200 A and a lateral resolution of 2 microns. This capability of simultaneously monitoring the interfacial displacement with high local depth resolution has yielded several novel results. We have found with lysozyme that (1) the normal growth rate is oscillatory, and (2) depending on the impurity content of the solution, the growth step density is either greater or lower at the periphery of a facet than in its center. The repartitioning of Na plus and Cl minus ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed, to interpret the large body of data in unified way. The results strongly suggest that (1) the ion to lysozyne ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter electron microscopy results appear to confirm this finding, which could have far-reaching consequences for x-ray diffraction studies. A computational model for diffusive-convective transport in protein crystallization has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies and our kinetics data for the growth of lysozyme. The results show that even in the small cell employed, protein concentration nonuniformities and gravity-driven solutal convection can be significant. The calculated convection velocities are of the same order to magnitude as those found in earlier experiments. As expected, convective transport, i.e., at Og, lysozyme crystal growth remains kinetically limited. The salt distribution in the crystal is predicted to be non-uniform at both 1g and 0g, as a consequence of protein depletion in the solution. Static and dynamic light scattering studies in undersaturated and supersaturated solutions have been performed. Diffusivities in undersaturated solutions, were found to vary with lysozyme concentrations. Depending on the salt concentration, the diffusivities either increase or decrease. Interestingly, the corresponding static scattering intensities behave oppositely, Our current analysis indicates that these changes are inconsistent with aggregation in undersaturated solutions. However, the data are compatible with concentration-dependent changes of the interactions between protein and salt.

Rosenberger, Franz↗
DOE OSTI2021

In-situ Investigations on Gold Nanoparticles Stabilization Mechanisms in Biological Environments Containing HSA

Nanoparticles (NPs) developments advance innovative biomedical applications. However, complex interactions and the low colloidal stability of NPs in biological media restrict their widespread utilization. The influence of NPs properties on the colloidal stability for gold NPs with 5 and 40 nm in diameter with two surface modifications, methoxy-polyethylene glycol-sulfhydryl (PEG) and citrate, in NaCl and human serum albumin (HSA) protein solution, is investigated. This study is based on small-angle X-ray scattering (SAXS) methods allowing the in-situ monitoring of interactions in physiological conditions. The PEG coating provides high colloidal stability for NPs of both sizes. For 5 nm NPs in NaCl solution, a stable 3D self-assembled body-centered cubic (BCC) arrangement is detected with an interparticle distance of 20.7 ± 0.1 nm. In protein solution, this distance increases to 21.9 ± 0.1 nm by protein penetration inside the ordered structure. For citrate-capped NPs, a different mechanism is observed. The protein particles attach to the NPs surfaces, and an appropriate concentration of proteins results in a stable suspension. Cryogenic transmission electron microscopy (Cryo-TEM), UV–visible spectroscopy, and dynamic light scattering (DLS) support the SAXS results. The findings will pave the way to design and synthesize NPs with controlled behaviors in biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Biomarkers for tau pathology

The aggregation of fibrils of hyperphosphorylated and C-terminally truncated microtubule-associated tau protein characterizes 80% of all dementia disorders, the most common neurodegenerative disorders. These so-called tauopathies are hitherto not curable and their diagnosis, especially at early disease stages, has traditionally proven difficult. A keystone in the diagnosis of tauopathies was the development of methods to assess levels of tau protein in vivo in cerebrospinal fluid, which has significantly improved our knowledge about these conditions. Tau proteins have also been measured in blood, but the importance of tau-related changes in blood is still unclear. The recent addition of positron emission tomography ligands to visualize, map and quantify tau pathology has further contributed with information about the temporal and spatial characteristics of tau accumulation in the living brain. Together, the measurement of tau with fluid biomarkers and positron emission tomography constitutes the basis for a highly active field of research. This review describes the current state of biomarkers for tau biomarkers derived from neuroimaging and from the analysis of bodily fluids and their roles in the detection, diagnosis and prognosis of tau-associated neurodegenerative disorders, as well as their associations with neuropathological findings, and aims to provide a perspective on how these biomarkers might be employed prospectively in research and clinical settings.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

A setup for simultaneous measurement of infrared spectra and light scattering signals: Watching amyloid fibrils grow from intact proteins

A setup for the simultaneous measurement of mid-infrared spectra and static light scattering is described that can be used for the analysis of the formation of nanoscale and microscopic aggregates from smaller molecules to biopolymers. It can be easily integrated into sample chambers of infrared spectrometers or combined with laser beams from tunable infrared lasers. Here, its use for the analysis of the formation of amyloid fibrils from intact proteins is demonstrated. The formation of amyloid fibrils or plaques from proteins is a widespread and pathogenetic relevant process, and a number of diseases are caused and correlated with the deposition of amyloid fibrils in cells and tissues. The molecular mechanisms of these transformations, however, are still unclear. We report here the simultaneous measurement of infrared spectra and static light scattering for the analysis of fibril formation from egg-white lysozyme. The transformation of the native form into non-native forms rich in β-sheet structure is measured by analysis of the amide I spectral region in the infrared spectra, which is sensitive for local structures. At the same time, light scattering signals at forward direction as well as the forward/backward ratio, which are sensitive for the number of scattering centers and their approximate sizes, respectively, are collected for the analysis of fibril growth. Thermodynamic and kinetic parameters as well as mechanistic information are deduced from the combination of the two complementary techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Mycobacterium tuberculosis SatS is a chaperone for the SecA2 protein export pathway

The SecA2 protein export system is critical for the virulence of Mycobacterium tuberculosis. However, the mechanism of this export pathway remains unclear. Through a screen for suppressors of a secA2 mutant, we identified a new player in the mycobacterial SecA2 pathway that we named SatS for SecA2 (two) Suppressor. In M. tuberculosis, SatS is required for the export of a subset of SecA2 substrates and for growth in macrophages. We further identify a role for SatS as a protein export chaperone. SatS exhibits multiple properties of a chaperone, including the ability to bind to and protect substrates from aggregation. Our structural studies of SatS reveal a distinct combination of a new fold and hydrophobic grooves resembling preprotein-binding sites of the SecB chaperone. These results are significant in better defining a molecular pathway for M. tuberculosis pathogenesis and in expanding our appreciation of the diversity among chaperones and protein export systems.

Miller, Brittany K.↗
DOE OSTI2015

Molecular Details of Olfactomedin Domains Provide Pathway to Structure-Function Studies

Olfactomedin (OLF) domains are found within extracellular, multidomain proteins in numerous tissues of multicellular organisms. Even though these proteins have been implicated in human disorders ranging from cancers to attention deficit disorder to glaucoma, little is known about their structure(s) and function(s). Here we biophysically, biochemically, and structurally characterize OLF domains from H. sapiens olfactomedin-1 (npoh-OLF, also called noelin, pancortin, OLFM1, and hOlfA), and M. musculus gliomedin (glio-OLF, also called collomin, collmin, and CRG-L2), and compare them with available structures of myocilin (myoc-OLF) recently reported by us and R. norvegicus glio-OLF and M. musculus latrophilin-3 (lat3-OLF) by others. Although the five-bladed β-propeller architecture remains unchanged, numerous physicochemical characteristics differ among these OLF domains. First, npoh-OLF and glio-OLF exhibit prominent, yet distinct, positive surface charges and copurify with polynucleotides. Second, whereas npoh-OLF and myoc-OLF exhibit thermal stabilities typical of human proteins near 55°C, and most myoc-OLF variants are destabilized and highly prone to aggregation, glio-OLF is nearly 20°C more stable and significantly more resistant to chemical denaturation. Phylogenetically, glio-OLF is most similar to primitive OLFs, and structurally, glio-OLF is missing distinguishing features seen in OLFs such as the disulfide bond formed by N- and C- terminal cysteines, the sequestered Ca 2+ ion within the propeller central hydrophilic cavity, and a key loop-stabilizing cation-π interaction on the top face of npoh-OLF and myoc-OLF. While deciphering the explicit biological functions, ligands, and binding partners for OLF domains will likely continue to be a challenging long-term experimental pursuit, we used structural insights gained here to generate a new antibody selective for myoc-OLF over npoh-OLF and glio-OLF as a first step in overcoming the impasse in detailed functional characterization of these biomedically important protein domains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Structural and functional characterization of a missense mutant of human γS-crystallin associated with dominant infantile cataracts

Highlights: • We present structural analysis of γS-Y67N, a cataract-active mutant of human γS-crystallin. • Mutation alters the energetics and kinetics of γS-Y67N underlying structural changes. • Aggregation is enhanced and Ca{sup 2+} binding is compromised as a result of mutation. • We explore the structure-function relationship in human γS-crystallin in health and disease. Infantile cataracts constitute one of the most important causes of childhood blindness worldwide. Human γS-crystallin is the most abundant protein in the eye lens cortex. A missense mutant of human γS-crystallin, Y67N (abbreviated hereafter as γS-Y67N) is recently reported to be associated with dominant infantile cataracts. To understand the structural basis for γS-Y67N to cause lens opacification, we constructed, expressed and purified γS-Y67N and its wild-type (abbreviated hereafter as γS-WT) and studied the structural and functional differences between them in solution using circular dichroism (CD), differential scanning calorimetry (DSC), fluorescence spectroscopy and extrinsic spectral probes. Extensive equilibrium characterization indicate that replacement of the highly conserved Tyr at 67th position by Asn distorts the conserved Tyr corner at the second Greek key motif in the N-terminal domain (NTD) and leads to substantial loss of structural stability. Our intrinsic fluorescence quenching results reveal differential in-vitro refolding kinetics identifying partially folded kinetic intermediates for both proteins. Extrinsic fluorescence studies further reveal loosening up of the compact structure of γS-crystallin upon mutation associated with enhanced aggregation. As Ca{sup 2+} homeostasis is a crucial regulator of lens transparency, we further investigated the Ca{sup 2+}-binding properties of γS-WT and γS-Y67N by isothermal titration calorimetry (ITC) to identify lens Ca{sup 2+} distribution in health and in disease. Overall, our results highlight the vital role of conserved Tyr corners in stabilizing Greek key motifs and provide useful structural and functional insights into the mechanism of cataract formation in humans.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Differential Misfolding Properties of Glaucoma-Associated Olfactomedin Domains from Humans and Mice

Mutations in myocilin, predominantly within its olfactomedin (OLF) domain, are causative for the heritable form of open angle glaucoma in humans. Surprisingly, mice expressing Tyr423His mutant myocilin, corresponding to a severe glaucoma-causing mutation (Tyr437His) in human subjects, exhibit a weak, if any, glaucoma phenotype. To address possible protein-level discrepancies between mouse and human OLFs, which might lead to this outcome, biophysical properties of mouse OLF were characterized for comparison with those of human OLF. Here, the 1.55 Å resolution crystal structure of mouse OLF reveals an asymmetric 5-bladed β-propeller that is nearly indistinguishable from previous structures of human OLF. Wild-type and selected mutant mouse OLFs mirror thermal stabilities of their human OLF counterparts, including characteristic stabilization in the presence of calcium. Mouse OLF forms thioflavin T-positive aggregates with a similar end-point morphology as human OLF, but amyloid aggregation kinetic rates of mouse OLF are faster than human OLF. Simulations and experiments support the interpretation that kinetics of mouse OLF are faster because of a decreased charge repulsion arising from more neutral surface electrostatics. Taken together, phenotypic differences observed in mouse and human studies of mutant myocilin could be a function of aggregation kinetics rates, which would alter the lifetime of putatively toxic protofibrillar intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The denatured state of HIV ‐1 protease under native conditions

Abstract The denatured state of several proteins has been shown to display transient structures that are relevant for folding, stability, and aggregation. To detect them by nuclear magnetic resonance (NMR) spectroscopy, the denatured state must be stabilized by chemical agents or changes in temperature. This makes the environment different from that experienced in biologically relevant processes. Using high‐resolution heteronuclear NMR spectroscopy, we have characterized several denatured states of a monomeric variant of HIV‐1 protease, which is natively structured in water, induced by different concentrations of urea, guanidinium chloride, and acetic acid. We have extrapolated the chemical shifts and the relaxation parameters to the denaturant‐free denatured state at native conditions, showing that they converge to the same values. Subsequently, we characterized the conformational properties of this biologically relevant denatured state under native conditions by advanced molecular dynamics simulations and validated the results by comparison to experimental data. We show that the denatured state of HIV‐1 protease under native conditions displays rich patterns of transient native and non‐native structures, which could be of relevance to its guidance through a complex folding process.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

High-speed atomic force microscopy reveals structural dynamics of α-synuclein monomers and dimers

We report that α-Synuclein (α-syn) is the major component of the intraneuronal inclusions called Lewy bodies, which are the pathological hallmark of Parkinson’s disease. α-Syn is capable of self-assembly into many different species, such as soluble oligomers and fibrils. Even though attempts to resolve the structures of the protein have been made, detailed understanding about the structures and their relationship with the different aggregation steps is lacking, which is of interest to provide insights into the pathogenic mechanism of Parkinson’s disease. Here we report the structural flexibility of α-syn monomers and dimers in an aqueous solution environment as probed by single-molecule time-lapse high-speed AFM. In addition, we present the molecular basis for the structural transitions using discrete molecular dynamics (DMD) simulations. α-Syn monomers assume a globular conformation, which is capable of forming tail-like protrusions over dozens of seconds. Importantly, a globular monomer can adopt fully extended conformations. Dimers, on the other hand, are less dynamic and show a dumbbell conformation that experiences morphological changes over time. DMD simulations revealed that the α-syn monomer consists of several tightly packed small helices. The tail-like protrusions are also helical with a small β-sheet, acting as a “hinge”. Monomers within dimers have a large interfacial interaction area and are stabilized by interactions in the non-amyloid central (NAC) regions. Additionally, the dimer NAC-region of each α-syn monomer forms a β-rich segment. Moreover, NAC-regions are located in the hydrophobic core of the dimer.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Leveraging Topological Events in Tracking Graphs for Understanding Particle Diffusion

Abstract Single particle tracking (SPT) of fluorescent molecules provides significant insights into the diffusion and relative motion of tagged proteins and other structures of interest in biology. However, despite the latest advances in high‐resolution microscopy, individual particles are typically not distinguished from clusters of particles. This lack of resolution obscures potential evidence for how merging and splitting of particles affect their diffusion and any implications on the biological environment. The particle tracks are typically decomposed into individual segments at observed merge and split events, and analysis is performed without knowing the true count of particles in the resulting segments. Here, we address the challenges in analyzing particle tracks in the context of cancer biology. In particular, we study the tracks of KRAS protein, which is implicated in nearly 20% of all human cancers, and whose clustering and aggregation have been linked to the signaling pathway leading to uncontrolled cell growth. We present a new analysis approach for particle tracks by representing them as tracking graphs and using topological events – merging and splitting, to disambiguate the tracks. Using this analysis, we infer a lower bound on the count of particles as they cluster and create conditional distributions of diffusion speeds before and after merge and split events. Using thousands of time‐steps of simulated and in‐vitro SPT data, we demonstrate the efficacy of our method, as it offers the biologists a new, detailed look into the relationship between KRAS clustering and diffusion speeds.

McDonald, T.↗
NASA NTRS1996

Lack of Evidence for Prenucleation Aggregate Formation in Lysozyme Crystal Growth Solutions

There have been numerous claims of large concentrations of prenucleation aggregates in supersaturated as well as undersaturated lysozyme solutions at high salt concentrations. The presence of these aggregates was derived from measurements of the light or neutron scattering intensity, ultracentrifugation and dialysis behavior, as well as over-simplified crystal growth kinetics considerations. In all these interpretations it has been assumed that lysozyme solutions are either ideal or that protein interactions are independent of salt concentration. Contrary to these presumptions, our static and dynamic light scattering experiments provide evidence that lysozyme forms highly non-ideal, strongly interacting solutions. At low salt concentrations, the scattering intensities fall well below the values expected for an ideal, monomeric solution at the same protein concentration, while diffusivities increase with increasing protein concentration. Upon increase in salt concentration, these trends are eventually reversed. This enhancement in scattering intensity and decrease in diffusivity was widely interpreted as sign of aggregate formation. Yet, a quantitative interpretation of the scattering behavior over the whole salt concentration range can only be given in terms of a transition from net repulsion to net attraction between lysozyme monomers. Increased salt screening of the electrostatic repulsion among the protein macro-ions, together with attractive protein interactions, such as van der Waals, hydrophobic and hydration forces, provide an unambiguous mechanism for the observed transition and a more physical interpretation of the various observations.

Muschol, Martin↗
DOE OSTI2020

Enhanced purification coupled with biophysical analyses shows cross-β structure as a core building block for Streptococcus mutans functional amyloids

Streptococcus mutans is an etiologic agent of human dental caries that forms dental plaque biofilms containing functional amyloids. Three amyloidogenic proteins, P1, WapA, and Smu_63c were previously identified. C123 and AgA are naturally occurring amyloid-forming fragments of P1 and WapA, respectively. We determined that four amyloidophilic dyes, ThT, CDy11, BD-oligo, and MK-H4, differentiate C123, AgA, and Smu_63c amyloid from monomers, but non-specific binding to bacterial cells in the absence of amyloid precludes their utility for identifying amyloid in biofilms. Congo redinduced birefringence is a more specific indicator of amyloid formation and differentiates biofilms formed by wild-type S. mutans from a triple ΔP1/WapA/Smu_63c mutant with reduced bioiflm forming capabilities. Amyloid accumulation is a late event, appearing in older S. mutans biofilms after 60hours of growth. Amyloid derived from pure preparations of all three proteins is visualized by electron microscopy as mat-like structures. Typical amyloid fibers become evident following protease digestion to eliminate non-specific aggregates and monomers. Amyloid mats, similar in appearance to those reported in S. mutans biofilm extracellular matrices, are reconstituted by co-incubation of monomers and amyloid fibers. X-ray fiber diffraction of amyloid mats and fibers from all three proteins demonstrate patterns reflective of a cross-β amyloid structure.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field

We have implemented the Martini force field within Lawrence Livermore National Laboratory’s molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein–lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗