Stocking of captive-bred fish can cause long-term population decline and gene pool replacement: predictions from a population dynamics model incorporating density-dependent mortality

Releasing captive-bred fish into natural environments (stocking) is common in fisheries worldwide. Although stocking is believed to have a positive effect on fish abundance over the short term, little is known about the long-term consequences of recurrent stocking and its influence on natural populations. In fact, there are growing concerns that genetically maladapted captive-bred fish can eventually reduce the abundance of natural population. In this study, we develop a simple model to quantitatively investigate the condition under which recurrent stocking has long-term effects on the natural population. Using a population dynamics model that takes into account a density-dependent recruitment, a gene responsible for the fitness difference between wild and captive-bred fish, and hybridization between them, we show that there is little or no contribution of recurrent stocking to the stock enhancement without a replacement of the wild gene pool by the captive-bred gene pool. The model further predicted that stocking of an intermediate level causes a reduction, rather than enhancement, of population size over the long term. The population decline due to stocking was attributed to the fitness disadvantage of captive-bred fish and strong overcompensation at recruitment stage. These results suggest that it would be difficult to simultaneously attain population size recovery and conservation of the local gene pool when captive-bred fish have fitness disadvantage in the wild, although caution is needed when applying the predictions from the simplified model to a specific species or population.

Supramolecular Organization of Heptapyrenotide Oligomers—An in Depth Investigation by Molecular Dynamics Simulations

We present a molecular modeling study based on ab initio and classical molecular dynamics calculations, for the investigation of the tridimensional structure and supramolecular assembly formation of heptapyrenotide oligomers in water solution. Our calculations show that free oligomers self-assemble in helical structures characterized by an inner core formed by π- stacked pyrene units, and external grooves formed by the linker moieties. The coiling of the linkers has high ordering, dominated by hydrogen-bond interactions among the phosphate and amide groups. Our models support a mechanism of longitudinal supramolecular oligomerization based on interstrand pyrene intercalation. Only a minimal number of pyrene units intercalate at one end, favoring formation of very extended longitudinal chains, as also detected by AFM experiment. Our results provide a structural explanation of the mechanism of chirality amplification in 1:1 mixtures of standard heptapyrenotides and modified oligomers with covalently linked deoxycytidine, based on selective molecular recognition and binding of the nucleotide to the groove of the left-wound helix.

Unraveling the Conformational Determinants of Peptide Dendrimers Using Molecular Dynamics Simulations

Peptide dendrimers are synthetic tree-like molecules composed of amino acids. There are at least two kinds of preferential structural behaviors exhibited by these molecules, which acquire either compact or noncompact shapes. However, the key structural determinants of such behaviors remained, until now, unstudied. Herein, we conduct a comprehensive investigation of the structural determinants of peptide dendrimers by employing long molecular dynamics simulations to characterize an extended set of third generation dendrimers. Our results clearly show that a trade-off between electrostatic effects and hydrogen bond formation controls structure acquisition in these systems. Moreover, by selectively changing the dendrimers charge we are able to manipulate the exhibited compactness. In contrast, the length of branching residues does not seem to be a major structural determinant. Our results are in accordance with the most recent experimental evidence and shed some light on the key molecular level interactions controlling structure acquisition in these systems. Thus, the results presented constitute valuable insights that can contribute to the development of truly tailor-made dendritic systems.

Mechanisms of recognition and binding of α-TTP to the plasma membrane by multi-scale molecular dynamics simulations

We used multiple sets of simulations both at the atomistic and coarse-grained level of resolution to investigate interaction and binding of α-tochoperol transfer protein (α-TTP) to phosphatidylinositol phosphate lipids (PIPs). Our calculations indicate that enrichment of membranes with such lipids facilitate membrane anchoring. Atomistic models suggest that PIP can be incorporated into the binding cavity of α-TTP and therefore confirm that such protein can work as lipid exchanger between the endosome and the plasma membrane. Comparison of the atomistic models of the α-TTP-PIPs complex with membrane-bound α-TTP revealed different roles for the various basic residues composing the basic patch that is key for the protein/ligand interaction. Such residues are of critical importance as several point mutations at their position lead to severe forms of ataxia with vitamin E deficiency (AVED) phenotypes. Specifically, R221 is main residue responsible for the stabilization of the complex. R68 and R192 exchange strong interactions in the protein or in the membrane complex only, suggesting that the two residues alternate contact formation, thus facilitating lipid flipping from the membrane into the protein cavity during the lipid exchange process. Finally, R59 shows weaker interactions with PIPs anyway with a clear preference for specific phosphorylation positions, hinting a role in early membrane selectivity for the protein. Altogether, our simulations reveal significant aspects at the atomistic scale of interactions of α-TTP with the plasma membrane and with PIP, providing clarifications on the mechanism of intracellular vitamin E trafficking and helping establishing the role of key residue for the functionality of α-TTP.

Single Molecular Conductance of Tolanes: Experimental and Theoretical Study on the Junction Evolution Dependent on the Anchoring Group

Employing a scanning tunneling microscopy based beak junction technique and mechanically controlled break junction experiments, we investigated tolane (diphenylacetylene)-type single molecular junctions having four different anchoring groups (SH, pyridyl (PY), NH2, and CN) at a solid/liquid interface. The combination of current–distance and current–voltage measurements and their quantitative statistical analysis revealed the following sequence for junction formation probability and stability: PY > SH > NH2 > CN. For all single molecular junctions investigated, we observed the evolution through multiple junction configurations, with a particularly well-defined binding geometry for PY. The comparison of density functional theory type model calculations and molecular dynamics simulations with the experimental results revealed structure and mechanistic details of the evolution of the different types of (single) molecular junctions upon stretching quantitatively.

Atomistic simulations of 2D bicomponent self-assembly: from molecular recognition to self-healing

Supramolecular two-dimensional engineering epitomizes the design of complex molecular architectures through recognition events in multicomponent self-assembly. Despite being the subject of in-depth experimental studies, such articulated phenomena have not been yet elucidated in time and space with atomic precision. Here we use atomistic molecular dynamics to simulate the recognition of complementary hydrogen-bonding modules forming 2D porous networks on graphite. We describe the transition path from the melt to the crystalline hexagonal phase and show that self-assembly proceeds through a series of intermediate states featuring a plethora of polygonal types. Finally, we design a novel bicomponent system possessing kinetically improved self-healing ability in silico, thus demonstrating that a priori engineering of 2D self-assembly is possible.

A novel computer simulation method for simulating the multiscale transduction dynamics of signal proteins

Signal proteins are able to adapt their response to a change in the environment, governing in this way a broad variety of important cellular processes in living systems. While conventional molecular-dynamics (MD) techniques can be used to explore the early signaling pathway of these protein systems at atomistic resolution, the high computational costs limit their usefulness for the elucidation of the multiscale transduction dynamics of most signaling processes, occurring on experimental timescales. To cope with the problem, we present in this paper a novel multiscale-modeling method, based on a combination of the kinetic Monte-Carlo- and MD-technique, and demonstrate its suitability for investigating the signaling behavior of the photoswitch light-oxygen-voltage-2-Jα domain from Avena Sativa (AsLOV2-Jα) and an AsLOV2-Jα-regulated photoactivable Rac1-GTPase (PA-Rac1), recently employed to control the motility of cancer cells through light stimulus. More specifically, we show that their signaling pathways begin with a residual re-arrangement and subsequent H-bond formation of amino acids near to the flavin-mononucleotide chromophore, causing a coupling between β-strands and subsequent detachment of a peripheral α-helix from the AsLOV2-domain. In the case of the PA-Rac1 system we find that this latter process induces the release of the AsLOV2-inhibitor from the switchII-activation site of the GTPase, enabling signal activation through effector-protein binding. These applications demonstrate that our approach reliably reproduces the signaling pathways of complex signal proteins, ranging from nanoseconds up to seconds at affordable computational costs.

IGHD II: A Novel GH-1 Gene Mutation (GH-L76P) Severely Affects GH Folding, Stability, and Secretion

CONTEXT The autosomal dominant form of GH deficiency (IGHD II) is characterized by markedly reduced GH secretion combined with low concentrations of IGF-1 leading to short stature. OBJECTIVE Structure-function analysis of a missense mutation in the GH-1 gene converting codon 76 from leucine (L) to proline (P) yielding a mutant GH-L76P peptide. DESIGN, SETTINGS, AND PATIENTS Heterozygosity for GH-L76P/wt-GH was identified in a nonconsanguineous Spanish family. The index patients, two siblings, a boy and a girl, were referred for assessment of their short stature (-3.2 and -3.8 SD). Their grandmother, father, and aunt were also carrying the same mutation and showed severe short stature; therefore, IGHD II was diagnosed. INTERVENTIONS AND RESULTS AtT-20 cells coexpressing both wt-GH and GH-L76P showed a reduced GH secretion (P < .001) after forskolin stimulation when compared with the cells expressing only wt-GH. In silico mutagenesis and molecular dynamics simulations presented alterations of correct folding and mutant stability compared with wt-GH. Therefore, further structural analysis of the GH-L76P mutant was performed using expressed and purified proteins in Escherichia coli by thermofluor assay and fast degradation proteolysis assay. Both assays revealed that the GH-L76P mutant is unstable and misfolded compared to wt-GH confirming the bioinformatic model prediction. CONCLUSIONS This is the first report of a family suffering from short stature caused by IGHD II, which severely affects intracellular GH folding and stability as well as secretion, highlighting the necessity of functional analysis of any GH variant for defining new mechanisms as a cause for IGHD II.