PubMed Health⌕ Search

Biomedical subjects

Riccardo Baron

Publications and source records attributed to Riccardo Baron.

8 recordsLinked to original sources

Configurational entropies of lipids in pure and mixed bilayers from atomic-level and coarse-grained molecular dynamics simulations.

Single-chain and single-fragment configurational entropies of lipid tails in hydrated lipid bilayers are evaluated from molecular dynamics simulations using the quasi-harmonic approximation. The entropy distribution along individual acyl tails is obtained and compared to that of corresponding hydrocarbon chains in the liquid phase. We consider pure dipalmitoylphosphatidylcholine and mixed dioleoylphosphatidylcholine/dioleoylphosphatidylethanolamine bilayers. The systems are modeled at different levels of spatial resolution: In an atomic-level (AL) model all (heavy) atoms are explicitly simulated; in a coarse-grained (CG) model particles (beads) representing groups of covalently bound atoms are used, which map approximately four non-hydrogen atoms to one interaction site. Single-chain and single-fragment entropies and correlations between the motions of (single) acyl chains are compared. A good correspondence is found between the flexibility of the AL and CG models. The loss in configurational entropy due to the reduction in the number of degrees of freedom upon coarse-graining of the model is estimated. The CG model shows about 4 times faster convergence of the chain entropies than the more detailed AL model. Corrections to the quasi-harmonic entropy estimates were found to be small for the CG model. For the AL model, the correction due to mode anharmonicities is small, but the correction due to pairwise (supralinear) mode correlations is sizable.

Computer Simulation↗

Biomolecular modeling: Goals, problems, perspectives.

Computation based on molecular models is playing an increasingly important role in biology, biological chemistry, and biophysics. Since only a very limited number of properties of biomolecular systems is actually accessible to measurement by experimental means, computer simulation can complement experiment by providing not only averages, but also distributions and time series of any definable quantity, for example, conformational distributions or interactions between parts of systems. Present day biomolecular modeling is limited in its application by four main problems: 1) the force-field problem, 2) the search (sampling) problem, 3) the ensemble (sampling) problem, and 4) the experimental problem. These four problems are discussed and illustrated by practical examples. Perspectives are also outlined for pushing forward the limitations of biomolecular modeling.

Computer Simulation↗

Configurational entropy change of netropsin and distamycin upon DNA minor-groove binding.

Binding of a small molecule to a macromolecular target reduces its conformational freedom, resulting in a negative entropy change that opposes the binding. The goal of this study is to estimate the configurational entropy change of two minor-groove-binding ligands, netropsin and distamycin, upon binding to the DNA duplex d(CGCGAAAAACGCG).d(CGCGTTTTTCGCG). Configurational entropy upper bounds based on 10-ns molecular dynamics simulations of netropsin and distamycin in solution and in complex with DNA in solution were estimated using the covariance matrix of atom-positional fluctuations. The results suggest that netropsin and distamycin lose a significant amount of configurational entropy upon binding to the DNA minor groove. The estimated changes in configurational entropy for netropsin and distamycin are -127 J K(-1) mol(-1) and -104 J K(-1) mol(-1), respectively. Estimates of the configurational entropy contributions of parts of the ligands are presented, showing that the loss of configurational entropy is comparatively more pronounced for the flexible tails than for the relatively rigid central body.

Binding Sites↗

Comparison of atomic-level and coarse-grained models for liquid hydrocarbons from molecular dynamics configurational entropy estimates.

Molecular liquids can be modeled at different levels of spatial resolution. In atomic-level (AL) models, all (heavy) atoms can be explicitly simulated. In coarse-grained (CG) models, particles (beads) that represent groups of covalently bound atoms are used as elementary units. Ideally, a CG model should reproduce the thermodynamic and structural properties of the corresponding AL model after mapping to the lower-resolution scale. In the present work, two such models are investigated: (i) the classical GROMOS atomic-level model; (ii) a CG model recently proposed by Marrink et al., which maps approximately four non-hydrogen atoms to one bead [J. Phys. Chem. B 2004, 108, 750]. The study is restricted to n-alkanes whose aliphatic fragments are abundantly found in lipids of biological interest. Additionally, cis-9-octadecene is included, as a template chain of the lipid dioleoylphosphatidylcholine (DOPC). The two representations of molecules in the liquid phase are compared in terms of average molecular structures, extent of configurational space sampled, and single-molecule entropies. An approximate method is used to estimate the rotational contributions to the absolute configurational entropy. Good correspondence between the AL and CG representations is found. The loss in configurational entropy due to the reduction in degrees of freedom upon coarse-graining of the model is estimated.

Algorithms↗

Conformational and dynamical properties of disaccharides in water: a molecular dynamics study.

Explicit-solvent molecular dynamics simulations (50 ns, 300 K) of the eight reducing glucose disaccharides (kojibiose, sophorose, nigerose, laminarabiose, maltose, cellobiose, isomaltose, and gentiobiose) have been carried out using the GROMOS 45A4 force field (including a recently reoptimized carbohydrate parameter set), to investigate and compare their conformational preferences, intramolecular hydrogen-bonding patterns, torsional dynamics, and configurational entropies. The calculated average values of the glycosidic torsional angles agree well with available experimental data, providing validation for the force field and simulation methodology employed in this study. These simulations show in particular that: 1) (1-->6)-linked disaccharides are characterized by an increased flexibility, the absence of any persistent intramolecular hydrogen bond and a significantly higher configurational entropy (compared to the other disaccharides); 2) cellobiose presents a highly persistent interresidue hydrogen bond and a significantly lower configurational entropy (compared to the other disaccharides); 3) persistent hydrogen bonds are observed for all disaccharides (except (1-->6)-linked) and typically involve a hydrogen donor in the reducing residue and an acceptor in the nonreducing one; 4) the probability distributions associated with the glycosidic dihedral angles and psi are essentially unimodal for all disaccharides, and full rotation around these angles occurs at most once or twice for (never for psi) on the 50-ns timescale; and 5) the timescales associated with torsional transitions (except around and psi) range from approximately 30 ps (rotation of hydroxyl groups) to the nanosecond range (rotation of the lactol and hydroxymethyl groups, and around the omega-glycosidic dihedral angle in (1-->6)-linked disaccharides).

Carbohydrate Sequence↗

The GROMOS software for biomolecular simulation: GROMOS05.

We present the latest version of the Groningen Molecular Simulation program package, GROMOS05. It has been developed for the dynamical modelling of (bio)molecules using the methods of molecular dynamics, stochastic dynamics, and energy minimization. An overview of GROMOS05 is given, highlighting features not present in the last major release, GROMOS96. The organization of the program package is outlined and the included analysis package GROMOS++ is described. Finally, some applications illustrating the various available functionalities are presented.

Algorithms↗

Principles of carbopeptoid folding: a molecular dynamics simulation study.

The conformational spaces of five oligomers of tetrahydrofuran-based carbopeptoids in chloroform and dimethyl sulfoxide were investigated through nine molecular dynamics simulations. Prompted by nuclear magnetic resonance experiments that indicated various stable folds for some but not all of these carbopeptoids, their folding behaviour was investigated as a function of stereochemistry, chain length and solvent. The conformational distributions of these molecules were analysed in terms of occurrence of hydrogen bonds, backbone torsional-angle distributions, conformational clustering and solute configurational entropy. While a cis-linkage across the tetrahydrofuran ring favours right-handed helical structures, a trans-linkage results in a larger conformational variability. Intra-solute hydrogen bonding is reduced with increasing chain length and with increasing solvent polarity. Solute configurational entropies confirm the picture obtained: they are smaller for cis- than for trans-linked peptides, for chloroform than for dimethyl sulfoxide as solvent and for shorter peptide chains. The simulations provide an atomic picture of molecular conformational variability that is consistent with the available experimental data.

Computer Simulation↗