Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Molecular Dynamic simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

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

At least 1,549 records · Page 86Linked to original sources
PubMed1998

Molecular dynamics simulation of a 13-mer duplex DNA: a PvuII substrate.

Parallel version of AMBER 4.1 was ported and optimised on the Indian parallel supercomputer PARAM OpenFrame built around Sun Ultra Sparc processors. This version of AMBER program was then used to carry out molecular dynamics (MD) simulations on 5'-TGACCAGCTGGTC-3', a substrate for PvuII enzyme. MD simulations in water are carried out under following conditions: (i) unconstrained at 300 K (230 ps); (ii) unconstrained at 283 K (500 ps); (iii) Watson-Crick basepair constrained at 283 K (1 ns); and (iv) Watson-Crick basepair constrained with ions at 283 K (1.2 ns). In all these simulation studies, the molecule was observed to be bending and maximum distortions in the double helix around was seen around the G7:C7' basepair, which is the phosphodiester bond that is cleaved by PvuII. Analysis of MD simulation with ions carried out for 1.2 ns also pointed out that the conformation of double helix alternates between a conformation close to B-form and close to A-form. It is argued that a bent non-standard conformation is recognised by the PvuII enzyme. The maximum bend occurs at the G7:C7' region, weakening the phosphodiester bond and allows His48 to get placed in such a fashion to permit the scission through a general base mechanism. The bending and distortion observed is a property of the sequence which acts as a substrate for PvuII enzyme. This is confirmed by carrying out MD studies on the Dickerson's sequence d(CGCGAATTCGCG)2 as a reference molecule, which practically does not bend or get deformed.

Base Composition↗
PubMed2005

Molecular dynamics simulations of the oxidation of aluminum nanoparticles.

The oxidation of aluminum nanoparticles is studied with classical molecular dynamics and the Streitz-Mintmire (Streitz, F. H.; Mintmire, J. W. Phys. Rev. B 1994, 50, 11996) electrostatic plus (ES+) potential that allows for the variation of electrostatic charge on all atoms in the simulation. The structure and charge distributions of bulk crystalline alpha-Al(2)O(3), a surface slab of alpha-Al(2)O(3) with an exposed (0001) basal plane, and an isolated Al(2)O(3) nanoparticle are studied. Constant NVT simulations of the oxidation of aluminum nanoparticles are also performed with different oxygen exposures. The calculations simulate a thermostated one-time exposure of an aluminum nanoparticle to different numbers of surface oxygen atoms. In the first set of oxidation studies, the overall approximate ratios of Al to O in the nanoparticle are 1:1 and 2:1. The nanoparticles are annealed to 3000 K and are then cooled to 500, 1000, or 1500 K. The atomic kinetic energy is scaled during the simulation to maintain the desired temperature. The structure and charge distributions in the oxidized nanoparticles differ from each other and from those of the bulk Al(2)O(3) phases. In the Al(1)O(1) simulation, an oxide shell forms that stabilizes the shape of the particle, and thus the original structure of the nanoparticle is approximately retained. In the case of Al(1)O(0.5), there is insufficient oxygen to form a complete oxide shell, and the oxidation results in particles of irregular shapes and rough surfaces. The particle surface is rough, and the nanoparticle is deformed.

Aluminum↗
PubMed2002

A large scale molecular dynamics simulation code using the fast multipole algorithm (FMD): performance and application.

We present the performance of the fast classical molecular dynamics (MD) code, fast molecular dynamics (FMD), designed for efficient, object-oriented, and scalable large scale simulations, and summarize its application to a liquid crystalline cluster. FMD uses an implementation of the three-dimensional fast multipole method, developed in our group. The fast multipole method offers an efficient way (order O(N)) to handle long range electrostatic interactions, thus, enabling more realistic simulations of large molecular systems. Performance testing was carried out on IBM SP2, SGI Origin 2000, and CRAY T3E massively parallel systems using the MPI massage passing library. The electrostatic forces were tested on models of up to 100,000 randomly placed charges, and on protein and liquid crystalline molecular systems of over 99,000 atoms. Tests on the stability of the method are presented, along with comparisons with direct calculations, the NAMD2 code, and the physical multipole-based cell-multipole method.

Algorithms↗
PubMed2002

Explicit solvent molecular dynamics simulation of duplex formed by the modified oligonucleotide with alternating phosphate/phosphonate internucleoside linkages and its natural counterpart.

Impact of the internucleoside linkage modification by inserting a methylene group on the ability of the modified oligonucleotide to hybridize with a natural DNA strand was studied by fully solvated molecular dynamics (MD) simulations. Three undecamer complexes were analyzed: natural dT(11).dA(11) duplex as a reference and two its analogs with alternating modified and natural linkages in the deoxyadenosine chain. The isopolar, non-isosteric modified linkages were of 5'-O-PO(2)-CH(2)-O-3' (5'PC3') or 5'-O-CH(2)-PO(2)-O-3' (5'CP3') type. Simulations were performed by using the AMBER 5.0 software package with the force field completed by a set of parameters needed to model the modified segments. Both modifications were found to lead to double helical complexes, in which the thymidine strand as well as deoxyriboses and unmodified linkages in the adenosine strand adopted conformations typical for the B-type structure. For each of the two conformational richer modified linkages two stable conformations were found at 300 K: the -ggt and ggt for the 5'PC3' and ggg, tgg for the 5'CP3', respectively. Both modified chains adopted helical conformations with heightened values of the inclination parameter but without affecting the Watson-Crick hydrogen bonds.

Computer Simulation↗
PubMed2003

Characteristic domain motion in the ribosome recycling factor revealed by 15N NMR relaxation experiments and molecular dynamics simulations.

The backbone dynamics of ribosome recycling factor (RRF) from Escherichia coli in water were characterized by (15)N NMR relaxation analysis and molecular dynamics (MD) simulation. RRF is composed of two domains connected by a joint region that consists of two peptide chains, such that the overall structure seems to mimic that of tRNA. MD trajectories indicated that the relative orientation of domains varies on the nanosecond time scale. We analyzed the observed (15)N T(1), T(2), and NOE using an extended model-free spectral density function in which the domain motions with a nanosecond time scale were considered. At 30 degrees C, the order parameters of slow motion () were determined to be approximately 0.9 for domain I and 0.7 for domain II, respectively. These values indicate that domain I is nearly fixed on the molecular diffusion frame, and domain II is wobbling in a cone for which the semi-angle is about 30 degrees.

Escherichia coli↗
PubMed2000

Dynamics of proteins predicted by molecular dynamics simulations and analytical approaches: application to alpha-amylase inhibitor.

The dynamics of alpha-amylase inhibitors has been investigated using molecular dynamics (MD) simulations and two analytical approaches, the Gaussian network model (GNM) and anisotropic network model (ANM). MD simulations use a full atomic approach with empirical force fields, while the analytical approaches are based on a coarse-grained single-site-per-residue model with a single-parameter harmonic potential between sufficiently close (r </= 7 A) residue pairs. The major difference between the GNM and the ANM is that no directional preferences can be obtained in the GNM, all residue fluctuations being theoretically isotropic, while ANM does incorporate directional preferences. The dominant modes of motions are identified by (i) the singular value decomposition (SVD) of the MD trajectory matrices, and (ii) the similarity transformation of the Kirchhoff matrices of inter-residue contacts in the GNM or ANM. The mean-square fluctuations of individual residues and the cross-correlations between domain movements retain the same characteristics, in all approaches-although the dispersion of modes and detailed amplitudes of motion obtained in the ANM conform more closely with MD results. The major weakness of the analytical approaches appears, on the other hand, to be their inadequacy to account for the anharmonic motions or multimeric transitions driven by the slowest collective mode observed in MD. Such motions usually suffer, however, from MD sampling inefficiencies, and multiple independent runs should be tested before making conclusions about their validity and detailed mechanisms. Overall this study invites attention to (i) the robustness of the average properties (mean-square fluctuations, cross-correlations) controlled by the low frequency motions, which are invariably reproduced in all approaches, and (ii) the utility and efficiency of the ANM, the computational time cost of which is of the order of "minutes" (real time), as opposed to "days" for MD simulations. Proteins 2000;40:512-524.

Computer Simulation↗
PubMed2004

How environment supports a state: molecular dynamics simulations of two states in bacteriorhodopsin suggest lipid and water compensation.

The light-driven proton pump bacteriorhodopsin (bR) is a transmembrane protein that uses large conformational changes for proton transfer from the cytoplasmic to the extracellular regions. Crystal structures, due to their solvent conditions, do not resolve the effect of lipid molecules on these protein conformational changes. To begin to understand the molecular details behind such large conformational changes, we simulated two conformations of wild-type bacteriorhodopsin, one of the dark-adapted state and the second of an intermediate (M(O)) state, each within an explicit dimyristoyl-phosphatidylcholine (DMPC) lipid bilayer. The simulations included all-hydrogen and all-atom representations of protein, lipid, and water and were performed for 20 ns. We investigate the equilibrium properties and the dynamic motions of the two conformations in the lipid setting. We note that the conformational state of the M(O) intermediate bR remains markedly different from the dark-adapted bR state in that the M(O) intermediate shows rearrangement of the cytoplasmic portions of helices C, F, and G, and nearby loops. This difference in the states remained throughout the simulations, and the results are stable on the molecular dynamics timescale and provide an illustration of the changes in both lipid and water that help to stabilize a particular state. Our analysis focuses on how the environment adjusts to these two states and on how the dynamics of the helices, loops, and water molecules can be related to the pump mechanism of bacteriorhodopsin. For example, water generally behaves in the same manner on the extracellular sides of both simulations but is decreased in the cytoplasmic region of the M(O) intermediate. We suspect that the different water behavior is closely related to the fluctuations of microcavities volume in the protein interior, which is strongly coupled to the collective motion of the protein. Our simulation result suggests that experimental observation can be useful to verify a decreased number of waters in the cytoplasmic regions of the late-intermediate stages by measuring the rate of water exchange with the interior of the protein.

Bacteriorhodopsins↗
PubMed2004

Molecular dynamics simulation of solvent-polymer interdiffusion: Fickian diffusion.

The interdiffusion of a solvent into a polymer melt has been studied using large scale molecular dynamics and Monte Carlo simulation techniques. The solvent concentration profile and weight gain by the polymer have been measured as a function of time. The weight gain is found to scale as t(1/2), which is expected for Fickian diffusion. The concentration profiles are fit very well assuming Fick's second law with a constant diffusivity. The diffusivity found from fitting Fick's second law is found to be independent of time and equal to the self-diffusion constant in the dilute solvent limit. We separately calculated the diffusivity as a function of concentration using the Darken equation and found that the diffusivity is essentially constant for the concentration range relevant for interdiffusion.

Journal Article↗
PubMed2026

Investigating the mechanisms of PhIP-induced colorectal cancer through network toxicology, machine learning, and molecular dynamics simulation.

BACKGROUND: Over the past few years, 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP)- a compound from grilled or processed meats-has emerged as a major player in cancer development, especially colorectal cancer (CRC). This work dives into its potential links to CRC and uncovers the key genes that bridge this connection. METHODS: We tapped into various databases to pinpoint target genes tied to PhIP and CRC, then ran protein-protein interaction (PPI) analyses for visualization. Next, we explored underlying mechanisms through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. To nail down predictions, we tested 107 machine learning pipelines and picked the best one, validating its accuracy and the core genes' prognostic value across datasets. Next, molecular docking and dynamics simulations probed the interactions between these genes and PhIP. Finally, cell proliferation was assessed using Cell Counting Kit-8 (CCK-8) and 5-ethynyl-2'-deoxyuridine (EdU) assays, and polymerase chain reaction (PCR) was performed to validate the expression levels of the hub genes. RESULTS: Our analysis identified 39 overlapping genes, from which a machine learning model (glmBoost + Enet) identified six candidate targets: CDK4, CEBPB, COMT, SOX9, TIMP1, and TOP2A. To prioritize these, a hierarchical screening framework was applied. Molecular docking and dynamics simulations identified CDK4, COMT, and TIMP1 as the most stable interactors with PhIP. Functional assays confirmed that PhIP treatment significantly enhanced the proliferation of CRC cells. Crucially, quantitative PCR (qPCR) validation in multiple CRC cell lines identified TIMP1 as the primary target, showing the most consistent and significant upregulation upon PhIP exposure. CONCLUSIONS: In essence, these genes drive PhIP is role in CRC, offering novel insights into its molecular pathways. This could reshape how we tackle food-related pollutants, paving the way for better prevention and targeted therapies.

Colorectal cancer (CRC)↗
PubMed2004

Nucleotide-dependent conformational changes in HisP: molecular dynamics simulations of an ABC transporter nucleotide-binding domain.

ATP-binding cassette (ABC) transporters mediate the movement of molecules across cell membranes in both prokaryotes and eukaryotes. In ABC transporters, solute translocation occurs after ATP is either bound or hydrolyzed at the intracellular nucleotide-binding domains (NBDs). Molecular dynamics (MD) simulations have been employed to study the interactions of nucleotide with NBD. The results of extended (approximately 20 ns) MD simulations of HisP (total simulation time approximately 80 ns), the NBD of the histidine transporter HisQMP2J from Salmonella typhimurium, are presented. Analysis of the MD trajectories reveals conformational changes within HisP that are dependent on the presence of ATP in the binding pocket of the protein, and are sensitive to the presence/absence of Mg ions bound to the ATP. These changes are predominantly confined to the alpha-helical subdomain of HisP. Specifically there is a rotation of three alpha-helices within the subdomain, and a movement of the signature sequence toward the bound nucleotide. In addition, there is considerable conformational flexibility in a conserved glutamine-containing loop, which is situated at the interface between the alpha-helical subdomain and the F1-like subdomain. These results support the mechanism for ATP-induced conformational transitions derived from the crystal structures of other NBDs.

ATP-Binding Cassette Transporters↗
PubMed1991

Solution structure of [d(GTATATAC)]2 via restrained molecular dynamics simulations with nuclear magnetic resonance constraints derived from relaxation matrix analysis of two-dimensional nuclear Overhauser effect experiments.

Two-dimensional nuclear Overhauser effect (2D NOE) spectra have been used as the experimental basis for determining the solution structure of the duplex [d(GTATATAC)]2 employing restrained molecular dynamics (rMD) simulations. The MARDIGRAS algorithm has been employed to construct a set of 233 interproton distance constraints via iterative complete relaxation matrix analysis utilizing the peak intensities from the 2D NOE spectra obtained for different mixing times and model structures. The upper and lower bounds for each of the constraints, defining size of a flat-well potential function term used in the rMD simulations, were conservatively chosen as the largest or smallest value calculated by MARDIGRAS. Three different starting models were utilized in several rMD calculations: energy-minimized A-DNA, B-DNA, and a structure containing wrinkled D-DNA in the interior. Considerable effort was made to define the appropriate force constants to be employed with the NOE terms in the AMBER force field, using as criteria the average constraints deviation, the constraints violation energy and the total energy. Of the 233 constraints, one was generated indirectly, but proved to be crucial in defining the structure: the cross-strand A5-H2 A5-H2 distance. As those two protons resonate isochronously for the self-complementary duplex, the distance cannot be determined directly. However, the general pattern of 2D NOE peak intensities, spin-lattice relaxation time (T1) values, and 31P nuclear magnetic resonance spectra lead to use of the A3-H2 A7-H2 distance for A5-H2 A5-H2 as well. Five rMD runs, with different random number seeds, were made for each of the three starting structures with the full distance constraint set. The average structure from all 15 runs and the five-structure averages from each starting structure were all quite similar. Two rMD runs for each starting structure were made with the A5-H2 A5-H2 constraint missing. The average of these six rMD runs revealed differences in structure, compared to that with the full set of constraints, primarily for the middle two base-pairs involving the missing cross-strand constraint but global deviations also were found. Conformational analysis of the resulting structures revealed that the inner four to six base-pairs differed in structure from the termini. Furthermore, an alternating structure was suggested with features alternating for the A-T and T-A steps.

Algorithms↗
PubMed2004

Thermophilicity of wild type and mutant cold shock proteins by molecular dynamics simulation.

In the attempt to clarify possible mechanisms underlying thermal stability of proteins, we study through molecular dynamics thermophile Bc-Csp, mesophile Bs-CspB, and selected mutants. These proteins have been extensively characterized experimentally; researchers showed that differential thermostability among the wild type proteins is fundamentally linked to one or two mutated amino acids, and that the nature of the effect is electrostatic. They also inferred an atomistic mechanism related to removal of unfavorable interactions, rather than to the formation of salt bridges. Molecular dynamics allows us to confirm and support both hypotheses. Several other collective parameters have also been monitored in relation to thermophilicity, such as global and local rigidity, permanence and number of hydrogen bonds, or of salt links. None of these clearly correlates with the thermal stability of the presently studied proteins.

Bacillus↗
PubMed2006

Molecular dynamics simulations on the inhibition of cyclin-dependent kinases 2 and 5 in the presence of activators.

Interests in CDK2 and CDK5 have stemmed mainly from their association with cancer and neuronal migration or differentiation related diseases and the need to design selective inhibitors for these kinases. Molecular dynamics (MD) simulations have not only become a viable approach to drug design because of advances in computer technology but are increasingly an integral part of drug discovery processes. It is common in MD simulations of inhibitor/CDK complexes to exclude the activator of the CDKs in the structural models to keep computational time tractable. In this paper, we present simulation results of CDK2 and CDK5 with roscovitine using models with and without their activators (cyclinA and p25). While p25 was found to induce slight changes in CDK5, the calculations support that cyclinA leads to significant conformational changes near the active site of CDK2. This suggests that detailed and structure-based inhibitor design targeted at these CDKs should employ activator-included models of the kinases. Comparisons between P/CDK2/cyclinA/roscovitine and CDK5/p25/roscovitine complexes reveal differences in the conformations of the glutamine around the active sites, which may be exploited to find highly selective inhibitors with respect to CDK2 and CDK5.

Cyclin-Dependent Kinase 2↗
PubMed2006

Molecular dynamics simulation and automated docking of the pro-apoptotic Bax protein and its complex with a peptide designed from the Bax-binding domain of anti-apoptotic Ku70.

Bax, a multi-domain protein belonging to the large family of Bcl-2 proteins, has a pivotal role for the initiation of the cytochrome c-mediated apoptosis, a vital physiologic process to eliminate damaged or unwanted cells. In response to specific stimuli Bax translocates from cytosol to mitochondria outer membrane where a process of oligomerization occurs with pore formation through which cytochrome c and other death molecules escape. The pro-death action of Bax is regulated by the interaction with other pro-survival proteins. However, the conformational changes and the structural details necessary for homo and hetero interaction with other regulating proteins are largely unknown. This article reports a combined investigation of molecular dynamics (MD) simulation and automated docking that evidence the molecular regions of Bax involved in the binding with anti-apoptotic exapeptide (Bip) designed from Ku70, a subunit of the protein complex essential for non-homologous DNA repair but that inhibits also the Bax translocation to mitochondria. Since Bip suppresses apoptosis induced by several anti-cancer drugs, it appears relevant to achieve a better understanding of the structural and dynamical aspects that characterize the Bip-Bax complex in view of potential therapeutic implications. The present results show that the Bax region with the highest affinity for Bip is located in proximity of BH3 homology domain of Bax and also involves the alpha-helices 1 and 8. Moreover, the comparison of essential motions of Bax at 300 and 400 K before and after the formation of the complex with Bip evidences how the binding with the exa-peptide affects the collective motions of specific molecular districts of Bax considered to have functional relevance.

Amino Acid Sequence↗
PubMed2002

Comparative molecular dynamics simulations of amphotericin B-cholesterol/ergosterol membrane channels.

Amphotericin B (AmB) is a very effective anti-fungal polyene macrolide antibiotic whose usage is limited by its toxicity. Lack of a complete understanding of AmB's molecular mechanism has impeded attempts to design less toxic AmB derivatives. The antibiotic is known to interact with sterols present in the cell membrane to form ion channels that disrupt membrane function. The slightly higher affinity of AmB toward ergosterol (dominant sterol in fungal cells) than cholesterol (mammalian sterol) is regarded as the most essential factor on which antifungal chemotherapy is based. To study these differences at the molecular level, two realistic model membrane channels containing molecules of AmB, sterol (cholesterol or ergosterol), phospholipid, and water were studied by molecular dynamics (MD) simulations. Comparative analysis of the simulation data revealed that the sterol type has noticeable effect on the properties of AmB membrane channels. In addition to having a larger size, the AmB channel in the ergosterol-containing membrane has a more pronounced pattern of intermolecular hydrogen bonds. The interaction between the antibiotic and ergosterol is more specific than between the antibiotic and cholesterol. These observed differences suggest that the channel in the ergosterol-containing membrane is more stable and, due to its larger size, would have a higher ion conductance. These observations are in agreement with experiments.

Amphotericin B↗
PubMed2006

Unfolding transition state and intermediates of the tumor suppressor p16INK4a investigated by molecular dynamics simulations.

The ankyrin repeat is one of the most common protein motifs and is involved in protein-protein interactions. It consists of 33 residues that assume a beta-hairpin helix-loop-helix fold. Mutagenesis and kinetic experiments (Phi-value analysis of the folding transition state) have shown that the tumor suppressor p16(INK4a), a four-repeat protein, unfolds sequentially starting from the two N-terminal repeats. Here, the flexibility of p16(INK4a) at room temperature and its unfolding mechanism at high temperature have been investigated by multiple molecular dynamics runs in explicit water for a total simulation time of 0.65 micros. The transition state ensemble (TSE) of p16(INK4a) was identified by monitoring both the deviation from the experimental Phi values and sudden conformational changes along the unfolding trajectories. Conformations in the TSE have a mainly unstructured second repeat whereas the other repeats are almost completely folded. A rigid-body displacement of the first repeat involving both a rotation and translation is observed in all molecular dynamics simulations at high temperature. The Trp(15), Pro(75), and Ala(76) side-chains are more buried in the TSE than the native state. The sequential unfolding starting at the second repeat is in agreement with the mutagenesis studies whereas the displacement of the first repeat and the presence of nonnative interactions at the TSE are simulation results which supplement the experimental data. Furthermore, the unfolding trajectories reveal the presence of two on-pathway intermediates with partial alpha-helical structure. Finally, on the basis of the available experimental and simulation results we suggest that in modular proteins the shift of the folding TSE toward the native structure upon reduction of the number of tandem repeats is consistent with the Hammond effect.

Binding Sites↗
PubMed2002

Structural analysis of human CCR2b and primate CCR2b by molecular modeling and molecular dynamics simulation.

CCR2b, a chemokine receptor for MCP-1, -2, -3, -4, plays an important role in a variety of diseases involving infection, inflammation, and/or injury, as well as being a coreceptor for HIV-1 infection. Two models of human CCR2b (hCCR2b) were generated by homology modeling and 1 ns restrained molecular dynamics (MD) simulation. In one only C113-C190 forms a disulfide bond (SS model); in another the potential C32-C277 disulfide bond was formed (2SS model). Analysis of the structures and averaged displacements of Calpha atoms of the N-terminal residues shows that the main differences between the SS and 2SS models lie in a region D25YDYGAPCHKFD36; in the extracellular part of the 2SS model the accessible surfaces of N12, F23, Y26, Y28 and F35 are obviously raised and a more stable H-bond net is formed. The potential energy of the 2SS-water assembly finally fluctuated around -43,020 kJ x mol(-1), which is about 302 kJ x mol(-1) lower than that of the SS-water assembly. All these results suggest that the 2SS model is more favorable. The CCR2b genes of 17 primates were sequenced and four CCR2b models for primates Ateles paniscus (A. pan), Hylobates leucogyneus(H. leu), Papio cynocephalus (P. cyn) and Trachypithecus francoist ( T. fra) were generated based on the 2SS model. A comparison of hCCR2b with primate CCR2b also supports the importance of the region D25YDYGAPCHKFD36. Electrostatic potential maps of human and primate CCR2b all display the dipolar characteristics of CCR2b with the negative pole located in the extracellular part and a strong positive pole in the cytoplasmic part. Based on the CCR2b model, we suggest that the main functional residues fall in the D25YDYGAPCHKFD36 region, and the negative electrostatic feature is a non-specific, but necessary, factor for ligands or gp120/CD4 binding.

Amino Acid Sequence↗
PubMed2000

Effect of the wall roughness on slip and rheological properties of hexadecane in molecular dynamics simulation of couette shear flow between two sinusoidal walls

Molecularly thin liquid films of alkanes in extreme conditions in a boundary lubrication regime have been investigated. The wall is modeled as a rough atomic sinusoidal wall. The effect on the boundary condition of the roughness characteristics, given by the period and amplitude of the sinusoidal wall, is studied here. The effect of the molecular length of the lubricating fluid is also examined here. The results show that the relative size of the fluid molecules and wall roughness determines the slip or nonslip boundary conditions. The effect of wall roughness characteristics on the rheological properties of the lubrication film is also studied.

Journal Article↗